Holuhraun Eruption Research Articles

A New Calibration of the OPAM Thermobarometer for Anhydrous and Hydrous Mafic Systems

Abstract Melt-based thermobarometers are essential tools to recover pre-eruptive magma storage conditions through their application to bulk rock and liquid chemistry. In active volcanic systems, thermobarometric results can be combined with independent geophysical data during or after an eruption to validate conceptual models. In this contribution, we revisit the thermobarometer for melts equilibrated with the mineral assemblage of olivine + plagioclase + augitic clinopyroxene (OPAM). We first demonstrate that the most widely applied OPAM thermobarometer suffers from both random and systematic uncertainty, even for anhydrous melts, and that the uncertainty increases proportionally with melt H2O. To address this issue, we use a modern compilation of anhydrous and hydrous OPAM-saturated experiments to regress a new empirical melt-based OPAM thermometer and barometer. Our new equations recover a validation dataset with a standard error estimate of ±1.14 kbar and ± 36 °C for pressure and temperature, respectively, as well as a low systematic uncertainty that does not depend on melt H2O. Additionally, we present a novel statistical approach to determine the probability that a given melt is OPAM-saturated, which can be used alongside rigorous petrographic and geochemical observations. Our thermobarometer and saturation test are presented as a user-friendly R script that reads from an input csv file to be populated with natural data. We benchmark the new calibrations on the products of the 2015 eruption of Wolf Volcano (Isabela Island, Galápagos archipelago) and the 2014–2015 Holuhraun eruption (Iceland), both of which have independent geophysical estimates of magma storage that agree well with our thermobarometric results.

Sulfide saturation and resorption modulates sulfur and metal availability during the 2014–15 Holuhraun eruption, Iceland

Mafic magmas may experience multiple stages of sulfide saturation and resorption during ascent and differentiation. Quenched tephra erupted during the 2014–15 Holuhraun eruption preserve abundant evidence for sulfide resorption, offering a rare opportunity to explore the sulfide life cycle from nucleation to resorption. Specifically, we combine detailed textural and chemical analyses of sulfides and silicate melts with geochemical models of sulfide saturation and degassing. This integrative approach demonstrates that sulfides began nucleating in melts with ~8 wt% MgO, persisted during fractionation to 6.5 wt% MgO, before resorbing heterogeneously in response to sulfur degassing. Sulfides are preserved preferentially in confined geometries within and between crystals, suggesting that kinetic effects impeded sulfur loss from the melt and maintained local sulfide saturation on eruption. The proportion of sulfides exhibiting breakdown textures increases throughout the eruption, coincident with decreasing magma discharge, indicating that sulfide resorption and degassing are kinetically limited. Sulfides likely modulate the emission of sulfur and chalcophile elements to the atmosphere and surface environment, with implications for assessing the environmental impacts and societal hazards of basaltic fissure eruptions.

SO2 emission rates and incorporation into the air pollution dispersion forecast during the 2021 eruption of Fagradalsfjall, Iceland

During the low-effusion rate Fagradalsfjall eruption (19 March – 18 September 2021), the emission of sulfur dioxide (SO2) was frequently measured using ground-based UV spectrometers. The total SO2 emitted during the entire eruption was 970 ± 540 kt, which is only about 6% of the SO2 emitted during the similar length Holuhraun eruption (2014–2015). The eruption was divided into five phases based on visual observations, including the number of active vents and the occurrence of lava fountaining. The SO2 emission rate ranged from 44 ± 19 kg/s in Phase 2 to 85 ± 29 kg/s in Phase 5, with an average of 64 ± 34 kg/s for the entire eruption. There was notable variability in SO2 on short timescales, with measurements on 11 August 2021 ranging from 17 to 78 kg/s.SO2 flux measurements were made using scanning DOAS instruments located at different distances from and orientations relative to the eruption site augmented by traverses. Four hundred and forty-four scan and traverse measurements met quality criteria and were used, along with plume height and meteorological data, to calculate SO2 fluxes while accounting for wind-related uncertainties.A tendency for stronger SO2 flux concurrent with higher amplitude seismic tremor and the occurrence of lava fountaining was observed during Phases 4 and 5 which were characterized by intermittent crater activity including observable effusion of lava and gas release interspersed with long repose times. This tendency was used to refine the calculation of the amount of SO2 emitted during variably vigorous activity. The continuous seismic tremor time series was used to quantify how long during these eruption phases strong/weak activity was exhibited to improve the calculated SO2 flux during these Phases.The total SO2 emissions derived from field measurements align closely with results obtained by combining melt inclusion and groundmass glass analyses with lava effusion rate measurements (910 ± 230 kt SO2). Specifically, utilizing the maximum S content found in evolved melt inclusions and the least remaining S content in accompanying quenched groundmasses provides an identical result between field measurements and the petrological calculations. This suggests that the maximum SO2 release calculated from petrological estimates should be preferentially used to initialize gas dispersion models for basaltic eruptions when other measurements are lacking.During the eruption, the CALPUFF dispersion model was used to forecast ground-level exposure to SO2. The SO2 emission rates measured by DOAS were used as input for the dispersion model, with updates made when a significant change was measured. A detailed analysis of one mid-distance station over the entire eruption shows that the model performed very well at predicting the presence of volcanic SO2 when it was measured. However, it frequently predicted the presence of SO2 that was not measured and the concentrations forecasted had no correlation with the concentrations measured. Various approaches to improve the model forecast were tested, including updating plume height and SO2 flux source terms based on measurements. These approaches did not unambiguously improve the model performance but suggest that improvements might be achieved in more-polluted conditions.

Degradation of the 2014–2015 Holuhraun vent-proximal edifice in Iceland

AbstractThe earliest stages of volcanic vent degradation are rarely measured, leaving a gap in the knowledge that informs landform degradation models of cinder cones and other monogenetic vent structures. We documented the initial degradation of a 500-m-long spatter rampart at the primary vent of the 2014–2015 Holuhraun eruption in northern Iceland with high-resolution topographic change maps derived from terrestrial laser scanning (TLS) and photogrammetric surveys using an unoccupied aircraft system (UAS). Topographic differencing shows a total negative volume change of 42,637 m3, and a total positive volume change (basal deposition) of 10,304 m3 (primarily as deposition at the base of steeply sloping surfaces). Two distinct styles of volume changes were observed on the interior and exterior of the spatter rampart. Material on the interior of the vent was removed from oversteepened slopes by discrete rockfalls, while diffusive processes were qualitatively evident on the exterior slopes. We propose a novel conceptual landform evolution model for spatter ramparts that combines rockfall processes on the interior walls, diffusive gravitational sliding on the exterior slopes, and incorporates cooling contraction and compaction over the entire edifice to describe the observed modes of topographic change during the onset of degradation. Potential hazards at fresh spatter ramparts are rockfalls at high slope areas of the vent interior walls where contacts between spatter clasts are prone to weakening by fumarolic activity, weathering, and settling. To capture such hazards, our data suggest a cadence for monitoring changes yearly for the first few years post-eruption, and at longer intervals thereafter.

How well are aerosol–cloud interactions represented in climate models? – Part 1: Understanding the sulfate aerosol production from the 2014–15 Holuhraun eruption

Abstract. For over 6 months, the 2014–2015 effusive eruption at Holuhraun, Iceland, injected considerable amounts of sulfur dioxide (SO2) into the lower troposphere with a daily rate of up to one-third of the global emission rate, causing extensive air pollution across Europe. The large injection of SO2, which oxidises to form sulfate aerosol (SO42-), provides a natural experiment offering an ideal opportunity to scrutinise state-of-the-art general circulation models' (GCMs) representation of aerosol–cloud interactions (ACIs). Here we present Part 1 of a two-part model inter-comparison using the Holuhraun eruption as a framework to analyse ACIs. We use SO2 retrievals from the Infrared Atmospheric Sounding Interferometer (IASI) instrument and ground-based measurements of SO2 and SO42- mass concentrations across Europe, in conjunction with a trajectory analysis using the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model, to assess the spatial and chemical evolution of the volcanic plume as simulated by five GCMs and a chemical transport model (CTM). IASI retrievals of plume altitude and SO2 column load reveal that the volcanic perturbation is largely contained within the lower troposphere. Compared to the satellite observations, the models capture the spatial evolution and vertical variability of the plume reasonably well, although the models often overestimate the plume altitude. HYSPLIT trajectories are used to attribute to Holuhraun emissions 111 instances of elevated sulfurous surface mass concentrations recorded at European Monitoring and Evaluation Programme (EMEP) stations during September and October 2014. Comparisons with the simulated concentrations show that the modelled ratio of SO2 to SO42- during these pollution episodes is often underestimated and overestimated for the young and mature plume, respectively. Models with finer vertical resolutions near the surface are found to better capture these elevated sulfurous ground-level concentrations. Using an exponential function to describe the decay of observed surface mass concentration ratios of SO2 to SO42- with plume age, the in-plume oxidation rate constant is estimated as 0.032 ± 0.002 h−1 (1.30 ± 0.08 d e-folding time), with a near-vent ratio of 25 ± 5 (µg m−3 of SO2 / µg m−3 of SO42-). The majority of the corresponding derived modelled oxidation rate constants are lower than the observed estimate. This suggests that the representation of the oxidation pathway/s in the simulated plumes is too slow. Overall, despite their coarse spatial resolutions, the six models show reasonable skill in capturing the spatial and chemical evolution of the Holuhraun plume. This capable representation of the underlying aerosol perturbation is essential to enable the investigation of the eruption's impact on ACIs in the second part of this study.

Morphological transitions between lobate resurfacing and distal breakout lava flows in flood basalts: insights from analog experiments

Continental flood basalts (CFBs) are dominated by two characteristic lava morphologies. The first type, referred to as ‘compound’ or ‘hummocky pāhoehoe,’ exhibits pillow-like lava flow lobes with cross-sections of ~ 0.5–2 m and thin chilled margins. The second type, referred to as ‘simple’ or ‘sheet lobes’ preserves more massive, inflated flow interiors that are laterally continuous on scales of 100s of meters to kilometers. Previous hypotheses suggest that two factors may contribute to stratigraphic changes in morphology from ‘compound’ to ‘simple’: 1) increased eruption duration or 2) increased extrusion rate. We test the hypothesis that a large increase in extrusion rate would result in flow morphology transitioning from multiple small lobes to inflated sheet lobes due to a shift in flow propagation from intraflow resurfacing-dominated to marginal breakout-dominated. Using polyethylene glycol (PEG) wax extruded into a circular water-filled tank 130 cm in diameter, we produced larger, more complex experiments than previous studies. Our efforts simulated more complex lava fields which change flow morphology with distance from the eruptive vent, characteristic of CFBs. Whereas previous PEG studies linked extrusion rate to near-source surface morphologies, our experiments evaluated how flow propagation mechanisms change with variable extrusion rate and distance from the source. Two flow propagation styles were identified: 1) resurfacing, in which molten material breaks through the surface of a flow and covers the older crust and 2) marginal breakouts, in which molten material extends beyond the crust at the active distal margin of the flow. Flows that propagated via marginal breakouts were found to have lower proportions of resurfaced area and vice versa. We show that significant resurfacing is needed to preserve internal chilled boundaries within a flow and a low-extrusion-rate surface morphology, whereas marginal breakout-dominated flows tend to inflate the pillow-like surface morphology preserving a massive interior at great distances from the vent. Higher and more steady extrusion rates tend to decrease the extent of resurfacing and increase the distance between the source and preserved low-extrusion-rate surface morphologies. We find that an extrusion rate increase equivalent to a jump in the extrusion rate scaling factor, Ψ value, from < 1 to > 5 would be necessary to ensure a switch from resurfacing-dominated lobate morphologies to marginal breakout-dominated propagation style. This amounts to a factor of 125 increase in effusion rate for fissure eruptions and a factor of 625 for point source eruptions, assuming no change in vent geometry. This would be equivalent to an effusion rate of 0.2 m3/s, as documented in 1987–1990 Kīlauea eruptions, increasing to 125 m3/s, which was commonly measured during the 2014 Holuhraun eruption in Iceland and the 2018 eruption at Leilani Estates in Hawai‘i. Thus, we propose that continental flood basalts do not require unusually large effusion rates, but instead were active for a longer and more consistent time period than smaller-volume eruptions.

High temperature spectral emissivity of glass and crystal-bearing basalts

Temperature is a key parameter controlling the rheology of lava flows. However, the unpredictable hazard of eruptions prevents direct measurements of hot volcanic bodies. Thus, the temperature of lava products is mostly retrieved by remote sensing techniques (ground- or satellite-based detectors) build on measuring infrared (IR) radiance. These well-established techniques are however subjected to important errors related to the poor knowledge of one of the most critical parameters, namely spectral emissivity. The aim of our study is to fill this gap by exploring the emissivity–temperature relationship of three different basalts. We performed in situ emissivity measurements at high temperatures (up to 1800 K) over a wide spectral range (350–8000 cm−1; 1.25–28.6 μm), using a non-contact IR apparatus. To unravel the complex radiative behavior of the samples with temperature, structural, chemical and textural analyses using Raman spectroscopy, XRD, DSC, SEM, EMPA and TEM were systematically performed. Our results show that spectral emissivity varies with temperature, wavenumber, and the sample degree of polymerization combined with total FeO content. Spectral emissivity is greatly affected by the crystallization of Fe-Mg-rich crystals at the micro-scale. Consequently, this study proves that spectral emissivity cannot be assumed constant for a single magmatic composition. Finally, our laboratory-measured values of spectral emissivity allowed refining the remote sensing temperature of the 2014–2015 Holuhraun eruption (Bárðarbunga, volcano) by 50 K in average. These new insights will ultimately reduce uncertainties in temperature estimates based on IR technology, information that are crucial to improve hazard assessment in volcanic crisis.

Linking lava morphologies to effusion rates for the 2014–2015 Holuhraun lava flow field, Iceland

AbstractDetermining the parameters that control fissure-fed lava morphologies is critical for reconstructing the complex emplacement histories of eruptions on Earth and other planetary bodies. We used a geomorphological map of the 2014–2015 Holuhraun lava flow field, in combination with new constraints on lava emplacement chronology and two independently derived time-averaged discharge rate (TADR) data sets, to analyze correlations between lava morphology and effusion rate. Results show that lava morphologies are dominantly controlled by effusion rate at the vent during the early phases of the eruption and by lava transport processes as the system evolves. Initially, TADR and its variance, which reflect pulsation in the lava supply rate from the vent, directly affect lava emplacement styles. However, as the eruption progresses, the lava transport system exerts a stronger control with channels and ponds that can either dampen variation in local effusion rate or create surges during sudden drainage events. The Holuhraun eruption predominantly produced rubbly lava in its earlier eruption phases and transitioned into the production of spiny lava toward the end of the eruption. However, a drop of TADR during the first phase of the eruption correlates with a decrease in rubbly lava formation and an increase in spiny lava production. This suggests that a change in effusion rate caused the observed transition in lava type. Our findings show that rubbly lava is formed under relatively high local effusion rates with pulsating supply conditions, whereas spiny lava is formed under lower local effusion rates and steadier supply.

Geomorphological characterization of the 2014–2015 Holuhraun lava flow-field in Iceland

The 2014–2015 Holuhraun eruption extruded >1 km3 of lava in a barren region of the Icelandic highlands. Due to its large volume and the abundance of data for this eruption, Holuhraun is an ideal site to investigate fissure-fed eruption products for comparison with other large lava flow-fields on Earth and other planetary bodies. To characterize lava morphologies associated with this event, we used 0.01–0.5 m/pixel image data, acquired from aerial surveys and small Unoccupied Aircraft Systems (sUAS) to create a 1:800-scale facies map that was ground-truthed using field observations during the summers of 2015–2019. Each facies region exhibits similar attributes in the remote sensing data, including albedo, surface texture, and geomorphology. However, at our mapping scale of 1:800, the facies typically include mixtures of lava types. Results show that transitional lava types (e.g., rubbly pāhoehoe, slabby pāhoehoe, and spiny pāhoehoe) dominate the 2014–2015 Holuhraun lava flow-field (83.82 km2), rather than the traditional end-members of ʻaʻā and classical pāhoehoe. At 1:800-scale, we distinguish the following eight facies (with the percentage of total flow-field area shown in parentheses): rubbly (57.35%), spiny (25.96%), undifferentiated rubbly–spiny (9.59%), shelly (5.58%), pāhoehoe (1.24%), flat-lying–knobby (0.58%), vent-proximal edifice (0.19%), and channel interior (0.16%). Field observations show that initial coherent pāhoehoe surfaces were episodically disrupted to produce slabby and rubbly textures that resemble ʻaʻā in remote sensing data. Our observations also show that continued solidification of the lava beneath brecciated surfaces can cause the surfaces of disrupted lobes to stabilize, or restabilize and undergo inflation. These factors complicate the use of surface texture as a direct indicator of emplacement style, which can change over the course of an eruption. This complexity has important implications for reconstructing the emplacement history of flow-fields on Earth and other planetary bodies.

Increased respiratory morbidity associated with exposure to a mature volcanic plume from a large Icelandic fissure eruption

The 2014–15 Holuhraun eruption in Iceland was the largest fissure eruption in over 200 years, emitting prodigious amounts of gas and particulate matter into the troposphere. Reykjavík, the capital area of Iceland (250 km from eruption site) was exposed to air pollution events from advection of (i) a relatively young and chemically primitive volcanic plume with a high sulphur dioxide gas (SO2) to sulphate PM (SO42−) ratio, and (ii) an older and chemically mature volcanic plume with a low SO2/SO42− ratio. Whereas the advection and air pollution caused by the primitive plume were successfully forecast and forewarned in public advisories, the mature plume was not. Here, we show that exposure to the mature plume is associated with an increase in register-measured health care utilisation for respiratory disease by 23% (95% CI 19.7–27.4%) and for asthma medication dispensing by 19.3% (95% CI 9.6–29.1%). Absence of public advisories is associated with increases in visits to primary care medical doctors and to the hospital emergency department. We recommend that operational response to volcanic air pollution considers both primitive and mature types of plumes.

Lava–water interaction and hydrothermal activity within the 2014–2015 Holuhraun Lava Flow Field, Iceland

Lava that erupted during the 2014–2015 Holuhraun eruption in Iceland flowed into a proglacial river system, resulting in aqueous cooling of the lava and an ephemeral hydrothermal system. We carried out a monitoring study of this system from 2015 to 2018 to document the cooling of the lava over this time, using thermocouple measurements and data-logging sensors. The heat loss rate from advection through this hydrothermal system in August 2015 was ~5.5 × 108 W; since eruption, aqueous cooling likely accounted for ~1% of the total heat loss from the lava. This estimate excludes steam losses from fumaroles as well as any groundwater that was not released to the surface, and thus is a lower bound. Near the terminus of the flow, advection of heat by flowing water may have locally accounted for tens of percent of the total cooling of that part of the flow. Our data quantify the importance of water cooling for this lava flow and can be compared with models to better understand lava–water interactions more generally. We also provide detailed methods for simple, low-cost monitoring of similar instances in the future.

Probabilistic hazard maps for operational use: the case of SO2 air pollution during the Holuhraun eruption (Bárðarbunga, Iceland) in 2014–2015

The Holuhraun fissure eruption (Iceland) in 2014–2015, which originated from the Barðarbunga volcanic system, was exceptional in several respects. It lasted 6 months and, throughout its duration, it released up to 9.6 Mt of SO2 in the atmosphere. The main recorded hazard affecting the entire country over the 6 months was the constant presence of a low-level gas cloud that led to recurrent air pollution episodes. The Icelandic Meteorological Office responded to this human health hazard by (1) setting up a forecasting system to anticipate the distribution of SO2 over Iceland and (2) preparing probabilistic hazard maps to support the decisions taken by the Icelandic Civil Protection in demarcating the accessible area around the eruption site. This paper introduces some technical aspects of the application of the CALPUFF numerical model to this eruption like the SO2 dispersal forecasting setup, the volcanic source numerical description, and the Monte Carlo procedure adopted for the creation of the probabilistic hazard maps. CALPUFF-based maps were created in January 2015, when the eruption was still ongoing, with the assumption that the eruption would be continuing with the same intensity. Maps for the entire country and for a smaller domain were produced, the latter showing the likelihood to exceeding an hourly concentration of 2600 μg/m3 (1 ppm) of SO2 for the spring season, a level chosen by the Icelandic Civil Protection for the delineation of the area of restricted access around the eruption site. As during the eruption there was no time for a rigorous evaluation of the model accuracy, we then undertook a retrospective analysis of CALPUFF model performance comparing the forecasted hourly SO2 concentration with real-time measurements at key-sites. The model did reproduce the hourly observations (the maximum within a 24-running window) with a level of agreement of 50.6% in Mývatn (85 km from the eruption site) and 50.4% in Reykjavik (258 km), when instances of null pairs have been removed. In Mývatn, the model overestimated the concentration more than 22% of the time. In Hofn (104 km), the model accuracy is 81.7% and occurrences of underestimation are higher than 11% of the time. In Reyðarfjorður (at 124 km), the model accuracy is assessed to be 82.7% and the model overestimates occurring 15.1% of the time. Possible explanations for the observed mismatch between model results and measurements include the spatial resolution of meteorological data field, the capability into reproducing chemical reactions of SO2 in atmosphere, and the reduced extension of the numerical domain. In addition, the model performance is strongly dependent on the source descriptors (e.g., strength of the SO2 flux and injection height) that, in this contribution, have been kept constant over long periods—neglecting, in this way, the natural variability of a dynamic emission of SO2. This consideration points toward the need of frequent and high-quality observation data for the initialization of dispersal numerical model. In light of this retrospective analysis, the probabilistic hazard maps possibly over-estimated the area exposed to high level of SO2 concentration. All the same, this paper reports on how quantitative probabilistic hazard mapping can be used for mitigating the health risk of volcanic SO2 emissions during a volcanic crisis to the benefit of operational hazard monitoring in support an effective crisis response.

The dynamics of a long-lasting effusive eruption modulated by Earth tides

Deciphering eruptive dynamics in near-real time is essential when dealing with hazard assessment and population evacuation. Recognition of specific patterns in time-series measured during volcanic activity may help decipher distinctive behavior at active volcanoes, providing insights into the underlying driving mechanisms. Multi-parameter data sets usually agree on the overall trend characterizing the temporal evolution of an eruption providing insights into the first-order eruptive dynamics. However, second-order variations detected in different data sets remain often poorly understood. The 2014-2015 Holuhraun eruption (Iceland) offers an excellent opportunity to investigate the dynamics of a long-lasting effusive eruption. We analyze the seismic tremor and the volcanic radiated power emitted by the lava field during the 6 months of the eruption using Singular Spectrum Analysis (SSA). In both geophysical time-series, we identify periods from ∼5 to ∼32 days coinciding with periods of Earth tides. Here we show that ∼50% of both signals are composed of tidal periods suggesting that magma movements follow frequencies imposed by lunisolar forces within the crust and at Earth's surface.

A Reappraisal of Seismicity Recorded During the 1996 Gjálp Eruption, Iceland, in Light of the 2014–2015 Bárðarbunga–Holuhraun Lateral Dike Intrusion

This study performs a reanalysis of the seismicity recorded during the 1996 Gjalp eruption that occurred at NW Vatnajokull, Iceland. The seismicity was recorded by the temporary HOTSPOT network consisting of 30 three-component broadband stations. In total 301 events were identified between 29 September and 12 October and their phases were manually picked. A velocity model was estimated from P-phase travel times by using VELEST. Events were first located using the algorithm NONLINLOC in order to obtain absolute locations. Precise relative locations were obtained with HYPODD by utilizing catalog and cross-correlation differential travel times. Results show that events clustered first along the SW rim of the Barðarbunga caldera and later along the Gjalp fissure, with most hypocentral depths located between 3 and 8 km. Waveforms of the 10 largest events that followed the Barðarbunga earthquake were inverted in order to obtain moment tensors. For all events we found that the deviatoric moment tensor fits the data better than pure double-couple or full moment tensor solutions. Events along the Barðarbunga caldera exhibited reverse focal mechanisms, while those at the Gjalp fissure exhibited mostly strike-slip faulting. Seismic velocity variations calculated using ambient noise interferometry, point to the possibility that a small subglacial eruption occurred at Barðarbunga before the main earthquake. This removed melt from the magma chamber causing its roof to collapse, and also resulted in the lateral migration of magma towards the Gjalp fissure. The 2014–2015 Barðarbunga–Holuhraun eruption shares common characteristics with the 1996 Gjalp eruption, although the size of the latter was much smaller.

Operational response and hazards assessment during the 2014–2015 volcanic crisis at Bárðarbunga volcano and associated eruption at Holuhraun, Iceland

On August 16, 2014, an intense seismic swarm started below the eastern part of Bárðarbunga caldera in the north-western part of Vatnajökull ice cap, Iceland, marking the onset of a major rifting event. The intense seismicity, corroborated by a complicated deformation pattern, migrated away from the Bárðarbunga caldera for >45 km, until a small effusive eruption started in Holuhraun, outside of the ice cap, on August 29 which lasted a few hours. On August 31 a second, more intense, eruption started. It lasted 180 days and is the largest effusive eruption in Iceland since the Laki event in 1783–1784. A lava field covered about 84 km2 in the middle of the highlands with a bulk volume of ~1.44 km3. Although the potential danger of ice melting and subsequent flooding coupled with tephra fallout was the main concern at the beginning of the crisis, the main associated hazard turned out to be an abundant and constant release of volcanic gases that polluted the country threatening human health and the environment for several months. Very little (or almost no) ash was produced. The impact of the Holuhraun eruption outside of Iceland was minor if compared with what happened in 2010 when Eyjafjallajökull erupted. However, on three occasions the dense gas cloud released by the fissure eruption reached some European countries where unusually high SO2 concentrations were measured at ground level. Managing and living with such a prolonged effusive event was a major effort for the Icelandic scientific community and civil protection authorities and it significantly affected the local population in the northern and eastern parts of the country.This paper retraces the main events that accompanied the volcanic crisis at Bárðarbunga and focuses on the scientific and operative response provided throughout its duration. Particular attention is given to the direct links between the observed phenomena, monitoring effort, hazard assessment and concrete actions taken by the authorities in charge of risk and mitigation assessment. The different roles of the involved institutions are clarified, and the wide and collaborative response on a national basis is explained.A major finding of this post-event analysis is that a strong internal collaboration between the scientific community and operative institutions is the basic element for proper crisis management. Iceland has established protocols that, while imperfect, ensured a well-coordinated response to the Holuhraun eruption and the preceding unrest, and will continue to provide a prompt and effective response for future events. Important realizations, including the need for a proper ground-based alert level system, have been raised and they will be a matter for further development.

Large-scale particulate air pollution and chemical fingerprint of volcanic sulfate aerosols from the 2014–2015 Holuhraun flood lava eruption of Bárðarbunga volcano (Iceland)

Abstract. Volcanic sulfate aerosols play a key role in air quality and climate. However, the rate of oxidation of sulfur dioxide (SO2) precursor gas to sulfate aerosols (SO42-) in volcanic plumes is poorly known, especially in the troposphere. Here we determine the chemical speciation as well as the intensity and temporal persistence of the impact on air quality of sulfate aerosols from the 2014–2015 Holuhraun flood lava eruption of Icelandic volcano Bárðarbunga. To do so, we jointly analyse a set of SO2 observations from satellite (OMPS and IASI) and ground-level measurements from air quality monitoring stations together with high temporal resolution mass spectrometry measurements of an Aerosol Chemical Speciation Monitor (ACSM) performed far from the volcanic source. We explore month/year long ACSM data in France from stations in contrasting environments, close and far from industrial sulfur-rich activities. We demonstrate that volcanic sulfate aerosols exhibit a distinct chemical signature in urban/rural conditions, with NO3:SO4 mass concentration ratios lower than for non-volcanic background aerosols. These results are supported by thermodynamic simulations of aerosol composition, using the ISORROPIA II model, which show that ammonium sulfate aerosols are preferentially formed at a high concentration of sulfate, leading to a decrease in the production of particulate ammonium nitrate. Such a chemical signature is however more difficult to identify at heavily polluted industrial sites due to a high level of background noise in sulfur. Nevertheless, aged volcanic sulfates can be distinguished from freshly emitted industrial sulfates according to their contrasting degree of anion neutralization. Combining AERONET (AErosol RObotic NETwork) sunphotometric data with ACSM observations, we also show a long persistence over weeks of pollution in volcanic sulfate aerosols, while SO2 pollution disappears in a few days at most. Finally, gathering 6-month long datasets from 27 sulfur monitoring stations of the EMEP (European Monitoring and Evaluation Programme) network allows us to demonstrate a much broader large-scale European pollution, in both SO2 and SO4, associated with the Holuhraun eruption, from Scandinavia to France. While widespread SO2 anomalies, with ground-level mass concentrations far exceeding background values, almost entirely result from the volcanic source, the origin of sulfate aerosols is more complex. Using a multi-site concentration-weighted trajectory analysis, emissions from the Holuhraun eruption are shown to be one of the main sources of SO4 at all EMEP sites across Europe and can be distinguished from anthropogenic emissions from eastern Europe but also from Great Britain. A wide variability in SO2:SO4 mass concentration ratios, ranging from 0.8 to 8.0, is shown at several stations geographically dispersed at thousands of kilometres from the eruption site. Despite this apparent spatial complexity, we demonstrate that these mass oxidation ratios can be explained by a simple linear dependency on the age of the plume, with a SO2-to-SO4 oxidation rate of 0.23 h−1. Most current studies generally focus on SO2, an unambiguous and more readily measured marker of the volcanic plume. However, the long persistence of the chemical fingerprint of volcanic sulfate aerosols at continental scale, as shown for the Holuhraun eruption here, casts light on the impact of tropospheric eruptions and passive degassing activities on air quality, health, atmospheric chemistry and climate.

Clinopyroxene–Liquid Equilibria and Geothermobarometry in Natural and Experimental Tholeiites: the 2014–2015 Holuhraun Eruption, Iceland

Abstract Clinopyroxene–liquid geothermobarometry is a widely used tool for estimating the conditions under which mafic magmas are stored before they erupt. However, redox variability, sector zoning and disequilibrium crystallization present major challenges to the robust estimation of magma storage conditions. Moreover, most recent studies seeking to address these challenges have focused on clinopyroxenes from alkalic systems and are thus of limited use for understanding clinopyroxenes from the tholeiitic systems that dominate global magma budgets. Here we combine observations on natural clinopyroxenes from the 2014–2015 Holuhraun lava in Iceland with observations on experimental clinopyroxenes synthesized during high-pressure, high-temperature experiments on the same lava in order to investigate clinopyroxene–liquid equilibria in tholeiitic systems and optimize of geothermobarometric strategies. Natural clinopyroxenes from the 2014–2015 Holuhraun lava are sector zoned, with {1-11} hourglass sectors being enriched in the enstatite–ferrosillite component at the expense of all other components with respect to {hk0} prism sectors. In contrast with observations on clinopyroxenes from alkalic systems, sector zoning in clinopyroxenes from the 2014–2015 Holuhraun lava is characterized by differences in Ca and Na contents as well as in Ti and Al contents. The products of crystallization experiments performed at 100–600 MPa and 1140–1220 °C on a powdered starting glass at two sets of melt H2O content–oxygen fugacity conditions (∼0·1 wt % H2O and close to the graphite-oxygen redox buffer, and 0·5–1·0 wt % H2O and approximately one and half log units above the quartz–fayalite–magnetite redox buffer) demonstrate that clinopyroxene crystals from nominally equilibrium experiments can preserve strongly disequilibrium compositions. The compositional systematics of experimental clinopyroxenes are consistent with the presence of sector zoning. Furthermore, the magnitude of compositional variability increases with decreasing melt H2O content and increasing deviations of experimental temperatures below clinopyroxene liquidus temperatures (i.e. degrees of undercooling sensu lato), indicating that kinetic processes play a key role in controlling clinopyroxene compositions, even under notionally equilibrium conditions. Few published analyses of experimental clinopyroxene crystals may thus represent truly equilibrium compositions. Stoichiometric calculations on natural and experimental clinopyroxenes show that Fe3+ is a major constituent of clinopyroxenes from tholeiitic magmas under naturally relevant oxygen fugacity conditions. They also show that Fe3+ is most likely incorporated as Ca- and Al- bearing Ca–Fe-Tschermak’s component rather than Na-bearing aegirine component at oxygen fugacities up to one and a half log units above the quartz–fayalite–magnetite buffer. Elevated oxygen fugacities are thus less likely to compromise clinopyroxene–liquid geothermobarometry than previously thought. Guided by our experimental results, we combined published descriptions of clinopyroxene–liquid equilibria with geothermobarometric equations to develop an internally consistent and widely applicable method for performing geothermobarometry on tholeiitic magmas that does not require equilibrium zones to be selected a priori. Applying this method to natural clinopyroxene crystals from the 2014–2015 Holuhraun lava that formed under low but variable degrees of undercooling (perhaps 25 °C or less) returns values in excellent agreement with those from independent methods (232 ± 86 MPa, 1161 ± 11 °C). Robust estimates of magma storage conditions can thus be obtained by performing clinopyroxene–liquid geothermobarometry on tholeiitic magmas when disequilibrium is suitably accounted for.

Landscape evolution associated with the 2014–2015 Holuhraun eruption in Iceland

The 2014–2015 Holuhraun eruption in Iceland developed between the outlet glacier Dyngjujökull and the Askja central volcano and extruded a bulk lava volume of over 1 km3 onto the floodplain of the Jökulsá á Fjöllum river, making it the largest effusive eruption in Iceland during the past 230 years. Time-series monitoring using a combination of traditional aerial imaging, unmanned aerial systems, and field-based geodetic surveys, established an unprecedented record of the hydrological response of the river system to this lava flow. We observed: (1) the formation of lava-dammed lakes and channels produced during dam-breaching events; (2) percolation of glacial meltwater into the porous and permeable lava, forming an ephemeral hydrothermal system that included hot pools and hot springs that emerged from the lava flow front; and (3) the formation of new seepage channels caused by upwelling of water around the periphery of the lava flow. The observations show that lava flows, like the one produced by the 2014–2015 Holuhraun eruption, can cause significant hydrological changes that continue for several years after the lava is emplaced. Documenting these processes is therefore crucial for our interpretation of volcanic landscapes and processes of lava–water interaction on both Earth and Mars.

Satellite-derived sulfur dioxide (SO&amp;lt;sub&amp;gt;2&amp;lt;/sub&amp;gt;) emissions from the 2014–2015 Holuhraun eruption (Iceland)

Abstract. The 6-month-long 2014–2015 Holuhraun eruption was the largest in Iceland for 200 years, emitting huge quantities of sulfur dioxide (SO2) into the troposphere, at times overwhelming European anthropogenic emissions. Weather, terrain and latitude made continuous ground-based or UV satellite sensor measurements challenging. Infrared Atmospheric Sounding Interferometer (IASI) data are used to derive the first time series of daily SO2 mass present in the atmosphere and its vertical distribution over the entire eruption period. A new optimal estimation scheme is used to calculate daily SO2 fluxes and average e-folding time every 12 h. For the 6 months studied, the SO2 flux was observed to be up to 200 kt day−1 and the minimum total SO2 erupted mass was 4.4±0.8 Tg. The average SO2 e-folding time was 2.4±0.6 days. Where comparisons are possible, these results broadly agree with ground-based near-source measurements, independent remote-sensing data and values obtained from model simulations from a previous paper. The results highlight the importance of using high-resolution time series data to accurately estimate volcanic SO2 emissions. The SO2 mass missed due to thermal contrast is estimated to be of the order of 3 % of the total emission when compared to measurements by the Ozone Monitoring Instrument. A statistical correction for cloud based on the AVHRR cloud-CCI data set suggested that the SO2 mass missed due to cloud cover could be significant, up to a factor of 2 for the plume within the first kilometre from the vent. Applying this correction results in a total erupted mass of 6.7±0.4 Tg and little change in average e-folding time. The data set derived can be used for comparisons to other ground- and satellite-based measurements and to petrological estimates of the SO2 flux. It could also be used to initialise climate model simulations, helping to better quantify the environmental and climatic impacts of future Icelandic fissure eruptions and simulations of past large-scale flood lava eruptions.

Respiratory health among professionals exposed to extreme SO2 levels from a volcanic eruption.

Objective The Holuhraun eruption of fall and winter 2014-15 produced large amounts of sulfur dioxide (SO 2). The aim of this study was to determine if exposure to extreme SO 2levels affected the health of individuals working at the eruption site. Methods During January‒March 2015, earth scientists, technicians, and law enforcement personnel who were about to work at the eruption site were invited to a respiratory health examination. Symptom reports and lung function measures, forced expiratory volume in one second (FEV 1) and forced vital capacity (FVC) were collected before and after an eruption site visit. Those with previous exposure (N=27) reported symptoms retrospectively. Results Altogether, 41 individuals were invited to participate, 32 underwent a clinical examination at a hospital respiratory health clinic (baseline); 27 reported symptoms during earlier visits to the eruption site (retrospective symptom reports), 17 were re-examined 1-6 days after visiting the eruption site (follow-up). All participants' lung function was within normal range both before and after exposure. At baseline, average FEV 1was 107.4% of predicted versus 106.6 at follow-up (P =0.82); average FVC was 107.0% of predicted at baseline versus 107.4% at follow-up (P=0.35). Eye and nasal irritation were more frequently reported during eruption site exposure by 24% versus 6% (P =0.37) for both. Conclusion Although "healthy-worker" effects cannot be excluded, our data indicate that SO 2exposure was associated with relatively mild and transient respiratory symptoms with no clinical signs of airway inflammation or airway obstruction.