On 31 August a new eruption began from the same fissure and is still ongoing at the time of writing. After 4 September the movement associated with the dyke was minor, suggesting an approximate equilibrium between inflow of magma into the dyke and magma flowing out of it feeding the eruption. Minor eruptions may have occurred under Vatnajškull; shallow ice depressions marked by circular crevasses (ice cauldrons) were discovered in the period 27/08-07/09, indicating leakage of magma or magmatic heat to the glacier causing basal melting ( Fig. 1 and 2b). On 5 September, aircraft radar profiling showed that the ice surface in the centre of the B ‡r!arbunga caldera had subsided 16 m relative to the surroundings, resulting in a 0.32±0.08 km 3 subsidence bowl ( can be compared to a 1 day interferogram over the ice surface spanning 27 -28 August (Fig. 1), that has maximum line-of-sight (LOS) increase of 57 cm, indicating 55-70 cm of subsidence, during 24 hours. From 24 August to 6 September 16 M≥5 earthquakes occurred on the caldera boundary.Over 22000 earthquakes were automatically detected 16/08-06/09 2014, 5000 of which have been manually checked. Four thousand of these have been relatively relocated, defining the dyke segments. Ground deformation in areas outside the Vatnajškull ice cap, and on nunataks within the ice cap, is well mapped by a combination of InSAR, continuously recording GPS sites, and campaign GPS measurements. The GPS observations and analysis give the temporal evolution of the three-dimensional displacements used in the modelling (Fig. 1). Interferometric analysis of synthetic aperture radar images from the COSMO-SkyMed, RADARSAT-2 and TerraSAR-X satellites was used to form 11 interferograms showing LOS change spanning different time intervals (Supplementary Fig. 2). The analysis of seismic and geodetic data is described in Methods.Initial modelling of the dyke, with no a priori constraints on position, strike or dip, show the deformation data require the dyke to be approximately vertical and line up with the seismicity (Extended Data item 4). We therefore fixed the dip to be vertical and the lateral position of the dyke to coincide with the earthquake locations.We modelled the dyke as a series of rectangular patches and estimated the opening and slip on each patch ( Fig. 3a; see Supplementary Figures 3-4 for slip and standard deviations of opening). We used a Markov-chain Monte Carlo approach to estimate 7 the multivariate probability distribution for all model parameters (Methods) on each day 16/08-06/09 2014 (Fig. 2d). The results suggest that most of the magma injected into the dyke is shallower than the seismicity, which mostly spans the depth range from 5 to 8 km below sea level (see Fig. 2c and Methods). While magma may extend to depths greater than 9 km near the centre of the ice cap, towards the edge of the ice cap where constraints from InSAR and GPS are much better, significant opening is all shallower than 5 km (Fig. 3a). The total volume intruded into the dyke by 28 August was 0.48-0...

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...

[1] Using a combination of body wave and surface wave data sets to reveal the mantle plume and plume head, this study presents a tomographic image of the mantle structure beneath Iceland to 400 km depth. Data comes primarily from the PASSCAL-HOTSPOT deployment of 30 broadband instruments over a period of 2 years, and is supplemented by data from the SIL and ICEMELT networks. Three sets of relative teleseismic body wave arrival times are generated through cross correlation: S and SKS arrivals at 0.03-0.1 Hz, and P and PKIKP arrivals at 0.03-0.1 and 0.8-2.0 Hz. Prior to inversion the crustal portion of the travel time anomalies is removed using the crustal model ICECRTb. This step has a significant effect on the mantle velocity variations imaged down to a depth of $250 km. Inversion of relative arrival times only provides information on lateral velocity variations. Surface waves are therefore used to provide absolute velocity information for the uppermost mantle beneath Iceland. The average wave number for the Love wave fundamental mode at 0.020 and 0.024 Hz is measured and used to invert for the average S velocity. Combination of the body wave and surface wave information reveals a predominantly horizontal low-velocity anomaly extending from the Moho down to $250 km depth, interpreted as a plume head. Below the plume head a near-cylindrical lowvelocity anomaly with a radius of $100 km and peak V P and V S anomalies of À2% and À4%, respectively, extends down to the maximum depth of resolution at 400 km. Within the plume head, in the uppermost mantle above the core of the plume, there is a relatively high velocity with a maximum V P and V S anomaly of +2%. This high-velocity anomaly may be the result of the extreme degree of melt extraction necessary to generate the thick (46 km) crust in central Iceland. Comparison of the plume volumetric flux implied by our images, the crustal generation rate, and the degree of melting suggested by rare earth element inversions, suggests that (1) mantle material must be flowing horizontally away from the plume core faster than the overlying lithosphere and (2)

[1] Through combination of surface wave and body wave constraints we derive a threedimensional (3-D) crustal S velocity model and Moho map for Iceland. It reveals a vast plumbing system feeding mantle plume melt into upper crustal magma chambers where crustal formation takes place. The method is based on the partitioned waveform inversion to which we add additional observations. Love waves from six local events recorded on the HOTSPOT-SIL networks are fitted, S n travel times from the same events measured, previous observations of crustal thickness are added, and all three sets of constraints simultaneously inverted for our 3-D model. In the upper crust (0-15 km) an elongated low-velocity region extends along the length of the Northern, Eastern and Western Neovolcanic Zones. The lowest velocities (À7%) are found at 5-10 km below the two most active volcanic complexes: Hekla and Bárdarbunga-Grímsvötn. In the lower crust (>15 km) the low-velocity region can be represented as a vertical cylinder beneath central Iceland. The low-velocity structure is interpreted as the thermal halo of pipe work which connects the region of melt generation in the uppermost mantle beneath central Iceland to active volcanoes along the neovolcanic zones. Crustal thickness in Iceland varies from 15-20 km beneath the Reykjanes Peninsula, Krafla and the extinct Snaefellsnes rift zone, to 46 km beneath central Iceland. The average crustal thickness is 29 km. The variations in thickness can be explained in terms of the temporal variation in plume productivity over the last $20 Myr, the Snaefellsnes rift zone being active during a minimum in plume productivity. Variations in crustal thickness do not depart significantly from an isostatically predicted crustal thickness. The best fit linear isostatic relation implies an average density jump of 4% across the Moho. Rare earth element inversions of basalt compositions on Iceland suggest a melt thickness (i.e., crustal thickness) of 15-20 km, given passive upwelling. The observed crustal thickness of up to 46 km implies active fluxing of source material through the melt zone by the mantle plume at up to 3 times the passive rate.

Summary We report the results of the highest‐resolution teleseismic tomography study yet performed of the upper mantle beneath Iceland. The experiment used data gathered by the Iceland Hotspot Project, which operated a 35‐station network of continuously recording, digital, broad‐band seismometers over all of Iceland 1996–1998. The structure of the upper mantle was determined using the ACH damped least‐squares method and involved 42 stations, 3159 P‐wave, and 1338 S‐wave arrival times, including the phases P, pP, sP, PP, SP, PcP, PKIKP, pPKIKP, S, sS, SS, SKS and Sdiff. Artefacts, both perceptual and parametric, were minimized by well‐tested smoothing techniques involving layer thinning and offset‐and‐averaging. Resolution is good beneath most of Iceland from ∼60 km depth to a maximum of ∼450 km depth and beneath the Tjornes Fracture Zone and near‐shore parts of the Reykjanes ridge. The results reveal a coherent, negative wave‐speed anomaly with a diameter of 200–250 km and anomalies in P‐wave speed, VP, as strong as −2.7 per cent and in S‐wave speed, VS, as strong as −4.9 per cent. The anomaly extends from the surface to the limit of good resolution at ∼450 km depth. In the upper ∼250 km it is centred beneath the eastern part of the Middle Volcanic Zone, coincident with the centre of the ∼100 mGal Bouguer gravity low over Iceland, and a lower crustal low‐velocity zone identified by receiver functions. This is probably the true centre of the Iceland hotspot. In the upper ∼200 km, the low‐wave‐speed body extends along the Reykjanes ridge but is sharply truncated beneath the Tjornes Fracture Zone. This suggests that material may flow unimpeded along the Reykjanes ridge from beneath Iceland but is blocked beneath the Tjornes Fracture Zone. The magnitudes of the VP, VS and VP/VS anomalies cannot be explained by elevated temperature alone, but favour a model of maximum temperature anomalies < 200 K, along with up to ∼2 per cent of partial melt in the depth range ∼100–300 km beneath east‐central Iceland. The anomalous body is approximately cylindrical in the top 250 km but tabular in shape at greater depth, elongated north–south and generally underlying the spreading plate boundary. Such a morphological change and its relationship to surface rift zones are predicted to occur in convective upwellings driven by basal heating, passive upwelling in response to plate separation and lateral temperature gradients. Although we cannot resolve structure deeper than ∼450 km, and do not detect a bottom to the anomaly, these models suggest that it extends no deeper than the mantle transition zone. Such models thus suggest a shallow origin for the Iceland hotspot rather than a deep mantle plume, and imply that the hotspot has been located on the spreading ridge in the centre of the north Atlantic for its entire history, and is not fixed relative to other Atlantic hotspots. The results are consistent with recent, regional full‐thickness mantle tomography and whole‐mantle tomography images that show a strong, low‐wave‐speed anomaly ben...