SUMMARY We present 2-D numerical simulations of convergence at a hyperextended passive margin with exhumed subcontinental mantle. We consider viscoelasto-plastic deformation, heat transfer and thermomechanical coupling by shear heating and associated thermal softening due to temperature dependent viscosity. The simulations show subduction initiation for convergence velocities of 2 cm yr−1, initial Moho temperatures of 525 °C and maximal deviatoric stresses of ca. 800 MPa, around the Moho, prior to localization. Subduction initiates in the region with thinned continental crust and is controlled by a thermally activated ductile shear zone in the mantle lithosphere. The shear zone temperature can be predicted with a recently published analytical expression. The criterion for subduction initiation is a temperature difference of at least 225 °C between predicted temperature and initial Moho temperature. The modelled forced subduction broadly agrees with geological data and reconstructions of subduction during closure of the Piemont-Liguria basin, caused by convergence of the European and Adriatic plates during the Alpine orogeny.
Abstract. Many plate tectonic processes, such as subduction initiation, are embedded in long-term (>100 Myr) geodynamic cycles often involving subsequent phases of extension, cooling without plate deformation and convergence. However, the impact of upper mantle convection on lithosphere dynamics during such long-term cycles is still poorly understood. We have designed two-dimensional upper-mantle-scale (down to a depth of 660 km) thermo-mechanical numerical models of coupled lithosphere–mantle deformation. We consider visco–elasto–plastic deformation including a combination of diffusion, dislocation and Peierls creep law mechanisms. Mantle densities are calculated from petrological phase diagrams (Perple_X) for a Hawaiian pyrolite. Our models exhibit realistic Rayleigh numbers between 106 and 107, and the model temperature, density and viscosity structures agree with geological and geophysical data and observations. We tested the impact of the viscosity structure in the asthenosphere on upper mantle convection and lithosphere dynamics. We also compare models in which mantle convection is explicitly modelled with models in which convection is parameterized by Nusselt number scaling of the mantle thermal conductivity. Further, we quantified the plate driving forces necessary for subduction initiation in 2D thermo-mechanical models of coupled lithosphere–mantle deformation. Our model generates a 120 Myr long geodynamic cycle of subsequent extension (30 Myr), cooling (70 Myr) and convergence (20 Myr) coupled to upper mantle convection in a single and continuous simulation. Fundamental features such as the formation of hyperextended margins, upper mantle convective flow and subduction initiation are captured by the simulations presented here. Compared to a strong asthenosphere, a weak asthenosphere leads to the following differences: smaller value of plate driving forces necessary for subduction initiation (15 TN m−1 instead of 22 TN m−1) and locally larger suction forces. The latter assists in establishing single-slab subduction rather than double-slab subduction. Subduction initiation is horizontally forced, occurs at the transition from the exhumed mantle to the hyperextended passive margin and is caused by thermal softening. Spontaneous subduction initiation due to negative buoyancy of the 400 km wide, cooled, exhumed mantle is not observed after 100 Myr in model history. Our models indicate that long-term lithosphere dynamics can be strongly impacted by sub-lithosphere dynamics. The first-order processes in the simulated geodynamic cycle are applicable to orogenies that resulted from the opening and closure of embryonic oceans bounded by magma-poor hyperextended rifted margins, which might have been the case for the Alpine orogeny.
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.
Abstract. Orogenic belts formed by collision are impressive manifestations of plate tectonics. Observations from orogenic belts, like the Western Alps, indicate an important involvement of the mantle lithosphere, significant burial and exhumation of continental and oceanic crustal rocks and the importance of the plate interface strength that can be modified, for example, by the presence of serpentinites. A popular model for the formation of such belts is the so-called orogenic wedge model. However, most wedge models consider crustal deformation only and do, hence, not consider subduction, the impact of related buoyancy forces arising from density differences between subducted crust and surrounding mantle and the effects of different plate interface strength. Here, we quantify the relative importance of buoyancy and shear forces in building collisional orogenic wedges. We leverage two-dimensional (2D) petrological-thermo-mechanical numerical simulations of a long-term (ca. 170 Myr) lithosphere deformation cycle involving subsequent hyperextension, cooling, convergence, subduction and collision. We compare simulations employing density fields calculated with linearized equations of state with simulations employing density fields calculated by phase equilibria models including metamorphic reactions. Further, we consider serpentinisation of the mantle material, exhumed in the hyperextended basin. Our models show that differences in density structure and in shear strength of serpentinites or upper crust have a strong impact on the evolution of orogenic wedges. Higher serpentinite strength causes a dominance of shear over buoyancy forces, resulting in either thrust-sheet dominated orogenic wedges, involving some diapiric exhumation at their base, or relamination of crustal material below the overriding plate. Lower serpentinite strength (equal importance of shear and buoyancy forces) generates orogenic wedges that are dominated by diapiric or channel-flow exhumation. Deep subduction (> 80 km) and subsequent surface exhumation of continental crust along the subduction interface occurs in these models. Employing phase equilibria density models decreases the average buoyancy contrasts, allows for deeper subduction of continental crust and reduces the average topography of the wedge by several kilometers. A decrease of upper crustal shear strength causes smaller maximal crustal burial depths. Progressive subduction of continental crust increases upward-directed buoyancy forces of the growing wedge and in turn increases horizontal driving forces. These driving forces eventually reach magnitudes (≈ 18 TN m−1) which were required to initiate subduction during convergence. We suggest that the evolving relation between shear and buoyancy forces and the increase of horizontal driving force related to the growing Alpine orogenic wedge has significantly slowed down (or choked) subduction of the European plate below the Adriatic one between 35 and 25 Ma. This buoyancy-related choking could have caused the reorganization of plate motion and the initiation of subduction of the Adriatic plate. We discuss potential applications and implications of our model results to the Pyrenean and Alpine orogenies.
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.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.