The transition between blueschist and eclogite plays an important role in subduction zones via dehydration and densification processes in descending oceanic slabs. There are a number of previous petrological studies describing potential mineral reactions taking place at the transition. An experimental determination of such reactions could help constrain the pressure–temperature conditions of the transition as well as the processes of dehydration. However, previous experimental contributions have focused on the stability of spontaneously formed hydrous minerals in basaltic compositions rather than on reactions among already formed blueschist facies minerals. Therefore, this study conducted three groups of experiments to explore the metamorphic reactions among blueschist facies minerals at conditions corresponding to warm subduction, where faster reaction rates are possible on the time scale of laboratory experiments. The first group of experiments was to establish experimental reversals of the reaction glaucophane+paragonite to jadeite+pyrope+quartz+H2O over the range of 2.2–3.5 GPa and 650–820°C. This reaction has long been treated as key to the blueschist–eclogite transition. However, only the growth of glaucophane+paragonite was observed at the intersectional stability field of both paragonite and jadeite+quartz, confirming thermodynamic calculations that the reaction is not stable in the system Na2O–MgO–Al2O3–SiO2–H2O. The second set of experiments involved unreversed experiments using glaucophane+zoisite ±quartz in low‐Fe and Ca‐rich systems and were run at 1.8–2.4 GPa and 600–780°C. These produced omphacite+paragonite/kyanite+H2O accompanied by compositional shifts in the sodium amphibole, glaucophane, towards sodium–calcium amphiboles such as winchite (☐(CaNa)(Mg4Al)Si8O22(OH)2) and barroisite (☐(CaNa)(Mg3Al2)(AlSi7)O22(OH)2). This suggests that a two‐step dehydration occurs, first involving the breakdown of glaucophane+zoisite towards a paragonite‐bearing assemblage, then the breakdown of paragonite to release H2O. It also indicates that sodium–calcium amphibole can coexist with eclogite phases, thereby extending the thermal stability of amphibole to greater subduction zone depths. The third set of experiments was an experimental investigation at 2.0–2.4 GPa and 630–850°C involving a high‐Fe (Fe#=Fetotal/(Fetotal+Mg)≈0.36) natural glaucophane, synthetic paragonite and their eclogite‐forming reaction products. The results indicated that garnet and omphacite grew over most of these pressure–temperature conditions, which demonstrates the importance of Fe‐rich glaucophane in forming the key eclogite assemblage of garnet+omphacite, even under warm subduction zone conditions. Based on the experiments of this study, reaction between glaucophane+zoisite is instrumental in controlling dehydration processes at the blueschist–eclogite transition during warm subduction.
Slab surface temperature is one of the key parameters that incur first-order changes in subduction dynamics. However, the current thermal models are based on empirical thermal parameters and do not accurately capture the complex pressure–temperature paths of the subducting slab, prompting significant uncertainties on slab temperature estimations. In this study, we investigate whether the dehydration-melting of glaucophane can be used to benchmark the temperature in the slab. We observe that dehydration and melting of glaucophane occur at relatively low temperatures compared to the principal hydrous phases in the slab and produce highly conductive Na-rich melt. The electrical properties of glaucophane and its dehydration products are notably different from the hydrous minerals and silicate melts. Hence, we conclude that the thermodynamic instability of glaucophane in the slab provides a unique petrological criterion for tracking temperature in the present-day subduction systems through magnetotelluric profiles.
more » « less- Award ID(s):
- 1753125
- NSF-PAR ID:
- 10307547
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Scientific Reports
- Volume:
- 11
- Issue:
- 1
- ISSN:
- 2045-2322
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract The thermal evolution of subducting slabs controls a range of subduction processes, yet we lack a robust understanding of how thermal structure develops over a subduction zone's lifetime. We investigate the time‐dependence of slab thermal structure using dynamically consistent, time evolving models. Pressure‐temperature (
P–T ) conditions along the slab Moho and slab top exhibit substantial variability throughout the various phases of subduction: initiation, free sinking, and mature subduction. This variability occurs in response to time‐dependent subduction properties (e.g., fast vs. slow convergence) and thermal structure inherited from previous phases (e.g., due to upper plate aging). At a given depth, the slab cools rapidly during initiation, after which slower cooling occurs. In the case of the Moho, additional cooling occurs during the free sinking phase. We explore the implications of time‐dependent thermal structure on exhumed rocks and slab dehydration. Modeled slab topP–T paths span much of theP–T space associated with exhumed rocks, suggesting a significant component of recorded variability may have dynamic origins. Coupling ourP–T profiles with thermodynamic models of oceanic lithosphere, we show that dehydrating ultramafic rocks at the slab Moho provide the bulk of hydrous fluid at subarc depths during the earliest phases. Over subsequent phases, these rocks carry fluids into the deeper mantle, and it is mafic crust along the slab top that releases water at subarc depths. We conclude that varying subduction conditions, and non‐steady‐state thermal structure, challenge the utility of kinematically driven models with constant subduction parameters, particularly for investigating thermal structure in the geological past. -
more » « less
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 A key feature of subduction zone geodynamics and thermal structure is the point at which the slab and mantle mechanically couple. This point defines the depth at which traction between slab and mantle begins to drive mantle wedge circulation and also corresponds with a rapid increase in temperature along the slab‐mantle interface. Here, we consider the effects of the backarc thermal structure and slab thermal parameter on coupling depth using two‐dimensional thermomechanical models of oceanic‐continental convergent margins. Coupling depth is strongly correlated with backarc lithospheric thickness, and weakly correlated with slab thermal parameter. Slab‐mantle coupling becomes significant where weak, hydrous antigorite reacts to form strong, anhydrous olivine and pyroxene along the slab‐mantle interface. Highly efficient (predominantly advective) heat transfer in the asthenospheric mantle wedge and inefficient (predominantly conductive) heat transfer in the lithospheric mantle wedge results in competing feedbacks that stabilize the antigorite‐out reaction at depths determined primarily by the mechanical thickness of the backarc lithosphere. For subduction zone segments where backarc lithospheric thickness can be inverted from surface heat flow, our results provide a regression model that can be applied with slab thermal parameter to predict coupling depth. Consistently high backarc heat flow in circum‐Pacific subduction zones suggests uniformly thin overriding plates likely regulated by lithospheric erosion caused by hydration and melting processes under volcanic arcs. This may also explain a common depth of slab‐mantle coupling globally.
-
more » « less
Abstract The nature and cause of deep earthquakes remain enduring unknowns in the field of seismology. We present new models of thermal structures of subducted slabs traced to mantle transition zone depths that permit a detailed comparison between slab pressure/temperature (
P/T ) paths and hydrated/carbonated mineral phase relations. We find a remarkable correlation between slabs capable of transporting water to transition zone depths in dense hydrous magnesium silicates with slabs that produce seismicity below ∼300‐km depth, primarily between 500 and 700 km. This depth range also coincides with theP/T conditions at which oceanic crustal lithologies in cold slabs are predicted to intersect the carbonate‐bearing basalt solidus to produce carbonatitic melts. Both forms of fluid evolution are well represented by sublithospheric diamonds whose inclusions record the existence of melts, fluids, or supercritical liquids derived from hydrated or carbonate‐bearing slabs at depths (∼300–700 km) generally coincident with deep‐focus earthquakes. We propose that the hydrous and carbonated fluids released from subducted slabs at these depths lead to fluid‐triggered seismicity, fluid migration, diamond precipitation, and inclusion crystallization. Deep focus earthquake hypocenters could track the general region of deep fluid release, migration, and diamond formation in the mantle. The thermal modeling of slabs in the mantle and the correlation between sublithospheric diamonds, deep focus earthquakes, and slabs at depth demonstrate a deep subduction pathway to the mantle transition zone for carbon and volatiles that bypasses shallower decarbonation and dehydration processes.
