Articles | Volume 17, issue 2
https://doi.org/10.5194/gmd-17-899-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-899-2024
© Author(s) 2024. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Development and technical paper
| 
02 Feb 2024
Development and technical paper |
| 02 Feb 2024
| 02 Feb 2024
Graphics-processing-unit-accelerated ice flow solver for unstructured meshes using the Shallow-Shelf Approximation (FastIceFlo v1.0.1)
Department of Mechanical Engineering, University of North Dakota, North Dakota, USA
Ludovic Räss
Laboratory of Hydraulics, Hydrology and Glaciology (VAW), ETH Zurich, Zurich, Switzerland
Swiss Federal Institute for Forest, Snow and Landscape Research (WSL), Birmensdorf, Switzerland
now at: Swiss Geocomputing Centre, Faculty of Geosciences and Environment, University of Lausanne, Lausanne, Switzerland
Mathieu Morlighem
Department of Earth Sciences, Dartmouth College, New Hampshire, USA
Related authors
No articles found.
Tim van den Akker, Mathieu Morlighem, and Nina Kirchner
EGUsphere, https://doi.org/10.5194/egusphere-2026-6156, https://doi.org/10.5194/egusphere-2026-6156, 2026
This preprint is open for discussion and under review for The Cryosphere (TC).
Short summary
Short summary
Ice sheet models initialised from observations often suffer an initial “shock” when starting future simulations. The usual fix, a relaxation period where the model settles freely, lets the ice drift from observed conditions and may never reach a steady state. We introduce nudging: gently steering the model during relaxation to remove drift, keep the ice close to observations, and start with the observed thinning rates. We demonstrate it on C.H. Ostenfeld Glacier, northern Greenland.
Adam J. Hepburn, Sammie Buzzard, Andrew J. Sole, Stephen J. Livingstone, Felix Ng, Mathieu Morlighem, Elizabeth A. Bagshaw, Caroline Clason, Lisa Craw, Christine F. Dow, Samuel Doyle, Jonathan Hawkins, Matthew Peacey, and Robert Storrar
The Cryosphere, 20, 5675–5696, https://doi.org/10.5194/tc-20-5675-2026, https://doi.org/10.5194/tc-20-5675-2026, 2026
Short summary
Short summary
Lakes at glacier margins can drain beneath the ice in devastating, floods called jökulhlaups. We present the first 2D model of jökulhlaup propagation, first applying it to a synthetic glacier, before then assessing its ability to reproduce an observed record of lake fill-drain cycles in Greenland. Our results match a 17-year record of flood timing but underpredicts peak discharge likely because we are missing certain physical processes. Our results hold promise for simulating jökulhlaups.
Jamie Barnett, Felicity A. Holmes, Sarah L. Greenwood, Mathieu Morlighem, Nina Kirchner, and Martin Jakobsson
The Cryosphere, 20, 3599–3617, https://doi.org/10.5194/tc-20-3599-2026, https://doi.org/10.5194/tc-20-3599-2026, 2026
Short summary
Short summary
Computer models used to predict future change of the Greenland Ice Sheet are uncertain, especially in how they represent iceberg calving. We compare several calving approaches by testing model results against satellite observations of changes at three unique floating ice shelves in Greenland. We then extend the simulations to the year 2300 to explore future ice loss, finding that warming of the atmosphere or ocean is more important than the choice of calving method.
Thibault Duretz, Albert de Montserrat, Rubén Sevilla, Ludovic Räss, Ivan Utkin, and Arne Spang
Geosci. Model Dev., 19, 5343–5362, https://doi.org/10.5194/gmd-19-5343-2026, https://doi.org/10.5194/gmd-19-5343-2026, 2026
Short summary
Short summary
Geodynamic modeling helps scientists understand how the Earth deforms. New computer methods make these simulations faster and more efficient, especially on powerful computers. They automatically adjust settings for better performance and can handle complex materials and flow types. This approach makes it easier to study large, detailed models of Earth processes.
Thibault Duretz, Albert de Montserrat, Rubén Sevilla, Ludovic Räss, Ivan Utkin, and Arne Spang
Geosci. Model Dev., 19, 5343–5362, https://doi.org/10.5194/gmd-19-5343-2026, https://doi.org/10.5194/gmd-19-5343-2026, 2026
Short summary
Short summary
Geodynamic modeling helps scientists understand how the Earth deforms. New computer methods make these simulations faster and more efficient, especially on powerful computers. They automatically adjust settings for better performance and can handle complex materials and flow types. This approach makes it easier to study large, detailed models of Earth processes.
Mansa Krishna, Gong Cheng, and Mathieu Morlighem
The Cryosphere, 20, 3533–3558, https://doi.org/10.5194/tc-20-3533-2026, https://doi.org/10.5194/tc-20-3533-2026, 2026
Short summary
Short summary
Estimates of the Greenland Ice Sheet’s contribution to sea level rise are affected by uncertainties in the bed topography. Traditional, physics-based methods for inferring the bed elevation are limited to fast-flowing areas of the ice sheet. We use machine learning models informed with two physical laws to infer the bed elevation for different regions in Greenland, showing that this method can be used to infer the bed elevation in slower-moving, sparsely surveyed regions of the ice sheet.
Daniel Abele, Thomas Kleiner, Yannic Fischler, Benjamin Uekermann, Gerasimos Chourdakis, Mathieu Morlighem, Achim Basermann, Christian Bischof, Hans-Joachim Bungartz, and Angelika Humbert
Geosci. Model Dev., 19, 5019–5039, https://doi.org/10.5194/gmd-19-5019-2026, https://doi.org/10.5194/gmd-19-5019-2026, 2026
Short summary
Short summary
Accurate simulations of ice-sheet evolution require a multi-physics approach that captures the interaction between glaciers and their subglacial environment. For this purpose, we have extended ice sheet and subglacial hydrology models using a coupling library. This simplifies coupled Earth systems simulations by allowing models to interact with minimal effort and computational expense. We have evaluated the solution regarding performance and choice of coupling library.
James R. Jordan, Frank Pattyn, Daniel Abele, Torsten Albrecht, Jorge Alvarez-Solas, Jowan M. Barnes, Tijn Berends, Jorge A. Bernales, Javier Blasco, Gong Cheng, Youngmin Choi, Stephen L. Cornford, Cruz Garcia-Molina, Fabien Gillet-Chaulet, G. Hilmar Gudmundsson, Angelika Humbert, Gunter R. Leguy, William H. Lipscomb, Marisa Montoya, Daniel Moreno-Parada, Mathieu Morlighem, Tyler Pelle, Alexander Robinson, Martin Rückamp, Hélène Seroussi, Yanmei Tian, Luisa Wagner, Roderick S. W. van de Wal, Liyun Zhao, and Thomas Zwinger
EGUsphere, https://doi.org/10.5194/egusphere-2026-1962, https://doi.org/10.5194/egusphere-2026-1962, 2026
Short summary
Short summary
Ice calving is a recent inclusion in ice sheet models and there has not been a systematic test of how well it has been implemented. We show that models can accurately represent a given rate of ice calving, properties at the ice front evolve smoothly during calving, and that there are no consistent differences in ice behaviour between various approaches to representing calving in numerical ice models. This gives us confidence in their ability for use in sea level rise prediction simulations.
Ankit Pramanik, Sarah L. Greenwood, Mathieu Morlighem, Jamie Barnett, Felicity A. Holmes, Richard Gyllencreutz, and Carl Regnéll
EGUsphere, https://doi.org/10.5194/egusphere-2026-172, https://doi.org/10.5194/egusphere-2026-172, 2026
Short summary
Short summary
Ice sheet/glacier retreat is driving expansion of glacial lakes impounded at the ice margin. Using an ice sheet model, we find that lake filling and drainage causes systematic changes to glacier thickness, velocity and mass loss. Rising lake levels dramatically increase mass loss as the ice margin becomes buoyant, while rapid lake drainage can drive crevassing and destabilise the margin. Steady lake drainage, in contrast, has a stabilising effect, and may mitigate mass loss and flood hazards.
Arne Spang, Marcel Thielmann, Casper Pranger, Albert de Montserrat, and Ludovic Räss
Geosci. Model Dev., 19, 369–388, https://doi.org/10.5194/gmd-19-369-2026, https://doi.org/10.5194/gmd-19-369-2026, 2026
Short summary
Short summary
Concentration of deformation is difficult to capture accurately in computer simulations. We present a number of challenges associated with concentrated viscous deformation and demonstrate strategies to overcome them. The strategies include automatic selection of appropriate time steps to react to rapid changes in model behavior, automatic rescaling to avoid rounding errors, and three methods to prevent model instability. This way, we are able to accurately capture very fast viscous deformation.
Johanna Beckmann, Ronja Reese, Felicity S. McCormack, Sue Cook, Lawrence Bird, Dawid Gwyther, Daniel Richards, Matthias Scheiter, Yu Wang, Hélène Seroussi, Ayako Abe‐Ouchi, Torsten Albrecht, Jorge Alvarez‐Solas, Xylar S. Asay‐Davis, Jean‐Baptiste Barre, Constantijn J. Berends, Jorge Bernales, Javier Blasco, Justine Caillet, David M. Chandler, Violaine Coulon, Richard Cullather, Christophe Dumas, Benjamin K. Galton‐Fenzi, Julius Garbe, Fabien Gillet‐Chaulet, Rupert Gladstone, Heiko Goelzer, Nicholas R. Golledge, Ralf Greve, G. Hilmar Gudmundsson, Holly Kyeore Han, Trevor R. Hillebrand, Matthew J. Hoffman, Philippe Huybrechts, Nicolas C. Jourdain, Ann Kristin Klose, Petra M. Langebroek, Gunter R. Leguy, William H. Lipscomb, Daniel P. Lowry, Pierre Mathiot, Marisa Montoya, Mathieu Morlighem, Sophie Nowicki, Frank Pattyn, Antony J. Payne, Tyler Pelle, Aurélien Quiquet, Alexander Robinson, Leopekka Saraste, Erika G. Simon, Sainan Sun, Jake P. Twarog, Luke D. Trusel, Benoit Urruty, Jonas Van Breedam, Roderik S. W. van de Wal, Chen Zhao, and Thomas Zwinger
EGUsphere, https://doi.org/10.5194/egusphere-2025-4069, https://doi.org/10.5194/egusphere-2025-4069, 2025
Short summary
Short summary
Antarctica holds enough ice to raise sea levels by many meters, but its future is uncertain. Warm ocean water melts ice shelves from below, letting inland ice flow faster into the sea. By 2300, Antarctica could add 0.6–4.4 m to sea levels. Our study identifies two key factors—how strongly shelves melt and how the ice responds. These explain much of the range, and refining them in models may improve future predictions.
Jamie Barnett, Felicity A. Holmes, Joshua Cuzzone, Henning Åkesson, Mathieu Morlighem, Matt O'Regan, Johan Nilsson, Nina Kirchner, and Martin Jakobsson
The Cryosphere, 19, 3631–3653, https://doi.org/10.5194/tc-19-3631-2025, https://doi.org/10.5194/tc-19-3631-2025, 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 the 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 is required to predict future changes and reliable sea level rise estimates.
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.
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.
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.
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.
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.
Ludovic Räss, Ivan Utkin, Thibault Duretz, Samuel Omlin, and Yuri Y. Podladchikov
Geosci. Model Dev., 15, 5757–5786, https://doi.org/10.5194/gmd-15-5757-2022, https://doi.org/10.5194/gmd-15-5757-2022, 2022
Short summary
Short summary
Continuum mechanics-based modelling of physical processes at large scale requires huge computational resources provided by massively parallel hardware such as graphical processing units. We present a suite of numerical algorithms, implemented using the Julia language, that efficiently leverages the parallelism. We demonstrate that our implementation is efficient, scalable and robust and showcase applications to various geophysical problems.
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.
Cited articles
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
Brædstrup, C. F., Damsgaard, A., and Egholm, D. L.: Ice-sheet modelling accelerated by graphics cards, Comput. Geosci., 72, 210–220, https://doi.org/10.1016/j.cageo.2014.07.019, 2014. a
Castleman, B. A., Schlegel, N.-J., Caron, L., Larour, E., and Khazendar, A.: Derivation of bedrock topography measurement requirements for the reduction of uncertainty in ice-sheet model projections of Thwaites Glacier, The Cryosphere, 16, 761–778, https://doi.org/10.5194/tc-16-761-2022, 2022. a
Chen, X., Zhang, X., Church, J. A., Watson, C. S., King, M. A., Monselesan, D., Legresy, B., and Harig, C.: The increasing rate of global mean sea-level rise during 1993–2014, Nat. Clim. Change, 7, 492–495, https://doi.org/10.1038/nclimate3325, 2017. a
Cornford, S. L., Martin, D. F., Graves, D. T., Ranken, D. F., Le Brocq, A. M., Gladstone, R. M., Payne, A. J., Ng, E. G., and Lipscomb, W. H.: Adaptive mesh, finite volume modeling of marine ice sheets, J. Comput. Phys., 232, 529–549, https://doi.org/10.1016/j.jcp.2012.08.037, 2013. a
Duretz, T., Räss, L., Podladchikov, Y., and Schmalholz, S.: Resolving thermomechanical coupling in two and three dimensions: spontaneous strain localization owing to shear heating, Geophys. J. Int., 216, 365–379, https://doi.org/10.1093/gji/ggy434, 2019. a, b
Gilbarg, D. and Trudinger, N. S.: Elliptic partial differential equations of second order, vol. 224, Springer, https://doi.org/10.1007/978-3-642-61798-0, 1977. a
Glen, J. W.: The creep of polycrystalline ice, P. Roy. Soc. Lond.-Ser. A, 228, 519–538, https://doi.org/10.1098/rspa.1955.0066, 1955. a
Goelzer, H., Robinson, A., Seroussi, H., and Van De Wal, R. S.: Recent progress in Greenland ice sheet modelling, Curr. Clim. Change Rep., 3, 291–302, https://doi.org/10.1007/s40641-017-0073-y, 2017. a
Habbal, F., Larour, E., Morlighem, M., Seroussi, H., Borstad, C. P., and Rignot, E.: Optimal numerical solvers for transient simulations of ice flow using the Ice Sheet System Model (ISSM versions 4.2.5 and 4.11), Geosci. Model Dev., 10, 155–168, https://doi.org/10.5194/gmd-10-155-2017, 2017. a
Häfner, D., Nuterman, R., and Jochum, M.: Fast, cheap, and turbulent – Global ocean modeling with GPU acceleration in python, J. Adv. Model. Earth Sy., 13, e2021MS002717, https://doi.org/10.1029/2021MS002717, 2021. a
Hinkel, J., Lincke, D., Vafeidis, A. T., Perrette, M., Nicholls, R. J., Tol, R. S., Marzeion, B., Fettweis, X., Ionescu, C., and Levermann, A.: Coastal flood damage and adaptation costs under 21st century sea-level rise, P. Natl. Acad. Sci. USA, 111, 3292–3297, https://doi.org/10.1073/pnas.1222469111, 2014. a
IPCC: Climate Change 2021 The Physical Science Basis, Working Group 1 (WG1) Contribution to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, Tech. rep., Intergovernmental Panel on Climate Change, 2021. a
Jouvet, G., Cordonnier, G., Kim, B., Lüthi, M., Vieli, A., and Aschwanden, A.: Deep learning speeds up ice flow modelling by several orders of magnitude, J. Glaciol., 68, 651–664, https://doi.org/10.1017/jog.2021.120, 2022. a
Kelley, C. and Liao, L.-Z.: Explicit pseudo-transient continuation, Computing, 15, 18, 2013. a
Kopp, R. E., Kemp, A. C., Bittermann, K., Horton, B. P., Donnelly, J. P., Gehrels, W. R., Hay, C. C., Mitrovica, J. X., Morrow, E. D., and Rahmstorf, S.: Temperature-driven global sea-level variability in the Common Era, P. Natl. Acad. Sci. USA, 113, E1434–E1441, https://doi.org/10.1073/pnas.1517056113, 2016. 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.-Earth, 117, 22, https://doi.org/10.1029/2011JF002140, 2012. a, b, c
MacAyeal, D. R.: Large-scale ice flow over a viscous basal sediment: Theory and application to ice stream B, Antarctica, J. Geophys. Res.-Sol. Ea., 94, 4071–4087, 1989. 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., Noel, 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 multibeam echo sounding combined with mass conservation, Geophys. Res. Lett., 44, 11–51, https://doi.org/10.1002/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., 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., Van Den Broeke, M. R., Van Ommen, T. D., Wessem, M. V., and Young, D. A.: Deep glacial troughs and stabilizing ridges unveiled beneath the margins of the Antarctic ice sheet, Nat. Geosci., 13, 132–137, https://doi.org/10.1038/s41561-019-0510-8, 2020. a
Omlin, S., Räss, L., and Podladchikov, Y. Y.: Simulation of three-dimensional viscoelastic deformation coupled to porous fluid flow, Tectonophysics, 746, 695–701, https://doi.org/10.1016/j.tecto.2017.08.012, 2018. a
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
Poliakov, A. N., Cundall, P. A., Podladchikov, Y. Y., and Lyakhovsky, V. A.: An explicit inertial method for the simulation of viscoelastic flow: an evaluation of elastic effects on diapiric flow in two- and three- layers models, Flow and creep in the solar system: observations, modeling and theory, in: Flow and Creep in the Solar System: Observations, Modeling and Theory, edited by: Stone, D. B. and Runcorn, S. K., NATO ASI Series, vol. 391, Springer, Dordrecht, 175–195, https://doi.org/10.1007/978-94-015-8206-3_12, 1993. a
Räss, L., Duretz, T., and Podladchikov, Y.: Resolving hydromechanical coupling in two and three dimensions: spontaneous channelling of porous fluids owing to decompaction weakening, Geophys. J. Int., 218, 1591–1616, https://doi.org/10.1093/gji/ggz239, 2019. a, b, c
Räss, L., Licul, A., Herman, F., Podladchikov, Y. Y., and Suckale, J.: Modelling thermomechanical ice deformation using an implicit pseudo-transient method (FastICE v1.0) based on graphical processing units (GPUs), Geosci. Model Dev., 13, 955–976, https://doi.org/10.5194/gmd-13-955-2020, 2020. a, b, c, d, e, f, g, h
Räss, L., Utkin, I., Duretz, T., Omlin, S., and Podladchikov, Y. Y.: Assessing the robustness and scalability of the accelerated pseudo-transient method, Geosci. Model Dev., 15, 5757–5786, https://doi.org/10.5194/gmd-15-5757-2022, 2022. a, b, c
Rietbroek, R., Brunnabend, S.-E., Kusche, J., Schröter, J., and Dahle, C.: Revisiting the contemporary sea-level budget on global and regional scales, P. Natl. Acad. Sci. USA, 113, 1504–1509, https://doi.org/10.1073/pnas.1519132113, 2016. a
Rignot, E., Mouginot, J., and Scheuchl, B.: Antarctic grounding line mapping from differential satellite radar interferometry, Geophys. Res. Lett., 38, 49, https://doi.org/10.1029/2011GL047109, 2011. a
Sandip, A., Morlighem, M., and Räss, L.: AnjaliSandip/FastIceFlo: FastIceFlo v1.0.1 (v1.0.1), Zenodo [code], https://doi.org/10.5281/zenodo.8356351, 2023. a
Tezaur, I. K., Perego, M., Salinger, A. G., Tuminaro, R. S., and Price, S. F.: Albany/FELIX: a parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis, Geosci. Model Dev., 8, 1197–1220, https://doi.org/10.5194/gmd-8-1197-2015, 2015. a
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.
We solve momentum balance for unstructured meshes to predict ice flow for real glaciers using a...
