Articles | Volume 13, issue 4
https://doi.org/10.5194/gmd-13-1845-2020
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/gmd-13-1845-2020
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Development and technical paper
| 
09 Apr 2020
Development and technical paper |
| 09 Apr 2020
| 09 Apr 2020
SICOPOLIS-AD v1: an open-source adjoint modeling framework for ice sheet simulation enabled by the algorithmic differentiation tool OpenAD
Liz C. Logan
Oden Institute for Computational Engineering and Sciences, University of Texas at Austin, 201 East 24th Street, Austin, TX 78712, USA
Sri Hari Krishna Narayanan
Mathematics and Computer Science Division, Argonne National Laboratory, Lemont, IL 60439, USA
Ralf Greve
Institute of Low Temperature Science, Hokkaido University, Kita-19, Nishi-8, Kita-ku, Sapporo 060-0819, Japan
Arctic Research Center, Hokkaido University, Kita-21, Nishi-11, Kita-ku, Sapporo 001-0021, Japan
Oden Institute for Computational Engineering and Sciences, University of Texas at Austin, 201 East 24th Street, Austin, TX 78712, USA
Jackson School of Geosciences, University of Texas at Austin, 201 East 24th Street, Austin, TX 78712, USA
Institute for Geophysics, University of Texas at Austin, J.J. Pickle Research Campus, Bldg. 196, 10100 Burnet Road (R2200), Room 2.236, Austin, TX 78758, USA
Related authors
Liz C. Logan, Luc L. Lavier, Eunseo Choi, Eh Tan, and Ginny A. Catania
The Cryosphere, 11, 117–132, https://doi.org/10.5194/tc-11-117-2017, https://doi.org/10.5194/tc-11-117-2017, 2017
Short summary
Short summary
Global sea level rise prediction is a pressing and unresolved problem, one whose solution depends upon glaciologists better predicting ice sheet shrinkage due to iceberg calving. We present a numerical model that is capable of simulating ice flow and breakage that leads to iceberg calving and find that a material property that captures both the fluid- and solid-like behavior of ice simultaneously is a necessary condition for studying areas of glaciers in contact with ocean water prone to calve.
Abdullah A. Fahad, Andrea Molod, Krzysztof Wargan, Dimitris Menemenlis, Patrick Heimbach, Atanas Trayanov, Ehud Strobach, and Lawrence Coy
EGUsphere, https://doi.org/10.21203/rs.3.rs-1892797/v2, https://doi.org/10.21203/rs.3.rs-1892797/v2, 2025
This preprint is open for discussion and under review for Atmospheric Chemistry and Physics (ACP).
Short summary
Short summary
This study used a 1-degree GEOS-MITgcm coupled GCM to analyze the Northern Hemisphere (NH) stratospheric temperature response to external forcing. Results show the NH polar stratospheric temperature increased from 1992 to 2000, contrary to the expectation of stratospheric cooling with rising CO2. However, from 2000 to 2020, the temperature decreased. The study concluded that changes in CO2 and Ozone drive the meridional eddy transport of heat, dictating polar stratospheric temperature behavior.
Andrew Porter and Patrick Heimbach
State Planet Discuss., https://doi.org/10.5194/sp-2024-32, https://doi.org/10.5194/sp-2024-32, 2024
Revised manuscript accepted for SP
Short summary
Short summary
Numerical ocean forecasting is a key part of accurate models of the earth system. However, they require powerful computing resources and the architectures of the necessary computers are evolving rapidly. Unfortunately, this is a disruptive change – an ocean model must be modified to enable it to make use of this new computing hardware. This paper reviews what has been done in this area and identifies solutions to enable operational ocean forecasts to make use of the new computing hardware.
Patrick Heimbach, Fearghal O'Donncha, Jose Maria Garcia-Valdecasas, Alain Arnaud, and Liying Wan
State Planet Discuss., https://doi.org/10.5194/sp-2024-18, https://doi.org/10.5194/sp-2024-18, 2024
Revised manuscript accepted for SP
Short summary
Short summary
Operational ocean prediction relies on computationally expensive numerical models and complex workflows known as data assimilation, in which models are fit to observations to produce optimal initial conditions for prediction. Machine learning has the potential to vastly accelerate ocean prediction, improve uncertainty quantification through massive surrogate model-based ensembles, and render simulations more accurate through better model calibration. We review developments and challenges.
Laurent Bertino, Patrick Heimbach, Ed Blockley, and Einar Ólason
State Planet Discuss., https://doi.org/10.5194/sp-2024-24, https://doi.org/10.5194/sp-2024-24, 2024
Revised manuscript accepted for SP
Short summary
Short summary
Forecasts of sea ice are in high demand in the polar regions, they are also quickly improving and becoming more easily accessible to non-experts. We provide here a brief status of the short-term forecasting services – typically 10 days ahead – and an outlook of their upcoming developments.
Matteo Willeit, Reinhard Calov, Stefanie Talento, Ralf Greve, Jorjo Bernales, Volker Klemann, Meike Bagge, and Andrey Ganopolski
Clim. Past, 20, 597–623, https://doi.org/10.5194/cp-20-597-2024, https://doi.org/10.5194/cp-20-597-2024, 2024
Short summary
Short summary
We present transient simulations of the last glacial inception with the coupled climate–ice sheet model CLIMBER-X showing a rapid increase in Northern Hemisphere ice sheet area and a sea level drop by ~ 35 m, with the vegetation feedback playing a key role. Overall, our simulations confirm and refine previous results showing that climate-vegetation–cryosphere–carbon cycle feedbacks play a fundamental role in the transition from interglacial to glacial states.
Hélène Seroussi, Vincent Verjans, Sophie Nowicki, Antony J. Payne, Heiko Goelzer, William H. Lipscomb, Ayako Abe-Ouchi, Cécile Agosta, Torsten Albrecht, Xylar Asay-Davis, Alice Barthel, Reinhard Calov, Richard Cullather, Christophe Dumas, Benjamin K. Galton-Fenzi, Rupert Gladstone, Nicholas R. Golledge, Jonathan M. Gregory, Ralf Greve, Tore Hattermann, Matthew J. Hoffman, Angelika Humbert, Philippe Huybrechts, Nicolas C. Jourdain, Thomas Kleiner, Eric Larour, Gunter R. Leguy, Daniel P. Lowry, Chistopher M. Little, Mathieu Morlighem, Frank Pattyn, Tyler Pelle, Stephen F. Price, Aurélien Quiquet, Ronja Reese, Nicole-Jeanne Schlegel, Andrew Shepherd, Erika Simon, Robin S. Smith, Fiammetta Straneo, Sainan Sun, Luke D. Trusel, Jonas Van Breedam, Peter Van Katwyk, Roderik S. W. van de Wal, Ricarda Winkelmann, Chen Zhao, Tong Zhang, and Thomas Zwinger
The Cryosphere, 17, 5197–5217, https://doi.org/10.5194/tc-17-5197-2023, https://doi.org/10.5194/tc-17-5197-2023, 2023
Short summary
Short summary
Mass loss from Antarctica is a key contributor to sea level rise over the 21st century, and the associated uncertainty dominates sea level projections. We highlight here the Antarctic glaciers showing the largest changes and quantify the main sources of uncertainty in their future evolution using an ensemble of ice flow models. We show that on top of Pine Island and Thwaites glaciers, Totten and Moscow University glaciers show rapid changes and a strong sensitivity to warmer ocean conditions.
Stefania A. Ciliberti, Enrique Alvarez Fanjul, Jay Pearlman, Kirsten Wilmer-Becker, Pierre Bahurel, Fabrice Ardhuin, Alain Arnaud, Mike Bell, Segolene Berthou, Laurent Bertino, Arthur Capet, Eric Chassignet, Stefano Ciavatta, Mauro Cirano, Emanuela Clementi, Gianpiero Cossarini, Gianpaolo Coro, Stuart Corney, Fraser Davidson, Marie Drevillon, Yann Drillet, Renaud Dussurget, Ghada El Serafy, Katja Fennel, Marcos Garcia Sotillo, Patrick Heimbach, Fabrice Hernandez, Patrick Hogan, Ibrahim Hoteit, Sudheer Joseph, Simon Josey, Pierre-Yves Le Traon, Simone Libralato, Marco Mancini, Pascal Matte, Angelique Melet, Yasumasa Miyazawa, Andrew M. Moore, Antonio Novellino, Andrew Porter, Heather Regan, Laia Romero, Andreas Schiller, John Siddorn, Joanna Staneva, Cecile Thomas-Courcoux, Marina Tonani, Jose Maria Garcia-Valdecasas, Jennifer Veitch, Karina von Schuckmann, Liying Wan, John Wilkin, and Romane Zufic
State Planet, 1-osr7, 2, https://doi.org/10.5194/sp-1-osr7-2-2023, https://doi.org/10.5194/sp-1-osr7-2-2023, 2023
Carl Wunsch, Sarah Williamson, and Patrick Heimbach
Ocean Sci., 19, 1253–1275, https://doi.org/10.5194/os-19-1253-2023, https://doi.org/10.5194/os-19-1253-2023, 2023
Short summary
Short summary
Data assimilation methods that couple observations with dynamical models are essential for understanding climate change. Here,
climateincludes all sub-elements (ocean, atmosphere, ice, etc.). A common form of combination arises from sequential estimation theory, a methodology susceptible to a variety of errors that can accumulate through time for long records. Using two simple analogs, examples of these errors are identified and discussed, along with suggestions for accommodating them.
David S. Trossman, Caitlin B. Whalen, Thomas W. N. Haine, Amy F. Waterhouse, An T. Nguyen, Arash Bigdeli, Matthew Mazloff, and Patrick Heimbach
Ocean Sci., 18, 729–759, https://doi.org/10.5194/os-18-729-2022, https://doi.org/10.5194/os-18-729-2022, 2022
Short summary
Short summary
How the ocean mixes is not yet adequately represented by models. There are many challenges with representing this mixing. A model that minimizes disagreements between observations and the model could be used to fill in the gaps from observations to better represent ocean mixing. But observations of ocean mixing have large uncertainties. Here, we show that ocean oxygen, which has relatively small uncertainties, and observations of ocean mixing provide information similar to the model.
Amy Solomon, Céline Heuzé, Benjamin Rabe, Sheldon Bacon, Laurent Bertino, Patrick Heimbach, Jun Inoue, Doroteaciro Iovino, Ruth Mottram, Xiangdong Zhang, Yevgeny Aksenov, Ronan McAdam, An Nguyen, Roshin P. Raj, and Han Tang
Ocean Sci., 17, 1081–1102, https://doi.org/10.5194/os-17-1081-2021, https://doi.org/10.5194/os-17-1081-2021, 2021
Short summary
Short summary
Freshwater in the Arctic Ocean plays a critical role in the global climate system by impacting ocean circulations, stratification, mixing, and emergent regimes. In this review paper we assess how Arctic Ocean freshwater changed in the 2010s relative to the 2000s. Estimates from observations and reanalyses show a qualitative stabilization in the 2010s due to a compensation between a freshening of the Beaufort Gyre and a reduction in freshwater in the Amerasian and Eurasian basins.
Matthias Scheiter, Marius Schaefer, Eduardo Flández, Deniz Bozkurt, and Ralf Greve
The Cryosphere, 15, 3637–3654, https://doi.org/10.5194/tc-15-3637-2021, https://doi.org/10.5194/tc-15-3637-2021, 2021
Short summary
Short summary
We simulate the current state and future evolution of the Mocho-Choshuenco ice cap in southern Chile (40°S, 72°W) with the ice-sheet model SICOPOLIS. Under different global warming scenarios, we project ice mass losses between 56 % and 97 % by the end of the 21st century. We quantify the uncertainties based on an ensemble of climate models and on the temperature dependence of the equilibrium line altitude. Our results suggest a considerable deglaciation in southern Chile in the next 80 years.
Christopher Chambers, Ralf Greve, Bas Altena, and Pierre-Marie Lefeuvre
The Cryosphere, 14, 3747–3759, https://doi.org/10.5194/tc-14-3747-2020, https://doi.org/10.5194/tc-14-3747-2020, 2020
Short summary
Short summary
The topography of the rock below the Greenland ice sheet is not well known. One long valley appears as a line of dips because of incomplete data. So we use ice model simulations that unblock this valley, and these create a watercourse that may represent a form of river over 1000 km long under the ice. When we melt ice at the bottom of the ice sheet only in the deep interior, water can flow down the valley only when the valley is unblocked. It may have developed while an ice sheet was present.
Heiko Goelzer, Sophie Nowicki, Anthony Payne, Eric Larour, Helene Seroussi, William H. Lipscomb, Jonathan Gregory, Ayako Abe-Ouchi, Andrew Shepherd, Erika Simon, Cécile Agosta, Patrick Alexander, Andy Aschwanden, Alice Barthel, Reinhard Calov, Christopher Chambers, Youngmin Choi, Joshua Cuzzone, Christophe Dumas, Tamsin Edwards, Denis Felikson, Xavier Fettweis, Nicholas R. Golledge, Ralf Greve, Angelika Humbert, Philippe Huybrechts, Sebastien Le clec'h, Victoria Lee, Gunter Leguy, Chris Little, Daniel P. Lowry, Mathieu Morlighem, Isabel Nias, Aurelien Quiquet, Martin Rückamp, Nicole-Jeanne Schlegel, Donald A. Slater, Robin S. Smith, Fiamma Straneo, Lev Tarasov, Roderik van de Wal, and Michiel van den Broeke
The Cryosphere, 14, 3071–3096, https://doi.org/10.5194/tc-14-3071-2020, https://doi.org/10.5194/tc-14-3071-2020, 2020
Short summary
Short summary
In this paper we use a large ensemble of Greenland ice sheet models forced by six different global climate models to project ice sheet changes and sea-level rise contributions over the 21st century.
The results for two different greenhouse gas concentration scenarios indicate that the Greenland ice sheet will continue to lose mass until 2100, with contributions to sea-level rise of 90 ± 50 mm and 32 ± 17 mm for the high (RCP8.5) and low (RCP2.6) scenario, respectively.
Hélène Seroussi, Sophie Nowicki, Antony J. Payne, Heiko Goelzer, William H. Lipscomb, Ayako Abe-Ouchi, Cécile Agosta, Torsten Albrecht, Xylar Asay-Davis, Alice Barthel, Reinhard Calov, Richard Cullather, Christophe Dumas, Benjamin K. Galton-Fenzi, Rupert Gladstone, Nicholas R. Golledge, Jonathan M. Gregory, Ralf Greve, Tore Hattermann, Matthew J. Hoffman, Angelika Humbert, Philippe Huybrechts, Nicolas C. Jourdain, Thomas Kleiner, Eric Larour, Gunter R. Leguy, Daniel P. Lowry, Chistopher M. Little, Mathieu Morlighem, Frank Pattyn, Tyler Pelle, Stephen F. Price, Aurélien Quiquet, Ronja Reese, Nicole-Jeanne Schlegel, Andrew Shepherd, Erika Simon, Robin S. Smith, Fiammetta Straneo, Sainan Sun, Luke D. Trusel, Jonas Van Breedam, Roderik S. W. van de Wal, Ricarda Winkelmann, Chen Zhao, Tong Zhang, and Thomas Zwinger
The Cryosphere, 14, 3033–3070, https://doi.org/10.5194/tc-14-3033-2020, https://doi.org/10.5194/tc-14-3033-2020, 2020
Short summary
Short summary
The Antarctic ice sheet has been losing mass over at least the past 3 decades in response to changes in atmospheric and oceanic conditions. This study presents an ensemble of model simulations of the Antarctic evolution over the 2015–2100 period based on various ice sheet models, climate forcings and emission scenarios. Results suggest that the West Antarctic ice sheet will continue losing a large amount of ice, while the East Antarctic ice sheet could experience increased snow accumulation.
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.
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.
Reinhard Calov, Sebastian Beyer, Ralf Greve, Johanna Beckmann, Matteo Willeit, Thomas Kleiner, Martin Rückamp, Angelika Humbert, and Andrey Ganopolski
The Cryosphere, 12, 3097–3121, https://doi.org/10.5194/tc-12-3097-2018, https://doi.org/10.5194/tc-12-3097-2018, 2018
Short summary
Short summary
We present RCP 4.5 and 8.5 projections for the Greenland glacial system with the new glacial system model IGLOO 1.0, which incorporates the ice sheet model SICOPOLIS 3.3, a model of basal hydrology and a parameterization of submarine melt of outlet glaciers. Surface temperature and mass balance anomalies from the MAR climate model serve as forcing delivering projections for the contribution of the Greenland ice sheet to sea level rise and submarine melt of Helheim and Store outlet glaciers.
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.
Nat Wilson, Fiammetta Straneo, and Patrick Heimbach
The Cryosphere, 11, 2773–2782, https://doi.org/10.5194/tc-11-2773-2017, https://doi.org/10.5194/tc-11-2773-2017, 2017
Short summary
Short summary
We estimate submarine melt rates from ice tongues in northern Greenland using WorldView satellite imagery. At Ryder Glacier, melt is strongly concentrated around regions where subglacier channels likely enter the fjord. At the 79 North Glacier, we find a large volume imbalance in which melting removes a greater quantity of ice than is replaced by inflow over the grounding line. This leads us to suggest that a reduction in the spatial extent of the ice tongue is possible over the coming decade.
Hakime Seddik, Ralf Greve, Thomas Zwinger, and Shin Sugiyama
The Cryosphere, 11, 2213–2229, https://doi.org/10.5194/tc-11-2213-2017, https://doi.org/10.5194/tc-11-2213-2017, 2017
Short summary
Short summary
The Shirase Glacier in Antarctica is studied by means of a computer model. This model implements two physical approaches to represent the glacier flow dynamics. This study finds that it is important to use the more precise and sophisticated method in order to better understand and predict the evolution of fast flowing glaciers. This may be important to more accurately predict the sea level change due to global warming.
Rupert Michael Gladstone, Roland Charles Warner, Benjamin Keith Galton-Fenzi, Olivier Gagliardini, Thomas Zwinger, and Ralf Greve
The Cryosphere, 11, 319–329, https://doi.org/10.5194/tc-11-319-2017, https://doi.org/10.5194/tc-11-319-2017, 2017
Short summary
Short summary
Computer models are used to simulate the behaviour of glaciers and ice sheets. It has been found that such models are required to be run at very high resolution (which means high computational expense) in order to accurately represent the evolution of marine ice sheets (ice sheets resting on bedrock below sea level), in certain situations which depend on sub-glacial physical processes.
Jorge Bernales, Irina Rogozhina, Ralf Greve, and Maik Thomas
The Cryosphere, 11, 247–265, https://doi.org/10.5194/tc-11-247-2017, https://doi.org/10.5194/tc-11-247-2017, 2017
Short summary
Short summary
This study offers a hard test to the models commonly used to simulate an ice sheet evolution over multimillenial timescales. Using an example of the Antarctic Ice Sheet, we evaluate the performance of such models against observations and highlight a strong impact of different approaches towards modeling rapidly flowing ice sectors. In particular, our results show that inferences of previous studies may need significant adjustments to be adopted by a different type of ice sheet models.
Liz C. Logan, Luc L. Lavier, Eunseo Choi, Eh Tan, and Ginny A. Catania
The Cryosphere, 11, 117–132, https://doi.org/10.5194/tc-11-117-2017, https://doi.org/10.5194/tc-11-117-2017, 2017
Short summary
Short summary
Global sea level rise prediction is a pressing and unresolved problem, one whose solution depends upon glaciologists better predicting ice sheet shrinkage due to iceberg calving. We present a numerical model that is capable of simulating ice flow and breakage that leads to iceberg calving and find that a material property that captures both the fluid- and solid-like behavior of ice simultaneously is a necessary condition for studying areas of glaciers in contact with ocean water prone to calve.
Daniel N. Goldberg, Sri Hari Krishna Narayanan, Laurent Hascoet, and Jean Utke
Geosci. Model Dev., 9, 1891–1904, https://doi.org/10.5194/gmd-9-1891-2016, https://doi.org/10.5194/gmd-9-1891-2016, 2016
Short summary
Short summary
Geophysical adjoint models are powerful tools, allowing sensitivity studies that are not possible otherwise, and enabling optimized fit of models to observing data sets. The complexity involved requires the use of algorithmic differentiation (AD) software, but AD adjoint calculation for ice models can be slow, with prohibitive memory requirements. In this paper, we present a method to improve the performance of ice model adjoint generation, in terms of timing, memory load, and accuracy.
T. Goelles, C. E. Bøggild, and R. Greve
The Cryosphere, 9, 1845–1856, https://doi.org/10.5194/tc-9-1845-2015, https://doi.org/10.5194/tc-9-1845-2015, 2015
Short summary
Short summary
Soot (black carbon) and dust particles darken the surface of ice sheets and glaciers as they accumulate. This causes more ice to melt, which releases more particles from within the ice. This positive feedback mechanism is studied with a new two-dimensional model, mimicking the conditions of Greenland, under different climate warming scenarios. In the warmest scenario, the additional ice sheet mass loss until the year 3000 is up to 7%.
A. Levermann, R. Winkelmann, S. Nowicki, J. L. Fastook, K. Frieler, R. Greve, H. H. Hellmer, M. A. Martin, M. Meinshausen, M. Mengel, A. J. Payne, D. Pollard, T. Sato, R. Timmermann, W. L. Wang, and R. A. Bindschadler
Earth Syst. Dynam., 5, 271–293, https://doi.org/10.5194/esd-5-271-2014, https://doi.org/10.5194/esd-5-271-2014, 2014
T. Sato, T. Shiraiwa, R. Greve, H. Seddik, E. Edelmann, and T. Zwinger
Clim. Past, 10, 393–404, https://doi.org/10.5194/cp-10-393-2014, https://doi.org/10.5194/cp-10-393-2014, 2014
O. Gagliardini, T. Zwinger, F. Gillet-Chaulet, G. Durand, L. Favier, B. de Fleurian, R. Greve, M. Malinen, C. Martín, P. Råback, J. Ruokolainen, M. Sacchettini, M. Schäfer, H. Seddik, and J. Thies
Geosci. Model Dev., 6, 1299–1318, https://doi.org/10.5194/gmd-6-1299-2013, https://doi.org/10.5194/gmd-6-1299-2013, 2013
F. Gillet-Chaulet, O. Gagliardini, H. Seddik, M. Nodet, G. Durand, C. Ritz, T. Zwinger, R. Greve, and D. G. Vaughan
The Cryosphere, 6, 1561–1576, https://doi.org/10.5194/tc-6-1561-2012, https://doi.org/10.5194/tc-6-1561-2012, 2012
Related subject area
Cryosphere
Tuning parameters of a sea ice model using machine learning
WRF-Chem simulations of snow nitrate and other physicochemical properties in northern China
Clustering simulated snow profiles to form avalanche forecast regions
SnowQM 1.0: a fast R package for bias-correcting spatial fields of snow water equivalent using quantile mapping
Simulation of snow albedo and solar irradiance profile with the Two-streAm Radiative TransfEr in Snow (TARTES) v2.0 model
Evaluation of MITgcm-based ocean reanalyses for the Southern Ocean
Improvements in the land surface configuration to better simulate seasonal snow cover in the European Alps with the CNRM-AROME (cycle 46) convection-permitting regional climate model
A three-stage model pipeline predicting regional avalanche danger in Switzerland (RAvaFcast v1.0.0): a decision-support tool for operational avalanche forecasting
A global–land snow scheme (GLASS) v1.0 for the GFDL Earth System Model: formulation and evaluation at instrumented sites
A gradient-boosted tree framework to model the ice thickness of the World's glaciers (IceBoost v1)
Design and performance of ELSA v2.0: an isochronal model for ice-sheet layer tracing
Southern Ocean Ice Prediction System version 1.0 (SOIPS v1.0): description of the system and evaluation of synoptic-scale sea ice forecasts
Lagrangian tracking of sea ice in Community Ice CodE (CICE; version 5)
openAMUNDSEN v1.0: an open-source snow-hydrological model for mountain regions
OpenFOAM-avalanche 2312: depth-integrated models beyond dense-flow avalanches
Refactoring the elastic–viscous–plastic solver from the sea ice model CICE v6.5.1 for improved performance
CMIP6 models overestimate sea ice melt, growth & conduction relative to ice mass balance buoy estimates
A new 3D full-Stokes calving algorithm within Elmer/Ice (v9.0)
Towards deep learning solutions for classification of automated snow height measurements (CleanSnow v1.0.0)
Quantitative Sub-Ice and Marine Tracing of Antarctic Sediment Provenance (TASP v1.0)
Coupling framework (1.0) for the Úa (2023b) ice sheet model and the FESOM-1.4 z-coordinate ocean model in an Antarctic domain
A novel numerical implementation for the surface energy budget of melting snowpacks and glaciers
SnowPappus v1.0, a blowing-snow model for large-scale applications of the Crocus snow scheme
A stochastic parameterization of ice sheet surface mass balance for the Stochastic Ice-Sheet and Sea-Level System Model (StISSM v1.0)
Graphics-processing-unit-accelerated ice flow solver for unstructured meshes using the Shallow-Shelf Approximation (FastIceFlo v1.0.1)
A finite-element framework to explore the numerical solution of the coupled problem of heat conduction, water vapor diffusion, and settlement in dry snow (IvoriFEM v0.1.0)
AvaFrame com1DFA (v1.3): a thickness-integrated computational avalanche module – theory, numerics, and testing
Universal differential equations for glacier ice flow modelling
A new model for supraglacial hydrology evolution and drainage for the Greenland Ice Sheet (SHED v1.0)
Modeling sensitivities of thermally and hydraulically driven ice stream surge cycling
A parallel implementation of the confined–unconfined aquifer system model for subglacial hydrology: design, verification, and performance analysis (CUAS-MPI v0.1.0)
Automatic snow type classification of snow micropenetrometer profiles with machine learning algorithms
An empirical model to calculate snow depth from daily snow water equivalent: SWE2HS 1.0
A wind-driven snow redistribution module for Alpine3D v3.3.0: adaptations designed for downscaling ice sheet surface mass balance
The CryoGrid community model (version 1.0) – a multi-physics toolbox for climate-driven simulations in the terrestrial cryosphere
Glacier Energy and Mass Balance (GEMB): a model of firn processes for cryosphere research
Sensitivity of NEMO4.0-SI3 model parameters on sea ice budgets in the Southern Ocean
Introducing CRYOWRF v1.0: multiscale atmospheric flow simulations with advanced snow cover modelling
SUHMO: an adaptive mesh refinement SUbglacial Hydrology MOdel v1.0
Improving snow albedo modeling in the E3SM land model (version 2.0) and assessing its impacts on snow and surface fluxes over the Tibetan Plateau
The Multiple Snow Data Assimilation System (MuSA v1.0)
The Stochastic Ice-Sheet and Sea-Level System Model v1.0 (StISSM v1.0)
Improved representation of the contemporary Greenland ice sheet firn layer by IMAU-FDM v1.2G
Modeling the small-scale deposition of snow onto structured Arctic sea ice during a MOSAiC storm using snowBedFoam 1.0.
Benchmarking the vertically integrated ice-sheet model IMAU-ICE (version 2.0)
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)
Anton Korosov, Yue Ying, and Einar Ólason
Geosci. Model Dev., 18, 885–904, https://doi.org/10.5194/gmd-18-885-2025, https://doi.org/10.5194/gmd-18-885-2025, 2025
Short summary
Short summary
We have developed a new method to improve the accuracy of sea ice models, which predict how ice moves and deforms due to wind and ocean currents. Traditional models use parameters that are often poorly defined. The new approach uses machine learning to fine-tune these parameters by comparing simulated ice drift with satellite data. The method identifies optimal settings for the model by analysing patterns in ice deformation. This results in more accurate simulations of sea ice drift forecasting.
Xia Wang, Tao Che, Xueyin Ruan, Shanna Yue, Jing Wang, Chun Zhao, and Lei Geng
Geosci. Model Dev., 18, 651–670, https://doi.org/10.5194/gmd-18-651-2025, https://doi.org/10.5194/gmd-18-651-2025, 2025
Short summary
Short summary
We employed the WRF-Chem model to parameterize atmospheric nitrate deposition in snow and evaluate its performance in simulating snow cover, snow depth, and concentrations of dust and nitrate using new observations from northern China. The results generally exhibit reasonable agreement with field observations in northern China, demonstrating the model's capability to simulate snow properties, including concentrations of reservoir species.
Simon Horton, Florian Herla, and Pascal Haegeli
Geosci. Model Dev., 18, 193–209, https://doi.org/10.5194/gmd-18-193-2025, https://doi.org/10.5194/gmd-18-193-2025, 2025
Short summary
Short summary
We present a method for avalanche forecasters to analyze patterns in snowpack model simulations. It uses fuzzy clustering to group small regions into larger forecast areas based on snow characteristics, locations, and temporal history. Tested in the Columbia Mountains in two winter seasons, it closely matched real forecast regions regions and identified major avalanche hazard patterns. This approach simplifies complex model outputs, helping forecasters make informed decisions.
Adrien Michel, Johannes Aschauer, Tobias Jonas, Stefanie Gubler, Sven Kotlarski, and Christoph Marty
Geosci. Model Dev., 17, 8969–8988, https://doi.org/10.5194/gmd-17-8969-2024, https://doi.org/10.5194/gmd-17-8969-2024, 2024
Short summary
Short summary
We present a method to correct snow cover maps (represented in terms of snow water equivalent) to match better-quality maps. The correction can then be extended backwards and forwards in time for periods when better-quality maps are not available. The method is fast and gives good results. It is then applied to obtain a climatology of the snow cover in Switzerland over the past 60 years at a resolution of 1 d and 1 km. This is the first time that such a dataset has been produced.
Ghislain Picard and Quentin Libois
Geosci. Model Dev., 17, 8927–8953, https://doi.org/10.5194/gmd-17-8927-2024, https://doi.org/10.5194/gmd-17-8927-2024, 2024
Short summary
Short summary
The Two-streAm Radiative TransfEr in Snow (TARTES) is a radiative transfer model to compute snow albedo in the solar domain and the profiles of light and energy absorption in a multi-layered snowpack whose physical properties are user defined. It uniquely considers snow grain shape flexibly, based on recent insights showing that snow does not behave as a collection of ice spheres but instead as a random medium. TARTES is user-friendly yet performs comparably to more complex models.
Yoshihiro Nakayama, Alena Malyarenko, Hong Zhang, Ou Wang, Matthis Auger, Yafei Nie, Ian Fenty, Matthew Mazloff, Armin Köhl, and Dimitris Menemenlis
Geosci. Model Dev., 17, 8613–8638, https://doi.org/10.5194/gmd-17-8613-2024, https://doi.org/10.5194/gmd-17-8613-2024, 2024
Short summary
Short summary
Global- and basin-scale ocean reanalyses are becoming easily accessible. However, such ocean reanalyses are optimized for their entire model domains and their ability to simulate the Southern Ocean requires evaluation. We conduct intercomparison analyses of Massachusetts Institute of Technology General Circulation Model (MITgcm)-based ocean reanalyses. They generally perform well for the open ocean, but open-ocean temporal variability and Antarctic continental shelves require improvements.
Diego Monteiro, Cécile Caillaud, Matthieu Lafaysse, Adrien Napoly, Mathieu Fructus, Antoinette Alias, and Samuel Morin
Geosci. Model Dev., 17, 7645–7677, https://doi.org/10.5194/gmd-17-7645-2024, https://doi.org/10.5194/gmd-17-7645-2024, 2024
Short summary
Short summary
Modeling snow cover in climate and weather forecasting models is a challenge even for high-resolution models. Recent simulations with CNRM-AROME have shown difficulties when representing snow in the European Alps. Using remote sensing data and in situ observations, we evaluate modifications of the land surface configuration in order to improve it. We propose a new surface configuration, enabling a more realistic simulation of snow cover, relevant for climate and weather forecasting applications.
Alessandro Maissen, Frank Techel, and Michele Volpi
Geosci. Model Dev., 17, 7569–7593, https://doi.org/10.5194/gmd-17-7569-2024, https://doi.org/10.5194/gmd-17-7569-2024, 2024
Short summary
Short summary
By harnessing AI models, this work enables processing large amounts of data, including weather conditions, snowpack characteristics, and historical avalanche data, to predict human-like avalanche forecasts in Switzerland. Our proposed model can significantly assist avalanche forecasters in their decision-making process, thereby facilitating more efficient and accurate predictions crucial for ensuring safety in Switzerland's avalanche-prone regions.
Enrico Zorzetto, Sergey Malyshev, Paul Ginoux, and Elena Shevliakova
Geosci. Model Dev., 17, 7219–7244, https://doi.org/10.5194/gmd-17-7219-2024, https://doi.org/10.5194/gmd-17-7219-2024, 2024
Short summary
Short summary
We describe a new snow scheme developed for use in global climate models, which simulates the interactions of snowpack with vegetation, atmosphere, and soil. We test the new snow model over a set of sites where in situ observations are available. We find that when compared to a simpler snow model, this model improves predictions of seasonal snow and of soil temperature under the snowpack, important variables for simulating both the hydrological cycle and the global climate system.
Niccolò Maffezzoli, Eric Rignot, Carlo Barbante, Troels Petersen, and Sebastiano Vascon
EGUsphere, https://doi.org/10.5194/egusphere-2024-2455, https://doi.org/10.5194/egusphere-2024-2455, 2024
Short summary
Short summary
In this work we introduces IceBoost, a machine learning framework to model the ice thickness distribution of all the World's glaciers, with greater accuracy than state-of-the-art methods. The model is trained on 3.7 million measurements globally available and provides skillful estimates across all regions. This advancement will help in better assessing future sea level changes, freshwater resources, with significance for both the scientific community and society at large.
Therese Rieckh, Andreas Born, Alexander Robinson, Robert Law, and Gerrit Gülle
Geosci. Model Dev., 17, 6987–7000, https://doi.org/10.5194/gmd-17-6987-2024, https://doi.org/10.5194/gmd-17-6987-2024, 2024
Short summary
Short summary
We present the open-source model ELSA, which simulates the internal age structure of large ice sheets. It creates layers of snow accumulation at fixed times during the simulation, which are used to model the internal stratification of the ice sheet. Together with reconstructed isochrones from radiostratigraphy data, ELSA can be used to assess ice sheet models and to improve their parameterization. ELSA can be used coupled to an ice sheet model or forced with its output.
Fu Zhao, Xi Liang, Zhongxiang Tian, Ming Li, Na Liu, and Chengyan Liu
Geosci. Model Dev., 17, 6867–6886, https://doi.org/10.5194/gmd-17-6867-2024, https://doi.org/10.5194/gmd-17-6867-2024, 2024
Short summary
Short summary
In this work, we introduce a newly developed Antarctic sea ice forecasting system, namely the Southern Ocean Ice Prediction System (SOIPS). The system is based on a regional sea ice‒ocean‒ice shelf coupled model and can assimilate sea ice concentration observations. By assessing the system's performance in sea ice forecasts, we find that the system can provide reliable Antarctic sea ice forecasts for the next 7 d and has the potential to guide ship navigation in the Antarctic sea ice zone.
Chenhui Ning, Shiming Xu, Yan Zhang, Xuantong Wang, Zhihao Fan, and Jiping Liu
Geosci. Model Dev., 17, 6847–6866, https://doi.org/10.5194/gmd-17-6847-2024, https://doi.org/10.5194/gmd-17-6847-2024, 2024
Short summary
Short summary
Sea ice models are mainly based on non-moving structured grids, which is different from buoy measurements that follow the ice drift. To facilitate Lagrangian analysis, we introduce online tracking of sea ice in Community Ice CodE (CICE). We validate the sea ice tracking with buoys and evaluate the sea ice deformation in high-resolution simulations, which show multi-fractal characteristics. The source code is openly available and can be used in various scientific and operational applications.
Ulrich Strasser, Michael Warscher, Erwin Rottler, and Florian Hanzer
Geosci. Model Dev., 17, 6775–6797, https://doi.org/10.5194/gmd-17-6775-2024, https://doi.org/10.5194/gmd-17-6775-2024, 2024
Short summary
Short summary
openAMUNDSEN is a fully distributed open-source snow-hydrological model for mountain catchments. It includes process representations of an empirical, semi-empirical, and physical nature. It uses temperature, precipitation, humidity, radiation, and wind speed as forcing data and is computationally efficient, of a modular nature, and easily extendible. The Python code is available on GitHub (https://github.com/openamundsen/openamundsen), including documentation (https://doc.openamundsen.org).
Matthias Rauter and Julia Kowalski
Geosci. Model Dev., 17, 6545–6569, https://doi.org/10.5194/gmd-17-6545-2024, https://doi.org/10.5194/gmd-17-6545-2024, 2024
Short summary
Short summary
Snow avalanches can form large powder clouds that substantially exceed the velocity and reach of the dense core. Only a few complex models exist to simulate this phenomenon, and the respective hazard is hard to predict. This work provides a novel flow model that focuses on simple relations while still encapsulating the significant behaviour. The model is applied to reconstruct two catastrophic powder snow avalanche events in Austria.
Till Andreas Soya Rasmussen, Jacob Poulsen, Mads Hvid Ribergaard, Ruchira Sasanka, Anthony P. Craig, Elizabeth C. Hunke, and Stefan Rethmeier
Geosci. Model Dev., 17, 6529–6544, https://doi.org/10.5194/gmd-17-6529-2024, https://doi.org/10.5194/gmd-17-6529-2024, 2024
Short summary
Short summary
Earth system models (ESMs) today strive for better quality based on improved resolutions and improved physics. A limiting factor is the supercomputers at hand and how best to utilize them. This study focuses on the refactorization of one part of a sea ice model (CICE), namely the dynamics. It shows that the performance can be significantly improved, which means that one can either run the same simulations much cheaper or advance the system according to what is needed.
Alex Edward West and Edward William Blockley
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-121, https://doi.org/10.5194/gmd-2024-121, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
This study uses ice mass balance buoys – temperature and height-measuring devices frozen into sea ice – to find how well climate models simulate the melt & growth of, and conduction of heat through, Arctic sea ice. This may help understand why models produce varying amounts of sea ice in the present day. We find models tend to show more melt, growth or conduction for a given ice thickness than the buoys, though the difference is smaller for models with more physically realistic thermodynamics.
Iain Wheel, Douglas I. Benn, Anna J. Crawford, Joe Todd, and Thomas Zwinger
Geosci. Model Dev., 17, 5759–5777, https://doi.org/10.5194/gmd-17-5759-2024, https://doi.org/10.5194/gmd-17-5759-2024, 2024
Short summary
Short summary
Calving, the detachment of large icebergs from glaciers, is one of the largest uncertainties in future sea level rise projections. This process is poorly understood, and there is an absence of detailed models capable of simulating calving. A new 3D calving model has been developed to better understand calving at glaciers where detailed modelling was previously limited. Importantly, the new model is very flexible. By allowing for unrestricted calving geometries, it can be applied at any location.
Jan Svoboda, Marc Ruesch, David Liechti, Corinne Jones, Michele Volpi, Michael Zehnder, and Jürg Schweizer
EGUsphere, https://doi.org/10.5194/egusphere-2024-1752, https://doi.org/10.5194/egusphere-2024-1752, 2024
Short summary
Short summary
Accurately measuring snow height is key for modeling approaches in climate sciences, snow hydrology and avalanche forecasting. Erroneous snow height measurements often occur when the snow height is low or changes, for instance, during a snowfall in the summer. We prepare a new benchmark dataset with annotated snow height data and demonstrate how to improve the measurement quality using modern deep learning approaches. Our approach can be easily implemented into a data pipeline for snow modeling.
Jim Marschalek, Edward Gasson, Tina van de Flierdt, Claus-Dieter Hillenbrand, Martin Siegert, and Liam Holder
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-104, https://doi.org/10.5194/gmd-2024-104, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
Ice sheet models can help predict how Antarctica's ice sheets respond to environmental change, and such models benefit from comparison to geological data. Here, we use an ice sheet model output, plus other data, to predict the erosion of debris and trace its transport to where it is deposited on the ocean floor. This allows the results of ice sheet modelling to be directly and quantitively compared to real-world data, helping to reduce uncertainty regarding Antarctic sea level contribution.
Ole Richter, Ralph Timmermann, G. Hilmar Gudmundsson, and Jan De Rydt
EGUsphere, https://doi.org/10.5194/egusphere-2024-648, https://doi.org/10.5194/egusphere-2024-648, 2024
Short summary
Short summary
The new coupled ice sheet-ocean model addresses challenges related to horizontal resolution through advanced mesh flexibility, enabled by the use of unstructured grids. We describe the new model, verify its functioning in an idealised setting and demonstrate its advantages in a global-ocean/Antarctic ice sheet domain. The results of this study comprise an important step towards improving predictions of the Antarctic contribution to sea level rise over centennial time scales.
Kévin Fourteau, Julien Brondex, Fanny Brun, and Marie Dumont
Geosci. Model Dev., 17, 1903–1929, https://doi.org/10.5194/gmd-17-1903-2024, https://doi.org/10.5194/gmd-17-1903-2024, 2024
Short summary
Short summary
In this paper, we provide a novel numerical implementation for solving the energy exchanges at the surface of snow and ice. By combining the strong points of previous models, our solution leads to more accurate and robust simulations of the energy exchanges, surface temperature, and melt while preserving a reasonable computation time.
Matthieu Baron, Ange Haddjeri, Matthieu Lafaysse, Louis Le Toumelin, Vincent Vionnet, and Mathieu Fructus
Geosci. Model Dev., 17, 1297–1326, https://doi.org/10.5194/gmd-17-1297-2024, https://doi.org/10.5194/gmd-17-1297-2024, 2024
Short summary
Short summary
Increasing the spatial resolution of numerical systems simulating snowpack evolution in mountain areas requires representing small-scale processes such as wind-induced snow transport. We present SnowPappus, a simple scheme coupled with the Crocus snow model to compute blowing-snow fluxes and redistribute snow among grid points at 250 m resolution. In terms of numerical cost, it is suitable for large-scale applications. We present point-scale evaluations of fluxes and snow transport occurrence.
Lizz Ultee, Alexander A. Robel, and Stefano Castruccio
Geosci. Model Dev., 17, 1041–1057, https://doi.org/10.5194/gmd-17-1041-2024, https://doi.org/10.5194/gmd-17-1041-2024, 2024
Short summary
Short summary
The surface mass balance (SMB) of an ice sheet describes the net gain or loss of mass from ice sheets (such as those in Greenland and Antarctica) through interaction with the atmosphere. We developed a statistical method to generate a wide range of SMB fields that reflect the best understanding of SMB processes. Efficiently sampling the variability of SMB will help us understand sources of uncertainty in ice sheet model projections.
Anjali Sandip, Ludovic Räss, and Mathieu Morlighem
Geosci. Model Dev., 17, 899–909, https://doi.org/10.5194/gmd-17-899-2024, https://doi.org/10.5194/gmd-17-899-2024, 2024
Short summary
Short summary
We solve momentum balance for unstructured meshes to predict ice flow for real glaciers using a pseudo-transient method on graphics processing units (GPUs) and compare it to a standard central processing unit (CPU) implementation. We justify the GPU implementation by applying the price-to-performance metric for up to million-grid-point spatial resolutions. This study represents a first step toward leveraging GPU processing power, enabling more accurate polar ice discharge predictions.
Julien Brondex, Kévin Fourteau, Marie Dumont, Pascal Hagenmuller, Neige Calonne, François Tuzet, and Henning Löwe
Geosci. Model Dev., 16, 7075–7106, https://doi.org/10.5194/gmd-16-7075-2023, https://doi.org/10.5194/gmd-16-7075-2023, 2023
Short summary
Short summary
Vapor diffusion is one of the main processes governing snowpack evolution, and it must be accounted for in models. Recent attempts to represent vapor diffusion in numerical models have faced several difficulties regarding computational cost and mass and energy conservation. Here, we develop our own finite-element software to explore numerical approaches and enable us to overcome these difficulties. We illustrate the capability of these approaches on established numerical benchmarks.
Matthias Tonnel, Anna Wirbel, Felix Oesterle, and Jan-Thomas Fischer
Geosci. Model Dev., 16, 7013–7035, https://doi.org/10.5194/gmd-16-7013-2023, https://doi.org/10.5194/gmd-16-7013-2023, 2023
Short summary
Short summary
Avaframe - the open avalanche framework - provides open-source tools to simulate and investigate snow avalanches. It is utilized for multiple purposes, the two main applications being hazard mapping and scientific research of snow processes. We present the theory, conversion to a computer model, and testing for one of the core modules used for simulations of a particular type of avalanche, the so-called dense-flow avalanches. Tests check and confirm the applicability of the utilized method.
Jordi Bolibar, Facundo Sapienza, Fabien Maussion, Redouane Lguensat, Bert Wouters, and Fernando Pérez
Geosci. Model Dev., 16, 6671–6687, https://doi.org/10.5194/gmd-16-6671-2023, https://doi.org/10.5194/gmd-16-6671-2023, 2023
Short summary
Short summary
We developed a new modelling framework combining numerical methods with machine learning. Using this approach, we focused on understanding how ice moves within glaciers, and we successfully learnt a prescribed law describing ice movement for 17 glaciers worldwide as a proof of concept. Our framework has the potential to discover important laws governing glacier processes, aiding our understanding of glacier physics and their contribution to water resources and sea-level rise.
Prateek Gantayat, Alison F. Banwell, Amber A. Leeson, James M. Lea, Dorthe Petersen, Noel Gourmelen, and Xavier Fettweis
Geosci. Model Dev., 16, 5803–5823, https://doi.org/10.5194/gmd-16-5803-2023, https://doi.org/10.5194/gmd-16-5803-2023, 2023
Short summary
Short summary
We developed a new supraglacial hydrology model for the Greenland Ice Sheet. This model simulates surface meltwater routing, meltwater drainage, supraglacial lake (SGL) overflow, and formation of lake ice. The model was able to reproduce 80 % of observed lake locations and provides a good match between the observed and modelled temporal evolution of SGLs.
Kevin Hank, Lev Tarasov, and Elisa Mantelli
Geosci. Model Dev., 16, 5627–5652, https://doi.org/10.5194/gmd-16-5627-2023, https://doi.org/10.5194/gmd-16-5627-2023, 2023
Short summary
Short summary
Physically meaningful modeling of geophysical system instabilities is numerically challenging, given the potential effects of purely numerical artifacts. Here we explore the sensitivity of ice stream surge activation to numerical and physical model aspects. We find that surge characteristics exhibit a resolution dependency but converge at higher horizontal grid resolutions and are significantly affected by the incorporation of bed thermal and sub-glacial hydrology models.
Yannic Fischler, Thomas Kleiner, Christian Bischof, Jeremie Schmiedel, Roiy Sayag, Raban Emunds, Lennart Frederik Oestreich, and Angelika Humbert
Geosci. Model Dev., 16, 5305–5322, https://doi.org/10.5194/gmd-16-5305-2023, https://doi.org/10.5194/gmd-16-5305-2023, 2023
Short summary
Short summary
Water underneath ice sheets affects the motion of glaciers. This study presents a newly developed code, CUAS-MPI, that simulates subglacial hydrology. It is designed for supercomputers and is hence a parallelized code. We measure the performance of this code for simulations of the entire Greenland Ice Sheet and find that the code works efficiently. Moreover, we validated the code to ensure the correctness of the solution. CUAS-MPI opens new possibilities for simulations of ice sheet hydrology.
Julia Kaltenborn, Amy R. Macfarlane, Viviane Clay, and Martin Schneebeli
Geosci. Model Dev., 16, 4521–4550, https://doi.org/10.5194/gmd-16-4521-2023, https://doi.org/10.5194/gmd-16-4521-2023, 2023
Short summary
Short summary
Snow layer segmentation and snow grain classification are essential diagnostic tasks for cryospheric applications. A SnowMicroPen (SMP) can be used to that end; however, the manual classification of its profiles becomes infeasible for large datasets. Here, we evaluate how well machine learning models automate this task. Of the 14 models trained on the MOSAiC SMP dataset, the long short-term memory model performed the best. The findings presented here facilitate and accelerate SMP data analysis.
Johannes Aschauer, Adrien Michel, Tobias Jonas, and Christoph Marty
Geosci. Model Dev., 16, 4063–4081, https://doi.org/10.5194/gmd-16-4063-2023, https://doi.org/10.5194/gmd-16-4063-2023, 2023
Short summary
Short summary
Snow water equivalent is the mass of water stored in a snowpack. Based on exponential settling functions, the empirical snow density model SWE2HS is presented to convert time series of daily snow water equivalent into snow depth. The model has been calibrated with data from Switzerland and validated with independent data from the European Alps. A reference implementation of SWE2HS is available as a Python package.
Eric Keenan, Nander Wever, Jan T. M. Lenaerts, and Brooke Medley
Geosci. Model Dev., 16, 3203–3219, https://doi.org/10.5194/gmd-16-3203-2023, https://doi.org/10.5194/gmd-16-3203-2023, 2023
Short summary
Short summary
Ice sheets gain mass via snowfall. However, snowfall is redistributed by the wind, resulting in accumulation differences of up to a factor of 5 over distances as short as 5 km. These differences complicate estimates of ice sheet contribution to sea level rise. For this reason, we have developed a new model for estimating wind-driven snow redistribution on ice sheets. We show that, over Pine Island Glacier in West Antarctica, the model improves estimates of snow accumulation variability.
Sebastian Westermann, Thomas Ingeman-Nielsen, Johanna Scheer, Kristoffer Aalstad, Juditha Aga, Nitin Chaudhary, Bernd Etzelmüller, Simon Filhol, Andreas Kääb, Cas Renette, Louise Steffensen Schmidt, Thomas Vikhamar Schuler, Robin B. Zweigel, Léo Martin, Sarah Morard, Matan Ben-Asher, Michael Angelopoulos, Julia Boike, Brian Groenke, Frederieke Miesner, Jan Nitzbon, Paul Overduin, Simone M. Stuenzi, and Moritz Langer
Geosci. Model Dev., 16, 2607–2647, https://doi.org/10.5194/gmd-16-2607-2023, https://doi.org/10.5194/gmd-16-2607-2023, 2023
Short summary
Short summary
The CryoGrid community model is a new tool for simulating ground temperatures and the water and ice balance in cold regions. It is a modular design, which makes it possible to test different schemes to simulate, for example, permafrost ground in an efficient way. The model contains tools to simulate frozen and unfrozen ground, snow, glaciers, and other massive ice bodies, as well as water bodies.
Alex S. Gardner, Nicole-Jeanne Schlegel, and Eric Larour
Geosci. Model Dev., 16, 2277–2302, https://doi.org/10.5194/gmd-16-2277-2023, https://doi.org/10.5194/gmd-16-2277-2023, 2023
Short summary
Short summary
This is the first description of the open-source Glacier Energy and Mass Balance (GEMB) model. GEMB models the ice sheet and glacier surface–atmospheric energy and mass exchange, as well as the firn state. The model is evaluated against the current state of the art and in situ observations and is shown to perform well.
Yafei Nie, Chengkun Li, Martin Vancoppenolle, Bin Cheng, Fabio Boeira Dias, Xianqing Lv, and Petteri Uotila
Geosci. Model Dev., 16, 1395–1425, https://doi.org/10.5194/gmd-16-1395-2023, https://doi.org/10.5194/gmd-16-1395-2023, 2023
Short summary
Short summary
State-of-the-art Earth system models simulate the observed sea ice extent relatively well, but this is often due to errors in the dynamic and other processes in the simulated sea ice changes cancelling each other out. We assessed the sensitivity of these processes simulated by the coupled ocean–sea ice model NEMO4.0-SI3 to 18 parameters. The performance of the model in simulating sea ice change processes was ultimately improved by adjusting the three identified key parameters.
Varun Sharma, Franziska Gerber, and Michael Lehning
Geosci. Model Dev., 16, 719–749, https://doi.org/10.5194/gmd-16-719-2023, https://doi.org/10.5194/gmd-16-719-2023, 2023
Short summary
Short summary
Most current generation climate and weather models have a relatively simplistic description of snow and snow–atmosphere interaction. One reason for this is the belief that including an advanced snow model would make the simulations too computationally demanding. In this study, we bring together two state-of-the-art models for atmosphere (WRF) and snow cover (SNOWPACK) and highlight both the feasibility and necessity of such coupled models to explore underexplored phenomena in the cryosphere.
Anne M. Felden, Daniel F. Martin, and Esmond G. Ng
Geosci. Model Dev., 16, 407–425, https://doi.org/10.5194/gmd-16-407-2023, https://doi.org/10.5194/gmd-16-407-2023, 2023
Short summary
Short summary
We present and validate a novel subglacial hydrology model, SUHMO, based on an adaptive mesh refinement framework. We propose the addition of a pseudo-diffusion to recover the wall melting in channels. Computational performance analysis demonstrates the efficiency of adaptive mesh refinement on large-scale hydrologic problems. The adaptive mesh refinement approach will eventually enable better ice bed boundary conditions for ice sheet simulations at a reasonable computational cost.
Dalei Hao, Gautam Bisht, Karl Rittger, Edward Bair, Cenlin He, Huilin Huang, Cheng Dang, Timbo Stillinger, Yu Gu, Hailong Wang, Yun Qian, and L. Ruby Leung
Geosci. Model Dev., 16, 75–94, https://doi.org/10.5194/gmd-16-75-2023, https://doi.org/10.5194/gmd-16-75-2023, 2023
Short summary
Short summary
Snow with the highest albedo of land surface plays a vital role in Earth’s surface energy budget and water cycle. This study accounts for the impacts of snow grain shape and mixing state of light-absorbing particles with snow on snow albedo in the E3SM land model. The findings advance our understanding of the role of snow grain shape and mixing state of LAP–snow in land surface processes and offer guidance for improving snow simulations and radiative forcing estimates in Earth system models.
Esteban Alonso-González, Kristoffer Aalstad, Mohamed Wassim Baba, Jesús Revuelto, Juan Ignacio López-Moreno, Joel Fiddes, Richard Essery, and Simon Gascoin
Geosci. Model Dev., 15, 9127–9155, https://doi.org/10.5194/gmd-15-9127-2022, https://doi.org/10.5194/gmd-15-9127-2022, 2022
Short summary
Short summary
Snow cover plays an important role in many processes, but its monitoring is a challenging task. The alternative is usually to simulate the snowpack, and to improve these simulations one of the most promising options is to fuse simulations with available observations (data assimilation). In this paper we present MuSA, a data assimilation tool which facilitates the implementation of snow monitoring initiatives, allowing the assimilation of a wide variety of remotely sensed snow cover information.
Vincent Verjans, Alexander A. Robel, Helene Seroussi, Lizz Ultee, and Andrew F. Thompson
Geosci. Model Dev., 15, 8269–8293, https://doi.org/10.5194/gmd-15-8269-2022, https://doi.org/10.5194/gmd-15-8269-2022, 2022
Short summary
Short summary
We describe the development of the first large-scale ice sheet model that accounts for stochasticity in a range of processes. Stochasticity allows the impacts of inherently uncertain processes on ice sheets to be represented. This includes climatic uncertainty, as the climate is inherently chaotic. Furthermore, stochastic capabilities also encompass poorly constrained glaciological processes that display strong variability at fine spatiotemporal scales. We present the model and test experiments.
Max Brils, Peter Kuipers Munneke, Willem Jan van de Berg, and Michiel van den Broeke
Geosci. Model Dev., 15, 7121–7138, https://doi.org/10.5194/gmd-15-7121-2022, https://doi.org/10.5194/gmd-15-7121-2022, 2022
Short summary
Short summary
Firn covers the Greenland ice sheet (GrIS) and can temporarily prevent mass loss. Here, we present the latest version of our firn model, IMAU-FDM, with an application to the GrIS. We improved the density of fallen snow, the firn densification rate and the firn's thermal conductivity. This leads to a higher air content and 10 m temperatures. Furthermore we investigate three case studies and find that the updated model shows greater variability and an increased sensitivity in surface elevation.
Océane Hames, Mahdi Jafari, David Nicholas Wagner, Ian Raphael, David Clemens-Sewall, Chris Polashenski, Matthew D. Shupe, Martin Schneebeli, and Michael Lehning
Geosci. Model Dev., 15, 6429–6449, https://doi.org/10.5194/gmd-15-6429-2022, https://doi.org/10.5194/gmd-15-6429-2022, 2022
Short summary
Short summary
This paper presents an Eulerian–Lagrangian snow transport model implemented in the fluid dynamics software OpenFOAM, which we call snowBedFoam 1.0. We apply this model to reproduce snow deposition on a piece of ridged Arctic sea ice, which was produced during the MOSAiC expedition through scan measurements. The model appears to successfully reproduce the enhanced snow accumulation and deposition patterns, although some quantitative uncertainties were shown.
Constantijn J. Berends, Heiko Goelzer, Thomas J. Reerink, Lennert B. Stap, and Roderik S. W. van de Wal
Geosci. Model Dev., 15, 5667–5688, https://doi.org/10.5194/gmd-15-5667-2022, https://doi.org/10.5194/gmd-15-5667-2022, 2022
Short summary
Short summary
The rate at which marine ice sheets such as the West Antarctic ice sheet will retreat in a warming climate and ocean is still uncertain. Numerical ice-sheet models, which solve the physical equations that describe the way glaciers and ice sheets deform and flow, have been substantially improved in recent years. Here we present the results of several years of work on IMAU-ICE, an ice-sheet model of intermediate complexity, which can be used to study ice sheets of both the past and the future.
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.
Cited articles
Arthern, R. J., Winebrenner, D. P., and Vaughan, D. G.: Antarctic snow
accumulation mapped using polarization of 4.3-cm wavelength microwave
emission, J. Geophys. Res.-Atmos., 111, D06107,
https://doi.org/10.1029/2004JD005667, 2006. a
Aschwanden, A., Bueler, E., Khroulev, C., and Blatter, H.: An enthalpy
formulation for glaciers and ice sheets, J. Glaciol., 58, 441–457,
https://doi.org/10.3189/2012JoG11J088, 2012. a
AWI FusionForge: Alfred Wegener Institute for Polar and Marine Research FusionForge, available at: https://swrepo1.awi.de, last access: 6 April 2020. a
Balmaseda, M. A.: Editorial for Ocean Reanalysis Intercomparison Special Issue,
Clim. Dynam., 49, 707–708, https://doi.org/10.1007/s00382-017-3813-8, 2017. a
Bamber, J. L., Layberry, R. L., and Gogenini, S. P.: A new ice thickness and
bed data set for the Greenland ice sheet 1. Measurement, data reduction,
and errors, J. Geophys. Res.-Atmos., 106,
33773–33780, https://doi.org/10.1029/2001JD900054, 2001. a
Bamber, J. L., Griggs, J. A., Hurkmans, R. T. W. L., Dowdeswell, J. A., Gogineni, S. P., Howat, I., Mouginot, J., Paden, J., Palmer, S., Rignot, E., and Steinhage, D.: A new bed elevation dataset for Greenland, The Cryosphere, 7, 499–510, https://doi.org/10.5194/tc-7-499-2013, 2013. a
Bernales, J., Rogozhina, I., Greve, R., and Thomas, M.: Comparison of hybrid schemes for the combination of shallow approximations in numerical simulations of the Antarctic Ice Sheet, The Cryosphere, 11, 247–265, https://doi.org/10.5194/tc-11-247-2017, 2017. a, b, c
Brinkerhoff, D. J. and Johnson, J. V.: Data assimilation and prognostic whole ice sheet modelling with the variationally derived, higher order, open source, and fully parallel ice sheet model VarGlaS, The Cryosphere, 7, 1161–1184, https://doi.org/10.5194/tc-7-1161-2013, 2013. a, b
Budd, W. F., Jenssen, D., and Smith, I. N.: A three-dimensional time-dependent
model of the West Antarctic ice sheet, Ann. Glaciol., 5, 29–36,
1984. a
Errico, R. M. and Vukicevic, T.: Sensitivity analysis using an adjoint of the
PSU-NCAR mesoseale model, Mon. Weather Rev., 120, 1644–1660,
https://doi.org/10.1175/1520-0493(1992)120<1644:SAUAAO>2.0.CO;2, 1992. a, b
Forget, G., Campin, J.-M., Heimbach, P., Hill, C. N., Ponte, R. M., and Wunsch, C.: ECCO version 4: an integrated framework for non-linear inverse modeling and global ocean state estimation, Geosci. Model Dev., 8, 3071–3104, https://doi.org/10.5194/gmd-8-3071-2015, 2015. a
Forth, S., Hovland, P., Phipps, E., Utke, J., and Walther, A. (Eds.): Recent
Advances in Algorithmic Differentiation, Lecture Notes in Computational
Science and Engineering, Springer Science & Business Media, Vol. 87,
https://doi.org/10.1007/978-3-642-30023-3, 2012. a
Fortuin, J. P. F. and Oerlemans, J.: Parameterization of the annual surface
temperature and mass balance of Antarctica, Ann. Glaciol., 14,
78–84, https://doi.org/10.3189/S0260305500008302, 1990. a
Fretwell, P., Pritchard, H. D., Vaughan, D. G., Bamber, J. L., Barrand, N. E., Bell, R., Bianchi, C., Bingham, R. G., Blankenship, D. D., Casassa, G., Catania, G., Callens, D., Conway, H., Cook, A. J., Corr, H. F. J., Damaske, D., Damm, V., Ferraccioli, F., Forsberg, R., Fujita, S., Gim, Y., Gogineni, P., Griggs, J. A., Hindmarsh, R. C. A., Holmlund, P., Holt, J. W., Jacobel, R. W., Jenkins, A., Jokat, W., Jordan, T., King, E. C., Kohler, J., Krabill, W., Riger-Kusk, M., Langley, K. A., Leitchenkov, G., Leuschen, C., Luyendyk, B. P., Matsuoka, K., Mouginot, J., Nitsche, F. O., Nogi, Y., Nost, O. A., Popov, S. V., Rignot, E., Rippin, D. M., Rivera, A., Roberts, J., Ross, N., Siegert, M. J., Smith, A. M., Steinhage, D., Studinger, M., Sun, B., Tinto, B. K., Welch, B. C., Wilson, D., Young, D. A., Xiangbin, C., and Zirizzotti, A.: Bedmap2: improved ice bed, surface and thickness datasets for Antarctica, The Cryosphere, 7, 375–393, https://doi.org/10.5194/tc-7-375-2013, 2013. a
Gelb, A. (Ed.): Applied Optimal Estimation, The MIT Press, 1974. a
Giering, R., Kaminski, T., and Slawig, T.: Generating efficient derivative code
with TAF: Adjoint and tangent linear Euler flow around an airfoil, Future
Gener. Comp. Sy., 21, 1345–1355,
https://doi.org/10.1016/j.future.2004.11.003, 2005. a, b
Giles, M. B.: Collected Matrix Derivative Results for Forward and Reverse Mode
Algorithmic Differentiation, in: Advances in Automatic Differentiation,
edited by: Bischof, C. H., Bücker, H. M., Hovland, P., Naumann, U., and
Utke, J., Springer Berlin Heidelberg, Berlin, Heidelberg, 35–44, 2008. a
Gillet-Chaulet, F., Gagliardini, O., Seddik, H., Nodet, M., Durand, G., Ritz, C., Zwinger, T., Greve, R., and Vaughan, D. G.: Greenland ice sheet contribution to sea-level rise from a new-generation ice-sheet model, The Cryosphere, 6, 1561–1576, https://doi.org/10.5194/tc-6-1561-2012, 2012. a, b
Glasser, N. F. and Scambos, T. A.: A structural glaciological analysis of the
2002 Larsen B ice-shelf collapse, J. Glaciol., 54, 3–16,
https://doi.org/10.3189/002214308784409017, 2008. a
Glen, J. W.: The creep of polycrystalline ice, P. R. Soc.
A, 228, 519–538, https://doi.org/10.1098/rspa.1955.0066, 1955. a
Goelzer, H., Nowicki, S., Edwards, T., Beckley, M., Abe-Ouchi, A., Aschwanden, A., Calov, R., Gagliardini, O., Gillet-Chaulet, F., Golledge, N. R., Gregory, J., Greve, R., Humbert, A., Huybrechts, P., Kennedy, J. H., Larour, E., Lipscomb, W. H., Le clec'h, S., Lee, V., Morlighem, M., Pattyn, F., Payne, A. J., Rodehacke, C., Rückamp, M., Saito, F., Schlegel, N., Seroussi, H., Shepherd, A., Sun, S., van de Wal, R., and Ziemen, F. A.: Design and results of the ice sheet model initialisation experiments initMIP-Greenland: an ISMIP6 intercomparison, The Cryosphere, 12, 1433–1460, https://doi.org/10.5194/tc-12-1433-2018, 2018. a, b
Goldberg, D. N. and Heimbach, P.: Parameter and state estimation with a time-dependent adjoint marine ice sheet model, The Cryosphere, 7, 1659–1678, https://doi.org/10.5194/tc-7-1659-2013, 2013. a, b, c
Goldberg, D. N. and Sergienko, O. V.: Data assimilation using a hybrid ice flow model, The Cryosphere, 5, 315–327, https://doi.org/10.5194/tc-5-315-2011, 2011. a, b
Goldberg, D. N., Heimbach, P., Joughin, I., and Smith, B.: Committed retreat of Smith, Pope, and Kohler Glaciers over the next 30 years inferred by transient model calibration, The Cryosphere, 9, 2429–2446, https://doi.org/10.5194/tc-9-2429-2015, 2015. a
Goldberg, D. N., Narayanan, S. H. K., Hascoet, L., and Utke, J.: An optimized treatment for algorithmic differentiation of an important glaciological fixed-point problem, Geosci. Model Dev., 9, 1891–1904, https://doi.org/10.5194/gmd-9-1891-2016, 2016. a, b
Greve, R.: A continuum-mechanical formulation for shallow polythermal ice
sheets, Philos. T. R. Soc. A, 355, 921–974,
https://doi.org/10.1098/rsta.1997.0050, 1997. a
Greve, R.: Geothermal heat flux distribution for the Greenland ice sheet,
derived by combining a global representation and information from deep ice
cores, Polar Data J., 3, 22–36, https://doi.org/10.20575/00000006, 2019. a
Greve, R. and Blatter, H.: Dynamics of Ice Sheets and Glaciers, Springer,
Berlin, Germany, https://doi.org/10.1007/978-3-642-03415-2, 2009. a, b
Greve, R. and Blatter, H.: Comparison of thermodynamics solvers in the
polythermal ice sheet model SICOPOLIS, Polar Sci., 10, 11–23,
https://doi.org/10.1016/j.polar.2015.12.004, 2016. a, b
Greve, R. and Calov, R.: Comparison of numerical schemes for the solution of
the ice-thickness equation in a dynamic/thermodynamic ice-sheet model,
J. Comput. Phys., 179, 649–664,
https://doi.org/10.1006/jcph.2002.7081, 2002. a
Griewank, A. and Walther, A.: Algorithm 799: revolve: an implementation of
checkpointing for the reverse or adjoint mode of computational
differentiation, ACM T. Math. Software, 26, 19–45,
https://doi.org/10.1145/347837.347846, 2000. a
Griewank, A. and Walther, A.: Evaluating Derivatives: Principles and Techniques
of Algorithmic Differentiation, Other Titles in Applied Mathematics, Society
for Industrial and Applied Mathematics, Philadelphia, PA, US, 2nd Edn.,
https://doi.org/10.1137/1.9780898717761, 2008. a
Hascoët, L. and Morlighem, M.: Source-to-source adjoint Algorithmic
Differentiation of an ice sheet model written in C, Optim. Method.
Softw., 33, 1–8, https://doi.org/10.1080/10556788.2017.1396600, 2018. a
Hascoët, L. and Utke, J.: Programming language features, usage patterns,
and the efficiency of generated adjoint code, Optim. Method.
Softw., 31, 885–903, https://doi.org/10.1080/10556788.2016.1146269, 2016. a
Heimbach, P., Hill, C., and Giering, R.: Automatic generation of efficient
adjoint code for a parallel Navier-Stokes solver, in: Computational Science
– ICCS 2002, edited by: Sloot, P. M. A., Hoekstra, A. G., Tan, C. J. K., and
Dongarra, J. J., Springer, Berlin, Heidelberg, 1019–1028, 2002. a
Heimbach, P., Hill, C., and Giering, R.: An efficient exact adjoint of the
parallel MIT general circulation model, generated via automatic
differentiation, Future Gener. Comp. Sy., 21, 1356–1371,
https://doi.org/10.1016/j.future.2004.11.010, 2005. a, b
Hoffman, M. J., Perego, M., Price, S. F., Lipscomb, W. H., Zhang, T., Jacobsen, D., Tezaur, I., Salinger, A. G., Tuminaro, R., and Bertagna, L.: MPAS-Albany Land Ice (MALI): a variable-resolution ice sheet model for Earth system modeling using Voronoi grids, Geosci. Model Dev., 11, 3747–3780, https://doi.org/10.5194/gmd-11-3747-2018, 2018. a
Hutter, K.: Theoretical Glaciology; Material Science of Ice and the Mechanics
of Glaciers and Ice Sheets, D. Reidel Publishing Company, Dordrecht, The
Netherlands, 1983. a
Isaac, T., Stadler, G., and Ghattas, O.: Solution of nonlinear Stokes
equations discretized by high-order finite elements on nonconforming and
anisotropic meshes, with application to ice sheet dynamics, SIAM J.
Sci. Comput., 37, B804–B833, https://doi.org/10.1137/140974407, 2015. a, b
Joughin, I., Tulaczyk, S., Bamber, J. L., Blankenship, D., Holt, J. W.,
Scambos, T., and Vaughan, D. G.: Basal conditions for Pine Island and
Thwaites Glaciers, West Antarctica, determined using satellite and
airborne data, J. Glaciol., 55, 245–257,
https://doi.org/10.3189/002214309788608705, 2009. a
Khazendar, A., Rignot, E., and Larour, E.: Larsen B Ice Shelf rheology
preceding its disintegration inferred by a control method, Geophys.
Res. Lett., 34, L19503, https://doi.org/10.1029/2007GL030980, 2007. a
Larour, E., Rignot, E., Joughin, I., and Aubry, D.: Rheology of the Ronne Ice
Shelf, Antarctica, inferred from satellite radar interferometry data using
an inverse control method, Geophys. Res. Lett., 32, L05503,
https://doi.org/10.1029/2004GL021693, 2005. a
Larour, E., Utke, J., Csatho, B., Schenk, A., Seroussi, H., Morlighem, M., Rignot, E., Schlegel, N., and Khazendar, A.: Inferred basal friction and surface mass balance of the Northeast Greenland Ice Stream using data assimilation of ICESat (Ice Cloud and land Elevation Satellite) surface altimetry and ISSM (Ice Sheet System Model), The Cryosphere, 8, 2335–2351, https://doi.org/10.5194/tc-8-2335-2014, 2014. a, b
Le Brocq, A. M., Payne, A. J., and Vieli, A.: An improved Antarctic dataset for high resolution numerical ice sheet models (ALBMAP v1), Earth Syst. Sci. Data, 2, 247–260, https://doi.org/10.5194/essd-2-247-2010, 2010. a
Liao, H.-J., Liu, J.-G., Wang, L., and Xiang, T.: Differentiable programming
tensor networks, Phys. Rev. X, 9, 031041,
https://doi.org/10.1103/PhysRevX.9.031041, 2019. a
Lliboutry, L. and Duval, P.: Various isotropic and anisotropic ices found in
glaciers and polar ice caps and their corresponding rheologies, Ann.
Geophys., 3, 207–224, 1985. a
Logan, L. C., Narayanan, S. H. K., Greve, R., and Heimbach, P.: SICOPOLIS-AD:
Quick-Start Manual, Technical Report ANL/MCS-TM-382 Rev 0.10, Argonne
National Laboratory, ANL, Argonne, IL, https://doi.org/10.2172/1499025, 2019. a, b, c
Logan, L. C., Narayanan, S. H. K., Greve, R., and Heimbach, P.: SICOPOLIS-AD v1 (Version 1), Zenodo, https://doi.org/10.5281/zenodo.3686393, 2020. a
MacAyeal, D. R.: Ice-shelf response to ice-stream discharge fluctuations: III.
The effects of ice-stream imbalance on the Ross Ice Shelf, Antarctica,
J. Glaciol., 35, 38–42, https://doi.org/10.3189/002214389793701545, 1989. a
MacAyeal, D. R., Firestone, J., and Waddington, E.: Paleothermometry by control
methods, J. Glaciol., 37, 326–338,
https://doi.org/10.3189/S0022143000005761, 1991. a
Marotzke, J., Giering, R., Zhang, K. Q., Stammer, D., Hill, C., and Lee, T.:
Construction of the adjoint MIT ocean general circulation model and
application to Atlantic heat transport variability, J. Geophys.
Res.-Oceans, 104, 29529–29547, https://doi.org/10.1029/1999JC900236, 1999. a
Marsiat, I.: Simulation of the northern hemisphere continental ice sheets over
the last glacial-interglacial cycle: Experiments with a latitude-longitude
vertically integrated ice sheet model coupled to zonally averaged climate
model, Paleoclimates, 1, 59–98, 1994. a
Meehl, G. A., Moss, R., Taylor, K. E., Eyring, V., Stouffer, R. J., Bony, S.,
and Stevens, B.: Climate model intercomparisons: preparing for the next
phase, Eos T. Am. Geophys. Un., 95, 77–78,
https://doi.org/10.1002/2014EO090001, 2014. a
Morland, L. W.: Thermomechanical balances of ice sheet flows, Geophys.
Astro. Fluid, 29, 237–266, https://doi.org/10.1080/03091928408248191,
1984. a
Morland, L. W.: Unconfined ice-shelf flow, in: Dynamics of the West Antarctic
Ice Sheet, edited by: van der Veen, C. J. and Oerlemans, J., D.
Reidel Publishing Company, Dordrecht, The Netherlands, 99–116, 1987. 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, L14502, 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, https://doi.org/10.1002/jgrf.20125, 2013. a, b
Mosbeux, C., Gillet-Chaulet, F., and Gagliardini, O.: Comparison of adjoint and nudging methods to initialise ice sheet model basal conditions, Geosci. Model Dev., 9, 2549–2562, https://doi.org/10.5194/gmd-9-2549-2016, 2016. a
Narayanan, S. H. K.: sriharikrishna/OpenAD: SICOPOLIS-AD v1 (Version SICOPOLIS-AD), Zenodo, https://doi.org/10.5281/zenodo.3361744, 2019. a
Naumann, U.: The Art of Differentiating Computer Programs: An Introduction to
Algorithmic Differentiation, Software, Environments and Tools, 333 pp.,
https://doi.org/10.1137/1.9781611972078, SIAM Press, 2012. a
Perego, M., Price, S., and Stadler, G.: Optimal initial conditions for coupling
ice sheet models to Earth system models, J. Geophys. Res.-Earth, 119, 1894–1917, https://doi.org/10.1002/2014JF003181, 2014. a
Petra, N., Zhu, H., Stadler, G., Hughes, T. J. R., and Ghattas, O.: An inexact
Gauss-Newton method for inversion of basal sliding and rheology parameters
in a nonlinear Stokes ice sheet model, J. Glaciol., 58, 889–903,
https://doi.org/10.3189/2012JoG11J182, 2012. a
Ritz, C.: EISMINT intercomparison experiment: comparison of existing
Greenland models, Laboratoire de Glaciologie et de Géophysique
de l’Environnement, Saint Martin d’Hères, France, 105
1997. a
Robinson, A., Calov, R., and Ganopolski, A.: An efficient regional energy-moisture balance model for simulation of the Greenland Ice Sheet response to climate change, The Cryosphere, 4, 129–144, https://doi.org/10.5194/tc-4-129-2010, 2010. a
Rückamp, M., Greve, R., and Humbert, A.: Comparative simulations of the
evolution of the Greenland ice sheet under simplified Paris Agreement
scenarios with the models SICOPOLIS and ISSM, Polar Sci., 21, 14–25,
https://doi.org/10.1016/j.polar.2018.12.003, 2019. a
Seroussi, H., Nowicki, S., Simon, E., Abe-Ouchi, A., Albrecht, T., Brondex, J., Cornford, S., Dumas, C., Gillet-Chaulet, F., Goelzer, H., Golledge, N. R., Gregory, J. M., Greve, R., Hoffman, M. J., Humbert, A., Huybrechts, P., Kleiner, T., Larour, E., Leguy, G., Lipscomb, W. H., Lowry, D., Mengel, M., Morlighem, M., Pattyn, F., Payne, A. J., Pollard, D., Price, S. F., Quiquet, A., Reerink, T. J., Reese, R., Rodehacke, C. B., Schlegel, N.-J., Shepherd, A., Sun, S., Sutter, J., Van Breedam, J., van de Wal, R. S. W., Winkelmann, R., and Zhang, T.: initMIP-Antarctica: an ice sheet model initialization experiment of ISMIP6, The Cryosphere, 13, 1441–1471, https://doi.org/10.5194/tc-13-1441-2019, 2019. a, b
SICOPOLIS.net: Ice sheet model SICOPOLIS, homepage, available at: http://www.sicopolis.net, last access: 6 April 2020. a
Streubel, T., Griewank, A., Radons, M., and Bernt, J.-U.: Representation and Analysis of Piecewise Linear Functions in Abs-Normal Form, in: System Modeling and Optimization, edited by: Pötzsche C., Heuberger C., Kaltenbacher B., Rendl F., CSMO 2013, IFIP Advances in Information and Communication Technology, 443, 327–336, https://doi.org/10.1007/978-3-662-45504-3_32, Springer, Berlin, Heidelberg, 1997. a
Talagrand, O. and Courtier, P.: Variational assimilation of meteorological
observations with the adjoint vorticity equation. I: Theory, Q.
J. Roy. Meteor. Soc., 113, 1311–1328,
https://doi.org/10.1002/qj.49711347812, 1987. a
Thacker, W. C. and Long, R. B.: Fitting dynamics to data, J.
Geophys. Res.-Oceans, 93, 1227–1240, https://doi.org/10.1029/JC093iC02p01227,
1988. a
Uppala, S. M., Kållberg, P. W., Simmons, A. J., Andrae, U.,
Da Costa Bechtold, V., Fiorino, M., Gibson, J. K., Haseler, J., Hernandez,
A., Kelly, G. A., Li, X., Onogi, K., Saarinen, S., Sokka, N., Allan, R. P.,
Andersson, E., Arpe, K., Balmaseda, M. A., Beljaars, A. C. M., Van De Berg,
L., Bidlot, J., Bormann, N., Caires, S., Chevallier, F., Dethof, A.,
Dragosavac, M., Fisher, M., Fuentes, M., Hagemann, S., Hólm, E., Hoskins,
B. J., Isaksen, L., Janssen, P. A. E. M., Jenne, R., Mcnally, A. P., Mahfouf,
J.-F., Morcrette, J.-J., Rayner, N. A., Saunders, R. W., Simon, P., Sterl,
A., Trenberth, K. E., Untch, A., Vasiljevic, D., Viterbo, P., and Woollen,
J.: The ERA-40 re-analysis, Q. J. Roy. Meteor.
Soc., 131, 2961–3012, https://doi.org/10.1256/qj.04.176, 2005.
a
Utke, J., Naumann, U., Fagan, M., Tallent, N., Strout, M., Heimbach, P., Hill,
C., Ozyurt, D., and Wunsch, C.: OpenAD/F: A modular open source tool for
automatic differentiation of Fortran codes, ACM T.
Math. Softw., 34, 18, https://doi.org/10.1145/1377596.1377598, 2008. a
Vieli, A. and Payne, A. J.: Application of control methods for modelling the
flow of Pine Island Glacier, West Antarctica, Ann. Glaciol., 36,
197–204, https://doi.org/10.3189/172756403781816338, 2003. a
Waddington, E. D., Neumann, T. A., Koutnik, M. R., Marshall, H.-P., and Morse,
D. L.: Inference of accumulation-rate patterns from deep layers in glaciers
and ice sheets, J. Glaciol., 53, 694–712,
https://doi.org/10.3189/002214307784409351, 2007. a
Weertman, J.: The theory of glacier sliding, J. Glaciol., 5,
287–303, https://doi.org/10.3189/S0022143000029038, 1964. a
Weis, M., Greve, R., and Hutter, K.: Theory of shallow ice shelves, Continuum
Mech. Therm., 11, 15–50, https://doi.org/10.1007/s001610050102, 1999. a
Wunsch, C. and Heimbach, P.: Practical global oceanic state estimation, Physica
D, 230, 197–208, https://doi.org/10.1016/j.physd.2006.09.040, 2007. a, b, c
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.
A new capability has been developed for the ice sheet model SICOPOLIS (SImulation COde for...
