An official website of the United States government
Abstract The dip angles of slabs are among the clearest characteristics of subduction zones, but the factors that control them remain obscure. Here, slab dip angles and subduction parameters, including subduction duration, the nature of the overriding plate, slab age, and convergence rate, are determined for 153 transects along subduction zones for the present day. We present a comprehensive tabulation of subduction duration based on isotopic ages of arc initiation and stratigraphic, structural, plate tectonic and seismic indicators of subduction initiation. We present two ages for subduction zones, a long‐term age and a reinitiation age. Using cross correlation and multivariate regression, we find that (1) subduction duration is the primary parameter controlling slab dips with slabs tending to have shallower dips at subduction zones that have been in existence longer; (2) the long‐term age of subduction duration better explains variation of shallow dip than reinitiation age; (3) overriding plate nature could influence shallow dip angle, where slabs below continents tend to have shallower dips; (4) slab age contributes to slab dip, with younger slabs having steeper shallow dips; and (5) the relations between slab dip and subduction parameters are depth dependent, where the ability of subduction duration and overriding plate nature to explain observed variation decreases with depth. The analysis emphasizes the importance of subduction history and the long‐term regional state of a subduction zone in determining slab dip and is consistent with mechanical models of subduction.
more »
« less
Differentiating induced versus spontaneous subduction initiation using thermomechanical models and metamorphic soles
Abstract Despite the critical role of subduction in plate tectonics, the dynamics of its initiation remains unclear. High-temperature low-pressure metamorphic soles are vestiges of subduction initiation, providing records of the pressure and temperature conditions along the subducting slab surface during subduction initiation that can possibly differentiate the two end-member subduction initiation modes: spontaneous and induced. Here, using numerical models, we show that the slab surface temperature reaches 800–900 °C at ~1 GPa over a wide range of parameter values for spontaneous subduction initiation whereas for induced subduction initiation, such conditions can be reached only if the age of the overriding plate is <5 Ma. These modeling results indicate that spontaneous subduction initiation would be more favorable for creating high-temperature conditions. However, the synthesis of our modeling results and geological observations indicate that the majority of the metamorphic soles likely formed during induced subduction initiation that involved a young overriding plate.
more »
« less
- Award ID(s):
- 1850683
- PAR ID:
- 10281780
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Nature Communications
- Volume:
- 12
- Issue:
- 1
- ISSN:
- 2041-1723
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Abstract. Retrograde metamorphic rocks provide key insights into the pressure–temperature (P–T) evolution of exhumed material, and resultant P–T constraints have direct implications for the mechanical and thermal conditions of subduction interfaces. However, constraining P–T conditions of retrograde metamorphic rocks has historically been challenging and has resulted in debate about the conditions experienced by these rocks. In this work, we combine elastic thermobarometry with oxygen isotope thermometry to quantify the P–T evolution of retrograde metamorphic rocks of the Cycladic Blueschist Unit (CBU), an exhumed subduction complex exposed on Syros, Greece. We employ quartz-in-garnet and quartz-in-epidote barometry to constrain pressures of garnet and epidote growth near peak subduction conditions and during exhumation, respectively. Oxygen isotope thermometry of quartz and calcite within boudin necks was used to estimate temperatures during exhumation and to refine pressure estimates. Three distinct pressure groups are related to different metamorphic events and fabrics: high-pressure garnet growth at ∼1.4–1.7 GPa between 500–550 ∘C, retrograde epidote growth at ∼1.3–1.5 GPa between 400–500 ∘C, and a second stage of retrograde epidote growth at ∼1.0 GPa and 400 ∘C. These results are consistent with different stages of deformation inferred from field and microstructural observations, recording prograde subduction to blueschist–eclogite facies and subsequent retrogression under blueschist–greenschist facies conditions. Our new results indicate that the CBU experienced cooling during decompression after reaching maximum high-pressure–low-temperature conditions. These P–T conditions and structural observations are consistent with exhumation and cooling within the subduction channel in proximity to the refrigerating subducting plate, prior to Miocene core-complex formation. This study also illustrates the potential of using elastic thermobarometry in combination with structural and microstructural constraints, to better understand the P–T-deformation conditions of retrograde mineral growth in high-pressure–low-temperature (HP/LT) metamorphic terranes.more » « less
-
Abstract During the subduction of an oceanic plate, fluids are released from metabasaltic crust, metasediment, and serpentinite under high‐pressure/low‐temperature conditions. Although some fluids may eventually leave the slab, some participate in metamorphic reactions within the slab during subduction and exhumation. To identify fluid sources and other controls influencing mineral composition, we report the in situ‐measured δ18O of lawsonite and garnet in blueschist‐ to eclogite‐facies rocks from 10 subduction zones that represent various field settings, including mélanges, structurally coherent terranes, and an eclogite xenolith derived from a subducted plate. Lawsonite records distinct δ18O depending on the host rock type and other rock types that were fluid sources during lawsonite growth. In general, lawsonite in metabasalt (7.6 ± 0.2–14.8 ± 1.1‰) is isotopically lighter than in metasediment (20.6 ± 1.4–24.1 ± 1.3‰) but heavier than in metagabbro (4.0 ± 0.4–7.9 ± 0.3‰). The extent of δ18O fractionation was evaluated for lawsonite–fluid and lawsonite–garnet pairs as a function of temperature (T). Results demonstrate that variations of >1.7‰ in lawsonite and >0.9‰ in garnet are not related to changingT. More likely, the relative contributions of fluids derived from isotopically heavier lithologies (e.g., sediments) versus lighter lithologies (e.g., ultramafic rocks) are the major control. Monte Carlo simulations were performed to investigate the sources of metasomatic fluids and the water/rock ratio that formed lawsonite‐bearing metasomatite. Results indicate that δ18OLwsand δ18OGrtrecord interactions with fluids sourced from diverse lithologies (sediment, serpentinite), further supporting that δ18OLwsis a useful indicator of subduction fluid‐rock interactions.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
