Articles | Volume 16, issue 2
https://doi.org/10.5194/tc-16-689-2022
© Author(s) 2022. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/tc-16-689-2022
© Author(s) 2022. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
A comparison of the stability and performance of depth-integrated ice-dynamics solvers
Department of Earth Science and Astrophysics, Complutense University of Madrid, Madrid, Spain
Geosciences Institute CSIC-UCM, Madrid, Spain
Potsdam Institute for Climate Impact Research, Potsdam, Germany
Climate and Global Dynamics Laboratory, National Center for Atmospheric Research, Boulder, CO 80305, USA
Daniel Goldberg
School of GeoSciences, University of Edinburgh, Edinburgh, UK
William H. Lipscomb
Climate and Global Dynamics Laboratory, National Center for Atmospheric Research, Boulder, CO 80305, USA
Related authors
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.
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.
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
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.
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.
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.
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.
Tim van den Akker, William H. Lipscomb, Gunter R. Leguy, Willem Jan van de Berg, and Roderik S. W. van de Wal
EGUsphere, https://doi.org/10.5194/egusphere-2025-3380, https://doi.org/10.5194/egusphere-2025-3380, 2025
This preprint is open for discussion and under review for The Cryosphere (TC).
Short summary
Short summary
The Antarctic Ice Sheet is currently thinning, especially at major outlet glaciers. Including present-day thinning rates in models is a modeller's choice and can affect future projections. This study quantifies the impact of current imbalance on forced future projections, revealing strong regional and short-term (up to 2100) effects when these mass change rates are included.
Patrick Schmitt, Fabien Maussion, Daniel N. Goldberg, and Philipp Gregor
EGUsphere, https://doi.org/10.5194/egusphere-2025-3401, https://doi.org/10.5194/egusphere-2025-3401, 2025
This preprint is open for discussion and under review for Geoscientific Model Development (GMD).
Short summary
Short summary
To improve large-scale understanding of glaciers, we developed a new data assimilation method that integrates available observations in a dynamically consistent way, while taking their timestamps into account. It is designed to flexibly include new glacier data as it becomes available. We tested the method with idealized experiments and found promising results in terms of accuracy and efficiency, showing strong potential for real-world applications.
Laure Moinat, Florian Franziskakis, Christian Vérard, Daniel N. Goldberg, and Maura Brunetti
EGUsphere, https://doi.org/10.5194/egusphere-2025-2946, https://doi.org/10.5194/egusphere-2025-2946, 2025
This preprint is open for discussion and under review for Geoscientific Model Development (GMD).
Short summary
Short summary
We describe a new tool, biogeodyn-MITgcmIS, that consistently reproduces the global-scale dynamics of the ocean, atmosphere, vegetation and ice on multimillennial timescales at low computational cost. Evaluated against observations and state-of-the-art Earth system models, it includes offline coupling to models of vegetation, hydrology and a newly developed global-scale ice sheet. Using arbitrary continental configurations, it enables studies of past and present climates on Earth or exoplanets.
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.
Claire K. Yung, Xylar S. Asay-Davis, Alistair Adcroft, Christopher Y. S. Bull, Jan De Rydt, Michael S. Dinniman, Benjamin K. Galton-Fenzi, Daniel Goldberg, David E. Gwyther, Robert Hallberg, Matthew Harrison, Tore Hattermann, David M. Holland, Denise Holland, Paul R. Holland, James R. Jordan, Nicolas C. Jourdain, Kazuya Kusahara, Gustavo Marques, Pierre Mathiot, Dimitris Menemenlis, Adele K. Morrison, Yoshihiro Nakayama, Olga Sergienko, Robin S. Smith, Alon Stern, Ralph Timmermann, and Qin Zhou
EGUsphere, https://doi.org/10.5194/egusphere-2025-1942, https://doi.org/10.5194/egusphere-2025-1942, 2025
Short summary
Short summary
ISOMIP+ compares 12 ocean models that simulate ice-ocean interactions in a common, idealised, static ice shelf cavity setup, aiming to assess and understand inter-model variability. Models simulate similar basal melt rate patterns, ocean profiles and circulation but differ in ice-ocean boundary layer properties and spatial distributions of melting. Ice-ocean boundary layer representation is a key area for future work, as are realistic-domain ice sheet-ocean model intercomparisons.
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.
Jowan M. Barnes, G. Hilmar Gudmundsson, Daniel N. Goldberg, and Sainan Sun
EGUsphere, https://doi.org/10.5194/egusphere-2025-328, https://doi.org/10.5194/egusphere-2025-328, 2025
Short summary
Short summary
Calving is where ice breaks off the front of glaciers. It has not been included widely in modelling as it is difficult to represent. We use our ice flow model to investigate the effects of calving floating ice shelves in West Antarctica. More calving leads to more ice loss and greater sea level rise, with local differences due to the shape of the bedrock. We find that ocean forcing and calving should be considered equally when trying to improve how models represent the real world.
Tim van den Akker, William H. Lipscomb, Gunter R. Leguy, Willem Jan van de Berg, and Roderik S. W. van de Wal
EGUsphere, https://doi.org/10.5194/egusphere-2025-441, https://doi.org/10.5194/egusphere-2025-441, 2025
Short summary
Short summary
Ice sheet models to simulate future sea level rise require parameterizations, like for the friction at the bedrock. Studies have quantified the effect of using different parameterizations, and some have concluded that projections are sensitive to the choice of the specific parameterization. In this study, we show that you can make an ice sheet model sensitive to the basal friction parameterization, and that for equally defendable modellers choices you can also make the model insensitive to this.
William H. Lipscomb, David Behar, and Monica Ainhorn Morrison
The Cryosphere, 19, 793–803, https://doi.org/10.5194/tc-19-793-2025, https://doi.org/10.5194/tc-19-793-2025, 2025
Short summary
Short summary
As communities try to adapt to climate change, they look for "actionable science" that can inform decision-making. There are risks in relying on novel results that are not yet accepted by the science community. We propose a practical criterion for determining which scientific claims are actionable. We show how premature acceptance of sea-level-rise predictions can lead to confusion and backtracking, and we suggest best practices for communication between scientists and adaptation planners.
Samar Minallah, William Lipscomb, Gunter Leguy, and Harry Zekollari
EGUsphere, https://doi.org/10.5194/egusphere-2024-4152, https://doi.org/10.5194/egusphere-2024-4152, 2025
Short summary
Short summary
We implemented a new modeling framework within an Earth system model to study the evolution of mountain glaciers under different climate scenarios and applied it to the European Alps. Alpine glaciers will lose a large volume fraction under current temperatures, with near complete ice loss under warmer scenarios. This is the first use of a 3D, higher-order ice flow model for regional-scale glacier simulations that will enable assessments of coupled land ice and Earth system processes.
Tim van den Akker, William H. Lipscomb, Gunter R. Leguy, Jorjo Bernales, Constantijn J. Berends, Willem Jan van de Berg, and Roderik S. W. van de Wal
The Cryosphere, 19, 283–301, https://doi.org/10.5194/tc-19-283-2025, https://doi.org/10.5194/tc-19-283-2025, 2025
Short summary
Short summary
In this study, we present an improved way of representing ice thickness change rates in an ice sheet model. We apply this method using two ice sheet models of the Antarctic Ice Sheet. We found that the two largest outlet glaciers on the Antarctic Ice Sheet, Thwaites Glacier and Pine Island Glacier, will collapse without further warming on a timescale of centuries. This would cause a sea level rise of about 1.2 m globally.
Michele Petrini, Meike D. W. Scherrenberg, Laura Muntjewerf, Miren Vizcaino, Raymond Sellevold, Gunter R. Leguy, William H. Lipscomb, and Heiko Goelzer
The Cryosphere, 19, 63–81, https://doi.org/10.5194/tc-19-63-2025, https://doi.org/10.5194/tc-19-63-2025, 2025
Short summary
Short summary
Anthropogenic warming is causing accelerated Greenland ice sheet melt. Here, we use a computer model to understand how prolonged warming and ice melt could threaten ice sheet stability. We find a threshold beyond which Greenland will lose more than 80 % of its ice over several thousand years, due to the interaction of surface and solid-Earth processes. Nearly complete Greenland ice sheet melt occurs when the ice margin disconnects from a region of high elevation in western Greenland.
Heiko Goelzer, Petra M. Langebroek, Andreas Born, Stefan Hofer, Konstanze Haubner, Michele Petrini, Gunter Leguy, William H. Lipscomb, and Katherine Thayer-Calder
EGUsphere, https://doi.org/10.5194/egusphere-2024-3045, https://doi.org/10.5194/egusphere-2024-3045, 2025
Short summary
Short summary
On the backdrop of observed accelerating ice sheet mass loss over the last few decades, there is growing interest in the role of ice sheet changes in global climate projections. In this regard, we have coupled an Earth system model with an ice sheet model and have produced an initial set of climate projections including an interactive coupling with a dynamic Greenland ice sheet.
Mira Berdahl, Gunter R. Leguy, William H. Lipscomb, Bette L. Otto-Bliesner, Esther C. Brady, Robert A. Tomas, Nathan M. Urban, Ian Miller, Harriet Morgan, and Eric J. Steig
Clim. Past, 20, 2349–2371, https://doi.org/10.5194/cp-20-2349-2024, https://doi.org/10.5194/cp-20-2349-2024, 2024
Short summary
Short summary
Studying climate conditions near the Antarctic ice sheet (AIS) during Earth’s past warm periods informs us about how global warming may influence AIS ice loss. Using a global climate model, we investigate climate conditions near the AIS during the Last Interglacial (129 to 116 kyr ago), a period with warmer global temperatures and higher sea level than today. We identify the orbital and freshwater forcings that could cause ice loss and probe the mechanisms that lead to warmer climate conditions.
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.
Jean-François Lemieux, William H. Lipscomb, Anthony Craig, David A. Bailey, Elizabeth C. Hunke, Philippe Blain, Till A. S. Rasmussen, Mats Bentsen, Frédéric Dupont, David Hebert, and Richard Allard
Geosci. Model Dev., 17, 6703–6724, https://doi.org/10.5194/gmd-17-6703-2024, https://doi.org/10.5194/gmd-17-6703-2024, 2024
Short summary
Short summary
We present the latest version of the CICE model. It solves equations that describe the dynamics and the growth and melt of sea ice. To do so, the domain is divided into grid cells and variables are positioned at specific locations in the cells. A new implementation (C-grid) is presented, with the velocity located on cell edges. Compared to the previous B-grid, the C-grid allows for a natural coupling with some oceanic and atmospheric models. It also allows for ice transport in narrow channels.
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.
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
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.
David T. Bett, Alexander T. Bradley, C. Rosie Williams, Paul R. Holland, Robert J. Arthern, and Daniel N. Goldberg
The Cryosphere, 18, 2653–2675, https://doi.org/10.5194/tc-18-2653-2024, https://doi.org/10.5194/tc-18-2653-2024, 2024
Short summary
Short summary
A new ice–ocean model simulates future ice sheet evolution in the Amundsen Sea sector of Antarctica. Substantial ice retreat is simulated in all scenarios, with some retreat still occurring even with no future ocean melting. The future of small "pinning points" (islands of ice that contact the seabed) is an important control on this retreat. Ocean melting is crucial in causing these features to go afloat, providing the link by which climate change may affect this sector's sea level contribution.
Tong Zhang, William Colgan, Agnes Wansing, Anja Løkkegaard, Gunter Leguy, William H. Lipscomb, and Cunde Xiao
The Cryosphere, 18, 387–402, https://doi.org/10.5194/tc-18-387-2024, https://doi.org/10.5194/tc-18-387-2024, 2024
Short summary
Short summary
The geothermal heat flux determines how much heat enters from beneath the ice sheet, and thus impacts the temperature and the flow of the ice sheet. In this study we investigate how much geothermal heat flux impacts the initialization of the Greenland ice sheet. We use the Community Ice Sheet Model with two different initialization methods. We find a non-trivial influence of the choice of heat flow boundary conditions on the ice sheet initializations for further designs of ice sheet modeling.
Hélène Seroussi, Vincent Verjans, Sophie Nowicki, Antony J. Payne, Heiko Goelzer, William H. Lipscomb, Ayako Abe-Ouchi, Cécile Agosta, Torsten Albrecht, Xylar Asay-Davis, Alice Barthel, Reinhard Calov, Richard Cullather, Christophe Dumas, Benjamin K. Galton-Fenzi, Rupert Gladstone, Nicholas R. Golledge, Jonathan M. Gregory, Ralf Greve, Tore Hattermann, Matthew J. Hoffman, Angelika Humbert, Philippe Huybrechts, Nicolas C. Jourdain, Thomas Kleiner, Eric Larour, Gunter R. Leguy, Daniel P. Lowry, Chistopher M. Little, Mathieu Morlighem, Frank Pattyn, Tyler Pelle, Stephen F. Price, Aurélien Quiquet, Ronja Reese, Nicole-Jeanne Schlegel, Andrew Shepherd, Erika Simon, Robin S. Smith, Fiammetta Straneo, Sainan Sun, Luke D. Trusel, Jonas Van Breedam, Peter Van Katwyk, Roderik S. W. van de Wal, Ricarda Winkelmann, Chen Zhao, Tong Zhang, and Thomas Zwinger
The Cryosphere, 17, 5197–5217, https://doi.org/10.5194/tc-17-5197-2023, https://doi.org/10.5194/tc-17-5197-2023, 2023
Short summary
Short summary
Mass loss from Antarctica is a key contributor to sea level rise over the 21st century, and the associated uncertainty dominates sea level projections. We highlight here the Antarctic glaciers showing the largest changes and quantify the main sources of uncertainty in their future evolution using an ensemble of ice flow models. We show that on top of Pine Island and Thwaites glaciers, Totten and Moscow University glaciers show rapid changes and a strong sensitivity to warmer ocean conditions.
Beatriz Recinos, Daniel Goldberg, James R. Maddison, and Joe Todd
The Cryosphere, 17, 4241–4266, https://doi.org/10.5194/tc-17-4241-2023, https://doi.org/10.5194/tc-17-4241-2023, 2023
Short summary
Short summary
Ice sheet models generate forecasts of ice sheet mass loss, a significant contributor to sea level rise; thus, capturing the complete range of possible projections of mass loss is of critical societal importance. Here we add to data assimilation techniques commonly used in ice sheet modelling (a Bayesian inference approach) and fully characterize calibration uncertainty. We successfully propagate this type of error onto sea level rise projections of three ice streams in West Antarctica.
René R. Wijngaard, Adam R. Herrington, William H. Lipscomb, Gunter R. Leguy, and Soon-Il An
The Cryosphere, 17, 3803–3828, https://doi.org/10.5194/tc-17-3803-2023, https://doi.org/10.5194/tc-17-3803-2023, 2023
Short summary
Short summary
We evaluate the ability of the Community Earth System Model (CESM2) to simulate cryospheric–hydrological variables, such as glacier surface mass balance (SMB), over High Mountain Asia (HMA) by using a global grid (~111 km) with regional refinement (~7 km) over HMA. Evaluations of two different simulations show that climatological biases are reduced, and glacier SMB is improved (but still too negative) by modifying the snow and glacier model and using an updated glacier cover dataset.
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.
Constantijn J. Berends, Roderik S. W. van de Wal, Tim van den Akker, and William H. Lipscomb
The Cryosphere, 17, 1585–1600, https://doi.org/10.5194/tc-17-1585-2023, https://doi.org/10.5194/tc-17-1585-2023, 2023
Short summary
Short summary
The rate at which the Antarctic ice sheet will melt because of anthropogenic climate change is uncertain. Part of this uncertainty stems from processes occurring beneath the ice, such as the way the ice slides over the underlying bedrock.
Inversion methodsattempt to use observations of the ice-sheet surface to calculate how these sliding processes work. We show that such methods cannot fully solve this problem, so a substantial uncertainty still remains in projections of sea-level rise.
Mira Berdahl, Gunter Leguy, William H. Lipscomb, Nathan M. Urban, and Matthew J. Hoffman
The Cryosphere, 17, 1513–1543, https://doi.org/10.5194/tc-17-1513-2023, https://doi.org/10.5194/tc-17-1513-2023, 2023
Short summary
Short summary
Contributions to future sea level from the Antarctic Ice Sheet remain poorly constrained. One reason is that ice sheet model initialization methods can have significant impacts on how the ice sheet responds to future forcings. We investigate the impacts of two key parameters used during model initialization. We find that these parameter choices alone can impact multi-century sea level rise by up to 2 m, emphasizing the need to carefully consider these choices for sea level rise predictions.
Helen Ockenden, Robert G. Bingham, Andrew Curtis, and Daniel Goldberg
The Cryosphere, 16, 3867–3887, https://doi.org/10.5194/tc-16-3867-2022, https://doi.org/10.5194/tc-16-3867-2022, 2022
Short summary
Short summary
Hills and valleys hidden under the ice of Thwaites Glacier have an impact on ice flow and future ice loss, but there are not many three-dimensional observations of their location or size. We apply a mathematical theory to new high-resolution observations of the ice surface to predict the bed topography beneath the ice. There is a good correlation with ice-penetrating radar observations. The method may be useful in areas with few direct observations or as a further constraint for other methods.
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.
Adrian K. Turner, William H. Lipscomb, Elizabeth C. Hunke, Douglas W. Jacobsen, Nicole Jeffery, Darren Engwirda, Todd D. Ringler, and Jonathan D. Wolfe
Geosci. Model Dev., 15, 3721–3751, https://doi.org/10.5194/gmd-15-3721-2022, https://doi.org/10.5194/gmd-15-3721-2022, 2022
Short summary
Short summary
We present the dynamical core of the MPAS-Seaice model, which uses a mesh consisting of a Voronoi tessellation with polygonal cells. Such a mesh allows variable mesh resolution in different parts of the domain and the focusing of computational resources in regions of interest. We describe the velocity solver and tracer transport schemes used and examine errors generated by the model in both idealized and realistic test cases and examine the computational efficiency of the model.
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.
Conrad P. Koziol, Joe A. Todd, Daniel N. Goldberg, and James R. Maddison
Geosci. Model Dev., 14, 5843–5861, https://doi.org/10.5194/gmd-14-5843-2021, https://doi.org/10.5194/gmd-14-5843-2021, 2021
Short summary
Short summary
Sea level change due to the loss of ice sheets presents great risk for coastal communities. Models are used to forecast ice loss, but their evolution depends strongly on properties which are hidden from observation and must be inferred from satellite observations. Common methods for doing so do not allow for quantification of the uncertainty inherent or how it will affect forecasts. We provide a framework for quantifying how this
initialization uncertaintyaffects ice loss forecasts.
Gunter R. Leguy, William H. Lipscomb, and Xylar S. Asay-Davis
The Cryosphere, 15, 3229–3253, https://doi.org/10.5194/tc-15-3229-2021, https://doi.org/10.5194/tc-15-3229-2021, 2021
Short summary
Short summary
We present numerical features of the Community Ice Sheet Model in representing ocean termini glaciers. Using idealized test cases, we show that applying melt in a partly grounded cell is beneficial, in contrast to recent studies. We confirm that parameterizing partly grounded cells yields accurate ice sheet representation at a grid resolution of ~2 km (arguably 4 km), allowing ice sheet simulations at a continental scale. The choice of basal friction law also influences the ice flow.
Mira Berdahl, Gunter Leguy, William H. Lipscomb, and Nathan M. Urban
The Cryosphere, 15, 2683–2699, https://doi.org/10.5194/tc-15-2683-2021, https://doi.org/10.5194/tc-15-2683-2021, 2021
Short summary
Short summary
Antarctic ice shelves are vulnerable to warming ocean temperatures and have already begun thinning in response to increased basal melt rates. Sea level is expected to rise due to Antarctic contributions, but uncertainties in rise amount and timing remain largely unquantified. To facilitate uncertainty quantification, we use a high-resolution ice sheet model to build, test, and validate an ice sheet emulator and generate probabilistic sea level rise estimates for 100 and 200 years in the future.
Jowan M. Barnes, Thiago Dias dos Santos, Daniel Goldberg, G. Hilmar Gudmundsson, Mathieu Morlighem, and Jan De Rydt
The Cryosphere, 15, 1975–2000, https://doi.org/10.5194/tc-15-1975-2021, https://doi.org/10.5194/tc-15-1975-2021, 2021
Short summary
Short summary
Some properties of ice flow models must be initialised using observed data before they can be used to produce reliable predictions of the future. Different models have different ways of doing this, and the process is generally seen as being specific to an individual model. We compare the methods used by three different models and show that they produce similar outputs. We also demonstrate that the outputs from one model can be used in other models without introducing large uncertainties.
William H. Lipscomb, Gunter R. Leguy, Nicolas C. Jourdain, Xylar Asay-Davis, Hélène Seroussi, and Sophie Nowicki
The Cryosphere, 15, 633–661, https://doi.org/10.5194/tc-15-633-2021, https://doi.org/10.5194/tc-15-633-2021, 2021
Short summary
Short summary
This paper describes Antarctic climate change experiments in which the Community Ice Sheet Model is forced with ocean warming predicted by global climate models. Generally, ice loss begins slowly, accelerates by 2100, and then continues unabated, with widespread retreat of the West Antarctic Ice Sheet. The mass loss by 2500 varies from about 150 to 1300 mm of equivalent sea level rise, based on the predicted ocean warming and assumptions about how this warming drives melting beneath ice shelves.
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.
Heiko Goelzer, Sophie Nowicki, Anthony Payne, Eric Larour, Helene Seroussi, William H. Lipscomb, Jonathan Gregory, Ayako Abe-Ouchi, Andrew Shepherd, Erika Simon, Cécile Agosta, Patrick Alexander, Andy Aschwanden, Alice Barthel, Reinhard Calov, Christopher Chambers, Youngmin Choi, Joshua Cuzzone, Christophe Dumas, Tamsin Edwards, Denis Felikson, Xavier Fettweis, Nicholas R. Golledge, Ralf Greve, Angelika Humbert, Philippe Huybrechts, Sebastien Le clec'h, Victoria Lee, Gunter Leguy, Chris Little, Daniel P. Lowry, Mathieu Morlighem, Isabel Nias, Aurelien Quiquet, Martin Rückamp, Nicole-Jeanne Schlegel, Donald A. Slater, Robin S. Smith, Fiamma Straneo, Lev Tarasov, Roderik van de Wal, and Michiel van den Broeke
The Cryosphere, 14, 3071–3096, https://doi.org/10.5194/tc-14-3071-2020, https://doi.org/10.5194/tc-14-3071-2020, 2020
Short summary
Short summary
In this paper we use a large ensemble of Greenland ice sheet models forced by six different global climate models to project ice sheet changes and sea-level rise contributions over the 21st century.
The results for two different greenhouse gas concentration scenarios indicate that the Greenland ice sheet will continue to lose mass until 2100, with contributions to sea-level rise of 90 ± 50 mm and 32 ± 17 mm for the high (RCP8.5) and low (RCP2.6) scenario, respectively.
Hélène Seroussi, Sophie Nowicki, Antony J. Payne, Heiko Goelzer, William H. Lipscomb, Ayako Abe-Ouchi, Cécile Agosta, Torsten Albrecht, Xylar Asay-Davis, Alice Barthel, Reinhard Calov, Richard Cullather, Christophe Dumas, Benjamin K. Galton-Fenzi, Rupert Gladstone, Nicholas R. Golledge, Jonathan M. Gregory, Ralf Greve, Tore Hattermann, Matthew J. Hoffman, Angelika Humbert, Philippe Huybrechts, Nicolas C. Jourdain, Thomas Kleiner, Eric Larour, Gunter R. Leguy, Daniel P. Lowry, Chistopher M. Little, Mathieu Morlighem, Frank Pattyn, Tyler Pelle, Stephen F. Price, Aurélien Quiquet, Ronja Reese, Nicole-Jeanne Schlegel, Andrew Shepherd, Erika Simon, Robin S. Smith, Fiammetta Straneo, Sainan Sun, Luke D. Trusel, Jonas Van Breedam, Roderik S. W. van de Wal, Ricarda Winkelmann, Chen Zhao, Tong Zhang, and Thomas Zwinger
The Cryosphere, 14, 3033–3070, https://doi.org/10.5194/tc-14-3033-2020, https://doi.org/10.5194/tc-14-3033-2020, 2020
Short summary
Short summary
The Antarctic ice sheet has been losing mass over at least the past 3 decades in response to changes in atmospheric and oceanic conditions. This study presents an ensemble of model simulations of the Antarctic evolution over the 2015–2100 period based on various ice sheet models, climate forcings and emission scenarios. Results suggest that the West Antarctic ice sheet will continue losing a large amount of ice, while the East Antarctic ice sheet could experience increased snow accumulation.
Cited articles
Arthern, R. J. and Williams, C. R.: The sensitivity of West Antarctica to the
submarine melting feedback, Geophys. Res. Lett., 44, 2352–2359,
https://doi.org/10.1002/2017GL072514, 2017. a
Arthern, R. J., Hindmarsh, R. C. A., and Williams, C. R.: Flow speed within the
Antarctic ice sheet and its controls inferred from satellite observations, J. Geophys. Res.-Earth, 120, 1171–1188, https://doi.org/10.1002/2014JF003239, 2015. a, b
Blatter, H.: Velocity and stress fields in grounded glaciers – a simple
algorithm for including deviatoric stress gradients, J. Glaciol., 41,
333–344, 1995. a
Bueler, E.: Lectures at Karthaus: Numerical modelling of ice sheets and ice
shelves,
https://glaciers.gi.alaska.edu/sites/default/files/Notes_icesheetmod_Bueler2014.pdf (last access: 23 February 2022), 2009. a
Bueler, E. and Brown, J.: Shallow shelf approximation as a “sliding law” in a
thermodynamically coupled ice sheet model, J. Geophys. Res., 114, F03008,
https://doi.org/10.1029/2008JF001179, 2009. a
Dukowicz, J. K.: Reformulating the full-Stokes ice sheet model for a more efficient computational solution, The Cryosphere, 6, 21–34, https://doi.org/10.5194/tc-6-21-2012, 2012. a
Glen, J. W.: The creep of polycrystalline ice, Proc. R. Soc. London A, 228,
519–538, 1955. a
Goelzer, H., Nowicki, S., Edwards, T., Beckley, M., Abe-Ouchi, A., Aschwanden, A., Calov, R., Gagliardini, O., Gillet-Chaulet, F., Golledge, N. R., Gregory, J., Greve, R., Humbert, A., Huybrechts, P., Kennedy, J. H., Larour, E., Lipscomb, W. H., Le clec'h, S., Lee, V., Morlighem, M., Pattyn, F., Payne, A. J., Rodehacke, C., Rückamp, M., Saito, F., Schlegel, N., Seroussi, H., Shepherd, A., Sun, S., van de Wal, R., and Ziemen, F. A.: Design and results of the ice sheet model initialisation experiments initMIP-Greenland: an ISMIP6 intercomparison, The Cryosphere, 12, 1433–1460, https://doi.org/10.5194/tc-12-1433-2018, 2018. a, b, c
Goelzer, H., Nowicki, S., Payne, A., Larour, E., Seroussi, H., Lipscomb, W. H., Gregory, J., Abe-Ouchi, A., Shepherd, A., Simon, E., Agosta, C., Alexander, P., Aschwanden, A., Barthel, A., Calov, R., Chambers, C., Choi, Y., Cuzzone, J., Dumas, C., Edwards, T., Felikson, D., Fettweis, X., Golledge, N. R., Greve, R., Humbert, A., Huybrechts, P., Le clec'h, S., Lee, V., Leguy, G., Little, C., Lowry, D. P., Morlighem, M., Nias, I., Quiquet, A., Rückamp, M., Schlegel, N.-J., Slater, D. A., Smith, R. S., Straneo, F., Tarasov, L., van de Wal, R., and van den Broeke, M.: The future sea-level contribution of the Greenland ice sheet: a multi-model ensemble study of ISMIP6, The Cryosphere, 14, 3071–3096, https://doi.org/10.5194/tc-14-3071-2020, 2020. a
Greve, R. and Blatter, H.: Dynamics of Ice Sheets and Glaciers, 1st edn.,
Springer-Verlag, Berlin, https://doi.org/10.1007/978-3-642-03415-2, 2009. a, b
Hindmarsh, R.: A numerical comparison of approximations to the Stokes equations
used in ice sheet and glacier modeling, J. Geophys. Res., 109, F01012,
https://doi.org/10.1029/2003JF000065, 2004. a, b
Hindmarsh, R.: The role of membrane-like stresses in determining the stability
and sensitivity of the Antarctic ice sheets: back pressure and grounding line
motion, Philos. T. R. Soc. A, 364, 1733–1767,
https://doi.org/10.1098/rsta.2006.1797, 2006. a, b
Hindmarsh, R. C.: Notes on basic glaciological computational methods and
algorithms, in: Continuum mechanics and applications in geophysics and the
environment, edited by Straughan, B., Greve, R., Ehrentraut, H., and Wang,
Y., 222–249, Springer, Berling,
http://nora.nerc.ac.uk/id/eprint/19612/ (last access: 23 February 2022), 2001. a, b
Hutter, K.: Theoretical Glaciology, Mathematical Approaches to Geophysics, 1st edn., edited by: Rikitaki, T. and Hazenwinkel, M., D.
Reidel Publishing Company, Dordrecht, Boston, Lancaster, https://doi.org/10.1007/978-94-015-1167-4, 1983. a
Isaacson, E. and Keller, H. B.: Analysis of numerical methods, 1st edn., Courier
Corporation, ISBN 9780486680293, 2012. a
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, c, d
Lipscomb, W. H., Leguy, G. R., Jourdain, N. C., Asay-Davis, X., Seroussi, H., and Nowicki, S.: ISMIP6-based projections of ocean-forced Antarctic Ice Sheet evolution using the Community Ice Sheet Model, The Cryosphere, 15, 633–661, https://doi.org/10.5194/tc-15-633-2021, 2021. a
Lipscomb, W., Sacks, B., Price, S., Hoffman, M., Martin, D. F., Ranken, D., Hagdorn, M., Salinger, A., Evans, K., Norman, M., Shannon, S., Thayer-Calder, K., Leguy, G., Tezaur, I. K., Kennedy, J. H., Goelzer, H., Kluzek, E., Rutt, I. C., and Iesulauro Barker, E.: ESCOMP/CISM: CISM Slab Tests (cism_slab_tests), Zenodo [code], https://doi.org/10.5281/zenodo.5889016, 2022. a
MacAyeal, D. R.: Large-scale ice flow over a viscous basal sediment – Theory
and application to Ice Stream B, Antarctica, J. Geophys. Res., 94,
4071–4087, 1989. a
Morland, L.: Unconfined ice-shelf flow, in: Dynamics of the West Antarctic Ice Sheet, edited by: Van der Veen, C. J. and Oerlemans, J., Glac. Quat. G.,
4, 99–116, Springer, https://doi.org/10.1007/978-94-009-3745-1_6, 1987. a
Morland, L. W.: Thermomechanical balances of ice sheet flows, Geophys. Astro. Fluid, 29, 237–266, https://doi.org/10.1080/03091928408248191,
1984. a
Pattyn, F.: A new three-dimensional higher-order thermomechanical ice sheet
model: Basic sensitivity, ice stream development, and ice flow across
subglacial lakes, J. Geophys. Res., 108, 2382, https://doi.org/10.1029/2002JB002329,
2003. a
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
Pattyn, F., Perichon, L., Aschwanden, A., Breuer, B., de Smedt, B., Gagliardini, O., Gudmundsson, G. H., Hindmarsh, R. C. A., Hubbard, A., Johnson, J. V., Kleiner, T., Konovalov, Y., Martin, C., Payne, A. J., Pollard, D., Price, S., Rückamp, M., Saito, F., Souček, O., Sugiyama, S., and Zwinger, T.: Benchmark experiments for higher-order and full-Stokes ice sheet models (ISMIP–HOM), The Cryosphere, 2, 95–108, https://doi.org/10.5194/tc-2-95-2008, 2008. a, b, c, d
Pollard, D. and DeConto, R. M.: Description of a hybrid ice sheet-shelf model, and application to Antarctica, Geosci. Model Dev., 5, 1273–1295, https://doi.org/10.5194/gmd-5-1273-2012, 2012. 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, b
Robinson, A.: palma-ice/yelmo: Solver Stability v1.0, Zenodo [code], https://doi.org/10.5281/zenodo.5791864, 2021. 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, d
Seroussi, H., Nowicki, S., Simon, E., Abe-Ouchi, A., Albrecht, T., Brondex, J., Cornford, S., Dumas, C., Gillet-Chaulet, F., Goelzer, H., Golledge, N. R., Gregory, J. M., Greve, R., Hoffman, M. J., Humbert, A., Huybrechts, P., Kleiner, T., Larour, E., Leguy, G., Lipscomb, W. H., Lowry, D., Mengel, M., Morlighem, M., Pattyn, F., Payne, A. J., Pollard, D., Price, S. F., Quiquet, A., Reerink, T. J., Reese, R., Rodehacke, C. B., Schlegel, N.-J., Shepherd, A., Sun, S., Sutter, J., Van Breedam, J., van de Wal, R. S. W., Winkelmann, R., and Zhang, T.: initMIP-Antarctica: an ice sheet model initialization experiment of ISMIP6, The Cryosphere, 13, 1441–1471, https://doi.org/10.5194/tc-13-1441-2019, 2019. a, b
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
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.
Here we investigate the numerical stability of several commonly used methods in order to...