Heterogeneity in mantle carbon content from CO2-undersaturated basalts

Deep mantle earthquakes linked to CO2 degassing at the mid-Atlantic ridge

Volatiles (CO2, H2O) play a fundamental role in mantle melting beneath ocean spreading centers, but what role they play during the melt migration remains unknown. Using seismological data recorded by ocean-bottom seismometers, here we report the presence of deep earthquakes at 10-20 km depth in the mantle along the Mid-Atlantic Ridge axis, much below the brittle-ductile boundary. Syntheses of regional basaltic rock samples and their geochemical analyses indicate the presence of an abnormally high quantity of CO2 (~0.4-3.0 wt%) in the primary melts. As the degassing of a high concentration of dissolved CO2 produces volume change, we suggest that deep earthquakes in the mantle result from the degassing of CO2. The large concentration of CO2 in the primitive melt will influence the presence of melt beneath the lithosphere-asthenosphere boundary at sub-solidus temperatures.

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.

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.

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.

High-pressure elastic properties of dolomite melt supporting carbonate-induced melting in deep upper mantle

Petrologic studies suggest that carbonate-rich melts are present in the Earth’s upper mantle and play an important role in the deep carbon cycle. However, seismic detection of these melts is difficult due to the lack of understanding of the elastic properties of carbonate melts. Here we determined the sound velocity and density of dolomite melt at upper mantle conditions using high-pressure experiments and theoretical simulations. The calculated velocities of carbonate melt-bearing mantle using these new elasticity data were compared with global seismic observations. The result suggests that ∼0.05% carbonate-rich melt can be pervasively present in the deep upper mantle, implying a global average mantle carbon concentration of 80-140 ppm.

Deeply subducted carbonates likely cause low-degree melting of the upper mantle and thus play an important role in the deep carbon cycle. However, direct seismic detection of carbonate-induced partial melts in the Earth’s interior is hindered by our poor knowledge on the elastic properties of carbonate melts. Here we report the first experimentally determined sound velocity and density data on dolomite melt up to 5.9 GPa and 2046 K by in-situ ultrasonic and sink-float techniques, respectively, as well as first-principles molecular dynamics simulations of dolomite melt up to 16 GPa and 3000 K. Using our new elasticity data, the calculated V_P/V_S ratio of the deep upper mantle (∼180–330 km) with a small amount of carbonate-rich melt provides a natural explanation for the elevated V_P/V_S ratio of the upper mantle from global seismic observations, supporting the pervasive presence of a low-degree carbonate-rich partial melt (∼0.05%) that is consistent with the volatile-induced or redox-regulated initial melting in the upper mantle as argued by petrologic studies. This carbonate-rich partial melt region implies a global average carbon (C) concentration of 80–140 ppm. by weight in the deep upper mantle source region, consistent with the mantle carbon content determined from geochemical studies.

Deep roots for mid-ocean-ridge volcanoes revealed by plagioclase-hosted melt inclusions

The global mid-ocean ridge system is the most extensive magmatic system on our planet and is the site of 75 per cent of Earth's volcanism¹. The vertical extent of mid-ocean-ridge magmatic systems has been considered to be restricted: even at the ultraslow-spreading Gakkel mid-ocean ridge under the Arctic Ocean, where the lithosphere is thickest, crystallization depths of magmas that feed eruptions are thought to be less than nine kilometres². These depths were determined using the volatile-element contents of melt inclusions, which are small volumes of magma that become trapped within crystallizing minerals. In studies of basaltic magmatic systems, olivine is the mineral of choice for this approach^2-6. However, pressures derived from olivine-hosted melt inclusions are at odds with pressures derived from basalt major-element barometers⁷ and geophysical measurements of lithospheric thickness⁸. Here we present a comparative study of olivine- and plagioclase-hosted melt inclusions from the Gakkel mid-ocean ridge. We show that the volatile contents of plagioclase-hosted melt inclusions correspond to much higher crystallization pressures (with a mean value of 270 megapascals) than olivine-hosted melt inclusions (with a mean value of 145 megapascals). The highest recorded pressure that we find equates to a depth 16.4 kilometres below the seafloor. Such higher depths are consistent with both the thickness of the Gakkel mid-ocean ridge lithosphere and with pressures reconstructed from glass compositions. In contrast to previous studies using olivine-hosted melt inclusions, our results demonstrate that mid-ocean-ridge volcanoes may have magmatic roots deep in the lithospheric mantle, at least at ultraslow-spreading ridges.

The mechanism of tidal triggering of earthquakes at mid-ocean ridges

The strong tidal triggering of mid-ocean ridge earthquakes has remained unexplained because the earthquakes occur preferentially during low tide, when normal faulting earthquakes should be inhibited. Using Axial Volcano on the Juan de Fuca ridge as an example, we show that the axial magma chamber inflates/deflates in response to tidal stresses, producing Coulomb stresses on the faults that are opposite in sign to those produced by the tides. When the magma chamber's bulk modulus is sufficiently low, the phase of tidal triggering is inverted. We find that the stress dependence of seismicity rate conforms to triggering theory over the entire tidal stress range. There is no triggering stress threshold and stress shadowing is just a continuous function of stress decrease. We find the viscous friction parameter A to be an order of magnitude smaller than laboratory measurements. The high tidal sensitivity at Axial Volcano results from the shallow earthquake depths.

CO2 flux emissions from the Earth's most actively degassing volcanoes, 2005-2015

The global carbon dioxide (CO₂) flux from subaerial volcanoes remains poorly quantified, limiting our understanding of the deep carbon cycle during geologic time and in modern Earth. Past attempts to extrapolate the global volcanic CO₂ flux have been biased by observations being available for a relatively small number of accessible volcanoes. Here, we propose that the strong, but yet unmeasured, CO₂ emissions from several remote degassing volcanoes worldwide can be predicted using regional/global relationships between the CO₂/S_T ratio of volcanic gases and whole-rock trace element compositions (e.g., Ba/La). From these globally linked gas/rock compositions, we predict the CO₂/S_T gas ratio of 34 top-degassing remote volcanoes with no available gas measurements. By scaling to volcanic SO₂ fluxes from a global catalogue, we estimate a cumulative “unmeasured” CO₂ output of 11.4 ± 1.1 Mt/yr (or 0.26 ± 0.02·10¹² mol/yr). In combination with the measured CO₂ output of 27.4 ± 3.6 Mt/yr (or 0.62 ± 0.08·10¹² mol/yr), our results constrain the time-averaged (2005–2015) cumulative CO₂ flux from the Earth’s 91 most actively degassing subaerial volcanoes at 38.7 ± 2.9 Mt/yr (or 0.88 ± 0.06·10¹² mol/yr).

Crystalline polymeric carbon dioxide stable at megabar pressures

Carbon dioxide is a widespread simple molecule in the Universe. In spite of its simplicity it has a very complex phase diagram, forming both amorphous and crystalline extended phases above 40 GPa. The stability range and nature of these phases are still debated, especially in view of their possible role within the deep carbon cycle. Here, we report static synchrotron X-ray diffraction and Raman high-pressure experiments in the megabar range providing evidence for the stability of the polymeric phase V at pressure-temperature conditions relevant to the Earth's lowermost mantle. The equation of state has been extended to 120 GPa and, contrary to earlier experimental findings, neither dissociation into diamond and ε-oxygen nor ionization was observed. Severe deviatoric stress and lattice deformation along with preferred orientation are removed on progressive annealing, thus suggesting CO₂-V as the stable structure also above one megabar.

The 2011 eruption of Nabro volcano, Eritrea: perspectives on magmatic processes from melt inclusions

The 2011 eruption of Nabro volcano, Eritrea, produced one of the largest volcanic sulphur inputs to the atmosphere since the 1991 eruption of Mt. Pinatubo, yet has received comparatively little scientific attention. Nabro forms part of an off-axis alignment, broadly perpendicular to the Afar Rift, and has a history of large-magnitude explosive silicic eruptions, as well as smaller more mafic ones. Here, we present and analyse extensive petrological data obtained from samples of trachybasaltic tephra erupted during the 2011 eruption to assess the pre-eruptive magma storage system and explain the large sulphur emission. We show that the eruption involved two texturally distinct batches of magma, one of which was more primitive and richer in sulphur than the other, which was higher in water (up to 2.5 wt%). Modelling of the degassing and crystallisation histories demonstrates that the more primitive magma rose rapidly from depth and experienced degassing crystallisation, while the other experienced isobaric cooling in the crust at around 5 km depth. Interaction between the two batches occurred shortly before the eruption. The eruption itself was likely triggered by recharge-induced destabilisation of vertically extensive mush zone under the volcano. This could potentially account for the large volume of sulphur released. Some of the melt inclusions are volatile undersaturated, and suggest that the original water content of the magma was around 1.3 wt%, which is relatively high for an intraplate setting, but consistent with seismic studies of the Afar plume. This eruption was smaller than some geological eruptions at Nabro, but provides important insights into the plumbing systems and dynamics of off-axis volcanoes in Afar.