Search

Increased rates of deformation and seismicity are well-established precursors to volcanic eruptions, and their interpretation forms the basis for eruption warnings worldwide. Rates of ground displacement and the number of earthquakes escalate before many eruptions 1-3 , as magma forces its way towards the surface. However, the pre-eruptive patterns of deformation and seismicity vary widely. Here we show how an eruption beginning on 19 March 2021 at Fagradalsfjall, Iceland, was preceded by a period of tectonic stress release ending with a decline in deformation and seismicity over several days preceding the eruption onset. High rates of deformation and seismicity occurred from 24 February to mid-March in relation to gradual emplacement of an approximately 9-km-long magma-filled dyke, between the surface and 8 km depth (volume approximately 34 × 10 6  m 3 ), as well as the triggering of strike-slip earthquakes up to magnitude M W 5.64. As stored tectonic stress was systematically released, there was less lateral migration of magma and a reduction in both the deformation rates and seismicity. Weaker crust near the surface may also have contributed to reduced seismicity, as the depth of active magma emplacement progressively shallowed. This demonstrates that the interaction between volcanoes and tectonic stress as well as crustal layering need to be fully considered when forecasting eruptions., (© 2022. The Author(s).)

Large volume effusive eruptions with relatively minor observed precursory signals are at odds with widely used models to interpret volcano deformation. Here we propose a new modelling framework that resolves this discrepancy by accounting for magma buoyancy, viscoelastic crustal properties, and sustained magma channels. At low magma accumulation rates, the stability of deep magma bodies is governed by the magma-host rock density contrast and the magma body thickness. During eruptions, inelastic processes including magma mush erosion and thermal effects, can form a sustained channel that supports magma flow, driven by the pressure difference between the magma body and surface vents. At failure onset, it may be difficult to forecast the final eruption volume; pressure in a magma body may drop well below the lithostatic load, create under-pressure and initiate a caldera collapse, despite only modest precursors.

Crust at many divergent plate boundaries forms primarily by the injection of vertical sheet-like dykes, some tens of kilometres long. Previous models of rifting events indicate either lateral dyke growth away from a feeding source, with propagation rates decreasing as the dyke lengthens, or magma flowing vertically into dykes from an underlying source, with the role of topography on the evolution of lateral dykes not clear. Here we show how a recent segmented dyke intrusion in the Bárðarbunga volcanic system grew laterally for more than 45 kilometres at a variable rate, with topography influencing the direction of propagation. Barriers at the ends of each segment were overcome by the build-up of pressure in the dyke end; then a new segment formed and dyke lengthening temporarily peaked. The dyke evolution, which occurred primarily over 14 days, was revealed by propagating seismicity, ground deformation mapped by Global Positioning System (GPS), interferometric analysis of satellite radar images (InSAR), and graben formation. The strike of the dyke segments varies from an initially radial direction away from the Bárðarbunga caldera, towards alignment with that expected from regional stress at the distal end. A model minimizing the combined strain and gravitational potential energy explains the propagation path. Dyke opening and seismicity focused at the most distal segment at any given time, and were simultaneous with magma source deflation and slow collapse at the Bárðarbunga caldera, accompanied by a series of magnitude M > 5 earthquakes. Dyke growth was slowed down by an effusive fissure eruption near the end of the dyke. Lateral dyke growth with segment barrier breaking by pressure build-up in the dyke distal end explains how focused upwelling of magma under central volcanoes is effectively redistributed over long distances to create new upper crust at divergent plate boundaries.

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.

Pressure influences both magma production and the failure of magma chambers. Changes in pressure interact with the local tectonic settings and can affect magmatic activity. Present-day reduction in ice load on subglacial volcanoes due to global warming is modifying pressure conditions in magmatic systems. The large pulse in volcanic production at the end of the last glaciation in Iceland suggests a link between unloading and volcanism, and models of that process can help to evaluate future scenarios. A viscoelastic model of glacio-isostatic adjustment that considers melt generation demonstrates how surface unloading may lead to a pulse in magmatic activity. Iceland's ice caps have been thinning since 1890 and glacial rebound at rates exceeding 20 mm yr(-1) is ongoing. Modelling predicts a significant amount of 'additional' magma generation under Iceland due to ice retreat. The unloading also influences stress conditions in shallow magma chambers, modifying their failure conditions in a manner that depends critically on ice retreat, the shape and depth of magma chambers as well as the compressibility of the magma. An annual cycle of land elevation in Iceland, due to seasonal variation of ice mass, indicates an annual modulation of failure conditions in subglacial magma chambers.