An official website of the United States government
Plate tectonic reconstructions of three of the best-defined Cenozoic subduction initiation (SI) events in the western Pacific, Izu-Bonin-Mariana, Vanuatu, and Puysegur subduction zones, show substantial components of strike-slip motion before and during the subduction initiation. Using computational models, we show that strike-slip motion has a large influence on the effective strength of incipient margins and the ease of subduction initiation. The parameter space associated with visco-elasto-plastic rheologies, plate weakening, and plate forces and kinematics is explored and we show that subduction initiates more easily with a higher force, a faster weakening, or greater strike-slip motion. With the analytical solution, we demonstrate that the effect of strike-slip motion can be equivalently represented by a modified weakening rate. Along transpressive margins, we show that a block of oceanic crust can become trapped between a new thrust fault and the antecedent strike-slip fault and is consistent with structural reconstructions and gravity models of the Puysegur margin. Together, models and observations suggest that subduction initiation can be triggered when margins become progressively weakened to the point that the resisting forces become smaller than the driving forces, and as the negative buoyancy builds up, the intraplate stress eventually turns from compressional into extensional. The analytical formulation of the initiation time,tSI, marking the moment when intraplate stress flips sign, is validated with a computational models. The analytical solution shows thattSIis dominated by convergence velocity, while the plate age, strike-slip velocity, and weakening rate all have a smaller but still important effect on the time scale of subduction initiation.
more »
« less
A simple force balance model of subduction initiation
SUMMARY The initiation and development of subduction zones are associated with substantial stress changes both within plates and at plate boundaries. We formulate a simple analytical model based on the force balance equation of a subduction zone, and validate it with numerical calculations of highly non-linear, coupled thermomechanical system. With two kinds of boundary conditions with either fixed velocity or fixed force in the far-field, we quantitatively analyse the role of each component in the force balance equation, including slab pull, interplate friction, plate bending and basal traction, on the kinematics and stress state of a subducting plate. Based on the numerical and analytical models, we discuss the evolution of plate curvature, the role of plastic yielding and elasticity, and how different factors affect the timing of subduction initiation. We demonstrate with the presence of plastic yielding for a plate of thickness, H, that the bending force is proportional to H2, instead of H3 as previously thought. Although elasticity increases the force required to start nucleating subduction it does not substantially change the total work required to initiate a subduction zone when the yielding stress is small. The analytical model provides an excellent fit to the total work and time to initiate subduction and the force and velocity as a function of convergence and time. Plate convergence and weakening rate during nucleation are the dominant factors influencing the force balance of the plate, and 200 km of plate convergence is typically required to bring a nascent subduction zone into a self-sustaining state. The closed-form solution now provides a framework to better interpret even more complex, time-dependent systems in three dimensions.
more »
« less
- Award ID(s):
- 2049086
- PAR ID:
- 10394676
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 232
- Issue:
- 1
- ISSN:
- 0956-540X
- Page Range / eLocation ID:
- 128 to 146
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Despite significant study, when and how plate tectonics initiated on Earth remains contentious. Geologic evidence from some of Earth's earliest cratons has been interpreted as reflecting the formation of initial continental blocks by non‐subduction processes, which then trigger subduction initiation at their margins. Numerical models of mantle convection with a plastic yield stress rheology have shown this scenario is plausible. However, whether continents can trigger subduction initiation has not been tested with other rheologies. We, therefore, use numerical models of mantle convection with an imposed continental block to test whether continents facilitate subduction initiation with a grain‐damage mechanism, where weak shear zones form by grain size reduction. Our results show that continents modestly enhance stresses in the lithosphere, but not enough to significantly impact lithospheric damage or subduction initiation: continents have minimal influence on lithospheric damage or plate speed, nor does subduction preferentially initiate at the continental margin. A new regime diagram that includes continental blocks shows only a small shift in the boundary between the mobile‐lid and stagnant‐lid regimes when continents are added. However, as we do find that stresses are modestly enhanced at the continental margin in our models, we develop a scaling law for this stress enhancement to more fully test whether continents could trigger subduction initiation on early Earth. We find that lithospheric stresses supplied by continents are not sufficient to initiate subduction on the early Earth on their own with grain‐damage rheology; instead, additional factors would be required.more » « less
-
Abstract Despite extensive modeling efforts, the dynamics of subduction initiation (SI), including the role of elasticity, are not fully understood. Using two‐dimensional thermomechanical models with visco‐plastic (VP) and visco‐elasto‐plastic (VEP) rheologies, we systematically investigate the role of elasticity in intraoceanic SI using two model setups: spontaneous initiation without imposed convergence and induced initiation with imposed convergence. In spontaneous models, the overriding plate age of <20 Ma and the subducting plate age of >50 Ma generally lead to vertically driven SI with either rheology, but for a given age contrast, SI is easier to occur with the VEP rheology. In induced model with either rheology, when the two plates are young and have a small age contrast, the resulting SI is horizontally driven, and elasticity does not affect SI significantly, regardless of the convergence rate. However, when the thermal age contrast is large and a convergence rate is relatively low, the SI in induced models is vertically driven and similar to that in the spontaneous models, and the VEP rheology leads to faster SI than the VP rheology. This effect of elasticity becomes smaller with increasing initial horizontal compressional stress but does not become fully negated by the initial stress of <∼50 MPa. Therefore, inclusion of elasticity with reasonable shear modulus and initial stress values results in a weaker slab, making it easier for vertically driven SI to occur when the age contrast is relatively large.more » « less
-
The accretion of future allochthonous terranes (e.g., microcontinents or oceanic plateaus) onto the southern margin of Asia occurred repeatedly during the evolution and closure of the Tethyan oceanic realm, but the specific geodynamic processes of this protracted convergence, successive accretion, and subduction zone initiation remain largely unknown. Here, we use numerical models to better understand the dynamics that govern multiple terrane accretions and the polarity of new subduction zone initiation. Our results show that the sediments surrounding the future terranes and the structural complexity of the overriding plate are important factors that affect accretion of multiple plates and guide subduction polarity. Wide (≥400 km) and buoyant terranes with sediments behind them and fast continental plate motions are favorable for multiple unidirectional subduction zone jumps, which are also referred to as subduction zone transference, and successive terrane-accretion events. The jumping times (∼3−20+ m.y.) are mainly determined by the convergence rates and rheology of the overriding complex plate with preceding terrane collisions, which increase with slower convergence rates and/or a greater number of preceding terrane collisions. Our work provides new insights into the key geodynamic conditions governing multiple subduction zone jumps induced by successive accretion and discusses Tethyan evolution at a macro level. More than 50 m.y. after India-Asia collision, subduction has yet to initiate along the southern Indian plate, which may be the joint result of slower plate convergence and partitioned deformation across southern Asia.more » « less
-
Subduction initiation remains one of the least understood global processes of plate tectonics. Prominent models have been cast in terms of two broad classes: spontaneous cases due to lithospheric gravitational instabilities and induced cases due to forced plate convergence. Yet gravitationally unstable lithosphere is old, strong, and difficult to begin to bend into a subduction zone and convergent forces necessary to begin subduction are often too large given the plates involved. These models also consider the asthenospheric mantle as passive, even though relative motion between slabs and the asthenosphere has long been regarded as a strong control on subduction dynamics. Here I propose that subduction-transform edge propagator (STEP) faults can initiate subduction depending on the absolute motion of the STEP fault with respect to the asthenosphere. STEP faults form where subduction zones end and the subducting plate tears forming a down flexed transcurrent plate boundary at the surface shearing against the adjacent rear arc lithospheric plate. However, STEP faults are not simple transcurrent boundaries. Absolute motion of the down flexed STEP fault edge with respect to the surrounding asthenosphere can produce a strong sea anchor force that either continues to bend the edge downward, initiating subduction, or opposes slab bending, inhibiting subduction. In the south Pacific, the southern end of the New Hebrides Trench and the northern end of the Tonga Trench are type-example STEP faults with opposite senses of dip but both moving northward with respect to the asthenosphere. The northward dipping New Hebrides STEP fault moves northward in a mantle reference frame creating a strong asthenospheric flow against the STEP fault edge, inducing active subduction at the Matthew-Hunter trench. In contrast, the Tonga STEP fault dips southward but also has a northward component of motion with respect to the mantle. Asthenosphere thus flows southward beneath the down flexed Tonga STEP fault edge opposing further bending. Subduction does not initiate at the Tonga STEP fault despite a ~100 Myr age contrast between the Pacific and north Fiji and Lau basin lithospheres. Since absolute plate motions reflect the sum of all forces acting on the entire lithospheric plate, a strong sea anchor mantle force may be generated at a STEP fault edge, initiating subduction (or inhibiting it), even where lithosphere is old, strong, and resists bending and without requiring large convergent forces between plates, overcoming these objections to previous models.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/gji/ggac332
