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.
Research article
| 
25 Feb 2022
Research article |
| 25 Feb 2022
| 25 Feb 2022
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
Therese Rieckh, Andreas Born, Alexander Robinson, Robert Law, and Gerrit Gülle
EGUsphere, https://doi.org/10.5194/egusphere-2023-3127, https://doi.org/10.5194/egusphere-2023-3127, 2024
Short summary
Short summary
We present the open-source model ELSA, which simulates the internal age structure of large ice sheets. ELSA is used coupled to a full ice sheet model and creates individual layers of accumulation with fixed time stamps during the simulation, modeling 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 parametrization.
Antonio Juarez-Martinez, Javier Blasco, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas
EGUsphere, https://doi.org/10.5194/egusphere-2023-2863, https://doi.org/10.5194/egusphere-2023-2863, 2024
Short summary
Short summary
We present sea-level projections for Antarctica in the context of ISMIP6-2300 with several forcings, but extending the simulations until 2500, and showing that more than 3 metres 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.
Jan Swierczek-Jereczek, Marisa Montoya, Konstantin Latychev, Alexander Robinson, Jorge Alvarez-Solas, and Jerry Mitrovica
EGUsphere, https://doi.org/10.5194/egusphere-2023-2869, https://doi.org/10.5194/egusphere-2023-2869, 2023
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 runs 100,000 simulation years in only minutes. Thus, the evolution of ice sheets is modelled with better accuracy, which is critical for a good mitigation of climate change and, in particular, sea-level rise.
Daniel Moreno-Parada, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas
EGUsphere, https://doi.org/10.5194/egusphere-2023-2690, https://doi.org/10.5194/egusphere-2023-2690, 2023
Short summary
Short summary
We introduce Nix, an ice-sheet model designed for understanding how large masses of ice behave. Nix as a computer program that simulates the movement and temperature changes in ice sheets. Nix helps us study how ice sheets respond to changes in the atmosphere and ocean. We found that how fast ice melts under the shelves and how heat is exchanged, play a role in determining the future of ice sheets. Nix is a useful tool for learning more about how climate change affects polar ice sheets.
Javier Blasco, Ilaria Tabone, Daniel Moreno-Parada, Alexander Robinson, Jorge Alvarez-Solas, Frank Pattyn, and Marisa Montoya
Clim. Past Discuss., https://doi.org/10.5194/cp-2023-76, https://doi.org/10.5194/cp-2023-76, 2023
Revised manuscript under review for CP
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.
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.
Daniel Moreno, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas
The Cryosphere Discuss., https://doi.org/10.5194/tc-2022-97, https://doi.org/10.5194/tc-2022-97, 2022
Preprint under review for TC
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.
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.
Alexander Robinson, Jorge Alvarez-Solas, Marisa Montoya, Heiko Goelzer, Ralf Greve, and Catherine Ritz
Geosci. Model Dev., 13, 2805–2823, https://doi.org/10.5194/gmd-13-2805-2020, https://doi.org/10.5194/gmd-13-2805-2020, 2020
Short summary
Short summary
Here we describe Yelmo v1.0, an intuitive and state-of-the-art hybrid ice sheet model. The model design and physics are described, and benchmark simulations are provided to validate its performance. Yelmo is a versatile ice sheet model that can be applied to a wide variety of problems.
Jorge Alvarez-Solas, Marisa Montoya, and Alexander Robinson
Clim. Past Discuss., https://doi.org/10.5194/cp-2019-96, https://doi.org/10.5194/cp-2019-96, 2019
Publication in CP not foreseen
Short summary
Short summary
Modelling the past abrupt climate changes often resorts to the use of freshwater flux (FWF) in the North Atlantic as an effective method to cause reorganizations of the Atlantic Meridional Overturning Circulation. This procedure has allowed to reproduce the timing of the events. However, the required FWF is inconsistent with reconstructions. Conversely, using a forcing derived from the sea-level record results in a poor fit with the data, highlighting the need of exploring other mechanisms.
Ilaria Tabone, Alexander Robinson, Jorge Alvarez-Solas, and Marisa Montoya
The Cryosphere, 13, 1911–1923, https://doi.org/10.5194/tc-13-1911-2019, https://doi.org/10.5194/tc-13-1911-2019, 2019
Short summary
Short summary
Recent reconstructions show that the North East Greenland Ice Stream (NEGIS) retreated away from its present-day position by 20–40 km during MIS-3. Atmospheric and external forcings were proposed as potential causes of this retreat, but the role of the ocean was not considered. Here, using a 3-D ice-sheet model, we suggest that oceanic warming is sufficient to induce a retreat of the NEGIS margin of many tens of kilometres during MIS-3, helping to explain this conundrum.
Jorge Alvarez-Solas, Rubén Banderas, Alexander Robinson, and Marisa Montoya
Clim. Past, 15, 957–979, https://doi.org/10.5194/cp-15-957-2019, https://doi.org/10.5194/cp-15-957-2019, 2019
Short summary
Short summary
The last glacial period was marked by the existence of of abrupt climatic changes; it is generally accepted that the presence of ice sheets played an important role in their occurrence. While an important effort has been made to investigate the dynamics and evolution of the Laurentide ice sheet during this period, the Eurasian ice sheet (EIS) has not received much attention. Here we investigate the response of the EIS to millennial-scale climate variability using a hybrid 3-D ice-sheet model.
Ilaria Tabone, Alexander Robinson, Jorge Alvarez-Solas, and Marisa Montoya
Clim. Past, 15, 593–609, https://doi.org/10.5194/cp-15-593-2019, https://doi.org/10.5194/cp-15-593-2019, 2019
Short summary
Short summary
By using a 3-D hybrid ice-sheet–shelf model, we investigate the impact of millennial-scale oceanic variability on the Greenland Ice Sheet (GrIS) evolution during the last glacial period (LGP). We show that the GrIS may have strongly reacted to oceanic temperature fluctuations associated with Dansgaard–Oeschger cycles, contributing to sea-level variations of more than 1 m. Our results open the chance for a non-negligible role of the GrIS in millennial-scale oceanic reorganisations during the LGP.
Javier Blasco, Ilaria Tabone, Jorge Alvarez-Solas, Alexander Robinson, and Marisa Montoya
Clim. Past, 15, 121–133, https://doi.org/10.5194/cp-15-121-2019, https://doi.org/10.5194/cp-15-121-2019, 2019
Short summary
Short summary
The LGP is a period punctuated by the presence of several abrupt climate events and sea-level variations of up to 20 m at millennial timescales. The origin of those fluctuations is attributed to NH paleo ice sheets, but a contribution from the AIS cannot be excluded. Here, for the first time, we investigate the response of the AIS to millennial climate variability using an ice sheet–shelf model. We shows that the AIS produces substantial sea-level rises and grounding line migrations.
Rubén Banderas, Jorge Alvarez-Solas, Alexander Robinson, and Marisa Montoya
Geosci. Model Dev., 11, 2299–2314, https://doi.org/10.5194/gmd-11-2299-2018, https://doi.org/10.5194/gmd-11-2299-2018, 2018
Short summary
Short summary
Here we present a new approach to force ice-sheet models offline, which accounts for a more realistic treatment of millennial-scale climate variability as compared to the existing methods. Our results reveal that an incorrect representation of the characteristic pattern of millennial-scale climate variability within the climate forcing not only affects NH ice-volume variations at millennial timescales but has consequences for glacial–interglacial ice-volume changes too.
Ilaria Tabone, Javier Blasco, Alexander Robinson, Jorge Alvarez-Solas, and Marisa Montoya
Clim. Past, 14, 455–472, https://doi.org/10.5194/cp-14-455-2018, https://doi.org/10.5194/cp-14-455-2018, 2018
Short summary
Short summary
The response of the Greenland Ice Sheet (GrIS) to palaeo-oceanic changes on a glacial–interglacial timescale is studied from a modelling perspective. A 3-D hybrid ice-sheet–shelf model which includes a parameterization of the basal melting rate at the GrIS marine margins is used. The results show that the oceanic forcing plays a key role in the GrIS evolution, not only by controlling the ice retreat during the deglaciation but also by driving the ice expansion in glacial periods.
Jorge Alvarez-Solas, Rubén Banderas, Alexander Robinson, and Marisa Montoya
Clim. Past Discuss., https://doi.org/10.5194/cp-2017-143, https://doi.org/10.5194/cp-2017-143, 2017
Revised manuscript not accepted
Short summary
Short summary
The last glacial period was marked by the existence of of abrupt climatic changes. It is generally accepted that the presence of ice sheets played an
important role in their occurrence. While an important effort has been made to investigate the dynamics and evolution of the Laurentide Ice Sheet during this period, the Eurasian Ice Sheet (EIS) has not received much attention. Here we investigate the response of the EIS to millennial-scale climate variability. We use a hybrid 3D ice-sheet model.
Mario Krapp, Alexander Robinson, and Andrey Ganopolski
The Cryosphere, 11, 1519–1535, https://doi.org/10.5194/tc-11-1519-2017, https://doi.org/10.5194/tc-11-1519-2017, 2017
Short summary
Short summary
We present the snowpack model SEMIC. It calculates snow height, surface temperature, surface albedo, and the surface mass balance of snow- and ice-covered surfaces while using meteorological data as input. In this paper we describe how SEMIC works and how well it compares with snowpack data of a more sophisticated regional climate model applied to the Greenland ice sheet. Because of its simplicity and efficiency, SEMIC can be used as a coupling interface between atmospheric and ice sheet models.
A. Robinson and M. Perrette
Geosci. Model Dev., 8, 1877–1883, https://doi.org/10.5194/gmd-8-1877-2015, https://doi.org/10.5194/gmd-8-1877-2015, 2015
Short summary
Short summary
Here we present a concise interface to the NetCDF library designed to simplify reading and writing tasks of up to 6-D arrays in Fortran programs.
R. Calov, A. Robinson, M. Perrette, and A. Ganopolski
The Cryosphere, 9, 179–196, https://doi.org/10.5194/tc-9-179-2015, https://doi.org/10.5194/tc-9-179-2015, 2015
Short summary
Short summary
Ice discharge into the ocean from outlet glaciers is an important
component of mass loss of the Greenland ice sheet. Here, we present a
simple parameterization of ice discharge for coarse resolution ice
sheet models, suitable for large ensembles or long-term palaeo
simulations. This parameterization reproduces in a good approximation
the present-day ice discharge compared with estimates, and the
simulation of the present-day ice sheet elevation is considerably
improved.
A. Robinson and H. Goelzer
The Cryosphere, 8, 1419–1428, https://doi.org/10.5194/tc-8-1419-2014, https://doi.org/10.5194/tc-8-1419-2014, 2014
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 Discuss., https://doi.org/10.5194/cp-2024-19, https://doi.org/10.5194/cp-2024-19, 2024
Preprint under review for CP
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.
William H. Lipscomb, David Behar, and Monica Ainhorn Morrison
EGUsphere, https://doi.org/10.5194/egusphere-2024-534, https://doi.org/10.5194/egusphere-2024-534, 2024
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.
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.
Jean-Francois Lemieux, William Lipscomb, Anthony Craig, David A. Bailey, Elizabeth Hunke, Philippe Blain, Till A. S. Rasmussen, Mats Bentsen, Frederic Dupont, David Hebert, and Richard Allard
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2023-239, https://doi.org/10.5194/gmd-2023-239, 2024
Preprint under review for GMD
Short summary
Short summary
We present the latest version of the CICE model. CICE solves equations that describe the dynamics and the growth/melt of sea ice. To do so, the domain is divided in grid cells and the 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 a natural coupling with some oceanic and atmospheric models. It also allows the transport of ice in narrow channels.
Therese Rieckh, Andreas Born, Alexander Robinson, Robert Law, and Gerrit Gülle
EGUsphere, https://doi.org/10.5194/egusphere-2023-3127, https://doi.org/10.5194/egusphere-2023-3127, 2024
Short summary
Short summary
We present the open-source model ELSA, which simulates the internal age structure of large ice sheets. ELSA is used coupled to a full ice sheet model and creates individual layers of accumulation with fixed time stamps during the simulation, modeling 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 parametrization.
Antonio Juarez-Martinez, Javier Blasco, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas
EGUsphere, https://doi.org/10.5194/egusphere-2023-2863, https://doi.org/10.5194/egusphere-2023-2863, 2024
Short summary
Short summary
We present sea-level projections for Antarctica in the context of ISMIP6-2300 with several forcings, but extending the simulations until 2500, and showing that more than 3 metres 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.
Jan Swierczek-Jereczek, Marisa Montoya, Konstantin Latychev, Alexander Robinson, Jorge Alvarez-Solas, and Jerry Mitrovica
EGUsphere, https://doi.org/10.5194/egusphere-2023-2869, https://doi.org/10.5194/egusphere-2023-2869, 2023
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 runs 100,000 simulation years in only minutes. Thus, the evolution of ice sheets is modelled with better accuracy, which is critical for a good mitigation of climate change and, in particular, sea-level rise.
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.
Daniel Moreno-Parada, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas
EGUsphere, https://doi.org/10.5194/egusphere-2023-2690, https://doi.org/10.5194/egusphere-2023-2690, 2023
Short summary
Short summary
We introduce Nix, an ice-sheet model designed for understanding how large masses of ice behave. Nix as a computer program that simulates the movement and temperature changes in ice sheets. Nix helps us study how ice sheets respond to changes in the atmosphere and ocean. We found that how fast ice melts under the shelves and how heat is exchanged, play a role in determining the future of ice sheets. Nix is a useful tool for learning more about how climate change affects polar ice sheets.
Michele Petrini, Meike Scherrenberg, Laura Muntjewerf, Miren Vizcaino, Raymond Sellevold, Gunter Leguy, William Lipscomb, and Heiko Goelzer
The Cryosphere Discuss., https://doi.org/10.5194/tc-2023-154, https://doi.org/10.5194/tc-2023-154, 2023
Preprint under review for TC
Short summary
Short summary
In this study, we investigate with a numerical model the stability of the Greenland ice-sheet under prolonged sustained warming and ice melt. We show that there is a threshold beyond which the ice-sheet will lose more than 80 % of its mass over tens of thousand of years. The point of no return is reached when the ice-sheet disconnects from a region of high topography in western Greenland. This threshold is determined by the interaction of surface and solid-Earth processes.
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.
Javier Blasco, Ilaria Tabone, Daniel Moreno-Parada, Alexander Robinson, Jorge Alvarez-Solas, Frank Pattyn, and Marisa Montoya
Clim. Past Discuss., https://doi.org/10.5194/cp-2023-76, https://doi.org/10.5194/cp-2023-76, 2023
Revised manuscript under review for CP
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.
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.
David T. Bett, Alexander T. Bradley, C. Rosie Williams, Paul R. Holland, Robert J. Arthern, and Daniel N. Goldberg
The Cryosphere Discuss., https://doi.org/10.5194/tc-2023-77, https://doi.org/10.5194/tc-2023-77, 2023
Revised manuscript accepted for TC
Short summary
Short summary
A new specialised ice/ocean coupled model simulates the future ice sheet evolution in the Amundsen Sea sector of Antarctica. The model predicts substantial ocean-driven ice retreat over the 125-year simulations. The future of small ‘pinning points’ (islands of grounded ice) are an important control on this evolution. Ocean melting is crucial in ungrounding these pinning points, providing the mechanism by which future climate change may affect the sea level contribution from this sector.
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.
Daniel Moreno, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas
The Cryosphere Discuss., https://doi.org/10.5194/tc-2022-97, https://doi.org/10.5194/tc-2022-97, 2022
Preprint under review for TC
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.
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.
Sophie Nowicki, Heiko Goelzer, Hélène Seroussi, Anthony J. Payne, William H. Lipscomb, Ayako Abe-Ouchi, Cécile Agosta, Patrick Alexander, Xylar S. Asay-Davis, Alice Barthel, Thomas J. Bracegirdle, Richard Cullather, Denis Felikson, Xavier Fettweis, Jonathan M. Gregory, Tore Hattermann, Nicolas C. Jourdain, Peter Kuipers Munneke, Eric Larour, Christopher M. Little, Mathieu Morlighem, Isabel Nias, Andrew Shepherd, Erika Simon, Donald Slater, Robin S. Smith, Fiammetta Straneo, Luke D. Trusel, Michiel R. van den Broeke, and Roderik van de Wal
The Cryosphere, 14, 2331–2368, https://doi.org/10.5194/tc-14-2331-2020, https://doi.org/10.5194/tc-14-2331-2020, 2020
Short summary
Short summary
This paper describes the experimental protocol for ice sheet models taking part in the Ice Sheet Model Intercomparion Project for CMIP6 (ISMIP6) and presents an overview of the atmospheric and oceanic datasets to be used for the simulations. The ISMIP6 framework allows for exploring the uncertainty in 21st century sea level change from the Greenland and Antarctic ice sheets.
Stephen L. Cornford, Helene Seroussi, Xylar S. Asay-Davis, G. Hilmar Gudmundsson, Rob Arthern, Chris Borstad, Julia Christmann, Thiago Dias dos Santos, Johannes Feldmann, Daniel Goldberg, Matthew J. Hoffman, Angelika Humbert, Thomas Kleiner, Gunter Leguy, William H. Lipscomb, Nacho Merino, Gaël Durand, Mathieu Morlighem, David Pollard, Martin Rückamp, C. Rosie Williams, and Hongju Yu
The Cryosphere, 14, 2283–2301, https://doi.org/10.5194/tc-14-2283-2020, https://doi.org/10.5194/tc-14-2283-2020, 2020
Short summary
Short summary
We present the results of the third Marine Ice Sheet Intercomparison Project (MISMIP+). MISMIP+ is one in a series of exercises that test numerical models of ice sheet flow in simple situations. This particular exercise concentrates on the response of ice sheet models to the thinning of their floating ice shelves, which is of interest because numerical models are currently used to model the response to contemporary and near-future thinning in Antarctic ice shelves.
Alexander Robinson, Jorge Alvarez-Solas, Marisa Montoya, Heiko Goelzer, Ralf Greve, and Catherine Ritz
Geosci. Model Dev., 13, 2805–2823, https://doi.org/10.5194/gmd-13-2805-2020, https://doi.org/10.5194/gmd-13-2805-2020, 2020
Short summary
Short summary
Here we describe Yelmo v1.0, an intuitive and state-of-the-art hybrid ice sheet model. The model design and physics are described, and benchmark simulations are provided to validate its performance. Yelmo is a versatile ice sheet model that can be applied to a wide variety of problems.
Heiko Goelzer, Brice P. Y. Noël, Tamsin L. Edwards, Xavier Fettweis, Jonathan M. Gregory, William H. Lipscomb, Roderik S. W. van de Wal, and Michiel R. van den Broeke
The Cryosphere, 14, 1747–1762, https://doi.org/10.5194/tc-14-1747-2020, https://doi.org/10.5194/tc-14-1747-2020, 2020
Short summary
Short summary
Future sea-level change projections with process-based ice sheet models are typically driven with surface mass balance forcing derived from climate models. In this work we address the problems arising from a mismatch of the modelled ice sheet geometry with the one used by the climate model. The proposed remapping method reproduces the original forcing data closely when applied to the original geometry and produces a physically meaningful forcing when applied to different modelled geometries.
Anders Levermann, Ricarda Winkelmann, Torsten Albrecht, Heiko Goelzer, Nicholas R. Golledge, Ralf Greve, Philippe Huybrechts, Jim Jordan, Gunter Leguy, Daniel Martin, Mathieu Morlighem, Frank Pattyn, David Pollard, Aurelien Quiquet, Christian Rodehacke, Helene Seroussi, Johannes Sutter, Tong Zhang, Jonas Van Breedam, Reinhard Calov, Robert DeConto, Christophe Dumas, Julius Garbe, G. Hilmar Gudmundsson, Matthew J. Hoffman, Angelika Humbert, Thomas Kleiner, William H. Lipscomb, Malte Meinshausen, Esmond Ng, Sophie M. J. Nowicki, Mauro Perego, Stephen F. Price, Fuyuki Saito, Nicole-Jeanne Schlegel, Sainan Sun, and Roderik S. W. van de Wal
Earth Syst. Dynam., 11, 35–76, https://doi.org/10.5194/esd-11-35-2020, https://doi.org/10.5194/esd-11-35-2020, 2020
Short summary
Short summary
We provide an estimate of the future sea level contribution of Antarctica from basal ice shelf melting up to the year 2100. The full uncertainty range in the warming-related forcing of basal melt is estimated and applied to 16 state-of-the-art ice sheet models using a linear response theory approach. The sea level contribution we obtain is very likely below 61 cm under unmitigated climate change until 2100 (RCP8.5) and very likely below 40 cm if the Paris Climate Agreement is kept.
Raymond Sellevold, Leonardus van Kampenhout, Jan T. M. Lenaerts, Brice Noël, William H. Lipscomb, and Miren Vizcaino
The Cryosphere, 13, 3193–3208, https://doi.org/10.5194/tc-13-3193-2019, https://doi.org/10.5194/tc-13-3193-2019, 2019
Short summary
Short summary
We evaluate a downscaling method to calculate ice sheet surface mass balance with global climate models, despite their coarse resolution. We compare it with high-resolution climate modeling. Despite absence of fine-scale simulation of individual energy and mass contributors, the method provides realistic vertical SMB gradients that can be used in forcing of ice sheet models, e.g., for sea level projections. Also, the climate model simulation is improved with the method implemented interactively.
Jorge Alvarez-Solas, Marisa Montoya, and Alexander Robinson
Clim. Past Discuss., https://doi.org/10.5194/cp-2019-96, https://doi.org/10.5194/cp-2019-96, 2019
Publication in CP not foreseen
Short summary
Short summary
Modelling the past abrupt climate changes often resorts to the use of freshwater flux (FWF) in the North Atlantic as an effective method to cause reorganizations of the Atlantic Meridional Overturning Circulation. This procedure has allowed to reproduce the timing of the events. However, the required FWF is inconsistent with reconstructions. Conversely, using a forcing derived from the sea-level record results in a poor fit with the data, highlighting the need of exploring other mechanisms.
Ilaria Tabone, Alexander Robinson, Jorge Alvarez-Solas, and Marisa Montoya
The Cryosphere, 13, 1911–1923, https://doi.org/10.5194/tc-13-1911-2019, https://doi.org/10.5194/tc-13-1911-2019, 2019
Short summary
Short summary
Recent reconstructions show that the North East Greenland Ice Stream (NEGIS) retreated away from its present-day position by 20–40 km during MIS-3. Atmospheric and external forcings were proposed as potential causes of this retreat, but the role of the ocean was not considered. Here, using a 3-D ice-sheet model, we suggest that oceanic warming is sufficient to induce a retreat of the NEGIS margin of many tens of kilometres during MIS-3, helping to explain this conundrum.
Jorge Alvarez-Solas, Rubén Banderas, Alexander Robinson, and Marisa Montoya
Clim. Past, 15, 957–979, https://doi.org/10.5194/cp-15-957-2019, https://doi.org/10.5194/cp-15-957-2019, 2019
Short summary
Short summary
The last glacial period was marked by the existence of of abrupt climatic changes; it is generally accepted that the presence of ice sheets played an important role in their occurrence. While an important effort has been made to investigate the dynamics and evolution of the Laurentide ice sheet during this period, the Eurasian ice sheet (EIS) has not received much attention. Here we investigate the response of the EIS to millennial-scale climate variability using a hybrid 3-D ice-sheet model.
Hélène Seroussi, Sophie Nowicki, Erika Simon, Ayako Abe-Ouchi, Torsten Albrecht, Julien Brondex, Stephen Cornford, Christophe Dumas, Fabien Gillet-Chaulet, Heiko Goelzer, Nicholas R. Golledge, Jonathan M. Gregory, Ralf Greve, Matthew J. Hoffman, Angelika Humbert, Philippe Huybrechts, Thomas Kleiner, Eric Larour, Gunter Leguy, William H. Lipscomb, Daniel Lowry, Matthias Mengel, Mathieu Morlighem, Frank Pattyn, Anthony J. Payne, David Pollard, Stephen F. Price, Aurélien Quiquet, Thomas J. Reerink, Ronja Reese, Christian B. Rodehacke, Nicole-Jeanne Schlegel, Andrew Shepherd, Sainan Sun, Johannes Sutter, Jonas Van Breedam, Roderik S. W. van de Wal, Ricarda Winkelmann, and Tong Zhang
The Cryosphere, 13, 1441–1471, https://doi.org/10.5194/tc-13-1441-2019, https://doi.org/10.5194/tc-13-1441-2019, 2019
Short summary
Short summary
We compare a wide range of Antarctic ice sheet simulations with varying initialization techniques and model parameters to understand the role they play on the projected evolution of this ice sheet under simple scenarios. Results are improved compared to previous assessments and show that continued improvements in the representation of the floating ice around Antarctica are critical to reduce the uncertainty in the future ice sheet contribution to sea level rise.
Ilaria Tabone, Alexander Robinson, Jorge Alvarez-Solas, and Marisa Montoya
Clim. Past, 15, 593–609, https://doi.org/10.5194/cp-15-593-2019, https://doi.org/10.5194/cp-15-593-2019, 2019
Short summary
Short summary
By using a 3-D hybrid ice-sheet–shelf model, we investigate the impact of millennial-scale oceanic variability on the Greenland Ice Sheet (GrIS) evolution during the last glacial period (LGP). We show that the GrIS may have strongly reacted to oceanic temperature fluctuations associated with Dansgaard–Oeschger cycles, contributing to sea-level variations of more than 1 m. Our results open the chance for a non-negligible role of the GrIS in millennial-scale oceanic reorganisations during the LGP.
William H. Lipscomb, Stephen F. Price, Matthew J. Hoffman, Gunter R. Leguy, Andrew R. Bennett, Sarah L. Bradley, Katherine J. Evans, Jeremy G. Fyke, Joseph H. Kennedy, Mauro Perego, Douglas M. Ranken, William J. Sacks, Andrew G. Salinger, Lauren J. Vargo, and Patrick H. Worley
Geosci. Model Dev., 12, 387–424, https://doi.org/10.5194/gmd-12-387-2019, https://doi.org/10.5194/gmd-12-387-2019, 2019
Short summary
Short summary
This paper describes the Community Ice Sheet Model (CISM) version 2.1. CISM solves equations for ice flow, heat conduction, surface melting, and other processes such as basal sliding and iceberg calving. It can be used for ice-sheet-only simulations or as the ice sheet component of the Community Earth System Model. Model solutions have been verified for standard test problems. CISM can efficiently simulate the whole Greenland ice sheet, with results that are broadly consistent with observations.
Javier Blasco, Ilaria Tabone, Jorge Alvarez-Solas, Alexander Robinson, and Marisa Montoya
Clim. Past, 15, 121–133, https://doi.org/10.5194/cp-15-121-2019, https://doi.org/10.5194/cp-15-121-2019, 2019
Short summary
Short summary
The LGP is a period punctuated by the presence of several abrupt climate events and sea-level variations of up to 20 m at millennial timescales. The origin of those fluctuations is attributed to NH paleo ice sheets, but a contribution from the AIS cannot be excluded. Here, for the first time, we investigate the response of the AIS to millennial climate variability using an ice sheet–shelf model. We shows that the AIS produces substantial sea-level rises and grounding line migrations.
Matthew J. Hoffman, Mauro Perego, Stephen F. Price, William H. Lipscomb, Tong Zhang, Douglas Jacobsen, Irina Tezaur, Andrew G. Salinger, Raymond Tuminaro, and Luca Bertagna
Geosci. Model Dev., 11, 3747–3780, https://doi.org/10.5194/gmd-11-3747-2018, https://doi.org/10.5194/gmd-11-3747-2018, 2018
Short summary
Short summary
MPAS-Albany Land Ice (MALI) is a new variable-resolution land ice model that uses unstructured grids on a plane or sphere. MALI is built for Earth system modeling on high-performance computing platforms using existing software libraries. MALI simulates the evolution of ice thickness, velocity, and temperature, and it includes schemes for simulating iceberg calving and the flow of water beneath ice sheets and its effect on ice sliding. The model is demonstrated for the Antarctic ice sheet.
Rubén Banderas, Jorge Alvarez-Solas, Alexander Robinson, and Marisa Montoya
Geosci. Model Dev., 11, 2299–2314, https://doi.org/10.5194/gmd-11-2299-2018, https://doi.org/10.5194/gmd-11-2299-2018, 2018
Short summary
Short summary
Here we present a new approach to force ice-sheet models offline, which accounts for a more realistic treatment of millennial-scale climate variability as compared to the existing methods. Our results reveal that an incorrect representation of the characteristic pattern of millennial-scale climate variability within the climate forcing not only affects NH ice-volume variations at millennial timescales but has consequences for glacial–interglacial ice-volume changes too.
Heiko Goelzer, Sophie Nowicki, Tamsin Edwards, Matthew Beckley, Ayako Abe-Ouchi, Andy Aschwanden, Reinhard Calov, Olivier Gagliardini, Fabien Gillet-Chaulet, Nicholas R. Golledge, Jonathan Gregory, Ralf Greve, Angelika Humbert, Philippe Huybrechts, Joseph H. Kennedy, Eric Larour, William H. Lipscomb, Sébastien Le clec'h, Victoria Lee, Mathieu Morlighem, Frank Pattyn, Antony J. Payne, Christian Rodehacke, Martin Rückamp, Fuyuki Saito, Nicole Schlegel, Helene Seroussi, Andrew Shepherd, Sainan Sun, Roderik van de Wal, and Florian A. Ziemen
The Cryosphere, 12, 1433–1460, https://doi.org/10.5194/tc-12-1433-2018, https://doi.org/10.5194/tc-12-1433-2018, 2018
Short summary
Short summary
We have compared a wide spectrum of different initialisation techniques used in the ice sheet modelling community to define the modelled present-day Greenland ice sheet state as a starting point for physically based future-sea-level-change projections. Compared to earlier community-wide comparisons, we find better agreement across different models, which implies overall improvement of our understanding of what is needed to produce such initial states.
Ilaria Tabone, Javier Blasco, Alexander Robinson, Jorge Alvarez-Solas, and Marisa Montoya
Clim. Past, 14, 455–472, https://doi.org/10.5194/cp-14-455-2018, https://doi.org/10.5194/cp-14-455-2018, 2018
Short summary
Short summary
The response of the Greenland Ice Sheet (GrIS) to palaeo-oceanic changes on a glacial–interglacial timescale is studied from a modelling perspective. A 3-D hybrid ice-sheet–shelf model which includes a parameterization of the basal melting rate at the GrIS marine margins is used. The results show that the oceanic forcing plays a key role in the GrIS evolution, not only by controlling the ice retreat during the deglaciation but also by driving the ice expansion in glacial periods.
Jorge Alvarez-Solas, Rubén Banderas, Alexander Robinson, and Marisa Montoya
Clim. Past Discuss., https://doi.org/10.5194/cp-2017-143, https://doi.org/10.5194/cp-2017-143, 2017
Revised manuscript not accepted
Short summary
Short summary
The last glacial period was marked by the existence of of abrupt climatic changes. It is generally accepted that the presence of ice sheets played an
important role in their occurrence. While an important effort has been made to investigate the dynamics and evolution of the Laurentide Ice Sheet during this period, the Eurasian Ice Sheet (EIS) has not received much attention. Here we investigate the response of the EIS to millennial-scale climate variability. We use a hybrid 3D ice-sheet model.
Mario Krapp, Alexander Robinson, and Andrey Ganopolski
The Cryosphere, 11, 1519–1535, https://doi.org/10.5194/tc-11-1519-2017, https://doi.org/10.5194/tc-11-1519-2017, 2017
Short summary
Short summary
We present the snowpack model SEMIC. It calculates snow height, surface temperature, surface albedo, and the surface mass balance of snow- and ice-covered surfaces while using meteorological data as input. In this paper we describe how SEMIC works and how well it compares with snowpack data of a more sophisticated regional climate model applied to the Greenland ice sheet. Because of its simplicity and efficiency, SEMIC can be used as a coupling interface between atmospheric and ice sheet models.
Stephen F. Price, Matthew J. Hoffman, Jennifer A. Bonin, Ian M. Howat, Thomas Neumann, Jack Saba, Irina Tezaur, Jeffrey Guerber, Don P. Chambers, Katherine J. Evans, Joseph H. Kennedy, Jan Lenaerts, William H. Lipscomb, Mauro Perego, Andrew G. Salinger, Raymond S. Tuminaro, Michiel R. van den Broeke, and Sophie M. J. Nowicki
Geosci. Model Dev., 10, 255–270, https://doi.org/10.5194/gmd-10-255-2017, https://doi.org/10.5194/gmd-10-255-2017, 2017
Short summary
Short summary
We introduce the Cryospheric Model Comparison Tool (CmCt) and propose qualitative and quantitative metrics for evaluating ice sheet model simulations against observations. Greenland simulations using the Community Ice Sheet Model are compared to gravimetry and altimetry observations from 2003 to 2013. We show that the CmCt can be used to score simulations of increasing complexity relative to observations of dynamic change in Greenland over the past decade.
Daniel N. Goldberg, Sri Hari Krishna Narayanan, Laurent Hascoet, and Jean Utke
Geosci. Model Dev., 9, 1891–1904, https://doi.org/10.5194/gmd-9-1891-2016, https://doi.org/10.5194/gmd-9-1891-2016, 2016
Short summary
Short summary
Geophysical adjoint models are powerful tools, allowing sensitivity studies that are not possible otherwise, and enabling optimized fit of models to observing data sets. The complexity involved requires the use of algorithmic differentiation (AD) software, but AD adjoint calculation for ice models can be slow, with prohibitive memory requirements. In this paper, we present a method to improve the performance of ice model adjoint generation, in terms of timing, memory load, and accuracy.
D. N. Goldberg, P. Heimbach, I. Joughin, and B. Smith
The Cryosphere, 9, 2429–2446, https://doi.org/10.5194/tc-9-2429-2015, https://doi.org/10.5194/tc-9-2429-2015, 2015
Short summary
Short summary
We calibrate a time-dependent ice model through optimal fit to transient observations of surface elevation and velocity, a novel procedure in glaciology and in particular for an ice stream with a dynamic grounding line. We show this procedure gives a level of confidence in model projections that cannot be achieved through more commonly used glaciological data assimilation methods. We show that Smith Glacier is in a state of retreat regardless of climatic forcing for the next several decades.
A. Robinson and M. Perrette
Geosci. Model Dev., 8, 1877–1883, https://doi.org/10.5194/gmd-8-1877-2015, https://doi.org/10.5194/gmd-8-1877-2015, 2015
Short summary
Short summary
Here we present a concise interface to the NetCDF library designed to simplify reading and writing tasks of up to 6-D arrays in Fortran programs.
R. Calov, A. Robinson, M. Perrette, and A. Ganopolski
The Cryosphere, 9, 179–196, https://doi.org/10.5194/tc-9-179-2015, https://doi.org/10.5194/tc-9-179-2015, 2015
Short summary
Short summary
Ice discharge into the ocean from outlet glaciers is an important
component of mass loss of the Greenland ice sheet. Here, we present a
simple parameterization of ice discharge for coarse resolution ice
sheet models, suitable for large ensembles or long-term palaeo
simulations. This parameterization reproduces in a good approximation
the present-day ice discharge compared with estimates, and the
simulation of the present-day ice sheet elevation is considerably
improved.
A. Robinson and H. Goelzer
The Cryosphere, 8, 1419–1428, https://doi.org/10.5194/tc-8-1419-2014, https://doi.org/10.5194/tc-8-1419-2014, 2014
G. R. Leguy, X. S. Asay-Davis, and W. H. Lipscomb
The Cryosphere, 8, 1239–1259, https://doi.org/10.5194/tc-8-1239-2014, https://doi.org/10.5194/tc-8-1239-2014, 2014
J. G. Fyke, W. J. Sacks, and W. H. Lipscomb
Geosci. Model Dev., 7, 1183–1195, https://doi.org/10.5194/gmd-7-1183-2014, https://doi.org/10.5194/gmd-7-1183-2014, 2014
D. N. Goldberg and P. Heimbach
The Cryosphere, 7, 1659–1678, https://doi.org/10.5194/tc-7-1659-2013, https://doi.org/10.5194/tc-7-1659-2013, 2013
S. H. Mernild, W. H. Lipscomb, D. B. Bahr, V. Radić, and M. Zemp
The Cryosphere, 7, 1565–1577, https://doi.org/10.5194/tc-7-1565-2013, https://doi.org/10.5194/tc-7-1565-2013, 2013
Related subject area
Discipline: Ice sheets | Subject: Numerical Modelling
Modeling the timing of Patagonian Ice Sheet retreat in the Chilean Lake District from 22–10 ka
Using specularity content to evaluate eight geothermal heat flow maps of Totten Glacier
Surging of a Hudson Strait-scale ice stream: subglacial hydrology matters but the process details mostly do not
Regularization and L-curves in ice sheet inverse models: a case study in the Filchner–Ronne catchment
Quantifying the uncertainty in the Eurasian ice-sheet geometry at the Penultimate Glacial Maximum (Marine Isotope Stage 6)
Biases in ice sheet models from missing noise-induced drift
The stability of present-day Antarctic grounding lines – Part 2: Onset of irreversible retreat of Amundsen Sea glaciers under current climate on centennial timescales cannot be excluded
The stability of present-day Antarctic grounding lines – Part 1: No indication of marine ice sheet instability in the current geometry
Geothermal heat flux is the dominant source of uncertainty in englacial-temperature-based dating of ice rise formation
Improving interpretation of sea-level projections through a machine-learning-based local explanation approach
Subglacial hydrology modulates basal sliding response of the Antarctic ice sheet to climate forcing
The predictive power of ice sheet models and the regional sensitivity of ice loss to basal sliding parameterisations: a case study of Pine Island and Thwaites glaciers, West Antarctica
Simulations of firn processes over the Greenland and Antarctic ice sheets: 1980–2021
Evaluation of six geothermal heat flux maps for the Antarctic Lambert–Amery glacial system
Impact of runoff temporal distribution on ice dynamics
Can changes in deformation regimes be inferred from crystallographic preferred orientations in polar ice?
Stabilizing effect of mélange buttressing on the marine ice-cliff instability of the West Antarctic Ice Sheet
Effective coefficient of diffusion and permeability of firn at Dome C and Lock In, Antarctica, and of various snow types – estimates over the 100–850 kg m−3 density range
The instantaneous impact of calving and thinning on the Larsen C Ice Shelf
Derivation of bedrock topography measurement requirements for the reduction of uncertainty in ice-sheet model projections of Thwaites Glacier
A new vertically integrated MOno-Layer Higher-Order (MOLHO) ice flow model
On the contribution of grain boundary sliding type creep to firn densification – an assessment using an optimization approach
Marine ice sheet experiments with the Community Ice Sheet Model
The transferability of adjoint inversion products between different ice flow models
Inferring the basal sliding coefficient field for the Stokes ice sheet model under rheological uncertainty
The tipping points and early warning indicators for Pine Island Glacier, West Antarctica
Sensitivity of ice sheet surface velocity and elevation to variations in basal friction and topography in the full Stokes and shallow-shelf approximation frameworks using adjoint equations
Quantifying the effect of ocean bed properties on ice sheet geometry over 40 000 years with a full-Stokes model
Bayesian calibration of firn densification models
A kinematic formalism for tracking ice–ocean mass exchange on the Earth's surface and estimating sea-level change
Results of the third Marine Ice Sheet Model Intercomparison Project (MISMIP+)
Ocean-forced evolution of the Amundsen Sea catchment, West Antarctica, by 2100
Parameter sensitivity analysis of dynamic ice sheet models – numerical computations
Simulated retreat of Jakobshavn Isbræ during the 21st century
Development of physically based liquid water schemes for Greenland firn-densification models
Regional grid refinement in an Earth system model: impacts on the simulated Greenland surface mass balance
initMIP-Antarctica: an ice sheet model initialization experiment of ISMIP6
Modeling the response of northwest Greenland to enhanced ocean thermal forcing and subglacial discharge
Assessment of the Greenland ice sheet–atmosphere feedbacks for the next century with a regional atmospheric model coupled to an ice sheet model
Sensitivity of centennial mass loss projections of the Amundsen basin to the friction law
Retreat of Thwaites Glacier, West Antarctica, over the next 100 years using various ice flow models, ice shelf melt scenarios and basal friction laws
Comparison of four calving laws to model Greenland outlet glaciers
Neutral equilibrium and forcing feedbacks in marine ice sheet modelling
Representation of basal melting at the grounding line in ice flow models
Stopping the flood: could we use targeted geoengineering to mitigate sea level rise?
Basal friction of Fleming Glacier, Antarctica – Part 1: Sensitivity of inversion to temperature and bedrock uncertainty
Basal friction of Fleming Glacier, Antarctica – Part 2: Evolution from 2008 to 2015
Simulated dynamic regrounding during marine ice sheet retreat
Dynamic response of Antarctic Peninsula Ice Sheet to potential collapse of Larsen C and George VI ice shelves
Simulated retreat of Jakobshavn Isbræ since the Little Ice Age controlled by geometry
Joshua Cuzzone, Matias Romero, and Shaun A. Marcott
The Cryosphere, 18, 1381–1398, https://doi.org/10.5194/tc-18-1381-2024, https://doi.org/10.5194/tc-18-1381-2024, 2024
Short summary
Short summary
We simulate the retreat history of the Patagonian Ice Sheet (PIS) across the Chilean Lake District from 22–10 ka. These results improve our understanding of the response of the PIS to deglacial warming and the patterns of deglacial ice margin retreat where gaps in the geologic record still exist, and they indicate that changes in large-scale precipitation during the last deglaciation played an important role in modulating the response of ice margin change across the PIS to deglacial warming.
Yan Huang, Liyun Zhao, Michael Wolovick, Yiliang Ma, and John C. Moore
The Cryosphere, 18, 103–119, https://doi.org/10.5194/tc-18-103-2024, https://doi.org/10.5194/tc-18-103-2024, 2024
Short summary
Short summary
Geothermal heat flux (GHF) is an important factor affecting the basal thermal environment of an ice sheet and crucial for its dynamics. But it is poorly defined for the Antarctic ice sheet. We simulate the basal temperature and basal melting rate with eight different GHF datasets. We use specularity content as a two-sided constraint to discriminate between local wet or dry basal conditions. Two medium-magnitude GHF distribution maps rank well, showing that most of the inland bed area is frozen.
Matthew Drew and Lev Tarasov
The Cryosphere, 17, 5391–5415, https://doi.org/10.5194/tc-17-5391-2023, https://doi.org/10.5194/tc-17-5391-2023, 2023
Short summary
Short summary
The interaction of fast-flowing regions of continental ice sheets with their beds governs how quickly they slide and therefore flow. The coupling of fast ice to its bed is controlled by the pressure of meltwater at its base. It is currently poorly understood how the physical details of these hydrologic systems affect ice speedup. Using numerical models we find, surprisingly, that they largely do not, except for the duration of the surge. This suggests that cheap models are sufficient.
Michael Wolovick, Angelika Humbert, Thomas Kleiner, and Martin Rückamp
The Cryosphere, 17, 5027–5060, https://doi.org/10.5194/tc-17-5027-2023, https://doi.org/10.5194/tc-17-5027-2023, 2023
Short summary
Short summary
The friction underneath ice sheets can be inferred from observed velocity at the top, but this inference requires smoothing. The selection of smoothing has been highly variable in the literature. Here we show how to rigorously select the best smoothing, and we show that the inferred friction converges towards the best knowable field as model resolution improves. We use this to learn about the best description of basal friction and to formulate recommended best practices for other modelers.
Oliver G. Pollard, Natasha L. M. Barlow, Lauren J. Gregoire, Natalya Gomez, Víctor Cartelle, Jeremy C. Ely, and Lachlan C. Astfalck
The Cryosphere, 17, 4751–4777, https://doi.org/10.5194/tc-17-4751-2023, https://doi.org/10.5194/tc-17-4751-2023, 2023
Short summary
Short summary
We use advanced statistical techniques and a simple ice-sheet model to produce an ensemble of plausible 3D shapes of the ice sheet that once stretched across northern Europe during the previous glacial maximum (140,000 years ago). This new reconstruction, equivalent in volume to 48 ± 8 m of global mean sea-level rise, will improve the interpretation of high sea levels recorded from the Last Interglacial period (120 000 years ago) that provide a useful perspective on the future.
Alexander Robel, Vincent Verjans, and Aminat Ambelorun
EGUsphere, https://doi.org/10.5194/egusphere-2023-2546, https://doi.org/10.5194/egusphere-2023-2546, 2023
Short summary
Short summary
The average size of many glaciers and ice sheets changes when noise is added to the system. The reasons for this "drift" in glacier state is intrinsic to the dynamics of how ice flows and the bumpiness of the Earth's surface. We argue that not including noise in projections of ice sheet evolution over coming decades and centuries is a pervasive source of bias in these computer models, and so realistic variability in glacier and climate processes must be included in models.
Ronja Reese, Julius Garbe, Emily A. Hill, Benoît Urruty, Kaitlin A. Naughten, Olivier Gagliardini, Gaël Durand, Fabien Gillet-Chaulet, G. Hilmar Gudmundsson, David Chandler, Petra M. Langebroek, and Ricarda Winkelmann
The Cryosphere, 17, 3761–3783, https://doi.org/10.5194/tc-17-3761-2023, https://doi.org/10.5194/tc-17-3761-2023, 2023
Short summary
Short summary
We use an ice sheet model to test where current climate conditions in Antarctica might lead. We find that present-day ocean and atmosphere conditions might commit an irreversible collapse of parts of West Antarctica which evolves over centuries to millennia. Importantly, this collapse is not irreversible yet.
Emily A. Hill, Benoît Urruty, Ronja Reese, Julius Garbe, Olivier Gagliardini, Gaël Durand, Fabien Gillet-Chaulet, G. Hilmar Gudmundsson, Ricarda Winkelmann, Mondher Chekki, David Chandler, and Petra M. Langebroek
The Cryosphere, 17, 3739–3759, https://doi.org/10.5194/tc-17-3739-2023, https://doi.org/10.5194/tc-17-3739-2023, 2023
Short summary
Short summary
The grounding lines of the Antarctic Ice Sheet could enter phases of irreversible retreat or advance. We use three ice sheet models to show that the present-day locations of Antarctic grounding lines are reversible with respect to a small perturbation away from their current position. This indicates that present-day retreat of the grounding lines is not yet irreversible or self-enhancing.
Aleksandr Montelli and Jonathan Kingslake
The Cryosphere, 17, 195–210, https://doi.org/10.5194/tc-17-195-2023, https://doi.org/10.5194/tc-17-195-2023, 2023
Short summary
Short summary
Thermal modelling and Bayesian inversion techniques are used to evaluate the uncertainties inherent in inferences of ice-sheet evolution from borehole temperature measurements. We show that the same temperature profiles may result from a range of parameters, of which geothermal heat flux through underlying bedrock plays a key role. Careful model parameterisation and evaluation of heat flux are essential for inferring past ice-sheet evolution from englacial borehole thermometry.
Jeremy Rohmer, Remi Thieblemont, Goneri Le Cozannet, Heiko Goelzer, and Gael Durand
The Cryosphere, 16, 4637–4657, https://doi.org/10.5194/tc-16-4637-2022, https://doi.org/10.5194/tc-16-4637-2022, 2022
Short summary
Short summary
To improve the interpretability of process-based projections of the sea-level contribution from land ice components, we apply the machine-learning-based
SHapley Additive exPlanationsapproach to a subset of a multi-model ensemble study for the Greenland ice sheet. This allows us to quantify the influence of particular modelling decisions (related to numerical implementation, initial conditions, or parametrisation of ice-sheet processes) directly in terms of sea-level change contribution.
Elise Kazmierczak, Sainan Sun, Violaine Coulon, and Frank Pattyn
The Cryosphere, 16, 4537–4552, https://doi.org/10.5194/tc-16-4537-2022, https://doi.org/10.5194/tc-16-4537-2022, 2022
Short summary
Short summary
The water at the interface between ice sheets and underlying bedrock leads to lubrication between the ice and the bed. Due to a lack of direct observations, subglacial conditions beneath the Antarctic ice sheet are poorly understood. Here, we compare different approaches in which the subglacial water could influence sliding on the underlying bedrock and suggest that it modulates the Antarctic ice sheet response and increases uncertainties, especially in the context of global warming.
Jowan M. Barnes and G. Hilmar Gudmundsson
The Cryosphere, 16, 4291–4304, https://doi.org/10.5194/tc-16-4291-2022, https://doi.org/10.5194/tc-16-4291-2022, 2022
Short summary
Short summary
Models must represent how glaciers slide along the bed, but there are many ways to do so. In this paper, several sliding laws are tested and found to affect different regions of the Antarctic Ice Sheet in different ways and at different speeds. However, the variability in ice volume loss due to sliding-law choices is low compared to other factors, so limited empirical knowledge of sliding does not prevent us from making predictions of how an ice sheet will evolve.
Brooke Medley, Thomas A. Neumann, H. Jay Zwally, Benjamin E. Smith, and C. Max Stevens
The Cryosphere, 16, 3971–4011, https://doi.org/10.5194/tc-16-3971-2022, https://doi.org/10.5194/tc-16-3971-2022, 2022
Short summary
Short summary
Satellite altimeters measure the height or volume change over Earth's ice sheets, but in order to understand how that change translates into ice mass, we must account for various processes at the surface. Specifically, snowfall events generate large, transient increases in surface height, yet snow fall has a relatively low density, which means much of that height change is composed of air. This air signal must be removed from the observed height changes before we can assess ice mass change.
Haoran Kang, Liyun Zhao, Michael Wolovick, and John C. Moore
The Cryosphere, 16, 3619–3633, https://doi.org/10.5194/tc-16-3619-2022, https://doi.org/10.5194/tc-16-3619-2022, 2022
Short summary
Short summary
Basal thermal conditions are important to ice dynamics and sensitive to geothermal heat flux (GHF). We estimate basal thermal conditions of the Lambert–Amery Glacier system with six GHF maps. Recent GHFs inverted from aerial geomagnetic observations produce a larger warm-based area and match the observed subglacial lakes better than the other GHFs. The modelled basal melt rate is 10 to hundreds of millimetres per year in fast-flowing glaciers feeding the Amery Ice Shelf and smaller inland.
Basile de Fleurian, Richard Davy, and Petra M. Langebroek
The Cryosphere, 16, 2265–2283, https://doi.org/10.5194/tc-16-2265-2022, https://doi.org/10.5194/tc-16-2265-2022, 2022
Short summary
Short summary
As temperature increases, more snow and ice melt at the surface of ice sheets. Here we use an ice dynamics and subglacial hydrology model with simplified geometry and climate forcing to study the impact of variations in meltwater on ice dynamics. We focus on the variations in length and intensity of the melt season. Our results show that a longer melt season leads to faster glaciers, but a more intense melt season reduces glaciers' seasonal velocities, albeit leading to higher peak velocities.
Maria-Gema Llorens, Albert Griera, Paul D. Bons, Ilka Weikusat, David J. Prior, Enrique Gomez-Rivas, Tamara de Riese, Ivone Jimenez-Munt, Daniel García-Castellanos, and Ricardo A. Lebensohn
The Cryosphere, 16, 2009–2024, https://doi.org/10.5194/tc-16-2009-2022, https://doi.org/10.5194/tc-16-2009-2022, 2022
Short summary
Short summary
Polar ice is formed by ice crystals, which form fabrics that are utilised to interpret how ice sheets flow. It is unclear whether fabrics result from the current flow regime or if they are inherited. To understand the extent to which ice crystals can be reoriented when ice flow conditions change, we simulate and evaluate multi-stage ice flow scenarios according to natural cases. We find that second deformation regimes normally overprint inherited fabrics, with a range of transitional fabrics.
Tanja Schlemm, Johannes Feldmann, Ricarda Winkelmann, and Anders Levermann
The Cryosphere, 16, 1979–1996, https://doi.org/10.5194/tc-16-1979-2022, https://doi.org/10.5194/tc-16-1979-2022, 2022
Short summary
Short summary
Marine cliff instability, if it exists, could dominate Antarctica's contribution to future sea-level rise. It is likely to speed up with ice thickness and thus would accelerate in most parts of Antarctica. Here, we investigate a possible mechanism that might stop cliff instability through cloaking by ice mélange. It is only a first step, but it shows that embayment geometry is, in principle, able to stop marine cliff instability in most parts of West Antarctica.
Neige Calonne, Alexis Burr, Armelle Philip, Frédéric Flin, and Christian Geindreau
The Cryosphere, 16, 967–980, https://doi.org/10.5194/tc-16-967-2022, https://doi.org/10.5194/tc-16-967-2022, 2022
Short summary
Short summary
Modeling gas transport in ice sheets from surface to close-off is key to interpreting climate archives. Estimates of the diffusion coefficient and permeability of snow and firn are required but remain a large source of uncertainty. We present a new dataset of diffusion coefficients and permeability from 20 to 120 m depth at two Antarctic sites. We suggest predictive formulas to estimate both properties over the entire 100–850 kg m3 density range, i.e., anywhere within the ice sheet column.
Tom Mitcham, G. Hilmar Gudmundsson, and Jonathan L. Bamber
The Cryosphere, 16, 883–901, https://doi.org/10.5194/tc-16-883-2022, https://doi.org/10.5194/tc-16-883-2022, 2022
Short summary
Short summary
We modelled the response of the Larsen C Ice Shelf (LCIS) and its tributary glaciers to the calving of the A68 iceberg and validated our results with observations. We found that the impact was limited, confirming that mostly passive ice was calved. Through further calving experiments we quantified the total buttressing provided by the LCIS and found that over 80 % of the buttressing capacity is generated in the first 5 km of the ice shelf downstream of the grounding line.
Blake A. Castleman, Nicole-Jeanne Schlegel, Lambert Caron, Eric Larour, and Ala Khazendar
The Cryosphere, 16, 761–778, https://doi.org/10.5194/tc-16-761-2022, https://doi.org/10.5194/tc-16-761-2022, 2022
Short summary
Short summary
In the described study, we derive an uncertainty range for global mean sea level rise (SLR) contribution from Thwaites Glacier in a 200-year period under an extreme ocean warming scenario. We derive the spatial and vertical resolutions needed for bedrock data acquisition missions in order to limit global mean SLR contribution from Thwaites Glacier to ±2 cm in a 200-year period. We conduct sensitivity experiments in order to present the locations of critical regions in need of accurate mapping.
Thiago Dias dos Santos, Mathieu Morlighem, and Douglas Brinkerhoff
The Cryosphere, 16, 179–195, https://doi.org/10.5194/tc-16-179-2022, https://doi.org/10.5194/tc-16-179-2022, 2022
Short summary
Short summary
Projecting the future evolution of Greenland and Antarctica and their potential contribution to sea level rise often relies on computer simulations carried out by numerical ice sheet models. Here we present a new vertically integrated ice sheet model and assess its performance using different benchmarks. The new model shows results comparable to a three-dimensional model at relatively lower computational cost, suggesting that it is an excellent alternative for long-term simulations.
Timm Schultz, Ralf Müller, Dietmar Gross, and Angelika Humbert
The Cryosphere, 16, 143–158, https://doi.org/10.5194/tc-16-143-2022, https://doi.org/10.5194/tc-16-143-2022, 2022
Short summary
Short summary
Firn is the interstage product between snow and ice. Simulations describing the process of firn densification are used in the context of estimating mass changes of the ice sheets and past climate reconstructions. The first stage of firn densification takes place in the upper few meters of the firn column. We investigate how well a material law describing the process of grain boundary sliding works for the numerical simulation of firn densification in this stage.
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.
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.
Olalekan Babaniyi, Ruanui Nicholson, Umberto Villa, and Noémi Petra
The Cryosphere, 15, 1731–1750, https://doi.org/10.5194/tc-15-1731-2021, https://doi.org/10.5194/tc-15-1731-2021, 2021
Short summary
Short summary
We consider the problem of inferring unknown parameter fields under additional uncertainty for an ice sheet model from synthetic surface ice flow velocity measurements. Our results indicate that accounting for model uncertainty stemming from the presence of nuisance parameters is crucial. Namely our findings suggest that using nominal values for these parameters, as is often done in practice, without taking into account the resulting modeling error can lead to overconfident and biased results.
Sebastian H. R. Rosier, Ronja Reese, Jonathan F. Donges, Jan De Rydt, G. Hilmar Gudmundsson, and Ricarda Winkelmann
The Cryosphere, 15, 1501–1516, https://doi.org/10.5194/tc-15-1501-2021, https://doi.org/10.5194/tc-15-1501-2021, 2021
Short summary
Short summary
Pine Island Glacier has contributed more to sea-level rise over the past decades than any other glacier in Antarctica. Ice-flow modelling studies have shown that it can undergo periods of rapid mass loss, but no study has shown that these future changes could cross a tipping point and therefore be effectively irreversible. Here, we assess the stability of Pine Island Glacier, quantifying the changes in ocean temperatures required to cross future tipping points using statistical methods.
Gong Cheng, Nina Kirchner, and Per Lötstedt
The Cryosphere, 15, 715–742, https://doi.org/10.5194/tc-15-715-2021, https://doi.org/10.5194/tc-15-715-2021, 2021
Short summary
Short summary
We present an inverse modeling approach to improve the understanding of spatiotemporally variable processes at the inaccessible base of an ice sheet by determining the sensitivity of direct surface observations to perturbations of basal conditions. Time dependency is proved to be important in these types of problems. The effect of perturbations is analyzed based on analytical and numerical solutions.
Clemens Schannwell, Reinhard Drews, Todd A. Ehlers, Olaf Eisen, Christoph Mayer, Mika Malinen, Emma C. Smith, and Hannes Eisermann
The Cryosphere, 14, 3917–3934, https://doi.org/10.5194/tc-14-3917-2020, https://doi.org/10.5194/tc-14-3917-2020, 2020
Short summary
Short summary
To reduce uncertainties associated with sea level rise projections, an accurate representation of ice flow is paramount. Most ice sheet models rely on simplified versions of the underlying ice flow equations. Due to the high computational costs, ice sheet models based on the complete ice flow equations have been restricted to < 1000 years. Here, we present a new model setup that extends the applicability of such models by an order of magnitude, permitting simulations of 40 000 years.
Vincent Verjans, Amber A. Leeson, Christopher Nemeth, C. Max Stevens, Peter Kuipers Munneke, Brice Noël, and Jan Melchior van Wessem
The Cryosphere, 14, 3017–3032, https://doi.org/10.5194/tc-14-3017-2020, https://doi.org/10.5194/tc-14-3017-2020, 2020
Short summary
Short summary
Ice sheets are covered by a firn layer, which is the transition stage between fresh snow and ice. Accurate modelling of firn density properties is important in many glaciological aspects. Current models show disagreements, are mostly calibrated to match specific observations of firn density and lack thorough uncertainty analysis. We use a novel calibration method for firn models based on a Bayesian statistical framework, which results in improved model accuracy and in uncertainty evaluation.
Surendra Adhikari, Erik R. Ivins, Eric Larour, Lambert Caron, and Helene Seroussi
The Cryosphere, 14, 2819–2833, https://doi.org/10.5194/tc-14-2819-2020, https://doi.org/10.5194/tc-14-2819-2020, 2020
Short summary
Short summary
The mathematical formalism presented in this paper aims at simplifying computational strategies for tracking ice–ocean mass exchange in the Earth system. To this end, we define a set of generic, and quite simple, descriptions of evolving land, ocean and ice interfaces and present a unified method to compute the sea-level contribution of evolving ice sheets. The formalism can be applied to arbitrary geometries and at all timescales.
Stephen L. Cornford, Helene Seroussi, Xylar S. Asay-Davis, G. Hilmar Gudmundsson, Rob Arthern, Chris Borstad, Julia Christmann, Thiago Dias dos Santos, Johannes Feldmann, Daniel Goldberg, Matthew J. Hoffman, Angelika Humbert, Thomas Kleiner, Gunter Leguy, William H. Lipscomb, Nacho Merino, Gaël Durand, Mathieu Morlighem, David Pollard, Martin Rückamp, C. Rosie Williams, and Hongju Yu
The Cryosphere, 14, 2283–2301, https://doi.org/10.5194/tc-14-2283-2020, https://doi.org/10.5194/tc-14-2283-2020, 2020
Short summary
Short summary
We present the results of the third Marine Ice Sheet Intercomparison Project (MISMIP+). MISMIP+ is one in a series of exercises that test numerical models of ice sheet flow in simple situations. This particular exercise concentrates on the response of ice sheet models to the thinning of their floating ice shelves, which is of interest because numerical models are currently used to model the response to contemporary and near-future thinning in Antarctic ice shelves.
Alanna V. Alevropoulos-Borrill, Isabel J. Nias, Antony J. Payne, Nicholas R. Golledge, and Rory J. Bingham
The Cryosphere, 14, 1245–1258, https://doi.org/10.5194/tc-14-1245-2020, https://doi.org/10.5194/tc-14-1245-2020, 2020
Gong Cheng and Per Lötstedt
The Cryosphere, 14, 673–691, https://doi.org/10.5194/tc-14-673-2020, https://doi.org/10.5194/tc-14-673-2020, 2020
Short summary
Short summary
We present a time-dependent inverse method for ice sheet modeling. By investigating the sensitivity of the observations of the velocity and the height at the surface to the basal conditions of the ice, we show that if the basal parameters are time dependent, then time cannot be ignored in the inversion. By looking at the numerical features, we conclude that adding the height information of an ice sheet in the velocity inversion procedure could improve the robustness of the inference.
Xiaoran Guo, Liyun Zhao, Rupert M. Gladstone, Sainan Sun, and John C. Moore
The Cryosphere, 13, 3139–3153, https://doi.org/10.5194/tc-13-3139-2019, https://doi.org/10.5194/tc-13-3139-2019, 2019
Vincent Verjans, Amber A. Leeson, C. Max Stevens, Michael MacFerrin, Brice Noël, and Michiel R. van den Broeke
The Cryosphere, 13, 1819–1842, https://doi.org/10.5194/tc-13-1819-2019, https://doi.org/10.5194/tc-13-1819-2019, 2019
Short summary
Short summary
Firn models rely on empirical approaches for representing the percolation and refreezing of meltwater through the firn column. We develop liquid water schemes of different levels of complexity for firn models and compare their performances with respect to observations of density profiles from Greenland. Our results demonstrate that physically advanced water schemes do not lead to better agreement with density observations. Uncertainties in other processes contribute more to model discrepancy.
Leonardus van Kampenhout, Alan M. Rhoades, Adam R. Herrington, Colin M. Zarzycki, Jan T. M. Lenaerts, William J. Sacks, and Michiel R. van den Broeke
The Cryosphere, 13, 1547–1564, https://doi.org/10.5194/tc-13-1547-2019, https://doi.org/10.5194/tc-13-1547-2019, 2019
Short summary
Short summary
A new tool is evaluated in which the climate and surface mass balance (SMB) of the Greenland ice sheet are resolved at 55 and 28 km resolution, while the rest of the globe is modelled at ~110 km. The local refinement of resolution leads to improved accumulation (SMB > 0) compared to observations; however ablation (SMB < 0) is deteriorated in some regions. This is attributed to changes in cloud cover and a reduced effectiveness of a model-specific vertical downscaling technique.
Hélène Seroussi, Sophie Nowicki, Erika Simon, Ayako Abe-Ouchi, Torsten Albrecht, Julien Brondex, Stephen Cornford, Christophe Dumas, Fabien Gillet-Chaulet, Heiko Goelzer, Nicholas R. Golledge, Jonathan M. Gregory, Ralf Greve, Matthew J. Hoffman, Angelika Humbert, Philippe Huybrechts, Thomas Kleiner, Eric Larour, Gunter Leguy, William H. Lipscomb, Daniel Lowry, Matthias Mengel, Mathieu Morlighem, Frank Pattyn, Anthony J. Payne, David Pollard, Stephen F. Price, Aurélien Quiquet, Thomas J. Reerink, Ronja Reese, Christian B. Rodehacke, Nicole-Jeanne Schlegel, Andrew Shepherd, Sainan Sun, Johannes Sutter, Jonas Van Breedam, Roderik S. W. van de Wal, Ricarda Winkelmann, and Tong Zhang
The Cryosphere, 13, 1441–1471, https://doi.org/10.5194/tc-13-1441-2019, https://doi.org/10.5194/tc-13-1441-2019, 2019
Short summary
Short summary
We compare a wide range of Antarctic ice sheet simulations with varying initialization techniques and model parameters to understand the role they play on the projected evolution of this ice sheet under simple scenarios. Results are improved compared to previous assessments and show that continued improvements in the representation of the floating ice around Antarctica are critical to reduce the uncertainty in the future ice sheet contribution to sea level rise.
Mathieu Morlighem, Michael Wood, Hélène Seroussi, Youngmin Choi, and Eric Rignot
The Cryosphere, 13, 723–734, https://doi.org/10.5194/tc-13-723-2019, https://doi.org/10.5194/tc-13-723-2019, 2019
Short summary
Short summary
Many glaciers along the coast of Greenland have been retreating. It has been suggested that this retreat is triggered by the presence of warm water in the fjords, and surface melt at the top of the ice sheet is exacerbating this problem. Here, we quantify the vulnerability of northwestern Greenland to further warming using a numerical model. We find that in current conditions, this sector alone will contribute more than 1 cm to sea rise level by 2100, and up to 3 cm in the most extreme scenario.
Sébastien Le clec'h, Sylvie Charbit, Aurélien Quiquet, Xavier Fettweis, Christophe Dumas, Masa Kageyama, Coraline Wyard, and Catherine Ritz
The Cryosphere, 13, 373–395, https://doi.org/10.5194/tc-13-373-2019, https://doi.org/10.5194/tc-13-373-2019, 2019
Short summary
Short summary
Quantifying the future contribution of the Greenland ice sheet (GrIS) to sea-level rise in response to atmospheric changes is important but remains challenging. For the first time a full representation of the feedbacks between a GrIS model and a regional atmospheric model was implemented. The authors highlight the fundamental need for representing the GrIS topography change feedbacks with respect to the atmospheric component face to the strong impact on the projected sea-level rise.
Julien Brondex, Fabien Gillet-Chaulet, and Olivier Gagliardini
The Cryosphere, 13, 177–195, https://doi.org/10.5194/tc-13-177-2019, https://doi.org/10.5194/tc-13-177-2019, 2019
Short summary
Short summary
Here, we apply a synthetic perturbation to the most active drainage basin of Antarctica and show that centennial mass loss projections obtained through ice flow models depend strongly on the implemented friction law, i.e. the mathematical relationship between basal drag and sliding velocities. In particular, the commonly used Weertman law considerably underestimates the sea-level contribution of this basin in comparison to two water pressure-dependent laws which rely on stronger physical bases.
Hongju Yu, Eric Rignot, Helene Seroussi, and Mathieu Morlighem
The Cryosphere, 12, 3861–3876, https://doi.org/10.5194/tc-12-3861-2018, https://doi.org/10.5194/tc-12-3861-2018, 2018
Short summary
Short summary
Thwaites Glacier, West Antarctica, has experienced rapid grounding line retreat and mass loss in the past decades. In this study, we simulate the evolution of Thwaites Glacier over the next century using different model configurations. Overall, we estimate a 5 mm contribution to global sea level rise from Thwaites Glacier in the next 30 years. However, a 300 % uncertainty is found over the next 100 years, ranging from 14 to 42 mm, depending on the model setup.
Youngmin Choi, Mathieu Morlighem, Michael Wood, and Johannes H. Bondzio
The Cryosphere, 12, 3735–3746, https://doi.org/10.5194/tc-12-3735-2018, https://doi.org/10.5194/tc-12-3735-2018, 2018
Short summary
Short summary
Calving is an important mechanism that controls the dynamics of Greenland outlet glaciers. We test and compare four calving laws and assess which calving law has better predictive abilities. Overall, the calving law based on von Mises stress is more satisfactory than other laws, but new parameterizations should be derived to better capture the detailed processes involved in calving.
Rupert M. Gladstone, Yuwei Xia, and John Moore
The Cryosphere, 12, 3605–3615, https://doi.org/10.5194/tc-12-3605-2018, https://doi.org/10.5194/tc-12-3605-2018, 2018
Short summary
Short summary
Computer models for the simulation of marine ice sheets (ice sheets resting on bedrock below sea level) historically show poor numerical convergence for grounding line (the boundary between grounded and floating parts of the ice sheet) movement. We have further characterised the nature of the numerical problems leading to poor convergence and highlighted implications for the design of computer experiments that test grounding line movement.
Hélène Seroussi and Mathieu Morlighem
The Cryosphere, 12, 3085–3096, https://doi.org/10.5194/tc-12-3085-2018, https://doi.org/10.5194/tc-12-3085-2018, 2018
Michael J. Wolovick and John C. Moore
The Cryosphere, 12, 2955–2967, https://doi.org/10.5194/tc-12-2955-2018, https://doi.org/10.5194/tc-12-2955-2018, 2018
Short summary
Short summary
In this paper, we explore the possibility of using locally targeted geoengineering to slow the rate of an ice sheet collapse. We find that an intervention as big as existing large civil engineering projects could have a 30 % probability of stopping an ice sheet collapse, while larger interventions have better odds of success. With more research to improve upon the simple designs we considered, it may be possible to perfect a design that was both achievable and had good odds of success.
Chen Zhao, Rupert M. Gladstone, Roland C. Warner, Matt A. King, Thomas Zwinger, and Mathieu Morlighem
The Cryosphere, 12, 2637–2652, https://doi.org/10.5194/tc-12-2637-2018, https://doi.org/10.5194/tc-12-2637-2018, 2018
Short summary
Short summary
A combination of computer modelling and observational data were used to infer the resistance to ice flow at the bed of the Fleming Glacier on the Antarctic Peninsula. The model was also used to simulate the distribution of temperature within the ice, which governs the rate at which the ice can deform. This is especially important for glaciers like the Fleming Glacier, which has both regions of rapid deformation and regions of rapid sliding at the bed.
Chen Zhao, Rupert M. Gladstone, Roland C. Warner, Matt A. King, Thomas Zwinger, and Mathieu Morlighem
The Cryosphere, 12, 2653–2666, https://doi.org/10.5194/tc-12-2653-2018, https://doi.org/10.5194/tc-12-2653-2018, 2018
Short summary
Short summary
A combination of computer modelling and observational data were used to infer the resistance to ice flow at the bed of the Fleming Glacier on the Antarctic Peninsula in both 2008 and 2015. The comparison suggests the grounding line retreated by ~ 9 km from 2008 to 2015. The retreat may be enhanced by a positive feedback between friction, melting and sliding at the glacier bed.
Lenneke M. Jong, Rupert M. Gladstone, Benjamin K. Galton-Fenzi, and Matt A. King
The Cryosphere, 12, 2425–2436, https://doi.org/10.5194/tc-12-2425-2018, https://doi.org/10.5194/tc-12-2425-2018, 2018
Short summary
Short summary
We used an ice sheet model to simulate temporary regrounding of a marine ice sheet retreating across a retrograde bedrock slope. We show that a sliding relation incorporating water-filled cavities and the ice overburden pressure at the base allows the temporary regrounding to occur. This suggests that choice of basal sliding relation can be important when modelling grounding line behaviour of regions where potential ice rises and pinning points are present and regrounding could occur.
Clemens Schannwell, Stephen Cornford, David Pollard, and Nicholas E. Barrand
The Cryosphere, 12, 2307–2326, https://doi.org/10.5194/tc-12-2307-2018, https://doi.org/10.5194/tc-12-2307-2018, 2018
Short summary
Short summary
Despite the speculation on the state and fate of Larsen C Ice Shelf, a key unknown factor remains: what would be the effects of ice-shelf collapse on upstream drainage basins and thus global sea levels? In our paper three state-of-the-art numerical ice-sheet models were used to simulate the volume evolution of the inland ice sheet to ice-shelf collapse at Larsen C and George VI ice shelves. Our results suggest sea-level rise of up to ~ 4 mm for Larsen C ice shelf and ~ 22 for George VI ice shelf.
Nadine Steiger, Kerim H. Nisancioglu, Henning Åkesson, Basile de Fleurian, and Faezeh M. Nick
The Cryosphere, 12, 2249–2266, https://doi.org/10.5194/tc-12-2249-2018, https://doi.org/10.5194/tc-12-2249-2018, 2018
Short summary
Short summary
We use an ice flow model to reconstruct the retreat of Jakobshavn Isbræ since 1850, forced by increased ocean warming and calving. Fjord geometry governs locations of rapid retreat: narrow and shallow areas act as intermittent pinning points for decades, followed by delayed rapid retreat without additional climate warming. These areas may be used to locate potential moraine buildup. Evidently, historic retreat and geometric influences have to be analyzed individually for each glacier system.
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...
