Juxtaposed slab dehydration, decarbonation and seismotectonic variation beneath the Philippine subduction zone based on 3-D modeling

Largescale volcanic eruptions and earthquakes are occurring frequently in the Philippines, and research has shown that slab metamorphism and diversity alter the impacts of subducted oceanic plates by changing water‒carbon productivity and interplate stability. Within the framework of the thermal evolution history of subducting slabs, the relationships between subduction zone seismicity characterized by both regular megathrust earthquakes and slow slip events of various magnitudes and long-term slab dehydration-decarbonation evolution in the Philippines remain poorly understood. Here, we constructed a comprehensive thermal model incorporating 3-D slab geometric data for the incoming plate and a 3-D subduction velocity field based on the MORVEL plate motion dataset for the Philippine subduction zone with high spatial and temporal resolutions. Our findings reveal that subduction seismicity and arc volcanism are prominent in belt-shaped regions with high thermal gradients (> 5 °C/km) and large-scale slab dehydration (> 0.05 wt%/km). Dehydration of serpentinite in ultramafic rocks in the subducting slab and decarbonation of carbonate minerals preferentially contribute to the generation and transport of fluids and carbonate melts, thus facilitating seismicity and carbon-rich magmatism. Our results suggest that slab geometry diversity-induced juxtaposed slab dehydration-decarbonation processes play a vital role in the generation of megathrust earthquakes below the forearc.

Numerical modeling the process of deep slab dehydration and magmatism

This study uses a 2D high-resolution thermo-mechanical coupled model to investigate the dynamic processes of deep plate hydration, dehydration, and subsequent magmatic activity in ocean-continent subduction zones. We reveal the pathways and temporal evolution of water transport to the deep mantle during the subduction process. Plate dehydration plays a critical role in triggering partial melting of the deep mantle and related magmatic activity. Our study shows significant differences in the volumes of melt produced at different depths, with dehydration reactions in deeper regions being weaker compared to shallower ones. It takes a longer time to reach the suitable P-T conditions for hydrous melting in the deep mantle. The results highlight the geophysical significance of water transport in deep subduction zones and its role in magmatic processes, particularly in the formation of magma chambers beneath continental plates.

Control of slab tears and slab flat wedging on volcanism in the Alaska subduction zone

Multistage plate subduction plays a crucial role in magmatism; however, the mechanisms by which deep geodynamic processes govern volcanism in the Alaska subduction zone remain controversial. Using numerous travel-time data from several seismic arrays, we constructed high-resolution tomographic models to investigate the velocity structure of the Pacific Plate and Yakutat slab. Our tomographic results revealed high-velocity anomalies in the Pacific Plate and Yakutat slab, while the low-velocity areas within the Pacific Plate were identified as slab tears. We suggest that the Pacific Plate transitioned from oblique subduction along the Aleutian volcano chain to lower-angle subduction beneath the Pacific-Yakutat Plate interaction zone, forming two slab tears that enhance hot asthenosphere materials upwelling. The partial melting of the mantle wedge induced by Pacific slab dehydration and the concurrent upwelling of mantle materials jointly drove volcanism in the transition zone. Conversely, the flat subduction of the Yakutat slab into the mantle wedge overlying the Pacific slab effectively hindered the upwelling of hot hybrid materials, cooling the Pacific mantle wedge. These results offer a new perspective on the influence of slab dynamics on volcanic and magmatic processes in the region and represent an advancement in our understanding compared to previous studies, which did not resolve the tears within the slab or their geodynamic implications at this level of detail.

Light δD apatites reveal deep origin water in North China Craton intracontinental granites and basalts

Water is essential to the formation of intracontinental granites, but its origin remains elusive. Here we address this scientific problem by analyzing D/H isotopes of apatites, hydrous minerals in Jurassic and Early Cretaceous granites and basalts from eastern North China Craton, where water was previously interpreted as derived from subducting slab. Results reveal extremely low δD values in pristine Early Cretaceous granitic (-203‰ to -127‰) and basaltic (-197‰ to -107‰) apatites, contrasting with relatively high δD values (-137‰ to -47‰) in Jurassic granites. Given the depth-dependent D/H isotopic fractionation during slab dehydration and high-water contents in coeval primitive mafic magmas, the Early Cretaceous magma water is attributed to the stagnant slab within the mantle transition zone. Secular change in the depth of water aligns with steepening of subducting Paleo-Pacific plate from Jurassic to Early Cretaceous, demonstrating the potential of apatite H isotopes in tracing water origin in granites and basalts.

Superionic iron hydride shapes ultralow-velocity zones at Earth's core-mantle boundary

Seismological studies have exposed numerous ultralow velocity zones (ULVZs) exhibiting extraordinary physical attributes at Earth's core-mantle boundary, yet their compositions and origins remain controversial. Water-iron reaction can generate unique phases under lowermost-mantle conditions and likely plays a crucial role in forming ULVZs. Through first-principles molecular dynamic simulations with machine learning techniques, we determine that iron hydride, the product of water-iron reaction, is stable as a superionic phase at the core-mantle boundary. This superionic iron hydride has much slower velocities and a higher density than the ambient mantle under lowermost-mantle conditions. Accumulation of iron hydride, created through either a chemical reaction between subducted water and iron or solidification of core material entrained in the lower mantle by convection, can explain the seismic observations of ULVZs particularly those associated with subduction. This work suggests that water may have a substantial role in creating seismic heterogeneities at the core-mantle boundary.

A role for subducting clays in the water transportation into the Earth's lower mantle

Subducting sedimentary layer typically contains water and hydrated clay minerals. The stability of clay minerals under such hydrous subduction environment would therefore constraint the lithology and physical properties of the subducting slab interface. Here we show that pyrophyllite (Al2Si4O10(OH)2), one of the representative clay minerals in the alumina-silica-water (Al2O3-SiO2-H2O, ASH) system, breakdowns to contain further hydrated minerals, gibbsite (Al(OH)3) and diaspore (AlO(OH)), when subducts along a water-saturated cold subduction geotherm. Such a hydration breakdown occurs at a depth of ~135 km to uptake water by ~1.8 wt%. Subsequently, dehydration breakdown occurs at ~185 km depth to release back the same amount of water, after which the net crystalline water content is preserved down to ~660 km depth, delivering a net amount of ~5.0 wt% H2O in a phase assemblage containing δ-AlOOH and phase Egg (AlSiO3(OH)). Our results thus demonstrate the importance of subducting clays to account the delivery of ~22% of water down to the lower mantle.

Early arc crust formation preserved in the Grenadines archipelago, southern Lesser Antilles arc

Intra-arc diversity in volcanic activity and composition is ubiquitous, but its underlying causes remain largely unresolved in many settings. In this work, we examine such variability in the Grenadines archipelago, southern Lesser Antilles arc. Here, juxtaposed volcanic centres exhibit eruptive longevities and chemistries distinct from northern counterparts in the same arc. Our goal is to explain this deviation by investigating variations in magmatic processes using petrological data from erupted crustal xenoliths and lavas, and interpreting these findings within the context of the archipelago's tectonic history and geophysical structure. Textural analyses of xenoliths reveal crystallization over a wide range of pressure-temperature-melt composition conditions in the crust. Mineral phases display discrete compositional trends pointing towards significant inter-island variability in underlying plumbing systems. The geochemical variety of erupted magmas is reminiscent of the entire arc. We speculate that the Grenadines represents the early onset of subduction forming the modern-day Lesser Antilles arc. Extrusive volcanism initiated as submarine activity. Subsequent uplift eroded the original topography of these volcanic centres following the eventual cessation of volcanism in the Neogene. The positioning of the Grenadines on an elevated platform provides rare modern insight into early arc crust formation not commonly preserved in established active arcs.

Sulfur disproportionation in deep COHS slab fluids drives mantle wedge oxidation

Sulfur degassed at volcanic arcs calls for dissolved sulfate ions (S6+) released by subduction-zone fluids, oxidizing (in association with carbon) the subarc mantle, but sulfur speciation in subduction fluids at subarc depths remains unclear. We apply electrolytic fluid thermodynamics to model the dissolution behavior of pyrite in metacarbonate sediments as a function of P, T and rock redox state up to 4.3 gigapascals and 730°C. At subarc depth and the redox conditions of the fayalite-magnetite-quartz oxygen buffer, pyrite dissolution releases oxidized sulfur in fluids by disproportionation into sulfate, bisulfite, and sulfide species. These findings indicate that oxidized, sulfur-rich carbon-oxygen-hydrogen-sulfur (COHS) fluids form within subducting slabs at depths greater than 100 kilometers independent from slab redox state and that sulfur can be more effective than the concomitantly dissolved carbon at oxidizing the mantle wedge, especially when carbonates are stable.

Understanding subduction infancy to mature subduction in Southwest Japan via the self-consistent formation of a weak slab interface

The weak slab interface controls long-term subduction dynamics. A weak hydrous layer at the slab interface promotes mechanical decoupling between the forearc mantle and the subducting slab and converts a hot forearc mantle to a cold mantle. Often referred to as a cold nose, the cold forearc mantle, plays a key role in the transition from subduction infancy to mature subduction. This study was the first to numerically demonstrate the self-consistent formation of a weak hydrous layer with permeability anisotropy based on the Southwest Japan subduction zone case, where transition-related geological features were present. Our models showed that mechanical decoupling by spontaneous downdip growth of the weak hydrous layer created a cold nose by converting a hot forearc mantle to a cold mantle. The emergence of a cold nose explained the migration of the forearc-to-arc volcanic front, expressed as the formation of mid-Miocene forearc high-magnesium andesite and Quaternary arc adakite. Furthermore, the weak hydrous layer providing a pathway for free-water transport toward the mantle wedge tip elucidates slab/mantle-derived geochemical components in deep groundwater as well as large S-wave delay times and non-volcanic seismic tremors in the forearc.

Tracing fluid infiltration into oceanic crust up to ultra-high-pressure conditions

Fluid-rock interaction within the altered oceanic crust and across the slab-mantle boundary during subduction facilitates element transfer, but the dynamics of fluid transport and fluid-rock exchange during upward fluid migration are still unclear. A study of metamorphic fluid-rock interaction within a section of subducted oceanic crust was carried out on eclogites and metasediments of the ultra-high-pressure Lago di Cignana Unit (NW Italian Alps). The P-T modeling of a quartzschist shows that garnet grew during the prograde and sporadically during the retrograde path and that phengite mainly records the peak to retrograde conditions. Microscale geochemical analysis of garnets has revealed a systematic evolution of oxygen isotopic composition with garnet major element zonation, with extreme within-sample core-rim variations in δ18O between 18 and 4‰ providing evidence for external fluid influx. Garnet in eclogites and calcschists, as well as garnet cores in quartz-rich lithologies, shows normal compositional zoning, as expected for prograde garnet growth, and a relatively constant oxygen isotopic composition. The outer garnet growth zones within a few metasediments show reverse or discontinuous zoning and progressively lower δ18O. Despite major element zoning, the isotopic composition of mica is homogeneous across chemical zoning in one eclogite and one quartzschist, but shows 6‰ variability in another quartzschist. In the underlying Zermatt-Saas serpentinites, antigorite from nine serpentinite samples shows some variation in δ18O, with average δ18O values for individual samples ranging from 1 to 6‰. These results provide evidence for two main stages of external fluid infiltration: (i) fluids from the dehydration of mafic lithologies entered the sequence at peak conditions around 3 GPa, as indicated by the oxygen composition of intermediate zones of mica and garnet, and (ii) low δ18O fluids from serpentinites infiltrated parts of the sedimentary package during exhumation prior to 1.5 GPa, as recorded by the 4‰ garnet outer rims. Samples recording external fluid infiltration are concentrated in the lower part of the sequence, indicating channelized fluid flow, suggesting focused fluid infiltration due to permeability contrasts between metasedimentary and eclogitic lithologies. Channelized fluid flow in the ultra-high-pressure metasediments of Lago di Cignana has not resulted in systematic decarbonation of the metasediments.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00410-023-02060-6.

Role of subduction dynamics on the unevenly distributed volcanism at the Middle American subduction system

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.

Silica-water superstructure and one-dimensional superionic conduit in Earth's mantle

Water in Earth's deep interior is predicted to be hydroxyl (OH-) stored in nominally anhydrous minerals, profoundly modulating both structure and dynamics of Earth's mantle. Here, we use a high-dimensional neuro-network potential and machine learning algorithm to investigate the weight percent water incorporation in stishovite, a main constituent of the subducted oceanic crust. We found that stishovite and water prefer forming medium- to long-range ordered superstructures, featuring one-dimensional (1D) water channels. Synthesizing single crystals of hydrous stishovite, we verified the ordering of OH- groups in the water channels through optical and nuclear magnetic resonance spectroscopy and found an average H-H distance of 2.05(3) Å, confirming simulation results. Upon heating, H atoms were predicted to behave fluid-like inside the channels, leading to an exotic 1D superionic state. Water-bearing stishovite could feature high ionic mobility and strong electrical anisotropy, manifesting as electrical heterogeneity in Earth's mantle.

Reestimation of slab dehydration fronts in Kuril-Kamchatka using updated global subduction zone thermal structures

Previous subduction thermal models are inconsistent with the values of forearc heat flow (50-140 mW/m²) and global P‒T conditions of exhumed rocks, both suggesting a shallow environment 200-300°C warmer than model predictions. Here, we revaluate these problems in Kuril-Kamchatka using 3D thermomechanical modeling that satisfies the observed subduction history and slab geometry, while our refined 3D slab thermal state is warmer than that predicted by previous 2D models and better matches observational constraints. We show that warmer slabs create hierarchical slab dehydration fronts at various forearc depths, causing fast and slow subduction earthquakes. We conclude that fast-to-slow subduction earthquakes all play a key role in balancing plate coupling energy release on megathrusts trenchward of high P-T volcanism.

Subducting volcaniclastic-rich upper crust supplies fluids for shallow megathrust and slow slip

Recurring slow slip along near-trench megathrust faults occurs at many subduction zones, but for unknown reasons, this process is not universal. Fluid overpressures are implicated in encouraging slow slip; however, links between slow slip, fluid content, and hydrogeology remain poorly known in natural systems. Three-dimensional seismic imaging and ocean drilling at the Hikurangi margin reveal a widespread and previously unknown fluid reservoir within the extensively hydrated (up to 47 vol % H2O) volcanic upper crust of the subducting Hikurangi Plateau large igneous province. This ~1.5 km thick volcaniclastic upper crust readily dewaters with subduction but retains half of its fluid content upon reaching regions with well-characterized slow slip. We suggest that volcaniclastic-rich upper crust at volcanic plateaus and seamounts is a major source of water that contributes to the fluid budget in subduction zones and may drive fluid overpressures along the megathrust that give rise to frequent shallow slow slip.

Tracking carbon from subduction to outgassing along the Aleutian-Alaska Volcanic Arc

Subduction transports volatiles between Earth's mantle, crust, and atmosphere, ultimately creating a habitable Earth. We use isotopes to track carbon from subduction to outgassing along the Aleutian-Alaska Arc. We find substantial along-strike variations in the isotopic composition of volcanic gases, explained by different recycling efficiencies of subducting carbon to the atmosphere via arc volcanism and modulated by subduction character. Fast and cool subduction facilitates recycling of ~43 to 61% sediment-derived organic carbon to the atmosphere through degassing of central Aleutian volcanoes, while slow and warm subduction favors forearc sediment removal, leading to recycling of ~6 to 9% altered oceanic crust carbon to the atmosphere through degassing of western Aleutian volcanoes. These results indicate that less carbon is returned to the deep mantle than previously thought and that subducting organic carbon is not a reliable atmospheric carbon sink over subduction time scales.

Water-sulfur-rich, oxidised adakite magmas are likely porphyry copper progenitors

The world's largest current Cu resource is volcanic arc-hosted, porphyry copper deposits. Whether unusual parental magmas or fortuitous combinations of processes accompanying emplacement of common parental arc magmas (e.g., basalt) is required for ore deposit formation, remains unclear. Spatial and tectonic associations of adakite (andesite with high La/Yb, Sr/Y) with porphyries exist, but genetic links are debated. Delayed saturation with Cu-bearing sulfides consequent to elevated redox state seems essential for late-stage exsolution of Cu-bearing hydrothermal fluids. Partial melting of igneous layers of subducted, hydrothermally altered oceanic crust in the eclogite stability field are invoked to account for andesitic compositions, residual garnet signatures, and the putative oxidised character of adakites. Alternative petrogeneses include partial melting of lower crustal, garnet-bearing sources and extensive intra-crustal amphibole fractionation. Here we demonstrate mineral-hosted, adakite glass (formerly melt) inclusions in lavas erupted subaqueously in the New Hebrides arc are oxidised relative to island arc (and mid-ocean ridge) basalts, are H₂O-S-Cl-rich, and moderately enriched in Cu. Polynomial fitting of chondrite-normalised, rare earth element abundance patterns shows the precursors of these erupted adakites were unequivocally derived from partial melting of subducted slab, and represent optimal porphyry copper progenitors.

Subduction thermal regime, petrological metamorphism and seismicity under the Mariana arc

Because of the steep subduction of a highly concave slab, researchers have characterized megathrusts under the Marianas as among the coldest and curviest plate coupling interfaces in various circum-Pacific subduction zones. Seismic tomography indicates that the heterogeneous underlying plate varies markedly in its subduction angle, velocity, and flexure along the strike and dip, while their effects on the thermal structure and intraslab earthquake occurrence remain enigmatic. By incorporating the 3-D MORVEL velocity and state-of-the-art slab geometry into thermomechanical modeling, we estimated the 3-D subduction thermal state and hydrothermal regime below the Marianas. We find that (1) the concave slab geometry and the complexity of the intraslab velocity variation in the Marianas are associated with a heterogeneous along-strike thermal regime and a cold mantle wedge beneath the central Marianas; (2) amphibolitization and eclogitization of subducted oceanic crust cause variations in fluid pressure and fluid release from the subduction interface, which may influence the distribution of interface seismicity in the Mariana system; (3) the concentration of active hydrothermal vents in the trench > 8 km deep is accompanied by a large temperature gradient and subsequent remarkable slab dehydration in the southern Marianas; and (4) slab dehydration (> 0.02 wt%/km) from 30 to 80 km indicates notable fluid release and potential fluid migration in subduction channels, which may correspond to the large water flux at depth beneath the Marianas.

Massive abiotic methane production in eclogite during cold subduction

Methane (CH₄) is a critical but overlooked component in the study of the deep carbon cycle. Abiotic CH₄ produced by serpentinization of ultramafic rocks has received extensive attention, but its formation and flux in mafic rocks during subduction remain poorly understood. Here, we report massive CH₄-rich fluid inclusions in well-zoned garnet from eclogites in Western Tianshan, China. Petrological characteristics and carbon-hydrogen isotopic compositions confirm the abiotic origin of this CH₄. Reconstructed P-T-fO₂-fluid trajectories and Deep Earth Water modeling imply that massive abiotic CH₄ was generated during cold subduction at depths of 50-120 km, whereas CO₂ was produced during exhumation. The massive production of abiotic CH₄ in eclogites may result from multiple mechanisms during prograde high pressure-ultrahigh pressure metamorphism. Our flux calculation proposes that abiotic CH₄ that has been formed in HP-UHP eclogites in cold subduction zones may represent one of the largest, yet overlooked, sources of abiotic CH₄ on Earth.

Elastic properties of Fe-bearing Akimotoite at mantle conditions: Implications for composition and temperature in lower mantle transition zone

The pyrolite model, which can reproduce the upper-mantle seismic velocity and density profiles, was suggested to have significantly lower velocities and density than seismic models in the lower mantle transition zone (MTZ). This argument has been taken as mineral-physics evidence for a compositionally distinct lower MTZ. However, previous studies only estimated the pyrolite velocities and density along a one-dimension (1D) geotherm and never considered the effect of lateral temperature heterogeneity. Because the majorite-perovskite-akimotoite triple point is close to the normal mantle geotherm in the lower MTZ, the lateral low-temperature anomaly can result in the presence of a significant fraction of akimotoite in pyrolitic lower MTZ. In this study, we reported the elastic properties of Fe-bearing akimotoite based on first-principles calculations. Combining with literature data, we found that the seismic velocities and density of the pyrolite model can match well those in the lower MTZ when the lateral temperature heterogeneity is modeled by a Gaussian distribution with a standard deviation of ∼100 K and an average temperature of dozens of K higher than the triple point of MgSiO3. We suggest that a harzburgite-rich lower MTZ is not required and the whole mantle convection is expected to be more favorable globally.

Relationship between tectonic tremors and 3-D distributions of thermal structure and dehydration in the Alaska subduction zone

The Alaska subduction zone is characterized by a subducting oceanic plateau, which is referred to as the Yakutat terrane. Tectonic tremors occur in this zone, and there are few volcanoes above the subducted Yakutat terrane. In this study, we performed a 3-D numerical simulation of a thermal structure associated with the simultaneous subduction of the Yakutat terrane and Pacific plate to elucidate the mechanism of tectonic tremors, which typically involve the presence of water. We calculated the water content distribution near the slab surface by using the thermal structure obtained from our simulation and phase diagrams of the hydrous minerals included in the slab. As a result, dehydration from the marine sedimentary layer and oceanic crust was observed near the area where tectonic tremors occurred. Tectonic tremors occur only in the Yakutat terrane because the marine sedimentary layer and oceanic crust are thicker there, and the total amount of water content in these layers is higher; therefore, the amount of dehydration is also higher there than in the Pacific plate. Additionally, there are few volcanoes above the subducted Yakutat terrane because little water remains within the slab beneath the volcanic chain where magma is produced.

Melting of subducted slab dictates trace element recycling in global arcs

Arc magma acquires continental crust-like trace element signatures through selective recycling of incompatible elements from the subducted slab. The long-standing model of element recycling through aqueous fluid from altered oceanic crust (AOC) and sediment melt has been challenged by the resurgence of mélange diapir (a mix of AOC, sediment, and serpentinite) and saline aqueous fluid models. Here, we present experimental data for near-solidus sediment melts and a framework for calculating trace element concentrations in subduction fluids from metamorphosed sediment and oceanic crust. We observe that variation of element ratios in global primitive arc basalts is comparable with that of sediment and/or oceanic crustal melt, rather than (saline) aqueous fluid or mélange melt. In particular, the systematic correlation of element ratios in arc basalt corresponds to element fractionation in slab melt with temperature and therefore follows a power function. Our findings suggest that slab melt is primarily responsible for element recycling to the arc.

Lithium, Oxygen and Magnesium Isotope Systematics of Volcanic Rocks in the Okinawa Trough: Implications for Plate Subduction Studies

Determining the influence of subduction input on back-arc basin magmatism is important for understanding material transfer and circulation in subduction zones. Although the mantle source of Okinawa Trough (OT) magmas is widely accepted to be modified by subducted components, the role of slab-derived fluids is poorly defined. Here, major element, trace element, and Li, O and Mg isotopic compositions of volcanic lavas from the middle OT (MOT) and southern OT (SOT) were analyzed. Compared with the MOT volcanic lavas, the T9-1 basaltic andesite from the SOT exhibited positive Pb anomalies, significantly lower Nd/Pb and Ce/Pb ratios, and higher Ba/La ratios, indicating that subducted sedimentary components affected SOT magma compositions. The δ7Li, δ18O, and δ26Mg values of the SOT basaltic andesite (−5.05‰ to 4.98‰, 4.83‰ to 5.80‰ and −0.16‰ to −0.09‰, respectively) differed from those of MOT volcanic lavas. Hence, the effect of the Philippine Sea Plate subduction component, (low δ7Li and δ18O and high δ26Mg) on magmas in the SOT was clearer than that in the MOT. This contrast likely appears because the amounts of fluids and/or melts derived from altered oceanic crust (AOC, lower δ18O) and/or subducted sediment (lower δ7Li, higher δ18O and δ26Mg) injected into magmas in the SOT are larger than those in the MOT and because the injection ratio between subducted AOC and sediment is always >1 in the OT. The distance between the subducting slab and overlying magma may play a significant role in controlling the differences in subduction components injected into magmas between the MOT and SOT.

Structural and Thermal Evolution of an Infant Subduction Shear Zone: Insights From Sub-Ophiolite Metamorphic Rocks Recovered From Oman Drilling Project Site BT-1B

Subduction interface thermal structure changes drastically within the first few million years of underthrusting (i.e., subduction infancy). Metamorphic soles beneath ophiolites record snapshots of dynamic conditions and mechanical coupling during subduction infancy. Beneath the Samail Ophiolite (Oman), the sole comprises structurally higher high-temperature (HT) and lower low-temperature (LT) units. This inverted metamorphic gradient has been attributed to evolving metamorphic Pressure-Temperature (P-T) conditions during infancy; however, peak P-T and timing of LT sole subduction are poorly constrained. Oman Drilling Project core BT-1B sampled the base of the ophiolite in a location lacking the HT sole. Metasedimentary and meta-mafic samples collected from 104 m of core reveal that the LT sole subducted to similar peak P as HT rocks preserved elsewhere in Oman, but experienced ∼300°C lower peak T. Prograde fabrics record Si-in-phengite and amphibole chemistries consistent with peak P-T of ∼7-10 kbar and ∼450-550°C in the epidote-amphibolite facies. Retrograde fabrics record a transition from near-pervasive ductile to localized brittle strain under greenschist facies conditions. Titanite U-Pb ages (n = 2) constrain timing of peak LT sole subduction to ∼91 Ma (post-dating initial HT sole subduction by ∼12-13 Myr) and dynamic retrogression through ∼90 Ma. Combined with existing geo/thermo-chronology, our results support a model of protracted subduction and accretion while the infant subduction zone experienced multi-phase, slow-fast-slow cooling. Temporal overlap of HT sole cooling (rehydration?) and ophiolite formation suggests that cooling may lead to interface weakening, facilitating upper-plate extension and spreading. The LT sole formed in a rapidly-refrigerating forearc after ophiolite formation and may reflect the transition to self-sustaining subduction.

Deep carbon cycle constrained by carbonate solubility

Earth's deep carbon cycle affects atmospheric CO₂, climate, and habitability. Owing to the extreme solubility of CaCO₃, aqueous fluids released from the subducting slab could extract all carbon from the slab. However, recycling efficiency is estimated at only around 40%. Data from carbonate inclusions, petrology, and Mg isotope systematics indicate Ca²⁺ in carbonates is replaced by Mg²⁺ and other cations during subduction. Here we determined the solubility of dolomite [CaMg(CO₃)₂] and rhodochrosite (MnCO₃), and put an upper limit on that of magnesite (MgCO₃) under subduction zone conditions. Solubility decreases at least two orders of magnitude as carbonates become Mg-rich. This decreased solubility, coupled with heterogeneity of carbon and water subduction, may explain discrepancies in carbon recycling estimates. Over a range of slab settings, we find aqueous dissolution responsible for mobilizing 10 to 92% of slab carbon. Globally, aqueous fluids mobilise [Formula: see text]% ([Formula: see text] Mt/yr) of subducted carbon from subducting slabs.

Fluid-rich subducting topography generates anomalous forearc porosity

The role of subducting topography on the mode of fault slip-particularly whether it hinders or facilitates large megathrust earthquakes-remains a controversial topic in subduction dynamics^1-5. Models have illustrated the potential for subducting topography to severely alter the structure, stress state and mechanics of subduction zones^4,6; however, direct geophysical imaging of the complex fracture networks proposed and the hydrology of both the subducting topography and the associated upper plate damage zones remains elusive. Here we use passive and controlled-source seafloor electromagnetic data collected at the northern Hikurangi Margin, New Zealand, to constrain electrical resistivity in a region of active seamount subduction. We show that a seamount on the incoming plate contains a thin, low-porosity basaltic cap that traps a conductive matrix of porous volcaniclastics and altered material over a resistive core, which allows 3.2 to 4.7 times more water to subduct, compared with normal, unfaulted oceanic lithosphere. In the forearc, we image a sediment-starved plate interface above a subducting seamount with similar electrical structure to the incoming plate seamount. A sharp resistive peak within the subducting seamount lies directly beneath a prominent upper plate conductive anomaly. The coincidence of this upper plate anomaly with the location of burst-type repeating earthquakes and seismicity associated with a recent slow slip event⁷ directly links subducting topography to the creation of fluid-rich damage zones in the forearc that alter the effective normal stress at the plate interface by modulating the fluid overpressure. In addition to severely modifying the structure and physical conditions of the upper plate, subducting seamounts represent an underappreciated mechanism for transporting a considerable flux of water to the forearc and deeper mantle.

Role of warm subduction in the seismological properties of the forearc mantle: An example from southwest Japan

A warm slab thermal structure plays an important role in controlling seismic properties of the slab and mantle wedge. Among warm subduction zones, most notably in southwest Japan, the spatial distribution of large S-wave delay times and deep nonvolcanic tremors in the forearc mantle indicate the presence of a serpentinite layer along the slab interface. However, the conditions under which such a layer is generated remains unclear. Using numerical models, we here show that a serpentinite layer begins to develop by the slab-derived fluids below the deeper end of the slab-mantle decoupling interface and grows toward the corner of the mantle wedge along the interface under warm subduction conditions only, explaining the large S-wave delay times in the forearc mantle. The serpentinite layer then allows continuous free-fluid flow toward the corner of the mantle wedge, presenting possible mechanisms for the deep nonvolcanic tremors in the forearc mantle.

Water transport to the core-mantle boundary

Water is transported to Earth's interior in lithospheric slabs at subduction zones. Shallow dehydration fuels hydrous island arc magmatism but some water is transported deeper in cool slab mantle. Further dehydration at ∼700 km may limit deeper transport but hydrated phases in slab crust have considerable capacity for transporting water to the core-mantle boundary. Quantifying how much remains the challenge.

What's down there? The structures, materials and environment of deep-seated slow slip and tremor

Deep-seated slow slip and tremor (SST), including slow slip events, episodic tremor and slip, and low-frequency earthquakes, occur downdip of the seismogenic zone of numerous subduction megathrusts and plate boundary strike-slip faults. These events represent a fascinating and perplexing mode of fault failure that has greatly broadened our view of earthquake dynamics. In this contribution, we review constraints on SST deformation processes from both geophysical observations of active subduction zones and geological observations of exhumed field analogues. We first provide an overview of what has been learned about the environment, kinematics and dynamics of SST from geodetic and seismologic data. We then describe the materials, deformation mechanisms, and metamorphic and fluid pressure conditions that characterize exhumed rocks from SST source depths. Both the geophysical and geological records strongly suggest the importance of a fluid-rich and high fluid pressure habitat for the SST source region. Additionally, transient deformation features preserved in the rock record, involving combined frictional-viscous shear in regions of mixed lithology and near-lithostatic fluid pressures, may scale with the tremor component of SST. While several open questions remain, it is clear that improved constraints on the materials, environment, structure, and conditions of the plate interface from geophysical imaging and geologic observations will enhance model representations of the boundary conditions and geometry of the SST deformation process. This article is part of a discussion meeting issue 'Understanding earthquakes using the geological record'.

The stability of subducted glaucophane with the Earth's secular cooling

The blueschist to eclogite transition is one of the major geochemical–metamorphic processes typifying the subduction zone, which releases fluids triggering earthquakes and arc volcanism. Although glaucophane is an index hydrous mineral for the blueschist facies, its stability at mantle depths in diverse subduction regimes of contemporary and early Earth has not been experimentally determined. Here, we show that the maximum depth of glaucophane stability increases with decreasing thermal gradients of the subduction system. Along cold subduction geotherm, glaucophane remains stable down ca. 240 km depth, whereas it dehydrates and breaks down at as shallow as ca. 40 km depth under warm subduction geotherm or the Proterozoic tectonic setting. Our results imply that secular cooling of the Earth has extended the stability of glaucophane and consequently enabled the transportation of water into deeper interior of the Earth, suppressing arc magmatism, volcanism, and seismic activities along subduction zones.

Along the cold subduction geotherm, glaucophane remains stable down to pressure and temperature (P–T) conditions of ca. 240 km depth, whereas under the warm subduction geotherm, it dehydrates and breaks down into pyroxenes and silica between ca. 50 and 100 km depths.

Field-Based Evidence for Intra-Slab High-Permeability Channel Formation at Eclogite-Facies Conditions During Subduction

Fluid release from subducting oceanic lithosphere is a key process for subduction zone geodynamics, from controlling arc volcanism to seismicity and tectonic exhumation. However, many fundamental details of fluid composition, flow pathways, and reactivity with slab-forming rocks remain to be thoroughly understood. In this study we investigate a multi-kilometer-long, high-pressure metasomatic system preserved in the lawsonite-eclogite metamorphic unit of Alpine Corsica, France. The fluid-mediated process was localized along a major intra-slab interface, which is the contact between basement and cover unit. Two distinct metasomatic stages are identified and discussed. We show that these two stages resulted from the infiltration of deep fluids that were derived from the same source and had the same slab-parallel, updip flow direction. By mass balance analysis, we quantify metasomatic mass changes along this fluid pathway and the time-integrated fluid fluxes responsible for them. In addition, we also assess carbon fluxes associated with these metasomatic events. The magnitude of the estimated fluid fluxes (104-105) indicates that major intra-slab interfaces such as lithological boundaries acted as fluid channels facilitating episodic pulses of fluid flow. We also show that when fluids are channelized, high time-integrated fluid fluxes lead to carbon fluxes several orders of magnitude higher than carbon fluxes generated by local dehydration reactions. Given the size and geologic features of the investigated metasomatic system, we propose that it represents the first reported natural analogue of the so-called high permeability channels predicted by numerical simulations.

Melting of subducted sediments reconciles geophysical images of subduction zones

Sediments play a key role in subduction. They help control the chemistry of arc volcanoes and the location of seismic hazards. Here, we present a new model describing the fate of subducted sediments that explains magnetotelluric models of subduction zones, which commonly show an enigmatic conductive anomaly at the trenchward side of volcanic arcs. In many subduction zones, sediments will melt trenchward of the source region for arc melts. High-pressure experiments show that these sediment melts will react with the overlying mantle wedge to produce electrically conductive phlogopite pyroxenites. Modelling of the Cascadia and Kyushu subduction zones shows that the products of sediment melting closely reproduce the magnetotelluric observations. Melting of subducted sediments can also explain K-rich volcanic rocks that are produced when the phlogopite pyroxenites melt during slab roll-back events. This process may also help constrain models for subduction zone seismicity. Since melts and phlogopite both have low frictional strength, damaging thrust earthquakes are unlikely to occur in the vicinity of the melting sediments, while increased fluid pressures may promote the occurrence of small magnitude earthquakes and episodic tremor and slip.

The role of the antigorite + brucite to olivine reaction in subducted serpentinites (Zermatt, Switzerland)

Metamorphic olivine formed by the reaction of antigorite + brucite is widespread in serpentinites that crop out in glacier-polished outcrops at the Unterer Theodulglacier, Zermatt. Olivine overgrows a relic magnetite mesh texture formed during ocean floor serpentinization. Serpentinization is associated with rodingitisation of mafic dykes. Metamorphic olivine coexists with magnetite, shows high Mg# of 94-97 and low trace element contents. A notable exception is 4 µg/g Boron (> 10 times primitive mantle), introduced during seafloor alteration and retained in metamorphic olivine. Olivine incorporated 100-140 µg/g H₂O in Si-vacancies, providing evidence for low SiO₂-activity imposed by brucite during olivine growth. No signs for hydrogen loss or major and minor element diffusional equilibration are observed. The occurrence of olivine in patches within the serpentinite mimics the former heterogeneous distribution of brucite, whereas the network of olivine-bearing veins and shear zones document the pathways of the escaping fluid produced by the olivine forming reaction. Relic Cr-spinels have a high Cr# of 0.5 and the serpentinites display little or no clinopyroxene, indicating that they derive from hydrated harzburgitic mantle that underwent significant melt depletion. The enrichment of Mg and depletion of Si results in the formation of brucite during seafloor alteration, a pre-requisite for later subduction-related olivine formation and fluid liberation. The comparison of calculated bulk rock brucite contents in the Zermatt-Saas with average IODP serpentinites suggests a large variation in fluid release during olivine formation. Between 3.4 and 7.2 wt% H₂O is released depending on the magnetite content in fully serpentinized harzburgites (average oceanic serpentinites). Thermodynamic modelling indicates that the fluid release in Zermatt occurred between 480 °C and 550 °C at 2-2.5 GPa with the Mg# of olivine varying from 68 to 95. However, the majority of the fluid released from this reaction was produced within a narrow temperature field of < 30 °C, at higher pressures 2.5 GPa and temperatures 550-600 °C than commonly thought. Fluids derived from the antigorite + brucite reaction might thus trigger eclogite facies equilibration in associated metabasalts, meta-gabbros, meta-rodingites and meta-sediments in the area. This focused fluid release has the potential to trigger intermediate depths earthquakes at 60-80 km in subducted oceanic lithosphere.

Contrasting volcano spacing along SW Japan arc caused by difference in age of subducting lithosphere

The SW Japan arc built by subduction of the Philippine Sea (PHS) plate exhibits uneven distribution of volcanoes: thirteen Quaternary composite volcanoes form in the western half of this arc, Kyushu Island, while only two in the eastern half, Chugoku district. Reconstruction of the PHS plate back to 14 Ma, together with examinations based on thermal structure models constrained by high-density heat flow data and a petrological model for dehydration reactions suggest that fluids are discharged actively at depths of 90–100 km in the hydrous layer at the top of the old (> 50 Ma), hence, cold lithosphere sinking beneath Kyushu Island. In contrast, the young (15–25 Ma) oceanic crust downgoing beneath Chugoku district releases fluids largely at shallower depths, i.e. beneath the non-volcanic forearc, to cause characteristic tectonic tremors and low-frequency earthquakes (LFEs) and be the source of specific brine springs. Much larger amounts of fluids supplied to the magma source region in the western SW Japan arc could build more densely-distributed volcanoes.

Variable water input controls evolution of the Lesser Antilles volcanic arc

Oceanic lithosphere carries volatiles, notably water, into the mantle through subduction at convergent plate boundaries. This subducted water exercises control on the production of magma, earthquakes, formation of continental crust and mineral resources. Identifying different potential fluid sources (sediments, crust and mantle lithosphere) and tracing fluids from their release to the surface has proved challenging¹. Atlantic subduction zones are a valuable endmember when studying this deep water cycle because hydration in Atlantic lithosphere, produced by slow spreading, is expected to be highly non-uniform². Here, as part of a multi-disciplinary project in the Lesser Antilles volcanic arc³, we studied boron trace element and isotopic fingerprints of melt inclusions. These reveal that serpentine-that is, hydrated mantle rather than crust or sediments-is a dominant supplier of subducted water to the central arc. This serpentine is most likely to reside in a set of major fracture zones subducted beneath the central arc over approximately the past ten million years. The current dehydration of these fracture zones coincides with the current locations of the highest rates of earthquakes and prominent low shear velocities, whereas the preceding history of dehydration is consistent with the locations of higher volcanic productivity and thicker arc crust. These combined geochemical and geophysical data indicate that the structure and hydration of the subducted plate are directly connected to the evolution of the arc and its associated seismic and volcanic hazards.

Strain-Induced Fabric Transition of Chlorite and Implications for Seismic Anisotropy in Subduction Zones

Seismic anisotropy of S-wave, trench-parallel or trench-normal polarization direction of fast S-wave, has been observed in the fore-arc and back-arc regions of subduction zones. Lattice preferred orientation (LPO) of elastically anisotropic chlorite has been suggested as one of the major causes of seismic anisotropy in subduction zones. However, there are two different LPOs of chlorite reported based on the previous studies of natural chlorite peridotites, which can produce different expression of seismic anisotropy. The mechanism for causing the two different LPOs of chlorite is not known. Therefore, we conducted deformation experiments of chlorite peridotite under high pressure–temperature conditions (P = 0.5–2.5 GPa, T = 540–720 °C). We found that two different chlorite LPOs were developed depending on the magnitude of shear strain. The type-1 chlorite LPO is characterized by the [001] axes aligned subnormal to the shear plane, and the type-2 chlorite LPO is characterized by a girdle distribution of the [001] axes subnormal to the shear direction. The type-1 chlorite LPO developed under low shear strain (γ ≤ 3.1 ± 0.3), producing trench-parallel seismic anisotropy. The type-2 chlorite LPO developed under high shear strain (γ ≥ 5.1 ± 1.5), producing trench-normal seismic anisotropy. The anisotropy of S-wave velocity (AVs) of chlorite was very strong up to AVs = 48.7% so that anomalous seismic anisotropy in subduction zones can be influenced by the chlorite LPOs.

Rapid subduction initiation and magmatism in the Western Pacific driven by internal vertical forces

Plate tectonics requires the formation of plate boundaries. Particularly important is the enigmatic initiation of subduction: the sliding of one plate below the other, and the primary driver of plate tectonics. A continuous, in situ record of subduction initiation was recovered by the International Ocean Discovery Program Expedition 352, which drilled a segment of the fore-arc of the Izu-Bonin-Mariana subduction system, revealing a distinct magmatic progression with a rapid timescale (approximately 1 million years). Here, using numerical models, we demonstrate that these observations cannot be produced by previously proposed horizontal external forcing. Instead a geodynamic evolution that is dominated by internal, vertical forces produces both the temporal and spatial distribution of magmatic products, and progresses to self-sustained subduction. Such a primarily internally driven initiation event is necessarily whole-plate scale and the rock sequence generated (also found along the Tethyan margin) may be considered as a smoking gun for this type of event.

The magmatic progression produced during the initiation of the Izu-Bonin-Marianas subduction zone took place rapidly over 1 million years, but it has been unclear why. Here, using numerical models, the authors show that subduction initiation was dominated by vertical forces, internal to the system itself, progressing to self-sustained subduction.

Intraplate volcanism originating from upwelling hydrous mantle transition zone

Most magmatism occurring on Earth is conventionally attributed to passive mantle upwelling at mid-ocean ridges, to slab devolatilization at subduction zones, or to mantle plumes. However, the widespread Cenozoic intraplate volcanism in northeast China^1-3 and the young petit-spot volcanoes^4-7 offshore of the Japan Trench cannot readily be associated with any of these mechanisms. In addition, the mantle beneath these types of volcanism is characterized by zones of anomalously low seismic velocity above and below the transition zone^8-12 (a mantle level located at depths between 410 and 660 kilometres). A comprehensive interpretation of these phenomena is lacking. Here we show that most (or possibly all) of the intraplate and petit-spot volcanism and low-velocity zones around the Japanese subduction zone can be explained by the Cenozoic interaction of the subducting Pacific slab with a hydrous mantle transition zone. Numerical modelling indicates that 0.2 to 0.3 weight per cent of water dissolved in mantle minerals that are driven out from the transition zone in response to subduction and retreat of a tectonic plate is sufficient to reproduce the observations. This suggests that a critical amount of water may have accumulated in the transition zone around this subduction zone, as well as in others of the Tethyan tectonic belt¹³ that are characterized by intraplate or petit-spot volcanism and low-velocity zones in the underlying mantle.

Evidence for the stability of ultrahydrous stishovite in Earth's lower mantle

The distribution and transportation of water in Earth's interior depends on the stability of water-bearing phases. The transition zone in Earth's mantle is generally accepted as an important potential water reservoir because its main constituents, wadsleyite and ringwoodite, can incorporate weight percent levels of H2O in their structures at mantle temperatures. The extent to which water can be transported beyond the transition zone deeper into the mantle depends on the water carrying capacity of minerals stable in subducted lithosphere. Stishovite is one of the major mineral components in subducting oceanic crust, yet the capacity of stishovite to incorporate water beyond at lower mantle conditions remains speculative. In this study, we combine in situ laser heating with synchrotron X-ray diffraction to show that the unit cell volume of stishovite synthesized under hydrous conditions is ∼2.3 to 5.0% greater than that of anhydrous stishovite at pressures of ∼27 to 58 GPa and temperatures of 1,240 to 1,835 K. Our results indicate that stishovite, even at temperatures along a mantle geotherm, can potentially incorporate weight percent levels of H2O in its crystal structure and has the potential to be a key phase for transporting and storing water in the lower mantle.

Hydrogen isotopic evidence for early oxidation of silicate Earth

The Moon-forming giant impact extensively melts and partially vaporizes the silicate Earth and delivers a substantial mass of metal to Earth's core. The subsequent evolution of the magma ocean and overlying atmosphere has been described by theoretical models but observable constraints on this epoch have proved elusive. Here, we report thermodynamic and climate calculations of the primordial atmosphere during the magma ocean and water ocean epochs respectively and forge new links with observations to gain insight into the behavior of volatiles on the Hadean Earth. As accretion wanes, Earth's magma ocean crystallizes, outgassing the bulk of its volatiles into the primordial atmosphere. The redox state of the magma ocean controls both the chemical composition of the outgassed volatiles and the hydrogen isotopic composition of water oceans that remain after hydrogen escape from the primordial atmosphere. The climate modeling indicates that multi-bar H₂-rich atmospheres generate sufficient greenhouse warming and rapid kinetics resulting in ocean-atmosphere H₂O-H₂ isotopic equilibration. Whereas water condenses and is mostly retained, molecular hydrogen does not condense and can escape, allowing large quantities (~10² bars) of hydrogen - if present - to be lost from the Earth in this epoch. Because the escaping inventory of H can be comparable to the hydrogen inventory in primordial water oceans, equilibrium deuterium enrichment can be large with a magnitude that depends on the initial atmospheric H₂ inventory. Under equilibrium partitioning, the water molecule concentrates deuterium and, to the extent that hydrogen in other forms (e.g., H₂) are significant species in the outgassed atmosphere, pronounced D/H enrichments (~1.5-2x) in the oceans are expected from equilibrium partitioning in this epoch. By contrast, the common view that terrestrial water has a carbonaceous chondritic source requires the oceans to preserve the isotopic composition of that source, undergoing minimal D-enrichment via equilibration with H₂ followed by hydrodynamic escape. Such minimal enrichment places upper limits on the amount of primordial atmospheric H₂ in contact with Hadean water oceans and implies oxidizing conditions (logfO₂>IW+1, H₂/H₂O<0.3) for outgassing from the magma ocean. Preservation of an approximate carbonaceous chondrite D/H signature in the oceans thus provides evidence that the observed oxidation of silicate Earth occurred before crystallization of the final magma ocean, yielding a new constraint on the timing of this critical event in Earth history. The seawater-carbonaceous chondrite "match" in D/H (to ~10-20%) further constrains the prior existence of an atmospheric H₂ inventory - of any origin - on post-giant-impact Earth to <20 bars, and suggests that the terrestrial mantle supplied the oxidant for the chemical resorption of metals during terrestrial late accretion.

Molybdenum systematics of subducted crust record reactive fluid flow from underlying slab serpentine dehydration

Fluids liberated from subducting slabs are critical in global geochemical cycles. We investigate the behaviour of Mo during slab dehydration using two suites of exhumed fragments of subducted, oceanic lithosphere. Our samples display a positive correlation of δ^98/95Mo_{NIST 3134} with Mo/Ce, from compositions close to typical mantle (−0.2‰ and 0.03, respectively) to very low values of both δ^98/95Mo_{NIST 3134} (−1‰) and Mo/Ce (0.002). Together with new, experimental data, we show that molybdenum isotopic fractionation is driven by preference of heavier Mo isotopes for a fluid phase over rutile, the dominant mineral host of Mo in eclogites. Moreover, the strongly perturbed δ^98/95Mo_{NIST 3134} and Mo/Ce of our samples requires that they experienced a large flux of oxidised fluid. This is consistent with channelised, reactive fluid flow through the subducted crust, following dehydration of the underlying, serpentinised slab mantle. The high δ^98/95Mo_{NIST 3134} of some arc lavas is the complement to this process.

Fluid liberation and migration from subducted oceanic slabs play a critical role in arc magmatism but the volume and origin of the released fluids is unclear and difficult to trace. Here, the authors use Molybdenum isotope ratios to tackle these problems.

Silicate Melt Inclusions in Diamonds of Eclogite Paragenesis from Placers on the Northeastern Siberian Craton

New findings of silicate-melt inclusions in two alluvial diamonds (from the Kholomolokh placer, northeastern Siberian Platform) are reported. Both diamonds exhibit a high degree of N aggregation state (60–70% B) suggesting their long residence in the mantle. Raman spectral analysis revealed that the composite inclusions consist of clinopyroxene and silicate glass. Hopper crystals of clinopyroxene were observed using scanning electron microscopy and energy-dispersive spectroscopic analyses; these are different in composition from the omphacite inclusions that co-exist in the same diamonds. The glasses in these inclusions contain relatively high SiO2, Al2O3, Na2O and, K2O. These composite inclusions are primary melt that partially crystallised at the cooling stage. Hopper crystals of clinopyroxene imply rapid cooling rates, likely related to the uplift of crystals in the kimberlite melt. The reconstructed composition of such primary melts suggests that they were formed as the product of metasomatised mantle. One of the most likely source of melts/fluids metasomatising the mantle could be a subducted slab.

Earthquakes track subduction fluids from slab source to mantle wedge sink

Subducting plates release fluids as they plunge into Earth's mantle and occasionally rupture to produce intraslab earthquakes. It is debated whether fluids and earthquakes are directly related. By combining seismic observations and geodynamic models from western Greece, and comparing across other subduction zones, we find that earthquakes effectively track the flow of fluids from their slab source at >80 km depth to their sink at shallow (<40 km) depth. Between source and sink, the fluids flow updip under a sealed plate interface, facilitating intraslab earthquakes. In some locations, the seal breaks and fluids escape through vents into the mantle wedge, thereby reducing the fluid supply and seismicity updip in the slab. The vents themselves may represent nucleation sites for larger damaging earthquakes.

Ferrous Iron Under Oxygen-Rich Conditions in the Deep Mantle

Recent experiments have demonstrated the existence of previously unknown iron oxides at high pressure and temperature including newly discovered pyrite-type FeO2 and FeO2Hx phases stable at deep terrestrial lower mantle pressures and temperatures. In the present study, we probed the iron oxidation state in high-pressure transformation products of Fe3+OOH goethite by in situ X-ray absorption spectroscopy in laser-heated diamond-anvil cell. At pressures and temperatures of ~91 GPa and 1,500-2,350 K, respectively, that is, in the previously reported stability field of FeO2Hx, a measured shift of -3.3 ± 0.1 eV of the Fe K-edge demonstrates that iron has turned from Fe3+ to Fe2+. We interpret this reductive valence change of iron by a concomitant oxidation of oxygen atoms from O2- to O-, in agreement with previous suggestions based on the structures of pyrite-type FeO2 and FeO2Hx phases. Such peculiar chemistry could drastically change our view of crystal chemistry in deep planetary interiors.

Water input into the Mariana subduction zone estimated from ocean-bottom seismic data

The water cycle at subduction zones remains poorly understood, although subduction is the only mechanism for water transport deep into Earth. Previous estimates of water flux^1-3 exhibit large variations in the amount of water that is subducted deeper than 100 kilometres. The main source of uncertainty in these calculations is the initial water content of the subducting uppermost mantle. Previous active-source seismic studies suggest that the subducting slab may be pervasively hydrated in the plate-bending region near the oceanic trench^4-7. However, these studies do not constrain the depth extent of hydration and most investigate young incoming plates, leaving subduction-zone water budgets for old subducting plates uncertain. Here we present seismic images of the crust and uppermost mantle around the central Mariana trench derived from Rayleigh-wave analysis of broadband ocean-bottom seismic data. These images show that the low mantle velocities that result from mantle hydration extend roughly 24 kilometres beneath the Moho discontinuity. Combined with estimates of subducting crustal water, these results indicate that at least 4.3 times more water subducts than previously calculated for this region³. If other old, cold subducting slabs contain correspondingly thick layers of hydrous mantle, as suggested by the similarity of incoming plate faulting across old, cold subducting slabs, then estimates of the global water flux into the mantle at depths greater than 100 kilometres must be increased by a factor of about three compared to previous estimates³. Because a long-term net influx of water to the deep interior of Earth is inconsistent with the geological record⁸, estimates of water expelled at volcanic arcs and backarc basins probably also need to be revised upwards⁹.

Shear heating reconciles thermal models with the metamorphic rock record of subduction

Some commonly referenced thermal-mechanical models of current subduction zones imply temperatures that are 100-500 °C colder at 30-80-km depth than pressure-temperature conditions determined thermobarometrically from exhumed metamorphic rocks. Accurately inferring subduction zone thermal structure, whether from models or rocks, is crucial for predicting metamorphic reactions and associated fluid release, subarc melting conditions, rheologies, and fault-slip phenomena. Here, we compile surface heat flow data from subduction zones worldwide and show that values are higher than can be explained for a frictionless subduction interface often assumed for modeling. An additional heat source--likely shear heating--is required to explain these forearc heat flow values. A friction coefficient of at least 0.03 and possibly as high as 0.1 in some cases explains these data, and we recommend a provisional average value of 0.05 ± 0.015 for modeling. Even small coefficients of friction can contribute several hundred degrees of heating at depths of 30-80 km. Adding such shear stresses to thermal models quantitatively reproduces the pressure-temperature conditions recorded by exhumed metamorphic rocks. Comparatively higher temperatures generally drive rock dehydration and densification, so, at a given depth, hotter rocks are denser than colder rocks, and harder to exhume through buoyancy mechanisms. Consequently--conversely to previous proposals--exhumed metamorphic rocks might overrepresent old-cold subduction where rocks at the slab interface are wetter and more buoyant than in young-hot subduction zones.

Xenon isotopic constraints on the history of volatile recycling into the mantle

The long-term exchange of volatile species (such as water, carbon, nitrogen and the noble gases) between deep Earth and surface reservoirs controls the habitability of the Earth's surface. The present-day volatile budget of the mantle reflects the integrated history of outgassing and retention of primordial volatiles delivered to the planet during accretion, volatile species generated by radiogenic ingrowth and volatiles transported into the mantle from surface reservoirs over time. Variations in the distribution of volatiles between deep Earth and surface reservoirs affect the viscosity, cooling rate and convective stress state of the solid Earth. Accordingly, constraints on the flux of surface volatiles transported into the deep Earth improve our understanding of mantle convection and plate tectonics. However, the history of surface volatile regassing into the mantle is not known. Here we use mantle xenon isotope systematics to constrain the age of initiation of volatile regassing into the deep Earth. Given evidence of prolonged evolution of the xenon isotopic composition of the atmosphere^1,2, we find that substantial recycling of atmospheric xenon into the deep Earth could not have occurred before 2.5 billion years ago. Xenon concentrations in downwellings remained low relative to ambient convecting mantle concentrations throughout the Archaean era, and the mantle shifted from a net degassing to a net regassing regime after 2.5 billion years ago. Because xenon is carried into the Earth's interior in hydrous mineral phases^3-5, our results indicate that downwellings were drier in the Archaean era relative to the present. Progressive drying of the Archean mantle would allow slower convection and decreased heat transport out of the mantle, suggesting non-monotonic thermal evolution of the Earth's interior.

Three-dimensional variations of the slab geometry correlate with earthquake distributions at the Cascadia subduction system

Significant along-strike variations of seismicity are observed at subduction zones, which are strongly influenced by physical properties of the plate interface and rheology of the crust and mantle lithosphere. However, the role of the oceanic side of the plate boundary on seismicity is poorly understood due to the lack of offshore instrumentations. Here tomographic results of the Cascadia subduction system, resolved with full-wave ambient noise simulation and inversion by integrating dense offshore and onshore seismic datasets, show significant variations of the oceanic lithosphere along strike and down dip from spreading centers to subduction. In central Cascadia, where seismicity is sparse, the slab is imaged as a large-scale low-velocity feature near the trench, which is attributed to a highly hydrated and strained oceanic lithosphere underlain by a layer of melts or fluids. The strong correlation suggests that the properties of the incoming oceanic plate play a significant role on seismicity.

Variations in seismicity are observed at subduction zones, but the oceanic sides remain poorly resolved. Here, the author presents tomographic results of the Cascadia subduction system demonstrating that there are significant variations of the oceanic lithosphere along the subduction zone.

Discovery of a hexagonal ultradense hydrous phase in (Fe,Al)OOH

A deep lower-mantle (DLM) water reservoir depends on availability of hydrous minerals which can store and transport water into the DLM without dehydration. Recent discoveries found hydrous phases AlOOH (Z = 2) with a CaCl₂-type structure and FeOOH (Z = 4) with a cubic pyrite-type structure stable under the high-pressure-temperature (P-T) conditions of the DLM. Our experiments at 107-136 GPa and 2,400 K have further demonstrated that (Fe,Al)OOH is stabilized in a hexagonal lattice. By combining powder X-ray-diffraction techniques with multigrain indexation, we are able to determine this hexagonal hydrous phase with a = 10.5803(6) Å and c = 2.5897(3) Å at 110 GPa. Hexagonal (Fe,Al)OOH can transform to the cubic pyrite structure at low T with the same density. The hexagonal phase can be formed when δ-AlOOH incorporates FeOOH produced by reaction between water and Fe, which may store a substantial quantity of water in the DLM.