Search for: All records

Hydration of the subduction zone forearc mantle wedge influences the downdip distribution of seismicity, the availability of fluids for arc magmatism, and Earth's long term water cycle. Reconstructions of present‐day subduction zone thermal structures using time‐invariant geodynamic models indicate relatively minor hydration, in contrast to many geophysical and geologic observations. We pair a dynamic, time‐evolving thermal model of subduction with phase equilibria modeling to investigate how variations in slab and forearc temperatures from subduction infancy through to maturity contribute to mantle wedge hydration. We find that thermal state during the intermediate period of subduction, as the slab freely descends through the upper mantle, promotes extensive forearc wedge hydration. In contrast, during early subduction the forearc is too hot to stabilize hydrous minerals in the mantle wedge, while during mature subduction, slab dehydration dominantly occurs beyond forearc depths. In our models, maximum wedge hydration during the intermediate phase is 60%–70% and falls to 20%–40% as quasi‐steady state conditions are approached during maturity. Comparison to global forearc H₂O capacities reveals that consideration of thermal evolution leads to an order of magnitude increase in estimates for current extents of wedge hydration and provides better agreement with geophysical observations. This suggests that hydration of the forearc mantle wedge represents a potential vast reservoir of H₂O, on the order of 3.4–5.9 × 10²¹ g globally. These results provide novel insights into the subduction zone water cycle, new constraints on the mantle wedge as a fluid reservoir and are useful to better understand geologic processes at plate margins. 

Talc‐rich rocks are common in exhumed subduction zone terranes and may explain geophysical observations of the subduction zone interface, particularly beneath Guerrero, Mexico, where the Cocos plate subducts horizontally beneath North America and episodic tremor and slow slip (ETS) occurs. We present petrologic models exploring (a) the degree of silica metasomatism required to produce talc in serpentinized peridotites at the pressure‐temperature conditions of the plate interface beneath Guerrero and (b) the amount of silica‐bearing water produced by rocks from the subducting Cocos plate and the location of fluid pulses. We estimate the volumes of talc produced by the advection of silica‐rich fluids into serpentinized peridotites at the plate interface over the history of the flat‐slab system. In the ETS‐hosting region, serpentinites must achieve ∼43 wt. % SiO₂to stabilize talc, but minor additions of silica beyond this produce large volumes of talc. Our models of Cocos plate dehydration predict that water flux into the interface averages 3.9 × 10⁴ kg m⁻² Myr⁻¹but suggest that only where subducting basalts undergo major dehydration reactions will sufficient amounts of silica‐rich fluids be produced to drive significant metasomatism. We suggest that talc produced by advective transport of aqueous silica alone cannot account for geophysical interpretations of km‐thick zones of talc‐rich rocks beneath Guerrero, although silica‐bearing fluids that migrate along the plate interface may promote broader metasomatism. Regions of predicted talc production do, however, overlap with the spatial occurrence of ETS, consistent with models of slow slip based on the frictional deformation of metasomatic lithologies. 

SUMMARY Tectonic plate motions predominantly result from a balance between the potential energy change of the subducting slab and viscous dissipation in the mantle, bending lithosphere and slab–upper plate interface. A wide range of observations from active subduction zones and exhumed rocks suggest that subduction interface shear zone rheology is sensitive to the composition of subducting crustal material—for example, sediments versus mafic igneous oceanic crust. Here we use 2-D numerical models of dynamically consistent subduction to systematically investigate how subduction interface viscosity influences large-scale subduction kinematics and dynamics. Our model consists of an oceanic slab subducting beneath an overriding continental plate. The slab includes an oceanic crustal/weak layer that controls the rheology of the interface. We implement a range of slab and interface strengths and explore how the kinematics respond for an initial upper mantle slab stage, and subsequent quasi-steady-state ponding near a viscosity jump at the 660-km-discontinuity. If material properties are suitably averaged, our results confirm the effect of interface strength on plate motions as based on simplified viscous dissipation analysis: a ∼2 order of magnitude increase in interface viscosity can decrease convergence speeds by ∼1 order of magnitude. However, the full dynamic solutions show a range of interesting behaviour including an interplay between interface strength and overriding plate topography and an end-member weak interface-weak slab case that results in slab break-off/tearing. Additionally, for models with a spatially limited, weak sediment strip embedded in regular interface material, as might be expected for the subduction of different types of oceanic materials through Earth’s history, the transient response of enhanced rollback and subduction velocity is different for strong and weak slabs. Our work substantiates earlier suggestions as to the importance of the plate interface, and expands the range of quantifiable links between plate reorganizations, the nature of the incoming and overriding plate and the potential geological record.