In August 2021, Askja volcano in Iceland returned to the spotlight after a sudden onset of rapid uplift followed decades of continuous subsidence. In this study the extended record of microgravity data from Askja between 1988 to 2017 is revisited, and new microgravity data from 2021 and 2022 are introduced, which were collected after the uplift had started. Askja caldera had been steadily subsiding since at least 1984 and was characterised by a net decrease in microgravity, potentially signalling the contraction of its magma chamber or eviction of magma either laterally or to deeper levels. The microgravity data indicate that despite ongoing subsidence between 2017 and early 2021, a significant gravity increase can be detected in the center of the caldera between 2017 and August 2021. This increase may be introduced during – or leading up to – the period of uplift. The new microgravity data also indicate that during the period of 40 cm uplift after August 2021 to fall 2022, gravity changes approach the free-air gradient, suggesting subsurface density decreases as a driving process. This process may relate to the vesiculation of magma previously emplaced in the volcano roots, a change in the hydrothermal system, or replacement of dense basaltic magma with less dense rhyolitic magma, or a combination of these processes. However, uncertainties for this period are elevated and may obscure a gravity signal expected from additional mass accumulation. The timing and high uncertainties of some campaigns make it challenging to be conclusive on the driving process behind the uplift, but future microgravity campaigns could help solve the ambiguity. The study also provides a description of potential pitfalls in microgravity campaigns and recommendations on how the reliability of microgravity data can be improved.

The 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 20km 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 4km depth closing by 2-4cm 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 (∼ 4km by ∼ 2km) while it is circular at Krafla (∼ 5km 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.5km 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.