Ármann Höskuldsson

Ash generation and distribution from the April-May 2010 eruption of Eyjafjallajökull, Iceland

The 39-day long eruption at the summit of Eyjafjallajökull volcano in April–May 2010 was of modest size but ash was widely dispersed. By combining data from ground surveys and remote sensing we show that the erupted material was 4.8±1.2·1011 kg (benmoreite and trachyte, dense rock equivalent volume 0.18±0.05 km3). About 20% was lava and water-transported tephra, 80% was airborne tephra (bulk volume 0.27 km3) transported by 3–10 km high plumes. The airborne tephra was mostly fine ash (diameter <1000 µm). At least 7·1010 kg (70 Tg) was very fine ash (<28 µm), several times more than previously estimated via satellite retrievals. About 50% of the tephra fell in Iceland with the remainder carried towards south and east, detected over ~7 million km2 in Europe and the North Atlantic. Of order 1010 kg (2%) are considered to have been transported longer than 600–700 km with <108 kg (<0.02%) reaching mainland Europe.

Tephra sedimentation during the 2010 Eyjafjallajökull eruption (Iceland) from deposit, radar, and satellite observations

The April-May 2010 eruption of Eyjafjallajokull volcano (Iceland) was characterized by a nearly continuous injection of tephra in the atmosphere that affected various economic sectors in Iceland and caused a global-wide interruption of air traffic. Eruptive activity during 4-8 May 2010 was characterized based on short-duration physical parameters in order to capture transient eruptive behavior of a long-lasting eruption (i.e., total grainsize distribution, erupted mass and mass eruption rate averaged over 30-minute activity). Resulting 30-minute total grainsize distribution based on both ground and MSG-SEVIRI satellite measurements is characterized by Mdphi of about 2 phi and a fine-ash content of about 30wt%. Accumulation rate varied by two orders of magnitude with an exponential decay away from the vent, whereas Mdphi shows a linear increase until about 18 km from vent reaching a plateau of about 4.5 phi between 20-56 km. Associated mass eruption rate is in between 0.6-1.2 x 10^5 kg s^-1. In-situ sampling showed how fine ash mainly fell as aggregates of various typologies. About 5 to 9 wt% of the erupted mass remained in the cloud up to 1000 km from the vent, suggesting that nearly half of the ash >7 phi settled as aggregates within the first 60 km. Particle sphericity and shape factor varied between 0.4 and 1 with no clear correlation with size and distance from vent. Our experiments also demonstrate how satellite retrievals and Doppler radar grainsize detection can provide real-time description of the source term but for a limited particle-size range.

MeMoVolc report on classification and dynamics of volcanic explosive eruptions

Classifications of volcanic eruptions were first introduced in the early twentieth century mostly based on qualitative observations of eruptive activity, and over time, they have gradually been developed to incorporate more quantitative descriptions of the eruptive products from both deposits and observations of active volcanoes. Progress in physical volcanology, and increased capability in monitoring, measuring and modelling of explosive eruptions, has highlighted shortcomings in the way we classify eruptions and triggered a debate around the need for eruption classification and the advantages and disadvantages of existing classification schemes. Here, we (i) review and assess existing classification schemes, focussing on subaerial eruptions; (ii) summarize the fundamental processes that drive and parameters that characterize explosive volcanism; (iii) identify and prioritize the main research that will improve the understanding, characterization and classification of volcanic eruptions and (iv) provide a roadmap for producing a rational and comprehensive classification scheme. In particular, classification schemes need to be objective-driven and simple enough to permit scientific exchange and promote transfer of knowledge beyond the scientific community. Schemes should be comprehensive and encompass a variety of products, eruptive styles and processes, including for example, lava flows, pyroclastic density currents, gas emissions and cinder cone or caldera formation. Open questions, processes and parameters that need to be addressed and better characterized in order to develop more comprehensive classification schemes and to advance our understanding of volcanic eruptions include conduit processes and dynamics, abrupt transitions in eruption regime, unsteadiness, eruption energy and energy balance.

Satellite detection, long-range transport, and air quality impacts of volcanic sulfur dioxide from the 2014-2015 flood lava eruption at Bárðarbunga (Iceland)

The 2014–2015 Barðarbunga-Veiðivotn fissure eruption at Holuhraun produced about 1.5 km3 of lava, making it the largest eruption in Iceland in more than 200 years. Over the course of the eruption, daily volcanic sulfur dioxide (SO2) emissions exceeded daily SO2 emissions from all anthropogenic sources in Europe in 2010 by at least a factor of 3. We present surface air quality observations from across Northern Europe together with satellite remote sensing data and model simulations of volcanic SO2 for September 2014. We show that volcanic SO2 was transported in the lowermost troposphere over long distances and detected by air quality monitoring stations up to 2750 km away from the source. Using retrievals from the Ozone Monitoring Instrument (OMI) and the Infrared Atmospheric Sounding Interferometer (IASI), we calculate an average daily SO2 mass burden of 99 ± 49 kilotons (kt) of SO2 from OMI and 61 ± 18 kt of SO2 from IASI for September 2014. This volcanic burden is at least a factor of 2 greater than the average SO2 mass burden between 2007 and 2009 due to anthropogenic emissions from the whole of Europe. Combining the observational data with model simulations using the United Kingdom Met Offices Numerical Atmospheric-dispersion Modelling Environment model, we are able to constrain SO2 emission rates to up to 120 kilotons per day (kt/d) during early September 2014, followed by a decrease to 20–60 kt/d between 6 and 22 September 2014, followed by a renewed increase to 60–120 kt/d until the end of September 2014. Based on these fluxes, we estimate that the eruption emitted a total of 2.0 ± 0.6 Tg of SO2 during September 2014, in good agreement with ground-based remote sensing and petrological estimates. Although satellite-derived and model-simulated vertical column densities of SO2 agree well, the model simulations are biased low by up to a factor of 8 when compared to surface observations of volcanic SO2 on 6–7 September 2014 in Ireland. These biases are mainly due to relatively small horizontal and vertical positional errors in the simulations of the volcanic plume occurring over transport distances of thousands of kilometers. Although the volcanic air pollution episodes were transient and lava-dominated volcanic eruptions are sporadic events, the observations suggest that (i) during an eruption, volcanic SO2 measurements should be assimilated for near real-time air quality forecasting and (ii) existing air quality monitoring networks should be retained or extended to monitor SO2 and other volcanic pollutants.

Late Pleistocene to Holocene eruptive activity of Pico de Orizaba, Eastern Mexico

The Late Pleistocene to Holocene eruptive history of Pico de Orizaba can be divided into 11 eurptive episodes. Each eruptive episode lasted several hundred years, the longest recorded being about 1000 years (the Xilomich episode). Intervals of dormancy range from millenia during the late Pleistocene to about 500 years, the shortest interval recorded in the Holocene. This difference could reflect either changes in the volcanos activity or that the older stratigraphic record is less complete than the younger. Eruptive mechanisms during the late Pleistocene were characterized by dome extrusions, lava flows and ash-and-scoria-flow generating eruptive columns. However, in Holocene time plinian activity became increasingly important. The increase in dacitic plinian eruptions over time is related to increased volumes of dacitic magma beneath Pico de Orizaba. We suggest that the magma reservoir under Pico de Orizaba is stratified. The last eruptive episode, which lasted from about 690 years bp until ad 1687, was initiated by a dacitic plinian eruption and was followed by effusive lava-forming eruptions. For the last 5,000 years the activity of the volcano has been gradually evolving towards such a trend, underlining the increasing importance of dacitic magma and stratification of the magma reservoir. Independent observations of Pico de Orizabas glacier early this century indicate that some increase in volcanic activity occurred between 1906 and 1947, and that it was probably fumarolic.

Geology of the Heimaey volcanic centre, south Iceland: early evolution of a central volcano in a propagating rift?

Geology of the Heimaey volcanic centre, south Iceland: early evolution of a central volcano in a propagating rift?

Magma–ice–sediment interactions and the origin of lava/hyaloclastite sequences in the Síða formation, South Iceland

Products of subglacial volcanism can illuminate reconstructions of paleo-environmental conditions on both local and regional scales. Competing interpretations of Pleistocene conditions in south Iceland have been proposed based on an extensive sequence of repeating lava-and-hyaloclastite deposits in the Síða district. We propose here a new eruptive model and refine the glacial environment during eruption based on field research and analytical data for the Síða district lava/hyaloclastite units. Field observations from this and previous studies reveal a repeating sequence of cogenetic lava and hyaloclastite deposits extending many kilometers from their presumed eruptive source. Glasses from lava selvages and unaltered hyaloclastites have very low H2O, S, and CO2 concentrations, indicating significant degassing at or close to atmospheric pressure prior to quenching. We also present a scenario that demonstrates virtual co-emplacement of the two eruptive products. Our data and model results suggest repeated eruptions under thin ice or partially subaerial conditions, rather than eruption under a thick ice sheet or subglacial conditions as previously proposed.

Contrasting mechanisms of magma fragmentation during coeval magmatic and hydromagmatic activity: the Hverfjall Fires fissure eruption, Iceland

Growing evidence for significant magmatic vesiculation prior to magma-water interaction (MWI) has brought into question the use of ‘diagnostic’ features, such as low vesicularities and blocky morphologies, to identify hydromagmatic pyroclasts. We address this question by quantifying co-variations in particle size, shape and texture in both magmatic and hydromagmatic deposits from the Hverfjall Fires fissure eruption, Iceland. Overlapping vesicularity and bubble number density distributions measured in rapidly quenched magmatic and hydromagmatic pyroclasts indicate a shared initial history of bubble nucleation and growth, with substantial vesiculation prior to MWI. Hydromagmatic fragmentation occurred principally by brittle mechanisms, where the length scale and geometry of fracturing was controlled by the bubble population. This suggests that the elevated fragmentation efficiency of hydromagmatic deposits is driven, at least in part, by brittle disintegration of vesicular pyroclasts due to high thermal stress generated during rapid cooling. In this way, the shape and size distributions of hydromagmatic pyroclasts, both critical input parameters for ash dispersion models, are strongly influenced by the dynamics of vesiculation prior to MWI. This result underlines the need to analyse multiple grain-size fractions to characterise the balance between magmatic and hydromagmatic processes. During the Hverfjall Fires eruption, the external water supply was sufficient to maintain MWI throughout the eruption, with no evidence for progressive exhaustion of a water reservoir. We suggest that both the longevity and the spatial distribution of MWI were determined by the pre-existing regional hydrology and represent continuous interaction between a propagating dike and a strong groundwater flow system hosted within permeable basalt lavas.

Multiple Coincident Eruptive Seismic Tremor Sources During the 2014-2015 Eruption at Holuhraun, Iceland

We analyze eruptive tremor during one of the largest effusive eruptions in historical times in Iceland (2014/15 Holuhraun eruption). Seismic array recordings are compared with effusion rates deduced from MODerate resolution Imaging Spectroradiometer (MODIS) recordings and ground video monitoring data and lead to the identification of three coexisting eruptive tremor sources. This contrasts other tremor studies that generally link eruptive tremor to only one source usually associated with the vent. The three sources are (i) a source that is stable in back azimuth and shows bursts with ramp-like decrease in amplitude at the beginning of the eruption. We link it to a process below the open vents where the bursts correlate with the opening of new vents and temporary increases in the lava fountaining height. (ii) a source moving by a few degrees per month while the tremor amplitude suddenly increases and decreases. Back azimuth and slowness correlate with the growing margins of the lava flow field whilst contact with a river led to fast increases of the tremor amplitude. (iii) a source moving by up to 25∘ southwards in 4 days that cannot be related to any observed surface activity and might be linked to intrusions. We therefore suggest that eruptive tremor amplitudes/ energies are used with caution when estimating eruptive volumes, effusion rates or the eruption explosivity as multiple sources can coexist during the eruption phase. Our results suggest that arrays can monitor both the growth of a lava flow field and the activity in the vents.

The 2014 Lake Askja rockslide‐induced tsunami: Optimization of numerical tsunami model using observed data

A large rockslide was released from the inner Askja caldera into Lake Askja, Iceland, on 21 July 2014. Upon entering the lake, it caused a large tsunami that traveled about ∼3 km across the lake and inundated the shore with vertical runup measuring up to 60–80 m. Following the event, comprehensive field data were collected, including GPS measurements of the inundation and multibeam echo soundings of the lake bathymetry. Using this exhaustive data set, numerical modeling of the tsunami has been conducted using both a nonlinear shallow water model and a Boussinesq-type model that includes frequency dispersion. To constrain unknown landslide parameters, a global optimization algorithm, Differential Evolution, was employed, resulting in a parameter set that minimized the deviation from measured inundation. The tsunami model of Lake Askja is the first example where we have been able to utilize field data to show that frequency dispersion is needed to explain the tsunami wave radiation pattern and that shallow water theory falls short. We were able to fit the trend in tsunami runup observations around the entire lake using the Boussinesq model. In contrast, the shallow water model gave a different runup pattern and produced pronounced offsets in certain areas. The well-documented Lake Askja tsunami thus provided a unique opportunity to explore and capture the essential physics of landslide tsunami generation and propagation through numerical modeling. Moreover, the study of the event is important because this dispersive nature is likely to occur for other subaerial impact tsunamis.