No abstract

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

Northern Hemisphere summer temperatures over the past 8000 years have been paced by the slow decrease in summer insolation resulting from the precession of the equinoxes. However, the causes of superposed century‐scale cold summer anomalies, of which the Little Ice Age (LIA) is the most extreme, remain debated, largely because the natural forcings are either weak or, in the case of volcanism, short lived. Here we present precisely dated records of ice‐cap growth from Arctic Canada and Iceland showing that LIA summer cold and ice growth began abruptly between 1275 and 1300 AD, followed by a substantial intensification 1430–1455 AD. Intervals of sudden ice growth coincide with two of the most volcanically perturbed half centuries of the past millennium. A transient climate model simulation shows that explosive volcanism produces abrupt summer cooling at these times, and that cold summers can be maintained by sea‐ice/ocean feedbacks long after volcanic aerosols are removed. Our results suggest that the onset of the LIA can be linked to an unusual 50‐year‐long episode with four large sulfur‐rich explosive eruptions, each with global sulfate loading >60 Tg. The persistence of cold summers is best explained by consequent sea‐ice/ocean feedbacks during a hemispheric summer insolation minimum; large changes in solar irradiance are not required.

Active volcanoes and hydrothermal systems underlie ice caps in Iceland. Glacier-volcano interactions produce meltwater that either drains toward the glacier margin or accumulates in subglacial lakes. Accumulated meltwater drains periodically in jökulhlaups from the subglacial lakes and occasionally during volcanic eruptions. The release of meltwater from glacial lakes can take place in two different mechanisms. Drainage can begin at pressures lower than the ice overburden in conduits that expand slowly due to melting of the ice walls by frictional and sensible heat in the water. Alternatively, the lake level rises until the ice dam is lifted and water pressure in excess of the ice overburden opens the waterways; the glacier is lifted along the flowpath to make space for the water. In this case, discharge rises faster than can be accommodated by melting of the conduits. Normally jökulhlaups do not lead to glacier surges but eruptions in ice-capped stratovolcanoes have caused rapid and extensive glacier sliding. Jökulhlaups from subglacial lakes may transport on the order of 10 7 tons of sediment per event but during violent volcanic eruptions, the sediment load has been 10 8 tons.

Surges are common in all the major ice caps in Iceland, and historical reports of surge occurrence go back several centuries. Data collection and regular observation over the last several decades have permitted a detailed description of several surges, from which it is possible to generalize on the nature of surging in Icelandic glaciers. Combining the historical records of glacier-front variations and recent field research, we summarize the geographic distribution of surging glaciers, their subglacial topography and geology, the frequency and duration of surges, changes in glacier surface geometry during the surge cycle, and measured velocity changes compared to calculated balance velocities. We note the indicators of surge onset and describe changes in ice, water and sediment fluxes during a surge. Surges accomplish a significant fraction of the total mass transport through the main outlet glaciers of ice caps in Iceland and have important implications for their hydrology. Our analysis of the data suggests that surge-type glaciers in Iceland are characterized by gently sloping surfaces and that they move too slowly to remain in balance given their accumulation rate. Surge frequency is neither regular nor clearly related to glacier size or mass balance. Steeply sloping glaciers, whether hard- or soft-bedded, seem to move sufficiently rapidly to keep in balance with the annual accumulation.