The glacial isostatic adjustment signal at present day in northern Europe and the British Isles estimated from geodetic observations and geophysical models

Abstract. The glacial isostatic adjustment (GIA) signal at present day is constrained via the joint inversion of geodetic observations and GIA models for a region encompassing northern Europe, the British Isles, and the Barents Sea. The constraining data are Global Positioning System (GPS) vertical crustal velocities and GRACE (Gravity Recovery and Climate Experiment) gravity data. When the data are inverted with a set of GIA models, the best-fit model for the vertical motion signal has a χ² value of approximately 1 and a maximum a posteriori uncertainty of 0.3–0.4 mm yr⁻¹. An elastic correction is applied to the vertical land motion rates that accounts for present-day changes to terrestrial hydrology as well as recent mass changes of ice sheets and glaciered regions. Throughout the study area, mass losses from Greenland dominate the elastic vertical signal and combine to give an elastic correction of up to +0.5 mm yr⁻¹ in central Scandinavia. Neglecting to use an elastic correction may thus introduce a small but persistent bias in model predictions of GIA vertical motion even in central Scandinavia where vertical motion is dominated by GIA due to past glaciations. The predicted gravity signal is generally less well-constrained than the vertical signal, in part due to uncertainties associated with the correction for contemporary ice mass loss in Svalbard and the Russian Arctic. The GRACE-derived gravity trend is corrected for present-day ice mass loss using estimates derived from the ICESat and CryoSat missions, although a difference in magnitude between GRACE-inferred and altimetry-inferred regional mass loss rates suggests the possibility of a non-negligible GIA response here either from millennial-scale or Little Ice Age GIA.