An official website of the United States government
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
Evolving Subduction Zone Thermal Structure Drives Extensive Forearc Mantle Wedge Hydration
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
- PAR ID:
- 10526727
- Publisher / Repository:
- American Geophysical Union
- Date Published:
- Journal Name:
- AGU Advances
- Volume:
- 5
- Issue:
- 4
- ISSN:
- 2576-604X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
In the current plate tectonic regime, thermal modeling, petrology, and seismology show that subsurface portions of cold slabs carry some of their volatiles into the deep upper mantle, mantle transition zone, and uppermost lower mantle avoiding the devolatilization occurring with normal arc and wedge subduction. Slab crustal remnants at these depths can melt by intersecting their carbonated solidus whereas slab mantle remnants can devolatilize by warming and metamorphosing to ‘dryer’ mineral assemblages. Since fluid release and earthquake production (“dehydration embrittlement”) operates down to ~300 km depths in all subduction zones, we propose, that deep-focus earthquakes trace the places of fluid release at deeper levels (350 to 750 km). Fluids in faults related to earthquake generation will become diamond-forming as they react with mantle rocks along the fault walls. Diamonds thus formed will record deformation produced by mantle convection and slab buckling during mantle storage. Lithospheric diamonds, stored in static ancient continental keels, lack the connection to this type of geodynamic regime that is evident for sublithospheric diamonds. However, a comparison between the two diamond types suggests a geologic model for lithospheric diamond formation in the ancient past. Lithospheric diamonds and sublithospheric diamonds both contain evidence for the recycling of sediments or surficial rocks that have equilibrated at low temperatures with seawater. The known way to inject these materials into diamond-forming regions is slab subduction. Hence both diamond types may have formed by variants of this same process that differ in depth and style over geologic time. Lithospheric diamonds are different from sublithospheric diamonds in critical ways: higher average N content, ages extending into the Paleoarchean, inclusion assemblages indicating formation at lower pressure, and lack of ubiquitous deformation features. Nitrogen content is the key to relating lithospheric diamonds to the subducting slab. Nitrogen occurs in clays and sediments at the slab surface or uppermost crust. Regardless of whether the slab is hot or cold during subduction, nitrogen will be removed into a mantle wedge if one exists. Additionally, diamonds will not survive in the melts/fluids generated in the wedge under oxidizing conditions. For sublithospheric diamonds, their low to non-existent nitrogen content occurs because they are derived from slab fluids/melts once nitrogen has been largely removed or from rocks deeper in the slab where nitrogen is scarce. The much higher nitrogen in lithospheric diamonds suggests that they formed from fluids/melts derived near the slab surface that contained N. In the Archean, such slabs must have subducted close to the nascent mantle keel with no mantle wedge so the fluids could be directly reduced by the mantle keel. We propose a gradual temporal change from shallow, keel-adjacent, mantle-wedge-poor subduction that produced lithospheric diamonds starting in the Paleoarchean to wedge-avoiding, cold and deep subduction that produced sublithospheric diamonds in the Paleozoic. This temporal change is consistent with many geologic features: an early stagnant lid and a buoyant Archean oceanic lithosphere; the slab-imbrication, advective thickening, and diamond-richness of portions of mantle keels; and anomalously diamond-rich ancient eclogites.more » « less
-
Abstract The dehydration and decarbonation in the subducting slab are intricately related and the knowledge of the physical properties of the resulting C–H–O fluid is crucial to interpret the petrological, geochemical, and geophysical processes associated with subduction zones. In this study, we investigate the C–H–O fluid released during the progressive devolatilization of carbonate-bearing serpentine-polymorph chrysotile, with in situ electrical conductivity measurements at high pressures and temperatures. The C–H–O fluid produced by carbonated chrysotile exhibits high electrical conductivity compared to carbon-free aqueous fluids and can be an excellent indicator of the migration of carbon in subduction zones. The crystallization of diamond and graphite indicates that the oxidized C–H–O fluids are responsible for the recycling of carbon in the wedge mantle. The carbonate and chrysotile bearing assemblages stabilize dolomite during the devolatilization process. This unique dolomite forming mechanism in chrysotile in subduction slabs may facilitate the transport of carbon into the deep mantle.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
-
Abstract Sulfur belongs among H2O, CO2, and Cl as one of the key volatiles in Earth’s chemical cycles. High oxygen fugacity, sulfur concentration, and δ34S values in volcanic arc rocks have been attributed to significant sulfate addition by slab fluids. However, sulfur speciation, flux, and isotope composition in slab-dehydrated fluids remain unclear. Here, we use high-pressure rocks and enclosed veins to provide direct constraints on subduction zone sulfur recycling for a typical oceanic lithosphere. Textural and thermodynamic evidence indicates the predominance of reduced sulfur species in slab fluids; those derived from metasediments, altered oceanic crust, and serpentinite have δ34S values of approximately −8‰, −1‰, and +8‰, respectively. Mass-balance calculations demonstrate that 6.4% (up to 20% maximum) of total subducted sulfur is released between 30–230 km depth, and the predominant sulfur loss takes place at 70–100 km with a net δ34S composition of −2.5 ± 3‰. We conclude that modest slab-to-wedge sulfur transport occurs, but that slab-derived fluids provide negligible sulfate to oxidize the sub-arc mantle and cannot deliver34S-enriched sulfur to produce the positive δ34S signature in arc settings. Most sulfur has negative δ34S and is subducted into the deep mantle, which could cause a long-term increase in the δ34S of Earth surface reservoirs.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1029/2023AV001121
