Articles | Volume 14, issue 5
https://doi.org/10.5194/gmd-14-2545-2021
© Author(s) 2021. 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-14-2545-2021
© Author(s) 2021. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Development and technical paper
| 
07 May 2021
Development and technical paper |
| 07 May 2021
| 07 May 2021
Assessment of numerical schemes for transient, finite-element ice flow models using ISSM v4.18
Thiago Dias dos Santos
CORRESPONDING AUTHOR
Department of Earth System Science, University of California, Irvine, CA, USA
Centro Polar e Climático, Universidade Federal do Rio Grande do Sul, Porto Alegre, RS, Brazil
Mathieu Morlighem
Department of Earth System Science, University of California, Irvine, CA, USA
Hélène Seroussi
Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, USA
Related authors
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.
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.
Stephen L. Cornford, Helene Seroussi, Xylar S. Asay-Davis, G. Hilmar Gudmundsson, Rob Arthern, Chris Borstad, Julia Christmann, Thiago Dias dos Santos, Johannes Feldmann, Daniel Goldberg, Matthew J. Hoffman, Angelika Humbert, Thomas Kleiner, Gunter Leguy, William H. Lipscomb, Nacho Merino, Gaël Durand, Mathieu Morlighem, David Pollard, Martin Rückamp, C. Rosie Williams, and Hongju Yu
The Cryosphere, 14, 2283–2301, https://doi.org/10.5194/tc-14-2283-2020, https://doi.org/10.5194/tc-14-2283-2020, 2020
Short summary
Short summary
We present the results of the third Marine Ice Sheet Intercomparison Project (MISMIP+). MISMIP+ is one in a series of exercises that test numerical models of ice sheet flow in simple situations. This particular exercise concentrates on the response of ice sheet models to the thinning of their floating ice shelves, which is of interest because numerical models are currently used to model the response to contemporary and near-future thinning in Antarctic ice shelves.
Thiago Dias dos Santos, Mathieu Morlighem, Hélène Seroussi, Philippe Remy Bernard Devloo, and Jefferson Cardia Simões
Geosci. Model Dev., 12, 215–232, https://doi.org/10.5194/gmd-12-215-2019, https://doi.org/10.5194/gmd-12-215-2019, 2019
Short summary
Short summary
The reduction of numerical errors in ice sheet modeling increases the results' accuracy reliability. We improve numerical accuracy by better capturing grounding line dynamics, while maintaining a low computational cost. We implement an adaptive mesh refinement (AMR) technique in the Ice Sheet System Model and compare AMR simulations with uniformly refined meshes. Our results show that the computational time with AMR is significantly shorter than for uniformly refined meshes for a given accuracy.
Ziad Rashed, Alexander A. Robel, and Hélène Seroussi
The Cryosphere, 19, 1775–1788, https://doi.org/10.5194/tc-19-1775-2025, https://doi.org/10.5194/tc-19-1775-2025, 2025
Short summary
Short summary
Sermeq Kujalleq, Greenland's largest glacier, has significantly retreated since the late 1990s in response to warming ocean temperatures. Using a large-ensemble approach, our simulations show that the retreat is mainly initiated by the arrival of warm water but sustained and accelerated by the glacier's position over deeper bed troughs and vigorous calving. We highlight the need for models of ice mélange to project glacier behavior under rapid calving regimes.
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.
Gong Cheng, Mansa Krishna, and Mathieu Morlighem
EGUsphere, https://doi.org/10.5194/egusphere-2025-1188, https://doi.org/10.5194/egusphere-2025-1188, 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, 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.
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.
Peter Van Katwyk, Baylor Fox-Kemper, Sophie Nowicki, Hélène Seroussi, and Karianne J. Bergen
EGUsphere, https://doi.org/10.5194/egusphere-2025-870, https://doi.org/10.5194/egusphere-2025-870, 2025
This preprint is open for discussion and under review for Geoscientific Model Development (GMD).
Short summary
Short summary
We present ISEFlow, a machine learning emulator that predicts how the melting of the Antarctic and Greenland ice sheets will contribute to sea level. ISEFlow is fast and accurate, allowing scientists to explore different climate scenarios with greater confidence. ISEFlow distinguishes between high and low emissions scenarios, helping us understand how today’s choices impact future sea levels. ISEFlow supports more reliable climate predictions and helps scientists study the future of ice sheets.
Benjamin Keith Galton-Fenzi, Richard Porter-Smith, Sue Cook, Eva Cougnon, David E. Gwyther, Wilma G. C. Huneke, Madelaine G. Rosevear, Xylar Asay-Davis, Fabio Boeira Dias, Michael S. Dinniman, David Holland, Kazuya Kusahara, Kaitlin A. Naughten, Keith W. Nicholls, Charles Pelletier, Ole Richter, Helene L. Seroussi, and Ralph Timmermann
EGUsphere, https://doi.org/10.5194/egusphere-2024-4047, https://doi.org/10.5194/egusphere-2024-4047, 2025
Short summary
Short summary
Melting beneath Antarctica’s floating ice shelves is key to future sea-level rise. We compare several different ocean simulations with satellite measurements, and provide the first multi-model average estimate of melting and refreezing driven by both ocean temperature and currents beneath ice shelves. The multi-model average can provide a useful tool for better understanding the role of ice shelf melting in present-day and future ice-sheet changes and informing coastal adaptation efforts.
James F. O'Neill, Tamsin L. Edwards, Daniel F. Martin, Courtney Shafer, Stephen L. Cornford, Hélène L. Seroussi, Sophie Nowicki, Mira Adhikari, and Lauren J. Gregoire
The Cryosphere, 19, 541–563, https://doi.org/10.5194/tc-19-541-2025, https://doi.org/10.5194/tc-19-541-2025, 2025
Short summary
Short summary
We use an ice sheet model to simulate the Antarctic contribution to sea level over the 21st century under a range of future climates and varying how sensitive the ice sheet is to different processes. We find that ocean temperatures increase and more snow falls on the ice sheet under stronger warming scenarios. When the ice sheet is sensitive to ocean warming, ocean melt-driven loss exceeds snowfall-driven gains, meaning that the sea level contribution is greater with more climate warming.
Luc Houriez, Eric Larour, Lambert Caron, Nicole-Jeanne Schlegel, Surendra Adhikari, Erik Ivins, Tyler Pelle, Hélène Seroussi, Eric Darve, and Martin Fischer
EGUsphere, https://doi.org/10.5194/egusphere-2024-4136, https://doi.org/10.5194/egusphere-2024-4136, 2025
Short summary
Short summary
We studied how interactions between the ice sheet and the Earth’s evolving surface affect the future of Thwaites Glacier in Antarctica. We find that small features in the bedrock play a major role in these interactions which can delay the glacier’s retreat by decades or even centuries. This can significantly reduce sea-level rise projections. Our work highlights resolution requirements for similar ice—earth models, and the importance of bedrock mapping efforts in Antarctica.
Vincent Verjans, Alexander A. Robel, Lizz Ultee, Helene Seroussi, Andrew F. Thompson, Lars Ackerman, Youngmin Choi, and Uta Krebs-Kanzow
EGUsphere, https://doi.org/10.5194/egusphere-2024-4067, https://doi.org/10.5194/egusphere-2024-4067, 2025
Short summary
Short summary
This study examines how random variations in climate may influence future ice loss from the Greenland Ice Sheet. We find that random climate variations are important for predicting future ice loss from the entire Greenland Ice Sheet over the next 20–30 years, but relatively unimportant after that period. Thus, uncertainty in sea level projections from the effect of climate variability on Greenland may play a role in coastal decision-making about sea level rise over the next few decades.
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.
Felicity Alice Holmes, Jamie Barnett, Henning Åkesson, Mathieu Morlighem, Johan Nilsson, Nina Kirchner, and Martin Jakobsson
EGUsphere, https://doi.org/10.5194/egusphere-2024-3839, https://doi.org/10.5194/egusphere-2024-3839, 2024
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 high emissions scenario. Very high levels of surface melt under a high emissions future leads to a sea-level rise contribution an order of magnitude higher than under a low emissions future.
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, and Javed Hassan
Earth Syst. Sci. Data Discuss., https://doi.org/10.5194/essd-2024-348, https://doi.org/10.5194/essd-2024-348, 2024
Revised manuscript accepted for ESSD
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 useful for GIS ice sheet modelling.
Jan De Rydt, Nicolas C. Jourdain, Yoshihiro Nakayama, Mathias van Caspel, Ralph Timmermann, Pierre Mathiot, Xylar S. Asay-Davis, Hélène Seroussi, Pierre Dutrieux, Ben Galton-Fenzi, David Holland, and Ronja Reese
Geosci. Model Dev., 17, 7105–7139, https://doi.org/10.5194/gmd-17-7105-2024, https://doi.org/10.5194/gmd-17-7105-2024, 2024
Short summary
Short summary
Global climate models do not reliably simulate sea-level change due to ice-sheet–ocean interactions. We propose a community modelling effort to conduct a series of well-defined experiments to compare models with observations and study how models respond to a range of perturbations in climate and ice-sheet geometry. The second Marine Ice Sheet–Ocean Model Intercomparison Project will continue to lay the groundwork for including ice-sheet–ocean interactions in global-scale IPCC-class models.
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.
Younghyun Koo, Gong Cheng, Mathieu Morlighem, and Maryam Rahnemoonfar
EGUsphere, https://doi.org/10.5194/egusphere-2024-1620, https://doi.org/10.5194/egusphere-2024-1620, 2024
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.
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.
Vincent Verjans, Alexander A. Robel, Helene Seroussi, Lizz Ultee, and Andrew F. Thompson
Geosci. Model Dev., 15, 8269–8293, https://doi.org/10.5194/gmd-15-8269-2022, https://doi.org/10.5194/gmd-15-8269-2022, 2022
Short summary
Short summary
We describe the development of the first large-scale ice sheet model that accounts for stochasticity in a range of processes. Stochasticity allows the impacts of inherently uncertain processes on ice sheets to be represented. This includes climatic uncertainty, as the climate is inherently chaotic. Furthermore, stochastic capabilities also encompass poorly constrained glaciological processes that display strong variability at fine spatiotemporal scales. We present the model and test experiments.
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.
Alexander A. Robel, Earle Wilson, and Helene Seroussi
The Cryosphere, 16, 451–469, https://doi.org/10.5194/tc-16-451-2022, https://doi.org/10.5194/tc-16-451-2022, 2022
Short summary
Short summary
Warm seawater may intrude as a thin layer below glaciers in contact with the ocean. Mathematical theory predicts that this intrusion may extend over distances of kilometers under realistic conditions. Computer models demonstrate that if this warm seawater causes melting of a glacier bottom, it can cause rates of glacier ice loss and sea level rise to be up to 2 times faster in response to potential future ocean warming.
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.
Jowan M. Barnes, Thiago Dias dos Santos, Daniel Goldberg, G. Hilmar Gudmundsson, Mathieu Morlighem, and Jan De Rydt
The Cryosphere, 15, 1975–2000, https://doi.org/10.5194/tc-15-1975-2021, https://doi.org/10.5194/tc-15-1975-2021, 2021
Short summary
Short summary
Some properties of ice flow models must be initialised using observed data before they can be used to produce reliable predictions of the future. Different models have different ways of doing this, and the process is generally seen as being specific to an individual model. We compare the methods used by three different models and show that they produce similar outputs. We also demonstrate that the outputs from one model can be used in other models without introducing large uncertainties.
William H. Lipscomb, Gunter R. Leguy, Nicolas C. Jourdain, Xylar Asay-Davis, Hélène Seroussi, and Sophie Nowicki
The Cryosphere, 15, 633–661, https://doi.org/10.5194/tc-15-633-2021, https://doi.org/10.5194/tc-15-633-2021, 2021
Short summary
Short summary
This paper describes Antarctic climate change experiments in which the Community Ice Sheet Model is forced with ocean warming predicted by global climate models. Generally, ice loss begins slowly, accelerates by 2100, and then continues unabated, with widespread retreat of the West Antarctic Ice Sheet. The mass loss by 2500 varies from about 150 to 1300 mm of equivalent sea level rise, based on the predicted ocean warming and assumptions about how this warming drives melting beneath ice shelves.
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.
Nicolas C. Jourdain, Xylar Asay-Davis, Tore Hattermann, Fiammetta Straneo, Hélène Seroussi, Christopher M. Little, and Sophie Nowicki
The Cryosphere, 14, 3111–3134, https://doi.org/10.5194/tc-14-3111-2020, https://doi.org/10.5194/tc-14-3111-2020, 2020
Short summary
Short summary
To predict the future Antarctic contribution to sea level rise, we need to use ice sheet models. The Ice Sheet Model Intercomparison Project for AR6 (ISMIP6) builds an ensemble of ice sheet projections constrained by atmosphere and ocean projections from the 6th Coupled Model Intercomparison Project (CMIP6). In this work, we present and assess a method to derive ice shelf basal melting in ISMIP6 from the CMIP6 ocean outputs, and we give examples of projected melt rates.
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.
Ronja Reese, Anders Levermann, Torsten Albrecht, Hélène Seroussi, and Ricarda Winkelmann
The Cryosphere, 14, 3097–3110, https://doi.org/10.5194/tc-14-3097-2020, https://doi.org/10.5194/tc-14-3097-2020, 2020
Short summary
Short summary
We compare 21st century projections of Antarctica's future sea-level contribution simulated with the Parallel Ice Sheet Model submitted to ISMIP6 with projections following the LARMIP-2 protocol based on the same model configuration. We find that (1) a preceding historic simulation increases mass loss by 5–50 % and that (2) the order of magnitude difference in the ice loss in our experiments following the two protocols can be explained by the translation of ocean forcing to sub-shelf melting.
Surendra Adhikari, Erik R. Ivins, Eric Larour, Lambert Caron, and Helene Seroussi
The Cryosphere, 14, 2819–2833, https://doi.org/10.5194/tc-14-2819-2020, https://doi.org/10.5194/tc-14-2819-2020, 2020
Short summary
Short summary
The mathematical formalism presented in this paper aims at simplifying computational strategies for tracking ice–ocean mass exchange in the Earth system. To this end, we define a set of generic, and quite simple, descriptions of evolving land, ocean and ice interfaces and present a unified method to compute the sea-level contribution of evolving ice sheets. The formalism can be applied to arbitrary geometries and at all timescales.
Sophie Nowicki, Heiko Goelzer, Hélène Seroussi, Anthony J. Payne, William H. Lipscomb, Ayako Abe-Ouchi, Cécile Agosta, Patrick Alexander, Xylar S. Asay-Davis, Alice Barthel, Thomas J. Bracegirdle, Richard Cullather, Denis Felikson, Xavier Fettweis, Jonathan M. Gregory, Tore Hattermann, Nicolas C. Jourdain, Peter Kuipers Munneke, Eric Larour, Christopher M. Little, Mathieu Morlighem, Isabel Nias, Andrew Shepherd, Erika Simon, Donald Slater, Robin S. Smith, Fiammetta Straneo, Luke D. Trusel, Michiel R. van den Broeke, and Roderik van de Wal
The Cryosphere, 14, 2331–2368, https://doi.org/10.5194/tc-14-2331-2020, https://doi.org/10.5194/tc-14-2331-2020, 2020
Short summary
Short summary
This paper describes the experimental protocol for ice sheet models taking part in the Ice Sheet Model Intercomparion Project for CMIP6 (ISMIP6) and presents an overview of the atmospheric and oceanic datasets to be used for the simulations. The ISMIP6 framework allows for exploring the uncertainty in 21st century sea level change from the Greenland and Antarctic ice sheets.
Stephen L. Cornford, Helene Seroussi, Xylar S. Asay-Davis, G. Hilmar Gudmundsson, Rob Arthern, Chris Borstad, Julia Christmann, Thiago Dias dos Santos, Johannes Feldmann, Daniel Goldberg, Matthew J. Hoffman, Angelika Humbert, Thomas Kleiner, Gunter Leguy, William H. Lipscomb, Nacho Merino, Gaël Durand, Mathieu Morlighem, David Pollard, Martin Rückamp, C. Rosie Williams, and Hongju Yu
The Cryosphere, 14, 2283–2301, https://doi.org/10.5194/tc-14-2283-2020, https://doi.org/10.5194/tc-14-2283-2020, 2020
Short summary
Short summary
We present the results of the third Marine Ice Sheet Intercomparison Project (MISMIP+). MISMIP+ is one in a series of exercises that test numerical models of ice sheet flow in simple situations. This particular exercise concentrates on the response of ice sheet models to the thinning of their floating ice shelves, which is of interest because numerical models are currently used to model the response to contemporary and near-future thinning in Antarctic ice shelves.
Donald A. Slater, Denis Felikson, Fiamma Straneo, Heiko Goelzer, Christopher M. Little, Mathieu Morlighem, Xavier Fettweis, and Sophie Nowicki
The Cryosphere, 14, 985–1008, https://doi.org/10.5194/tc-14-985-2020, https://doi.org/10.5194/tc-14-985-2020, 2020
Short summary
Short summary
Changes in the ocean around Greenland play an important role in determining how much the ice sheet will contribute to global sea level over the coming century. However, capturing these links in models is very challenging. This paper presents a strategy enabling an ensemble of ice sheet models to feel the effect of the ocean for the first time and should therefore result in a significant improvement in projections of the Greenland ice sheet's contribution to future sea level change.
Alice Barthel, Cécile Agosta, Christopher M. Little, Tore Hattermann, Nicolas C. Jourdain, Heiko Goelzer, Sophie Nowicki, Helene Seroussi, Fiammetta Straneo, and Thomas J. Bracegirdle
The Cryosphere, 14, 855–879, https://doi.org/10.5194/tc-14-855-2020, https://doi.org/10.5194/tc-14-855-2020, 2020
Short summary
Short summary
We compare existing coupled climate models to select a total of six models to provide forcing to the Greenland and Antarctic ice sheet simulations of the Ice Sheet Model Intercomparison Project (ISMIP6). We select models based on (i) their representation of current climate near Antarctica and Greenland relative to observations and (ii) their ability to sample a diversity of projected atmosphere and ocean changes over the 21st century.
Silje Smith-Johnsen, Basile de Fleurian, Nicole Schlegel, Helene Seroussi, and Kerim Nisancioglu
The Cryosphere, 14, 841–854, https://doi.org/10.5194/tc-14-841-2020, https://doi.org/10.5194/tc-14-841-2020, 2020
Short summary
Short summary
The Northeast Greenland Ice Stream (NEGIS) drains a large part of Greenland and displays fast flow far inland. However, the flow pattern is not well represented in ice sheet models. The fast flow has been explained by abnormally high geothermal heat flux. The heat melts the base of the ice sheet and the water produced may lubricate the bed and induce fast flow. By including high geothermal heat flux and a hydrology model, we successfully reproduce NEGIS flow pattern in an ice sheet model.
Anders Levermann, Ricarda Winkelmann, Torsten Albrecht, Heiko Goelzer, Nicholas R. Golledge, Ralf Greve, Philippe Huybrechts, Jim Jordan, Gunter Leguy, Daniel Martin, Mathieu Morlighem, Frank Pattyn, David Pollard, Aurelien Quiquet, Christian Rodehacke, Helene Seroussi, Johannes Sutter, Tong Zhang, Jonas Van Breedam, Reinhard Calov, Robert DeConto, Christophe Dumas, Julius Garbe, G. Hilmar Gudmundsson, Matthew J. Hoffman, Angelika Humbert, Thomas Kleiner, William H. Lipscomb, Malte Meinshausen, Esmond Ng, Sophie M. J. Nowicki, Mauro Perego, Stephen F. Price, Fuyuki Saito, Nicole-Jeanne Schlegel, Sainan Sun, and Roderik S. W. van de Wal
Earth Syst. Dynam., 11, 35–76, https://doi.org/10.5194/esd-11-35-2020, https://doi.org/10.5194/esd-11-35-2020, 2020
Short summary
Short summary
We provide an estimate of the future sea level contribution of Antarctica from basal ice shelf melting up to the year 2100. The full uncertainty range in the warming-related forcing of basal melt is estimated and applied to 16 state-of-the-art ice sheet models using a linear response theory approach. The sea level contribution we obtain is very likely below 61 cm under unmitigated climate change until 2100 (RCP8.5) and very likely below 40 cm if the Paris Climate Agreement is kept.
Hélène Seroussi, Sophie Nowicki, Erika Simon, Ayako Abe-Ouchi, Torsten Albrecht, Julien Brondex, Stephen Cornford, Christophe Dumas, Fabien Gillet-Chaulet, Heiko Goelzer, Nicholas R. Golledge, Jonathan M. Gregory, Ralf Greve, Matthew J. Hoffman, Angelika Humbert, Philippe Huybrechts, Thomas Kleiner, Eric Larour, Gunter Leguy, William H. Lipscomb, Daniel Lowry, Matthias Mengel, Mathieu Morlighem, Frank Pattyn, Anthony J. Payne, David Pollard, Stephen F. Price, Aurélien Quiquet, Thomas J. Reerink, Ronja Reese, Christian B. Rodehacke, Nicole-Jeanne Schlegel, Andrew Shepherd, Sainan Sun, Johannes Sutter, Jonas Van Breedam, Roderik S. W. van de Wal, Ricarda Winkelmann, and Tong Zhang
The Cryosphere, 13, 1441–1471, https://doi.org/10.5194/tc-13-1441-2019, https://doi.org/10.5194/tc-13-1441-2019, 2019
Short summary
Short summary
We compare a wide range of Antarctic ice sheet simulations with varying initialization techniques and model parameters to understand the role they play on the projected evolution of this ice sheet under simple scenarios. Results are improved compared to previous assessments and show that continued improvements in the representation of the floating ice around Antarctica are critical to reduce the uncertainty in the future ice sheet contribution to sea level rise.
Tyler Pelle, Mathieu Morlighem, and Johannes H. Bondzio
The Cryosphere, 13, 1043–1049, https://doi.org/10.5194/tc-13-1043-2019, https://doi.org/10.5194/tc-13-1043-2019, 2019
Short summary
Short summary
How ocean-induced melt under floating ice shelves will change as ocean currents evolve remains a big uncertainty in projections of sea level rise. In this study, we combine two of the most recently developed melt models to form PICOP, which overcomes the limitations of past models and produces accurate ice shelf melt rates. We find that our model is easy to set up and computationally efficient, providing researchers an important tool to improve the accuracy of their future glacial projections.
Joshua K. Cuzzone, Nicole-Jeanne Schlegel, Mathieu Morlighem, Eric Larour, Jason P. Briner, Helene Seroussi, and Lambert Caron
The Cryosphere, 13, 879–893, https://doi.org/10.5194/tc-13-879-2019, https://doi.org/10.5194/tc-13-879-2019, 2019
Short summary
Short summary
We present ice sheet modeling results of ice retreat over southwestern Greenland during the last 12 000 years, and we also test the impact that model horizontal resolution has on differences in the simulated spatial retreat and its associated rate. Results indicate that model resolution plays a minor role in simulated retreat in areas where bed topography is not complex but plays an important role in areas where bed topography is complex (such as fjords).
Mathieu Morlighem, Michael Wood, Hélène Seroussi, Youngmin Choi, and Eric Rignot
The Cryosphere, 13, 723–734, https://doi.org/10.5194/tc-13-723-2019, https://doi.org/10.5194/tc-13-723-2019, 2019
Short summary
Short summary
Many glaciers along the coast of Greenland have been retreating. It has been suggested that this retreat is triggered by the presence of warm water in the fjords, and surface melt at the top of the ice sheet is exacerbating this problem. Here, we quantify the vulnerability of northwestern Greenland to further warming using a numerical model. We find that in current conditions, this sector alone will contribute more than 1 cm to sea rise level by 2100, and up to 3 cm in the most extreme scenario.
Thiago Dias dos Santos, Mathieu Morlighem, Hélène Seroussi, Philippe Remy Bernard Devloo, and Jefferson Cardia Simões
Geosci. Model Dev., 12, 215–232, https://doi.org/10.5194/gmd-12-215-2019, https://doi.org/10.5194/gmd-12-215-2019, 2019
Short summary
Short summary
The reduction of numerical errors in ice sheet modeling increases the results' accuracy reliability. We improve numerical accuracy by better capturing grounding line dynamics, while maintaining a low computational cost. We implement an adaptive mesh refinement (AMR) technique in the Ice Sheet System Model and compare AMR simulations with uniformly refined meshes. Our results show that the computational time with AMR is significantly shorter than for uniformly refined meshes for a given accuracy.
Hongju Yu, Eric Rignot, Helene Seroussi, and Mathieu Morlighem
The Cryosphere, 12, 3861–3876, https://doi.org/10.5194/tc-12-3861-2018, https://doi.org/10.5194/tc-12-3861-2018, 2018
Short summary
Short summary
Thwaites Glacier, West Antarctica, has experienced rapid grounding line retreat and mass loss in the past decades. In this study, we simulate the evolution of Thwaites Glacier over the next century using different model configurations. Overall, we estimate a 5 mm contribution to global sea level rise from Thwaites Glacier in the next 30 years. However, a 300 % uncertainty is found over the next 100 years, ranging from 14 to 42 mm, depending on the model setup.
Youngmin Choi, Mathieu Morlighem, Michael Wood, and Johannes H. Bondzio
The Cryosphere, 12, 3735–3746, https://doi.org/10.5194/tc-12-3735-2018, https://doi.org/10.5194/tc-12-3735-2018, 2018
Short summary
Short summary
Calving is an important mechanism that controls the dynamics of Greenland outlet glaciers. We test and compare four calving laws and assess which calving law has better predictive abilities. Overall, the calving law based on von Mises stress is more satisfactory than other laws, but new parameterizations should be derived to better capture the detailed processes involved in calving.
Nicole-Jeanne Schlegel, Helene Seroussi, Michael P. Schodlok, Eric Y. Larour, Carmen Boening, Daniel Limonadi, Michael M. Watkins, Mathieu Morlighem, and Michiel R. van den Broeke
The Cryosphere, 12, 3511–3534, https://doi.org/10.5194/tc-12-3511-2018, https://doi.org/10.5194/tc-12-3511-2018, 2018
Short summary
Short summary
Using NASA supercomputers and a novel framework, in which Sandia National Laboratories' statistical software is embedded in the Jet Propulsion Laboratory's ice sheet model, we run a range of 100-year warming scenarios for Antarctica. We find that 1.2 m of sea level contribution is achievable, but not likely. Also, we find that bedrock topography beneath the ice drives potential for regional sea level contribution, highlighting the need for accurate bedrock mapping of the ice sheet interior.
Veronika Emetc, Paul Tregoning, Mathieu Morlighem, Chris Borstad, and Malcolm Sambridge
The Cryosphere, 12, 3187–3213, https://doi.org/10.5194/tc-12-3187-2018, https://doi.org/10.5194/tc-12-3187-2018, 2018
Short summary
Short summary
The paper includes a model that can be used to predict zones of fracture formation in both floating and grounded ice in Antarctica. We used observations and a statistics-based model to predict fractures in most ice shelves in Antarctica as an alternative to the damage-based approach. We can predict the location of observed fractures with an average success rate of 84% for grounded ice and 61% for floating ice and mean overestimation error of 26% and 20%, respectively.
Hélène Seroussi and Mathieu Morlighem
The Cryosphere, 12, 3085–3096, https://doi.org/10.5194/tc-12-3085-2018, https://doi.org/10.5194/tc-12-3085-2018, 2018
Chen Zhao, Rupert M. Gladstone, Roland C. Warner, Matt A. King, Thomas Zwinger, and Mathieu Morlighem
The Cryosphere, 12, 2637–2652, https://doi.org/10.5194/tc-12-2637-2018, https://doi.org/10.5194/tc-12-2637-2018, 2018
Short summary
Short summary
A combination of computer modelling and observational data were used to infer the resistance to ice flow at the bed of the Fleming Glacier on the Antarctic Peninsula. The model was also used to simulate the distribution of temperature within the ice, which governs the rate at which the ice can deform. This is especially important for glaciers like the Fleming Glacier, which has both regions of rapid deformation and regions of rapid sliding at the bed.
Chen Zhao, Rupert M. Gladstone, Roland C. Warner, Matt A. King, Thomas Zwinger, and Mathieu Morlighem
The Cryosphere, 12, 2653–2666, https://doi.org/10.5194/tc-12-2653-2018, https://doi.org/10.5194/tc-12-2653-2018, 2018
Short summary
Short summary
A combination of computer modelling and observational data were used to infer the resistance to ice flow at the bed of the Fleming Glacier on the Antarctic Peninsula in both 2008 and 2015. The comparison suggests the grounding line retreated by ~ 9 km from 2008 to 2015. The retreat may be enhanced by a positive feedback between friction, melting and sliding at the glacier bed.
Aleah Sommers, Harihar Rajaram, and Mathieu Morlighem
Geosci. Model Dev., 11, 2955–2974, https://doi.org/10.5194/gmd-11-2955-2018, https://doi.org/10.5194/gmd-11-2955-2018, 2018
Short summary
Short summary
Meltwater drainage beneath glaciers and ice sheets influences how fast they move and is complicated and constantly changing. Most models distinguish between
fastand
slowdrainage with different equations for each system. The SHAKTI model allows for the ice–water drainage arrangement to transition naturally between different types of flow. This model can be used to understand how drainage affects glacier speeds and the associated ice loss to further inform predictions of sea level rise.
Joshua K. Cuzzone, Mathieu Morlighem, Eric Larour, Nicole Schlegel, and Helene Seroussi
Geosci. Model Dev., 11, 1683–1694, https://doi.org/10.5194/gmd-11-1683-2018, https://doi.org/10.5194/gmd-11-1683-2018, 2018
Short summary
Short summary
This paper details the implementation of higher-order vertical finite elements in the Ice Sheet System Model (ISSM). When using higher-order vertical finite elements, fewer vertical layers are needed to accurately capture the thermal structure in an ice sheet versus a conventional linear vertical interpolation, therefore greatly improving model runtime speeds, particularly in higher-order stress balance ice sheet models. The implications for paleoclimate ice sheet simulations are discussed.
Konstanze Haubner, Jason E. Box, Nicole J. Schlegel, Eric Y. Larour, Mathieu Morlighem, Anne M. Solgaard, Kristian K. Kjeldsen, Signe H. Larsen, Eric Rignot, Todd K. Dupont, and Kurt H. Kjær
The Cryosphere, 12, 1511–1522, https://doi.org/10.5194/tc-12-1511-2018, https://doi.org/10.5194/tc-12-1511-2018, 2018
Short summary
Short summary
We investigate the effect of neglecting calving on Upernavik Isstrøm, West Greenland, between 1849 and 2012.
Our simulation is forced with observed terminus positions in discrete time steps and is responsive to the prescribed ice front changes.
Simulated frontal retreat is needed to obtain a realistic ice surface elevation and velocity evolution of Upernavik.
Using the prescribed terminus position change we gain insight to mass loss partitioning during different time periods.
Heiko Goelzer, Sophie Nowicki, Tamsin Edwards, Matthew Beckley, Ayako Abe-Ouchi, Andy Aschwanden, Reinhard Calov, Olivier Gagliardini, Fabien Gillet-Chaulet, Nicholas R. Golledge, Jonathan Gregory, Ralf Greve, Angelika Humbert, Philippe Huybrechts, Joseph H. Kennedy, Eric Larour, William H. Lipscomb, Sébastien Le clec'h, Victoria Lee, Mathieu Morlighem, Frank Pattyn, Antony J. Payne, Christian Rodehacke, Martin Rückamp, Fuyuki Saito, Nicole Schlegel, Helene Seroussi, Andrew Shepherd, Sainan Sun, Roderik van de Wal, and Florian A. Ziemen
The Cryosphere, 12, 1433–1460, https://doi.org/10.5194/tc-12-1433-2018, https://doi.org/10.5194/tc-12-1433-2018, 2018
Short summary
Short summary
We have compared a wide spectrum of different initialisation techniques used in the ice sheet modelling community to define the modelled present-day Greenland ice sheet state as a starting point for physically based future-sea-level-change projections. Compared to earlier community-wide comparisons, we find better agreement across different models, which implies overall improvement of our understanding of what is needed to produce such initial states.
Hafeez Jeofry, Neil Ross, Hugh F. J. Corr, Jilu Li, Mathieu Morlighem, Prasad Gogineni, and Martin J. Siegert
Earth Syst. Sci. Data, 10, 711–725, https://doi.org/10.5194/essd-10-711-2018, https://doi.org/10.5194/essd-10-711-2018, 2018
Short summary
Short summary
Accurately characterizing the complexities of the ice-sheet dynamic specifically close to the grounding line across the Weddell Sea (WS) sector in the ice-sheet models provides challenges to the scientific community. Our main objective is to comprehend these complexities, adding accuracy to the projection of future ice-sheet dynamics. Therefore, we have developed a new bed elevation digital elevation model across the WS sector, which will be of value to ice-sheet modelling experiments.
Felicity S. Graham, Mathieu Morlighem, Roland C. Warner, and Adam Treverrow
The Cryosphere, 12, 1047–1067, https://doi.org/10.5194/tc-12-1047-2018, https://doi.org/10.5194/tc-12-1047-2018, 2018
Short summary
Short summary
Ice sheet flow is anisotropic, depending on the nature of the stress applied. However, most large-scale ice sheet models rely on the Glen flow relation, which ignores anisotropic effects. We implement a flow relation (ESTAR) for anisotropic ice in a large-scale ice sheet model. In ice shelf simulations, the Glen flow relation overestimates velocities by up to 17 % compared with ESTAR. Our results have implications for ice sheet model simulations of paleo-ice extent and sea level rise prediction.
Eric Larour, Daniel Cheng, Gilberto Perez, Justin Quinn, Mathieu Morlighem, Bao Duong, Lan Nguyen, Kit Petrie, Silva Harounian, Daria Halkides, and Wayne Hayes
Geosci. Model Dev., 10, 4393–4403, https://doi.org/10.5194/gmd-10-4393-2017, https://doi.org/10.5194/gmd-10-4393-2017, 2017
Short summary
Short summary
This work presents a new way of carrying out simulations using the C++ based Ice Sheet System Model (ISSM) within a web page. This allows for a new generation of websites that can rely on the entire code of a climate model, without compromising or simplifying the physics implemented in such a model. We believe this approach will enable better education/outreach websites as well as improve access to complex climate models without compromising their integrity.
Hongju Yu, Eric Rignot, Mathieu Morlighem, and Helene Seroussi
The Cryosphere, 11, 1283–1296, https://doi.org/10.5194/tc-11-1283-2017, https://doi.org/10.5194/tc-11-1283-2017, 2017
Short summary
Short summary
We combine 2-D ice flow model with linear elastic fracture mechanics (LEFM) to model the calving behavior of Thwaites Glacier, West Antarctica. We find the combination of full-Stokes (FS) model and LEFM produces crevasses that are consistent with observations. We also find that calving is enhanced with pre-existing surface crevasses, shorter ice shelves or undercut at the ice shelf front. We conclude that the FS/LEFM combination is capable of constraining crevasse formation and iceberg calving.
Daniel Farinotti, Douglas J. Brinkerhoff, Garry K. C. Clarke, Johannes J. Fürst, Holger Frey, Prateek Gantayat, Fabien Gillet-Chaulet, Claire Girard, Matthias Huss, Paul W. Leclercq, Andreas Linsbauer, Horst Machguth, Carlos Martin, Fabien Maussion, Mathieu Morlighem, Cyrille Mosbeux, Ankur Pandit, Andrea Portmann, Antoine Rabatel, RAAJ Ramsankaran, Thomas J. Reerink, Olivier Sanchez, Peter A. Stentoft, Sangita Singh Kumari, Ward J. J. van Pelt, Brian Anderson, Toby Benham, Daniel Binder, Julian A. Dowdeswell, Andrea Fischer, Kay Helfricht, Stanislav Kutuzov, Ivan Lavrentiev, Robert McNabb, G. Hilmar Gudmundsson, Huilin Li, and Liss M. Andreassen
The Cryosphere, 11, 949–970, https://doi.org/10.5194/tc-11-949-2017, https://doi.org/10.5194/tc-11-949-2017, 2017
Short summary
Short summary
ITMIX – the Ice Thickness Models Intercomparison eXperiment – was the first coordinated performance assessment for models inferring glacier ice thickness from surface characteristics. Considering 17 different models and 21 different test cases, we show that although solutions of individual models can differ considerably, an ensemble average can yield uncertainties in the order of 10 ± 24 % the mean ice thickness. Ways forward for improving such estimates are sketched.
Henning Åkesson, Kerim H. Nisancioglu, Rianne H. Giesen, and Mathieu Morlighem
The Cryosphere, 11, 281–302, https://doi.org/10.5194/tc-11-281-2017, https://doi.org/10.5194/tc-11-281-2017, 2017
Short summary
Short summary
We present simulations of the history of Hardangerjøkulen ice cap in southern Norway using a dynamical ice sheet model. From mid-Holocene ice-free conditions 4000 years ago, Hardangerjøkulen grows nonlinearly in response to a linear climate forcing, reaching maximum extent during the Little Ice Age (~ 1750 AD). The ice cap exhibits spatially asymmetric growth and retreat and is highly sensitive to climate change. Our results call for reassessment of glacier reconstructions from proxy records.
Feras Habbal, Eric Larour, Mathieu Morlighem, Helene Seroussi, Christopher P. Borstad, and Eric Rignot
Geosci. Model Dev., 10, 155–168, https://doi.org/10.5194/gmd-10-155-2017, https://doi.org/10.5194/gmd-10-155-2017, 2017
Short summary
Short summary
This work presents the results from testing a suite of numerical solvers on a standard ice sheet benchmark test. We note the relevance of this test to practical simulations and identify the fastest solvers for the transient simulation. The highlighted solvers show significant speed-ups in relation to the default solver (~1.5–100 times faster) and enable a new capability for solving massive, high-resolution models that are critical for improving projections of ice sheets and sea-level change.
Sophie M. J. Nowicki, Anthony Payne, Eric Larour, Helene Seroussi, Heiko Goelzer, William Lipscomb, Jonathan Gregory, Ayako Abe-Ouchi, and Andrew Shepherd
Geosci. Model Dev., 9, 4521–4545, https://doi.org/10.5194/gmd-9-4521-2016, https://doi.org/10.5194/gmd-9-4521-2016, 2016
Short summary
Short summary
This paper describes an experimental protocol designed to quantify and understand the global sea level that arises due to past, present, and future changes in the Greenland and Antarctic ice sheets, along with investigating ice sheet–climate feedbacks. The Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) protocol includes targeted experiments, and a set of output diagnostic related to ice sheets, that are part of the 6th phase of the Coupled Model Intercomparison Project (CMIP6).
Eric Larour, Jean Utke, Anton Bovin, Mathieu Morlighem, and Gilberto Perez
Geosci. Model Dev., 9, 3907–3918, https://doi.org/10.5194/gmd-9-3907-2016, https://doi.org/10.5194/gmd-9-3907-2016, 2016
Short summary
Short summary
We present an approach to derive the adjoint state of the C++ coded Ice Sheet System Model. The approach enables data assimilation of observations to improve projections of polar ice sheet mass balance and contribution to sea-level rise. It is applicable to other Earth science frameworks relying on C++ and parallel computing, is non-intrusive, and enables computation of transient adjoints for any type of physics, hence providing insights into the sensitivities of any model to its inputs.
Janin Schaffer, Ralph Timmermann, Jan Erik Arndt, Steen Savstrup Kristensen, Christoph Mayer, Mathieu Morlighem, and Daniel Steinhage
Earth Syst. Sci. Data, 8, 543–557, https://doi.org/10.5194/essd-8-543-2016, https://doi.org/10.5194/essd-8-543-2016, 2016
Short summary
Short summary
The RTopo-2 data set provides consistent maps of global ocean bathymetry and ice surface topographies for Greenland and Antarctica at 30 arcsec grid spacing. We corrected data from earlier products in the areas of Petermann, Hagen Bræ, and Helheim glaciers, incorporated original data for the floating ice tongue of Nioghalvfjerdsfjorden Glacier, and applied corrections for the geometry of Getz, Abbot, and Fimbul ice shelf cavities. The data set is available from the PANGAEA database.
Xylar S. Asay-Davis, Stephen L. Cornford, Gaël Durand, Benjamin K. Galton-Fenzi, Rupert M. Gladstone, G. Hilmar Gudmundsson, Tore Hattermann, David M. Holland, Denise Holland, Paul R. Holland, Daniel F. Martin, Pierre Mathiot, Frank Pattyn, and Hélène Seroussi
Geosci. Model Dev., 9, 2471–2497, https://doi.org/10.5194/gmd-9-2471-2016, https://doi.org/10.5194/gmd-9-2471-2016, 2016
Short summary
Short summary
Coupled ice sheet–ocean models capable of simulating moving grounding lines are just becoming available. Such models have a broad range of potential applications in studying the dynamics of ice sheets and glaciers, including assessing their contributions to sea level change. Here we describe the idealized experiments that make up three interrelated Model Intercomparison Projects (MIPs) for marine ice sheet models and regional ocean circulation models incorporating ice shelf cavities.
Hongju Yu, Eric Rignot, Mathieu Morlighem, and Helene Seroussi
The Cryosphere Discuss., https://doi.org/10.5194/tc-2016-101, https://doi.org/10.5194/tc-2016-101, 2016
Revised manuscript not accepted
Short summary
Short summary
We performed a 2D Full-Stokes (FS) modeling study of grounding line dynamics and calving of Thwaites Glacier, West Antarctica. We compare FS with simplified models on grounding line migration and we combine FS with Linear Elastic Fracture Mechanics to simulate crevasse propagation. We find that only FS is able to provide reliable grounding line migration and to explain observed crevasse. We conclude that it may be essential to employ FS in the grounding line region for 2D simulations.
Pierre Rampal, Sylvain Bouillon, Einar Ólason, and Mathieu Morlighem
The Cryosphere, 10, 1055–1073, https://doi.org/10.5194/tc-10-1055-2016, https://doi.org/10.5194/tc-10-1055-2016, 2016
Short summary
Short summary
The Arctic sea ice cover has changed drastically over the last decades and undergone a shift in its dynamical regime, as seen by the increase of extreme fracturing events and the acceleration of sea ice drift. In this paper we present a new sea ice model, neXtSIM, that is capable of simulating both sea ice drift and deformation as observed from satellites, with similar spatial and temporal scaling properties. At the same time, the model reproduces sea ice area, extent, and volume correctly.
Johannes H. Bondzio, Hélène Seroussi, Mathieu Morlighem, Thomas Kleiner, Martin Rückamp, Angelika Humbert, and Eric Y. Larour
The Cryosphere, 10, 497–510, https://doi.org/10.5194/tc-10-497-2016, https://doi.org/10.5194/tc-10-497-2016, 2016
Short summary
Short summary
We implemented a level-set method in the ice sheet system model. This method allows us to dynamically evolve a calving front subject to user-defined calving rates. We apply the method to Jakobshavn Isbræ, West Greenland, and study its response to calving rate perturbations. We find its behaviour strongly dependent on the calving rate, which was to be expected. Both reduced basal drag and rheological shear margin weakening sustain the acceleration of this dynamic outlet glacier.
Laura A. Stevens, Fiamma Straneo, Sarah B. Das, Albert J. Plueddemann, Amy L. Kukulya, and Mathieu Morlighem
The Cryosphere, 10, 417–432, https://doi.org/10.5194/tc-10-417-2016, https://doi.org/10.5194/tc-10-417-2016, 2016
Short summary
Short summary
Here we pair detailed hydrographic measurements collected with an autonomous underwater vehicle as close as 150 m from the ice–ocean interface of the Saqqarliup sermia–Sarqardleq Fjord system, West Greenland, with modeled and observed subglacial discharge locations and magnitudes. We find evidence of two main types of subsurface glacially modified water localized in space that are consistent with runoff discharged at two locations along the grounding line.
K. Le Morzadec, L. Tarasov, M. Morlighem, and H. Seroussi
Geosci. Model Dev., 8, 3199–3213, https://doi.org/10.5194/gmd-8-3199-2015, https://doi.org/10.5194/gmd-8-3199-2015, 2015
Short summary
Short summary
A long-term challenge for any model of complex large-scale processes
is accounting for the impact of unresolved sub-grid (SG) processes.
We quantify the impact of SG mass-balance and ice fluxes on glacial
cycle ensemble results for North America. We find no easy solutions to
accurately capture these impacts. We show that SG process
representation and associated parametric uncertainties can have
significant impact on coarse resolution model results for glacial
cycle ice sheet evolution.
E. Larour, J. Utke, B. Csatho, A. Schenk, H. Seroussi, M. Morlighem, E. Rignot, N. Schlegel, and A. Khazendar
The Cryosphere, 8, 2335–2351, https://doi.org/10.5194/tc-8-2335-2014, https://doi.org/10.5194/tc-8-2335-2014, 2014
Short summary
Short summary
We present a temporal inversion of surface mass balance and basal friction for the Northeast Greenland Ice Sheet between 2003 and 2009, using the altimetry record from ICESat. The inversion relies on automatic differentiation of ISSM and demonstrates the feasibility of assimilating altimetry records into reconstructions of the Greenland Ice Sheet. The boundary conditions provide a snapshot of the state of the ice for this period and can be used for further process studies.
H. Seroussi, M. Morlighem, E. Larour, E. Rignot, and A. Khazendar
The Cryosphere, 8, 2075–2087, https://doi.org/10.5194/tc-8-2075-2014, https://doi.org/10.5194/tc-8-2075-2014, 2014
J. M. Lea, D. W. F. Mair, F. M. Nick, B. R. Rea, D. van As, M. Morlighem, P. W. Nienow, and A. Weidick
The Cryosphere, 8, 2031–2045, https://doi.org/10.5194/tc-8-2031-2014, https://doi.org/10.5194/tc-8-2031-2014, 2014
H. Seroussi, M. Morlighem, E. Rignot, J. Mouginot, E. Larour, M. Schodlok, and A. Khazendar
The Cryosphere, 8, 1699–1710, https://doi.org/10.5194/tc-8-1699-2014, https://doi.org/10.5194/tc-8-1699-2014, 2014
S. Adhikari, E. R. Ivins, E. Larour, H. Seroussi, M. Morlighem, and S. Nowicki
Solid Earth, 5, 569–584, https://doi.org/10.5194/se-5-569-2014, https://doi.org/10.5194/se-5-569-2014, 2014
Related subject area
Cryosphere
CMIP6 models overestimate sea ice melt, growth and conduction relative to ice mass balance buoy estimates
Coupling framework (1.0) for the Úa (2023b) ice sheet model and the FESOM-1.4 z-coordinate ocean model in an Antarctic domain
A gradient-boosted tree framework to model the ice thickness of the world's glaciers (IceBoost v1.1)
Towards deep-learning solutions for classification of automated snow height measurements (CleanSnow v1.0.2)
Quantitative sub-ice and marine tracing of Antarctic sediment provenance (TASP v1.0)
Tuning parameters of a sea ice model using machine learning
WRF-Chem simulations of snow nitrate and other physicochemical properties in northern China
Clustering simulated snow profiles to form avalanche forecast regions
SnowQM 1.0: a fast R package for bias-correcting spatial fields of snow water equivalent using quantile mapping
Simulation of snow albedo and solar irradiance profile with the Two-streAm Radiative TransfEr in Snow (TARTES) v2.0 model
Computationally efficient subglacial drainage modelling using Gaussian Process emulators: GlaDS-GP v1.0
Evaluation of MITgcm-based ocean reanalyses for the Southern Ocean
Improvements in the land surface configuration to better simulate seasonal snow cover in the European Alps with the CNRM-AROME (cycle 46) convection-permitting regional climate model
A three-stage model pipeline predicting regional avalanche danger in Switzerland (RAvaFcast v1.0.0): a decision-support tool for operational avalanche forecasting
A Flexible Snow Model (FSM 2.1.0) including a forest canopy
Anisotropic metric-based mesh adaptation for ice flow modelling in Firedrake
A global–land snow scheme (GLASS) v1.0 for the GFDL Earth System Model: formulation and evaluation at instrumented sites
Design and performance of ELSA v2.0: an isochronal model for ice-sheet layer tracing
Southern Ocean Ice Prediction System version 1.0 (SOIPS v1.0): description of the system and evaluation of synoptic-scale sea ice forecasts
Lagrangian tracking of sea ice in Community Ice CodE (CICE; version 5)
openAMUNDSEN v1.0: an open-source snow-hydrological model for mountain regions
OpenFOAM-avalanche 2312: depth-integrated models beyond dense-flow avalanches
Refactoring the elastic–viscous–plastic solver from the sea ice model CICE v6.5.1 for improved performance
A new 3D full-Stokes calving algorithm within Elmer/Ice (v9.0)
The Utrecht Finite Volume Ice-Sheet Model (UFEMISM version 2.0) – part 1: description and idealised experiments
A novel numerical implementation for the surface energy budget of melting snowpacks and glaciers
SnowPappus v1.0, a blowing-snow model for large-scale applications of the Crocus snow scheme
A stochastic parameterization of ice sheet surface mass balance for the Stochastic Ice-Sheet and Sea-Level System Model (StISSM v1.0)
Graphics-processing-unit-accelerated ice flow solver for unstructured meshes using the Shallow-Shelf Approximation (FastIceFlo v1.0.1)
A finite-element framework to explore the numerical solution of the coupled problem of heat conduction, water vapor diffusion, and settlement in dry snow (IvoriFEM v0.1.0)
Description and validation of the ice sheet model Nix v1.0
AvaFrame com1DFA (v1.3): a thickness-integrated computational avalanche module – theory, numerics, and testing
Universal differential equations for glacier ice flow modelling
A new model for supraglacial hydrology evolution and drainage for the Greenland Ice Sheet (SHED v1.0)
Modeling sensitivities of thermally and hydraulically driven ice stream surge cycling
A parallel implementation of the confined–unconfined aquifer system model for subglacial hydrology: design, verification, and performance analysis (CUAS-MPI v0.1.0)
Automatic snow type classification of snow micropenetrometer profiles with machine learning algorithms
An empirical model to calculate snow depth from daily snow water equivalent: SWE2HS 1.0
A wind-driven snow redistribution module for Alpine3D v3.3.0: adaptations designed for downscaling ice sheet surface mass balance
The CryoGrid community model (version 1.0) – a multi-physics toolbox for climate-driven simulations in the terrestrial cryosphere
Glacier Energy and Mass Balance (GEMB): a model of firn processes for cryosphere research
Sensitivity of NEMO4.0-SI3 model parameters on sea ice budgets in the Southern Ocean
Introducing CRYOWRF v1.0: multiscale atmospheric flow simulations with advanced snow cover modelling
SUHMO: an adaptive mesh refinement SUbglacial Hydrology MOdel v1.0
Improving snow albedo modeling in the E3SM land model (version 2.0) and assessing its impacts on snow and surface fluxes over the Tibetan Plateau
The Multiple Snow Data Assimilation System (MuSA v1.0)
The Stochastic Ice-Sheet and Sea-Level System Model v1.0 (StISSM v1.0)
Improved representation of the contemporary Greenland ice sheet firn layer by IMAU-FDM v1.2G
Modeling the small-scale deposition of snow onto structured Arctic sea ice during a MOSAiC storm using snowBedFoam 1.0.
Benchmarking the vertically integrated ice-sheet model IMAU-ICE (version 2.0)
Alex E. West and Edward W. Blockley
Geosci. Model Dev., 18, 3041–3064, https://doi.org/10.5194/gmd-18-3041-2025, https://doi.org/10.5194/gmd-18-3041-2025, 2025
Short summary
Short summary
This study uses ice mass balance buoys – temperature- and height-measuring devices frozen into sea ice – to find how well climate models simulate (1) melt and growth of Arctic sea ice and (2) conduction of heat through Arctic sea ice. This may help understand why models produce varying amounts of sea ice in the present day. We find that models tend to show more melt, growth or conduction for a given ice thickness than the buoys, although the difference is smaller for models with more physically realistic thermodynamics.
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.
Niccolò Maffezzoli, Eric Rignot, Carlo Barbante, Troels Petersen, and Sebastiano Vascon
Geosci. Model Dev., 18, 2545–2568, https://doi.org/10.5194/gmd-18-2545-2025, https://doi.org/10.5194/gmd-18-2545-2025, 2025
Short summary
Short summary
In this work we introduce IceBoost, a machine learning framework to model the ice thickness distribution of all the world's glaciers with greater accuracy than state-of-the-art methods. The model is trained on 3.7 million measurements globally available and provides skilful estimates across all regions. This advancement will help in better assessing future sea level changes and freshwater resources, with significance for both the scientific community and society at large.
Jan Svoboda, Marc Ruesch, David Liechti, Corinne Jones, Michele Volpi, Michael Zehnder, and Jürg Schweizer
Geosci. Model Dev., 18, 1829–1849, https://doi.org/10.5194/gmd-18-1829-2025, https://doi.org/10.5194/gmd-18-1829-2025, 2025
Short summary
Short summary
Accurately measuring snow height is key for modeling approaches in climate science, snow hydrology, and avalanche forecasting. Erroneous snow height measurements often occur when snow height is low or changes, for instance during snowfall in summer. We prepare a new benchmark dataset with annotated snow height data and demonstrate how to improve the measurement quality using modern deep-learning approaches. Our approach can be easily implemented in a data pipeline for snow modeling.
James W. Marschalek, Edward Gasson, Tina van de Flierdt, Claus-Dieter Hillenbrand, Martin J. Siegert, and Liam Holder
Geosci. Model Dev., 18, 1673–1708, https://doi.org/10.5194/gmd-18-1673-2025, https://doi.org/10.5194/gmd-18-1673-2025, 2025
Short summary
Short summary
Ice sheet models can help predict how Antarctica's ice sheets respond to environmental change, and such models benefit from comparison to geological data. Here, we use an ice sheet model output and other data to predict the erosion of debris and trace its transport to where it is deposited on the ocean floor. This allows the results of ice sheet modelling to be directly and quantitively compared to real-world data, helping to reduce uncertainty regarding Antarctic sea level contribution.
Anton Korosov, Yue Ying, and Einar Ólason
Geosci. Model Dev., 18, 885–904, https://doi.org/10.5194/gmd-18-885-2025, https://doi.org/10.5194/gmd-18-885-2025, 2025
Short summary
Short summary
We have developed a new method to improve the accuracy of sea ice models, which predict how ice moves and deforms due to wind and ocean currents. Traditional models use parameters that are often poorly defined. The new approach uses machine learning to fine-tune these parameters by comparing simulated ice drift with satellite data. The method identifies optimal settings for the model by analysing patterns in ice deformation. This results in more accurate simulations of sea ice drift forecasting.
Xia Wang, Tao Che, Xueyin Ruan, Shanna Yue, Jing Wang, Chun Zhao, and Lei Geng
Geosci. Model Dev., 18, 651–670, https://doi.org/10.5194/gmd-18-651-2025, https://doi.org/10.5194/gmd-18-651-2025, 2025
Short summary
Short summary
We employed the WRF-Chem model to parameterize atmospheric nitrate deposition in snow and evaluate its performance in simulating snow cover, snow depth, and concentrations of dust and nitrate using new observations from northern China. The results generally exhibit reasonable agreement with field observations in northern China, demonstrating the model's capability to simulate snow properties, including concentrations of reservoir species.
Simon Horton, Florian Herla, and Pascal Haegeli
Geosci. Model Dev., 18, 193–209, https://doi.org/10.5194/gmd-18-193-2025, https://doi.org/10.5194/gmd-18-193-2025, 2025
Short summary
Short summary
We present a method for avalanche forecasters to analyze patterns in snowpack model simulations. It uses fuzzy clustering to group small regions into larger forecast areas based on snow characteristics, locations, and temporal history. Tested in the Columbia Mountains in two winter seasons, it closely matched real forecast regions regions and identified major avalanche hazard patterns. This approach simplifies complex model outputs, helping forecasters make informed decisions.
Adrien Michel, Johannes Aschauer, Tobias Jonas, Stefanie Gubler, Sven Kotlarski, and Christoph Marty
Geosci. Model Dev., 17, 8969–8988, https://doi.org/10.5194/gmd-17-8969-2024, https://doi.org/10.5194/gmd-17-8969-2024, 2024
Short summary
Short summary
We present a method to correct snow cover maps (represented in terms of snow water equivalent) to match better-quality maps. The correction can then be extended backwards and forwards in time for periods when better-quality maps are not available. The method is fast and gives good results. It is then applied to obtain a climatology of the snow cover in Switzerland over the past 60 years at a resolution of 1 d and 1 km. This is the first time that such a dataset has been produced.
Ghislain Picard and Quentin Libois
Geosci. Model Dev., 17, 8927–8953, https://doi.org/10.5194/gmd-17-8927-2024, https://doi.org/10.5194/gmd-17-8927-2024, 2024
Short summary
Short summary
The Two-streAm Radiative TransfEr in Snow (TARTES) is a radiative transfer model to compute snow albedo in the solar domain and the profiles of light and energy absorption in a multi-layered snowpack whose physical properties are user defined. It uniquely considers snow grain shape flexibly, based on recent insights showing that snow does not behave as a collection of ice spheres but instead as a random medium. TARTES is user-friendly yet performs comparably to more complex models.
Tim Hill, Derek Bingham, Gwenn E. Flowers, and Matthew J. Hoffman
EGUsphere, https://doi.org/10.22541/essoar.172736254.41350153/v2, https://doi.org/10.22541/essoar.172736254.41350153/v2, 2024
Short summary
Short summary
Subglacial drainage models represent water flow beneath glaciers and ice sheets. Here, we train fast statistical models called Gaussian Process emulators to accelerate subglacial drainage modelling by ~1000 times. We use the fast emulator predictions to show that three of the model parameters are responsible for >90 % of the variance in model outputs. The fast GP emulators will enable future uncertainty quantification and calibration of these models.
Yoshihiro Nakayama, Alena Malyarenko, Hong Zhang, Ou Wang, Matthis Auger, Yafei Nie, Ian Fenty, Matthew Mazloff, Armin Köhl, and Dimitris Menemenlis
Geosci. Model Dev., 17, 8613–8638, https://doi.org/10.5194/gmd-17-8613-2024, https://doi.org/10.5194/gmd-17-8613-2024, 2024
Short summary
Short summary
Global- and basin-scale ocean reanalyses are becoming easily accessible. However, such ocean reanalyses are optimized for their entire model domains and their ability to simulate the Southern Ocean requires evaluation. We conduct intercomparison analyses of Massachusetts Institute of Technology General Circulation Model (MITgcm)-based ocean reanalyses. They generally perform well for the open ocean, but open-ocean temporal variability and Antarctic continental shelves require improvements.
Diego Monteiro, Cécile Caillaud, Matthieu Lafaysse, Adrien Napoly, Mathieu Fructus, Antoinette Alias, and Samuel Morin
Geosci. Model Dev., 17, 7645–7677, https://doi.org/10.5194/gmd-17-7645-2024, https://doi.org/10.5194/gmd-17-7645-2024, 2024
Short summary
Short summary
Modeling snow cover in climate and weather forecasting models is a challenge even for high-resolution models. Recent simulations with CNRM-AROME have shown difficulties when representing snow in the European Alps. Using remote sensing data and in situ observations, we evaluate modifications of the land surface configuration in order to improve it. We propose a new surface configuration, enabling a more realistic simulation of snow cover, relevant for climate and weather forecasting applications.
Alessandro Maissen, Frank Techel, and Michele Volpi
Geosci. Model Dev., 17, 7569–7593, https://doi.org/10.5194/gmd-17-7569-2024, https://doi.org/10.5194/gmd-17-7569-2024, 2024
Short summary
Short summary
By harnessing AI models, this work enables processing large amounts of data, including weather conditions, snowpack characteristics, and historical avalanche data, to predict human-like avalanche forecasts in Switzerland. Our proposed model can significantly assist avalanche forecasters in their decision-making process, thereby facilitating more efficient and accurate predictions crucial for ensuring safety in Switzerland's avalanche-prone regions.
Richard Essery, Giulia Mazzotti, Sarah Barr, Tobias Jonas, Tristan Quaife, and Nick Rutter
EGUsphere, https://doi.org/10.5194/egusphere-2024-2546, https://doi.org/10.5194/egusphere-2024-2546, 2024
Short summary
Short summary
How forests influence accumulation and melt of snow on the ground is of long-standing interest, but uncertainty remains in how best to model forest snow processes. We developed the Flexible Snow Model version 2 to quantify these uncertainties. In a first model demonstration, how unloading of intercepted snow from the forest canopy is represented is responsible for the largest uncertainty. Global mapping of forest distribution is also likely to be a large source of uncertainty in existing models.
Davor Dundovic, Joseph G. Wallwork, Stephan C. Kramer, Fabien Gillet-Chaulet, Regine Hock, and Matthew D. Piggott
EGUsphere, https://doi.org/10.5194/egusphere-2024-2649, https://doi.org/10.5194/egusphere-2024-2649, 2024
Short summary
Short summary
Accurate numerical studies of glaciers often require high-resolution simulations, which often prove too demanding even for modern computers. In this paper we develop a method that identifies whether different parts of a glacier require high or low resolution based on its physical features, such as its thickness and velocity. We show that by doing so we can achieve a more optimal simulation accuracy for the available computing resources compared to uniform resolution simulations.
Enrico Zorzetto, Sergey Malyshev, Paul Ginoux, and Elena Shevliakova
Geosci. Model Dev., 17, 7219–7244, https://doi.org/10.5194/gmd-17-7219-2024, https://doi.org/10.5194/gmd-17-7219-2024, 2024
Short summary
Short summary
We describe a new snow scheme developed for use in global climate models, which simulates the interactions of snowpack with vegetation, atmosphere, and soil. We test the new snow model over a set of sites where in situ observations are available. We find that when compared to a simpler snow model, this model improves predictions of seasonal snow and of soil temperature under the snowpack, important variables for simulating both the hydrological cycle and the global climate system.
Therese Rieckh, Andreas Born, Alexander Robinson, Robert Law, and Gerrit Gülle
Geosci. Model Dev., 17, 6987–7000, https://doi.org/10.5194/gmd-17-6987-2024, https://doi.org/10.5194/gmd-17-6987-2024, 2024
Short summary
Short summary
We present the open-source model ELSA, which simulates the internal age structure of large ice sheets. It creates layers of snow accumulation at fixed times during the simulation, which are used to model the internal stratification of the ice sheet. Together with reconstructed isochrones from radiostratigraphy data, ELSA can be used to assess ice sheet models and to improve their parameterization. ELSA can be used coupled to an ice sheet model or forced with its output.
Fu Zhao, Xi Liang, Zhongxiang Tian, Ming Li, Na Liu, and Chengyan Liu
Geosci. Model Dev., 17, 6867–6886, https://doi.org/10.5194/gmd-17-6867-2024, https://doi.org/10.5194/gmd-17-6867-2024, 2024
Short summary
Short summary
In this work, we introduce a newly developed Antarctic sea ice forecasting system, namely the Southern Ocean Ice Prediction System (SOIPS). The system is based on a regional sea ice‒ocean‒ice shelf coupled model and can assimilate sea ice concentration observations. By assessing the system's performance in sea ice forecasts, we find that the system can provide reliable Antarctic sea ice forecasts for the next 7 d and has the potential to guide ship navigation in the Antarctic sea ice zone.
Chenhui Ning, Shiming Xu, Yan Zhang, Xuantong Wang, Zhihao Fan, and Jiping Liu
Geosci. Model Dev., 17, 6847–6866, https://doi.org/10.5194/gmd-17-6847-2024, https://doi.org/10.5194/gmd-17-6847-2024, 2024
Short summary
Short summary
Sea ice models are mainly based on non-moving structured grids, which is different from buoy measurements that follow the ice drift. To facilitate Lagrangian analysis, we introduce online tracking of sea ice in Community Ice CodE (CICE). We validate the sea ice tracking with buoys and evaluate the sea ice deformation in high-resolution simulations, which show multi-fractal characteristics. The source code is openly available and can be used in various scientific and operational applications.
Ulrich Strasser, Michael Warscher, Erwin Rottler, and Florian Hanzer
Geosci. Model Dev., 17, 6775–6797, https://doi.org/10.5194/gmd-17-6775-2024, https://doi.org/10.5194/gmd-17-6775-2024, 2024
Short summary
Short summary
openAMUNDSEN is a fully distributed open-source snow-hydrological model for mountain catchments. It includes process representations of an empirical, semi-empirical, and physical nature. It uses temperature, precipitation, humidity, radiation, and wind speed as forcing data and is computationally efficient, of a modular nature, and easily extendible. The Python code is available on GitHub (https://github.com/openamundsen/openamundsen), including documentation (https://doc.openamundsen.org).
Matthias Rauter and Julia Kowalski
Geosci. Model Dev., 17, 6545–6569, https://doi.org/10.5194/gmd-17-6545-2024, https://doi.org/10.5194/gmd-17-6545-2024, 2024
Short summary
Short summary
Snow avalanches can form large powder clouds that substantially exceed the velocity and reach of the dense core. Only a few complex models exist to simulate this phenomenon, and the respective hazard is hard to predict. This work provides a novel flow model that focuses on simple relations while still encapsulating the significant behaviour. The model is applied to reconstruct two catastrophic powder snow avalanche events in Austria.
Till Andreas Soya Rasmussen, Jacob Poulsen, Mads Hvid Ribergaard, Ruchira Sasanka, Anthony P. Craig, Elizabeth C. Hunke, and Stefan Rethmeier
Geosci. Model Dev., 17, 6529–6544, https://doi.org/10.5194/gmd-17-6529-2024, https://doi.org/10.5194/gmd-17-6529-2024, 2024
Short summary
Short summary
Earth system models (ESMs) today strive for better quality based on improved resolutions and improved physics. A limiting factor is the supercomputers at hand and how best to utilize them. This study focuses on the refactorization of one part of a sea ice model (CICE), namely the dynamics. It shows that the performance can be significantly improved, which means that one can either run the same simulations much cheaper or advance the system according to what is needed.
Iain Wheel, Douglas I. Benn, Anna J. Crawford, Joe Todd, and Thomas Zwinger
Geosci. Model Dev., 17, 5759–5777, https://doi.org/10.5194/gmd-17-5759-2024, https://doi.org/10.5194/gmd-17-5759-2024, 2024
Short summary
Short summary
Calving, the detachment of large icebergs from glaciers, is one of the largest uncertainties in future sea level rise projections. This process is poorly understood, and there is an absence of detailed models capable of simulating calving. A new 3D calving model has been developed to better understand calving at glaciers where detailed modelling was previously limited. Importantly, the new model is very flexible. By allowing for unrestricted calving geometries, it can be applied at any location.
Constantijn J. Berends, Victor Azizi, Jorge Bernales, and Roderik S. W. van de Wal
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-5, https://doi.org/10.5194/gmd-2024-5, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
Ice-sheet models are computer programs that can simulate how the Greenland and Antarctic ice sheets will evolve in the future. The accuracy of these models depends on their resolution: how small the details are that the model can resolve. We have created a model with a variable resolution, which can resolve a lot of detail in areas where lots of changes happen in the ice, and less detail in areas where the ice does not move so much. This makes the model both accurate and fast.
Kévin Fourteau, Julien Brondex, Fanny Brun, and Marie Dumont
Geosci. Model Dev., 17, 1903–1929, https://doi.org/10.5194/gmd-17-1903-2024, https://doi.org/10.5194/gmd-17-1903-2024, 2024
Short summary
Short summary
In this paper, we provide a novel numerical implementation for solving the energy exchanges at the surface of snow and ice. By combining the strong points of previous models, our solution leads to more accurate and robust simulations of the energy exchanges, surface temperature, and melt while preserving a reasonable computation time.
Matthieu Baron, Ange Haddjeri, Matthieu Lafaysse, Louis Le Toumelin, Vincent Vionnet, and Mathieu Fructus
Geosci. Model Dev., 17, 1297–1326, https://doi.org/10.5194/gmd-17-1297-2024, https://doi.org/10.5194/gmd-17-1297-2024, 2024
Short summary
Short summary
Increasing the spatial resolution of numerical systems simulating snowpack evolution in mountain areas requires representing small-scale processes such as wind-induced snow transport. We present SnowPappus, a simple scheme coupled with the Crocus snow model to compute blowing-snow fluxes and redistribute snow among grid points at 250 m resolution. In terms of numerical cost, it is suitable for large-scale applications. We present point-scale evaluations of fluxes and snow transport occurrence.
Lizz Ultee, Alexander A. Robel, and Stefano Castruccio
Geosci. Model Dev., 17, 1041–1057, https://doi.org/10.5194/gmd-17-1041-2024, https://doi.org/10.5194/gmd-17-1041-2024, 2024
Short summary
Short summary
The surface mass balance (SMB) of an ice sheet describes the net gain or loss of mass from ice sheets (such as those in Greenland and Antarctica) through interaction with the atmosphere. We developed a statistical method to generate a wide range of SMB fields that reflect the best understanding of SMB processes. Efficiently sampling the variability of SMB will help us understand sources of uncertainty in ice sheet model projections.
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.
Julien Brondex, Kévin Fourteau, Marie Dumont, Pascal Hagenmuller, Neige Calonne, François Tuzet, and Henning Löwe
Geosci. Model Dev., 16, 7075–7106, https://doi.org/10.5194/gmd-16-7075-2023, https://doi.org/10.5194/gmd-16-7075-2023, 2023
Short summary
Short summary
Vapor diffusion is one of the main processes governing snowpack evolution, and it must be accounted for in models. Recent attempts to represent vapor diffusion in numerical models have faced several difficulties regarding computational cost and mass and energy conservation. Here, we develop our own finite-element software to explore numerical approaches and enable us to overcome these difficulties. We illustrate the capability of these approaches on established numerical benchmarks.
Daniel Moreno-Parada, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas
EGUsphere, https://doi.org/10.5194/egusphere-2023-2690, https://doi.org/10.5194/egusphere-2023-2690, 2023
Short summary
Short summary
We introduce Nix, an ice-sheet model designed for understanding how large masses of ice behave. Nix as a computer program that simulates the movement and temperature changes in ice sheets. Nix helps us study how ice sheets respond to changes in the atmosphere and ocean. We found that how fast ice melts under the shelves and how heat is exchanged, play a role in determining the future of ice sheets. Nix is a useful tool for learning more about how climate change affects polar ice sheets.
Matthias Tonnel, Anna Wirbel, Felix Oesterle, and Jan-Thomas Fischer
Geosci. Model Dev., 16, 7013–7035, https://doi.org/10.5194/gmd-16-7013-2023, https://doi.org/10.5194/gmd-16-7013-2023, 2023
Short summary
Short summary
Avaframe - the open avalanche framework - provides open-source tools to simulate and investigate snow avalanches. It is utilized for multiple purposes, the two main applications being hazard mapping and scientific research of snow processes. We present the theory, conversion to a computer model, and testing for one of the core modules used for simulations of a particular type of avalanche, the so-called dense-flow avalanches. Tests check and confirm the applicability of the utilized method.
Jordi Bolibar, Facundo Sapienza, Fabien Maussion, Redouane Lguensat, Bert Wouters, and Fernando Pérez
Geosci. Model Dev., 16, 6671–6687, https://doi.org/10.5194/gmd-16-6671-2023, https://doi.org/10.5194/gmd-16-6671-2023, 2023
Short summary
Short summary
We developed a new modelling framework combining numerical methods with machine learning. Using this approach, we focused on understanding how ice moves within glaciers, and we successfully learnt a prescribed law describing ice movement for 17 glaciers worldwide as a proof of concept. Our framework has the potential to discover important laws governing glacier processes, aiding our understanding of glacier physics and their contribution to water resources and sea-level rise.
Prateek Gantayat, Alison F. Banwell, Amber A. Leeson, James M. Lea, Dorthe Petersen, Noel Gourmelen, and Xavier Fettweis
Geosci. Model Dev., 16, 5803–5823, https://doi.org/10.5194/gmd-16-5803-2023, https://doi.org/10.5194/gmd-16-5803-2023, 2023
Short summary
Short summary
We developed a new supraglacial hydrology model for the Greenland Ice Sheet. This model simulates surface meltwater routing, meltwater drainage, supraglacial lake (SGL) overflow, and formation of lake ice. The model was able to reproduce 80 % of observed lake locations and provides a good match between the observed and modelled temporal evolution of SGLs.
Kevin Hank, Lev Tarasov, and Elisa Mantelli
Geosci. Model Dev., 16, 5627–5652, https://doi.org/10.5194/gmd-16-5627-2023, https://doi.org/10.5194/gmd-16-5627-2023, 2023
Short summary
Short summary
Physically meaningful modeling of geophysical system instabilities is numerically challenging, given the potential effects of purely numerical artifacts. Here we explore the sensitivity of ice stream surge activation to numerical and physical model aspects. We find that surge characteristics exhibit a resolution dependency but converge at higher horizontal grid resolutions and are significantly affected by the incorporation of bed thermal and sub-glacial hydrology models.
Yannic Fischler, Thomas Kleiner, Christian Bischof, Jeremie Schmiedel, Roiy Sayag, Raban Emunds, Lennart Frederik Oestreich, and Angelika Humbert
Geosci. Model Dev., 16, 5305–5322, https://doi.org/10.5194/gmd-16-5305-2023, https://doi.org/10.5194/gmd-16-5305-2023, 2023
Short summary
Short summary
Water underneath ice sheets affects the motion of glaciers. This study presents a newly developed code, CUAS-MPI, that simulates subglacial hydrology. It is designed for supercomputers and is hence a parallelized code. We measure the performance of this code for simulations of the entire Greenland Ice Sheet and find that the code works efficiently. Moreover, we validated the code to ensure the correctness of the solution. CUAS-MPI opens new possibilities for simulations of ice sheet hydrology.
Julia Kaltenborn, Amy R. Macfarlane, Viviane Clay, and Martin Schneebeli
Geosci. Model Dev., 16, 4521–4550, https://doi.org/10.5194/gmd-16-4521-2023, https://doi.org/10.5194/gmd-16-4521-2023, 2023
Short summary
Short summary
Snow layer segmentation and snow grain classification are essential diagnostic tasks for cryospheric applications. A SnowMicroPen (SMP) can be used to that end; however, the manual classification of its profiles becomes infeasible for large datasets. Here, we evaluate how well machine learning models automate this task. Of the 14 models trained on the MOSAiC SMP dataset, the long short-term memory model performed the best. The findings presented here facilitate and accelerate SMP data analysis.
Johannes Aschauer, Adrien Michel, Tobias Jonas, and Christoph Marty
Geosci. Model Dev., 16, 4063–4081, https://doi.org/10.5194/gmd-16-4063-2023, https://doi.org/10.5194/gmd-16-4063-2023, 2023
Short summary
Short summary
Snow water equivalent is the mass of water stored in a snowpack. Based on exponential settling functions, the empirical snow density model SWE2HS is presented to convert time series of daily snow water equivalent into snow depth. The model has been calibrated with data from Switzerland and validated with independent data from the European Alps. A reference implementation of SWE2HS is available as a Python package.
Eric Keenan, Nander Wever, Jan T. M. Lenaerts, and Brooke Medley
Geosci. Model Dev., 16, 3203–3219, https://doi.org/10.5194/gmd-16-3203-2023, https://doi.org/10.5194/gmd-16-3203-2023, 2023
Short summary
Short summary
Ice sheets gain mass via snowfall. However, snowfall is redistributed by the wind, resulting in accumulation differences of up to a factor of 5 over distances as short as 5 km. These differences complicate estimates of ice sheet contribution to sea level rise. For this reason, we have developed a new model for estimating wind-driven snow redistribution on ice sheets. We show that, over Pine Island Glacier in West Antarctica, the model improves estimates of snow accumulation variability.
Sebastian Westermann, Thomas Ingeman-Nielsen, Johanna Scheer, Kristoffer Aalstad, Juditha Aga, Nitin Chaudhary, Bernd Etzelmüller, Simon Filhol, Andreas Kääb, Cas Renette, Louise Steffensen Schmidt, Thomas Vikhamar Schuler, Robin B. Zweigel, Léo Martin, Sarah Morard, Matan Ben-Asher, Michael Angelopoulos, Julia Boike, Brian Groenke, Frederieke Miesner, Jan Nitzbon, Paul Overduin, Simone M. Stuenzi, and Moritz Langer
Geosci. Model Dev., 16, 2607–2647, https://doi.org/10.5194/gmd-16-2607-2023, https://doi.org/10.5194/gmd-16-2607-2023, 2023
Short summary
Short summary
The CryoGrid community model is a new tool for simulating ground temperatures and the water and ice balance in cold regions. It is a modular design, which makes it possible to test different schemes to simulate, for example, permafrost ground in an efficient way. The model contains tools to simulate frozen and unfrozen ground, snow, glaciers, and other massive ice bodies, as well as water bodies.
Alex S. Gardner, Nicole-Jeanne Schlegel, and Eric Larour
Geosci. Model Dev., 16, 2277–2302, https://doi.org/10.5194/gmd-16-2277-2023, https://doi.org/10.5194/gmd-16-2277-2023, 2023
Short summary
Short summary
This is the first description of the open-source Glacier Energy and Mass Balance (GEMB) model. GEMB models the ice sheet and glacier surface–atmospheric energy and mass exchange, as well as the firn state. The model is evaluated against the current state of the art and in situ observations and is shown to perform well.
Yafei Nie, Chengkun Li, Martin Vancoppenolle, Bin Cheng, Fabio Boeira Dias, Xianqing Lv, and Petteri Uotila
Geosci. Model Dev., 16, 1395–1425, https://doi.org/10.5194/gmd-16-1395-2023, https://doi.org/10.5194/gmd-16-1395-2023, 2023
Short summary
Short summary
State-of-the-art Earth system models simulate the observed sea ice extent relatively well, but this is often due to errors in the dynamic and other processes in the simulated sea ice changes cancelling each other out. We assessed the sensitivity of these processes simulated by the coupled ocean–sea ice model NEMO4.0-SI3 to 18 parameters. The performance of the model in simulating sea ice change processes was ultimately improved by adjusting the three identified key parameters.
Varun Sharma, Franziska Gerber, and Michael Lehning
Geosci. Model Dev., 16, 719–749, https://doi.org/10.5194/gmd-16-719-2023, https://doi.org/10.5194/gmd-16-719-2023, 2023
Short summary
Short summary
Most current generation climate and weather models have a relatively simplistic description of snow and snow–atmosphere interaction. One reason for this is the belief that including an advanced snow model would make the simulations too computationally demanding. In this study, we bring together two state-of-the-art models for atmosphere (WRF) and snow cover (SNOWPACK) and highlight both the feasibility and necessity of such coupled models to explore underexplored phenomena in the cryosphere.
Anne M. Felden, Daniel F. Martin, and Esmond G. Ng
Geosci. Model Dev., 16, 407–425, https://doi.org/10.5194/gmd-16-407-2023, https://doi.org/10.5194/gmd-16-407-2023, 2023
Short summary
Short summary
We present and validate a novel subglacial hydrology model, SUHMO, based on an adaptive mesh refinement framework. We propose the addition of a pseudo-diffusion to recover the wall melting in channels. Computational performance analysis demonstrates the efficiency of adaptive mesh refinement on large-scale hydrologic problems. The adaptive mesh refinement approach will eventually enable better ice bed boundary conditions for ice sheet simulations at a reasonable computational cost.
Dalei Hao, Gautam Bisht, Karl Rittger, Edward Bair, Cenlin He, Huilin Huang, Cheng Dang, Timbo Stillinger, Yu Gu, Hailong Wang, Yun Qian, and L. Ruby Leung
Geosci. Model Dev., 16, 75–94, https://doi.org/10.5194/gmd-16-75-2023, https://doi.org/10.5194/gmd-16-75-2023, 2023
Short summary
Short summary
Snow with the highest albedo of land surface plays a vital role in Earth’s surface energy budget and water cycle. This study accounts for the impacts of snow grain shape and mixing state of light-absorbing particles with snow on snow albedo in the E3SM land model. The findings advance our understanding of the role of snow grain shape and mixing state of LAP–snow in land surface processes and offer guidance for improving snow simulations and radiative forcing estimates in Earth system models.
Esteban Alonso-González, Kristoffer Aalstad, Mohamed Wassim Baba, Jesús Revuelto, Juan Ignacio López-Moreno, Joel Fiddes, Richard Essery, and Simon Gascoin
Geosci. Model Dev., 15, 9127–9155, https://doi.org/10.5194/gmd-15-9127-2022, https://doi.org/10.5194/gmd-15-9127-2022, 2022
Short summary
Short summary
Snow cover plays an important role in many processes, but its monitoring is a challenging task. The alternative is usually to simulate the snowpack, and to improve these simulations one of the most promising options is to fuse simulations with available observations (data assimilation). In this paper we present MuSA, a data assimilation tool which facilitates the implementation of snow monitoring initiatives, allowing the assimilation of a wide variety of remotely sensed snow cover information.
Vincent Verjans, Alexander A. Robel, Helene Seroussi, Lizz Ultee, and Andrew F. Thompson
Geosci. Model Dev., 15, 8269–8293, https://doi.org/10.5194/gmd-15-8269-2022, https://doi.org/10.5194/gmd-15-8269-2022, 2022
Short summary
Short summary
We describe the development of the first large-scale ice sheet model that accounts for stochasticity in a range of processes. Stochasticity allows the impacts of inherently uncertain processes on ice sheets to be represented. This includes climatic uncertainty, as the climate is inherently chaotic. Furthermore, stochastic capabilities also encompass poorly constrained glaciological processes that display strong variability at fine spatiotemporal scales. We present the model and test experiments.
Max Brils, Peter Kuipers Munneke, Willem Jan van de Berg, and Michiel van den Broeke
Geosci. Model Dev., 15, 7121–7138, https://doi.org/10.5194/gmd-15-7121-2022, https://doi.org/10.5194/gmd-15-7121-2022, 2022
Short summary
Short summary
Firn covers the Greenland ice sheet (GrIS) and can temporarily prevent mass loss. Here, we present the latest version of our firn model, IMAU-FDM, with an application to the GrIS. We improved the density of fallen snow, the firn densification rate and the firn's thermal conductivity. This leads to a higher air content and 10 m temperatures. Furthermore we investigate three case studies and find that the updated model shows greater variability and an increased sensitivity in surface elevation.
Océane Hames, Mahdi Jafari, David Nicholas Wagner, Ian Raphael, David Clemens-Sewall, Chris Polashenski, Matthew D. Shupe, Martin Schneebeli, and Michael Lehning
Geosci. Model Dev., 15, 6429–6449, https://doi.org/10.5194/gmd-15-6429-2022, https://doi.org/10.5194/gmd-15-6429-2022, 2022
Short summary
Short summary
This paper presents an Eulerian–Lagrangian snow transport model implemented in the fluid dynamics software OpenFOAM, which we call snowBedFoam 1.0. We apply this model to reproduce snow deposition on a piece of ridged Arctic sea ice, which was produced during the MOSAiC expedition through scan measurements. The model appears to successfully reproduce the enhanced snow accumulation and deposition patterns, although some quantitative uncertainties were shown.
Constantijn J. Berends, Heiko Goelzer, Thomas J. Reerink, Lennert B. Stap, and Roderik S. W. van de Wal
Geosci. Model Dev., 15, 5667–5688, https://doi.org/10.5194/gmd-15-5667-2022, https://doi.org/10.5194/gmd-15-5667-2022, 2022
Short summary
Short summary
The rate at which marine ice sheets such as the West Antarctic ice sheet will retreat in a warming climate and ocean is still uncertain. Numerical ice-sheet models, which solve the physical equations that describe the way glaciers and ice sheets deform and flow, have been substantially improved in recent years. Here we present the results of several years of work on IMAU-ICE, an ice-sheet model of intermediate complexity, which can be used to study ice sheets of both the past and the future.
Cited articles
Akin, J. and Tezduyar, T. E.: Calculation of the advective limit of the SUPG stabilization parameter for linear and higher-order elements, Comput. Method. Appl. M., 193, 1909–1922, https://doi.org/10.1016/j.cma.2003.12.050, 2004. a
Almeida, R. C. and Silva, R. S.: A stable Petrov–Galerkin method for convection-dominated problems, Comput. Method. Appl. M., 140, 291–304, https://doi.org/10.1016/S0045-7825(96)01108-5, 1997. a
Arnold, D. N., Brezzi, F., Cockburn, B., and Marini, L. D.: Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems, SIAM J. Numer. Anal., 39, 1749–1779, https://doi.org/10.1137/S0036142901384162, 2002. a
Aschwanden, A., Fahnestock, M. A., Truffer, M., Brinkerhoff, D. J., Hock, R., Khroulev, C., Mottram, R., and Khan, S. A.: Contribution of the Greenland Ice Sheet to sea level over the next millennium, Science Advances, 5, eaav9396, https://doi.org/10.1126/sciadv.aav9396, 2019. a
Babuška, I., Baumann, C., and Oden, J.: A discontinuous hp finite element method for diffusion problems: 1-D analysis, Comput. Math. Appl., 37, 103–122, https://doi.org/10.1016/S0898-1221(99)00117-0, 1999. a
Barnes, J. M., dos Santos, T. D., Goldberg, D., Gudmundsson, G. H., Morlighem, M., and De Rydt, J.: The transferability of adjoint inversion products between different ice flow models, The Cryosphere Discuss. [preprint], https://doi.org/10.5194/tc-2020-235, in review, 2020. a
Blank, H., Rudgyard, M., and Wathen, A.: Stabilised finite element methods for steady incompressible flow, Comput. Method. Appl. M., 174, 91–105, https://doi.org/10.1016/S0045-7825(98)00279-5, 1999. a
Blatter, H.: Velocity and stress-fields in grounded glaciers: A simple algorithm for including deviatoric stress gradients, J. Glaciol., 41, 333–344, 1995. a
Bochev, P. B., Gunzburger, M. D., and Shadid, J. N.: Stability of the SUPG finite element method for transient advection–diffusion problems, Comput. Method. Appl. M., 193, 2301–2323, https://doi.org/10.1016/j.cma.2004.01.026, 2004. a
Boris, J. P. and Book, D. L.: Flux-corrected transport. I. SHASTA, a fluid transport algorithm that works, J. Comput. Phys., 11, 38–69, https://doi.org/10.1016/0021-9991(73)90147-2, 1973. a
Brezzi, F., Manzini, G., Marini, D., Pietra, P., and Russo, A.: Discontinuous Galerkin approximations for elliptic problems, Numer. Meth. Part. D. E., 16, 365–378, https://doi.org/10.1002/1098-2426(200007)16:4<365::AID-NUM2>3.0.CO;2-Y, 2000. a
Brezzi, F., Marini, L. D., and Süli, E.: Discontinuous Galerkin methods for first-order hyperbolic problems, Math. Mod. Meth. Appl. S., 14, 1893–1903, https://doi.org/10.1142/S0218202504003866, 2004. a, b
Burman, E.: Consistent SUPG-method for transient transport problems: Stability and convergence, Comput. Method. Appl. M., 199, 1114–1123, https://doi.org/10.1016/j.cma.2009.11.023, 2010. a
Calle, J. L., Devloo, P. R., and Gomes, S. M.: Stabilized discontinuous Galerkin method for hyperbolic equations, Comput. Method. Appl. M., 194, 1861–1874, https://doi.org/10.1016/j.cma.2004.06.036, 2005. a
Christie, I., Griffiths, D. F., Mitchell, A. R., and Zienkiewicz, O. C.: Finite element methods for second order differential equations with significant first derivatives, Int. J. Numer. Meth. Eng., 10, 1389–1396, https://doi.org/10.1002/nme.1620100617, 1976. a
Cockburn, B.: An introduction to the Discontinuous Galerkin method for
convection-dominated problems, Springer, Berlin, Heidelberg, Germany, 150–268, https://doi.org/10.1007/BFb0096353, 1998. a, b, c, d
Cockburn, B.: Discontinuous Galerkin Methods for Convection-Dominated
Problems, Springer, Berlin, Heidelberg, Germany, 69–224, https://doi.org/10.1007/978-3-662-03882-6_2, 1999. a
Cockburn, B.: Discontinuous Galerkin methods, ZAMM-Z. Angew. Math. Me., 83, 731–754, https://doi.org/10.1002/zamm.200310088, 2003. a, b, c
Cockburn, B. and Shu, C.-W.: The Runge–Kutta local projection P1-discontinuous-Galerkin finite element method for scalar conservation laws, ESAIM: M2AN, 25, 337–361, https://doi.org/10.1051/m2an/1991250303371, 1991. a
Codina, R.: Comparison of some finite element methods for solving the diffusion-convection-reaction equation, Comput. Method. Appl. M., 156, 185–210, https://doi.org/10.1016/S0045-7825(97)00206-5, 1998. a, b
Codina, R.: On stabilized finite element methods for linear systems of convection–diffusion-reaction equations, Comput. Method. Appl. M., 188, 61–82, https://doi.org/10.1016/S0045-7825(00)00177-8, 2000. a
Cornford, S. L., Martin, D. F., Payne, A. J., Ng, E. G., Le Brocq, A. M., Gladstone, R. M., Edwards, T. L., Shannon, S. R., Agosta, C., van den Broeke, M. R., Hellmer, H. H., Krinner, G., Ligtenberg, S. R. M., Timmermann, R., and Vaughan, D. G.: Century-scale simulations of the response of the West Antarctic Ice Sheet to a warming climate, The Cryosphere, 9, 1579–1600, https://doi.org/10.5194/tc-9-1579-2015, 2015. a
Cornford, S. L., Martin, D. F., Lee, V., Payne, A. J., and Ng, E. G.: Adaptive mesh refinement versus subgrid friction interpolation in simulations of Antarctic ice dynamics, Ann. Glaciol., 57, 1–9, https://doi.org/10.1017/aog.2016.13, 2016. a
de Vahl Davis, G. and Mallinson, G.: An evaluation of upwind and central difference approximations by a study of recirculating flow, Comput. Fluids, 4, 29–43, https://doi.org/10.1016/0045-7930(76)90010-4, 1976. a
DeConto, R. M. and Pollard, D.: Contribution of Antarctica to past and future sea-level rise, Nature, 531, 591–597, https://doi.org/10.1038/nature17145, 2016. a
Devloo, P., Forti, T., and Gomes, S.: A combined continuous-discontinuous finite element method for convection-diffusion problems, Lat. Am. J. Solids Stru., 4, 229–246, 2007. a
Diddens, C.: Detailed finite element method modeling of evaporating multi-component droplets, J. Comput. Phys., 340, 670–687, https://doi.org/10.1016/j.jcp.2017.03.049, 2017. a
Docquier, D., Pollard, D., and Pattyn, F.: Thwaites Glacier grounding-line retreat: influence of width and buttressing parameterizations, J. Glaciol., 60, 305–313, https://doi.org/10.3189/2014JoG13J117, 2014. a
Donea, J.: A Taylor–Galerkin method for convective transport problems, Int. J. Numer. Meth. Eng., 20, 101–119, https://doi.org/10.1002/nme.1620200108, 1984a. a
Donea, J.: Recent advances in computational methods for steady and transient transport problems, Nuclear Engineering and Design, 80, 141–162, https://doi.org/10.1016/0029-5493(84)90163-8, 1984b. a
Donea, J., Giuliani, S., Laval, H., and Quartapelle, L.: Time-accurate solution of advection-diffusion problems by finite elements, Comput. Method. Appl. M., 45, 123–145, https://doi.org/10.1016/0045-7825(84)90153-1, 1984. a
Dupont, T. K. and Alley, R. B.: Assessment of the importance of ice-shelf buttressing to ice-sheet flow, Geophys. Res. Lett., 32, L04503, https://doi.org/10.1029/2004GL022024, 2005. a
Durand, G., Gagliardini, O., Zwinger, T., Le Meur, E., and Hindmarsh, R. C. A.: Full Stokes modeling of marine ice sheets: influence of the grid size, Ann. Glaciol., 50, 109–114, https://doi.org/10.3189/172756409789624283, 2009. a
Durand, G., Gagliardini, O., Favier, L., Zwinger, T., and Le Meur, E.: Impact of bedrock description on modeling ice sheet dynamics, Geophys. Res. Lett., 38, 1–6, https://doi.org/10.1029/2011GL048892, 2011. a
Favier, L., Durand, G., Cornford, S. L., Gudmundsson, G. H., Gagliardini, O., Gillet-Chaulet, F., Zwinger, T., Payne, A. J., and Le Brocq, A. M.: Retreat of Pine Island Glacier controlled by marine ice-sheet instability, Nat. Clim. Change, 4, 117–121, https://doi.org/10.1038/nclimate2094, 2014. a
Feldmann, J., Albrecht, T., Khroulev, C., Pattyn, F., and Levermann, A.: Resolution-dependent performance of grounding line motion in a shallow model compared with a full-Stokes model according to the MISMIP3d intercomparison, J. Glaciol., 60, 353–360, https://doi.org/10.3189/2014JoG13J093, 2014. a, b, c, d, e
Franca, L. P., Hauke, G., and Masud, A.: Revisiting stabilized finite element methods for the advective–diffusive equation, Comput. Method. Appl. M., 195, 1560–1572, https://doi.org/10.1016/j.cma.2005.05.028, 2006. a, b
Gagliardini, O., Zwinger, T., Gillet-Chaulet, F., Durand, G., Favier, L., de Fleurian, B., Greve, R., Malinen, M., Martín, C., Råback, P., Ruokolainen, J., Sacchettini, M., Schäfer, M., Seddik, H., and Thies, J.: Capabilities and performance of Elmer/Ice, a new-generation ice sheet model, Geosci. Model Dev., 6, 1299–1318, https://doi.org/10.5194/gmd-6-1299-2013, 2013. a, b, c
Goelzer, H., Nowicki, S., Payne, A., Larour, E., Seroussi, H., Lipscomb, W. H., Gregory, J., Abe-Ouchi, A., Shepherd, A., Simon, E., Agosta, C., Alexander, P., Aschwanden, A., Barthel, A., Calov, R., Chambers, C., Choi, Y., Cuzzone, J., Dumas, C., Edwards, T., Felikson, D., Fettweis, X., Golledge, N. R., Greve, R., Humbert, A., Huybrechts, P., Le clec'h, S., Lee, V., Leguy, G., Little, C., Lowry, D. P., Morlighem, M., Nias, I., Quiquet, A., Rückamp, M., Schlegel, N.-J., Slater, D. A., Smith, R. S., Straneo, F., Tarasov, L., van de Wal, R., and van den Broeke, M.: The future sea-level contribution of the Greenland ice sheet: a multi-model ensemble study of ISMIP6, The Cryosphere, 14, 3071–3096, https://doi.org/10.5194/tc-14-3071-2020, 2020. a
Goldberg, D., Holland, D. M., and Schoof, C.: Grounding line movement and ice shelf buttressing in marine ice sheets, J. Geophys. Res.-Earth, 114, 1–23, https://doi.org/10.1029/2008JF001227, 2009. a
Gresho, P. M. and Lee, R. L.: Don't suppress the wiggles – they're telling
you something, Tech. rep., United States, uCRL–82979, available at:
http://inis.iaea.org/search/search.aspx?orig_q=RN:11526189 (last access: 25 November 2020),
1979. a
Greve, R. and Blatter, H.: Dynamics of Ice Sheets and Glaciers, Advances in
Geophysical and Environmental Mechanics and Mathematics, 1st Edn., Springer,
Berlin, Heidelberg, Germany, 2009. a
Griffiths, D. and Lorenz, J.: An analysis of the petrov – galerkin finite element method, Comput. Method. Appl. M., 14, 39–64, https://doi.org/10.1016/0045-7825(78)90012-9, 1978. a
Gudmundsson, G. H.: GHilmarG/UaSource: Ua2019b (Version v2019b), Zenodo [code], https://doi.org/10.5281/zenodo.3706623, 2020. a, b, c, d
Gudmundsson, G. H., Krug, J., Durand, G., Favier, L., and Gagliardini, O.: The stability of grounding lines on retrograde slopes, The Cryosphere, 6, 1497–1505, https://doi.org/10.5194/tc-6-1497-2012, 2012. 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
Hansen, K. B., Arzani, A., and Shadden, S. C.: Finite element modeling of near-wall mass transport in cardiovascular flows, Int. J. Numer. Meth. Bio., 35, e3148, https://doi.org/10.1002/cnm.3148, e3148 cnm.3148, 2019. a
Heinrich, J. C., Huyakorn, P. S., Zienkiewicz, O. C., and Mitchell, A. R.: An `upwind' finite element scheme for two-dimensional convective transport equation, Int. J. Numer. Meth. Eng., 11, 131–143, https://doi.org/10.1002/nme.1620110113, 1977. a
Jameson, A.: Analysis and design of numerical schemes for gas dynamics, 1: Artificial diffusion, upwind biasing, limiters and their effect on accuracy and multigrid convergence, Int. J. Comput. Fluid D., 4, 171–218, https://doi.org/10.1080/10618569508904524, 1995. a
John, V., Knobloch, P., and Novo, J.: Finite elements for scalar convection-dominated equations and incompressible flow problems: a never ending story?, Comput. Visual. Sci., 19, 47–63, https://doi.org/10.1007/s00791-018-0290-5, 2018. a, b
Joughin, I., Smith, B. E., and Medley, B.: Marine Ice Sheet Collapse Potentially Under Way for the Thwaites Glacier Basin, West Antarctica, Science, 344, 735–738, https://doi.org/10.1126/science.1249055, 2014. a
Knopp, T., Lube, G., and Rapin, G.: Stabilized finite element methods with shock capturing for advection–diffusion problems, Comput. Method. Appl. M., 191, 2997–3013, https://doi.org/10.1016/S0045-7825(02)00222-0, 2002. a
Kuzmin, D. and Turek, S.: Flux Correction Tools for Finite Elements, J. Comput. Phys., 175, 525–558, https://doi.org/10.1006/jcph.2001.6955, 2002. a, b, c
Kuzmin, D., Möller, M., and Turek, S.: Multidimensional FEM-FCT schemes for arbitrary time stepping, Int. J. Numer. Meth. Fl., 42, 265–295, https://doi.org/10.1002/fld.493, 2003. a, b
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, 1–20, https://doi.org/10.1029/2011JF002140, 2012. a
Larour, E., Morlighem, M., and Seroussi, H.: Ice-Sheet and Sea-Level System
Model, svn repository, available at: https://issm.ess.uci.edu/svn/issm/issm/trunk (last access: 20 November 2020), 2020. a
LeVeque, R. J.: Numerical Methods for Conservation Laws, Lectures in mathematics ETH Zürich, 2nd edn., Birkhäuser, Basel, 1992. a
Löhner, R., Morgan, K., Peraire, J., and Vahdati, M.: Finite element flux-corrected transport (FEM–FCT) for the euler and Navier–Stokes equations, Int. J. Numer. Meth. Fl., 7, 1093–1109, https://doi.org/10.1002/fld.1650071007, 1987. a
MacAyeal, D.: 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, https://doi.org/10.1029/JB094iB04p04071, 1989. a, b, c, d
Martin, D. F., Cornford, S. L., and Payne, A. J.: Millennial-Scale Vulnerability of the Antarctic Ice Sheet to Regional Ice Shelf Collapse, Geophys. Res. Lett., 46, 1467–1475, https://doi.org/10.1029/2018GL081229, 2019. a
Möller, M., Kuzmin, D., and Turek, S.: Implicit Flux-Corrected Transport
Algorithm for Finite Element Simulation of the Compressible Euler Equations,
Springer, Berlin, Heidelberg, Germany, 325–354, https://doi.org/10.1007/978-3-642-18560-1_20, 2004. a
Morlighem, M.: MEaSUREs BedMachine Antarctica, Version 1, NASA National Snow and Ice Data Center Distributed Active Archive Center,
https://doi.org/10.5067/C2GFER6PTOS4, 2019. 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.-Earth, 118, 1746–1753, 2013. 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., Broeke, M. R. v. d.,
Ommen, T. D. v., 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
Mouginot, J., Rignot, E., and Scheuchl, B.: Continent-Wide, Interferometric SAR Phase, Mapping of Antarctic Ice Velocity, Geophys. Res. Lett., 46, 9710–9718, https://doi.org/10.1029/2019GL083826, 2019a. a
Mouginot, J., Rignot, E., and Scheuchl, B.: MEaSUREs Phase-Based Antarctica Ice Velocity Map, Version 1, NASA National Snow and Ice Data Center Distributed Active Archive Center, https://doi.org/10.5067/PZ3NJ5RXRH10, 2019b. a
Nakayama, Y., Manucharyan, G., Zhang, H., Dutrieux, P., Torres, H. S., Klein, P., Seroussi, H., Schodlok, M., Rignot, E., and Menemenlis, D.: Pathways of ocean heat towards Pine Island and Thwaites grounding lines, Sci. Rep.-UK, 9, 16649, https://doi.org/10.1038/s41598-019-53190-6, 2019. a
Ngo, A., Bastian, P., and Ippisch, O.: Numerical solution of steady-state groundwater flow and solute transport problems: Discontinuous Galerkin based methods compared to the Streamline Diffusion approach, Comput. Method. Appl. M., 294, 331–358, https://doi.org/10.1016/j.cma.2015.06.008, 2015. a
Pattyn, F., Schoof, C., Perichon, L., Hindmarsh, R. C. A., Bueler, E., de Fleurian, B., Durand, G., Gagliardini, O., Gladstone, R., Goldberg, D., Gudmundsson, G. H., Huybrechts, P., Lee, V., Nick, F. M., Payne, A. J., Pollard, D., Rybak, O., Saito, F., and Vieli, A.: Results of the Marine Ice Sheet Model Intercomparison Project, MISMIP, The Cryosphere, 6, 573–588, https://doi.org/10.5194/tc-6-573-2012, 2012. a
Pattyn, F., Perichon, L., Durand, G., Favier, L., Gagliardini, O., Hindmarsh, R. C. A., Zwinger, T., Albrecht, T., Cornford, S., Docquier, D., Fürst, J. J., Goldberg, D., Gudmundsson, G. H., Humbert, A., Hütten, M., Huybrechts, P., Jouvet, G., Kleiner, T., Larour, E., Martin, D., Morlighem, M., Payne, A. J., Pollard, D., Rückamp, M., Rybak, O., Seroussi, H., Thoma, M., and Wilkens, N.: Grounding-line migration in plan-view marine ice-sheet models: results of the ice2sea MISMIP3d intercomparison, J. Glaciol., 59, 410–422, https://doi.org/10.3189/2013JoG12J129, 2013. a, b, c, d
Ritz, C., Edwards, T. L., Durand, G., Payne, A. J., Peyaud, V., and Hindmarsh, R. C. A.: Potential sea-level rise from Antarctic Ice-Sheet instability constrained by observations, Nature, 528, 115–118, https://doi.org/10.1038/nature16147, https://doi.org/10.1038/nature16147, 2015. a
Schoof, C.: Marine ice-sheet dynamics. Part 1. The case of rapid sliding, J. Fluid Mech., 573, 27–55, https://doi.org/10.1017/S0022112006003570, 2007a. a, b
Schoof, C.: Ice sheet grounding line dynamics: Steady states, stability, and hysteresis, J. Geophys. Res.-Earth, 112, 1–19, https://doi.org/10.1029/2006JF000664, 2007b. a, b, c
Selmin, V.: The node-centred finite volume approach: Bridge between finite differences and finite elements, Comput. Method. Appl. M., 102, 107–138, https://doi.org/10.1016/0045-7825(93)90143-L, 1993. a, b
Seroussi, H., Morlighem, M., Rignot, E., Mouginot, J., Larour, E., Schodlok, M., and Khazendar, A.: Sensitivity of the dynamics of Pine Island Glacier, West Antarctica, to climate forcing for the next 50 years, The Cryosphere, 8, 1699–1710, https://doi.org/10.5194/tc-8-1699-2014, 2014b. a
Seroussi, H., Nakayama, Y., Larour, E., Menemenlis, D., Morlighem, M., Rignot, E., and Khazendar, A.: Continued retreat of Thwaites Glacier, West Antarctica, controlled by bed topography and ocean circulation, Geophys. Res. Lett., 44, 6191–6199, https://doi.org/10.1002/2017GL072910, 2017. 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
Shewchuk, J. R.: Triangle: Engineering a 2D Quality Mesh Generator and
Delaunay Triangulator, in: Applied Computational Geometry: Towards Geometric
Engineering, edited by: Lin, M. C. and Manocha, D., vol. 1148 of Lecture
Notes in Computer Science, Springer, 203–222, 1996. a
Van Wessem, J., Reijmer, C., Morlighem, M., Mouginot, J., Rignot, E., Medley, B., Joughin, I., Wouters, B., Depoorter, M., Bamber, J., Lenaerts, J. T. M., Van De Berg, W. J., Van Den Broeke, M. R., and Van Meijgaard, E.: Improved representation of East Antarctic surface mass balance in a regional atmospheric climate model, J. Glaciol., 60, 761–770, https://doi.org/10.3189/2014JoG14J051, 2014. a
Watanabe, N. and Kolditz, O.: Numerical stability analysis of two-dimensional solute transport along a discrete fracture in a porous rock matrix, Water Resour. Res., 51, 5855–5868, https://doi.org/10.1002/2015WR017164, 2015. a
Weertman, J.: Stability of the junction of an ice sheet and an ice shelf, J. Glaciol., 13, 3–11, https://doi.org/10.3189/S0022143000023327, 1974.
a
Zalesak, S. T.: Fully multidimensional flux-corrected transport algorithms for fluids, J. Comput. Phys., 31, 335–362, https://doi.org/10.1016/0021-9991(79)90051-2, 1979. a, b
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.
Numerical models are routinely used to understand the past and future behavior of ice sheets in...
