A benchmark study of numerical implementations of the sea level equation in GIA modelling

Abstract

 Users may download and print one copy of any publication from the public portal for the purpose of private study or research.  You may not further distribute the material or use it for any profit-making activity or commercial gain  You may freely distribute the URL identifying the publication in the public portal If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim.

SELEN<sup>4</sup> (SELEN version 4.0): a Fortran program for solving the gravitationally and topographically self-consistent Sea Level Equation in Glacial Isostatic Adjustment modeling

“…Last, SELEN 4 comes with a User guide and with a fully detailed theory background, which is particularly meant to illustrate the basic concepts of GIA to young scientists or colleagues and to allow transparency and reproducibility. 25 For simplified surface loads, recently a preliminary version of the new program has been successfully tested against other independently developed, but not yet publicly available, SLE solvers (see Martinec et al, 2018). After it has been progressively tion, Geophysical Journal International, 194, 45-60, 2013.…”

“…We have also ascertained that misfits of the order of 1 mm yr −1 are not uncommon in other highlatitude sites of both hemispheres. Disclosing the origin of the discrepancies in the two sets of GIA predictions (and therefore in the whole set of the GIA geodetic fingerprints considered in Figure 8) is not easy at this stage, and would demand a detailed model inter-comparison study like those performed in the GIA community by Spada et al (2011) and Martinec et al (2018). 10 We can however guess that the misfit between the two sets of GIA predictions stems from the different discretisations of the ice time-histories, from the effects of mantle compressibility, and possibly from the different rotation theories adopted.…”

“…v) SELEN 4 no longer calls the post glacial rebound solver TABOO as an internal subroutine to compute the viscoelastic loading and tidal Love numbers, which are instead supplied by the user through a data file. In this way, any set of Love numbers can be used in 15 SELEN 4 , possibly overcoming some of the intrinsic limitations of existing Love numbers calculators like TABOO, vi) Recently, a prototype version of SELEN 4 has been successfully validated in a community benchmark of independently developed GIA codes (Martinec et al, 2018).…”

Reversal of the Direction of Horizontal Velocities Induced by GIA as a Function of Mantle Viscosity

“…The colatitude resolution and maximum spherical harmonic degree are the same as in the FE model. The Love number computation was benchmarked in Spada et al (2011), and the combination with surface loads in Martinec et al (2018).…”

SELEN<sup>4</sup> (SELEN version 4.0): a Fortran program for solving the gravitationally and topographically self-consistent sea-level equation in glacial isostatic adjustment modeling

“…The SLE solver ISSM-SESAW v1.0 (Ice Sheet System Model -Solid Earth and Sea-level Adjustment Workbench; Adhikari et al, 2016), being oriented to short-term cryosphere and climate changes, only accounts for the elastic deformation of the Earth. The open-source SLE solver giapy (Kachuck, 2017;Martinec et al, 2018), available from https://github. com/skachuck/giapy (last access: 26 November 2019), can deal with complex ice models and viscoelastic rheology but does not take rotational effects into account.…”