Decadal geodetic variations in Ny-Ålesund (Svalbard): role of past and present ice-mass changes

The geodetic rates for the gravity variation and vertical uplift in polar regions subject to past and present-day ice-mass changes (PDIMCs) provide important insight into the rheological structure of the Earth. We provide an update of the rates observed at Ny-Ålesund, Svalbard. To do so, we extract and remove the significant seasonal content from the observations. The rate of gravity variations, derived from absolute and relative gravity measurements, is −1.39 ± 0.11 μGal yr⁻¹. The rate of vertical displacements is estimated using GPS and tide gauge measurements. We obtain 7.94 ± 0.21 and 8.29 ± 1.60 mm yr⁻¹, respectively. We compare the extracted signal with that predicted by GLDAS/Noah and ERA-interim hydrology models. We find that the seasonal gravity variations are well-represented by local hydrology changes contained in the ERA-interim model. The phase of seasonal vertical displacements are due to non-local continental hydrology and non-tidal ocean loading. However, a large part of the amplitude of the seasonal vertical displacements remains unexplained. The geodetic rates are used to investigate the asthenosphere viscosity and lithosphere/asthenosphere thicknesses. We first correct the updated geodetic rates for those induced by PDIMCs in Svalbard, using published results, and the sea level change due to the melting of the major ice reservoirs. We show that the latter are at the level of the geodetic rate uncertainties and are responsible for rates of gravity variations and vertical displacements of −0.29 ± 0.03 μGal yr⁻¹ and 1.11 ± 0.10 mm yr⁻¹, respectively. To account for the late Pleistocene deglaciation, we use the global ice evolution model ICE-3G. The Little Ice Age (LIA) deglaciation in Svalbard is modelled using a disc load model with a simple linear temporal evolution. The geodetic rates at Ny-Ålesund induced by the past deglaciations depend on the viscosity structure of the Earth. We find that viscous relaxation time due to the LIA deglaciation in Svalbard is more than 60 times shorter than that due to the Pleistocene deglaciation. We also find that the response to past and PDIMCs of an Earth model with asthenosphere viscosities ranging between 1.0 and 5.5 × 10¹⁸ Pa s and lithosphere (resp. asthenosphere) thicknesses ranging between 50 and 100 km (resp. 120 and 170 km) can explain the rates derived from geodetic observations.