Temporal Gravity Variations near Shrinking Vatnajökull Ice Cap, Iceland

Repeated gravity measurements were carried out from 1991 until 1999 at sites SE of Vatnajökull, Iceland, to estimate the mass flow and deformation accompanying the shrinking of the ice cap. Published GPS data show an uplift of about 13 ± 5 mm/a near the ice margin. A gravity decrease of –2 ± 1 μGal/...

Full description

Bibliographic Details
Published in:Pure and Applied Geophysics
Main Authors: Jacoby, W., Hartmann, O., Wallner, H., Smilde, P., Bürger, S., Sjöberg, L., Erlingsson, S., Wolf, D., Klemann, V., Sasgen, I.
Other Authors: 1.3 Earth System Modelling, 1.0 Geodesy and Remote Sensing, Departments, GFZ Publication Database, Deutsches GeoForschungsZentrum
Format: Article in Journal/Newspaper
Language:unknown
Published: 2009
Subjects:
Online Access:https://gfzpublic.gfz-potsdam.de/pubman/item/item_239307
Description
Summary:Repeated gravity measurements were carried out from 1991 until 1999 at sites SE of Vatnajökull, Iceland, to estimate the mass flow and deformation accompanying the shrinking of the ice cap. Published GPS data show an uplift of about 13 ± 5 mm/a near the ice margin. A gravity decrease of –2 ± 1 μGal/a relative to the Höfn base station, was observed for the same sites. Control measurements at the Höfn station showed a gravity decrease of –2 ± 0.5 µGal/a relative to the station RVIK 5473 at Reykjavík (about 250 km from Höfn). This is compatible, as a Bouguer effect, with a 10 ± 3 mm/a uplift rate of the IGS point at Höfn and an uplift rate of ~20 mm/a near the ice margin. Although the derived gravity change rates at individual sites have large uncertainties, the ensemble of the rates varies systematically and significantly with distance from the ice. The relationship between gravity and elevation changes and the shrinking ice mass is modelled as response to the loading history. The GPS data can be explained by 1-D modelling (i.e., an earth model with a 15-km thick elastic lithosphere and a 7·10^17 Pa·s asthenosphere viscosity), but not the gravity data. Based on 2-D modelling, the gravity data favour a low-viscosity plume in the form of a cylinder of 80 km radius and 10^17 to 10^18 Pa·s viscosity below a 6 km-thick elastic lid, embedded in a layered PREM-type earth, although the elevation data are less well explained by this model. Strain-porosity-hydrology effects are likely to enhance the magnitude of the gravity changes, but need verification by drilling. More accurate data may resolve the discrepancies or suggest improved models.