Articles | Volume 14, issue 5
Geosci. Model Dev., 14, 2545–2573, 2021
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.
Geosci. Model Dev., 14, 2545–2573, 2021
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
Assessment of numerical schemes for transient, finite-element ice flow models using ISSM v4.18
Thiago Dias dos Santos et al.
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.
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.
Vincent Verjans, Alexander Robel, Helene Seroussi, Lizz Ultee, and Andrew Thompson
EGUsphere, https://doi.org/10.5194/egusphere-2022-699, https://doi.org/10.5194/egusphere-2022-699, 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 to represent the impacts of inherently uncertain processes on ice sheets. 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.
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.
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.
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.
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
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)
SnowClim v1.0: high-resolution snow model and data for the western United States
Snow Multidata Mapping and Modeling (S3M) 5.1: a distributed cryospheric model with dry and wet snow, data assimilation, glacier mass balance, and debris-driven melt
MPAS-Seaice (v1.0.0): sea-ice dynamics on unstructured Voronoi meshes
Explicitly modelling microtopography in permafrost landscapes in a land surface model (JULES vn5.4_microtopography)
Geometric remapping of particle distributions in the Discrete Element Model for Sea Ice (DEMSI v0.0)
Mapping high-resolution basal topography of West Antarctica from radar data using non-stationary multiple-point geostatistics (MPS-BedMappingV1)
NEMO-Bohai 1.0: a high-resolution ocean and sea ice modelling system for the Bohai Sea, China
An improved regional coupled modeling system for Arctic sea ice simulation and prediction: a case study for 2018
WIFF1.0: a hybrid machine-learning-based parameterization of wave-induced sea ice floe fracture
The Whole Antarctic Ocean Model (WAOM v1.0): development and evaluation
SNICAR-ADv3: a community tool for modeling spectral snow albedo
STEMMUS-UEB v1.0.0: integrated modeling of snowpack and soil water and energy transfer with three complexity levels of soil physical processes
A versatile method for computing optimized snow albedo from spectrally fixed radiative variables: VALHALLA v1.0
Ice Algae Model Intercomparison Project phase 2 (IAMIP2)
Improved representation of the contemporary Greenland ice sheet firn layer by IMAU-FDM v1.2G
A Gaussian process emulator for simulating ice sheet–climate interactions on a multi-million-year timescale: CLISEMv1.0
SITool (v1.0) – a new evaluation tool for large-scale sea ice simulations: application to CMIP6 OMIP
fenics_ice 1.0: a framework for quantifying initialization uncertainty for time-dependent ice sheet models
Development of adjoint-based ocean state estimation for the Amundsen and Bellingshausen seas and ice shelf cavities using MITgcm–ECCO (66j)
Sensitivity of Northern Hemisphere climate to ice–ocean interface heat flux parameterizations
icepack: a new glacier flow modeling package in Python, version 1.0
Benefits of sea ice initialization for the interannual-to-decadal climate prediction skill in the Arctic in EC-Earth3
Coupling framework (1.0) for the PISM (1.1.4) ice sheet model and the MOM5 (5.1.0) ocean model via the PICO ice shelf cavity model in an Antarctic domain
Performance of MAR (v3.11) in simulating the drifting-snow climate and surface mass balance of Adélie Land, East Antarctica
The Utrecht Finite Volume Ice-Sheet Model: UFEMISM (version 1.0)
PERICLIMv1.0: a model deriving palaeo-air temperatures from thaw depth in past permafrost regions
Assessing the simulated soil hydrothermal regime of the active layer from the Noah-MP land surface model (v1.1) in the permafrost regions of the Qinghai–Tibet Plateau
CrocO_v1.0: a particle filter to assimilate snowpack observations in a spatialised framework
A fully coupled Arctic sea-ice–ocean–atmosphere model (ArcIOAM v1.0) based on C-Coupler2: model description and preliminary results
The Framework For Ice Sheet–Ocean Coupling (FISOC) V1.1
Comparison of sea ice kinematics at different resolutions modeled with a grid hierarchy in the Community Earth System Model (version 1.2.1)
Snow profile alignment and similarity assessment for aggregating, clustering, and evaluating snowpack model output for avalanche forecasting
Improvements in one-dimensional grounding-line parameterizations in an ice-sheet model with lateral variations (PSUICE3D v2.1)
Implementation of the RCIP scheme and its performance for 1-D age computations in ice-sheet models
COSIPY v1.3 – an open-source coupled snowpack and ice surface energy and mass balance model
Using Arctic ice mass balance buoys for evaluation of modelled ice energy fluxes
Impact of the ice thickness distribution discretization on the sea ice concentration variability in the NEMO3.6–LIM3 global ocean–sea ice model
Simulating the Early Holocene demise of the Laurentide Ice Sheet with BISICLES (public trunk revision 3298)
Extended enthalpy formulations in the Ice-sheet and Sea-level System Model (ISSM) version 4.17: discontinuous conductivity and anisotropic streamline upwind Petrov–Galerkin (SUPG) method
The Community Firn Model (CFM) v1.0
Description and validation of the ice-sheet model Yelmo (version 1.0)
Evaluating integrated surface/subsurface permafrost thermal hydrology models in ATS (v0.88) against observations from a polygonal tundra site
SICOPOLIS-AD v1: an open-source adjoint modeling framework for ice sheet simulation enabled by the algorithmic differentiation tool OpenAD
On the calculation of normalized viscous–plastic sea ice stresses
Modelling thermomechanical ice deformation using an implicit pseudo-transient method (FastICE v1.0) based on graphical processing units (GPUs)
Version 1 of a sea ice module for the physics-based, detailed, multi-layer SNOWPACK model
A module to convert spectral to narrowband snow albedo for use in climate models: SNOWBAL v1.2
On the discretization of the ice thickness distribution in the NEMO3.6-LIM3 global ocean–sea ice model
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.
Abby C. Lute, John Abatzoglou, and Timothy Link
Geosci. Model Dev., 15, 5045–5071, https://doi.org/10.5194/gmd-15-5045-2022, https://doi.org/10.5194/gmd-15-5045-2022, 2022
Short summary
Short summary
We developed a snow model that can be used to quantify snowpack over large areas with a high degree of spatial detail. We ran the model over the western United States, creating a snow and climate dataset for three time periods. Compared to observations of snowpack, the model captured the key aspects of snow across time and space. The model and dataset will be useful in understanding historical and future changes in snowpack, with relevance to water resources, agriculture, and ecosystems.
Francesco Avanzi, Simone Gabellani, Fabio Delogu, Francesco Silvestro, Edoardo Cremonese, Umberto Morra di Cella, Sara Ratto, and Hervé Stevenin
Geosci. Model Dev., 15, 4853–4879, https://doi.org/10.5194/gmd-15-4853-2022, https://doi.org/10.5194/gmd-15-4853-2022, 2022
Short summary
Short summary
Knowing in real time how much snow and glacier ice has accumulated across the landscape has significant implications for water-resource management and flood control. This paper presents a computer model – S3M – allowing scientists and decision makers to predict snow and ice accumulation during winter and the subsequent melt during spring and summer. S3M has been employed for real-world flood forecasting since the early 2000s but is here being made open source for the first time.
Adrian K. Turner, William H. Lipscomb, Elizabeth C. Hunke, Douglas W. Jacobsen, Nicole Jeffery, Darren Engwirda, Todd D. Ringler, and Jonathan D. Wolfe
Geosci. Model Dev., 15, 3721–3751, https://doi.org/10.5194/gmd-15-3721-2022, https://doi.org/10.5194/gmd-15-3721-2022, 2022
Short summary
Short summary
We present the dynamical core of the MPAS-Seaice model, which uses a mesh consisting of a Voronoi tessellation with polygonal cells. Such a mesh allows variable mesh resolution in different parts of the domain and the focusing of computational resources in regions of interest. We describe the velocity solver and tracer transport schemes used and examine errors generated by the model in both idealized and realistic test cases and examine the computational efficiency of the model.
Noah D. Smith, Eleanor J. Burke, Kjetil Schanke Aas, Inge H. J. Althuizen, Julia Boike, Casper Tai Christiansen, Bernd Etzelmüller, Thomas Friborg, Hanna Lee, Heather Rumbold, Rachael H. Turton, Sebastian Westermann, and Sarah E. Chadburn
Geosci. Model Dev., 15, 3603–3639, https://doi.org/10.5194/gmd-15-3603-2022, https://doi.org/10.5194/gmd-15-3603-2022, 2022
Short summary
Short summary
The Arctic has large areas of small mounds that are caused by ice lifting up the soil. Snow blown by wind gathers in hollows next to these mounds, insulating them in winter. The hollows tend to be wetter, and thus the soil absorbs more heat in summer. The warm wet soil in the hollows decomposes, releasing methane. We have made a model of this, and we have tested how it behaves and whether it looks like sites in Scandinavia and Siberia. Sometimes we get more methane than a model without mounds.
Adrian K. Turner, Kara J. Peterson, and Dan Bolintineanu
Geosci. Model Dev., 15, 1953–1970, https://doi.org/10.5194/gmd-15-1953-2022, https://doi.org/10.5194/gmd-15-1953-2022, 2022
Short summary
Short summary
We developed a technique to remap sea ice tracer quantities between circular discrete element distributions. This is needed for a global discrete element method sea ice model being developed jointly by Los Alamos National Laboratory and Sandia National Laboratories that has the potential to better utilize newer supercomputers with graphics processing units and better represent sea ice dynamics. This new remapping technique ameliorates the effect of element distortion created by sea ice ridging.
Zhen Yin, Chen Zuo, Emma J. MacKie, and Jef Caers
Geosci. Model Dev., 15, 1477–1497, https://doi.org/10.5194/gmd-15-1477-2022, https://doi.org/10.5194/gmd-15-1477-2022, 2022
Short summary
Short summary
We provide a multiple-point geostatistics approach to probabilistically learn from training images to fill large-scale irregular geophysical data gaps. With a repository of global topographic training images, our approach models high-resolution basal topography and quantifies the geospatial uncertainty. It generated high-resolution topographic realizations to investigate the impact of basal topographic uncertainty on critical subglacial hydrological flow patterns associated with ice velocity.
Yu Yan, Wei Gu, Andrea M. U. Gierisch, Yingjun Xu, and Petteri Uotila
Geosci. Model Dev., 15, 1269–1288, https://doi.org/10.5194/gmd-15-1269-2022, https://doi.org/10.5194/gmd-15-1269-2022, 2022
Short summary
Short summary
In this study, we developed NEMO-Bohai, an ocean–ice model for the Bohai Sea, China. This study presented the scientific design and technical choices of the parameterizations for the NEMO-Bohai model. The model was calibrated and evaluated with in situ and satellite observations of ocean and sea ice. NEMO-Bohai is intended to be a valuable tool for long-term ocean and ice simulations and climate change studies.
Chao-Yuan Yang, Jiping Liu, and Dake Chen
Geosci. Model Dev., 15, 1155–1176, https://doi.org/10.5194/gmd-15-1155-2022, https://doi.org/10.5194/gmd-15-1155-2022, 2022
Short summary
Short summary
We present an improved coupled modeling system for Arctic sea ice prediction. We perform Arctic sea ice prediction experiments with improved/updated physical parameterizations, which show better skill in predicting sea ice state as well as atmospheric and oceanic state in the Arctic compared with its predecessor. The improved model also shows extended predictive skill of Arctic sea ice after the summer season. This provides an added value of this prediction system for decision-making.
Christopher Horvat and Lettie A. Roach
Geosci. Model Dev., 15, 803–814, https://doi.org/10.5194/gmd-15-803-2022, https://doi.org/10.5194/gmd-15-803-2022, 2022
Short summary
Short summary
Sea ice is a composite of individual pieces, called floes, ranging in horizontal size from meters to kilometers. Variations in sea ice geometry are often forced by ocean waves, a process that is an important target of global climate models as it affects the rate of sea ice melting. Yet directly simulating these interactions is computationally expensive. We present a neural-network-based model of wave–ice fracture that allows models to incorporate their effect without added computational cost.
Ole Richter, David E. Gwyther, Benjamin K. Galton-Fenzi, and Kaitlin A. Naughten
Geosci. Model Dev., 15, 617–647, https://doi.org/10.5194/gmd-15-617-2022, https://doi.org/10.5194/gmd-15-617-2022, 2022
Short summary
Short summary
Here we present an improved model of the Antarctic continental shelf ocean and demonstrate that it is capable of reproducing present-day conditions. The improvements are fundamental and regard the inclusion of tides and ocean eddies. We conclude that the model is well suited to gain new insights into processes that are important for Antarctic ice sheet retreat and global ocean changes. Hence, the model will ultimately help to improve projections of sea level rise and climate change.
Mark G. Flanner, Julian B. Arnheim, Joseph M. Cook, Cheng Dang, Cenlin He, Xianglei Huang, Deepak Singh, S. McKenzie Skiles, Chloe A. Whicker, and Charles S. Zender
Geosci. Model Dev., 14, 7673–7704, https://doi.org/10.5194/gmd-14-7673-2021, https://doi.org/10.5194/gmd-14-7673-2021, 2021
Short summary
Short summary
We present the technical formulation and evaluation of a publicly available code and web-based model to simulate the spectral albedo of snow. Our model accounts for numerous features of the snow state and ambient conditions, including the the presence of light-absorbing matter like black and brown carbon, mineral dust, volcanic ash, and snow algae. Carbon dioxide snow, found on Mars, is also represented. The model accurately reproduces spectral measurements of clean and contaminated snow.
Lianyu Yu, Yijian Zeng, and Zhongbo Su
Geosci. Model Dev., 14, 7345–7376, https://doi.org/10.5194/gmd-14-7345-2021, https://doi.org/10.5194/gmd-14-7345-2021, 2021
Short summary
Short summary
We developed an integrated soil–snow–atmosphere model (STEMMUS-UEB) dedicated to the physical description of snow and soil processes with various complexities. With STEMMUS-UEB, we demonstrated that the snowpack affects not only the soil surface moisture conditions (in the liquid and ice phase) and energy-related states (albedo, LE) but also the subsurface soil water and vapor transfer, which contributes to a better understanding of the hydrothermal implications of the snowpack in cold regions.
Florent Veillon, Marie Dumont, Charles Amory, and Mathieu Fructus
Geosci. Model Dev., 14, 7329–7343, https://doi.org/10.5194/gmd-14-7329-2021, https://doi.org/10.5194/gmd-14-7329-2021, 2021
Short summary
Short summary
In climate models, the snow albedo scheme generally calculates only a narrowband or broadband albedo. Therefore, we have developed the VALHALLA method to optimize snow spectral albedo calculations through the determination of spectrally fixed radiative variables. The development of VALHALLA v1.0 with the use of the snow albedo model TARTES and the spectral irradiance model SBDART indicates a considerable reduction in calculation time while maintaining an adequate accuracy of albedo values.
Hakase Hayashida, Meibing Jin, Nadja S. Steiner, Neil C. Swart, Eiji Watanabe, Russell Fiedler, Andrew McC. Hogg, Andrew E. Kiss, Richard J. Matear, and Peter G. Strutton
Geosci. Model Dev., 14, 6847–6861, https://doi.org/10.5194/gmd-14-6847-2021, https://doi.org/10.5194/gmd-14-6847-2021, 2021
Short summary
Short summary
Ice algae are tiny plants like phytoplankton but they grow within sea ice. In polar regions, both phytoplankton and ice algae are the foundation of marine ecosystems and play an important role in taking up carbon dioxide in the atmosphere. However, state-of-the-art climate models typically do not include ice algae, and therefore their role in the climate system remains unclear. This project aims to address this knowledge gap by coordinating a set of experiments using sea-ice–ocean models.
Max Brils, Peter Kuipers Munneke, Willem Jan van de Berg, and Michiel van den Broeke
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2021-303, https://doi.org/10.5194/gmd-2021-303, 2021
Revised manuscript accepted for GMD
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 its 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.
Jonas Van Breedam, Philippe Huybrechts, and Michel Crucifix
Geosci. Model Dev., 14, 6373–6401, https://doi.org/10.5194/gmd-14-6373-2021, https://doi.org/10.5194/gmd-14-6373-2021, 2021
Short summary
Short summary
Ice sheets are an important component of the climate system and interact with the atmosphere through albedo variations and changes in the surface height. On very long timescales, it is impossible to directly couple ice sheet models with climate models and other techniques have to be used. Here we present a novel coupling method between ice sheets and the atmosphere by making use of an emulator to simulate ice sheet–climate interactions for several million years.
Xia Lin, François Massonnet, Thierry Fichefet, and Martin Vancoppenolle
Geosci. Model Dev., 14, 6331–6354, https://doi.org/10.5194/gmd-14-6331-2021, https://doi.org/10.5194/gmd-14-6331-2021, 2021
Short summary
Short summary
This study introduces a new Sea Ice Evaluation Tool (SITool) to evaluate the model skills on the bipolar sea ice simulations by providing performance metrics and diagnostics. SITool is applied to evaluate the CMIP6 OMIP simulations. By changing the atmospheric forcing from CORE-II to JRA55-do data, many aspects of sea ice simulations are improved. SITool will be useful for helping teams managing various versions of a sea ice model or tracking the time evolution of model performance.
Conrad P. Koziol, Joe A. Todd, Daniel N. Goldberg, and James R. Maddison
Geosci. Model Dev., 14, 5843–5861, https://doi.org/10.5194/gmd-14-5843-2021, https://doi.org/10.5194/gmd-14-5843-2021, 2021
Short summary
Short summary
Sea level change due to the loss of ice sheets presents great risk for coastal communities. Models are used to forecast ice loss, but their evolution depends strongly on properties which are hidden from observation and must be inferred from satellite observations. Common methods for doing so do not allow for quantification of the uncertainty inherent or how it will affect forecasts. We provide a framework for quantifying how this
initialization uncertaintyaffects ice loss forecasts.
Yoshihiro Nakayama, Dimitris Menemenlis, Ou Wang, Hong Zhang, Ian Fenty, and An T. Nguyen
Geosci. Model Dev., 14, 4909–4924, https://doi.org/10.5194/gmd-14-4909-2021, https://doi.org/10.5194/gmd-14-4909-2021, 2021
Short summary
Short summary
High ice shelf melting in the Amundsen Sea has attracted many observational campaigns in the past decade. One method to combine observations with numerical models is the adjoint method. After 20 iterations, the cost function, defined as a sum of the weighted model–data difference, is reduced by 65 % by adjusting initial conditions, atmospheric forcing, and vertical diffusivity. This study demonstrates adjoint-method optimization with explicit representation of ice shelf cavity circulation.
Xiaoxu Shi, Dirk Notz, Jiping Liu, Hu Yang, and Gerrit Lohmann
Geosci. Model Dev., 14, 4891–4908, https://doi.org/10.5194/gmd-14-4891-2021, https://doi.org/10.5194/gmd-14-4891-2021, 2021
Short summary
Short summary
The ice–ocean heat flux is one of the key elements controlling sea ice changes. It motivates our study, which aims to examine the responses of modeled climate to three ice–ocean heat flux parameterizations, including two old approaches that assume one-way heat transport and a new one describing a double-diffusive ice–ocean heat exchange. The results show pronounced differences in the modeled sea ice, ocean, and atmosphere states for the latter as compared to the former two parameterizations.
Daniel R. Shapero, Jessica A. Badgeley, Andrew O. Hoffman, and Ian R. Joughin
Geosci. Model Dev., 14, 4593–4616, https://doi.org/10.5194/gmd-14-4593-2021, https://doi.org/10.5194/gmd-14-4593-2021, 2021
Short summary
Short summary
This paper describes a new software package called "icepack" for modeling the flow of ice sheets and glaciers. Glaciologists use tools like icepack to better understand how ice sheets flow, what role they have played in shaping Earth's climate, and how much sea level rise we can expect in the coming decades to centuries. The icepack package includes several innovations to help researchers describe and solve interesting glaciological problems and to experiment with the underlying model physics.
Tian Tian, Shuting Yang, Mehdi Pasha Karami, François Massonnet, Tim Kruschke, and Torben Koenigk
Geosci. Model Dev., 14, 4283–4305, https://doi.org/10.5194/gmd-14-4283-2021, https://doi.org/10.5194/gmd-14-4283-2021, 2021
Short summary
Short summary
Three decadal prediction experiments with EC-Earth3 are performed to investigate the impact of ocean, sea ice concentration and thickness initialization, respectively. We find that the persistence of perennial thick ice in the central Arctic can affect the sea ice predictability in its adjacent waters via advection process or wind, despite those regions being seasonally ice free during two recent decades. This has implications for the coming decades as the thinning of Arctic sea ice continues.
Moritz Kreuzer, Ronja Reese, Willem Nicholas Huiskamp, Stefan Petri, Torsten Albrecht, Georg Feulner, and Ricarda Winkelmann
Geosci. Model Dev., 14, 3697–3714, https://doi.org/10.5194/gmd-14-3697-2021, https://doi.org/10.5194/gmd-14-3697-2021, 2021
Short summary
Short summary
We present the technical implementation of a coarse-resolution coupling between an ice sheet model and an ocean model that allows one to simulate ice–ocean interactions at timescales from centuries to millennia. As ice shelf cavities cannot be resolved in the ocean model at coarse resolution, we bridge the gap using an sub-shelf cavity module. It is shown that the framework is computationally efficient, conserves mass and energy, and can produce a stable coupled state under present-day forcing.
Charles Amory, Christoph Kittel, Louis Le Toumelin, Cécile Agosta, Alison Delhasse, Vincent Favier, and Xavier Fettweis
Geosci. Model Dev., 14, 3487–3510, https://doi.org/10.5194/gmd-14-3487-2021, https://doi.org/10.5194/gmd-14-3487-2021, 2021
Short summary
Short summary
This paper presents recent developments in the drifting-snow scheme of the regional climate model MAR and its application to simulate drifting snow and the surface mass balance of Adélie Land in East Antarctica. The model is extensively described and evaluated against a multi-year drifting-snow dataset and surface mass balance estimates available in the area. The model sensitivity to input parameters and improvements over a previously published version are also assessed.
Constantijn J. Berends, Heiko Goelzer, and Roderik S. W. van de Wal
Geosci. Model Dev., 14, 2443–2470, https://doi.org/10.5194/gmd-14-2443-2021, https://doi.org/10.5194/gmd-14-2443-2021, 2021
Short summary
Short summary
The largest uncertainty in projections of sea-level rise comes from ice-sheet retreat. To better understand how these ice sheets respond to the changing climate, ice-sheet models are used, which must be able to reproduce both their present and past evolution. We have created a model that is fast enough to simulate an ice sheet at a high resolution over the course of an entire 120 000-year glacial cycle. This allows us to study processes that cannot be captured by lower-resolution models.
Tomáš Uxa, Marek Křížek, and Filip Hrbáček
Geosci. Model Dev., 14, 1865–1884, https://doi.org/10.5194/gmd-14-1865-2021, https://doi.org/10.5194/gmd-14-1865-2021, 2021
Short summary
Short summary
We present a simple model that derives palaeo-air temperature characteristics related to the palaeo-active-layer thickness, which can be recognized using many relict periglacial features found in past permafrost regions. Its evaluation against modern temperature records and an experimental palaeo-air temperature reconstruction showed relatively high model accuracy, which suggests that it could become a useful tool for reconstructing Quaternary palaeo-environments.
Xiangfei Li, Tonghua Wu, Xiaodong Wu, Jie Chen, Xiaofan Zhu, Guojie Hu, Ren Li, Yongping Qiao, Cheng Yang, Junming Hao, Jie Ni, and Wensi Ma
Geosci. Model Dev., 14, 1753–1771, https://doi.org/10.5194/gmd-14-1753-2021, https://doi.org/10.5194/gmd-14-1753-2021, 2021
Short summary
Short summary
In this study, an ensemble simulation of 55296 scheme combinations for at a typical permafrost site on the Qinghai–Tibet Plateau (QTP) was conducted. The general performance of the Noah-MP model for snow cover events (SCEs), soil temperature (ST) and soil liquid water content (SLW) was assessed, and the sensitivities of parameterization schemes at different depths were investigated. We show that Noah-MP tends to overestimate SCEs and underestimate ST and topsoil SLW on the QTP.
Bertrand Cluzet, Matthieu Lafaysse, Emmanuel Cosme, Clément Albergel, Louis-François Meunier, and Marie Dumont
Geosci. Model Dev., 14, 1595–1614, https://doi.org/10.5194/gmd-14-1595-2021, https://doi.org/10.5194/gmd-14-1595-2021, 2021
Short summary
Short summary
In the mountains, the combination of large model error and observation sparseness is a challenge for data assimilation. Here, we develop two variants of the particle filter (PF) in order to propagate the information content of observations into unobserved areas. By adjusting observation errors or exploiting background correlation patterns, we demonstrate the potential for partial observations of snow depth and surface reflectance to improve model accuracy with the PF in an idealised setting.
Shihe Ren, Xi Liang, Qizhen Sun, Hao Yu, L. Bruno Tremblay, Bo Lin, Xiaoping Mai, Fu Zhao, Ming Li, Na Liu, Zhikun Chen, and Yunfei Zhang
Geosci. Model Dev., 14, 1101–1124, https://doi.org/10.5194/gmd-14-1101-2021, https://doi.org/10.5194/gmd-14-1101-2021, 2021
Short summary
Short summary
Sea ice plays a crucial role in global energy and water budgets. To get a better simulation of sea ice, we coupled a sea ice model with an atmospheric and ocean model to form a fully coupled system. The sea ice simulation results of this coupled system demonstrated that a two-way coupled model has better performance in terms of sea ice, especially in summer. This indicates that sea-ice–ocean–atmosphere interaction plays a crucial role in controlling Arctic summertime sea ice distribution.
Rupert Gladstone, Benjamin Galton-Fenzi, David Gwyther, Qin Zhou, Tore Hattermann, Chen Zhao, Lenneke Jong, Yuwei Xia, Xiaoran Guo, Konstantinos Petrakopoulos, Thomas Zwinger, Daniel Shapero, and John Moore
Geosci. Model Dev., 14, 889–905, https://doi.org/10.5194/gmd-14-889-2021, https://doi.org/10.5194/gmd-14-889-2021, 2021
Short summary
Short summary
Retreat of the Antarctic ice sheet, and hence its contribution to sea level rise, is highly sensitive to melting of its floating ice shelves. This melt is caused by warm ocean currents coming into contact with the ice. Computer models used for future ice sheet projections are not able to realistically evolve these melt rates. We describe a new coupling framework to enable ice sheet and ocean computer models to interact, allowing projection of the evolution of melt and its impact on sea level.
Shiming Xu, Jialiang Ma, Lu Zhou, Yan Zhang, Jiping Liu, and Bin Wang
Geosci. Model Dev., 14, 603–628, https://doi.org/10.5194/gmd-14-603-2021, https://doi.org/10.5194/gmd-14-603-2021, 2021
Short summary
Short summary
A multi-resolution tripolar grid hierarchy is constructed and integrated in CESM (version 1.2.1). The resolution range includes 0.45, 0.15, and 0.05°. Based on atmospherically forced sea ice experiments, the model simulates reasonable sea ice kinematics and scaling properties. Landfast ice thickness can also be systematically shifted due to non-convergent solutions to an
elastic–viscous–plastic (EVP) model. This work is a framework for multi-scale modeling of the ocean and sea ice with CESM.
Florian Herla, Simon Horton, Patrick Mair, and Pascal Haegeli
Geosci. Model Dev., 14, 239–258, https://doi.org/10.5194/gmd-14-239-2021, https://doi.org/10.5194/gmd-14-239-2021, 2021
Short summary
Short summary
The adoption of snowpack models in support of avalanche forecasting has been limited. To promote their operational application, we present a numerical method for processing multivariate snow stratigraphy profiles of mixed data types. Our algorithm enables applications like dynamical grouping and summarizing of model simulations, model evaluation, and data assimilation. By emulating the human analysis process, our approach will allow forecasters to familiarly interact with snowpack simulations.
David Pollard and Robert M. DeConto
Geosci. Model Dev., 13, 6481–6500, https://doi.org/10.5194/gmd-13-6481-2020, https://doi.org/10.5194/gmd-13-6481-2020, 2020
Short summary
Short summary
Buttressing by floating ice shelves at ice-sheet grounding lines is an
important process that affects ice retreat and whether structural failure
occurs in deep bathymetry. Here, we use a simple algorithm to better
represent 2-D grounding-line curvature in an ice-sheet model. Along with other
enhancements, this improves the performance in idealized-fjord intercomparisons
and enables better diagnosis of potential structural failure at future
retreating Antarctic grounding lines.
Fuyuki Saito, Takashi Obase, and Ayako Abe-Ouchi
Geosci. Model Dev., 13, 5875–5896, https://doi.org/10.5194/gmd-13-5875-2020, https://doi.org/10.5194/gmd-13-5875-2020, 2020
Short summary
Short summary
The present study introduces the rational function-based constrained interpolation profile (RCIP) method for use in 1 d dating computations in ice sheets and demonstrates the performance of the scheme. Comparisons are examined among the RCIP schemes and the first- and second-order upwind schemes. The results show that, in particular, the RCIP scheme preserves the pattern of input histories, in terms of the profile of internal annual layer thickness, better than the other schemes.
Tobias Sauter, Anselm Arndt, and Christoph Schneider
Geosci. Model Dev., 13, 5645–5662, https://doi.org/10.5194/gmd-13-5645-2020, https://doi.org/10.5194/gmd-13-5645-2020, 2020
Short summary
Short summary
Glacial changes play a key role from a socioeconomic, political, and scientific point of view. Here, we present the open-source coupled snowpack and ice surface energy and mass balance model, which provides a lean, flexible, and user-friendly framework for modeling distributed snow and glacier mass changes. The model provides a suitable platform for sensitivity, detection, and attribution analyses for glacier changes and a tool for quantifying inherent uncertainties.
Alex West, Mat Collins, and Ed Blockley
Geosci. Model Dev., 13, 4845–4868, https://doi.org/10.5194/gmd-13-4845-2020, https://doi.org/10.5194/gmd-13-4845-2020, 2020
Short summary
Short summary
This study calculates sea ice energy fluxes from data produced by ice mass balance buoys (devices measuring ice elevation and temperature). It is shown how the resulting dataset can be used to evaluate a coupled climate model (HadGEM2-ES), with biases in the energy fluxes seen to be consistent with biases in the sea ice state and surface radiation. This method has potential to improve sea ice model evaluation, so as to better understand spread in model simulations of sea ice state.
Eduardo Moreno-Chamarro, Pablo Ortega, and François Massonnet
Geosci. Model Dev., 13, 4773–4787, https://doi.org/10.5194/gmd-13-4773-2020, https://doi.org/10.5194/gmd-13-4773-2020, 2020
Short summary
Short summary
Climate models need to capture sea ice complexity to represent it realistically. Here we assess how distributing sea ice in discrete thickness categories impacts how sea ice variability is simulated in the NEMO3.6–LIM3 model. Simulations and satellite observations are compared by using k-means clustering of sea ice concentration in winter and summer between 1979 and 2014 at both poles. Little improvements in the modeled sea ice lead us to recommend using the standard number of five categories.
Ilkka S. O. Matero, Lauren J. Gregoire, and Ruza F. Ivanovic
Geosci. Model Dev., 13, 4555–4577, https://doi.org/10.5194/gmd-13-4555-2020, https://doi.org/10.5194/gmd-13-4555-2020, 2020
Short summary
Short summary
The Northern Hemisphere cooled by several degrees for a century 8000 years ago due to the collapse of an ice sheet in North America that released large amounts of meltwater into the North Atlantic and slowed down its circulation. We numerically model the ice sheet to understand its evolution during this event. Our results match data thanks to good ice dynamics but depend mostly on surface melt and snowfall. Further work will help us understand how past and future ice melt affects climate.
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.
C. Max Stevens, Vincent Verjans, Jessica M. D. Lundin, Emma C. Kahle, Annika N. Horlings, Brita I. Horlings, and Edwin D. Waddington
Geosci. Model Dev., 13, 4355–4377, https://doi.org/10.5194/gmd-13-4355-2020, https://doi.org/10.5194/gmd-13-4355-2020, 2020
Short summary
Short summary
Understanding processes in snow (firn), including compaction and airflow, is important for calculating how much mass the ice sheets are losing and for interpreting climate records from ice cores. We have developed the open-source Community Firn Model to simulate these processes. We used it to compare 13 different firn compaction equations and found that they do not agree within 10 %. We also show that including firn compaction in a firn-air model improves the match with data from ice cores.
Alexander Robinson, Jorge Alvarez-Solas, Marisa Montoya, Heiko Goelzer, Ralf Greve, and Catherine Ritz
Geosci. Model Dev., 13, 2805–2823, https://doi.org/10.5194/gmd-13-2805-2020, https://doi.org/10.5194/gmd-13-2805-2020, 2020
Short summary
Short summary
Here we describe Yelmo v1.0, an intuitive and state-of-the-art hybrid ice sheet model. The model design and physics are described, and benchmark simulations are provided to validate its performance. Yelmo is a versatile ice sheet model that can be applied to a wide variety of problems.
Ahmad Jan, Ethan T. Coon, and Scott L. Painter
Geosci. Model Dev., 13, 2259–2276, https://doi.org/10.5194/gmd-13-2259-2020, https://doi.org/10.5194/gmd-13-2259-2020, 2020
Short summary
Short summary
Computer simulations are important tools for understanding the response of Arctic permafrost to a warming climate. To build confidence in an emerging class of permafrost simulators, we evaluated the Advanced Terrestrial Simulator against field observations from a frozen tundra site near Utqiaġvik (formerly Barrow), Alaska. The 3-year simulations agree well with observations of snow depth, summer water table, soil temperature at multiple locations, and spatially averaged evaporation.
Liz C. Logan, Sri Hari Krishna Narayanan, Ralf Greve, and Patrick Heimbach
Geosci. Model Dev., 13, 1845–1864, https://doi.org/10.5194/gmd-13-1845-2020, https://doi.org/10.5194/gmd-13-1845-2020, 2020
Short summary
Short summary
A new capability has been developed for the ice sheet model SICOPOLIS (SImulation COde for POLythermal Ice Sheets) that enables the generation of derivative code, such as tangent linear or adjoint models, by means of algorithmic differentiation. It relies on the source transformation algorithmic (AD) differentiation tool OpenAD. The reverse mode of AD provides the adjoint model, SICOPOLIS-AD, which may be applied for comprehensive sensitivity analyses as well as gradient-based optimization.
Jean-François Lemieux and Frédéric Dupont
Geosci. Model Dev., 13, 1763–1769, https://doi.org/10.5194/gmd-13-1763-2020, https://doi.org/10.5194/gmd-13-1763-2020, 2020
Short summary
Short summary
Sea ice dynamics plays an important role in shaping the sea cover in polar regions. Winds and ocean currents exert large stresses on the sea ice cover. This can lead to the formation of long cracks and ridges, which strongly impact the exchange of heat, momentum and moisture between the atmosphere and the ocean. It is therefore crucial for a sea ice model to be able to represent these features. This article describes how internal sea ice stresses should be diagnosed from model simulations.
Ludovic Räss, Aleksandar Licul, Frédéric Herman, Yury Y. Podladchikov, and Jenny Suckale
Geosci. Model Dev., 13, 955–976, https://doi.org/10.5194/gmd-13-955-2020, https://doi.org/10.5194/gmd-13-955-2020, 2020
Short summary
Short summary
Accurate predictions of future sea level rise require numerical models that predict rapidly deforming ice. Localised ice deformation can be captured numerically only with high temporal and spatial resolution. This paper’s goal is to propose a parallel FastICE solver for modelling ice deformation. Our model is particularly useful for improving our process-based understanding of localised ice deformation. Our solver reaches a parallel efficiency of 99 % on GPU-based supercomputers.
Nander Wever, Leonard Rossmann, Nina Maaß, Katherine C. Leonard, Lars Kaleschke, Marcel Nicolaus, and Michael Lehning
Geosci. Model Dev., 13, 99–119, https://doi.org/10.5194/gmd-13-99-2020, https://doi.org/10.5194/gmd-13-99-2020, 2020
Short summary
Short summary
Sea ice is an important component of the global climate system. The presence of a snow layer covering sea ice can impact ice mass and energy budgets. The detailed, physics-based, multi-layer snow model SNOWPACK was modified to simulate the snow–sea-ice system, providing simulations of the snow microstructure, water percolation and flooding, and superimposed ice formation. The model is applied to in situ measurements from snow and ice mass-balance buoys installed in the Antarctic Weddell Sea.
Christiaan T. van Dalum, Willem Jan van de Berg, Quentin Libois, Ghislain Picard, and Michiel R. van den Broeke
Geosci. Model Dev., 12, 5157–5175, https://doi.org/10.5194/gmd-12-5157-2019, https://doi.org/10.5194/gmd-12-5157-2019, 2019
Short summary
Short summary
Climate models are often limited to relatively simple snow albedo schemes. Therefore, we have developed the SNOWBAL module to couple a climate model with a physically based wavelength dependent snow albedo model. Using SNOWBAL v1.2 to couple the snow albedo model TARTES with the regional climate model RACMO2 indicates a potential performance gain for the Greenland ice sheet.
François Massonnet, Antoine Barthélemy, Koffi Worou, Thierry Fichefet, Martin Vancoppenolle, Clément Rousset, and Eduardo Moreno-Chamarro
Geosci. Model Dev., 12, 3745–3758, https://doi.org/10.5194/gmd-12-3745-2019, https://doi.org/10.5194/gmd-12-3745-2019, 2019
Short summary
Short summary
Sea ice thickness varies considerably on spatial scales of several meters. However, contemporary climate models cannot resolve such scales yet. This is why sea ice models used in climate models include an ice thickness distribution (ITD) to account for this unresolved variability. Here, we explore with the ocean–sea ice model NEMO3.6-LIM3 the sensitivity of simulated mean Arctic and Antarctic sea ice states to the way the ITD is discretized.
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...
