An approach for constraining mantle viscosities through assimilation of palaeo sea level data into a glacial isostatic adjustment model

Abstract

Abstract. Glacial isostatic adjustment is largely governed by the rheological properties of the Earth's mantle. Large mass redistributions in the ocean–cryosphere system and the subsequent response of the viscoelastic Earth have led to dramatic sea level changes in the past. This process is ongoing, and in order to understand and predict current and future sea level changes, the knowledge of mantle properties such as viscosity is essential. In this study, we present a method to obtain estimates of mantle viscosities by the assimilation of relative sea level rates of change into a viscoelastic model of the lithosphere and mantle. We set up a particle filter with probabilistic resampling. In an identical twin experiment, we show that mantle viscosities can be recovered in a glacial isostatic adjustment model of a simple three-layer Earth structure consisting of an elastic lithosphere and two mantle layers of different viscosity. We investigate the ensemble behaviour on different parameters in the following three set-ups: (1) global observations data set since last glacial maximum with different ensemble initialisations and observation uncertainties, (2) regional observations from Fennoscandia or Laurentide/Greenland only, and (3) limiting the observation period to 10 ka until the present. We show that the recovery is successful in all cases if the target parameter values are properly sampled by the initial ensemble probability distribution. This even includes cases in which the target viscosity values are located far in the tail of the initial ensemble probability distribution. Experiments show that the method is successful if enough near-field observations are available. This makes it work best for a period after substantial deglaciation until the present when the number of sea level indicators is relatively high.

Influence of GIA Uncertainty on Climate Model Evaluation With GRACE/GRACE‐FO Satellite Gravimetry Data

“…Our approach is closely related to the work of Pan et al (2022) who similarly utilized an ensemble of GIA models with horizontally and vertically variable Earth structures to assess in particular variations in global mean sea-level since the last glacial maximum. All 56 GIA simulations considered in this study are performed with the VILMA model, that has been applied in many GIA-related studies in the past (Bagge et al, 2021;Dobslaw et al, 2020;Hoening et al, 2023;Huang et al, 2023;Klemann et al, 2008Klemann et al, , 2015Martinec et al, 2018;Schachtschneider et al, 2022). The numerical code solves the sea level equation in the spherical domain and accounts for continuum-mechanics, viscoelasticity, rotational feedback, self-gravitation, mass conservation, and the migration of coastlines.…”

EAT v0.9.6: a 1D testbed for physical-biogeochemical data assimilation in natural waters

“…The data assimilation methods are implemented in a generic way, allowing PDAF to be applied in various modelling applications like ocean physics (Brüning et al, 2021) and sea ice (Mu et al, 2019), ocean biogeochemistry (Goodliff et al, 2019;Pradhan et al, 2019), hydrology (Kurtz et al, 2016), geodynamo 150 (Fournier et al, 2013) and geodynamics (Schachtschneider et al, 2022) or transport dynamics in the atmosphere (Pardini et al, 2020). PDAF is implemented in Fortran 95 with some functions of Fortran 2003 and uses parallelization with the Message Passing Interface (MPI, Gropp et al, 1994) and OpenMP.…”