40Large volcanic eruptions on Earth commonly occur with collapse of the roof of a crustal magma 41 reservoir, forming a caldera. Only a few such collapses occur per century and lack of detailed 42 observations has obscured insight on mechanical interplay between collapse and eruption. We use Calderas are 1 -100 km diameter depressions found in volcanic regions of Earth and other planets. basaltic andesite) intrusive activity and eruptions (2,(9)(10)(11)(12). 59The consensus from field and modelling studies is that caldera collapse progresses from initial 60 surface downsag to fault-controlled subsidence (1, 8, 13, 14). The pre-collapse topography is obtained by subtracting the subsidence observed at the surface. As we recorded the caldera subsidence mainly on the ice (Fig. 1, Fig. S1), we made corrections and (Fig. 3A). We therefore conclude that suggestions of a large increase in ice flow out of the caldera 147 during these events (25) cannot be fitted with our data. 148Bedrock subsidence exceeding 1 m occurred within an area of 110 km 2 that extended beyond the 149 pre-existing caldera (Fig. 1, Fig. S1). After termination of collapse the total subsidence at the pre-150 existing caldera rims amounted to 3 to 11 meters ( Fig. 1D and 1E). Using subglacial radio-echo GPS station in the center of the caldera (Fig. 1A), including the rate of vertical rate of ice surface Cumulative number of M>4 caldera earthquakes, with magnitude evolution colored in red, blue and 176 grey representing clusters on the southern rim, the northern rim and smaller clusters, respectively 177 (see Fig. S5). E) Cumulative seismic moment for M>4 caldera earthquakes. from analysis of subaerial gas measurements (Fig. 4). This depth concurs with our regional on FTIR and Multi-GAS measurements (24). 194Seismicity and subsurface structure 195 We used seismic data and Distinct Element Method (DEM) numerical modelling (24), to 196 characterize the deeper collapse structure as the reactivation of a steeply-inclined ring fault (Fig. 5). 197We mostly observed seismicity at depths of 0-9 km beneath the northern and southern caldera rims 198( Fig. 5B), with earthquakes being more numerous on the northern rim. This spatial pattern of 199 seismicity is consistent with fracturing above a deflating magma reservoir that was elliptical in (Fig. 5C, D). Our best fitting models had preexisting faults dipping out at 80-85¡ from the caldera 207 center on the north side and at 85-90¡ toward the caldera center on the south side. The modeled pre- 208existing faults lay at 1-2 km below the surface on the north side and 3-4 km on the south side. 209Modeling of a more complex fault geometry or the inclusion of greater material heterogeneity may 210 further improve the data fit, but presently lacks robust geophysical constraints. components of the observed earthquakes at B ‡rdarbunga. We, however, narrowed down on 222 plausible solutions by using the micro-earthquakes (Fig. 5A). The moment tensor solutions are well 223 constrained, but the inferred d...

The effusive six months long 2014-2015 Bárðarbunga eruption (31 August-27 February) was the largest in Iceland for more than 200 years, producing 1.6 ± 0.3 km 3 of lava. The total SO 2 emission was 11 ± 5 Mt, more than the amount emitted from Europe in 2011. The ground level concentration of SO 2 exceeded the 350 µg m −3 hourly average health limit over much of Iceland for days to weeks. Anomalously high SO 2 concentrations were also measured at several locations in Europe in September. The lowest pH of fresh snowmelt at the eruption site was 3.3, and 3.2 in precipitation 105 km away from the source. Elevated dissolved H 2 SO 4 , HCl, HF, and metal concentrations were measured in snow and precipitation. Environmental pressures from the eruption and impacts on populated areas were reduced by its remoteness, timing, and the weather. The anticipated primary environmental pressure is on the surface waters, soils, and vegetation of Iceland.

The 6-month long eruption at Holuhraun (August 2014-February 2015) in the Bárðarbunga-Veiðivötn volcanic system was the largest effusive eruption in Iceland since the 1783-1784 CE Laki eruption. The lava flow field covered~84 km 2 and has an estimated bulk (i.e., including vesicles) volume of~1.44 km 3. The eruption had an average discharge rate of~90 m 3 /s making it the longest effusive eruption in modern times to sustain such high average flux. The first phase of the eruption (August 31, 2014 to mid-October 2014) had a discharge rate of~350 to 100 m 3 /s and was typified by lava transport via open channels and the formation of four lava flows, no. 1-4, which were emplaced side by side. The eruption began on a 1.8 km long fissure, feeding partly incandescent sheets of slabby pāhoehoe up to 500 m wide. By the following day the lava transport got confined to open channels and the dominant lava morphology changed to rubbly pāhoehoe and 'a'ā. The latter became the dominating morphology of lava flows no. 1-8. The second phase of the eruption (Mid-October to end November) had a discharge of~100-50 m 3 /s. During this time the lava transport system changed, via the formation of a b 1 km 2 lava pond~1 km east of the vent. The pond most likely formed in a topographical low created by a the pre-existing Holuhraun and the new Holuhraun lava flow fields. This pond became the main point of lava distribution, controlling the emplacement of subsequent flows (i.e. no. 5-8). Towards the end of this phase inflation plateaus developed in lava flow no. 1. These inflation plateaus were the surface manifestation of a growing lava tube system, which formed as lava ponded in the open lava channels creating sufficient lavastatic pressure in the fluid lava to lift the roof of the lava channels. This allowed new lava into the previously active lava channel lifting the channel roof via inflation. The final (third) phase, lasting from December to end-February 2015 had a mean discharge rate of~50 m 3 /s. In this phase the lava transport was mainly confined to lava tubes within lava flows no. 1-2, which fed breakouts that resurfaced N 19 km 2 of the flow field. The primary lava morphology from this phase was spiny pāhoehoe, which superimposed on the 'a'ā lava flows no. 1-3 and extended the entire length of the flow field (i.e. 17 km). This made the 2014-2015 Holuhraun a paired flow field, where both lava morphologies had similar length. We suggest that the similar length is a consequence of the pāhoehoe is fed from the tube system utilizing the existing 'a'ā lava channels, and thereby are controlled by the initial length of the 'a'ā flows.

Explosive volcanic eruptions are characterized by the rapid fragmentation of a magmatic melt into ash particles. In order to describe the energy dissipation during fragmentation it is important to understand the mechanism of material failure. A quantitative description of fragmentation is only possible under controlled laboratory conditions. Industrial silicate glasses have a high structural affinity with magmatic melts and have the advantage of being transparent, which allows the study of the evolution of fractures by optical methods on a time scale relevant for explosive volcanism. With this aim, a series of low speed edge-on hammer impact experiments on silicate glass targets has been conducted, leading to the generation of fragments in the grain-size spectra of volcanic ash. In order to verify the general transferability of the experimentally generated fragmentation dynamics to volcanic processes, the resulting products were compared, by means of statistical particle-shape analyses, to particles produced by standardized magma fragmentation experiments and to natural ash particles coming from deposits of basaltic and rhyolitic compositions from the 2004 Grimsvötn and the Quaternary Tepexitl tuff-ring eruptions, respectively. Natural ash particles from both Grimsvötn and Tepexitl show significant similarities with experimental fragments of thermally pre-stressed float glasses, indicating a dominant influence of preexisting stresses on particle shape and suggesting analogous fragmentation processes within the studied materials.

Fissure eruptions are commonly linked to magma dikes at depth and are associated with elastic and inelastic surface deformation. Elastic deformation is well described by subsidence occurring above the dike plane and uplift and lateral widening occurring perpendicular to the dike plane. Inelastic deformation is associated with the formation of a graben, which is bordered by graben parallel faults that might express as sets of fractures at the surface. Additionally, secondary structures, such as push-ups, bends and step overs, yield information about the deforming domain. However, once these structures are formed during fissure eruptions, they are rarely preserved in nature, due to the effects of rapid erosion, sediment coverage or overprinting by other faulting events. Therefore, simple normal fault displacements are commonly assumed at dikes. At the 2014/2015 Holuhraun eruption sites (Iceland), increasing evidence suggests that developing fractures exhibited variations in their displacement modes. In an attempt to investigate these variations, a fieldwork mapping project combining Terrestrial Laser Scanning (TLS) and Unmanned Aerial Vehicle (UAV)-based aerophoto analysis was undertaken. Using these data, we generated local high-resolution Digital Elevation Models (DEMs) and a structural map that facilitated the identification of kinematic indicators and the assessment of the observed structures. We identified 315 fracture segments from these satellite data. We measured the strike directions of single segments, including the amount of opening and opening angles, which indicate that many of the measured fractures show transtensional dislocations. Of these, ∼81% exhibit a significant left-lateral component and only ∼17% exhibit a right-lateral component. Here, we demonstrate that the local complexities in these fracture traces and geometries are closely related to variations in their transtensional opening directions. Moreover, we identified local changes in fracture azimuths and offsets close to eruption sites, which we speculate are associated with geometric changes in the magma feeder itself. The results highlight that the opening of fractures associated with an erupting fissure Müller et al. Structures at Holuhraun may record transtensional modes with both, left-lateral and right-lateral components. These results further highlight the value of using UAV-based high-resolution data to contribute to the integrity of the observations of the structural complexities produced by local geologic events.

The source conditions of volcanic plumes and collapsing fountains are investigated by means of large-scale experiments. In the experiments, gas-particle jets issuing from a cylindrical conduit are forced into the atmosphere at different mass flow rates. Dense jets (high particle volumetric concentration, e.g., C 0 > 0.01) generate collapsing fountains, whose height scales with the squared exit velocity. This is consistent with Bernoulli's equation, which is a good approximation if air entrainment is negligible. In this case, kinetic energy is transformed into potential energy without any significant loss by friction with the atmosphere. The dense collapsing fountain, on hitting the ground, generates an intense shear flow similar to a pyroclastic density current. Dilute hot jets (low particle volumetric concentration, e.g., C 0 < 0.01) dissipate their initial kinetic energy at much smaller heights than those predicted by Bernoulli's equation. This is an indication that part of the total mechanical energy is lost by friction with the atmosphere. Significant air entrainment results in this case, leading to the formation of a buoyant column (plume) from which particles settle similarly to pyroclastic fallout. The direct measurement of entrainment coefficient in the experiments suggests that dense collapsing fountains form only when air entrainment is not significant. This is a consequence of the large density difference between the jet and the atmosphere. Cold dilute experiments result in an entrainment coefficient of about 0.06, which is typical of pure jets of fluid dynamics. Hot dilute experiments result in an entrainment coefficient of about 0.11, which is typical of thermally buoyant plumes. The entrainment coefficients obtained by experiments were used as input data in numerical simulations of fountains and plumes. A numerical model was used to solve the classic top-hat system of governing equations, which averages the field variables (e.g., column velocity and density) across the column. The maximum heights calculated with the model agree well with those observed experimentally, showing that our entrainment coefficients are compatible with a top-hat model. Dimensional analysis of the experimental data shows that a value of 3 for the source densimetric Froude number characterizes the transition between dense collapsing fountains and dilute plumes. This value delimits the source conditions (exit velocity, conduit radius, and particle volumetric concentration) for pyroclastic flow (<3) and fallout (>3). © 2014 Springer-Verlag Berlin Heidelberg

The majority of Earth's volcanic eruptions occur beneath the sea, but few direct observations and samples limit our understanding of these unseen events. Subaerial eruptions lend some insights, but direct extrapolation from subaerial to deep-sea is precluded by the great differences in pressure, thermal conditions, density, rheology, and the interplay among them. Here we present laboratory fragmentation experiments that mimic deep-sea explosive eruptions and compare our laboratory observations with those from the kilometre-deep submarine eruption of Havre volcano, Kermadec arc, New Zealand in 2012. We find that the Havre eruption involved explosive fragmentation of magma by a pressure-insensitive interaction between cool water and

The 2010 eruption of Eyjafjallajökull volcano was characterized by pulsating activity. Discrete ash bursts merged at higher altitude and formed a sustained quasi-continuous eruption column. High-resolution near-field videos were recorded on 8-10 May, during the second explosive phase of the eruption, and supplemented by contemporary aerial observations. In the observed period, pulses occurred at intervals of 0.8 to 23.4 s (average, 4.2 s). On the basis of video analysis, the pulse volume and the velocity of the reversely buoyant jets that initiated each pulse were determined. The expansion history of jets was tracked until the pulses reached the height of transition from a negatively buoyant jet to a convective buoyant plume about 100 m above the vent. Based on the assumption that the density of the gas-solid mixture making up the pulse approximates that of the surrounding air at the level of transition from the jet to the plume, a mass flux ranging between 2.2 and 3.5 · 10 4 kg/s was calculated. This mass eruption rate is in good agreement with results obtained with simple models relating plume height with mass discharge at the vent. Our findings indicate that near-field measurements of eruption source parameters in a pulsating eruption may prove to be an effective monitoring tool. A comparison of the observed pulses with those generated in calibrated large-scale experiments reveals very similar characteristics and suggests that the analysis of near-field sensors could in the future help to constrain the triggering mechanism of explosive eruptions.