The isotopic composition of lavas associated with mantle plumes has previously been interpreted in the light of core–mantle interaction, suggesting that mantle plumes may transport core material to Earth’s surface 1–5 . However, a definitive fingerprint of Earth’s core in the mantle remains unconfirmed. Precious metals, such as ruthenium (Ru), are highly concentrated in the metallic core but extremely depleted in the silicate mantle. Recently discovered mass-independent Ru isotope variations ( ε 100 Ru) in ancient rocks show that the Ru isotope composition of accreted material changed during later stages of Earth’s growth 6 , indicating that the core and mantle must have different Ru isotope compositions. This illustrates the potential of Ru isotopes as a new tracer for core–mantle interaction. Here we report Ru isotope anomalies for ocean island basalts. Basalts from Hawaii have higher ε 100 Ru than the ambient mantle. Combined with unradiogenic tungsten (W) isotope ratios, this is diagnostic of a core contribution to their mantle sources. The combined Ru and W isotope systematics of Hawaiian basalts are best explained by simple core entrainment but addition of core-derived oxide minerals at the core–mantle boundary is a possibility.

The lavas associated with mantle plumes may sample domains throughout Earth’s mantle and probe its dynamics. However, plume studies are often only able to take snapshots in time, usually of the most recent plume activity, leaving the chemical and geodynamic evolution of major convective upwellings in Earth’s mantle poorly constrained. Here, we report the geodynamically key information of how the lithology and density of a plume change from plume head phase to tail. We use iron stable isotopes and thermodynamic modeling to show that the Galápagos plume has contained small, nearly constant, amounts of dense recycled crust over its 90-million-year history. Despite a temporal evolution in the amount of recycled crust-derived melt in Galápagos-related lavas, we show that this can be explained by plume cooling alone, without associated changes in the plume’s mantle source; results are also consistent with a plume rooted in a lower mantle low-velocity zone also sampling primordial components.

Recently reported major element, trace element, and volatile compositions of dacites erupted from Hekla Volcano, Iceland are inconsistent with an anatectic origin. In contrast, all of the data are consistent with dacitic melts being derived by modest amounts of crystallization of an anhydrous assemblage from an andesitic melt. Models that apply straightforward Gd, Yb, Th, and U concentrations, instead of ratios, to the problem show that the crystallization hypothesis remains valid, especially when uncertainties are propagated.

Using recent advancements in high-performance computing data assimilation to combine satellite InSAR data with numerical models, the prolonged unrest of the Sierra Negra volcano in the Galápagos was tracked to provide a fortuitous, but successful, forecast 5 months in advance of the 26 June 2018 eruption. Subsequent numerical simulations reveal that the evolution of the stress state in the host rock surrounding the Sierra Negra magma system likely controlled eruption timing. While changes in magma reservoir pressure remained modest (<15 MPa), modeled widespread Mohr-Coulomb failure is coincident with the timing of the 26 June 2018 moment magnitude 5.4 earthquake and subsequent eruption. Coulomb stress transfer models suggest that the faulting event triggered the 2018 eruption by encouraging tensile failure along the northern portion of the caldera. These findings provide a critical framework for understanding Sierra Negra’s eruption cycles and evaluating the potential and timing of future eruptions.

Hekla is an elongate volcano that lies at the intersection of the South Iceland Seismic Zone and the Eastern Volcanic Zone. We report major and trace element, oxygen isotopic, and H 2 O analyses on rocks, glass, melt inclusions, and minerals from almost all of the historical lavas and tephra deposits. This new dataset confirms the remarkable observation that not only are many eruptions compositionally zoned from felsic to mafic, but the extent of zoning relates directly to the length of repose since the previous eruption. Compositional data are consistent with the origin of the basaltic andesites and andesites by fractional crystallization, with no measurable crustal interaction once basaltic andesite has been produced. Although the 1104 CE Plinian rhyolite and 1158 CE dacite are also created by fractional crystallization, uranium–thorium isotopic disequilibria measured by others require that they evolved in a separate body, where magma is stored in a molten state for >10 4 years. Consistent trace element trends and ratios, as well as oxygenisotopic data, preclude significant crustal input into the evolving magma. The phenocryst assemblages are dominated by crystals that formed from their host melt; an exception is the 1158 CE dacite, which contains abundant crystals that formed from the 1104 CE rhyolite melt. A suite of thermobarometers indicates that most crystals formed in the lower crust at temperatures ranging from ∼1010 to 850 °C. Hekla's unique and systematic petrological time series and geophysical activity are attributed to the unusual geometry of the magma body, which we propose to be a tabular, vertically elongate macrodike, extending from the lower to the upper crust. The vertical body is recharged with basaltic andesite magma at the end of each eruption, which then undergoes cooling and crystallization until the subsequent eruption. The entire system is supplied by a lower-crustal body of basaltic andesite, which is produced by fractional crystallization of basaltic magma in a reservoir that is thermochemically buffered to ∼1010 °C. Cooling and crystallization of recharged basaltic andesite magma in a background geothermal gradient from the lower to the shallow crust accounts for the systematic relationship between repose and composition.

Many volcanoes erupt compositionally homogeneous magmas over timescales ranging from decades to millennia. This monotonous activity is thought to reflect a high degree of chemical homogeneity in their magmatic systems, leading to predictable eruptive behaviour. We combine petrological analyses of erupted crystals with new thermodynamic models to characterise the diversity of melts in magmatic systems beneath monotonous shield volcanoes in the Galápagos Archipelago (Wolf and Fernandina). In contrast with the uniform basaltic magmas erupted at the surface over long timescales, we find that the sub-volcanic systems contain extreme heterogeneity, with melts extending to rhyolitic compositions. Evolved melts are in low abundance and large volumes of basalt flushing through the crust from depth overprint their chemical signatures. This process will only maintain monotonous activity while the volume of melt entering the crust is high, raising the possibility of transitions to more silicic activity given a decrease in the crustal melt flux.

Ensemble based data assimilation approaches, such as the Ensemble Kalman Filter (EnKF), have been widely and successfully implemented to combine observations with dynamic models to investigate the evolution of a system’s state. Such inversions are powerful tools for providing forecasts as well as “hindcasting” events such as volcanic eruptions to investigate source parameters and triggering mechanisms. In this study, a high performance computing (HPC) adaptation of the EnKF is used to assimilate ground deformation observations from interferometric synthetic-aperture radar (InSAR) into high-fidelity, multiphysics finite element models to evaluate the prolonged unrest and June 26, 2018 eruption of Sierra Negra volcano, Galápagos. The stability of the Sierra Negra magma system is evaluated at each time step by estimating variations in reservoir overpressure, Mohr-Coulomb failure in the host rock, and tensile stress and failure along the reservoir boundary. The deformation of Sierra Negra is tracked over a decade, during which almost 5 meters of surface uplift has been recorded. The EnKF reveals that the evolution of the stress state in the host rock surrounding the Sierra Negra magma reservoir likely controlled the timing of the eruption. While increases in magma reservoir overpressure remained modest (< 10 MPa) throughout the data assimilation time period, significant Mohr-Coulomb failure is indicated in the lead up to the eruption coincident with increased seismicity along both trapdoor faults within Sierra Negra’s caldera and along the caldera’s ring faults. During the final stages of pre-eruptive unrest, the EnKF models indicate limited tensile failure, with no tensile failure along the northern portion of the magma system where the eruption commenced. Most strikingly, model calculations of significant through-going Mohr-Coulomb failure correspond in space and time with a Mw 5.4 earthquake recorded in the hours preceding the 2018 eruption. Subsequent stress modeling implicates the Mw 5.4 earthquake along the southern intra-caldera trapdoor fault as the potential catalyst for tensile failure and dike initiation along the reservoir to the north. In conclusion, the volcano EnKF approach successfully tracked the evolving stability of Sierra Negra, indicating great potential for future forecasting efforts.

Ensemble based data assimilation approaches, such as the Ensemble Kalman Filter (EnKF), have been widely and successfully implemented to combine observations with dynamic models to investigate the evolution of a system’s state. Such inversions are powerful tools for providing forecasts as well as “hindcasting” events such as volcanic eruptions to investigate source parameters and triggering mechanisms. In this study, a high performance computing (HPC) adaptation of the EnKF is used to assimilate ground deformation observations from interferometric synthetic-aperture radar (InSAR) into high-fidelity, multiphysics finite element models to evaluate the prolonged unrest and June 26, 2018 eruption of Sierra Negra volcano, Galápagos. The stability of the Sierra Negra magma system is evaluated at each time step by estimating variations in reservoir overpressure, Mohr-Coulomb failure in the host rock, and tensile stress and failure along the reservoir boundary. The deformation of Sierra Negra is tracked over a decade, during which almost 5 meters of surface uplift has been recorded. The EnKF reveals that the evolution of the stress state in the host rock surrounding the Sierra Negra magma reservoir likely controlled the timing of the eruption. While increases in magma reservoir overpressure remained modest (< 10 MPa) throughout the data assimilation time period, significant Mohr-Coulomb failure is indicated in the lead up to the eruption coincident with increased seismicity along both trapdoor faults within Sierra Negra’s caldera and along the caldera’s ring faults. During the final stages of pre-eruptive unrest, the EnKF models indicate limited tensile failure, with no tensile failure along the northern portion of the magma system where the eruption commenced. Most strikingly, model calculations of significant through-going Mohr-Coulomb failure correspond in space and time with a Mw 5.4 earthquake recorded in the hours preceding the 2018 eruption. Subsequent stress modeling implicates the Mw 5.4 earthquake along the southern intra-caldera trapdoor fault as the potential catalyst for tensile failure and dike initiation along the reservoir to the north. In conclusion, the volcano EnKF approach successfully tracked the evolving stability of Sierra Negra, indicating great potential for future forecasting efforts.

This former West Point graduate and coffee grower transformed igneous petrology and volcanology.

There are clear microstructural differences between mafic plutonic rocks that formed in a dynamic liquid-rich environment, in which crystals can be moved and re-arranged by magmatic currents, and those in which crystal nucleation and growth are essentially in situ and static. Crystalline enclaves, derived from deep crustal mushy zones and erupted in many volcanic settings, afford a unique opportunity to use the understanding of microstructural development, established from the study of intrusive plutons, to place constraints on the architecture of sub-volcanic systems. Here, we review the relevant microstructural literature, before applying these techniques to interrogate the crystallization environments of enclaves from the Kameni Islands of Santorini and Rábida Volcano in the Galápagos. Crystals in samples of deep-sourced material from both case studies preserve evidence of at least some time spent in a liquid-rich environment. The Kameni enclaves appear to record an early stage of crystallization during which crystals were free to move, with the bulk of crystallization occurring in a static, mushy environment. By contrast, the Rábida enclaves were sourced from an environment in which hydrodynamic sorting and re-arrangement by magmatic currents were common, consistent with a liquid-rich magma chamber. While presently active volcanoes are thought to be underlain by extensive regions rich in crystal mush, these examples preserve robust evidence for the presence of liquid-rich magma chambers in the geological record. This article is part of the Theo Murphy meeting issue ‘Magma reservoir architecture and dynamics'.

Complex magmatic systems beneath monotonous basaltic shield volcanoes in the Galápagos Archipelago MICHAEL J. STOCK1, DENNIS GEIST2, DAVID A. NEAVE3, MARCO BAGNARDI4, JOHN MACLENNAN1, BENJAMIN BERNARD5, IRIS BUISMAN1, MATTHEW L.M. GLEESON1 1Department of Earth Sciences, University of Cambridge, UK, ms2368@cam.ac.uk 2Department of Geology, Colgate University, USA 3School of Earth and Environmental Science, University of Manchester, UK 4JPL, California Institute of Technology, USA 5Instituto Geofísico, Escuela Politécnica Nacional, Ecuador

Extensive vertical deformation (>4.5 m) observed at Sierra Negra volcano Galapagos, Ecuador, between 1992 and the 2005 eruption led scientists to hypothesize that repeated faulting events relieved magma chamber overpressure and prevented eruption. To better understand the catalyst of the 2005 eruption, thermomechanical models are used to track the stress state and stability of the magma storage system during the 1992-2005 inflation events. Numerical experiments indicate that the host rock surrounding the Sierra Negra reservoir remained in compression with minimal changes in overpressure (similar to 10 MPa) leading up to the 2005 eruption. The lack of tensile failure and minimal overpressure accumulation likely inhibited dike initiation and accommodated the significant inflation without the need for pressure relief through shallow trapdoor faulting events. The models indicate that static stress transfer due to the M-w 5.4 earthquake 3 hr prior to the eruption most likely triggered tensile failure and catalyzed the 2005 eruption. Plain Language Summary Tracking the stability of a magma system in the lead up to a volcanic eruption requires investigating both the pressure state of the magma reservoir and stress accumulation in the host rock. New coupled conduit flow-magma reservoir pressurization models are used to evaluate the evolution of the magma reservoir of Sierra Negra volcano, Galapagos, in the lead up to its 2005 eruption. Stress calculations indicate that the magma reservoir was stable prior to the 2005 eruption and that the eruption was likely triggered by a M-w 5.4 earthquake that occurred 3hours prior to the event. The new modeling approach has important implications for tracking the stress evolution of magma systems to evaluate future unrest and eruption triggering mechanisms at volcanoes worldwide.

The hydromagmatic eruption that immediately preceded the 1968 caldera collapse of Fernandina Volcano, Galapagos, which had a volcano explosivity index (VEI) of 4, offers a case study of powerful eruptions where basaltic magma interacts with caldera-ponded water. The 4-d-long hydromagmatic eruption sequence records an early stage and a small fraction of the volume of magmatic withdrawal that led the caldera floor to lower 350 m over the next 10 d. Erupted tephra was lithic-rich. The small proportion of juvenile basaltic glass included blocky fragments, Pele's tears, and Pele's hair. Pyroclastic density currents swept across the western summit plateau 600-700 m above the vent and deposited dunes, cross-bedded and rubbly breccia deposits, imbricated lag blocks, and ash plasters, and toppled trees. Blocks ejected out of the caldera formed impact craters on the volcano's flank >600 m higher and >1 km away. Ejected blocks are mostly basalt but include cumulate olivine gabbro. The vent area enlarged by 300 x 10(6) m(3) during the eruption. A small adjacent fault-bounded block subsided after the eruption. Lake water and groundwater confined within the caldera by ring dikes were available to interact with hot rocks and magma. In our interpretation, this water helped to trigger and feed the eruption by interacting with rocks above a lowering magma column. Ecosystems recovered rapidly on the tephra. Eruptions have not diminished the island's biodiversity despite Fernandina's high rate of volcanic activity, including the massive resurfacing in 1968. Stratigraphic evidence suggests that the 1968 eruption may be only the latest in a series of explosive eruptions from the caldera.

The 2015 eruption of Wolf volcano was one of the largest eruptions in the Galápagos Islands since the onset of routine satellite‐based volcano monitoring. It therefore provides an excellent opportunity to combine geophysical and petrological data, to place detailed constraints on the architecture and dynamics of subvolcanic systems in the western archipelago. We present new geodetic models that show that pre‐eruptive inflation at Wolf was caused by magma accumulation in a shallow flat‐topped reservoir at ~1.1 km, whereas edifice‐scale deformation during the eruption was related to a deflationary source at 6.1–8.8 km. Petrological observations suggest that the erupted material was derived from both a subvolcanic mush and a liquid‐rich magma body. Using a combination of olivine‐plagioclase‐augite‐melt (OPAM) and clinopyroxene‐melt barometry, we show that the majority of magma equilibration, crystallization, and mush entrainment occurred at a depth equal to or greater than the deep geodetic source, with little petrological evidence of material sourced from shallower levels. Hence, our multidisciplinary study does not support a fully transcrustal magmatic system beneath Wolf volcano before the 2015 eruption but instead indicates two discrete storage regions, with a small magma lens at shallow levels and the major zone of magma storage in the lower crust, from which most of the erupted material was sourced. A predominance of lower crustal magma storage has previously been thought typical of subvolcanic systems in the eastern Galápagos Archipelago, but our new data suggest that this may also occur beneath the more active volcanoes of the western archipelago.

The older eastern Galapagos are different in almost every way from the historically active western Galapagos volcanoes. The western Galapagos volcanoes have steep upper slopes and are topped by large calderas, whereas none of the older islands has a caldera, an observation that is supported by recent gravity measurements. Moreover, the eastern islands tend to have been constructed by linear fissure systems and many are cut by faults. Most of the western volcanoes erupt evolved basalts with an exceedingly small range of Mg#, Lan/Smn, and Smn/Ybn. This is attributed to homogenization in a crustal-scale magmatic mush column, which is maintained in a thermochemical steady state, owing to high magma supply directly over the Galapagos mantle plume. The exceptions are volcanoes at the leading edge of the hotspot, which have yet to develop mush columns, and volcanoes that are waning in activity, because they are being carried away from the plume. In contrast, the eastern volcanoes erupt relatively primitive magmas, with a large range in Mg#, Lan/Smn, and Smn/Ybn. This is attributed to isolated, ephemeral magmatic plumbing systems supplied by smaller magmatic fluxes throughout their histories. Consequently, each batch of magma follows an independent course of evolution, owing to the low volume of hypersolidus material beneath these volcanoes. The magmatic flux to Galapagos volcanoes negatively correlates with the distance to the Galapagos Spreading Center (GSC). When the ridge was close to the plume, most of the plume-derived magma was directed to the ridge. Currently, the active volcanoes are much farther from the GSC, thus most of the plume-derived magma erupts on the Nazca Plate and can be focused beneath the large young shields. We define an intermediate sub-province comprising Rabida, Santiago and Pinzon volcanoes, which were most active about 1 Ma. They have all erupted dacites, rhyolites, and trachytes, similar to the dying stage of the western volcanoes, indicating that there was a relatively large volume of mush beneath them. Morphologically, however, they are more like the eastern volcanoes, and have erupted lavas with a large range in composition.