Articles | Volume 17, issue 13
https://doi.org/10.5194/gmd-17-5263-2024
© Author(s) 2024. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/gmd-17-5263-2024
© Author(s) 2024. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Development and technical paper
| 
10 Jul 2024
Development and technical paper |
| 10 Jul 2024
| 10 Jul 2024
FastIsostasy v1.0 – a regional, accelerated 2D glacial isostatic adjustment (GIA) model accounting for the lateral variability of the solid Earth
Department of Earth Physics and Astrophysics, Complutense University of Madrid, Madrid, Spain
Geosciences Institute, CSIC–UCM, Madrid, Spain
Marisa Montoya
Department of Earth Physics and Astrophysics, Complutense University of Madrid, Madrid, Spain
Geosciences Institute, CSIC–UCM, Madrid, Spain
Konstantin Latychev
Seakon, Toronto, Canada
Alexander Robinson
Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research, Potsdam, Germany
Jorge Alvarez-Solas
Department of Earth Physics and Astrophysics, Complutense University of Madrid, Madrid, Spain
Geosciences Institute, CSIC–UCM, Madrid, Spain
Jerry Mitrovica
Department of Earth and Planetary Sciences, Harvard University, Cambridge, Massachusetts, USA
Related authors
Lucía Gutiérrez-González, Alexander Robinson, Jorge Alvarez-Solas, Ilaria Tabone, Jan Swierczek-Jereczek, Daniel Moreno-Parada, and Marisa Montoya
EGUsphere, https://doi.org/10.5194/egusphere-2025-2616, https://doi.org/10.5194/egusphere-2025-2616, 2025
This preprint is open for discussion and under review for The Cryosphere (TC).
Short summary
Short summary
The Greenland ice sheet is considered a tipping element: if temperatures exceed its threshold, it would transition to a virtually ice-free state. We analyze its stability across the full range of glacial-interglacial temperatures, as well as those expected in the coming centuries. We find a future critical threshold between +1.5 and +2 K of global warming, another under colder climates, and persistent hysteresis across the full range of study.
Sergio Pérez-Montero, Jorge Alvarez-Solas, Jan Swierczek-Jereczek, Daniel Moreno-Parada, Alexander Robinson, and Marisa Montoya
Earth Syst. Dynam., 16, 915–937, https://doi.org/10.5194/esd-16-915-2025, https://doi.org/10.5194/esd-16-915-2025, 2025
Short summary
Short summary
The climate of the last 3 Myr has varied between cold and warm periods. Numerous independent mechanisms have been proposed to explain this; however, no effort has been made to study their competing effects. Here we present a simple but physically motivated model that includes these mechanisms in a modular way. We identify ice-sheet dynamics and lithosphere displacement as main triggers, but reproducing the climate records additionally requires the natural darkening of ice.
Sergio Pérez-Montero, Jorge Alvarez-Solas, Jan Swierczek-Jereczek, Daniel Moreno-Parada, Alexander Robinson, and Marisa Montoya
EGUsphere, https://doi.org/10.5194/egusphere-2025-2467, https://doi.org/10.5194/egusphere-2025-2467, 2025
Short summary
Short summary
Almost 3 million years ago, the planet began to experience a succession of cold and warm periods every 40,000 years. However, about 1 million years ago, they began to occur every 100,000 years. In this paper we explore how the change in the basal velocity of the ice sheets could have produced this change in behavior. On the other hand, we also see that in our model, decreasing in time the sensitivity of snowfall to temperature is also an effective mechanism with which to reproduce the records.
Nils Bochow, Philipp Hess, and Alexander Robinson
EGUsphere, https://doi.org/https://doi.org/10.48550/arXiv.2507.22485, https://doi.org/https://doi.org/10.48550/arXiv.2507.22485, 2025
This preprint is open for discussion and under review for The Cryosphere (TC).
Short summary
Short summary
This study presents a fast, physics-guided machine-learning method that downscales coarse climate fields to fine resolution while enforcing conservation of large-scale totals. Trained on regional climate simulations and driven by Earth system model output, it handles extremes and outperforms linear interpolation, providing realistic, high-resolution forcing for ice-sheet models and improving projections of Greenland’s sea-level contribution.
Lucía Gutiérrez-González, Alexander Robinson, Jorge Alvarez-Solas, Ilaria Tabone, Jan Swierczek-Jereczek, Daniel Moreno-Parada, and Marisa Montoya
EGUsphere, https://doi.org/10.5194/egusphere-2025-2616, https://doi.org/10.5194/egusphere-2025-2616, 2025
This preprint is open for discussion and under review for The Cryosphere (TC).
Short summary
Short summary
The Greenland ice sheet is considered a tipping element: if temperatures exceed its threshold, it would transition to a virtually ice-free state. We analyze its stability across the full range of glacial-interglacial temperatures, as well as those expected in the coming centuries. We find a future critical threshold between +1.5 and +2 K of global warming, another under colder climates, and persistent hysteresis across the full range of study.
Daniel Moreno-Parada, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas
Geosci. Model Dev., 18, 3895–3919, https://doi.org/10.5194/gmd-18-3895-2025, https://doi.org/10.5194/gmd-18-3895-2025, 2025
Short summary
Short summary
We introduce Nix, an ice-sheet model designed for understanding how large masses of ice behave. Nix is a computer programme that simulates the movement and temperature evolution in ice sheets. It helps us study how ice sheets respond to changes in the atmosphere and ocean. We found that ice temperatures play an essential role in determining the motion and stability of ice sheets. Nix is a useful tool for learning how climate change affects polar ice sheets.
Sergio Pérez-Montero, Jorge Alvarez-Solas, Jan Swierczek-Jereczek, Daniel Moreno-Parada, Alexander Robinson, and Marisa Montoya
Earth Syst. Dynam., 16, 915–937, https://doi.org/10.5194/esd-16-915-2025, https://doi.org/10.5194/esd-16-915-2025, 2025
Short summary
Short summary
The climate of the last 3 Myr has varied between cold and warm periods. Numerous independent mechanisms have been proposed to explain this; however, no effort has been made to study their competing effects. Here we present a simple but physically motivated model that includes these mechanisms in a modular way. We identify ice-sheet dynamics and lithosphere displacement as main triggers, but reproducing the climate records additionally requires the natural darkening of ice.
Ricarda Winkelmann, Donovan P. Dennis, Jonathan F. Donges, Sina Loriani, Ann Kristin Klose, Jesse F. Abrams, Jorge Alvarez-Solas, Torsten Albrecht, David Armstrong McKay, Sebastian Bathiany, Javier Blasco Navarro, Victor Brovkin, Eleanor Burke, Gokhan Danabasoglu, Reik V. Donner, Markus Drüke, Goran Georgievski, Heiko Goelzer, Anna B. Harper, Gabriele Hegerl, Marina Hirota, Aixue Hu, Laura C. Jackson, Colin Jones, Hyungjun Kim, Torben Koenigk, Peter Lawrence, Timothy M. Lenton, Hannah Liddy, José Licón-Saláiz, Maxence Menthon, Marisa Montoya, Jan Nitzbon, Sophie Nowicki, Bette Otto-Bliesner, Francesco Pausata, Stefan Rahmstorf, Karoline Ramin, Alexander Robinson, Johan Rockström, Anastasia Romanou, Boris Sakschewski, Christina Schädel, Steven Sherwood, Robin S. Smith, Norman J. Steinert, Didier Swingedouw, Matteo Willeit, Wilbert Weijer, Richard Wood, Klaus Wyser, and Shuting Yang
EGUsphere, https://doi.org/10.5194/egusphere-2025-1899, https://doi.org/10.5194/egusphere-2025-1899, 2025
This preprint is open for discussion and under review for Earth System Dynamics (ESD).
Short summary
Short summary
The Tipping Points Modelling Intercomparison Project (TIPMIP) is an international collaborative effort to systematically assess tipping point risks in the Earth system using state-of-the-art coupled and stand-alone domain models. TIPMIP will provide a first global atlas of potential tipping dynamics, respective critical thresholds and key uncertainties, generating an important building block towards a comprehensive scientific basis for policy- and decision-making.
Sergio Pérez-Montero, Jorge Alvarez-Solas, Jan Swierczek-Jereczek, Daniel Moreno-Parada, Alexander Robinson, and Marisa Montoya
EGUsphere, https://doi.org/10.5194/egusphere-2025-2467, https://doi.org/10.5194/egusphere-2025-2467, 2025
Short summary
Short summary
Almost 3 million years ago, the planet began to experience a succession of cold and warm periods every 40,000 years. However, about 1 million years ago, they began to occur every 100,000 years. In this paper we explore how the change in the basal velocity of the ice sheets could have produced this change in behavior. On the other hand, we also see that in our model, decreasing in time the sensitivity of snowfall to temperature is also an effective mechanism with which to reproduce the records.
Meghan E. King, Jessica R. Creveling, and Jerry X. Mitrovica
Clim. Past, 21, 53–65, https://doi.org/10.5194/cp-21-53-2025, https://doi.org/10.5194/cp-21-53-2025, 2025
Short summary
Short summary
In this study, we compute glacial–interglacial sea-level changes across the mid-Pliocene Warm Period (MPWP; 3.264–3.025 Ma) resulting from ice mass loss from different ice sheets. Our results quantify the relationship between changes in local sea level and global mean sea level (GMSL) and highlight the level of consistency in this mapping across different ice melt scenarios. These predictions can help guide site selection in efforts to constrain the sources and magnitude of MPWP GMSL change.
Ryan Love, Glenn A. Milne, Parviz Ajourlou, Soran Parang, Lev Tarasov, and Konstantin Latychev
Geosci. Model Dev., 17, 8535–8551, https://doi.org/10.5194/gmd-17-8535-2024, https://doi.org/10.5194/gmd-17-8535-2024, 2024
Short summary
Short summary
A relatively recent advance in glacial isostatic adjustment modeling has been the development of models that include 3D Earth structure, as opposed to 1D structure. However, a major limitation is the computational expense. We have developed a method using artificial neural networks to emulate the influence of 3D Earth models to affordably constrain the viscosity parameter space. Our results indicate that the misfits are of a scale such that useful predictions of relative sea level can be made.
Antonio Juarez-Martinez, Javier Blasco, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas
The Cryosphere, 18, 4257–4283, https://doi.org/10.5194/tc-18-4257-2024, https://doi.org/10.5194/tc-18-4257-2024, 2024
Short summary
Short summary
We present sea level projections for Antarctica in the context of ISMIP6-2300 with several forcings but extend the simulations to 2500, showing that more than 3 m of sea level contribution could be reached. We also test the sensitivity on a basal melting parameter and determine the timing of the loss of ice in the west region. All the simulations were carried out with the ice sheet model Yelmo.
Therese Rieckh, Andreas Born, Alexander Robinson, Robert Law, and Gerrit Gülle
Geosci. Model Dev., 17, 6987–7000, https://doi.org/10.5194/gmd-17-6987-2024, https://doi.org/10.5194/gmd-17-6987-2024, 2024
Short summary
Short summary
We present the open-source model ELSA, which simulates the internal age structure of large ice sheets. It creates layers of snow accumulation at fixed times during the simulation, which are used to model the internal stratification of the ice sheet. Together with reconstructed isochrones from radiostratigraphy data, ELSA can be used to assess ice sheet models and to improve their parameterization. ELSA can be used coupled to an ice sheet model or forced with its output.
Daniel Moreno-Parada, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas
The Cryosphere, 18, 4215–4232, https://doi.org/10.5194/tc-18-4215-2024, https://doi.org/10.5194/tc-18-4215-2024, 2024
Short summary
Short summary
Our study tries to understand how the ice temperature evolves in a large mass as in the case of Antarctica. We found a relation that tells us the ice temperature at any point. These results are important because they also determine how the ice moves. In general, ice moves due to slow deformation (as if pouring honey from a jar). Nevertheless, in some regions the ice base warms enough and melts. The liquid water then serves as lubricant and the ice slides and its velocity increases rapidly.
Javier Blasco, Ilaria Tabone, Daniel Moreno-Parada, Alexander Robinson, Jorge Alvarez-Solas, Frank Pattyn, and Marisa Montoya
Clim. Past, 20, 1919–1938, https://doi.org/10.5194/cp-20-1919-2024, https://doi.org/10.5194/cp-20-1919-2024, 2024
Short summary
Short summary
In this study, we assess Antarctic tipping points which may had been crossed during the mid-Pliocene Warm Period. For this, we use data from the PlioMIP2 ensemble. Additionally, we investigate various sources of uncertainty, like ice dynamics and bedrock configuration. Our research significantly enhances our comprehension of Antarctica's response to a warming climate, shedding light on potential future tipping points that may be surpassed.
Sönke Dangendorf, Qiang Sun, Thomas Wahl, Philip Thompson, Jerry X. Mitrovica, and Ben Hamlington
Earth Syst. Sci. Data, 16, 3471–3494, https://doi.org/10.5194/essd-16-3471-2024, https://doi.org/10.5194/essd-16-3471-2024, 2024
Short summary
Short summary
Sea-level information from the global ocean is sparse in time and space, with comprehensive data being limited to the period since 2005. Here we provide a novel reconstruction of sea level and its contributing causes, as determined by a Kalman smoother approach applied to tide gauge records over the period 1900–2021. The new reconstruction shows a continuing acceleration in global mean sea-level rise since 1970 that is dominated by melting land ice. Contributors vary significantly by region.
Natasha Valencic, Linda Pan, Konstantin Latychev, Natalya Gomez, Evelyn Powell, and Jerry X. Mitrovica
The Cryosphere, 18, 2969–2978, https://doi.org/10.5194/tc-18-2969-2024, https://doi.org/10.5194/tc-18-2969-2024, 2024
Short summary
Short summary
We quantify the effect of ongoing Antarctic bedrock uplift due to Ice Age or modern ice mass changes on estimates of ice thickness changes obtained from satellite-based ice height measurements. We find that variations in the Ice Age signal introduce an uncertainty in estimates of total Antarctic ice change of up to ~10%. Moreover, the usual assumption that the mapping between modern ice height and thickness changes is uniform systematically underestimates net Antarctic ice volume changes.
Daniel Moreno-Parada, Jorge Alvarez-Solas, Javier Blasco, Marisa Montoya, and Alexander Robinson
The Cryosphere, 17, 2139–2156, https://doi.org/10.5194/tc-17-2139-2023, https://doi.org/10.5194/tc-17-2139-2023, 2023
Short summary
Short summary
We have reconstructed the Laurentide Ice Sheet, located in North America during the Last Glacial Maximum (21 000 years ago). The absence of direct measurements raises a number of uncertainties. Here we study the impact of different physical laws that describe the friction as the ice slides over its base. We found that the Laurentide Ice Sheet is closest to prior reconstructions when the basal friction takes into account whether the base is frozen or thawed during its motion.
Matteo Willeit, Andrey Ganopolski, Alexander Robinson, and Neil R. Edwards
Geosci. Model Dev., 15, 5905–5948, https://doi.org/10.5194/gmd-15-5905-2022, https://doi.org/10.5194/gmd-15-5905-2022, 2022
Short summary
Short summary
In this paper we present the climate component of the newly developed fast Earth system model CLIMBER-X. It has a horizontal resolution of 5°x5° and is designed to simulate the evolution of the Earth system on temporal scales ranging from decades to >100 000 years. CLIMBER-X is available as open-source code and is expected to be a useful tool for studying past climate changes and for the investigation of the long-term future evolution of the climate.
Alexander Robinson, Daniel Goldberg, and William H. Lipscomb
The Cryosphere, 16, 689–709, https://doi.org/10.5194/tc-16-689-2022, https://doi.org/10.5194/tc-16-689-2022, 2022
Short summary
Short summary
Here we investigate the numerical stability of several commonly used methods in order to determine which of them are capable of resolving the complex physics of the ice flow and are also computationally efficient. We find that the so-called DIVA solver outperforms the others. Its representation of the physics is consistent with more complex methods, while it remains computationally efficient at high resolution.
Andreas Born and Alexander Robinson
The Cryosphere, 15, 4539–4556, https://doi.org/10.5194/tc-15-4539-2021, https://doi.org/10.5194/tc-15-4539-2021, 2021
Short summary
Short summary
Ice penetrating radar reflections from the Greenland ice sheet are the best available record of past accumulation and how these layers have been deformed over time by the flow of ice. Direct simulations of this archive hold great promise for improving our models and for uncovering details of ice sheet dynamics that neither models nor data could achieve alone. We present the first three-dimensional ice sheet model that explicitly simulates individual layers of accumulation and how they deform.
Javier Blasco, Jorge Alvarez-Solas, Alexander Robinson, and Marisa Montoya
The Cryosphere, 15, 215–231, https://doi.org/10.5194/tc-15-215-2021, https://doi.org/10.5194/tc-15-215-2021, 2021
Short summary
Short summary
During the Last Glacial Maximum the Antarctic Ice Sheet was larger and more extended than at present. However, neither its exact position nor the total ice volume are well constrained. Here we investigate how the different climatic boundary conditions, as well as basal friction configurations, affect the size and extent of the Antarctic Ice Sheet and discuss its potential implications.
Cited articles
A, G., Wahr, J., and Zhong, S.: Computations of the viscoelastic response of a 3-D compressible Earth to surface loading: an application to Glacial Isostatic Adjustment in Antarctica and Canada, Geophys. J. Int., 192, 557–572, https://doi.org/10.1093/gji/ggs030, 2013. a
Accardo, N. J., Wiens, D. A., Hernandez, S., Aster, R. C., Nyblade, A., Huerta, A., Anandakrishnan, S., Wilson, T., Heeszel, D. S., and Dalziel, I. W. D.: Upper mantle seismic anisotropy beneath the West Antarctic Rift System and surrounding region from shear wave splitting analysis, Geophys. J. Int., 198, 414–429, https://doi.org/10.1093/gji/ggu117, 2014. a
Adhikari, S., Ivins, E. R., Larour, E., Seroussi, H., Morlighem, M., and Nowicki, S.: Future Antarctic bed topography and its implications for ice sheet dynamics, Solid Earth, 5, 569–584, https://doi.org/10.5194/se-5-569-2014, 2014. a
Adhikari, S., Ivins, E. R., Larour, E., Caron, L., and Seroussi, H.: A kinematic formalism for tracking ice–ocean mass exchange on the Earth's surface and estimating sea-level change, The Cryosphere, 14, 2819–2833, https://doi.org/10.5194/tc-14-2819-2020, 2020. a
Amante, C. and Eakins, B. E.: ETOPO1 1 arc-minute global relief model: procedure, data sources and analysis, Tech. rep., NOAA, https://doi.org/10.7289/V5C8276M, 2009. a
Amelung, F. and Wolf, D.: Viscoelastic perturbations of the earth: significance of the incremental gravitational force in models of glacial isostasy, Geophys. J. Int., 117, 864–879, https://doi.org/10.1111/j.1365-246X.1994.tb02476.x, 1994. a
Argus, D. F., Peltier, W. R., Drummond, R., and Moore, A. W.: The Antarctica component of postglacial rebound model ICE-6G_C (VM5a) based on GPS positioning, exposure age dating of ice thicknesses, and relative sea level histories, Geophys. J. Int., 198, 537–563, https://doi.org/10.1093/gji/ggu140, 2014. a
Austermann, J., Hoggard, M. J., Latychev, K., Richards, F. D., and Mitrovica, J. X.: The effect of lateral variations in Earth structure on Last Interglacial sea level, Geophys. J. Int., 227, 1938–1960, https://doi.org/10.1093/gji/ggab289, 2021. a
Bagge, M., Klemann, V., Steinberger, B., Latinović, M., and Thomas, M.: Glacial‐Isostatic Adjustment Models Using Geodynamically Constrained 3D Earth Structures, Geochem. Geophy. Geosy., 22, e2021GC009853, https://doi.org/10.1029/2021GC009853, 2021.Please provide article number or page range. a
Barletta, V. R., Bevis, M., Smith, B. E., Wilson, T., Brown, A., Bordoni, A., Willis, M., Khan, S. A., Rovira-Navarro, M., Dalziel, I., Smalley, R., Kendrick, E., Konfal, S., Caccamise, D. J., Aster, R. C., Nyblade, A., and Wiens, D. A.: Observed rapid bedrock uplift in Amundsen Sea Embayment promotes ice-sheet stability, Science, 360, 1335–1339, https://doi.org/10.1126/science.aao1447, 2018. a, b
Beghein, C., Trampert, J., and van Heijst, H. J.: Radial anisotropy in seismic reference models of the mantle: ANISOTROPY IN MANTLE MODELS, J. Geophys. Res.-Sol. Ea., 111, B02303, https://doi.org/10.1029/2005JB003728, 2006. a
Behrendt, J. C.: Crustal and lithospheric structure of the West Antarctic Rift System from geophysical investigations – a review, Global Planet. Change, 23, 25–44, https://doi.org/10.1016/S0921-8181(99)00049-1, 1999. a
Besard, T., Foket, C., and De Sutter, B.: Effective Extensible Programming: Unleashing Julia on GPUs, IEEE T. Parall. Distr., 30, 827–841, https://doi.org/10.1109/TPDS.2018.2872064, 2019. a
Bezanson, J., Edelman, A., Karpinski, S., and Shah, V. B.: Julia: A Fresh Approach to Numerical Computing, SIAM Rev., 59, 65–98, https://doi.org/10.1137/141000671, 2017. a
Bogacki, P. and Shampine, L.: An efficient Runge-Kutta (4,5) pair, Comput. Math. Appl., 32, 15–28, https://doi.org/10.1016/0898-1221(96)00141-1, 1996. a
Book, C., Hoffman, M. J., Kachuck, S. B., Hillebrand, T. R., Price, S. F., Perego, M., and Bassis, J. N.: Stabilizing effect of bedrock uplift on retreat of Thwaites Glacier, Antarctica, at centennial timescales, Earth Planet. Sc. Lett., 597, 117798, https://doi.org/10.1016/j.epsl.2022.117798, 2022. a, b, c
Caron, L., Métivier, L., Greff-Lefftz, M., Fleitout, L., and Rouby, H.: Inverting Glacial Isostatic Adjustment signal using Bayesian framework and two linearly relaxing rheologies, Geophys. J. Int., 209, 1126–1147, https://doi.org/10.1093/gji/ggx083, 2017. a
Cogley, J. G.: Area of the Ocean, Mar. Geod., 35, 379–388, https://doi.org/10.1080/01490419.2012.709476, 2012. a
Coulon, V., Bulthuis, K., Whitehouse, P. L., Sun, S., Haubner, K., Zipf, L., and Pattyn, F.: Contrasting Response of West and East Antarctic Ice Sheets to Glacial Isostatic Adjustment, J. Geophys. Res.-Earth, 126, e2020JF006003, https://doi.org/10.1029/2020JF006003, 2021. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s
De Boer, B., Stocchi, P., Whitehouse, P. L., and Van De Wal, R. S.: Current state and future perspectives on coupled ice-sheet – sea-level modelling, Quaternary Sci. Rev., 169, 13–28, https://doi.org/10.1016/j.quascirev.2017.05.013, 2017. a
DeConto, R. M. and Pollard, D.: Contribution of Antarctica to past and future sea-level rise, Nature, 531, 591–597, https://doi.org/10.1038/nature17145, 2016. a
Farrell, W. E.: Deformation of the Earth by surface loads, Rev. Geophys., 10, 761, https://doi.org/10.1029/RG010i003p00761, 1972. a
Farrell, W. E. and Clark, J. A.: On Postglacial Sea Level, Geophys. J. Roy. Astron. Soc., 46, 647–667, https://doi.org/10.1111/j.1365-246X.1976.tb01252.x, 1976. a
Fretwell, P., Pritchard, H. D., Vaughan, D. G., Bamber, J. L., Barrand, N. E., Bell, R., Bianchi, C., Bingham, R. G., Blankenship, D. D., Casassa, G., Catania, G., Callens, D., Conway, H., Cook, A. J., Corr, H. F. J., Damaske, D., Damm, V., Ferraccioli, F., Forsberg, R., Fujita, S., Gim, Y., Gogineni, P., Griggs, J. A., Hindmarsh, R. C. A., Holmlund, P., Holt, J. W., Jacobel, R. W., Jenkins, A., Jokat, W., Jordan, T., King, E. C., Kohler, J., Krabill, W., Riger-Kusk, M., Langley, K. A., Leitchenkov, G., Leuschen, C., Luyendyk, B. P., Matsuoka, K., Mouginot, J., Nitsche, F. O., Nogi, Y., Nost, O. A., Popov, S. V., Rignot, E., Rippin, D. M., Rivera, A., Roberts, J., Ross, N., Siegert, M. J., Smith, A. M., Steinhage, D., Studinger, M., Sun, B., Tinto, B. K., Welch, B. C., Wilson, D., Young, D. A., Xiangbin, C., and Zirizzotti, A.: Bedmap2: improved ice bed, surface and thickness datasets for Antarctica, The Cryosphere, 7, 375–393, https://doi.org/10.5194/tc-7-375-2013, 2013. a
Frigo, M. and Johnson, S.: The Design and Implementation of FFTW3, P. IEEE, 93, 216–231, https://doi.org/10.1109/JPROC.2004.840301, 2005. a
Garbe, J., Albrecht, T., Levermann, A., Donges, J. F., and Winkelmann, R.: The hysteresis of the Antarctic Ice Sheet, Nature, 585, 538–544, https://doi.org/10.1038/s41586-020-2727-5, 2020. a
Gasperini, P., Dal Forno, G., and Boschi, E.: Linear or non-linear rheology in the Earth's mantle: the prevalence of power-law creep in the postglacial isostatic readjustment of Laurentia, Geophys. J. Int., 157, 1297–1302, https://doi.org/10.1111/j.1365-246X.2004.02319.x, 2004. a
Goelzer, H., Coulon, V., Pattyn, F., de Boer, B., and van de Wal, R.: Brief communication: On calculating the sea-level contribution in marine ice-sheet models, The Cryosphere, 14, 833–840, https://doi.org/10.5194/tc-14-833-2020, 2020. a, b, c, d
Gomez, N., Mitrovica, J. X., Huybers, P., and Clark, P. U.: Sea level as a stabilizing factor for marine-ice-sheet grounding lines, Nat. Geosci., 3, 850–853, https://doi.org/10.1038/ngeo1012, 2010. a, b, c
Gomez, N., Pollard, D., Mitrovica, J. X., Huybers, P., and Clark, P. U.: Evolution of a coupled marine ice sheet-sea level model, J. Geophys. Res.-Earth, 117, F01013, https://doi.org/10.1029/2011JF002128, 2012. a, b
Gomez, N., Pollard, D., and Holland, D.: Sea-level feedback lowers projections of future Antarctic Ice-Sheet mass loss, Nat. Commun., 6, 8798, https://doi.org/10.1038/ncomms9798, 2015. a, b
Gregory, J. M., Griffies, S. M., Hughes, C. W., Lowe, J. A., Church, J. A., Fukimori, I., Gomez, N., Kopp, R. E., Landerer, F., Cozannet, G. L., Ponte, R. M., Stammer, D., Tamisiea, M. E., and Van De Wal, R. S. W.: Concepts and Terminology for Sea Level: Mean, Variability and Change, Both Local and Global, Surv. Geophys., 40, 1251–1289, https://doi.org/10.1007/s10712-019-09525-z, 2019. a
Gudmundsson, G. H., Krug, J., Durand, G., Favier, L., and Gagliardini, O.: The stability of grounding lines on retrograde slopes, The Cryosphere, 6, 1497–1505, https://doi.org/10.5194/tc-6-1497-2012, 2012. a
Heeszel, D. S., Wiens, D. A., Anandakrishnan, S., Aster, R. C., Dalziel, I. W. D., Huerta, A. D., Nyblade, A. A., Wilson, T. J., and Winberry, J. P.: Upper mantle structure of central and West Antarctica from array analysis of Rayleigh wave phase velocities, J. Geophys. Res.-Sol. Ea., 121, 1758–1775, https://doi.org/10.1002/2015JB012616, 2016. a, b
Huang, P., Steffen, R., Steffen, H., Klemann, V., Wu, P., van der Wal, W., Martinec, Z., and Tanaka, Y.: A commercial finite element approach to modelling Glacial Isostatic Adjustment on spherical self-gravitating compressible earth models, Geophys. J. Int., 235, 2231–2256, https://doi.org/10.1093/gji/ggad354, 2023. a, b
Ivins, E. R., Caron, L., Adhikari, S., and Larour, E.: Notes on a compressible extended Burgers model of rheology, Geophys. J. Int., 228, 1975–1991, https://doi.org/10.1093/gji/ggab452, 2021. a
Ivins, E. R., van der Wal, W., Wiens, D. A., Lloyd, A. J., and Caron, L.: Antarctic upper mantle rheology, Geological Society, London, Memoirs, 56, M56–2020–19, https://doi.org/10.1144/M56-2020-19, 2022. a, b, c
Kachuck, S. B.: giapy: Glacial Isostatic Adjustment in PYthon (1.0.0), https://github.com/skachuck/giapy/ (last access: 4 July 2024), 2017. a
Kachuck, S. B., Martin, D. F., Bassis, J. N., and Price, S. F.: Rapid Viscoelastic Deformation Slows Marine Ice Sheet Instability at Pine Island Glacier, Geophys. Res. Lett., 47, e2019GL086446, https://doi.org/10.1029/2019GL086446, 2020. a, b, c
Kang, K., Zhong, S., Geruo, A., and Mao, W.: The effects of non-Newtonian rheology in the upper mantle on relative sea level change and geodetic observables induced by glacial isostatic adjustment process, Geophys. J. Int., 228, 1887–1906, https://doi.org/10.1093/gji/ggab428, 2021. a
Kendall, R. A., Mitrovica, J. X., and Milne, G. A.: On post-glacial sea level – II. Numerical formulation and comparative results on spherically symmetric models, Geophys. J. Int., 161, 679–706, https://doi.org/10.1111/j.1365-246X.2005.02553.x, 2005. a, b
Klemann, V., Martinec, Z., and Ivins, E. R.: Glacial isostasy and plate motion, J. Geodynam., 46, 95–103, https://doi.org/10.1016/j.jog.2008.04.005, 2008. a
Konrad, H., Thoma, M., Sasgen, I., Klemann, V., Grosfeld, K., Barbi, D., and Martinec, Z.: The Deformational Response of a Viscoelastic Solid Earth Model Coupled to a Thermomechanical Ice Sheet Model, Surv. Geophys., 35, 1441–1458, https://doi.org/10.1007/s10712-013-9257-8, 2014. a, b, c
Konrad, H., Sasgen, I., Pollard, D., and Klemann, V.: Potential of the solid-Earth response for limiting long-term West Antarctic Ice Sheet retreat in a warming climate, Earth Planet. Sc. Lett., 432, 254–264, https://doi.org/10.1016/j.epsl.2015.10.008, 2015. a
Konrad, H., Sasgen, I., Klemann, V., Thoma, M., Grosfeld, K., and Martinec, Z.: Sensitivity of grounding-line dynamics to viscoelastic deformation of the solid-earth in an idealized scenario, Polarforschung, Bremerhaven, Alfred Wegener Institute for Polar and Marine Research & German Society of Polar Research, 85, 89–99, https://doi.org/10.2312/POLFOR.2016.005, 2016. a
Kreuzer, M., Albrecht, T., Nicola, L., Reese, R., and Winkelmann, R.: Oceanic gateways in Antarctica – Impact of relative sea-level change on sub-shelf melt, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2023-2737, 2023. a
Kulp, S. A. and Strauss, B. H.: New elevation data triple estimates of global vulnerability to sea-level rise and coastal flooding, Nat. Commun., 10, 4844, https://doi.org/10.1038/s41467-019-12808-z, 2019. a
Lingle, C. S. and Clark, J. A.: A numerical model of interactions between a marine ice sheet and the solid earth: Application to a West Antarctic ice stream, J. Geophys. Res., 90, 1100, https://doi.org/10.1029/JC090iC01p01100, 1985. a, b, c, d
Lipscomb, W. H., Price, S. F., Hoffman, M. J., Leguy, G. R., Bennett, A. R., Bradley, S. L., Evans, K. J., Fyke, J. G., Kennedy, J. H., Perego, M., Ranken, D. M., Sacks, W. J., Salinger, A. G., Vargo, L. J., and Worley, P. H.: Description and evaluation of the Community Ice Sheet Model (CISM) v2.1, Geosci. Model Dev., 12, 387–424, https://doi.org/10.5194/gmd-12-387-2019, 2019. a, b
Lloyd, A. J., Wiens, D. A., Nyblade, A. A., Anandakrishnan, S., Aster, R. C., Huerta, A. D., Wilson, T. J., Dalziel, I. W. D., Shore, P. J., and Zhao, D.: A seismic transect across West Antarctica: Evidence for mantle thermal anomalies beneath the Bentley Subglacial Trench and the Marie Byrd Land Dome, J. Geophys. Res.-Sol. Ea., 120, 8439–8460, https://doi.org/10.1002/2015JB012455, 2015. a
Lloyd, A. J., Wiens, D. A., Zhu, H., Tromp, J., Nyblade, A. A., Aster, R. C., Hansen, S. E., Dalziel, I. W. D., Wilson, T., Ivins, E. R., and O’Donnell, J. P.: Seismic Structure of the Antarctic Upper Mantle Imaged with Adjoint Tomography, J. Geophys. Res.-Sol. Ea., 125, 2019JB017823, https://doi.org/10.1029/2019JB017823, 2020. a
Martinec, Z.: Spectral-finite element approach to three-dimensional viscoelastic relaxation in a spherical earth, Geophys. J. Int., 142, 117–141, https://doi.org/10.1046/j.1365-246x.2000.00138.x, 2000. a
Martinec, Z., Klemann, V., van der Wal, W., Riva, R. E. M., Spada, G., Sun, Y., Melini, D., Kachuck, S. B., Barletta, V., Simon, K., A, G., and James, T. S.: A benchmark study of numerical implementations of the sea level equation in GIA modelling, Geophys. J. Int., 215, 389–414, https://doi.org/10.1093/gji/ggy280, 2018. a
Mitrovica, J. X. and Milne, G. A.: On post-glacial sea level: I. General theory, Geophys. J. Int., 154, 253–267, https://doi.org/10.1046/j.1365-246X.2003.01942.x, 2003. a, b
Mitrovica, J. X., Milne, G. A., and Davis, J. L.: Glacial isostatic adjustment on a rotating earth, Geophys. J. Int., 147, 562–578, https://doi.org/10.1046/j.1365-246x.2001.01550.x, 2001. a
Mitrovica, J. X., Gomez, N., and Clark, P. U.: The Sea-Level Fingerprint of West Antarctic Collapse, Science, 323, 753–753, https://doi.org/10.1126/science.1166510, 2009. a
Morelli, A. and Danesi, S.: Seismological imaging of the Antarctic continental lithosphere: a review, Global Planet. Change, 42, 155–165, https://doi.org/10.1016/j.gloplacha.2003.12.005, 2004. a
Morlighem, M., Rignot, E., Binder, T., Blankenship, D., Drews, R., Eagles, G., Eisen, O., Ferraccioli, F., Forsberg, R., Fretwell, P., Goel, V., Greenbaum, J. S., Gudmundsson, H., Guo, J., Helm, V., Hofstede, C., Howat, I., Humbert, A., Jokat, W., Karlsson, N. B., Lee, W. S., Matsuoka, K., Millan, R., Mouginot, J., Paden, J., Pattyn, F., Roberts, J., Rosier, S., Ruppel, A., Seroussi, H., Smith, E. C., Steinhage, D., Sun, B., Broeke, M. R. V. D., Ommen, T. D. V., Wessem, M. V., and Young, D. A.: Deep glacial troughs and stabilizing ridges unveiled beneath the margins of the Antarctic ice sheet, Nat. Geosci., 13, 132–137, https://doi.org/10.1038/s41561-019-0510-8, 2020. a
Nield, G. A., Barletta, V. R., Bordoni, A., King, M. A., Whitehouse, P. L., Clarke, P. J., Domack, E., Scambos, T. A., and Berthier, E.: Rapid bedrock uplift in the Antarctic Peninsula explained by viscoelastic response to recent ice unloading, Earth Planet. Sc. Lett., 397, 32–41, https://doi.org/10.1016/j.epsl.2014.04.019, 2014. a, b
Nield, G. A., Whitehouse, P. L., van der Wal, W., Blank, B., O'Donnell, J. P., and Stuart, G. W.: The impact of lateral variations in lithospheric thickness on glacial isostatic adjustment in West Antarctica, Geophys. J. Int., 214, 811–824, https://doi.org/10.1093/gji/ggy158, 2018. a, b, c
Pan, L., Powell, E. M., Latychev, K., Mitrovica, J. X., Creveling, Gomez, N., Hoggard, M. J., and Clark, P. J.: Rapid postglacial rebound amplifies global sea level rise following West Antarctic Ice Sheet collapse, Sci. Adv., 7, eabf7787, https://doi.org/10.1126/sciadv.abf7787, 2021. a
Pan, L., Milne, G. A., Latychev, K., Goldberg, S. L., Austermann, J., Hoggard, M. J., and Mitrovica, J. X.: The influence of lateral Earth structure on inferences of global ice volume during the Last Glacial Maximum, Quaternary Sci. Rev., 290, 107644, https://doi.org/10.1016/j.quascirev.2022.107644, 2022. a, b, c, d, e
Pattyn, F.: Sea-level response to melting of Antarctic ice shelves on multi-centennial timescales with the fast Elementary Thermomechanical Ice Sheet model (f.ETISh v1.0), The Cryosphere, 11, 1851–1878, https://doi.org/10.5194/tc-11-1851-2017, 2017. a
Peltier, W. R., Argus, D. F., and Drummond, R.: Space geodesy constrains ice age terminal deglaciation: The global ICE‐6G_C (VM5a) model, J. Geophys. Res.-Sol. Ea., 120, 450–487, https://doi.org/10.1002/2014JB011176, 2015. a
Peltier, W. R., Argus, D. F., and Drummond, R.: Comment on “An Assessment of the ICE‐6G_C (VM5a) Glacial Isostatic Adjustment Model” by Purcell et al., J. Geophys. Res.-Sol. Ea., 123, 2019–2028, https://doi.org/10.1002/2016JB013844, 2018. a, b, c
Purcell, A.: The significance of pre-stress advection and internal buoyancy in the flat-Earth formulation, in: Dynamics of the Ice Age Earth – A Modern Perspective, Trans Tech Publications LTD, 24, 105–122, ISBN 978-0-87849-810-9, 1998. a
Purcell, A., Tregoning, P., and Dehecq, A.: An assessment of the ICE6G_C(VM5a) glacial isostatic adjustment model, J. Geophys. Res.-Sol. Ea., 121, 3939–3950, https://doi.org/10.1002/2015JB012742, 2016. a
Quiquet, A., Dumas, C., Ritz, C., Peyaud, V., and Roche, D. M.: The GRISLI ice sheet model (version 2.0): calibration and validation for multi-millennial changes of the Antarctic ice sheet, Geosci. Model Dev., 11, 5003–5025, https://doi.org/10.5194/gmd-11-5003-2018, 2018. a
Rackauckas, C. and Nie, Q.: DifferentialEquations.jl – A Performant and Feature-Rich Ecosystem for Solving Differential Equations in Julia, J. Open Res. Softw., 5, 15, https://doi.org/10.5334/jors.151, 2017. a
Robinson, A., Alvarez-Solas, J., Montoya, M., Goelzer, H., Greve, R., and Ritz, C.: Description and validation of the ice-sheet model Yelmo (version 1.0), Geosci. Model Dev., 13, 2805–2823, https://doi.org/10.5194/gmd-13-2805-2020, 2020. a, b, c
Rückamp, M., Greve, R., and Humbert, A.: Comparative simulations of the evolution of the Greenland ice sheet under simplified Paris Agreement scenarios with the models SICOPOLIS and ISSM, Polar Sci., 21, 14–25, https://doi.org/10.1016/j.polar.2018.12.003, 2019. a
Sasgen, I., Martín-Español, A., Horvath, A., Klemann, V., Petrie, E. J., Wouters, B., Horwath, M., Pail, R., Bamber, J. L., Clarke, P. J., Konrad, H., Wilson, T., and Drinkwater, M. R.: Altimetry, gravimetry, GPS and viscoelastic modeling data for the joint inversion for glacial isostatic adjustment in Antarctica (ESA STSE Project REGINA), Earth Syst. Sci. Data, 10, 493–523, https://doi.org/10.5194/essd-10-493-2018, 2018. a
Seroussi, H., Nowicki, S., Payne, A. J., Goelzer, H., Lipscomb, W. H., Abe-Ouchi, A., Agosta, C., Albrecht, T., Asay-Davis, X., Barthel, A., Calov, R., Cullather, R., Dumas, C., Galton-Fenzi, B. K., Gladstone, R., Golledge, N. R., Gregory, J. M., Greve, R., Hattermann, T., Hoffman, M. J., Humbert, A., Huybrechts, P., Jourdain, N. C., Kleiner, T., Larour, E., Leguy, G. R., Lowry, D. P., Little, C. M., Morlighem, M., Pattyn, F., Pelle, T., Price, S. F., Quiquet, A., Reese, R., Schlegel, N.-J., Shepherd, A., Simon, E., Smith, R. S., Straneo, F., Sun, S., Trusel, L. D., Van Breedam, J., van de Wal, R. S. W., Winkelmann, R., Zhao, C., Zhang, T., and Zwinger, T.: ISMIP6 Antarctica: a multi-model ensemble of the Antarctic ice sheet evolution over the 21st century, The Cryosphere, 14, 3033–3070, https://doi.org/10.5194/tc-14-3033-2020, 2020. a
Spada, G. and Melini, D.: SELEN4 (SELEN version 4.0): a Fortran program for solving the gravitationally and topographically self-consistent sea-level equation in glacial isostatic adjustment modeling, Geosci. Model Dev., 12, 5055–5075, https://doi.org/10.5194/gmd-12-5055-2019, 2019. a, b, c, d
Spada, G., Barletta, V. R., Klemann, V., Riva, R. E., Martinec, Z., Gasperini, P., Lund, B., Wolf, D., Vermeersen, L. L. A., and King, M. A.: A benchmark study for glacial isostatic adjustment codes: A GIA benchmark study, Geophys. J. Int., 185, 106–132, https://doi.org/10.1111/j.1365-246X.2011.04952.x, 2011. a, b, c, d, e, f, g, h, i
Swierczek-Jereczek, J.: JanJereczek/FastIsostasy.jl: v1.0.0, Zenodo [code], https://doi.org/10.5281/zenodo.11277347, 2024a. a
Swierczek-Jereczek, J.: JanJereczek/isostasy_data: Version 0.2 (v0.2.0). Zenodo [data set], https://doi.org/10.5281/zenodo.11175418, 2024b. a
Swierczek-Jereczek, J.: Simulating the GIA response of Antarctica during the last glacial cycle with FastIsosatsy, GitHub [video], https://github.com/JanJereczek/FastIsostasy.j, last access: 4 July 2024c. a
van Calcar, C. J., van de Wal, R. S. W., Blank, B., de Boer, B., and van der Wal, W.: Simulation of a fully coupled 3D glacial isostatic adjustment – ice sheet model for the Antarctic ice sheet over a glacial cycle, Geosci. Model Dev., 16, 5473–5492, https://doi.org/10.5194/gmd-16-5473-2023, 2023. a, b, c
van der Wal, W., Wu, P., Wang, H., and Sideris, M. G.: Sea levels and uplift rate from composite rheology in glacial isostatic adjustment modeling, J. Geodynam., 50, 38–48, https://doi.org/10.1016/j.jog.2010.01.006, 2010. a
van der Wal, W., Whitehouse, P. L., and Schrama, E. J. O.: Effect of GIA models with 3D composite mantle viscosity on GRACE mass balance estimates for Antarctica, Earth Planet. Sc. Lett., 414, 134–143, https://doi.org/10.1016/j.epsl.2015.01.001, 2015. a, b
Weerdesteijn, M. F. M., Naliboff, J. B., Conrad, C. P., Reusen, J. M., Steffen, R., Heister, T., and Zhang, J.: Modeling Viscoelastic Solid Earth Deformation Due To Ice Age and Contemporary Glacial Mass Changes in ASPECT, Geochem. Geophy. Geosy., 24, e2022GC010813, https://doi.org/10.1029/2022GC010813, 2023. a, b, c, d, e
Whitehouse, P. L., Gomez, N., King, M. A., and Wiens, D. A.: Solid Earth change and the evolution of the Antarctic Ice Sheet, Nat.e Commun., 10, 503, https://doi.org/10.1038/s41467-018-08068-y, 2019. a, b, c, d
Wiens, D. A., Shen, W., and Lloyd, A. J.: The seismic structure of the Antarctic upper mantle, Geological Society, London, Memoirs, 56, M56–2020–18, https://doi.org/10.1144/M56-2020-18, 2022. a
Winkelmann, R., Martin, M. A., Haseloff, M., Albrecht, T., Bueler, E., Khroulev, C., and Levermann, A.: The Potsdam Parallel Ice Sheet Model (PISM-PIK) – Part 1: Model description, The Cryosphere, 5, 715–726, https://doi.org/10.5194/tc-5-715-2011, 2011. a
Wu, P.: Using commercial finite element packages for the study of earth deformations, sea levels and the state of stress, Geophys. J. Int., 158, 401–408, https://doi.org/10.1111/j.1365-246X.2004.02338.x, 2004. a
Wu, P. and Wang, H.: Effects of mode coupling and location of rotational axis on glacial induced rotational deformation in a laterally heterogeneous viscoelastic earth, Geophys. J. Int., 167, 853–859, https://doi.org/10.1111/j.1365-246X.2006.03103.x, 2006. a
Zhong, S., Kang, K., A, G., and Qin, C.: CitcomSVE: A Three‐Dimensional Finite Element Software Package for Modeling Planetary Mantle’s Viscoelastic Deformation in Response to Surface and Tidal Loads, Geochem. Geophy. Geosy., 23, e2022GC010359, https://doi.org/10.1029/2022GC010359, 2022. a, b, c, d, e, f, g
Short summary
Ice sheets present a thickness of a few kilometres, leading to a vertical deformation of the crust of up to a kilometre. This process depends on properties of the solid Earth, which can be regionally very different. We propose a model that accounts for this often-ignored heterogeneity and run 100 000 simulation years in minutes. Thus, the evolution of ice sheets is modeled with better accuracy, which is critical for a good mitigation of climate change and, in particular, sea-level rise.
Ice sheets present a thickness of a few kilometres, leading to a vertical deformation of the...
