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
Direct observations of transient weakening during phase transformations in quartz and olivine
Phase transformations are widely invoked as a source of rheological weakening during subduction, continental collision, mantle convection and various other geodynamic phenomena. However, despite more than half a century of research, the likelihood and magnitude of such weakening in nature remain poorly constrained. Here we use experiments performed on a synchrotron beamline to reveal transient weakening of up to three orders of magnitude during the polymorphic quartz to coesite (SiO2) and olivine to ringwoodite (Fe2SiO4) phase transitions. Weakening becomes increasingly prominent as the transformation outpaces deformation. We suggest that this behaviour is broadly applicable among silicate minerals undergoing first-order phase transitions and examine the likelihood of weakening due to the olivine-spinel, (Mg,Fe)2SiO4, transformation during subduction. Modelling suggests that cold, wet slabs are most susceptible to transformational weakening, consistent with geophysical observations of slab stagnation in the mantle transition zone beneath the western Pacific. Our study highlights the importance of incorporating transformational weakening into geodynamic simulations and provides a quantitative basis for doing so.
more »
« less
- PAR ID:
- 10634328
- Publisher / Repository:
- Nature
- Date Published:
- Journal Name:
- Nature Geoscience
- Volume:
- 18
- Issue:
- 6
- ISSN:
- 1752-0894
- Page Range / eLocation ID:
- 548 to 554
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract One hypothesized mechanism that triggers deep‐focus earthquakes in oceanic subducting slabs below ∼300 km depth is transformational faulting due to the olivine‐to‐spinel phase transition. This study uses finite element modeling to investigate phase transformation‐induced stress redistribution and material weakening in olivine. A thermodynamically consistent constitutive model is developed to capture the evolution of phase transformation in olivine under different pressure and temperature conditions. The overall numerical model enables considering multiscale material features, including the polycrystalline structure, mesoscale heterogeneity, and various phases or variants of phases at the microscopic level, and accounts for viscoplastic behaviors with thermo‐mechanical coupling effects. The model is validated with several benchmarks, including a phase diagram of phase transformation from olivine to spinel. The validated model is used to study the interactive behaviors between defects (heterogeneity) and phase transformation. The simulation results reveal that spinel formation under pressure initiates near inclusions and along the grain boundaries, consistent with experimental observations. At lower temperatures, the transformation leads to the formation of thin conjugate bands of spinel diagonal to the compression loading direction. Local stress analysis along these bands also suggests the initiation of faulting. In contrast, the numerical results at higher transformation rates show that significant spinel formation occurs over a larger area at elevated temperatures, leading to ductile behavior, which agrees with experimental findings. Numerical simulation of multiple inclusions under confined pressure also shows the formation of a network of spinel bands resembling phase‐transformation patterns observed in the laboratory experiments. Additionally, stress softening patterns due to phase transformation are similar to experimental observations.more » « less
-
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
-
Olivine lattice preferred orientation (LPO), or texture, forms in relation to deformation mechanisms such as dislocation creep and can be observed in the upper mantle as seismic anisotropy. Olivine is also mechanically anisotropic, meaning that it responds to stresses differently depending on the direction of the stress. Understanding the interplay between anisotropic viscosity (AV) and LPO, and their role in deformation, is necessary for relating seismic anisotropy to mantle flow patterns. In this study, we employ three methods to predict olivine texture (D-Rex, MDM, and MDM+AV) in a shear box model and a subduction model. D-Rex and MDM are two representative texture development methods that have been compared before, and our results are in line with previous studies showing that textures computed by D-Rex develop faster and are stronger and more point-like than textures calculated with MDM. MDM+AV uses the same isotropic mantle stresses and particle paths as D-Rex and MDM but includes the effect of AV for texture predictions. MDM+AV predicts a texture similar to MDM with a distinct girdle-like orientation for simple shear deformation or at low strain in the subduction model. At larger strains, MDM+AV’s textures are more point-like and stronger compared to the other two methods. The effective viscosity for MDM+AV drops by up to 60% in a shear box model and can be either strengthened or weakened relative to isotropic viscosity for different regions of the subduction model experiencing different patterns of deformation. Our results emphasize the significant role of AV in olivine texture development, which could substantially affect geodynamic processes in the upper mantle.more » « less
-
Deep-focus earthquakes that occur at 350–660 km, where pressures p =12-23 GPa and temperature T =1800-2000 K, are generally assumed to be caused by olivine→spinel phase transformation, see pioneering works [1–10]. However, there are many existing puzzles: (a) What are the mechanisms for jump from geological 10−17−10−15 s−1 to seismic 10−103s−1(see [3]) strain rates? Is it possible without phase transformation? (b) How does metastable olivine, which does not completely transform to spinel at high temperature and deeply in the region of stability of spinel for over the million years, suddenly transforms during seconds and generates seismic strain rates 10−103s−1 that produce strong seismic waves? (c) How to connect deviatorically dominated seismic signals with volume-change dominated transformation strain during phase transformations [9,11]? Here we introduce a combination of several novel concepts that allow us to resolve the above puzzles quantitatively. We treat the transformation in olivine like plastic strain-induced (instead of pressure/stress-induced) and find an analytical 3D solution for coupled deformation-transformation-heating processes in a shear band. This solution predicts conditions for severe (singular) transformation-induced plasticity (TRIP) and self-blown-up deformation-transformation-heating process due to positive thermomechanochemical feedback between TRIP and strain-induced transformation. In nature, this process leads to temperature in a band exceeding the unstable stationary temperature, above which the self-blown-up shear-heating process in the shear band occurs after finishing the phase transformation. Without phase transformation and TRIP, significant temperature and strain rate increase is impossible. Due to the much smaller band thickness in the laboratory, heating within the band does not occur, and plastic flow after the transformation is very limited. Our findings change the main concepts in studying the initiation of the deep-focus earthquakes and phase transformations during plastic flow in geophysics in general. The latter may change the interpretation of different geological phenomena, e.g., the possibility of the appearance of microdiamond directly in the cold Earth crust within shear-bands [12] during tectonic activities without subduction to the mantle and uplifting. Developed theory of the self-blown-up transformation-TRIP-heating process is applicable outside geophysics for various processes in materials under pressure and shear, e.g., for new routes of material synthesis [12,13], friction and wear, surface treatment, penetration of the projectiles and meteorites, and severe plastic deformation and mechanochemical technologies.more » « less
