Geodynamic processes in Antarctica such as glacial isostatic adjustment (GIA) and post-seismic deformation are measured by geo-detic observations such as global navigation satellite systems (GNSS) and satellite gravimetry. GNSS measurements have comprised both continuous measurements and episodic measurements since the mid-1990s. The estimated velocities typically reach an accuracy of 1 mm a⁻¹ for horizontal velocities and 2 mm a⁻¹ for vertical velocities. However, the elastic deformation due to present-day ice-load change needs to be considered accordingly. Space gravimetry derives mass changes from small variations in the inter-satellite distance of a pair of satellites, starting with the GRACE (Gravity Recovery and Climate Experiment) satellite mission in 2002 and continuing with the GRACE-FO (GRACE Follow-On) mission launched in 2018. The spatial resolution of the measurements is low (about 300 km) but the measurement error is homogeneous across Ant-arctica. The estimated trends contain signals from ice-mass change, and local and global GIA signals. To combine the strengths of the individual datasets, statistical combinations of GNSS, GRACE and satellite altimetry data have been developed. These combinations rely on realistic error estimates and assumptions of snow density. Nevertheless, they capture signals that are missing from geodynamic forward models such as the large uplift in the Amundsen Sea sector caused by a low-viscous response to century-scale ice-mass changes. ... Read more

To fully exploit data from the Gravity Recovery and Climate Experiment (GRACE), we separate geophysical signals observed by GRACE in Antarctica by deriving high-spatial resolution maps for present-day glacial isostatic adjustment (GIA) and ice-mass changes with the least possible noise level. For this, we simultaneously (i) improve the postprocessing of gravity data and (ii) consistently combine them with high-resolution data from Ice Cloud and land Elevation Satellite altimeter (ICESat) and Regional Atmospheric Climate Model 2.3 (RACMO). We use GPS observations to discriminate between various candidate spatial patterns of vertical motions caused by GIA. The ICESat-RACMO combination determines the spatial resolution of estimated ice-mass changes. The results suggest the capability of the developed approach to retrieve the complex spatial pattern of present-day GIA, such as a pronounced subsidence in the proximity of the Kamb Ice Stream and pronounced uplift in the Amundsen Sea Sector. ... Read more

The main goal of this thesis involves the development of a refined methodology to
separate the mass change signals associated with glacial isostatic adjustment (GIA)
from those of surface ice/firn by exploiting the strengths of independent data sets,
such as those from gravimetry, altimetry, climate data, and others. To achieve this,
various research efforts were conducted addressing specific aspects of the methodology and subsequent data processing. This led to a number of new contributions to the topic, ... Read more

We present a high resolution model of the linear trend in the Earth’s mass variations based on DMT-2 (Delft Mass Transport model, release 2). DMT-2 was produced primarily from K-Band Ranging (KBR) data of the Gravity Recovery And Climate Experiment (GRACE). It comprises a time series of monthly solutions complete to spherical harmonic degree 120. A novel feature in its production was the accurate computation and incorporation of stochastic properties of coloured noise when processing KBR data. The unconstrained DMT-2 monthly solutions are used to estimate the linear trend together with a bias, as well as annual and semi-annual sinusoidal terms. The linear term is further processed with an anisotropic Wiener filter, which uses full noise and signal covariance matrices. Given the fact that noise in an unconstrained model of the trend is reduced substantially as compared to monthly solutions, the Wiener filter associated with the trend is much less aggressive compared to a Wiener filter applied to monthly solutions. Consequently, the trend estimate shows an enhanced spatial resolution. It allows signals in relatively small water bodies, such as Aral sea and Ladoga lake, to be detected. Over the ice sheets, it allows for a clear identification of signals associated with some outlet glaciers or their groups. We compare the obtained trend estimate with the ones from the CSR-RL05 model using (i) the same approach based on monthly noise covariance matrices and (ii) a commonly-used approach based on the DDK-filtered monthly solutions. We use satellite altimetry data as independent control data. The comparison demonstrates a high spatial resolution of the DMT-2 linear trend. We link this to the usage of high-accuracy monthly noise covariance matrices, which is due to an accurate computation and incorporation of coloured noise when processing KBR data. A preliminary comparison of the linear trend based on DMT-2 with that computed from GSFC global mascons v01 reveals, among other, a high concentration of the signal along the coast for both models in areas like the ice sheets, Gulf of Alaska, and Iceland. ... Read more

There has been considerable research in the literature focused on computing and forecasting sea-level changes in terms of constant trends or rates. The Antarctic ice sheet is one of the main contributors to sea-level change with highly uncertain rates of glacial thinning and accumulation. Geodetic observing systems such as the Gravity Recovery and Climate Experiment (GRACE) and the Global Positioning System (GPS) are routinely used to estimate these trends. In an effort to improve the accuracy and reliability of these trends, this study investigates a technique that allows the estimated rates, along with co-estimated seasonal components, to vary in time. For this, state space models are defined and then solved by a Kalman filter (KF). The reliable estimation of noise parameters is one of the main problems encountered when using a KF approach, which is solved by numerically optimizing likelihood. Since the optimization problem is non-convex, it is challenging to find an optimal solution. To address this issue, we limited the parameter search space using classical least-squares adjustment (LSA). In this context, we also tested the usage of inequality constraints by directly verifying whether they are supported by the data. The suggested technique for time-series analysis is expanded ... Read more