An official website of the United States government
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 H2O 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 H2O, on the order of 3.4–5.9 × 1021 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.
more »
« less
This content will become publicly available on June 1, 2026
Constraining Solid Dynamics, Interface Rheology, and Slab Hydration in the Hikurangi Subduction Zone Using 3‐Dimensional Fully Dynamic Models
Abstract Simulating present‐day solid Earth deformation and volatile cycling requires integrating diverse geophysical data sets and advanced numerical techniques to model complex multiphysics processes at high resolutions. Subduction zone modeling is particularly challenging due to the large geographic extent, localized deformation zones, and the strong feedbacks between reactive fluid transport and solid deformation. Here, we develop new workflows for simulating 3‐dimensional thermal‐mechanical subduction and patterns of volatile dehydration at convergent margins, adaptable to include reactive fluid transport. We apply these workflows to the Hikurangi margin, where recent geophysical investigations have offered unprecedented insight into the structure and processes coupling fluid transport and solid deformation across broad spatiotemporal scales. Geophysical data sets constraining the downgoing and overriding plate structure are collated with the Geodynamic World Builder, which provides the initial conditions for forward simulations using the open‐source geodynamic modeling software code ASPECT. We systematically examine how plate interface rheology and hydration of the downgoing plate and upper mantle influence Pacific–Australian convergence and seismic anisotropy. Models prescribing a plate boundary viscosity of Pa s best reproduce observed plate velocities, and changing the configuration of the Pacific–Australia plate boundary directly influences the modeled plate motions. Models considering hydrated olivine fabrics best reproduce observations of seismic anisotropy. Predicted patterns of slab dehydration and mantle melting correlate well with observations of seismic attenuation and arc volcanism. These results suggest that hydration‐related rheological heterogeneity and related fluid weakening may strongly influence slab dynamics. Future investigations integrating coupled fluid transport and global mantle flow will provide insight into the feedbacks between subduction dynamics, fluid pathways, and arc volcanism.
more »
« less
- Award ID(s):
- 1925677
- PAR ID:
- 10630693
- Publisher / Repository:
- Wiley
- Date Published:
- Journal Name:
- Geochemistry, Geophysics, Geosystems
- Volume:
- 26
- Issue:
- 6
- ISSN:
- 1525-2027
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Low‐angle subduction has been shown to have a profound impact on subduction processes. However, the mechanisms that initiate, drive, and sustain flat‐slab subduction are debated. Within all subduction zone systems, metamorphic dehydration reactions within the down‐going slab have been hypothesized to produce seismicity, and to produce water that fluxes melting of the asthenospheric wedge leading to arc magmatism. In this work, we examine the role hydration plays in influencing slab buoyancy and the geometry of the downgoing oceanic plate. When water is introduced to the oceanic lithosphere, it is incorporated into hydrous phases, which results in lowered rock densities. The net effect of this process is an increase in the buoyancy of the downgoing oceanic lithosphere. To better understand the role of water in low‐angle subduction settings, we model flat‐slab subduction in Alaska, where the thickened oceanic lithosphere of the Yakutat oceanic plateau is subducting beneath the continental lithosphere. In this work, we calculate the thermal conditions and stable mineral assemblages in the slab crust and mantle in order to assess the role that water plays in altering the density of the subducting slab. Our slab density results show that a moderate amount of hydration (1–1.5 wt% H2O) in the subducting crust and upper lithospheric mantle reduces slab density by 0.5%–0.8% relative to an anhydrous slab, and is sufficient to maintain slab buoyancy to 300–400 km from the trench. These models show that water is a viable factor in influencing the subduction geometry in Alaska, and is likely important globally.more » « less
-
Abstract The Alaska Peninsula has a long history of plate subduction with along‐arc variations in volcanic eruption styles and geochemistry. However, the sub‐arc melting processes that feed these volcanoes are unclear. The Alaska slab morphology below 200 km depth remains debated due to limited seismic data and thus low tomography resolution in this region. Here we utilize the newly available regional and teleseismic data to build 3‐D high‐resolutionVPandVSmodels to 660 km depth. We find that the high‐velocity Pacific Plate subducts to the bottom of the mantle transition zone (MTZ) with complex deformation and gaps. In the southwest, we observe a wide gap in the high‐velocity slab at 200–500 km depths. Toward the northeast, the slab becomes more continuous extending to the MTZ with a few holes below 200 km. We interpret these gaps as a slab tear that coincides with the subducted ancient Kula‐Pacific Ridge. We also invert for 3‐DVPandVP/VSmodels to 200 km depth with higher resolution and find strong along‐strike changes in slab dehydration and sub‐arc melting, indicated by lowVPand highVP/VSanomalies. Slab dehydration and sub‐arc melting are most extensive below the Pavlof and Shumagin segments in the southwest, becoming limited below the Chignik and Chirikof segments in the northeast, and extensive again beneath the Kodiak segment further to the northeast. We propose that the variations of slab hydration at the outer rise significantly influence slab dehydration at greater depths and further control sub‐arc melting beneath the Alaska Peninsula.more » « less
-
A typical subduction of an oceanic plate beneath a continent is expected to be accompanied by arc volcanoes along the convergent margin. However, subduction of the Cocos plate at the Middle American subduction system has resulted in an uneven distribution of magmatism/volcanism along strike. Here we construct a new three-dimensional shear-wave velocity model of the entire Middle American subduction system, using full-wave ambient noise tomography. Our model reveals significant variations of the oceanic plates along strike and down dip, in correspondence with either weakened or broken slabs after subduction. The northern and southern segments of the Cocos plate, including the Mexican flat slab subduction, are well imaged as high-velocity features, where a low-velocity mantle wedge exists and demonstrate a strong correlation with the arc volcanoes. Subduction of the central Cocos plate encounters a thick high-velocity feature beneath North America, which hinders the formation of a typical low-velocity mantle wedge and arc volcanoes. We suggest that the presence of slab tearing at both edges of the Mexican flat slab has been modifying the mantle flows, resulting in the unusual arc volcanism.more » « less
-
Abstract Whether and how subduction increases the oxidation state of Earth's mantle are two of the most important unresolved questions in solid Earth geochemistry. Using data from the southern Cascade arc (California, USA), we show quantitatively for the first time that increases in arc magma oxidation state are fundamentally linked to mass transfer of isotopically heavy sulfate from the subducted plate into the mantle wedge. We investigate multiple hypotheses related to plate dehydration and melting and the rise and reaction of slab melts with mantle peridotite in the wedge, focusing on electron balance between redox-sensitive iron and sulfur during these processes. These results show that unless slab-derived silicic melts contain much higher dissolved sulfur than is indicated by currently available experimental data, arc magma generation by mantle wedge melting must involve multiple stages of mantle metasomatism by slab-derived oxidized and sulfur-bearing hydrous components.more » « less
