4 research outputs found
Deformation at Krafla and Bjarnarflag geothermal areas, Northern Volcanic Zone of Iceland, 1993-2015
Get PDFThe Krafla volcanic system has geothermal areas within the Krafla caldera and at Bjarnarflag in the Krafla fissure swarm, 9-km south of the Krafla caldera. Arrays of boreholes extract geothermal fluids for power plants in both areas. We collected and analyzed InSAR, GPS, and leveling data spanning 1993–2015 in order to investigate crustal deformation in these areas. The volcanic zone hosting the geothermal areas is also subject to large scale regional deformation processes, including plate spreading and deflation of the Krafla volcanic system. These deformation processes have to be taken into account in order to isolate the geothermal deformation signal. Plate spreading produces the largest horizontal displacements, but the regional deformation pattern also suggests readjustment of the Krafla system at depth after the 1975–1984 Krafla rifting episode. Observed deformation can be fit by an inflation source at about 20 km depth north of Krafla and a deflation source at similar depth directly below the Krafla caldera. Deflation signal along the fissure swarm can be reproduced by a 1-km wide sill at 4 km depth closing by 2–4 cm per year. These sources are considered to approximate the combined effects of vertical deformation associated with plate spreading and post-rifting response. Local deformation at the geothermal areas is well resolved in addition to these signals. InSAR shows that deformation at Bjarnarflag is elongated along the direction of the Krafla fissure swarm (∼4 km by ∼2 km) while it is circular at Krafla (∼5 km diameter). Rates of deflation at Krafla and Bjarnarflag geothermal areas have been relatively steady. Average volume decrease of about 6.6 ×10⁵ m³/yr for Krafla and 3.9 ×10⁵ m³/yr for Bjanarflag are found at sources located at ∼1.5 km depth, when interpreted by a spherical point source of pressure. This volume change represents about 8 ×10 −3 m³/ton of the mass of geothermal fluid extracted per year, indicating important renewal of the geothermal reservoir by water flow
Post-eruptive volcano inflation following major magma drainage: Interplay between models of viscoelastic response influence and models of magma inflow at Bárðarbunga caldera, Iceland, 2015-2018
Get PDF
<p>Unrest at Bárðarbunga after a caldera collapse in 2014-2015 includes elevated seismicity beginning about six months after the eruption ended, including nine Mw>4.5 earthquakes. The earthquakes occurred mostly on the northern and southern parts of a caldera ring fault. Global Navigation Satellite System (GNSS, in particular, Global Positioning System; GPS) and Interferometric Synthetic Aperture Radar (InSAR) geodesy are applied to evaluate the spatial and temporal pattern of ground deformation around Bárðarbunga caldera outside the icecap, in 2015-2018, when deformation rates were relatively steady. The aim is to study the role of viscoelastic relaxation following major magma drainage versus renewed magma inflow as an explanation for the ongoing unrest.</p><p>The largest horizontal velocity is measured at GPS station KISA (3 km from caldera rim), 141 mm/yr in direction N47<sup>o</sup>E relative to the Eurasian plate in 2015-2018. GPS and InSAR observations show that the velocities decay rapidly outward from the caldera. We correct our observations for Glacial Isostatic Adjustment and plate spreading to extract the deformation related to volcanic activity. After this correction, some GPS sites show subsidence.</p><p>We use a reference Earth model to initially evaluate the contribution of viscoelastic processes to the observed deformation field. We model the deformation within a half-space composed of a 7-km thick elastic layer on top of a viscoelastic layer with a viscosity of 5 x 10<sup>18</sup> Pa s, considering two co-eruptive contributors to the viscoelastic relaxation: “non-piston” magma withdrawal at 10 km depth (modelled as pressure drop in a spherical source) and caldera collapse (modelled as surface unloading). The other model we test is the magma inflow in an elastic half-space. Both the viscoelastic relaxation and magma inflow create horizontal outward movements around the caldera, and uplift at the surface projection of the source center in 2015-2018. Viscoelastic response due to magma withdrawal results in subsidence in the area outside the icecap. Magma inflow creates rapid surface velocity decay as observed.</p><p>We explore further two parameters in the viscoelastic reference model: the viscosity and the "non-piston" magma withdrawal volume. Our comparison between the corrected InSAR velocities and viscoelastic models suggests a viscosity of 2.6×10<sup>18</sup> Pa s and 0.36 km<sup>3</sup> of “non-piston” magma withdrawal volume, given by the optimal reduced Chi-squared statistic. When the deformation is explained using only magma inflow into a single spherical source (and no viscoelastic response), the optimal model suggests an inflow rate at 1×10<sup>7</sup> m<sup>3</sup>/yr at 700 m depth. A magma inflow model with more model parameters is also a possible explanation, including sill inflation at 10 km together with slip on caldera ring faults. Our reference Earth model and the two end-member models suggest that there is a trade-off between the viscoelastic relaxation and the magma inflow, since they produce similar deformation signals outside the icecap. However, to reproduce details of the observed deformation, both processes are required. A viscoelastic-only model cannot fully explain the fast velocity decay away from the caldera, whereas a magma inflow-only model cannot explain the subsidence observed at several locations.</p>
</jats:p
Global link between deformation and volcanic eruption quantified by satellite imagery
Get PDFA key challenge for volcanological science and hazard management is that few of the world’s volcanoes are effectively monitored. Satellite imagery covers volcanoes globally throughout their eruptive cycles, independent of ground-based monitoring, providing a multidecadal archive suitable for probabilistic analysis linking deformation with eruption. Here we show that, of the 198 volcanoes systematically observed for the past 18 years, 54 deformed, of which 25 also erupted. For assessing eruption potential, this high proportion of deforming volcanoes that also erupted (46%), together with the proportion of non-deforming volcanoes that did not erupt (94%), jointly represent indicators with ‘strong’ evidential worth. Using a larger catalogue of 540 volcanoes observed for 3 years, we demonstrate how this eruption–deformation relationship is influenced by tectonic, petrological and volcanic factors. Satellite technology is rapidly evolving and routine monitoring of the deformation status of all volcanoes from space is anticipated, meaning probabilistic approaches will increasingly inform hazard decisions and strategic development
Micrometre-scale deformation observations reveal fundamental controls on geological rifting
Get PDFMany of the world’s largest volcanic eruptions are associated with geological rifting where major fractures open at the Earth’s surface, yet fundamental controls on the near-surface response to the rifting process are lacking. New high resolution observations gleaned from seismometer data during the 2014 Bárðarbunga basaltic dyke intrusion in Iceland allow us unprecedented access to the associated graben formation process on both sub-second and micrometre scales. We find that what appears as quasi steady-state near-surface rifting on lower resolution GPS observation comprises discrete staccato-like deformation steps as the upper crust unzips through repetitive low magnitude (M(W) < 0) failures on fracture patches estimated between 300 m(2) and 1200 m(2) in size. Stress drops for these events are one to two orders of magnitude smaller than expected for tectonic earthquakes, demonstrating that the uppermost crust in the rift zone is exceptionally weak
