Articles | Volume 18, issue 16
https://doi.org/10.5194/gmd-18-5311-2025
© Author(s) 2025. 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-18-5311-2025
© Author(s) 2025. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Model description paper
| 
26 Aug 2025
Model description paper |
| 26 Aug 2025
| 26 Aug 2025
A Python library for solving ice sheet modeling problems using physics-informed neural networks, PINNICLE v1.0
Department of Earth Sciences, Dartmouth College, Hanover, NH 03755, USA
Mansa Krishna
Department of Earth Sciences, Dartmouth College, Hanover, NH 03755, USA
Mathieu Morlighem
Department of Earth Sciences, Dartmouth College, Hanover, NH 03755, USA
Related authors
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.
Gong Cheng, Mathieu Morlighem, and G. Hilmar Gudmundsson
Geosci. Model Dev., 17, 6227–6247, https://doi.org/10.5194/gmd-17-6227-2024, https://doi.org/10.5194/gmd-17-6227-2024, 2024
Short summary
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.
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.
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.
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.
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.
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.
Gong Cheng, Mathieu Morlighem, and G. Hilmar Gudmundsson
Geosci. Model Dev., 17, 6227–6247, https://doi.org/10.5194/gmd-17-6227-2024, https://doi.org/10.5194/gmd-17-6227-2024, 2024
Short summary
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.
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.
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.
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.
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.
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.
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.
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
Abadi, M., Agarwal, A., Barham, P., Brevdo, E., Chen, Z., Citro, C., Corrado, G. S., Davis, A., Dean, J., Devin, M., Ghemawat, S., Goodfellow, I., Harp, A., Irving, G., Isard, M., Jia, Y., Jozefowicz, R., Kaiser, L., Kudlur, M., Levenberg, J., Mané, D., Monga, R., Moore, S., Murray, D., Olah, C., Schuster, M., Shlens, J., Steiner, B., Sutskever, I., Talwar, K., Tucker, P., Vanhoucke, V., Vasudevan, V., Viégas, F., Vinyals, O., Warden, P., Wattenberg, M., Wicke, M., Yu, Y., and Zheng, X.: TensorFlow: Large-scale machine learning on heterogeneous systems, https://www.tensorflow.org (last access: 14 August 2025), 2015. a
Aschwanden, A., Fahnestock, M. A., Truffer, M., Brinkerhoff, D. J., Hock, R., Khroulev, K., Mottram, R., and Khan, S. A.: Contribution of the Greenland Ice Sheet to sea level over the next millennium, Sci. Adv., 5, eaav9396, https://doi.org/10.1126/sciadv.aav9396, 2019. a
Aschwanden, A., Bartholomaus, T. C., Brinkerhoff, D. J., and Truffer, M.: Brief communication: A roadmap towards credible projections of ice sheet contribution to sea level, The Cryosphere, 15, 5705–5715, https://doi.org/10.5194/tc-15-5705-2021, 2021. a
Bauer, P., Thorpe, A., and Brunet, G.: The Quiet Revolution of Numerical Weather Prediction, Nature, 525, 47–55, https://doi.org/10.1038/nature14956, 2015. a
Bolibar, J., Sapienza, F., Maussion, F., Lguensat, R., Wouters, B., and Pérez, F.: Universal differential equations for glacier ice flow modelling, Geosci. Model Dev., 16, 6671–6687, https://doi.org/10.5194/gmd-16-6671-2023, 2023. a, b
Bradbury, J., Frostig, R., Hawkins, P., Johnson, M. J., Leary, C., Maclaurin, D., Necula, G., Paszke, A., VanderPlas, J., Wanderman-Milne, S., and Zhang, Q.: JAX: composable transformations of Python+NumPy programs, GitHub [code], http://github.com/jax-ml/jax (last access: 14 August 2025), 2018. a
Brinkerhoff, D. J.: Variational inference at glacier scale, J. Comput. Phys., 459, 111095, https://doi.org/10.1016/j.jcp.2022.111095, 2022. a
Budd, W. F., Keage, P. L., and Blundy, N. A.: Empirical studies of ice sliding, J. Glaciol., 23, 157–170, 1979. a
Cheng, G. and Lotstedt, P.: Parameter sensitivity analysis of dynamic ice sheet models – numerical computations, The Cryosphere, 14, 673–691, https://doi.org/10.5194/tc-14-673-2020, 2020. 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, b
Cheng, G., Morlighem, M., and Francis, S.: Forward and Inverse Modeling of Ice Sheet Flow Using Physics-Informed Neural Networks: Application to Helheim Glacier, Greenland, J. Geophys. Res. Mach. Learn. Comput., 1, e2024JH000169, https://doi.org/10.1029/2024JH000169, 2024. a, b, c, d, e, f, g, h, i, j, k, l, m
Cheng, G., Krishna, M., and Francis, S.: ISSMteam/PINNICLE: Official release version of PINNICLE, Zenodo [code], https://doi.org/10.5281/zenodo.15643042, 2025. a
Cho, J., Nam, S., Yang, H., Yun, S.-B., Hong, Y., and Park, E.: Separable Physics-Informed Neural Networks, arXiv [preprint], https://arxiv.org/abs/2306.15969 (last access: 14 August 2025), 2023. a
Colgan, W., MacGregor, J. A., Mankoff, K. D., Haagenson, R., Rajaram, H., Martos, Y. M., Morlighem, M., Fahnestock, M. A., and Kjeldsen, K. K.: Topographic Correction of Geothermal Heat Flux in Greenland and Antarctica, J. Geophys. Res., 126, e2020JF005598, https://doi.org/10.1029/2020JF005598, 2021. a
dos Santos, T. D., Morlighem, M., and Brinkerhoff, D.: A new vertically integrated, MOno-Layer Higher-Order ice flow model (MOLHO), The Cryosphere, 16, 179–195, https://doi.org/10.5194/tc-16-179-2022, 2022. a
Edwards, T. L., Brandon, M. A., Durand, G., Edwards, N. R., Golledge, N. R., Holden, P. B., Nias, I. J., Payne, A. J., Ritz, C., and Wernecke, A.: Revisiting Antarctic ice loss due to marine ice-cliff instability, Nature, 566, 58–64, https://doi.org/10.1038/s41586-019-0901-4, 2019. a
Fowler, A. C.: A theoretical treatment of the sliding of glaciers in the absence of cavitation, Philos. T. Roy. Soc. A, 298, 637–685, https://doi.org/10.1098/rsta.1981.0003, 1981. a
Furst, J. J., Durand, G., Gillet-Chaulet, F., Merino, N., Tavard, L., Mouginot, J., Gourmelen, N., and Gagliardini, O.: Assimilation of Antarctic velocity observations provides evidence for uncharted pinning points, The Cryosphere, 9, 1427–1443, https://doi.org/10.5194/tc-9-1427-2015, 2015. a
Gagliardini, O., Cohen, D., Raback, P., and Zwinger, T.: Finite-element modeling of subglacial cavities and related friction law, J. Geophys. Res.-Earth, 112, 1–11, https://doi.org/10.1029/2006JF000576, 2007. a, b
Glen, J. W.: The flow law of ice: A discussion of the assumptions made in glacier theory, their experimental foundations and consequences, IASH Publ., 47, 171–183, 1958. 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
Goldberg, D. N. and Heimbach, P.: Parameter and state estimation with a time-dependent adjoint marine ice sheet model, The Cryosphere, 7, 1659–1678, https://doi.org/10.5194/tc-7-1659-2013, 2013. a
Goldberg, D. N. and Sergienko, O. V.: Data assimilation using a hybrid ice flow model, The Cryosphere, 5, 315–327, https://doi.org/10.5194/tc-5-315-2011, 2011. a
He, Q., Perego, M., Howard, A. A., Karniadakis, G. E., and Stinis, P.: A hybrid deep neural operator/finite element method for ice-sheet modeling, J. Comput. Phys., 492, 112428, https://doi.org/10.1016/j.jcp.2023.112428, 2023. a, b
Iwasaki, Y. and Lai, C.-Y.: One-dimensional ice shelf hardness inversion: Clustering behavior and collocation resampling in physics-informed neural networks, J. Comput. Phys., 492, 112435, https://doi.org/10.1016/j.jcp.2023.112435, 2023. a, b, c, d
Joughin, I., Shapero, D., Dutrieux, P., and Smith, B.: Ocean-induced melt volume directly paces ice loss from Pine Island Glacier, Sci. Adv., 7, eabi5738, https://doi.org/10.1126/sciadv.abi5738, 2021. a
Jouvet, G. and Cordonnier, G.: Ice-flow model emulator based on physics-informed deep learning, J. Glaciol., 69, 1941–1955, https://doi.org/10.1017/jog.2023.73, 2023. a, b, c
Karniadakis, G. E., Kevrekidis, I. G., Lu, L., Perdikaris, P., Wang, S., and Yang, L.: Physics-informed machine learning, Nat. Rev. Phys., 3, 422–440, https://doi.org/10.1038/s42254-021-00314-5, 2021. a, b, c
Koo, Y., Cheng, G., Morlighem, M., and Rahnemoonfar, M.: Calibrating calving parameterizations using graph neural network emulators: Application to Helheim Glacier, East Greenland, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2024-1620, 2024. a
Lai, C.-Y., Hassanzadeh, P., Sheshadri, A., Sonnewald, M., Ferrari, R., and Balaji, V.: Machine Learning for Climate Physics and Simulations, Annu. Rev. Condens. Matt. Phys., 16, 343–365, https://doi.org/10.1146/annurev-conmatphys-043024-114758, 2024. 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, b
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
MacAyeal, D. R.: A tutorial on the use of control methods in ice-sheet modeling, J. Glaciol., 39, 91–98, 1993. a
Morlighem, M., Rignot, E., Seroussi, H., Larour, E., Ben Dhia, H., and Aubry, D.: Spatial patterns of basal drag inferred using control methods from a full-Stokes and simpler models for Pine Island Glacier, West Antarctica, Geophys. Res. Lett., 37, 1–6, https://doi.org/10.1029/2010GL043853, 2010. a
Morlighem, M., Seroussi, H., Larour, E., and Rignot, E.: Inversion of basal friction in Antarctica using exact and incomplete adjoints of a higher-order model, J. Geophys. Res., 118, 1746–1753, https://doi.org/10.1002/jgrf.20125, 2013. a
Morlighem, M., Williams, C. N., Rignot, E., An, L., Arndt, J. E., Bamber, J. L., Catania, G., Chauché, N., Dowdeswell, J. A., Dorschel, B., Fenty, I., Hogan, K., Howat, I., Hubbard, A., Jakobsson, M., Jordan, T. M., Kjeldsen, K. K., Millan, R., Mayer, L., Mouginot, J., Noël, B. P. Y., O'Cofaigh, C., Palmer, S., Rysgaard, S., Seroussi, H., Siegert, M. J., Slabon, P., Straneo, F., van den Broeke, M. R., Weinrebe, W., Wood, M., and Zinglersen, K. B.: BedMachine v3: Complete bed topography and ocean bathymetry mapping of Greenland from multi-beam echo sounding combined with mass conservation, Geophys. Res. Lett., 44, 11,051–11,061, https://doi.org/10.1002/2017GL074954, 2017GL074954, 2017. a
Morlighem, M., Rignot, E., Binder, T., Blankenship, D., Drews, R., Eagles, G., Eisen, O., Ferraccioli, F., Forsberg, R., Fretwell, P., Goel, V., Greenbaum, J. S., Gudmundsson, H., Guo, J., Helm, V., Hofstede, C., Howat, I., Humbert, A., Jokat, W., Karlsson, N. B., Lee, W. S., Matsuoka, K., Millan, R., Mouginot, J., Paden, J., Pattyn, F., Roberts, J., Rosier, S., Ruppel, A., Seroussi, H., Smith, E. C., Steinhage, D., Sun, B., v. d. Broeke, M. R., v. Ommen, T. D., v. Wessem, M., and Young, D. A.: Deep glacial troughs and stabilizing ridges unveiled beneath the margins of the Antarctic ice sheet, Nat. Geosci., 13, 132–137, https://doi.org/10.1038/s41561-019-0510-8, 2020. a
Morlighem, M., Goldberg, D., Barnes, J. M., Bassis, J. N., Benn, D. I., Crawford, A. J., Gudmundsson, G. H., and Seroussi, H.: The West Antarctic Ice Sheet may not be vulnerable to marine ice cliff instability during the 21st century, Sci. Adv., 10, eado7794, https://doi.org/10.1126/sciadv.ado7794, 2024. a
Paszke, A., Gross, S., Massa, F., Lerer, A., Bradbury, J., Chanan, G., Killeen, T., Lin, Z., Gimelshein, N., Antiga, L., Desmaison, A., Kopf, A., Yang, E., DeVito, Z., Raison, M., Tejani, A., Chilamkurthy, S., Steiner, B., Fang, L., Bai, J., and Chintala, S.: PyTorch: An Imperative Style, High-Performance Deep Learning Library, arXiv [preprint], https://doi.org/10.48550/arXiv.1912.01703, 2019. a
Raissi, M., Perdikaris, P., and Karniadakis, G. E.: Physics-informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations, J. Comput. Phys., 378, 686–707, https://doi.org/10.1016/j.jcp.2018.10.045, 2019. a, b, c
Raissi, M., Yazdani, A., and Karniadakis, G. E.: Hidden fluid mechanics: Learning velocity and pressure fields from flow visualizations, Science, 367, 6481, 1026–1030, https://doi.org/10.1126/science.aaw4741, 2020. a
Riel, B., Minchew, B., and Bischoff, T.: Data-driven inference of the sechanics of slip slong glacier beds using physics-informed neural networks: case study on Rutford Ice Stream, Antarctica, J. Adv. Model. Earth Syst., 13, e2021MS002621, https://doi.org/10.1029/2021MS002621, 2021. a, b, c
Robel, A. A., Seroussi, H., and Roe, G. H.: Marine ice sheet instability amplifies and skews uncertainty in projections of future sea-level rise, P. Natl. Acad. Sci. USA, 116, 14887–14892, https://doi.org/10.1073/pnas.1904822116, 2019. a
Schoof, C.: Ice sheet grounding line dynamics: Steady states, stability, and hysteresis, J. Geophys. Res., 112, 1–19, https://doi.org/10.1029/2006JF000664, 2007. a
Seo, J.: Solving real-world optimization tasks using physics-informed neural computing, Sci. Rep., 14, 202, https://doi.org/10.1038/s41598-023-49977-3, 2024. 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
Tancik, M., Srinivasan, P., Mildenhall, B., Fridovich-Keil, S., Raghavan, N., Singhal, U., Ramamoorthi, R., Barron, J., and Ng, R.: Fourier features let networks learn high frequency functions in low dimensional domains, Adv. Neural Inf. Process. Syst., 33, 7537–7547, 2020. a
Teisberg, T. O., Schroeder, D. M., and MacKie, E. J.: A machine learning approach to mass-conserving ice thickness interpolation, in: 2021 IEEE International Geoscience and Remote Sensing Symposium IGARSS, Brussels, Belgium, 2021, 8664–8667, https://doi.org/10.1109/IGARSS47720.2021.9555002, 2021. a
Tikhonov, A. N.: On the stability of inverse problems, C. R. (Dokl.) Acad. Sci. SSSR (N.S.), 39, 176–179, 1943. a
Wang, S., Teng, Y., and Perdikaris, P.: Understanding and mitigating gradient flow pathologies in physics-informed neural networks, SIAM J. Sci.Comput., 43, A3055–A3081, 2021. a
Wang, Y., Lai, C.-Y., Prior, D. J., and Cowen-Breen, C.: Deep learning the flow law of Antarctic ice shelves, Science, 387, 1219–1224, https://doi.org/10.1126/science.adp3300, 2025. a, b, c, d
Wernecke, A., Edwards, T. L., Holden, P. B., Edwards, N. R., and Cornford, S. L.: Quantifying the Impact of Bedrock Topography Uncertainty in Pine Island Glacier Projections for This Century, Geophys. Res. Lett., 49, e2021GL096589, https://doi.org/10.1029/2021GL096589, 2023. a
Wong, T.-T., Luk, W.-S., and Heng, P.-A.: Sampling with Hammersley and Halton points, J. Graph. Tools, 2, 9–24, 1997. a
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.
Predicting ice sheet contributions to sea level rise is challenging due to limited data and...
