Preprints
https://doi.org/10.5194/egusphere-2023-2869
https://doi.org/10.5194/egusphere-2023-2869
14 Dec 2023
 | 14 Dec 2023

FastIsostasy v1.0 – An accelerated regional GIA model accounting for the lateral variability of the solid Earth

Jan Swierczek-Jereczek, Marisa Montoya, Konstantin Latychev, Alexander Robinson, Jorge Alvarez-Solas, and Jerry Mitrovica

Abstract. The vast majority of ice-sheet modelling studies rely on simplified representations of the Glacial Isostatic Adjustment (GIA), which, among other limitations, do not account for lateral variations of the lithospheric thickness and upper-mantle viscosity. In studies using 3D GIA models, this has however been shown to have major impacts on the dynamics of marine-based sectors of Antarctica, which are likely to be the greatest contributors to sea-level rise in the coming centuries. This gap in comprehensiveness is explained by the fact that 3D GIA models are computationally expensive, seldomly open-source and require the implementation of an iterative coupling scheme to converge with the history of the ice-sheet model. To close this gap between "best" and "tractable" GIA models, we here propose FastIsostasy, a regional GIA model capturing lateral variations of the lithospheric thickness and mantle viscosity. By means of Fast-Fourier transforms and a hybrid collocation scheme to solve its underlying partial differential equation, FastIsostasy can simulate 100,000 years of high-resolution bedrock displacement in only minutes of single-CPU computation, including the changes in sea-surface height due to mass redistribution. Despite its 2D grid, FastIsostasy parametrises the depth-dependent viscosity in a physically meaningful way and therefore represents the depth dimension to a certain extent. FastIsostasy is here benchmarked against analytical, 1D and 3D GIA solutions and shows very good agreement with them. It is fully open-source, documented with many examples and provides a straight-forward interface for coupling to an ice-sheet model. The model is benchmarked here based on its implementation in Julia, while a Fortran version is also provided to allow for compatibility with most existing ice-sheet models. The Julia version provides additional features, including a vast library of time-stepping methods and GPU support.

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this preprint. The responsibility to include appropriate place names lies with the authors.
Share

Journal article(s) based on this preprint

10 Jul 2024
FastIsostasy v1.0 – a regional, accelerated 2D glacial isostatic adjustment (GIA) model accounting for the lateral variability of the solid Earth
Jan Swierczek-Jereczek, Marisa Montoya, Konstantin Latychev, Alexander Robinson, Jorge Alvarez-Solas, and Jerry Mitrovica
Geosci. Model Dev., 17, 5263–5290, https://doi.org/10.5194/gmd-17-5263-2024,https://doi.org/10.5194/gmd-17-5263-2024, 2024
Short summary
Download

The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.

Short summary
Ice sheets present a thickness of a few kilometres, leading to a vertical deformation of the...
Share