Rapid shifting of a deep magmatic source at Fagradalsfjall volcano, Iceland.

Recent Icelandic rifting events have illuminated the roles of centralized crustal magma reservoirs and lateral magma transport1-4, important characteristics of mid-ocean ridge magmatism1,5. A consequence of such shallow crustal processing of magmas4,5 is the overprinting of signatures that trace the origin, evolution and transport of melts in the uppermost mantle and lowermost crust6,7. Here we present unique insights into processes occurring in this zone from integrated petrologic and geochemical studies of the 2021 Fagradalsfjall eruption on the Reykjanes Peninsula in Iceland. Geochemical analyses of basalts erupted during the first 50 days of the eruption, combined with associated gas emissions, reveal direct sourcing from a near-Moho magma storage zone. Geochemical proxies, which signify different mantle compositions and melting conditions, changed at a rate unparalleled for individual basaltic eruptions globally. Initially, the erupted lava was dominated by melts sourced from the shallowest mantle but over the following three weeks became increasingly dominated by magmas generated at a greater depth. This exceptionally rapid trend in erupted compositions provides an unprecedented temporal record of magma mixing that filters the mantle signal, consistent with processing in near-Moho melt lenses containing 107-108 m3 of basaltic magma. Exposing previously inaccessible parts of this key magma processing zone to near-real-time investigations provides new insights into the timescales and operational mode of basaltic magma systems.

High 3He/4He in central Panama reveals a distal connection to the Galápagos plume.

It is well established that mantle plumes are the main conduits for upwelling geochemically enriched material from Earth's deep interior. The fashion and extent to which lateral flow processes at shallow depths may disperse enriched mantle material far (>1,000 km) from vertical plume conduits, however, remain poorly constrained. Here, we report He and C isotope data from 65 hydrothermal fluids from the southern Central America Margin (CAM) which reveal strikingly high 3He/4He (up to 8.9RA) in low-temperature (≤50 °C) geothermal springs of central Panama that are not associated with active volcanism. Following radiogenic correction, these data imply a mantle source 3He/4He >10.3RA (and potentially up to 26RA, similar to Galápagos hotspot lavas) markedly greater than the upper mantle range (8 ± 1RA). Lava geochemistry (Pb isotopes, Nb/U, and Ce/Pb) and geophysical constraints show that high 3He/4He values in central Panama are likely derived from the infiltration of a Galápagos plume-like mantle through a slab window that opened ∼8 Mya. Two potential transport mechanisms can explain the connection between the Galápagos plume and the slab window: 1) sublithospheric transport of Galápagos plume material channeled by lithosphere thinning along the Panama Fracture Zone or 2) active upwelling of Galápagos plume material blown by a "mantle wind" toward the CAM. We present a model of global mantle flow that supports the second mechanism, whereby most of the eastward transport of Galápagos plume material occurs in the shallow asthenosphere. These findings underscore the potential for lateral mantle flow to transport mantle geochemical heterogeneities thousands of kilometers away from plume conduits.

Ancient helium and tungsten isotopic signatures preserved in mantle domains least modified by crustal recycling.

Rare high-3He/4He signatures in ocean island basalts (OIB) erupted at volcanic hotspots derive from deep-seated domains preserved in Earth's interior. Only high-3He/4He OIB exhibit anomalous 182W-an isotopic signature inherited during the earliest history of Earth-supporting an ancient origin of high 3He/4He. However, it is not understood why some OIB host anomalous 182W while others do not. We provide geochemical data for the highest-3He/4He lavas from Iceland (up to 42.9 times atmospheric) with anomalous 182W and examine how Sr-Nd-Hf-Pb isotopic variations-useful for tracing subducted, recycled crust-relate to high 3He/4He and anomalous 182W. These data, together with data on global OIB, show that the highest-3He/4He and the largest-magnitude 182W anomalies are found only in geochemically depleted mantle domains-with high 143Nd/144Nd and low 206Pb/204Pb-lacking strong signatures of recycled materials. In contrast, OIB with the strongest signatures associated with recycled materials have low 3He/4He and lack anomalous 182W. These observations provide important clues regarding the survival of the ancient He and W signatures in Earth's mantle. We show that high-3He/4He mantle domains with anomalous 182W have low W and 4He concentrations compared to recycled materials and are therefore highly susceptible to being overprinted with low 3He/4He and normal (not anomalous) 182W characteristic of subducted crust. Thus, high 3He/4He and anomalous 182W are preserved exclusively in mantle domains least modified by recycled crust. This model places the long-term preservation of ancient high 3He/4He and anomalous 182W in the geodynamic context of crustal subduction and recycling and informs on survival of other early-formed heterogeneities in Earth's interior.

Melt inclusion constraints on petrogenesis of the 2014-2015 Holuhraun eruption, Iceland.

The 2014-2015 Holuhraun eruption, on the Bárðarbunga volcanic system in central Iceland, was one of the best-monitored basaltic fissure eruptions that has ever occurred, and presents a unique opportunity to link petrological and geochemical data with geophysical observations during a major rifting episode. We present major and trace element analyses of melt inclusions and matrix glasses from a suite of ten samples collected over the course of the Holuhraun eruption. The diversity of trace element ratios such as La/Yb in Holuhraun melt inclusions reveals that the magma evolved via concurrent mixing and crystallization of diverse primary melts in the mid-crust. Using olivine-plagioclase-augite-melt (OPAM) barometry, we calculate that the Holuhraun carrier melt equilibrated at 2.1 ± 0.7 kbar (7.5 ± 2.5 km), which is in agreement with the depths of earthquakes (6 ± 1 km) between Bárðarbunga central volcano and the eruption site in the days preceding eruption onset. Using the same approach, melt inclusions equilibrated at pressures between 0.5 and 8.0 kbar, with the most probable pressure being 3.2 kbar. Diffusion chronometry reveals minimum residence timescales of 1-12 days for melt inclusion-bearing macrocrysts in the Holuhraun carrier melt. By combining timescales of diffusive dehydration of melt inclusions with the calculated pressure of H2O saturation for the Holuhraun magma, we calculate indicative magma ascent rates of 0.12-0.29 m s-1. Our petrological and geochemical data are consistent with lateral magma transport from Bárðarbunga volcano to the eruption site in a shallow- to mid-crustal dyke, as has been suggested on the basis of seismic and geodetic datasets. This result is a significant step forward in reconciling petrological and geophysical interpretations of magma transport during volcano-tectonic episodes, and provides a critical framework for the interpretation of premonitory seismic and geodetic data in volcanically active regions.

A dynamic mid-crustal magma domain revealed by the 2023 to 2024 Sundhnúksgígar eruptions in Iceland.

Mid-crustal magma domains are the source of many basaltic eruptions. Lavas from individual eruptions are often chemically homogeneous, suggesting that they derive from single, well-mixed magma reservoirs. The 2023 to 2024 eruptions at Sundhnúksgígar in the Svartsengi volcanic system in Iceland provide an opportunity to observe the behavior of a mid-crustal magma domain at high spatial and temporal resolution by detailed sampling and geochemical characterization. We observed substantial mantle-derived geochemical variability in the products erupted in the first hours of the December 2023 and January, February, and March to May 2024 eruptions, indicating that the eruptions derived from multiple magma reservoirs, which mineral-melt equilibration pressures place in the mid crust. The unusual presence of geochemical heterogeneity in the mid-crustal magma domain provides insights into how dynamic and complex mid-crustal magma domains can be.

Ultrahigh-precision noble gas isotope analyses reveal pervasive subsurface fractionation in hydrothermal systems.

Mantle-derived noble gases in volcanic gases are powerful tracers of terrestrial volatile evolution, as they contain mixtures of both primordial (from Earth's accretion) and secondary (e.g., radiogenic) isotope signals that characterize the composition of deep Earth. However, volcanic gases emitted through subaerial hydrothermal systems also contain contributions from shallow reservoirs (groundwater, crust, atmosphere). Deconvolving deep and shallow source signals is critical for robust interpretations of mantle-derived signals. Here, we use a novel dynamic mass spectrometry technique to measure argon, krypton, and xenon isotopes in volcanic gas with ultrahigh precision. Data from Iceland, Germany, United States (Yellowstone, Salton Sea), Costa Rica, and Chile show that subsurface isotope fractionation within hydrothermal systems is a globally pervasive and previously unrecognized process causing substantial nonradiogenic Ar-Kr-Xe isotope variations. Quantitatively accounting for this process is vital for accurately interpreting mantle-derived volatile (e.g., noble gas and nitrogen) signals, with profound implications for our understanding of terrestrial volatile evolution.

Complex organic matter degradation by secondary consumers in chemolithoautotrophy-based subsurface geothermal ecosystems.

Microbial communities in terrestrial geothermal systems often contain chemolithoautotrophs with well-characterized distributions and metabolic capabilities. However, the extent to which organic matter produced by these chemolithoautotrophs supports heterotrophs remains largely unknown. Here we compared the abundance and activity of peptidases and carbohydrate active enzymes (CAZymes) that are predicted to be extracellular identified in metagenomic assemblies from 63 springs in the Central American and the Andean convergent margin (Argentinian backarc of the Central Volcanic Zone), as well as the plume-influenced spreading center in Iceland. All assemblies contain two orders of magnitude more peptidases than CAZymes, suggesting that the microorganisms more often use proteins for their carbon and/or nitrogen acquisition instead of complex sugars. The CAZy families in highest abundance are GH23 and CBM50, and the most abundant peptidase families are M23 and C26, all four of which degrade peptidoglycan found in bacterial cells. This implies that the heterotrophic community relies on autochthonous dead cell biomass, rather than allochthonous plant matter, for organic material. Enzymes involved in the degradation of cyanobacterial- and algal-derived compounds are in lower abundance at every site, with volcanic sites having more enzymes degrading cyanobacterial compounds and non-volcanic sites having more enzymes degrading algal compounds. Activity assays showed that many of these enzyme classes are active in these samples. High temperature sites (> 80°C) had similar extracellular carbon-degrading enzymes regardless of their province, suggesting a less well-developed population of secondary consumers at these sites, possibly connected with the limited extent of the subsurface biosphere in these high temperature sites. We conclude that in < 80°C springs, chemolithoautotrophic production supports heterotrophs capable of degrading a wide range of organic compounds that do not vary by geological province, even though the taxonomic and respiratory repertoire of chemolithoautotrophs and heterotrophs differ greatly across these regions.

Nitrogen Incorporation in Potassic and Micro- and Meso-Porous Minerals: Potential Biogeochemical Records and Targets for Mars Sampling.

We measured the N concentrations and isotopic compositions of 44 samples of terrestrial potassic and micro- and meso-porous minerals and a small number of whole-rocks to determine the extent to which N is incorporated and stored during weathering and low-temperature hydrothermal alteration in Mars surface/near-surface environments. The selection of these minerals and other materials was partly guided by the study of altered volcanic glass from Antarctica and Iceland, in which the incorporation of N as NH4+ in phyllosilicates is indicated by correlated concentrations of N and the LILEs (i.e., K, Ba, Rb, Cs), with scatter likely related to the presence of exchanged, occluded/trapped, or encapsulated organic/inorganic N occurring within structural cavities (e.g., in zeolites). The phyllosilicates, zeolites, and sulfates analyzed in this study contain between 0 and 99,120 ppm N and have δ15Nair values of -34‰ to +65‰. Most of these minerals, and the few siliceous hydrothermal deposits that were analyzed, have δ15N consistent with the incorporation of biologically processed N during low-temperature hydrothermal or weathering processes. Secondary ion mass spectrometry on altered hyaloclastites demonstrates the residency of N in smectites and zeolites, and silica. We suggest that geological materials known on Earth to incorporate and store N and known to be abundant at, or near, the surface of Mars should be considered targets for upcoming Mars sample return with the intent to identify any signs of ancient or modern life.

Linking deeply-sourced volatile emissions to plateau growth dynamics in southeastern Tibetan Plateau.

The episodic growth of high-elevation orogenic plateaux is controlled by a series of geodynamic processes. However, determining the underlying mechanisms that drive plateau growth dynamics over geological history and constraining the depths at which growth originates, remains challenging. Here we present He-CO2-N2 systematics of hydrothermal fluids that reveal the existence of a lithospheric-scale fault system in the southeastern Tibetan Plateau, whereby multi-stage plateau growth occurred in the geological past and continues to the present. He isotopes provide unambiguous evidence for the involvement of mantle-scale dynamics in lateral expansion and localized surface uplift of the Tibetan Plateau. The excellent correlation between 3He/4He values and strain rates, along the strike of Indian indentation into Asia, suggests non-uniform distribution of stresses between the plateau boundary and interior, which modulate southeastward growth of the Tibetan Plateau within the context of India-Asia convergence. Our results demonstrate that deeply-sourced volatile geochemistry can be used to constrain deep dynamic processes involved in orogenic plateau growth.

Unexpected large eruptions from buoyant magma bodies within viscoelastic crust.

Large volume effusive eruptions with relatively minor observed precursory signals are at odds with widely used models to interpret volcano deformation. Here we propose a new modelling framework that resolves this discrepancy by accounting for magma buoyancy, viscoelastic crustal properties, and sustained magma channels. At low magma accumulation rates, the stability of deep magma bodies is governed by the magma-host rock density contrast and the magma body thickness. During eruptions, inelastic processes including magma mush erosion and thermal effects, can form a sustained channel that supports magma flow, driven by the pressure difference between the magma body and surface vents. At failure onset, it may be difficult to forecast the final eruption volume; pressure in a magma body may drop well below the lithostatic load, create under-pressure and initiate a caldera collapse, despite only modest precursors.

Magma reservoir dynamics at Toba caldera, Indonesia, recorded by oxygen isotope zoning in quartz.

Quartz is a common phase in high-silica igneous rocks and is resistant to post-eruptive alteration, thus offering a reliable record of magmatic processes in silicic magma systems. Here we employ the 75 ka Toba super-eruption as a case study to show that quartz can resolve late-stage temporal changes in magmatic δ18O values. Overall, Toba quartz crystals exhibit comparatively high δ18O values, up to 10.2‰, due to magma residence within, and assimilation of, local granite basement. However, some 40% of the analysed quartz crystals display a decrease in δ18O values in outermost growth zones compared to their cores, with values as low as 6.7‰ (maximum ∆core-rim = 1.8‰). These lower values are consistent with the limited zircon record available for Toba, and the crystallisation history of Toba quartz traces an influx of a low-δ18O component into the magma reservoir just prior to eruption. Here we argue that this late-stage low-δ18O component is derived from hydrothermally-altered roof material. Our study demonstrates that quartz isotope stratigraphy can resolve magmatic events that may remain undetected by whole-rock or zircon isotope studies, and that assimilation of altered roof material may represent a viable eruption trigger in large Toba-style magmatic systems.

Gradual caldera collapse at Bárdarbunga volcano, Iceland, regulated by lateral magma outflow.

Large volcanic eruptions on Earth commonly occur with a collapse of the roof of a crustal magma reservoir, forming a caldera. Only a few such collapses occur per century, and the lack of detailed observations has obscured insight into the mechanical interplay between collapse and eruption. We use multiparameter geophysical and geochemical data to show that the 110-square-kilometer and 65-meter-deep collapse of Bárdarbunga caldera in 2014-2015 was initiated through withdrawal of magma, and lateral migration through a 48-kilometers-long dike, from a 12-kilometers deep reservoir. Interaction between the pressure exerted by the subsiding reservoir roof and the physical properties of the subsurface flow path explain the gradual, near-exponential decline of both collapse rate and the intensity of the 180-day-long eruption.

Evidence for primordial water in Earth's deep mantle.

The hydrogen-isotope [deuterium/hydrogen (D/H)] ratio of Earth can be used to constrain the origin of its water. However, the most accessible reservoir, Earth's oceans, may no longer represent the original (primordial) D/H ratio, owing to changes caused by water cycling between the surface and the interior. Thus, a reservoir completely isolated from surface processes is required to define Earth's original D/H signature. Here we present data for Baffin Island and Icelandic lavas, which suggest that the deep mantle has a low D/H ratio (δD more negative than -218 per mil). Such strongly negative values indicate the existence of a component within Earth's interior that inherited its D/H ratio directly from the protosolar nebula.