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 2014-2015 Holuhraun eruption, on the Bárðarbunga volcanic system in central Iceland, was one of the best-monitored basaltic fissure eruptions that has ever occurred, and presents a unique opportunity to link petrological and geochemical data with geophysical observations during a major rifting episode. We present major and trace element analyses of melt inclusions and matrix glasses from a suite of ten samples collected over the course of the Holuhraun eruption. The diversity of trace element ratios such as La/Yb in Holuhraun melt inclusions reveals that the magma evolved via concurrent mixing and crystallization of diverse primary melts in the mid-crust. Using olivine-plagioclase-augite-melt (OPAM) barometry, we calculate that the Holuhraun carrier melt equilibrated at 2.1 ± 0.7 kbar (7.5 ± 2.5 km), which is in agreement with the depths of earthquakes (6 ± 1 km) between Bárðarbunga central volcano and the eruption site in the days preceding eruption onset. Using the same approach, melt inclusions equilibrated at pressures between 0.5 and 8.0 kbar, with the most probable pressure being 3.2 kbar. Diffusion chronometry reveals minimum residence timescales of 1-12 days for melt inclusionbearing macrocrysts in the Holuhraun carrier melt. By combining timescales of diffusive dehydration of melt inclusions with the calculated pressure of H 2 O saturation for the Holuhraun magma, we calculate indicative magma ascent rates of 0.12-0.29 m s −1 . Our petrological and geochemical data are consistent with lateral magma transport from Bárðarbunga volcano to the eruption site in a shallow-to mid-crustal dyke, as has been suggested on the basis of seismic and geodetic datasets. This result is a significant step forward in reconciling petrological and geophysical interpretations of magma transport during volcano-tectonic episodes, and provides a critical framework for the interpretation of premonitory seismic and geodetic data in volcanically active regions.

The hydrogen-isotope [deuterium/hydrogen (D/H)] ratio of Earth can be used to constrain the origin of its water. However, the most accessible reservoir, Earth's oceans, may no longer represent the original (primordial) D/H ratio, owing to changes caused by water cycling between the surface and the interior. Thus, a reservoir completely isolated from surface processes is required to define Earth's original D/H signature. Here we present data for Baffin Island and Icelandic lavas, which suggest that the deep mantle has a low D/H ratio (δD more negative than -218 per mil). Such strongly negative values indicate the existence of a component within Earth's interior that inherited its D/H ratio directly from the protosolar nebula.

No abstract