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.
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.
Gong Cheng, Mansa Krishna, and Mathieu Morlighem
Geosci. Model Dev., 18, 5311–5327, https://doi.org/10.5194/gmd-18-5311-2025, https://doi.org/10.5194/gmd-18-5311-2025, 2025
Short summary
Short summary
Predicting ice sheet contributions to sea level rise is challenging due to limited data and uncertainties in key processes. Traditional models require complex methods that lack flexibility. We developed PINNICLE (Physics-Informed Neural Networks for Ice and CLimatE), an open-source Python library that integrates machine learning with physical laws to improve ice sheet modeling. By combining data and physics, PINNICLE enhances predictions and adaptability, providing a powerful tool for climate research and sea level rise projections.
Felicity A. Holmes, Jamie Barnett, Henning Åkesson, Mathieu Morlighem, Johan Nilsson, Nina Kirchner, and Martin Jakobsson
The Cryosphere, 19, 2695–2714, https://doi.org/10.5194/tc-19-2695-2025, https://doi.org/10.5194/tc-19-2695-2025, 2025
Short summary
Short summary
Northern Greenland contains some of the ice sheet's last remaining glaciers with floating ice tongues. One of these is Ryder Glacier, which has been relatively stable in recent decades, in contrast to nearby glaciers. Here, we use a computer model to simulate Ryder Glacier until 2300 under both a low- and a high-emissions scenario. Very high levels of surface melt under a high-emissions future lead to a sea level rise contribution that is an order of magnitude higher than under a low-emissions future.
Daniel Abele, Thomas Kleiner, Yannic Fischler, Benjamin Uekermann, Gerasimos Chourdakis, Mathieu Morlighem, Achim Basermann, Christian Bischof, Hans-Joachim Bungartz, and Angelika Humbert
EGUsphere, https://doi.org/10.5194/egusphere-2025-3345, https://doi.org/10.5194/egusphere-2025-3345, 2025
This preprint is open for discussion and under review for Geoscientific Model Development (GMD).
Short summary
Short summary
For accurate projections of the evolution of continental ice sheets in Greenland and Antartica, interactions between the ice and its environment must be included in simulations. For this purpose, we have implemented adapters for the ice sheet model ISSM and subglacial hydrology model CUAS-MPI for the coupling library preCICE. This simplifies the study of earth systems by allowing the models to interact with each other as well as with models of the oceans or atmosphere with very little effort.
Younghyun Koo, Gong Cheng, Mathieu Morlighem, and Maryam Rahnemoonfar
The Cryosphere, 19, 2583–2599, https://doi.org/10.5194/tc-19-2583-2025, https://doi.org/10.5194/tc-19-2583-2025, 2025
Short summary
Short summary
Calving, the breaking of ice bodies from the terminus of a glacier, plays an important role in the mass losses of Greenland ice sheets. However, calving parameters have been poorly understood because of the intensive computational demands of traditional numerical models. To address this issue and find the optimal calving parameter that best represents real observations, we develop deep-learning emulators based on graph neural network architectures.
Shfaqat A. Khan, Helene Seroussi, Mathieu Morlighem, William Colgan, Veit Helm, Gong Cheng, Danjal Berg, Valentina R. Barletta, Nicolaj K. Larsen, William Kochtitzky, Michiel van den Broeke, Kurt H. Kjær, Andy Aschwanden, Brice Noël, Jason E. Box, Joseph A. MacGregor, Robert S. Fausto, Kenneth D. Mankoff, Ian M. Howat, Kuba Oniszk, Dominik Fahrner, Anja Løkkegaard, Eigil Y. H. Lippert, Alicia Bråtner, and Javed Hassan
Earth Syst. Sci. Data, 17, 3047–3071, https://doi.org/10.5194/essd-17-3047-2025, https://doi.org/10.5194/essd-17-3047-2025, 2025
Short summary
Short summary
The surface elevation of the Greenland Ice Sheet is changing due to surface mass balance processes and ice dynamics, each exhibiting distinct spatiotemporal patterns. Here, we employ satellite and airborne altimetry data with fine spatial (1 km) and temporal (monthly) resolutions to document this spatiotemporal evolution from 2003 to 2023. This dataset of fine-resolution altimetry data in both space and time will support studies of ice mass loss and be useful for GIS ice sheet modeling.
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.
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.
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.
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.
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...