Deep carbon cycle constrained by carbonate solubility

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.

The effects of local variations in conditions on carbon storage and release in the continental mantle

Recent advances indicate that the amount of carbon released by gradual degassing from the mantle needs to be revised upwards, whereas the carbon supplied by plumes may have been overestimated in the past. Variations in rock types and oxidation state may be very local and exert strong influences on carbon storage and release mechanisms. Deep subduction may be prevented by diapirism in thick sedimentary packages, whereas carbonates in thinner sequences may be subducted. Carbonates stored in the mantle transition zone will melt when they heat up, recognized by coupled stable isotope systems (e.g. Mg, Zn, Ca). There is no single 'mantle oxygen fugacity', particularly in the thermal boundary layer (TBL) and lowermost lithosphere, where very local mixtures of rock types coexist. Carbonate-rich melts from either subduction or melting of the uppermost asthenosphere trap carbon by redox freezing or as carbonate-rich dykes in this zone. Deeply derived, reduced melts may form further diamond reservoirs, recognized as polycrystalline diamonds associated with websteritic silicate minerals. Carbon is released by either edge-driven convection, which tears sections of the TBL and lower lithosphere down so that they melt by a mixture of heating and oxidation, or by lateral advection of solids beneath rifts. Both mechanisms operate at steps in lithosphere thickness and result in carbonate-rich melts, explaining the spatial association of craton edges and carbonate-rich magmatism. High-pressure experiments on individual rock types, and increasingly on reactions between rocks and melts, are fine-tuning our understanding of processes and turning up unexpected results that are not seen in studies of single rocks. Future research should concentrate on elucidating local variations and integrating these with the interpretation of geophysical signals. Global concepts such as average sediment compositions and a uniform mantle oxidation state are not appropriate for small-scale processes; an increased focus on local variations will help to refine carbon budget models.

Earthquake-enhanced dissolved carbon cycles in ultra-deep ocean sediments

Hadal trenches are unique geological and ecological systems located along subduction zones. Earthquake-triggered turbidites act as efficient transport pathways of organic carbon (OC), yet remineralization and transformation of OC in these systems are not comprehensively understood. Here we measure concentrations and stable- and radiocarbon isotope signatures of dissolved organic and inorganic carbon (DOC, DIC) in the subsurface sediment interstitial water along the Japan Trench axis collected during the IODP Expedition 386. We find accumulation and aging of DOC and DIC in the subsurface sediments, which we interpret as enhanced production of labile dissolved carbon owing to earthquake-triggered turbidites, which supports intensive microbial methanogenesis in the trench sediments. The residual dissolved carbon accumulates in deep subsurface sediments and may continue to fuel the deep biosphere. Tectonic events can therefore enhance carbon accumulation and stimulate carbon transformation in plate convergent trench systems, which may accelerate carbon export into the subduction zones.

Carbonate-rich crust subduction drives the deep carbon and chlorine cycles

The flux balances of carbon and chlorine between subduction into the deep mantle and volcanic emissions into the atmosphere are crucial for the habitability of our planet1,2. However, pervasive loss of fluids from subducting slabs has been thought to cut off the delivery of both carbon and chlorine to the deep mantle owing to their high mobility under hydrous conditions3,4. Our new high-pressure experiments show that most carbonates (>75 wt%) in carbonate-rich crustal rocks-one of the main subducting carbon reservoirs-survive devolatilization and hydrous melting in cold and warm subduction zones, indicating that their subduction has driven the deep carbon cycle since the Mesoproterozoic. We found that KCl and NaCl, respectively, become stable phases crystallizing from hydrous carbonatite melts with low chlorine solubility in warm and hot subduction zones, resulting in the sequestration of chlorine in the solid residue in downwelling slabs. Accordingly, the subduction of carbonate-rich rocks facilitated highly effective recycling of both chlorine and carbon into the deep mantle at intermediate stages of Earth's history and led to declining atmospheric pCO2 and the formation of carbon-rich and chlorine-rich mantle reservoirs since the Mesoproterozoic. This period of optimal carbon and chlorine subduction may explain the ages of eclogitic diamonds and the formation of the HIMU mantle source.

Structural water in amorphous carbonate minerals: ab initio molecular dynamics simulations of X-ray pair distribution experiments

Water is known to play a controlling role in directing mineralization pathways and stabilizing metastable amorphous intermediates in hydrous carbonate mineral MCO₃·nH₂O systems, where M²⁺ is a divalent metal cation. Despite this recognition, the nature of the controls on crystallization is poorly understood, largely owing to the difficulty in characterizing the dynamically disordered structures of amorphous intermediates at the atomic scale. Here, we present a series of atomistic models, derived from ab initio molecular dynamics simulation, across a range of experimentally relevant cations (M = Ca, Mg, Sr) and hydration levels (0 ≤ n ≤ 2). Theoretical simulations of the dependence of the X-ray pair distribution function on the hydration level n show good agreement with available experimental data and thus provide further evidence for a lack of significant nanoscale structure in amorphous carbonates. Upon dehydration, the metal coordination number does not change significantly, but the relative extent of water dissociation increases, indicating that a thermodynamic driving force exists for water dissociation to accompany dehydration. Mg strongly favors monodentate conformation of carbonate ligands and shows a marked preference to exchange monodentate carbonate O for water O upon hydration, whereas Ca and Sr exchange mono- and bidentate carbonate ligands with comparable frequency. Water forms an extensive hydrogen bond network among both water and carbonate groups that exhibits frequent proton transfers for all three cations considered suggesting that proton mobility is likely predominantly due to water dissociation and proton transfer reactions rather than molecular water diffusion.

Evolution of Earth's tectonic carbon conveyor belt

Concealed deep beneath the oceans is a carbon conveyor belt, propelled by plate tectonics. Our understanding of its modern functioning is underpinned by direct observations, but its variability through time has been poorly quantified. Here we reconstruct oceanic plate carbon reservoirs and track the fate of subducted carbon using thermodynamic modelling. In the Mesozoic era, 250 to 66 million years ago, plate tectonic processes had a pivotal role in driving climate change. Triassic-Jurassic period cooling correlates with a reduction in solid Earth outgassing, whereas Cretaceous period greenhouse conditions can be linked to a doubling in outgassing, driven by high-speed plate tectonics. The associated 'carbon subduction superflux' into the subcontinental mantle may have sparked North American diamond formation. In the Cenozoic era, continental collisions slowed seafloor spreading, reducing tectonically driven outgassing, while deep-sea carbonate sediments emerged as the Earth's largest carbon sink. Subduction and devolatilization of this reservoir beneath volcanic arcs led to a Cenozoic increase in carbon outgassing, surpassing mid-ocean ridges as the dominant source of carbon emissions 20 million years ago. An increase in solid Earth carbon emissions during Cenozoic cooling requires an increase in continental silicate weathering flux to draw down atmospheric carbon dioxide, challenging previous views and providing boundary conditions for future carbon cycle models.