2026-06-01T12:26:28Zhttps://docta.ucm.es/rest/oai/request

oai:docta.ucm.es:20.500.14352/1161792025-01-28T01:14:01Zcom_20.500.14352_14col_20.500.14352_15

FastIsostasy v1.0 – a regional, accelerated 2D glacial isostatic adjustment (GIA) model accounting for the lateral variability of the solid Earth Swierczek-Jerecze, Jan Montoya Redondo, María Luisa Latychev, Konstantin Robinson, Alexander James Álvarez Solas, Jorge Mitrovica, Jerry 550.3 Marine-ice-sheet Finite-element approach Postglacial sea-level Rapid bedrock uplift Upper-mantle Antarctica Deformation Rheology Surface Validation Geofísica 2507 Geofísica 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 in the lithospheric thickness and upper-mantle viscosity. In studies of the last glacial cycle 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, rarely open-source and require a complex coupling scheme. To close this gap between “best” and “tractable” GIA models, we propose FastIsostasy here, a regional GIA model capturing lateral variations in 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 parameterises the depth-dependent viscosity and therefore represents the depth dimension to a certain extent. FastIsostasy is benchmarked here against analytical, as well as 1D and 3D numerical solutions, and shows good agreement with them. For a simulation of the last glacial cycle, its mean and maximal error over time and space respectively yield less than 5 % and 16 % compared to a 3D GIA model over the regional solution domain. FastIsostasy is open-source, is documented with many examples and provides a straightforward 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 adaptive time-stepping methods and GPU support. CriticalEarth Unión Europea Ministerio de Ciencia, Innovación y Universidades (España) Depto. de Física de la Tierra y Astrofísica Fac. de Ciencias Físicas TRUE pub 2025-01-27T09:18:58Z 2025-01-27T09:18:58Z 2024-07-10 journal article VoR https://hdl.handle.net/20.500.14352/116179 1991-959X 10.5194/gmd-17-5263-2024 1991-9603 eng info:eu-repo/grantAgreement/EC/H2020/956170/EU info:eu-repo/grantAgreement/EU//101044247/EU info:eu-repo/grantAgreement/AEI/Plan Estatal de Investigación Científica y Técnica y de Innovación 2017-2020/PID2020-117768RB-I00/ES/MODELIZACION DE LAS INTERACCIONES MARINAS Y DE LA PERDIDA DINAMICA DE HIELO EN LA ANTARTIDA/ Swierczek-Jereczek, J., Montoya, M., Latychev, K., Robinson, A., Alvarez-Solas, J., and Mitrovica, J.: FastIsostasy v1.0 – a regional, accelerated 2D glacial isostatic adjustment (GIA) model accounting for the lateral variability of the solid Earth, Geosci. Model Dev., 17, 5263–5290, https://doi.org/10.5194/gmd-17-5263-2024, 2024. Attribution 4.0 International http://creativecommons.org/licenses/by/4.0/ open access application/pdf Copernicus Publications