Deep magma storage during the 2021 La Palma eruption.

The 2021 La Palma eruption provided an unpreceded opportunity to test the relationship between earthquake hypocenters and the location of magma reservoirs. We performed density measurements on CO2-rich fluid inclusions (FIs) hosted in olivine crystals that are highly sensitive to pressure via calibrated Raman spectroscopy. This technique can revolutionize our knowledge of magma storage and transport during an ongoing eruption, given that it can produce precise magma storage depth constraints in near real time with minimal sample preparation. Our FIs have CO2 recorded densities from 0.73 to 0.98 g/cm3, translating into depths of 15 to 27 km, which falls within the reported deep seismic zone recording the main melt storage reservoir.

Hidden mechanical weaknesses within lava domes provided by buried high-porosity hydrothermal alteration zones.

Catastrophic lava dome collapse is considered an unpredictable volcanic hazard because the physical properties, stress conditions, and internal structure of lava domes are not well understood and can change rapidly through time. To explain the locations of dome instabilities at Merapi volcano, Indonesia, we combined geochemical and mineralogical analyses, rock physical property measurements, drone-based photogrammetry, and geoinformatics. We show that a horseshoe-shaped alteration zone that formed in 2014 was subsequently buried by renewed lava extrusion in 2018. Drone data, as well as geomechanical, mineralogical, and oxygen isotope data suggest that this zone is characterized by high-porosity hydrothermally altered materials that are mechanically weak. We additionally show that the new lava dome is currently collapsing along this now-hidden weak alteration zone, highlighting that a detailed understanding of dome architecture, made possible using the monitoring techniques employed here, is essential for assessing hazards associated with dome and edifice failure at volcanoes worldwide.

Sunda arc mantle source δ18O value revealed by intracrystal isotope analysis.

Magma plumbing systems underlying subduction zone volcanoes extend from the mantle through the overlying crust and facilitate protracted fractional crystallisation, assimilation, and mixing, which frequently obscures a clear view of mantle source compositions. In order to see through this crustal noise, we present intracrystal Secondary Ion Mass Spectrometry (SIMS) δ18O values in clinopyroxene from Merapi, Kelut, Batur, and Agung volcanoes in the Sunda arc, Indonesia, under which the thickness of the crust decreases from ca. 30 km at Merapi to ≤20 km at Agung. Here we show that mean clinopyroxene δ18O values decrease concomitantly with crustal thickness and that lavas from Agung possess mantle-like He-Sr-Nd-Pb isotope ratios and clinopyroxene mean equilibrium melt δ18O values of 5.7 ‰ (±0.2 1 SD) indistinguishable from the δ18O range for Mid Ocean Ridge Basalt (MORB). The oxygen isotope composition of the mantle underlying the East Sunda Arc is therefore largely unaffected by subduction-driven metasomatism and may thus represent a sediment-poor arc end-member.

Constraining the sub-arc, parental magma composition for the giant Altiplano-Puna Volcanic Complex, northern Chile.

The Andean continental arc is built upon the thickest crust on Earth, whose eruption products reflect varying degrees of crustal assimilation. In order to robustly model magma evolution and assimilation at subduction zones such as the Andes, the compositions of parental magmas feeding crustal magma reservoirs need to be defined. Here we present new olivine and clinopyroxene oxygen isotope data from rare mafic volcanic rocks erupted at the margins of the giant Altiplano-Puna Magma Body (APMB) of the Altiplano-Puna Volcanic Complex, Central Andes. Existing olivine and pyroxene δ18O values for the Central Andes are highly variable and potentially not representative of sub-arc parental compositions. However, new olivine (n = 6) and clinopyroxene (n = 12) δ18O values of six Central Andean volcanoes presented here display a narrow range, with averages at 6.0‰ ± 0.2 (2σ S.D.) and 6.7‰ ± 0.3 (2σ S.D.), consistent with a common history for the investigated minerals. These data allow us to estimate the δ18O values of sub-arc, parental melts to ca. 7.0‰ ± 0.2 (2σ S.D.). Parental melts feeding the APMB and associated volcanic centres are postulated to form in the felsic continental crust following assimilation of up to 28% high-δ18O basement rocks by mantle-derived magmas.

Hydrothermal alteration of andesitic lava domes can lead to explosive volcanic behaviour.

Dome-forming volcanoes are among the most hazardous volcanoes on Earth. Magmatic outgassing can be hindered if the permeability of a lava dome is reduced, promoting pore pressure augmentation and explosive behaviour. Laboratory data show that acid-sulphate alteration, common to volcanoes worldwide, can reduce the permeability on the sample lengthscale by up to four orders of magnitude and is the result of pore- and microfracture-filling mineral precipitation. Calculations using these data demonstrate that intense alteration can reduce the equivalent permeability of a dome by two orders of magnitude, which we show using numerical modelling to be sufficient to increase pore pressure. The fragmentation criterion shows that the predicted pore pressure increase is capable of fragmenting the majority of dome-forming materials, thus promoting explosive volcanism. It is crucial that hydrothermal alteration, which develops over months to years, is monitored at dome-forming volcanoes and is incorporated into real-time hazard assessments.

Crustal CO2 contribution to subduction zone degassing recorded through calc-silicate xenoliths in arc lavas.

Interaction between magma and crustal carbonate at active arc volcanoes has recently been proposed as a source of atmospheric CO2, in addition to CO2 released from the mantle and subducted oceanic crust. However, quantitative constraints on efficiency and timing of these processes are poorly established. Here, we present the first in situ carbon and oxygen isotope data of texturally distinct calcite in calc-silicate xenoliths from arc volcanics in a case study from Merapi volcano (Indonesia). Textures and C-O isotopic data provide unique evidence for decarbonation, magma-fluid interaction, and the generation of carbonate melts. We report extremely light δ13CPDB values down to -29.3‰ which are among the lowest reported in magmatic systems so far. Combined with the general paucity of relict calcite, these extremely low values demonstrate highly efficient remobilisation of crustal CO2 over geologically short timescales of thousands of years or less. This rapid release of large volumes of crustal CO2 may impact global carbon cycling.

Multi-level magma plumbing at Agung and Batur volcanoes increases risk of hazardous eruptions.

The island of Bali in Indonesia is home to two active stratovolcanoes, Agung and Batur, but relatively little is known of their underlying magma plumbing systems. Here we define magma storage depths and isotopic evolution of the 1963 and 1974 eruptions using mineral-melt equilibrium thermobarometry and oxygen and helium isotopes in mineral separates. Olivine crystallised from a primitive magma and has average δ18O values of 4.8‰. Clinopyroxene records magma storage at the crust-mantle boundary, and displays mantle-like isotope values for Helium (8.62 RA) and δ18O (5.0-5.8‰). Plagioclase reveals crystallisation in upper crustal storage reservoirs and shows δ18O values of 5.5-6.4‰. Our new thermobarometry and isotope data thus corroborate earlier seismic and InSAR studies that inferred upper crustal magma storage in the region. This type of multi-level plumbing architecture could drive replenishing magma to rapid volatile saturation, thus increasing the likelihood of explosive eruptions and the consequent hazard potential for the population of Bali.

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.

Boron isotope fractionation in magma via crustal carbonate dissolution.

Carbon dioxide released by arc volcanoes is widely considered to originate from the mantle and from subducted sediments. Fluids released from upper arc carbonates, however, have recently been proposed to help modulate arc CO2 fluxes. Here we use boron as a tracer, which substitutes for carbon in limestone, to further investigate crustal carbonate degassing in volcanic arcs. We performed laboratory experiments replicating limestone assimilation into magma at crustal pressure-temperature conditions and analysed boron isotope ratios in the resulting experimental glasses. Limestone dissolution and assimilation generates CaO-enriched glass near the reaction site and a CO2-dominated vapour phase. The CaO-rich glasses have extremely low δ(11)B values down to -41.5‰, reflecting preferential partitioning of (10)B into the assimilating melt. Loss of (11)B from the reaction site occurs via the CO2 vapour phase generated during carbonate dissolution, which transports (11)B away from the reaction site as a boron-rich fluid phase. Our results demonstrate the efficacy of boron isotope fractionation during crustal carbonate assimilation and suggest that low δ(11)B melt values in arc magmas could flag shallow-level additions to the subduction cycle.

Nannofossils in 2011 El Hierro eruptive products reinstate plume model for Canary Islands.

The origin and life cycle of ocean islands have been debated since the early days of Geology. In the case of the Canary archipelago, its proximity to the Atlas orogen led to initial fracture-controlled models for island genesis, while later workers cited a Miocene-Quaternary east-west age-progression to support an underlying mantle-plume. The recent discovery of submarine Cretaceous volcanic rocks near the westernmost island of El Hierro now questions this systematic age-progression within the archipelago. If a mantle-plume is indeed responsible for the Canaries, the onshore volcanic age-progression should be complemented by progressively younger pre-island sedimentary strata towards the west, however, direct age constraints for the westernmost pre-island sediments are lacking. Here we report on new age data obtained from calcareous nannofossils in sedimentary xenoliths erupted during the 2011 El Hierro events, which date the sub-island sedimentary rocks to between late Cretaceous and Pliocene in age. This age-range includes substantially younger pre-volcanic sedimentary rocks than the Jurassic to Miocene strata known from the older eastern islands and now reinstate the mantle-plume hypothesis as the most plausible explanation for Canary volcanism. The recently discovered Cretaceous submarine volcanic rocks in the region are, in turn, part of an older, fracture-related tectonic episode.