Articles | Volume 17, issue 16
https://doi.org/10.5194/gmd-17-6227-2024
© Author(s) 2024. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/gmd-17-6227-2024
© Author(s) 2024. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Numerical stabilization methods for level-set-based ice front migration
Department of Earth Sciences, Dartmouth College, Hanover, NH 03755, USA
Mathieu Morlighem
Department of Earth Sciences, Dartmouth College, Hanover, NH 03755, USA
G. Hilmar Gudmundsson
Department of Geography and Environmental Sciences, Northumbria University, Newcastle upon Tyne, NE1 8ST, UK
Related authors
Gong Cheng, Mansa Krishna, and Mathieu Morlighem
Geosci. Model Dev., 18, 5311–5327, https://doi.org/10.5194/gmd-18-5311-2025, https://doi.org/10.5194/gmd-18-5311-2025, 2025
Short summary
Short summary
Predicting ice sheet contributions to sea level rise is challenging due to limited data and uncertainties in key processes. Traditional models require complex methods that lack flexibility. We developed PINNICLE (Physics-Informed Neural Networks for Ice and CLimatE), an open-source Python library that integrates machine learning with physical laws to improve ice sheet modeling. By combining data and physics, PINNICLE enhances predictions and adaptability, providing a powerful tool for climate research and sea level rise projections.
Martin Rückamp, Gong Cheng, Karlheinz Gutjahr, Marco Möller, Petri K. E. Pellikka, and Christoph Mayer
EGUsphere, https://doi.org/10.5194/egusphere-2025-3150, https://doi.org/10.5194/egusphere-2025-3150, 2025
This preprint is open for discussion and under review for The Cryosphere (TC).
Short summary
Short summary
The study simulates the 21st-century evolution of Great Aletsch Glacier and Hintereisferner using full-Stokes ice dynamics and surface mass balance under different emission scenarios. Results show significant ice loss, with Hintereisferner expected to disappear by mid-century. Great Aletsch Glacier vanish by the end of the century under high-emission scenarios, but persist under lower-emission scenarios. These trends agree with large-scale models except some variability.
Younghyun Koo, Gong Cheng, Mathieu Morlighem, and Maryam Rahnemoonfar
The Cryosphere, 19, 2583–2599, https://doi.org/10.5194/tc-19-2583-2025, https://doi.org/10.5194/tc-19-2583-2025, 2025
Short summary
Short summary
Calving, the breaking of ice bodies from the terminus of a glacier, plays an important role in the mass losses of Greenland ice sheets. However, calving parameters have been poorly understood because of the intensive computational demands of traditional numerical models. To address this issue and find the optimal calving parameter that best represents real observations, we develop deep-learning emulators based on graph neural network architectures.
Shfaqat A. Khan, Helene Seroussi, Mathieu Morlighem, William Colgan, Veit Helm, Gong Cheng, Danjal Berg, Valentina R. Barletta, Nicolaj K. Larsen, William Kochtitzky, Michiel van den Broeke, Kurt H. Kjær, Andy Aschwanden, Brice Noël, Jason E. Box, Joseph A. MacGregor, Robert S. Fausto, Kenneth D. Mankoff, Ian M. Howat, Kuba Oniszk, Dominik Fahrner, Anja Løkkegaard, Eigil Y. H. Lippert, Alicia Bråtner, and Javed Hassan
Earth Syst. Sci. Data, 17, 3047–3071, https://doi.org/10.5194/essd-17-3047-2025, https://doi.org/10.5194/essd-17-3047-2025, 2025
Short summary
Short summary
The surface elevation of the Greenland Ice Sheet is changing due to surface mass balance processes and ice dynamics, each exhibiting distinct spatiotemporal patterns. Here, we employ satellite and airborne altimetry data with fine spatial (1 km) and temporal (monthly) resolutions to document this spatiotemporal evolution from 2003 to 2023. This dataset of fine-resolution altimetry data in both space and time will support studies of ice mass loss and be useful for GIS ice sheet modeling.
Joel A. Wilner, Mathieu Morlighem, and Gong Cheng
The Cryosphere, 17, 4889–4901, https://doi.org/10.5194/tc-17-4889-2023, https://doi.org/10.5194/tc-17-4889-2023, 2023
Short summary
Short summary
We use numerical modeling to study iceberg calving off of ice shelves in Antarctica. We examine four widely used mathematical descriptions of calving (
calving laws), under the assumption that Antarctic ice shelf front positions should be in steady state under the current climate forcing. We quantify how well each of these calving laws replicates the observed front positions. Our results suggest that the eigencalving and von Mises laws are most suitable for Antarctic ice shelves.
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.
Gong Cheng, Mansa Krishna, and Mathieu Morlighem
Geosci. Model Dev., 18, 5311–5327, https://doi.org/10.5194/gmd-18-5311-2025, https://doi.org/10.5194/gmd-18-5311-2025, 2025
Short summary
Short summary
Predicting ice sheet contributions to sea level rise is challenging due to limited data and uncertainties in key processes. Traditional models require complex methods that lack flexibility. We developed PINNICLE (Physics-Informed Neural Networks for Ice and CLimatE), an open-source Python library that integrates machine learning with physical laws to improve ice sheet modeling. By combining data and physics, PINNICLE enhances predictions and adaptability, providing a powerful tool for climate research and sea level rise projections.
Martin Rückamp, Gong Cheng, Karlheinz Gutjahr, Marco Möller, Petri K. E. Pellikka, and Christoph Mayer
EGUsphere, https://doi.org/10.5194/egusphere-2025-3150, https://doi.org/10.5194/egusphere-2025-3150, 2025
This preprint is open for discussion and under review for The Cryosphere (TC).
Short summary
Short summary
The study simulates the 21st-century evolution of Great Aletsch Glacier and Hintereisferner using full-Stokes ice dynamics and surface mass balance under different emission scenarios. Results show significant ice loss, with Hintereisferner expected to disappear by mid-century. Great Aletsch Glacier vanish by the end of the century under high-emission scenarios, but persist under lower-emission scenarios. These trends agree with large-scale models except some variability.
Felicity A. Holmes, Jamie Barnett, Henning Åkesson, Mathieu Morlighem, Johan Nilsson, Nina Kirchner, and Martin Jakobsson
The Cryosphere, 19, 2695–2714, https://doi.org/10.5194/tc-19-2695-2025, https://doi.org/10.5194/tc-19-2695-2025, 2025
Short summary
Short summary
Northern Greenland contains some of the ice sheet's last remaining glaciers with floating ice tongues. One of these is Ryder Glacier, which has been relatively stable in recent decades, in contrast to nearby glaciers. Here, we use a computer model to simulate Ryder Glacier until 2300 under both a low- and a high-emissions scenario. Very high levels of surface melt under a high-emissions future lead to a sea level rise contribution that is an order of magnitude higher than under a low-emissions future.
Daniel Abele, Thomas Kleiner, Yannic Fischler, Benjamin Uekermann, Gerasimos Chourdakis, Mathieu Morlighem, Achim Basermann, Christian Bischof, Hans-Joachim Bungartz, and Angelika Humbert
EGUsphere, https://doi.org/10.5194/egusphere-2025-3345, https://doi.org/10.5194/egusphere-2025-3345, 2025
This preprint is open for discussion and under review for Geoscientific Model Development (GMD).
Short summary
Short summary
For accurate projections of the evolution of continental ice sheets in Greenland and Antartica, interactions between the ice and its environment must be included in simulations. For this purpose, we have implemented adapters for the ice sheet model ISSM and subglacial hydrology model CUAS-MPI for the coupling library preCICE. This simplifies the study of earth systems by allowing the models to interact with each other as well as with models of the oceans or atmosphere with very little effort.
Younghyun Koo, Gong Cheng, Mathieu Morlighem, and Maryam Rahnemoonfar
The Cryosphere, 19, 2583–2599, https://doi.org/10.5194/tc-19-2583-2025, https://doi.org/10.5194/tc-19-2583-2025, 2025
Short summary
Short summary
Calving, the breaking of ice bodies from the terminus of a glacier, plays an important role in the mass losses of Greenland ice sheets. However, calving parameters have been poorly understood because of the intensive computational demands of traditional numerical models. To address this issue and find the optimal calving parameter that best represents real observations, we develop deep-learning emulators based on graph neural network architectures.
Sebastian H. R. Rosier, G. Hilmar Gudmundsson, Adrian Jenkins, and Kaitlin A. Naughten
The Cryosphere, 19, 2527–2557, https://doi.org/10.5194/tc-19-2527-2025, https://doi.org/10.5194/tc-19-2527-2025, 2025
Short summary
Short summary
Glaciers in the Amundsen Sea region of Antarctica have been retreating and losing mass, but their future contribution to global sea level rise remains highly uncertain. We use an ice sheet model and uncertainty quantification methods to evaluate the probable range of mass loss from this region for two future climate scenarios. We find that the rate of ice loss until 2100 will likely remain similar to present-day observations, with little sensitivity to climate scenario over this short time frame.
Shfaqat A. Khan, Helene Seroussi, Mathieu Morlighem, William Colgan, Veit Helm, Gong Cheng, Danjal Berg, Valentina R. Barletta, Nicolaj K. Larsen, William Kochtitzky, Michiel van den Broeke, Kurt H. Kjær, Andy Aschwanden, Brice Noël, Jason E. Box, Joseph A. MacGregor, Robert S. Fausto, Kenneth D. Mankoff, Ian M. Howat, Kuba Oniszk, Dominik Fahrner, Anja Løkkegaard, Eigil Y. H. Lippert, Alicia Bråtner, and Javed Hassan
Earth Syst. Sci. Data, 17, 3047–3071, https://doi.org/10.5194/essd-17-3047-2025, https://doi.org/10.5194/essd-17-3047-2025, 2025
Short summary
Short summary
The surface elevation of the Greenland Ice Sheet is changing due to surface mass balance processes and ice dynamics, each exhibiting distinct spatiotemporal patterns. Here, we employ satellite and airborne altimetry data with fine spatial (1 km) and temporal (monthly) resolutions to document this spatiotemporal evolution from 2003 to 2023. This dataset of fine-resolution altimetry data in both space and time will support studies of ice mass loss and be useful for GIS ice sheet modeling.
Ole Richter, Ralph Timmermann, G. Hilmar Gudmundsson, and Jan De Rydt
Geosci. Model Dev., 18, 2945–2960, https://doi.org/10.5194/gmd-18-2945-2025, https://doi.org/10.5194/gmd-18-2945-2025, 2025
Short summary
Short summary
The new coupled ice sheet–ocean model addresses challenges related to horizontal resolution through advanced mesh flexibility, enabled by the use of unstructured grids. We describe the new model, verify its functioning in an idealised setting and demonstrate its advantages in a global-ocean–Antarctic ice sheet domain. The results of this study comprise an important step towards improving predictions of the Antarctic contribution to sea level rise over centennial timescales.
Joshua K. Cuzzone, Aaron Barth, Kelsey Barker, and Mathieu Morlighem
The Cryosphere, 19, 1559–1575, https://doi.org/10.5194/tc-19-1559-2025, https://doi.org/10.5194/tc-19-1559-2025, 2025
Short summary
Short summary
We use an ice sheet model to simulate the Last Glacial Maximum conditions of the Laurentide Ice Sheet (LIS) across the northeastern United States. A complex thermal history existed for the LIS that caused high erosion across most of the NE USA but prevented erosion across high-elevation mountain peaks and areas where ice flow was slow. This has implications for geologic studies which rely on the erosional nature of the LIS to reconstruct its glacial history and landscape evolution.
Jowan M. Barnes, G. Hilmar Gudmundsson, Daniel N. Goldberg, and Sainan Sun
EGUsphere, https://doi.org/10.5194/egusphere-2025-328, https://doi.org/10.5194/egusphere-2025-328, 2025
Short summary
Short summary
Calving is where ice breaks off the front of glaciers. It has not been included widely in modelling as it is difficult to represent. We use our ice flow model to investigate the effects of calving floating ice shelves in West Antarctica. More calving leads to more ice loss and greater sea level rise, with local differences due to the shape of the bedrock. We find that ocean forcing and calving should be considered equally when trying to improve how models represent the real world.
Jamie Barnett, Felicity Alice Holmes, Joshua Cuzzone, Henning Åkesson, Mathieu Morlighem, Matt O'Regan, Johan Nilsson, Nina Kirchner, and Martin Jakobsson
EGUsphere, https://doi.org/10.5194/egusphere-2025-653, https://doi.org/10.5194/egusphere-2025-653, 2025
Short summary
Short summary
Understanding how ice sheets have changed in the past can allow us to make better predictions for the future. By running a state-of-the-art model of Ryder Glacier, North Greenland, over the past 12,000 years we find that both a warming atmosphere and ocean play a key role in the evolution of the Glacier. Our conclusions stress that accurately quantifying the ice sheet’s interactions with the ocean are required to predict future changes and reliable sea level rise estimates.
Francesca Baldacchino, Nicholas R. Golledge, Mathieu Morlighem, Huw Horgan, Alanna V. Alevropoulos-Borrill, Alena Malyarenko, Alexandra Gossart, Daniel P. Lowry, and Laurine van Haastrecht
The Cryosphere, 19, 107–127, https://doi.org/10.5194/tc-19-107-2025, https://doi.org/10.5194/tc-19-107-2025, 2025
Short summary
Short summary
Understanding how the Ross Ice Shelf flow is changing in a warming world is important for predicting ice sheet change. Field measurements show clear intra-annual variations in ice flow; however, it is unclear what mechanisms drive this variability. We show that local perturbations in basal melt can have a significant impact on ice flow speed, but a combination of forcings is likely driving the observed variability in ice flow.
Richard Parsons, Sainan Sun, G. Hilmar Gudmundsson, Jan Wuite, and Thomas Nagler
The Cryosphere, 18, 5789–5801, https://doi.org/10.5194/tc-18-5789-2024, https://doi.org/10.5194/tc-18-5789-2024, 2024
Short summary
Short summary
In 2022, multi-year landfast sea ice in Antarctica's Larsen B embayment disintegrated, after which time an increase in the rate at which Crane Glacier discharged ice into the ocean was observed. As the fast ice was joined to the glacier terminus, it could provide resistance against the glacier's flow, slowing down the rate of ice discharge. We used numerical modelling to quantify this resistive stress and found that the fast ice provided significant support to Crane prior to its disintegration.
Xianwei Wang, Hilmar Gudmundsson, and David Holland
EGUsphere, https://doi.org/10.5194/egusphere-2024-2790, https://doi.org/10.5194/egusphere-2024-2790, 2024
Short summary
Short summary
Understanding why iceberg calved during drifting in the ocean is important to understand the life cycle and the influence on the surrounding ocean of an iceberg. This study explains why iceberg A68a calved when approaching the South Georgia Island in late 2020 during its drifting in the Southern Ocean using satellite observation and modeling, which was caused by collision with seamount.
Brad Reed, J. A. Mattias Green, Adrian Jenkins, and G. Hilmar Gudmundsson
The Cryosphere, 18, 4567–4587, https://doi.org/10.5194/tc-18-4567-2024, https://doi.org/10.5194/tc-18-4567-2024, 2024
Short summary
Short summary
We use a numerical ice-flow model to simulate the response of a 1940s Pine Island Glacier to changes in melting beneath its ice shelf. A decadal period of warm forcing is sufficient to push the glacier into an unstable, irreversible retreat from its long-term position on a subglacial ridge to an upstream ice plain. This retreat can only be stopped when unrealistic cold forcing is applied. These results show that short warm anomalies can lead to quick and substantial increases in ice flux.
J. Rachel Carr, Emily A. Hill, and G. Hilmar Gudmundsson
The Cryosphere, 18, 2719–2737, https://doi.org/10.5194/tc-18-2719-2024, https://doi.org/10.5194/tc-18-2719-2024, 2024
Short summary
Short summary
The Greenland Ice Sheet is one of the world's largest glaciers and is melting quickly in response to climate change. It contains fast-flowing channels of ice that move ice from Greenland's centre to its coasts and allow Greenland to react quickly to climate warming. As a result, we want to predict how these glaciers will behave in the future, but there are lots of uncertainties. Here we assess the impacts of two main sources of uncertainties in glacier models.
Cristina Gerli, Sebastian Rosier, G. Hilmar Gudmundsson, and Sainan Sun
The Cryosphere, 18, 2677–2689, https://doi.org/10.5194/tc-18-2677-2024, https://doi.org/10.5194/tc-18-2677-2024, 2024
Short summary
Short summary
Recent efforts have focused on using AI and satellite imagery to track crevasses for assessing ice shelf damage and informing ice flow models. Our study reveals a weak connection between these observed products and damage maps inferred from ice flow models. While there is some improvement in crevasse-dense regions, this association remains limited. Directly mapping ice damage from satellite observations may not significantly improve the representation of these processes within ice flow models.
In-Woo Park, Emilia Kyung Jin, Mathieu Morlighem, and Kang-Kun Lee
The Cryosphere, 18, 1139–1155, https://doi.org/10.5194/tc-18-1139-2024, https://doi.org/10.5194/tc-18-1139-2024, 2024
Short summary
Short summary
This study conducted 3D thermodynamic ice sheet model experiments, and modeled temperatures were compared with 15 observed borehole temperature profiles. We found that using incompressibility of ice without sliding agrees well with observed temperature profiles in slow-flow regions, while incorporating sliding in fast-flow regions captures observed temperature profiles. Also, the choice of vertical velocity scheme has a greater impact on the shape of the modeled temperature profile.
Anjali Sandip, Ludovic Räss, and Mathieu Morlighem
Geosci. Model Dev., 17, 899–909, https://doi.org/10.5194/gmd-17-899-2024, https://doi.org/10.5194/gmd-17-899-2024, 2024
Short summary
Short summary
We solve momentum balance for unstructured meshes to predict ice flow for real glaciers using a pseudo-transient method on graphics processing units (GPUs) and compare it to a standard central processing unit (CPU) implementation. We justify the GPU implementation by applying the price-to-performance metric for up to million-grid-point spatial resolutions. This study represents a first step toward leveraging GPU processing power, enabling more accurate polar ice discharge predictions.
Youngmin Choi, Helene Seroussi, Mathieu Morlighem, Nicole-Jeanne Schlegel, and Alex Gardner
The Cryosphere, 17, 5499–5517, https://doi.org/10.5194/tc-17-5499-2023, https://doi.org/10.5194/tc-17-5499-2023, 2023
Short summary
Short summary
Ice sheet models are often initialized using snapshot observations of present-day conditions, but this approach has limitations in capturing the transient evolution of the system. To more accurately represent the accelerating changes in glaciers, we employed time-dependent data assimilation. We found that models calibrated with the transient data better capture past trends and more accurately reproduce changes after the calibration period, even with limited observations.
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.
Joel A. Wilner, Mathieu Morlighem, and Gong Cheng
The Cryosphere, 17, 4889–4901, https://doi.org/10.5194/tc-17-4889-2023, https://doi.org/10.5194/tc-17-4889-2023, 2023
Short summary
Short summary
We use numerical modeling to study iceberg calving off of ice shelves in Antarctica. We examine four widely used mathematical descriptions of calving (
calving laws), under the assumption that Antarctic ice shelf front positions should be in steady state under the current climate forcing. We quantify how well each of these calving laws replicates the observed front positions. Our results suggest that the eigencalving and von Mises laws are most suitable for Antarctic ice shelves.
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.
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.
Alice C. Frémand, Peter Fretwell, Julien A. Bodart, Hamish D. Pritchard, Alan Aitken, Jonathan L. Bamber, Robin Bell, Cesidio Bianchi, Robert G. Bingham, Donald D. Blankenship, Gino Casassa, Ginny Catania, Knut Christianson, Howard Conway, Hugh F. J. Corr, Xiangbin Cui, Detlef Damaske, Volkmar Damm, Reinhard Drews, Graeme Eagles, Olaf Eisen, Hannes Eisermann, Fausto Ferraccioli, Elena Field, René Forsberg, Steven Franke, Shuji Fujita, Yonggyu Gim, Vikram Goel, Siva Prasad Gogineni, Jamin Greenbaum, Benjamin Hills, Richard C. A. Hindmarsh, Andrew O. Hoffman, Per Holmlund, Nicholas Holschuh, John W. Holt, Annika N. Horlings, Angelika Humbert, Robert W. Jacobel, Daniela Jansen, Adrian Jenkins, Wilfried Jokat, Tom Jordan, Edward King, Jack Kohler, William Krabill, Mette Kusk Gillespie, Kirsty Langley, Joohan Lee, German Leitchenkov, Carlton Leuschen, Bruce Luyendyk, Joseph MacGregor, Emma MacKie, Kenichi Matsuoka, Mathieu Morlighem, Jérémie Mouginot, Frank O. Nitsche, Yoshifumi Nogi, Ole A. Nost, John Paden, Frank Pattyn, Sergey V. Popov, Eric Rignot, David M. Rippin, Andrés Rivera, Jason Roberts, Neil Ross, Anotonia Ruppel, Dustin M. Schroeder, Martin J. Siegert, Andrew M. Smith, Daniel Steinhage, Michael Studinger, Bo Sun, Ignazio Tabacco, Kirsty Tinto, Stefano Urbini, David Vaughan, Brian C. Welch, Douglas S. Wilson, Duncan A. Young, and Achille Zirizzotti
Earth Syst. Sci. Data, 15, 2695–2710, https://doi.org/10.5194/essd-15-2695-2023, https://doi.org/10.5194/essd-15-2695-2023, 2023
Short summary
Short summary
This paper presents the release of over 60 years of ice thickness, bed elevation, and surface elevation data acquired over Antarctica by the international community. These data are a crucial component of the Antarctic Bedmap initiative which aims to produce a new map and datasets of Antarctic ice thickness and bed topography for the international glaciology and geophysical community.
Sebastian H. R. Rosier, Christopher Y. S. Bull, Wai L. Woo, and G. Hilmar Gudmundsson
The Cryosphere, 17, 499–518, https://doi.org/10.5194/tc-17-499-2023, https://doi.org/10.5194/tc-17-499-2023, 2023
Short summary
Short summary
Future ice loss from Antarctica could raise sea levels by several metres, and key to this is the rate at which the ocean melts the ice sheet from below. Existing methods for modelling this process are either computationally expensive or very simplified. We present a new approach using machine learning to mimic the melt rates calculated by an ocean model but in a fraction of the time. This approach may provide a powerful alternative to existing methods, without compromising on accuracy or speed.
Bertie W. J. Miles, Chris R. Stokes, Adrian Jenkins, Jim R. Jordan, Stewart S. R. Jamieson, and G. Hilmar Gudmundsson
The Cryosphere, 17, 445–456, https://doi.org/10.5194/tc-17-445-2023, https://doi.org/10.5194/tc-17-445-2023, 2023
Short summary
Short summary
Satellite observations have shown that the Shirase Glacier catchment in East Antarctica has been gaining mass over the past 2 decades, a trend largely attributed to increased snowfall. Our multi-decadal observations of Shirase Glacier show that ocean forcing has also contributed to some of this recent mass gain. This has been caused by strengthening easterly winds reducing the inflow of warm water underneath the Shirase ice tongue, causing the glacier to slow down and thicken.
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.
Francesca Baldacchino, Mathieu Morlighem, Nicholas R. Golledge, Huw Horgan, and Alena Malyarenko
The Cryosphere, 16, 3723–3738, https://doi.org/10.5194/tc-16-3723-2022, https://doi.org/10.5194/tc-16-3723-2022, 2022
Short summary
Short summary
Understanding how the Ross Ice Shelf will evolve in a warming world is important to the future stability of Antarctica. It remains unclear what changes could drive the largest mass loss in the future and where places are most likely to trigger larger mass losses. Sensitivity maps are modelled showing that the RIS is sensitive to changes in environmental and glaciological controls at regions which are currently experiencing changes. These regions need to be monitored in a warming world.
Joshua K. Cuzzone, Nicolás E. Young, Mathieu Morlighem, Jason P. Briner, and Nicole-Jeanne Schlegel
The Cryosphere, 16, 2355–2372, https://doi.org/10.5194/tc-16-2355-2022, https://doi.org/10.5194/tc-16-2355-2022, 2022
Short summary
Short summary
We use an ice sheet model to determine what influenced the Greenland Ice Sheet to retreat across a portion of southwestern Greenland during the Holocene (about the last 12 000 years). Our simulations, constrained by observations from geologic markers, show that atmospheric warming and ice melt primarily caused the ice sheet to retreat rapidly across this domain. We find, however, that iceberg calving at the interface where the ice meets the ocean significantly influenced ice mass change.
Yannic Fischler, Martin Rückamp, Christian Bischof, Vadym Aizinger, Mathieu Morlighem, and Angelika Humbert
Geosci. Model Dev., 15, 3753–3771, https://doi.org/10.5194/gmd-15-3753-2022, https://doi.org/10.5194/gmd-15-3753-2022, 2022
Short summary
Short summary
Ice sheet models are used to simulate the changes of ice sheets in future but are currently often run in coarse resolution and/or with neglecting important physics to make them affordable in terms of computational costs. We conducted a study simulating the Greenland Ice Sheet in high resolution and adequate physics to test where the ISSM ice sheet code is using most time and what could be done to improve its performance for future computer architectures that allow massive parallel computing.
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.
Thomas Frank, Henning Åkesson, Basile de Fleurian, Mathieu Morlighem, and Kerim H. Nisancioglu
The Cryosphere, 16, 581–601, https://doi.org/10.5194/tc-16-581-2022, https://doi.org/10.5194/tc-16-581-2022, 2022
Short summary
Short summary
The shape of a fjord can promote or inhibit glacier retreat in response to climate change. We conduct experiments with a synthetic setup under idealized conditions in a numerical model to study and quantify the processes involved. We find that friction between ice and fjord is the most important factor and that it is possible to directly link ice discharge and grounding line retreat to fjord topography in a quantitative way.
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.
Emily A. Hill, Sebastian H. R. Rosier, G. Hilmar Gudmundsson, and Matthew Collins
The Cryosphere, 15, 4675–4702, https://doi.org/10.5194/tc-15-4675-2021, https://doi.org/10.5194/tc-15-4675-2021, 2021
Short summary
Short summary
Using an ice flow model and uncertainty quantification methods, we provide probabilistic projections of future sea level rise from the Filchner–Ronne region of Antarctica. We find that it is most likely that this region will contribute negatively to sea level rise over the next 300 years, largely as a result of increased surface mass balance. We identify parameters controlling ice shelf melt and snowfall contribute most to uncertainties in projections.
Matt O'Regan, Thomas M. Cronin, Brendan Reilly, Aage Kristian Olsen Alstrup, Laura Gemery, Anna Golub, Larry A. Mayer, Mathieu Morlighem, Matthias Moros, Ole L. Munk, Johan Nilsson, Christof Pearce, Henrieka Detlef, Christian Stranne, Flor Vermassen, Gabriel West, and Martin Jakobsson
The Cryosphere, 15, 4073–4097, https://doi.org/10.5194/tc-15-4073-2021, https://doi.org/10.5194/tc-15-4073-2021, 2021
Short summary
Short summary
Ryder Glacier is a marine-terminating glacier in north Greenland discharging ice into the Lincoln Sea. Here we use marine sediment cores to reconstruct its retreat and advance behavior through the Holocene. We show that while Sherard Osborn Fjord has a physiography conducive to glacier and ice tongue stability, Ryder still retreated more than 40 km inland from its current position by the Middle Holocene. This highlights the sensitivity of north Greenland's marine glaciers to climate change.
Thiago Dias dos Santos, Mathieu Morlighem, and Hélène Seroussi
Geosci. Model Dev., 14, 2545–2573, https://doi.org/10.5194/gmd-14-2545-2021, https://doi.org/10.5194/gmd-14-2545-2021, 2021
Short summary
Short summary
Numerical models are routinely used to understand the past and future behavior of ice sheets in response to climate evolution. As is always the case with numerical modeling, one needs to minimize biases and numerical artifacts due to the choice of numerical scheme employed in such models. Here, we assess different numerical schemes in time-dependent simulations of ice sheets. We also introduce a new parameterization for the driving stress, the force that drives the ice sheet 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.
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.
Bertie W. J. Miles, Jim R. Jordan, Chris R. Stokes, Stewart S. R. Jamieson, G. Hilmar Gudmundsson, and Adrian Jenkins
The Cryosphere, 15, 663–676, https://doi.org/10.5194/tc-15-663-2021, https://doi.org/10.5194/tc-15-663-2021, 2021
Short summary
Short summary
We provide a historical overview of changes in Denman Glacier's flow speed, structure and calving events since the 1960s. Based on these observations, we perform a series of numerical modelling experiments to determine the likely cause of Denman's acceleration since the 1970s. We show that grounding line retreat, ice shelf thinning and the detachment of Denman's ice tongue from a pinning point are the most likely causes of the observed acceleration.
Jan De Rydt, Ronja Reese, Fernando S. Paolo, and G. Hilmar Gudmundsson
The Cryosphere, 15, 113–132, https://doi.org/10.5194/tc-15-113-2021, https://doi.org/10.5194/tc-15-113-2021, 2021
Short summary
Short summary
We used satellite observations and numerical simulations of Pine Island Glacier, West Antarctica, between 1996 and 2016 to show that the recent increase in its flow speed can only be reproduced by computer models if stringent assumptions are made about the material properties of the ice and its underlying bed. These assumptions are not commonly adopted in ice flow modelling, and our results therefore have implications for future simulations of Antarctic ice flow and sea level projections.
Kate Winter, Emily A. Hill, G. Hilmar Gudmundsson, and John Woodward
Earth Syst. Sci. Data, 12, 3453–3467, https://doi.org/10.5194/essd-12-3453-2020, https://doi.org/10.5194/essd-12-3453-2020, 2020
Short summary
Short summary
Satellite measurements of the English Coast in the Antarctic Peninsula reveal that glaciers are thinning and losing mass, but ice thickness data are required to assess these changes, in terms of ice flux and sea level contribution. Our ice-penetrating radar measurements reveal that low-elevation subglacial channels control fast-flowing ice streams, which release over 39 Gt of ice per year to floating ice shelves. This topography could make ice flows susceptible to future instability.
Xiangbin Cui, Hafeez Jeofry, Jamin S. Greenbaum, Jingxue Guo, Lin Li, Laura E. Lindzey, Feras A. Habbal, Wei Wei, Duncan A. Young, Neil Ross, Mathieu Morlighem, Lenneke M. Jong, Jason L. Roberts, Donald D. Blankenship, Sun Bo, and Martin J. Siegert
Earth Syst. Sci. Data, 12, 2765–2774, https://doi.org/10.5194/essd-12-2765-2020, https://doi.org/10.5194/essd-12-2765-2020, 2020
Short summary
Short summary
We present a topographic digital elevation model (DEM) for Princess Elizabeth Land (PEL), East Antarctica. The DEM covers an area of approximately 900 000 km2 and was built from radio-echo sounding data collected in four campaigns since 2015. Previously, to generate the Bedmap2 topographic product, PEL’s bed was characterised from low-resolution satellite gravity data across an otherwise large (>200 km wide) data-free zone.
Eric Larour, Lambert Caron, Mathieu Morlighem, Surendra Adhikari, Thomas Frederikse, Nicole-Jeanne Schlegel, Erik Ivins, Benjamin Hamlington, Robert Kopp, and Sophie Nowicki
Geosci. Model Dev., 13, 4925–4941, https://doi.org/10.5194/gmd-13-4925-2020, https://doi.org/10.5194/gmd-13-4925-2020, 2020
Short summary
Short summary
ISSM-SLPS is a new projection system for future sea level that increases the resolution and accuracy of current projection systems and improves the way uncertainty is treated in such projections. This will pave the way for better inclusion of state-of-the-art results from existing intercomparison efforts carried out by the scientific community, such as GlacierMIP2 or ISMIP6, into sea-level projections.
Martin Rückamp, Angelika Humbert, Thomas Kleiner, Mathieu Morlighem, and Helene Seroussi
Geosci. Model Dev., 13, 4491–4501, https://doi.org/10.5194/gmd-13-4491-2020, https://doi.org/10.5194/gmd-13-4491-2020, 2020
Short summary
Short summary
We present enthalpy formulations within the Ice-Sheet and Sea-Level System model that show better performance than earlier implementations. A first experiment indicates that the treatment of discontinuous conductivities of the solid–fluid system with a geometric mean produce accurate results when applied to coarse vertical resolutions. In a second experiment, we propose a novel stabilization formulation that avoids the problem of thin elements. This method provides accurate and stable results.
Heiko Goelzer, Sophie Nowicki, Anthony Payne, Eric Larour, Helene Seroussi, William H. Lipscomb, Jonathan Gregory, Ayako Abe-Ouchi, Andrew Shepherd, Erika Simon, Cécile Agosta, Patrick Alexander, Andy Aschwanden, Alice Barthel, Reinhard Calov, Christopher Chambers, Youngmin Choi, Joshua Cuzzone, Christophe Dumas, Tamsin Edwards, Denis Felikson, Xavier Fettweis, Nicholas R. Golledge, Ralf Greve, Angelika Humbert, Philippe Huybrechts, Sebastien Le clec'h, Victoria Lee, Gunter Leguy, Chris Little, Daniel P. Lowry, Mathieu Morlighem, Isabel Nias, Aurelien Quiquet, Martin Rückamp, Nicole-Jeanne Schlegel, Donald A. Slater, Robin S. Smith, Fiamma Straneo, Lev Tarasov, Roderik van de Wal, and Michiel van den Broeke
The Cryosphere, 14, 3071–3096, https://doi.org/10.5194/tc-14-3071-2020, https://doi.org/10.5194/tc-14-3071-2020, 2020
Short summary
Short summary
In this paper we use a large ensemble of Greenland ice sheet models forced by six different global climate models to project ice sheet changes and sea-level rise contributions over the 21st century.
The results for two different greenhouse gas concentration scenarios indicate that the Greenland ice sheet will continue to lose mass until 2100, with contributions to sea-level rise of 90 ± 50 mm and 32 ± 17 mm for the high (RCP8.5) and low (RCP2.6) scenario, respectively.
Hélène Seroussi, Sophie Nowicki, Antony J. Payne, Heiko Goelzer, William H. Lipscomb, Ayako Abe-Ouchi, Cécile Agosta, Torsten Albrecht, Xylar Asay-Davis, Alice Barthel, Reinhard Calov, Richard Cullather, Christophe Dumas, Benjamin K. Galton-Fenzi, Rupert Gladstone, Nicholas R. Golledge, Jonathan M. Gregory, Ralf Greve, Tore Hattermann, Matthew J. Hoffman, Angelika Humbert, Philippe Huybrechts, Nicolas C. Jourdain, Thomas Kleiner, Eric Larour, Gunter R. Leguy, Daniel P. Lowry, Chistopher M. Little, Mathieu Morlighem, Frank Pattyn, Tyler Pelle, Stephen F. Price, Aurélien Quiquet, Ronja Reese, Nicole-Jeanne Schlegel, Andrew Shepherd, Erika Simon, Robin S. Smith, Fiammetta Straneo, Sainan Sun, Luke D. Trusel, Jonas Van Breedam, Roderik S. W. van de Wal, Ricarda Winkelmann, Chen Zhao, Tong Zhang, and Thomas Zwinger
The Cryosphere, 14, 3033–3070, https://doi.org/10.5194/tc-14-3033-2020, https://doi.org/10.5194/tc-14-3033-2020, 2020
Short summary
Short summary
The Antarctic ice sheet has been losing mass over at least the past 3 decades in response to changes in atmospheric and oceanic conditions. This study presents an ensemble of model simulations of the Antarctic evolution over the 2015–2100 period based on various ice sheet models, climate forcings and emission scenarios. Results suggest that the West Antarctic ice sheet will continue losing a large amount of ice, while the East Antarctic ice sheet could experience increased snow accumulation.
Cited articles
Benn, D. I. and Astrom, J. A.: Calving glaciers and ice shelves, Adv. Phys.-X, 3, 1513819, https://doi.org/10.1080/23746149.2018.1513819, 2018. a
Biswas, R., Devine, K. D., and Flahert, J. E.: Parallel, adaptive finite element methods for conservation laws, Appl. Numer. Math., 14, 255–283, https://doi.org/10.1016/0168-9274(94)90029-9, 1994. a
Black, T. E. and Joughin, I.: Weekly to monthly terminus variability of Greenland's marine-terminating outlet glaciers, The Cryosphere, 17, 1–13, https://doi.org/10.5194/tc-17-1-2023, 2023. a, b
Bondzio, J., Morlighem, M., Seroussi, H., Kleiner, T., Ruckamp, M., Mouginot, J., Moon, T., Larour, E., and Humbert, A.: The mechanisms behind Jakobshavn Isbræ's acceleration and mass loss: A 3-D thermomechanical model study, Geophys. Res. Lett., 44, 6252–6260, https://doi.org/10.1002/2017GL073309, 2017. a
Bondzio, J. H., Seroussi, H., Morlighem, M., Kleiner, T., Rückamp, M., Humbert, A., and Larour, E. Y.: Modelling calving front dynamics using a level-set method: application to Jakobshavn Isbræ, West Greenland, The Cryosphere, 10, 497–510, https://doi.org/10.5194/tc-10-497-2016, 2016. a, b, c
Cheng, D., Hayes, W., Larour, E., Mohajerani, Y., Wood, M., Velicogna, I., and Rignot, E.: Calving Front Machine (CALFIN): glacial termini dataset and automated deep learning extraction method for Greenland, 1972–2019, The Cryosphere, 15, 1663–1675, https://doi.org/10.5194/tc-15-1663-2021, 2021. a
Cheng, G.: Animations of the evolution of misfits, Zenodo [video], https://doi.org/10.5281/zenodo.10454554, 2024. a
Cheng, G., Morlighem, M., Mouginot, J., and Cheng, D.: Helheim Glacier's Terminus Position Controls Its Seasonal and Inter-Annual Ice Flow Variability, Geophys. Res. Lett., 49, e2021GL097085, https://doi.org/10.1029/2021GL097085, 2022. a
Choi, Y., Morlighem, M., Wood, M., and Bondzio, J. H.: Comparison of four calving laws to model Greenland outlet glaciers, The Cryosphere, 12, 3735–3746, https://doi.org/10.5194/tc-12-3735-2018, 2018. a
Choi, Y., Morlighem, M., Rignot, E., and Wood, M.: Ice dynamics will remain a primary driver of Greenland ice sheet mass loss over the next century, Nature Commun. Earth Environ., 2, 26, https://doi.org/10.1038/s43247-021-00092-z, 2021. a, b
Crawford, A. J., Benn, D. I., Todd, J., Astrom, J. A., Bassis, J. N., and Zwinger, T.: Marine ice-cliff instability modeling shows mixed- mode ice-cliff failure and yields calving rate parameterization, Nat. Commun., 12, 2701, https://doi.org/10.1038/s41467-021-23070-7, 2021. a
dos Santos, T. D., Morlighem, M., and Seroussi, H.: Assessment of numerical schemes for transient, finite-element ice flow models using ISSM v4.18, Geosci. Model Dev., 14, 2545–2573, https://doi.org/10.5194/gmd-14-2545-2021, 2021. a, b, c
Eriksson, K., Estep, D., Hansbo, P., and Johnson, C.: Computational differential equations, Cambridge University Press, ISBN 9780521567381, 1996. a
Gibou, F., Fedkiw, R., and Osher, S.: A review of level-set methods and some recent applications, J. Comput. Phys., 353, 82–109, https://doi.org/10.1016/j.jcp.2017.10.006, 2018. a
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
Gong, C.: enigne/Levelset: Level-set Stabilization (v1.1), Zenodo [code and data set], https://doi.org/10.5281/zenodo.10454657, 2024. a
Greve, R. and Blatter, H.: Dynamics of Ice Sheets and Glaciers, Advances in Geophysical and Environmental Mechanics and Mathematics, Springer Science & Business Media, https://doi.org/10.1007/978-3-642-03415-2, 2009. a
Gross, S. and Reusken, A.: Numerical Methods for Two-phase Incompressible Flows, vol. 40 of Springer Series in Computational Mathematics, 1–472, https://doi.org/10.1007/978-3-642-19686-7, 2011. a
Gudmundsson, G. H., Paolo, F. S., Adusumilli, S., and Fricker, H. A.: Instantaneous Antarctic ice sheet mass loss driven by thinning ice shelves, Geophys. Res. Lett., 46, 13903–13909, https://doi.org/10.1029/2019GL085027, 2019. a, b
Gudmundsson, H.: GHilmarG/UaSource: Ua2019b, Zenodo [code], https://doi.org/10.5281/zenodo.3706624, 2020. a, b, c
Hartmann, D., Meinke, M., and Schroeder, W.: The constrained reinitialization equation for level set methods, J. Comput. Phys., 229, 1514–1535, https://doi.org/10.1016/j.jcp.2009.10.042, 2010. a, b
Henri, F., Coquerelle, M., and Lubin, P.: Geometrical level set reinitialization using closest point method and kink detection for thin filaments, topology changes and two-phase flows, J. Comput. Phys., 448, 110704, https://doi.org/10.1016/j.jcp.2021.110704, 2022. a, b
IPCC: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, vol. In Press, Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, https://doi.org/10.1017/9781009157896, 2021. a
ISSM Team: Ice-sheet and Sea-level System Model source code, v4.23 r27696, Zenodo [code], https://doi.org/10.5281/zenodo.7850841, 2023. a, b
King, M. D., Howat, I. M., Candela, S. G., Noh, M. J., Jeong, S., Noel, B. P. Y., van den Broeke, M. R., Wouters, B., and Negrete, A.: Dynamic ice loss from the Greenland Ice Sheet driven by sustained glacier retreat, Commun. Earth Environ., 1, 1, https://doi.org/10.1038/s43247-020-0001-2, 2020. a
Larour, E., Seroussi, H., Morlighem, M., and Rignot, E.: Continental scale, high order, high spatial resolution, ice sheet modeling using the Ice Sheet System Model (ISSM), J. Geophys. Res., 117, 1–20, https://doi.org/10.1029/2011JF002140, 2012. a
Larson, M. G. and Bengzon, F.: The Finite Element Method: Theory, Implementation, and Applications, Springer Publishing Company, Incorporated, https://doi.org/10.1007/978-3-642-33287-6, 2013. a, b, c, d
Li, C., Xu, C., Gui, C., and Fox, M. D.: Level set evolution without re-initialization: A new variational formulation, in: 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition, Vol 1, Proceedings, edited by: Schmid, C., Soatto, S., and Tomasi, C., PROCEEDINGS – IEEE Computer Society Conference On Computer Vision And Pattern Recognition, 430–436, IEEE Comp Soc, conference on Computer Vision and Pattern Recognition, 20–25 June 2005, San Diego, CA, 2005. 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
Mouginot, J., Rignot, E., Bjørk, A. A., van den Broeke, M., Millan, R., Morlighem, M., Noël, B., Scheuchl, B., and Wood, M.: Forty-six years of Greenland Ice Sheet mass balance from 1972 to 2018, P. Natl. Acad. Sci. USA., 116, 9239–9244, https://doi.org/10.1073/pnas.1904242116, 2019. a, b
Pattyn, F. and Morlighem, M.: The uncertain future of the Antarctic Ice Sheet, Science, 367, 1331–1335, https://doi.org/10.1126/science.aaz5487, 2020. a
Pollard, D., DeConto, R. M., and Alley, R. B.: Potential Antarctic Ice Sheet retreat driven by hydrofracturing and ice cliff failure, Earth Planet Sc. Lett., 412, 112–121, https://doi.org/10.1016/j.epsl.2014.12.035, 2015. a
Seroussi, H., Nowicki, S., Payne, A. J., Goelzer, H., Lipscomb, W. H., Abe-Ouchi, A., Agosta, C., Albrecht, T., Asay-Davis, X., Barthel, A., Calov, R., Cullather, R., Dumas, C., Galton-Fenzi, B. K., Gladstone, R., Golledge, N. R., Gregory, J. M., Greve, R., Hattermann, T., Hoffman, M. J., Humbert, A., Huybrechts, P., Jourdain, N. C., Kleiner, T., Larour, E., Leguy, G. R., Lowry, D. P., Little, C. M., Morlighem, M., Pattyn, F., Pelle, T., Price, S. F., Quiquet, A., Reese, R., Schlegel, N.-J., Shepherd, A., Simon, E., Smith, R. S., Straneo, F., Sun, S., Trusel, L. D., Van Breedam, J., van de Wal, R. S. W., Winkelmann, R., Zhao, C., Zhang, T., and Zwinger, T.: ISMIP6 Antarctica: a multi-model ensemble of the Antarctic ice sheet evolution over the 21st century, The Cryosphere, 14, 3033–3070, https://doi.org/10.5194/tc-14-3033-2020, 2020. a
Sussman, M., Smereka, P., and Osher, S.: A Level set Aspproach for Computing Solutions to Incompressible Two-Phase Flow, J. Comput. Phys., 114, 146–159, https://doi.org/10.1006/jcph.1994.1155, 1994. a
Toure, M. K. and Soulaimani, A.: Stabilized finite element methods for solving the level set equation without reinitialization, Comput. Math. Appl., 71, 1602–1623, https://doi.org/10.1016/j.camwa.2016.02.028, 2016. a, b
Wood, M., Rignot, E., Fenty, I., An, L., Bjørk, A., van den Broeke, M., Cai, C., Kane, E., Menemenlis, D., Millan, R., Morlighem, M., Mouginot, J., Noël, B., Scheuchl, B., Velicogna, I., Willis, J. K., and Zhang, H.: Ocean forcing drives glacier retreat in Greenland, Sci. Adv., 7, 1315–1332, https://doi.org/10.1126/sciadv.aba7282, 2021. a
Short summary
We conducted a comprehensive analysis of the stabilization and reinitialization techniques currently employed in ISSM and Úa for solving level-set equations, specifically those related to the dynamic representation of moving ice fronts within numerical ice sheet models. Our results demonstrate that the streamline upwind Petrov–Galerkin (SUPG) method outperforms the other approaches. We found that excessively frequent reinitialization can lead to exceptionally high errors in simulations.
We conducted a comprehensive analysis of the stabilization and reinitialization techniques...