The basalts of the 2021 Fagradalsfjall eruption were the first erupted on the Reykjanes Peninsula in 781 years and offer a unique opportunity to determine the composition of the mantle underlying Iceland, in particular its oxygen isotope composition (δ¹⁸O values). The basalts show compositional variations in Zr/Y, Nb/Zr and Nb/Y values that span roughly half of the previously described range for Icelandic basaltic magmas and signal involvement of Icelandic plume (OIB) and Enriched Mid-Ocean Ridge Basalt (EMORB) in magma genesis. Here we show that Fagradalsfjall δ¹⁸O values are invariable (mean δ¹⁸O = 5.4 ± 0.3‰ 2 SD, N = 47) and indistinguishable from “normal” upper mantle, in contrast to significantly lower δ¹⁸O values reported for erupted materials elsewhere in Iceland (e.g., the 2014–2015 eruption at Holuhraun, Central Iceland). Thus, despite differing trace element characteristics, the melts that supplied the Fagradalsfjall eruption show no evidence for ¹⁸O-depleted mantle or interaction with low-δ¹⁸O crust and may therefore represent a useful mantle reference value in this part of the Icelandic plume system.

The Eldgjá eruption is the largest basalt lava flood of the Common Era. It has been linked to a major ice‐core sulfur (S) spike in 939–940 CE and Northern Hemisphere summer cooling in 940 CE. Despite its magnitude and potential climate impacts, uncertainties remain concerning the eruption timeline, atmospheric dispersal of emitted volatiles, and coincident volcanism in Iceland and elsewhere. Here, we present a comprehensive analysis of Greenland ice‐cores from 936 to 943 CE, revealing a complex volatile record and cryptotephra with numerous geochemical populations. Transitional alkali basalt tephra matching Eldgjá are found in 939–940 CE, while tholeiitic basalt shards present in 936/937 CE and 940/941 CE are compatible with contemporaneous Icelandic eruptions from Grímsvötn and Bárðarbunga‐Veiðivötn systems (including V‐Sv tephra). We also find four silicic tephra populations, one of which we link to the Jala Pumice of Ceboruco (Mexico) at 941 ± 1 CE. Triple S isotopes, Δ³³S, spanning 936–940 CE are indicative of upper tropospheric/lower stratospheric transport of aerosol sourced from the Icelandic fissure eruptions. However, anomalous Δ³³S (down to −0.4‰) in 940–941 CE evidence stratospheric aerosol transport consistent with summer surface cooling revealed by tree‐ring reconstructions. Tephra associated with the anomalous Δ³³S have a variety of compositions, complicating the attribution of climate cooling to Eldgjá alone. Nevertheless, our study confirms a major S emission from Eldgjá in 939–940 CE and implicates Eldgjá and a cluster of eruptions as triggers of summer cooling, severe winters, and privations in ∼940 CE.

The purpose of the EUNADICS-AV (European Natural Airborne Disaster Information and Coordination System for Aviation) prototype early warning system (EWS) is to develop the combined use of harmonised data products from satellite, ground-based and in situ instruments to produce alerts of airborne hazards (volcanic, dust, smoke and radionuclide clouds), satisfying the requirement of aviation air traffic management (ATM) stakeholders (https://cordis.europa.eu/project/id/723986, last access: 5 November 2021). The alert products developed by the EUNADICS-AV EWS, i.e. near-real-time (NRT) observations, email notifications and netCDF (Network Common Data Form) alert data products (called NCAP files), have shown significant interest in using selective detection of natural airborne hazards from polar-orbiting satellites. The combination of several sensors inside a single global system demonstrates the advantage of using a triggered approach to obtain selective detection from observations, which cannot initially discriminate the different aerosol types. Satellite products from hyperspectral ultraviolet–visible (UV–vis) and infrared (IR) sensors (e.g. TROPOMI – TROPOspheric Monitoring Instrument – and IASI – Infrared Atmospheric Sounding Interferometer) and a broadband geostationary imager (Spinning Enhanced Visible and InfraRed Imager; SEVIRI) and retrievals from ground-based networks (e.g. EARLINET – European Aerosol Research Lidar Network, E-PROFILE and the regional network from volcano observatories) are combined by our system to create tailored alert products (e.g. selective ash detection, SO₂ column and plume height, dust cloud, and smoke from wildfires). A total of 23 different alert products are implemented, using 1 geostationary and 13 polar-orbiting satellite platforms, 3 external existing service, and 2 EU and 2 regional ground-based networks. This allows for the identification and the tracking of extreme events. The EUNADICS-AV EWS has also shown the need to implement a future relay of radiological data (gamma dose rate and radionuclides concentrations in ground-level air) in the case of a nuclear accident. This highlights the interest of operating early warnings with the use of a homogenised dataset. For the four types of airborne hazard, the EUNADICS-AV EWS has demonstrated its capability to provide NRT alert data products to trigger data assimilation and dispersion modelling providing forecasts and inverse modelling for source term estimate. Not all of our alert data products (NCAP files) are publicly disseminated. Access to our alert products is currently restricted to key users (i.e. Volcanic Ash Advisory Centres, national meteorological services, the World Meteorological Organization, governments, volcano observatories and research collaborators), as these are considered pre-decisional products. On the other hand, thanks to the EUNADICS-AV–SACS (Support to Aviation Control Service) web interface (https://sacs.aeronomie.be, last access: 5 November 2021), the main part of the satellite observations used by the EUNADICS-AV EWS is shown in NRT, with public email notification of volcanic emission and delivery of tailored images and NCAP files. All of the ATM stakeholders (e.g. pilots, airlines and passengers) can access these alert products through this free channel.

The 10th century Eldgjá fissure eruption is the largest in Iceland in historical time. It erupted 21.0 km3 of magma, with 1.3 km3 as tephra in at least 16 explosive episodes from subaerial and subglacial vents, producing magmatic and phreatomagmatic deposits respectively. Grain-size distributions for these end-members show distinct differences at comparable distances from source: the former are coarser and unimodal; the latter are finer and bimodal. These distributions appear to record different primary fragmentation histories. In contrast, the vesicle-size distributions of pyroclasts from each type of deposit show the pyroclasts underwent similar vesicle nucleation and growth prior to fragmentation. This indicates that the role of glacial water was comparatively late-stage, re-fragmenting an already disrupting magma by quench granulation. The presence of microlite-rich domains within clasts reveals a history of complex conduit evolution, during the transition from a continuous dyke to focussed, discrete vents.

Grabens, or valleys formed during extensional tectonic events, are common but rarely observed during formation. In November 2023, inelastic surface deformation formed abruptly along Iceland's plate boundary in Grindavík. We documented graben formation in real‐time through satellite mapping (InSAR), seismicity, GNSS data, repeated lidar surveys, and field mapping. Five normal faults and ∼12 fissures ruptured the surface delineating two grabens separated by a horst, a context not present in other contemporary case studies. The graben normal faults slipped rapidly (over hours) and maximum surface motions coincided with the occurrence of turbulent seismic swarms in both space and time. Although 3 eruptions took place ∼15 km northeast of Grindavík from 2021 to 2023, attributed to magma intrusions (i.e., dikes), none of these also formed grabens. Thus, the Grindavík grabens shows evidence for tectonic origins. Real‐time monitoring of these phenomena provide insight into graben formation on Earth and potentially on other planets.

The interaction of magma with external water commonly enhances magma fragmentation through the conversion of thermal to mechanical energy and results in an increased production of fine‐grained volcanic tephra. Magma‐water interaction is thus of importance for hazard mitigation on both a local and a regional scales. The relative proportion of water that interacts with magma, quantified as the water‐to‐melt mass ratio, is thought to determine the efficiency of thermal to mechanical energy conversion, termed the fragmentation efficiency. Here, we analyze the pyroclast size distributions from the 10th century Eldgjá fissure eruption in Iceland, where parts of the fissure erupted subglacially and other erupted subaerially. The subglacially erupted magma passed through a column of glacial meltwater, resulting in a larger proportion of finer pyroclast sizes relative to the subaerially erupted, purely magmatic tephra. This finer grain size distribution has been attributed to quench‐granulation induced by enhanced cooling upon interaction with external water. We hypothesize that the additional fragmentation (surface) energy required to produce the finer grained hydromagmatic deposits is due to the conversion of thermal to mechanical energy associated with the entrainment of water into the volcanic jet, as it passed through a column of subglacial melt water. Based on field and granulometry data, we estimate that the interaction of the volcanic jet with the meltwater provided an additional fragmentation energy of approximately 3–14 kJ per kg of pyroclasts. We numerically model the hydrofragmentation energy within a jet that passes through a layer of meltwater. We find that the water‐to‐melt mass ratio of entrained water required to produce the additional fragmentation energy is in the range of 1–2, which requires a minimum ice melting rate of 10⁴ m³ s⁻¹. Our simulation results show that the water‐to‐melt ratio is an important parameter that controls the ascent of plume in the atmosphere.