Abstract:Abstract. We apply three-dimensional (3-D) thermo-mechanical numerical simulations of the shortening of the upper crustal region of a passive margin in order to investigate the control of 3-D laterally variable inherited structures on fold-and-thrust belt evolution and associated nappe formation. We consider tectonic inheritance by employing an initial model configuration with basement horst and graben structures having laterally variable geometry and with sedimentary layers having different mechanical strengt… Show more
“…This way, the generated passive margin, the marine basin system and the subsequent convergence are modelled in an internally consistent manner. Subduction initiation is horizontally forced (Stern, 2004;Stern and Gerya, 2018;Crameri et al, 2020), and a major lithospheric shear zone forms around the transition from the distal to the proximal margin (see also the discussion in Candioti et al, 2020). The ad hoc parameterized layer of serpentinite is not relevant for subduction initiation, as subduction is initiated also in the reference model without serpentinite (see Fig.…”
Section: Subduction Initiationmentioning
confidence: 99%
“…Such wedge models have also been used to study the formation of viscous fold nappes during fold and thrust belt evolution (e.g. Bauville and Schmalholz, 2015;Spitz et al, 2020) or the impact of surface processes on wedge formation (e.g. Willett, 1999).…”
Abstract. The dynamics of growing collisional orogens are mainly controlled by buoyancy and shear forces. However, the relative importance of these forces, their temporal evolution and their impact on the tectonic style of orogenic wedges remain elusive.
Here, we quantify buoyancy and shear forces during collisional orogeny and investigate their impact on orogenic wedge formation and exhumation of crustal rocks.
We leverage two-dimensional petrological–thermomechanical numerical simulations of a long-term (ca. 170 Myr) lithosphere deformation cycle involving subsequent hyperextension, cooling, convergence, subduction and collision.
Hyperextension generates a basin with exhumed continental mantle bounded by asymmetric passive margins.
Before convergence, we replace the top few kilometres of the exhumed mantle with serpentinite to investigate its role during subduction and collision. We study the impact of three parameters: (1) shear resistance, or strength, of serpentinites, controlling the strength of the evolving subduction interface; (2) strength of the continental upper crust; and (3) density structure of the subducted material.
Densities are determined by linearized equations of state or by petrological-phase equilibria calculations.
The three parameters control the evolution of the ratio of upward-directed buoyancy force to horizontal driving force, FB/FD=ArF, which controls the mode of orogenic wedge formation: ArF≈0.5 causes thrust-sheet-dominated wedges, ArF≈0.75 causes minor wedge formation due to relamination of subducted crust below the upper plate, and ArF≈1 causes buoyancy-flow- or diapir-dominated wedges involving exhumation of crustal material from great depth (>80 km).
Furthermore, employing phase equilibria density models reduces the average topography of wedges by several kilometres. We suggest that during the formation of the Pyrenees ArF⪅0.5 due to the absence of high-grade metamorphic rocks, whereas for the Alps ArF≈1 during exhumation of high-grade rocks and ArF⪅0.5 during the post-collisional stage.
In the models, FD increases during wedge growth and subduction and eventually reaches magnitudes (≈18 TN m−1) which are required to initiate subduction.
Such an increase in the horizontal force, required to continue driving subduction, might have “choked” the subduction of the European plate below the Adriatic one between 35 and 25 Ma and could have caused the reorganization of plate motion and subduction initiation of the Adriatic plate.
“…This way, the generated passive margin, the marine basin system and the subsequent convergence are modelled in an internally consistent manner. Subduction initiation is horizontally forced (Stern, 2004;Stern and Gerya, 2018;Crameri et al, 2020), and a major lithospheric shear zone forms around the transition from the distal to the proximal margin (see also the discussion in Candioti et al, 2020). The ad hoc parameterized layer of serpentinite is not relevant for subduction initiation, as subduction is initiated also in the reference model without serpentinite (see Fig.…”
Section: Subduction Initiationmentioning
confidence: 99%
“…Such wedge models have also been used to study the formation of viscous fold nappes during fold and thrust belt evolution (e.g. Bauville and Schmalholz, 2015;Spitz et al, 2020) or the impact of surface processes on wedge formation (e.g. Willett, 1999).…”
Abstract. The dynamics of growing collisional orogens are mainly controlled by buoyancy and shear forces. However, the relative importance of these forces, their temporal evolution and their impact on the tectonic style of orogenic wedges remain elusive.
Here, we quantify buoyancy and shear forces during collisional orogeny and investigate their impact on orogenic wedge formation and exhumation of crustal rocks.
We leverage two-dimensional petrological–thermomechanical numerical simulations of a long-term (ca. 170 Myr) lithosphere deformation cycle involving subsequent hyperextension, cooling, convergence, subduction and collision.
Hyperextension generates a basin with exhumed continental mantle bounded by asymmetric passive margins.
Before convergence, we replace the top few kilometres of the exhumed mantle with serpentinite to investigate its role during subduction and collision. We study the impact of three parameters: (1) shear resistance, or strength, of serpentinites, controlling the strength of the evolving subduction interface; (2) strength of the continental upper crust; and (3) density structure of the subducted material.
Densities are determined by linearized equations of state or by petrological-phase equilibria calculations.
The three parameters control the evolution of the ratio of upward-directed buoyancy force to horizontal driving force, FB/FD=ArF, which controls the mode of orogenic wedge formation: ArF≈0.5 causes thrust-sheet-dominated wedges, ArF≈0.75 causes minor wedge formation due to relamination of subducted crust below the upper plate, and ArF≈1 causes buoyancy-flow- or diapir-dominated wedges involving exhumation of crustal material from great depth (>80 km).
Furthermore, employing phase equilibria density models reduces the average topography of wedges by several kilometres. We suggest that during the formation of the Pyrenees ArF⪅0.5 due to the absence of high-grade metamorphic rocks, whereas for the Alps ArF≈1 during exhumation of high-grade rocks and ArF⪅0.5 during the post-collisional stage.
In the models, FD increases during wedge growth and subduction and eventually reaches magnitudes (≈18 TN m−1) which are required to initiate subduction.
Such an increase in the horizontal force, required to continue driving subduction, might have “choked” the subduction of the European plate below the Adriatic one between 35 and 25 Ma and could have caused the reorganization of plate motion and subduction initiation of the Adriatic plate.
“…First, high-resolution 3-D modelling requires assembling large matrices that are challenging to solve even with the sparse matrix format (Burov et al 2014;May et al 2014;Morra 2020). As a result, the maximum degree of freedom is usually limited to <10 9 in published 3-D models (Zhong et al 2000;Settari & Walters 2001;Tackley 2008;May et al 2014;Moresi et al 2014;Dannberg & Heister 2016;Spitz et al 2020;Rinaldi et al 2021). Secondly, existing porous flow models in deformable rocks commonly use simplified equations that might give inaccurate results (Settari & Walters 2001;Prevost 2014).…”
Section: Numerical Methods In Geosciencesmentioning
Summary
Two-phase flow equations that couple solid deformation and fluid migration have opened new research trends in geodynamical simulations and modelling of subsurface engineering. Physical nonlinearity of fluid-rock systems and strong coupling between flow and deformation in such equations lead to interesting predictions such as spontaneous formation of focused fluid flow in ductile/plastic rocks. However, numerical implementation of two-phase flow equations and their application to realistic geological environments with complex geometries and multiple stratigraphic layers is challenging. This study documents an efficient pseudo-transient solver for two-phase flow equations and describes the numerical theory and physical rationale. We provide a simple explanation for all steps involved in the development of a pseudo-transient numerical scheme for various types of equations. Two different constitutive models are used in our formulations: a bilinear viscous model with decompaction weakening and a viscoplastic model that allows decompaction weakening at positive effective pressures. The resulting numerical models are used to study fluid leakage from high porosity reservoirs into less porous overlying rocks. The interplay between time-dependent rock deformation and the buoyancy of ascending fluids leads to the formation of localized channels. The role of material parameters, reservoir topology, geological heterogeneity and porosity is investigated. Our results show that material parameters control the propagation speed of channels while the geometry of the reservoir controls their locations. Geological layers present in the overburden do not stop the propagation of the localized channels but rather modify their width, permeability, and growth speed.
“…These field observations are key to characterize the deformation within the orogenic wedge, and thus to discuss the outcome of rheological models (e.g. Spitz et al, 2020;Bauville et Schmalholz, 2015). Helvetic thrust, Alpine Sole Thrust, and late thrusts.…”
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.