Search for: All records

Abstract During subduction, the downgoing oceanic crust is exposed to high temperatures in the mantle wedge, causing volatile‐bearing minerals to break down and release hydrous fluids into the forearc. These fluids percolate upwards, reacting with the mantle wedge to form hydrated ultramafic lithologies, including serpentinite. To accurately track the fate and impact of water on the forearc, we develop time‐dependent models that self‐consistently capture both serpentinite ingrowth and the associated rheological weakening of the plate interface. Unlike many subduction models that investigate forearc serpentinization and prescribe plate velocities, geometries, or steady‐state conditions, our approach allows plates to evolve dynamically without predefined velocities or geometries. During subduction infancy, serpentinite accumulates rapidly. As subduction matures, serpentinite ingrowth decreases, and more serpentinite is also dragged downward by the slab. To elucidate the links between subduction dynamics and serpentinization, we consider variations in serpentinite strength and hydration state of the incoming plate. Subducting fully water‐saturated sediments yield ∼3× greater forearc serpentinite than within the moderately hydrated reference case. The water‐saturated case produces a weaker interface and, in turn, subduction zone convergence rates ∼40% higher than in an endmember case with anhydrous sediment. A lower serpentinite strength also produces higher convergence rates despite more downdragging of serpentinite from the forearc. We find that hydrous sediments not only lubricate the interface directly by weakening it, as previously suggested, but also by dehydrating and releasing water that produces weak serpentinite in the mantle wedge, with such feedback only able to be captured within fully coupled dynamic models. 

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. 

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.