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

Large earthquakes alter the stress in the surrounding crust, leading to triggered earthquakes and aftershocks. A number of time-dependent processes, including afterslip, pore-fluid flow and viscous relaxation of the lower crust and upper mantle, further modify the stress and pore pressure near the fault, and hence the tendency for triggered earthquakes. It has proved difficult, however, to distinguish between these processes on the basis of direct field observations, despite considerable effort. Here we present a unique combination of measurements consisting of satellite radar interferograms and water-level changes in geothermal wells following two magnitude-6.5 earthquakes in the south Iceland seismic zone. The deformation recorded in the interferograms cannot be explained by either afterslip or visco-elastic relaxation, but is consistent with rebound of a porous elastic material in the first 1-2 months following the earthquakes. This interpretation is confirmed by direct measurements which show rapid (1-2-month) recovery of the earthquake-induced water-level changes. In contrast, the duration of the aftershock sequence is projected to be approximately 3.5 years, suggesting that pore-fluid flow does not control aftershock duration. But because the surface strains are dominated by pore-pressure changes in the shallow crust, we cannot rule out a longer pore-pressure transient at the depth of the aftershocks. The aftershock duration is consistent with models of seismicity rate variations based on rate- and state-dependent friction laws.

Gradual inflation of magma chambers often precedes eruptions at highly active volcanoes. During such eruptions, rapid deflation occurs as magma flows out and pressure is reduced. Less is known about the deformation style at moderately active volcanoes, such as Eyjafjallajökull, Iceland, where an explosive summit eruption of trachyandesite beginning on 14 April 2010 caused exceptional disruption to air traffic, closing airspace over much of Europe for days. This eruption was preceded by an effusive flank eruption of basalt from 20 March to 12 April 2010. The 2010 eruptions are the culmination of 18 years of intermittent volcanic unrest. Here we show that deformation associated with the eruptions was unusual because it did not relate to pressure changes within a single magma chamber. Deformation was rapid before the first eruption (>5 mm per day after 4 March), but negligible during it. Lack of distinct co-eruptive deflation indicates that the net volume of magma drained from shallow depth during this eruption was small; rather, magma flowed from considerable depth. Before the eruption, a ∼0.05 km(3) magmatic intrusion grew over a period of three months, in a temporally and spatially complex manner, as revealed by GPS (Global Positioning System) geodetic measurements and interferometric analysis of satellite radar images. The second eruption occurred within the ice-capped caldera of the volcano, with explosivity amplified by magma-ice interaction. Gradual contraction of a source, distinct from the pre-eruptive inflation sources, is evident from geodetic data. Eyjafjallajökull's behaviour can be attributed to its off-rift setting with a 'cold' subsurface structure and limited magma at shallow depth, as may be typical for moderately active volcanoes. Clear signs of volcanic unrest signals over years to weeks may indicate reawakening of such volcanoes, whereas immediate short-term eruption precursors may be subtle and difficult to detect.

Most of the Earth's crust is created along 60,000 km of mid-ocean ridge system. Here, tectonic plates spread apart and, in doing so, gradually build up stress. This stress is released during rifting episodes, when bursts of magmatic activity lead to the injection of vertical sheets of magma -termed dykes -into the crust. Only 2% of the global mid-ocean ridge system is above sea level, so making direct observations of the rifting process is difficult. However, geodetic and seismic observations exist from spreading centres in Afar (East Africa) and Iceland that are exposed at the land surface. Rifting episodes are rare, but the few that have been well observed at these sites have operated with remarkably similar mechanisms. Specifically, magma is supplied to the crust in an intermittent manner, and is stored at multiple positions and depths. It then laterally intrudes in dykes within the brittle upper crust. Depending on the availability of magma, multiple magma centres can interact during one rifting episode. If we are to forecast large eruptions at spreading centres, rifting-cycle models will need to fully incorporate realistic crust and mantle properties, as well as the dynamic transport of magma.

No abstract

No abstract

The 18th historic eruption of Hekla started on 26 February, 2000. It was a short-lived but intense event, emitting basaltic andesitic (55.5 wt% SiO 2 ) pyroclastic fragments and lava. During the course of the eruption, monitoring was done by both instruments and direct observations, together providing unique insight into the current activity of Hekla. During the 12-day eruption, a total of 0.189 km 3 DRE of magma was emitted. The eruptive fissure split into five segments. The segments at the highest altitude were active during the first hours, while the segments at lower altitude continued throughout the eruption. The eruption started in a highly explosive manner giving rise to a Subplinian eruptive column and consequent basaltic pyroclastic flows fed by column collapses. After the explosive phase reached its maximum, the eruption went through three more phases, namely fire-fountaining, Strombolian bursts and lava effusion. In this paper, we describe the course of events of the eruption of Hekla and the origin of its magma, and then show that the discharge rate can be linked to different style of eruptive activity, which are controlled by fissure geometry. We also show that the eruption phases observed at Hekla can be linked with inferred magma chamber overpressure prior to the eruption.

Dynamically triggered seismicity followed shortly after a Ms 6.6 earthquake in Iceland on June 17, 2000. Smaller earthquakes occurred on the Reykjanes Peninsula up to 100 km from the mainshock rupture. Using interferometric analysis of Synthetic Aperture Radar images (InSAR), we measure crustal deformation associated with three triggered deformation events. The largest of these occurred at Lake Kleifarvatn, 85 km west of the mainshock epicenter. Modeling of the InSAR data reveals strikeslip on a north‐striking fault, with a geodetic moment of 6.2 × 1017 Nm, equivalent to magnitude Mw 5.8 earthquake. A seismological estimate of the moment is not yet available, because the seismic signature of this event is partly hidden by the mainshock waveform. The paucity of aftershocks on the triggered rupture plane suggests some aseismic slip there, compatible with a thin seismogenic crust, high heat‐flow, hydrothermal alteration and the presence of fluids in the area.