Articles | Volume 16, issue 19
https://doi.org/10.5194/gmd-16-5627-2023
© Author(s) 2023. 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-16-5627-2023
© Author(s) 2023. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Development and technical paper
| 
10 Oct 2023
Development and technical paper |
| 10 Oct 2023
| 10 Oct 2023
Modeling sensitivities of thermally and hydraulically driven ice stream surge cycling
Department of Physics and Physical Oceanography, Memorial University of Newfoundland, St. John’s, NL, A1B 3X7, Canada
Lev Tarasov
Department of Physics and Physical Oceanography, Memorial University of Newfoundland, St. John’s, NL, A1B 3X7, Canada
Elisa Mantelli
Department of Earth and Environmental Sciences, Ludwig-Maximillians-Universitaet Munich, Theresienstr. 41, 80333 Munich, Germany
Alfred Wegener Institute for Polar and Marine Research, Am Alten Hafen 26, 27568 Bremerhaven, Germany
Institute for Marine and Antarctic Studies, University of Tasmania, 20 Castray Esplanade, Battery Point TAS 7004, Australia
Related authors
Kevin Hank, Robert J. Arthern, C. Rosie Williams, Alex M. Brisbourne, Andrew M. Smith, James A. Smith, Anna Wåhlin, and Sridhar Anandakrishnan
EGUsphere, https://doi.org/10.5194/egusphere-2025-764, https://doi.org/10.5194/egusphere-2025-764, 2025
This preprint is open for discussion and under review for The Cryosphere (TC).
Short summary
Short summary
The slipperiness beneath ice sheets is a key source of uncertainty in sea level rise projections. Using both observations and model output, we infer the most probable representation of basal slipperiness in ice sheet models, enabling more accurate projections. For Pine Island Glacier, our results provide support for a Coulomb-type sliding law and widespread low effective pressures, potentially increasing sliding velocities in prognostic simulations and, hence, sea level rise projections.
Lev Tarasov, Benoit S. Lecavalier, Kevin Hank, and David Pollard
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-175, https://doi.org/10.5194/gmd-2024-175, 2025
Preprint under review for GMD
Short summary
Short summary
We document the glacial system model (GSM), a 3D glaciological ice sheet systems model specifically designed for large ensemble modelling in glacial cycle contexts. The model is distinguished by the breadth of relevant processes represented for this context. This ranges from meltwater surface drainage with proglacial lake formation to state-of-the-art subglacial sediment production/transport/deposition. The other key distinguishing design feature is attention to addressing process uncertainties.
Kevin Hank and Lev Tarasov
Clim. Past, 20, 2499–2524, https://doi.org/10.5194/cp-20-2499-2024, https://doi.org/10.5194/cp-20-2499-2024, 2024
Short summary
Short summary
The ice-rafted debris signature of Heinrich events in marine sedimentary cores is usually attributed to massive ice discharge from the Laurentide Ice Sheet. However, the driving mechanism of this pulsed discharge remains unclear. We compare three previously proposed hypotheses and examine the role of relevant system processes. We find ice stream surge cycling is the most likely mechanism, but its character is sensitive to both the geothermal heat flux and the form of the basal drag law.
Kevin Hank, Robert J. Arthern, C. Rosie Williams, Alex M. Brisbourne, Andrew M. Smith, James A. Smith, Anna Wåhlin, and Sridhar Anandakrishnan
EGUsphere, https://doi.org/10.5194/egusphere-2025-764, https://doi.org/10.5194/egusphere-2025-764, 2025
This preprint is open for discussion and under review for The Cryosphere (TC).
Short summary
Short summary
The slipperiness beneath ice sheets is a key source of uncertainty in sea level rise projections. Using both observations and model output, we infer the most probable representation of basal slipperiness in ice sheet models, enabling more accurate projections. For Pine Island Glacier, our results provide support for a Coulomb-type sliding law and widespread low effective pressures, potentially increasing sliding velocities in prognostic simulations and, hence, sea level rise projections.
Benoit S. Lecavalier and Lev Tarasov
The Cryosphere, 19, 919–953, https://doi.org/10.5194/tc-19-919-2025, https://doi.org/10.5194/tc-19-919-2025, 2025
Short summary
Short summary
We present the evolution of the Antarctic Ice Sheet (AIS) over the last 200 kyr by means of a history-matching analysis where an updated observational database constrained ~ 10 000 model simulations. During peak glaciation at the Last Glacial Maximum (LGM), the best-fitting sub-ensemble of AIS simulations reached an excess grounded ice volume relative to the present of 9.2 to 26.5 m equivalent sea level relative to the present. The LGM AIS volume can help resolve the LGM missing-ice problem.
Marilena Sophie Geng, Lev Tarasov, and April Sue Dalton
EGUsphere, https://doi.org/10.5194/egusphere-2025-495, https://doi.org/10.5194/egusphere-2025-495, 2025
This preprint is open for discussion and under review for Climate of the Past (CP).
Short summary
Short summary
We used a fully coupled ice-climate model to simulate the last two glacial inceptions, and compare the ensemble simulated ice sheet evolution to limited geological data. Our results show that Northern Hemisphere ice sheets grew rapidly, sometimes merging in ways not previously assumed and that capturing one glacial inception does not guarantee capturing another. These findings improve our understanding of ice-age dynamics and highlight challenges in predicting past and future climate evolution.
Lev Tarasov, Benoit S. Lecavalier, Kevin Hank, and David Pollard
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-175, https://doi.org/10.5194/gmd-2024-175, 2025
Preprint under review for GMD
Short summary
Short summary
We document the glacial system model (GSM), a 3D glaciological ice sheet systems model specifically designed for large ensemble modelling in glacial cycle contexts. The model is distinguished by the breadth of relevant processes represented for this context. This ranges from meltwater surface drainage with proglacial lake formation to state-of-the-art subglacial sediment production/transport/deposition. The other key distinguishing design feature is attention to addressing process uncertainties.
Benoit S. Lecavalier and Lev Tarasov
EGUsphere, https://doi.org/10.5194/egusphere-2024-3268, https://doi.org/10.5194/egusphere-2024-3268, 2024
Short summary
Short summary
To simulate the past evolution of the Antarctic ice sheet (AIS) during past warm and cold periods, a modelling analysis was performed that compared thousands of AIS simulations to a large collection of field observations. As the AIS changes, so does the surface load which leads to crustal deformation, gravitational and sea-level change. The present-day rate of bedrock deformation due to past AIS changes is used with satellite observations to infer AIS changes due to contemporary climate change.
Ryan Love, Glenn A. Milne, Parviz Ajourlou, Soran Parang, Lev Tarasov, and Konstantin Latychev
Geosci. Model Dev., 17, 8535–8551, https://doi.org/10.5194/gmd-17-8535-2024, https://doi.org/10.5194/gmd-17-8535-2024, 2024
Short summary
Short summary
A relatively recent advance in glacial isostatic adjustment modeling has been the development of models that include 3D Earth structure, as opposed to 1D structure. However, a major limitation is the computational expense. We have developed a method using artificial neural networks to emulate the influence of 3D Earth models to affordably constrain the viscosity parameter space. Our results indicate that the misfits are of a scale such that useful predictions of relative sea level can be made.
Kevin Hank and Lev Tarasov
Clim. Past, 20, 2499–2524, https://doi.org/10.5194/cp-20-2499-2024, https://doi.org/10.5194/cp-20-2499-2024, 2024
Short summary
Short summary
The ice-rafted debris signature of Heinrich events in marine sedimentary cores is usually attributed to massive ice discharge from the Laurentide Ice Sheet. However, the driving mechanism of this pulsed discharge remains unclear. We compare three previously proposed hypotheses and examine the role of relevant system processes. We find ice stream surge cycling is the most likely mechanism, but its character is sensitive to both the geothermal heat flux and the form of the basal drag law.
Daniel H. Richards, Elisa Mantelli, Samuel S. Pegler, and Sandra Piazolo
EGUsphere, https://doi.org/10.5194/egusphere-2024-3067, https://doi.org/10.5194/egusphere-2024-3067, 2024
Short summary
Short summary
Ice behaves differently depending on its crystal orientation, but how this affects its flow is unclear. We combine a range of previous models into a common equation to better understand crystal alignment. We tested a range of previous models on ice streams and divides, discovering that the best fit to observations comes from a) assuming neighbouring crystals have the same stress, and b) through describing the effect of crystal orientation on the flow in a way that allows directional variation.
Robert G. Bingham, Julien A. Bodart, Marie G. P. Cavitte, Ailsa Chung, Rebecca J. Sanderson, Johannes C. R. Sutter, Olaf Eisen, Nanna B. Karlsson, Joseph A. MacGregor, Neil Ross, Duncan A. Young, David W. Ashmore, Andreas Born, Winnie Chu, Xiangbin Cui, Reinhard Drews, Steven Franke, Vikram Goel, John W. Goodge, A. Clara J. Henry, Antoine Hermant, Benjamin H. Hills, Nicholas Holschuh, Michelle R. Koutnik, Gwendolyn J.-M. C. Leysinger Vieli, Emma J. Mackie, Elisa Mantelli, Carlos Martín, Felix S. L. Ng, Falk M. Oraschewski, Felipe Napoleoni, Frédéric Parrenin, Sergey V. Popov, Therese Rieckh, Rebecca Schlegel, Dustin M. Schroeder, Martin J. Siegert, Xueyuan Tang, Thomas O. Teisberg, Kate Winter, Shuai Yan, Harry Davis, Christine F. Dow, Tyler J. Fudge, Tom A. Jordan, Bernd Kulessa, Kenichi Matsuoka, Clara J. Nyqvist, Maryam Rahnemoonfar, Matthew R. Siegfried, Shivangini Singh, Verjan Višnjević, Rodrigo Zamora, and Alexandra Zuhr
EGUsphere, https://doi.org/10.5194/egusphere-2024-2593, https://doi.org/10.5194/egusphere-2024-2593, 2024
Short summary
Short summary
The ice sheets covering Antarctica have built up over millenia through successive snowfall events which become buried and preserved as internal surfaces of equal age detectable with ice-penetrating radar. This paper describes an international initiative to work together on this archival data to build a comprehensive 3-D picture of how old the ice is everywhere across Antarctica, and how this will be used to reconstruct past and predict future ice and climate behaviour.
Matthew Drew and Lev Tarasov
EGUsphere, https://doi.org/10.5194/egusphere-2024-620, https://doi.org/10.5194/egusphere-2024-620, 2024
Short summary
Short summary
We model the sediment-ice-climate system over North America for the last 2.58 Myr showing that ice sheets are capable of excavating features the size of the Hudson bay. This work provides a basis for reconstructing past landscapes important to climate modelling efforts, helping us to understand past earth system change.
Brian R. Crow, Lev Tarasov, Michael Schulz, and Matthias Prange
Clim. Past, 20, 281–296, https://doi.org/10.5194/cp-20-281-2024, https://doi.org/10.5194/cp-20-281-2024, 2024
Short summary
Short summary
An abnormally warm period around 400,000 years ago is thought to have resulted in a large melt event for the Greenland Ice Sheet. Using a sequence of climate model simulations connected to an ice model, we estimate a 50 % melt of Greenland compared to today. Importantly, we explore how the exact methodology of connecting the temperatures and precipitation from the climate model to the ice sheet model can influence these results and show that common methods could introduce errors.
Matthew Drew and Lev Tarasov
The Cryosphere, 17, 5391–5415, https://doi.org/10.5194/tc-17-5391-2023, https://doi.org/10.5194/tc-17-5391-2023, 2023
Short summary
Short summary
The interaction of fast-flowing regions of continental ice sheets with their beds governs how quickly they slide and therefore flow. The coupling of fast ice to its bed is controlled by the pressure of meltwater at its base. It is currently poorly understood how the physical details of these hydrologic systems affect ice speedup. Using numerical models we find, surprisingly, that they largely do not, except for the duration of the surge. This suggests that cheap models are sufficient.
Ryan Love, Lev Tarasov, Heather Andres, Alan Condron, Xu Zhang, and Gerrit Lohmann
EGUsphere, https://doi.org/10.5194/egusphere-2023-2225, https://doi.org/10.5194/egusphere-2023-2225, 2023
Preprint archived
Short summary
Short summary
Freshwater injection into bands across the North Atlantic are a mainstay of climate modelling when investigating topics such as climate change or the role of glacial runoff in the glacial climate system. However, this approach is unrealistic and results in a systematic bias in the climate response to a given flux of freshwater. We evaluate the magnitude of this bias by comparison to two other approaches for introducing freshwater into a coupled climate model setup for glacial conditions.
Benoit S. Lecavalier, Lev Tarasov, Greg Balco, Perry Spector, Claus-Dieter Hillenbrand, Christo Buizert, Catherine Ritz, Marion Leduc-Leballeur, Robert Mulvaney, Pippa L. Whitehouse, Michael J. Bentley, and Jonathan Bamber
Earth Syst. Sci. Data, 15, 3573–3596, https://doi.org/10.5194/essd-15-3573-2023, https://doi.org/10.5194/essd-15-3573-2023, 2023
Short summary
Short summary
The Antarctic Ice Sheet Evolution constraint database version 2 (AntICE2) consists of a large variety of observations that constrain the evolution of the Antarctic Ice Sheet over the last glacial cycle. This includes observations of past ice sheet extent, past ice thickness, past relative sea level, borehole temperature profiles, and present-day bedrock displacement rates. The database is intended to improve our understanding of past Antarctic changes and for ice sheet model calibrations.
Lev Tarasov and Michael Goldstein
EGUsphere, https://doi.org/10.5194/egusphere-2022-1410, https://doi.org/10.5194/egusphere-2022-1410, 2023
Preprint archived
Short summary
Short summary
This overview: 1. Illustrates how current climate and/or ice sheet model-based inferences about the past tend to have little interpretable value about the real world given inadequate accounting of uncertainties. 2. Explains Bayesian inference to a non-statistical community. 3. Sketches out some tractable inferential steps for computationally expensive models in a way that meaningfully accounts for uncertainties. 4. Lays out some steps for the community to move forward.
Lev Tarasov and Michael Goldstein
Clim. Past Discuss., https://doi.org/10.5194/cp-2021-145, https://doi.org/10.5194/cp-2021-145, 2021
Revised manuscript not accepted
Short summary
Short summary
This review: 1. Illustrates how current climate and/or ice sheet model-based inferences about the past tend to have limited interpretable value about the real world given inadequate accounting of uncertainties. 2. Explains Bayesian inference to a non-statistical community. 3. Sketches out tractable Bayesian inference for computationally expensive models in a way that meaningfully accounts for uncertainties. 4. Lays out some steps for the community to move forward.
Ryan Love, Heather J. Andres, Alan Condron, and Lev Tarasov
Clim. Past, 17, 2327–2341, https://doi.org/10.5194/cp-17-2327-2021, https://doi.org/10.5194/cp-17-2327-2021, 2021
Short summary
Short summary
Freshwater, in the form of glacial runoff, is hypothesized to play a critical role in centennial- to millennial-scale climate variability and climate transitions. We track the routing of glaciologically constrained freshwater volumes in glacial ocean simulations. Our simulations capture important generally not well-represented small-scale features (boundary currents, eddies). We show that the dilution of freshwater as it is transported to key climate regions reduces the freshening to 20 %–60 %.
Taimaz Bahadory, Lev Tarasov, and Heather Andres
Clim. Past, 17, 397–418, https://doi.org/10.5194/cp-17-397-2021, https://doi.org/10.5194/cp-17-397-2021, 2021
Short summary
Short summary
We present an ensemble of last glacial inception simulations using a fully coupled ice–climate model for the Northern Hemisphere. The ensemble largely captures inferred ice volume changes within proxy uncertainties. Notable features include an ice bridge across Davis Strait and between Greenland and Iceland. Via an equilibrium climate response experiment, we also demonstrate the potential value of fully coupled ice–climate modelling of last glacial inception to constrain future climate change.
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.
Heather J. Andres and Lev Tarasov
Clim. Past, 15, 1621–1646, https://doi.org/10.5194/cp-15-1621-2019, https://doi.org/10.5194/cp-15-1621-2019, 2019
Short summary
Short summary
Abrupt climate shifts of large magnitudes were common during glacial states, with explanations centred on the oceans. However, winds drive ocean surface currents so shifts in mean wind conditions could also have played a critical role. In a small ensemble of transient deglacial simulations, we find abrupt shifts in both jet stream location and variability over the North Atlantic. We show that the eastern North American ice sheet margin strongly constrains regional jet characteristics.
Laurie Menviel, Emilie Capron, Aline Govin, Andrea Dutton, Lev Tarasov, Ayako Abe-Ouchi, Russell N. Drysdale, Philip L. Gibbard, Lauren Gregoire, Feng He, Ruza F. Ivanovic, Masa Kageyama, Kenji Kawamura, Amaelle Landais, Bette L. Otto-Bliesner, Ikumi Oyabu, Polychronis C. Tzedakis, Eric Wolff, and Xu Zhang
Geosci. Model Dev., 12, 3649–3685, https://doi.org/10.5194/gmd-12-3649-2019, https://doi.org/10.5194/gmd-12-3649-2019, 2019
Short summary
Short summary
As part of the Past Global Changes (PAGES) working group on Quaternary Interglacials, we propose a protocol to perform transient simulations of the penultimate deglaciation for the Paleoclimate Modelling Intercomparison Project (PMIP4). This design includes time-varying changes in orbital forcing, greenhouse gas concentrations, continental ice sheets as well as freshwater input from the disintegration of continental ice sheets. Key paleo-records for model-data comparison are also included.
Taimaz Bahadory and Lev Tarasov
Geosci. Model Dev., 11, 3883–3902, https://doi.org/10.5194/gmd-11-3883-2018, https://doi.org/10.5194/gmd-11-3883-2018, 2018
Short summary
Short summary
We describe a new coupling between the Glacial Systems Model and the
LOVECLIM intermediate complexity climate model. The coupling is
distinguished from that of previous studies by greater completeness
and accuracy, with the intent of capturing the major feedbacks between
ice sheets and climate on glacial cycle timescales. The fully coupled
model will be used to examine the ice/climate phase space of past
glacial cycles.
Laurie Menviel, Emilie Capron, Aline Govin, Andrea Dutton, Lev Tarasov, Ayako Abe-Ouchi, Russell Drysdale, Philip Gibbard, Lauren Gregoire, Feng He, Ruza Ivanovic, Masa Kageyama, Kenji Kawamura, Amaelle Landais, Bette L. Otto-Bliesner, Ikumi Oyabu, Polychronis Tzedakis, Eric Wolff, and Xu Zhang
Clim. Past Discuss., https://doi.org/10.5194/cp-2018-106, https://doi.org/10.5194/cp-2018-106, 2018
Preprint withdrawn
Short summary
Short summary
The penultimate deglaciation (~ 138–128 ka), which represents the transition into the Last Interglacial period, provides a framework to investigate the climate and environmental response to large changes in boundary conditions. Here, as part of the PAGES-PMIP working group on Quaternary Interglacials, we propose a protocol to perform transient simulations of the penultimate deglaciation as well as a selection of paleo records for upcoming model-data comparisons.
Mark Kavanagh and Lev Tarasov
Geosci. Model Dev., 11, 3497–3513, https://doi.org/10.5194/gmd-11-3497-2018, https://doi.org/10.5194/gmd-11-3497-2018, 2018
Short summary
Short summary
We present and validate BrAHMs (BAsal Hydrology Model): a new
physically based basal hydrology model, which captures the two main
types of subglacial drainage systems (high-pressure distributed systems and
low-pressure channelized systems). BrAHMs is designed for continental
glacial cycle scale contexts, for which computational speed is
essential. This speed is accomplished, in part, by numerical methods
novel to basal hydrology contexts.
Masa Kageyama, Pascale Braconnot, Sandy P. Harrison, Alan M. Haywood, Johann H. Jungclaus, Bette L. Otto-Bliesner, Jean-Yves Peterschmitt, Ayako Abe-Ouchi, Samuel Albani, Patrick J. Bartlein, Chris Brierley, Michel Crucifix, Aisling Dolan, Laura Fernandez-Donado, Hubertus Fischer, Peter O. Hopcroft, Ruza F. Ivanovic, Fabrice Lambert, Daniel J. Lunt, Natalie M. Mahowald, W. Richard Peltier, Steven J. Phipps, Didier M. Roche, Gavin A. Schmidt, Lev Tarasov, Paul J. Valdes, Qiong Zhang, and Tianjun Zhou
Geosci. Model Dev., 11, 1033–1057, https://doi.org/10.5194/gmd-11-1033-2018, https://doi.org/10.5194/gmd-11-1033-2018, 2018
Short summary
Short summary
The Paleoclimate Modelling Intercomparison Project (PMIP) takes advantage of the existence of past climate states radically different from the recent past to test climate models used for climate projections and to better understand these climates. This paper describes the PMIP contribution to CMIP6 (Coupled Model Intercomparison Project, 6th phase) and possible analyses based on PMIP results, as well as on other CMIP6 projects.
Mohammad Hizbul Bahar Arif, Lev Tarasov, and Tristan Hauser
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2017-276, https://doi.org/10.5194/gmd-2017-276, 2018
Revised manuscript has not been submitted
Short summary
Short summary
This study is a first step answer to the following question: Can you
use emulators (machine learning techniques) to make the output of fast
simple climate models (a 2-D energy balance model in this test case)
indistinguishable from that of a much more computationally expensive
General Circulation climate model (GCM) within the uncertainties of
GCMs? Our preliminary test of this concept for large spatio-temporal
contexts gives a positive answer.
Masa Kageyama, Samuel Albani, Pascale Braconnot, Sandy P. Harrison, Peter O. Hopcroft, Ruza F. Ivanovic, Fabrice Lambert, Olivier Marti, W. Richard Peltier, Jean-Yves Peterschmitt, Didier M. Roche, Lev Tarasov, Xu Zhang, Esther C. Brady, Alan M. Haywood, Allegra N. LeGrande, Daniel J. Lunt, Natalie M. Mahowald, Uwe Mikolajewicz, Kerim H. Nisancioglu, Bette L. Otto-Bliesner, Hans Renssen, Robert A. Tomas, Qiong Zhang, Ayako Abe-Ouchi, Patrick J. Bartlein, Jian Cao, Qiang Li, Gerrit Lohmann, Rumi Ohgaito, Xiaoxu Shi, Evgeny Volodin, Kohei Yoshida, Xiao Zhang, and Weipeng Zheng
Geosci. Model Dev., 10, 4035–4055, https://doi.org/10.5194/gmd-10-4035-2017, https://doi.org/10.5194/gmd-10-4035-2017, 2017
Short summary
Short summary
The Last Glacial Maximum (LGM, 21000 years ago) is an interval when global ice volume was at a maximum, eustatic sea level close to a minimum, greenhouse gas concentrations were lower, atmospheric aerosol loadings were higher than today, and vegetation and land-surface characteristics were different from today. This paper describes the implementation of the LGM numerical experiment for the PMIP4-CMIP6 modelling intercomparison projects and the associated sensitivity experiments.
Ruza F. Ivanovic, Lauren J. Gregoire, Masa Kageyama, Didier M. Roche, Paul J. Valdes, Andrea Burke, Rosemarie Drummond, W. Richard Peltier, and Lev Tarasov
Geosci. Model Dev., 9, 2563–2587, https://doi.org/10.5194/gmd-9-2563-2016, https://doi.org/10.5194/gmd-9-2563-2016, 2016
Short summary
Short summary
This manuscript presents the experiment design for the PMIP4 Last Deglaciation Core experiment: a transient simulation of the last deglaciation, 21–9 ka. Specified model boundary conditions include time-varying orbital parameters, greenhouse gases, ice sheets, ice meltwater fluxes and other geographical changes (provided for 26–0 ka). The context of the experiment and the choices for the boundary conditions are explained, along with the future direction of the working group.
André Düsterhus, Alessio Rovere, Anders E. Carlson, Benjamin P. Horton, Volker Klemann, Lev Tarasov, Natasha L. M. Barlow, Tom Bradwell, Jorie Clark, Andrea Dutton, W. Roland Gehrels, Fiona D. Hibbert, Marc P. Hijma, Nicole Khan, Robert E. Kopp, Dorit Sivan, and Torbjörn E. Törnqvist
Clim. Past, 12, 911–921, https://doi.org/10.5194/cp-12-911-2016, https://doi.org/10.5194/cp-12-911-2016, 2016
Short summary
Short summary
This review/position paper addresses problems in creating new interdisciplinary databases for palaeo-climatological sea-level and ice-sheet data and gives an overview on new advances to tackle them. The focus therein is to define and explain strategies and highlight their importance to allow further progress in these fields. It also offers important insights into the general problem of designing competitive databases which are also applicable to other communities within the palaeo-environment.
A. Abe-Ouchi, F. Saito, M. Kageyama, P. Braconnot, S. P. Harrison, K. Lambeck, B. L. Otto-Bliesner, W. R. Peltier, L. Tarasov, J.-Y. Peterschmitt, and K. Takahashi
Geosci. Model Dev., 8, 3621–3637, https://doi.org/10.5194/gmd-8-3621-2015, https://doi.org/10.5194/gmd-8-3621-2015, 2015
Short summary
Short summary
We describe the creation of boundary conditions related to the presence of ice sheets, including ice-sheet extent and height, ice-shelf extent, and the distribution and altitude of ice-free land, at the Last Glacial Maximum (LGM), for use in LGM experiments conducted as part of the Coupled Modelling Intercomparison Project (CMIP5) and Palaeoclimate Modelling Intercomparison Project (PMIP3). The difference in the ice sheet boundary conditions as well as the climate response to them are discussed.
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.
H. Beltrami, G. S. Matharoo, L. Tarasov, V. Rath, and J. E. Smerdon
Clim. Past, 10, 1693–1706, https://doi.org/10.5194/cp-10-1693-2014, https://doi.org/10.5194/cp-10-1693-2014, 2014
R. Briggs, D. Pollard, and L. Tarasov
The Cryosphere, 7, 1949–1970, https://doi.org/10.5194/tc-7-1949-2013, https://doi.org/10.5194/tc-7-1949-2013, 2013
Related subject area
Cryosphere
Towards deep-learning solutions for classification of automated snow height measurements (CleanSnow v1.0.2)
Quantitative sub-ice and marine tracing of Antarctic sediment provenance (TASP v1.0)
Tuning parameters of a sea ice model using machine learning
WRF-Chem simulations of snow nitrate and other physicochemical properties in northern China
Clustering simulated snow profiles to form avalanche forecast regions
SnowQM 1.0: a fast R package for bias-correcting spatial fields of snow water equivalent using quantile mapping
Simulation of snow albedo and solar irradiance profile with the Two-streAm Radiative TransfEr in Snow (TARTES) v2.0 model
Computationally efficient subglacial drainage modelling using Gaussian Process emulators: GlaDS-GP v1.0
Evaluation of MITgcm-based ocean reanalyses for the Southern Ocean
Improvements in the land surface configuration to better simulate seasonal snow cover in the European Alps with the CNRM-AROME (cycle 46) convection-permitting regional climate model
A three-stage model pipeline predicting regional avalanche danger in Switzerland (RAvaFcast v1.0.0): a decision-support tool for operational avalanche forecasting
A 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)
The Utrecht Finite Volume Ice-Sheet Model (UFEMISM version 2.0) – part 1: description and idealised experiments
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)
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)
Jan Svoboda, Marc Ruesch, David Liechti, Corinne Jones, Michele Volpi, Michael Zehnder, and Jürg Schweizer
Geosci. Model Dev., 18, 1829–1849, https://doi.org/10.5194/gmd-18-1829-2025, https://doi.org/10.5194/gmd-18-1829-2025, 2025
Short summary
Short summary
Accurately measuring snow height is key for modeling approaches in climate science, snow hydrology, and avalanche forecasting. Erroneous snow height measurements often occur when snow height is low or changes, for instance during snowfall in summer. We prepare a new benchmark dataset with annotated snow height data and demonstrate how to improve the measurement quality using modern deep-learning approaches. Our approach can be easily implemented in a data pipeline for snow modeling.
James W. Marschalek, Edward Gasson, Tina van de Flierdt, Claus-Dieter Hillenbrand, Martin J. Siegert, and Liam Holder
Geosci. Model Dev., 18, 1673–1708, https://doi.org/10.5194/gmd-18-1673-2025, https://doi.org/10.5194/gmd-18-1673-2025, 2025
Short summary
Short summary
Ice sheet models can help predict how Antarctica's ice sheets respond to environmental change, and such models benefit from comparison to geological data. Here, we use an ice sheet model output and other data to predict the erosion of debris and trace its transport to where it is deposited on the ocean floor. This allows the results of ice sheet modelling to be directly and quantitively compared to real-world data, helping to reduce uncertainty regarding Antarctic sea level contribution.
Anton Korosov, Yue Ying, and Einar Ólason
Geosci. Model Dev., 18, 885–904, https://doi.org/10.5194/gmd-18-885-2025, https://doi.org/10.5194/gmd-18-885-2025, 2025
Short summary
Short summary
We have developed a new method to improve the accuracy of sea ice models, which predict how ice moves and deforms due to wind and ocean currents. Traditional models use parameters that are often poorly defined. The new approach uses machine learning to fine-tune these parameters by comparing simulated ice drift with satellite data. The method identifies optimal settings for the model by analysing patterns in ice deformation. This results in more accurate simulations of sea ice drift forecasting.
Xia Wang, Tao Che, Xueyin Ruan, Shanna Yue, Jing Wang, Chun Zhao, and Lei Geng
Geosci. Model Dev., 18, 651–670, https://doi.org/10.5194/gmd-18-651-2025, https://doi.org/10.5194/gmd-18-651-2025, 2025
Short summary
Short summary
We employed the WRF-Chem model to parameterize atmospheric nitrate deposition in snow and evaluate its performance in simulating snow cover, snow depth, and concentrations of dust and nitrate using new observations from northern China. The results generally exhibit reasonable agreement with field observations in northern China, demonstrating the model's capability to simulate snow properties, including concentrations of reservoir species.
Simon Horton, Florian Herla, and Pascal Haegeli
Geosci. Model Dev., 18, 193–209, https://doi.org/10.5194/gmd-18-193-2025, https://doi.org/10.5194/gmd-18-193-2025, 2025
Short summary
Short summary
We present a method for avalanche forecasters to analyze patterns in snowpack model simulations. It uses fuzzy clustering to group small regions into larger forecast areas based on snow characteristics, locations, and temporal history. Tested in the Columbia Mountains in two winter seasons, it closely matched real forecast regions regions and identified major avalanche hazard patterns. This approach simplifies complex model outputs, helping forecasters make informed decisions.
Adrien Michel, Johannes Aschauer, Tobias Jonas, Stefanie Gubler, Sven Kotlarski, and Christoph Marty
Geosci. Model Dev., 17, 8969–8988, https://doi.org/10.5194/gmd-17-8969-2024, https://doi.org/10.5194/gmd-17-8969-2024, 2024
Short summary
Short summary
We present a method to correct snow cover maps (represented in terms of snow water equivalent) to match better-quality maps. The correction can then be extended backwards and forwards in time for periods when better-quality maps are not available. The method is fast and gives good results. It is then applied to obtain a climatology of the snow cover in Switzerland over the past 60 years at a resolution of 1 d and 1 km. This is the first time that such a dataset has been produced.
Ghislain Picard and Quentin Libois
Geosci. Model Dev., 17, 8927–8953, https://doi.org/10.5194/gmd-17-8927-2024, https://doi.org/10.5194/gmd-17-8927-2024, 2024
Short summary
Short summary
The Two-streAm Radiative TransfEr in Snow (TARTES) is a radiative transfer model to compute snow albedo in the solar domain and the profiles of light and energy absorption in a multi-layered snowpack whose physical properties are user defined. It uniquely considers snow grain shape flexibly, based on recent insights showing that snow does not behave as a collection of ice spheres but instead as a random medium. TARTES is user-friendly yet performs comparably to more complex models.
Tim Hill, Derek Bingham, Gwenn E. Flowers, and Matthew J. Hoffman
EGUsphere, https://doi.org/10.22541/essoar.172736254.41350153/v2, https://doi.org/10.22541/essoar.172736254.41350153/v2, 2024
Short summary
Short summary
Subglacial drainage models represent water flow beneath glaciers and ice sheets. Here, we train fast statistical models called Gaussian Process emulators to accelerate subglacial drainage modelling by ~1000 times. We use the fast emulator predictions to show that three of the model parameters are responsible for >90 % of the variance in model outputs. The fast GP emulators will enable future uncertainty quantification and calibration of these models.
Yoshihiro Nakayama, Alena Malyarenko, Hong Zhang, Ou Wang, Matthis Auger, Yafei Nie, Ian Fenty, Matthew Mazloff, Armin Köhl, and Dimitris Menemenlis
Geosci. Model Dev., 17, 8613–8638, https://doi.org/10.5194/gmd-17-8613-2024, https://doi.org/10.5194/gmd-17-8613-2024, 2024
Short summary
Short summary
Global- and basin-scale ocean reanalyses are becoming easily accessible. However, such ocean reanalyses are optimized for their entire model domains and their ability to simulate the Southern Ocean requires evaluation. We conduct intercomparison analyses of Massachusetts Institute of Technology General Circulation Model (MITgcm)-based ocean reanalyses. They generally perform well for the open ocean, but open-ocean temporal variability and Antarctic continental shelves require improvements.
Diego Monteiro, Cécile Caillaud, Matthieu Lafaysse, Adrien Napoly, Mathieu Fructus, Antoinette Alias, and Samuel Morin
Geosci. Model Dev., 17, 7645–7677, https://doi.org/10.5194/gmd-17-7645-2024, https://doi.org/10.5194/gmd-17-7645-2024, 2024
Short summary
Short summary
Modeling snow cover in climate and weather forecasting models is a challenge even for high-resolution models. Recent simulations with CNRM-AROME have shown difficulties when representing snow in the European Alps. Using remote sensing data and in situ observations, we evaluate modifications of the land surface configuration in order to improve it. We propose a new surface configuration, enabling a more realistic simulation of snow cover, relevant for climate and weather forecasting applications.
Alessandro Maissen, Frank Techel, and Michele Volpi
Geosci. Model Dev., 17, 7569–7593, https://doi.org/10.5194/gmd-17-7569-2024, https://doi.org/10.5194/gmd-17-7569-2024, 2024
Short summary
Short summary
By harnessing AI models, this work enables processing large amounts of data, including weather conditions, snowpack characteristics, and historical avalanche data, to predict human-like avalanche forecasts in Switzerland. Our proposed model can significantly assist avalanche forecasters in their decision-making process, thereby facilitating more efficient and accurate predictions crucial for ensuring safety in Switzerland's avalanche-prone regions.
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.
Constantijn J. Berends, Victor Azizi, Jorge Bernales, and Roderik S. W. van de Wal
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-5, https://doi.org/10.5194/gmd-2024-5, 2024
Preprint under review for GMD
Short summary
Short summary
Ice-sheet models are computer programs that can simulate how the Greenland and Antarctic ice sheets will evolve in the future. The accuracy of these models depends on their resolution: how small the details are that the model can resolve. We have created a model with a variable resolution, which can resolve a lot of detail in areas where lots of changes happen in the ice, and less detail in areas where the ice does not move so much. This makes the model both accurate and fast.
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.
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.
Cited articles
Andrews, J. T. and MacLean, B.: Hudson Strait ice streams: A review of
stratigraphy, chronology and links with North Atlantic Heinrich events,
Boreas, 32, 4–17, https://doi.org/10.1080/03009480310001010, 2003. a
Arakawa, A. and Lamb, V. R.: Computational Design of the Basic Dynamical
Processes of the UCLA General Circulation Model, in: General Circulation
Models of the Atmosphere, edited by CHANG, J., vol. 17 of Methods in
Computational Physics: Advances in Research and Applications,
Elsevier, 173–265, https://doi.org/10.1016/B978-0-12-460817-7.50009-4,
1977. a, b
Bahadory, T. and Tarasov, L.: LCice 1.0 – a generalized Ice Sheet System Model coupler for LOVECLIM version 1.3: description, sensitivities, and validation with the Glacial Systems Model (GSM version D2017.aug17), Geosci. Model Dev., 11, 3883–3902, https://doi.org/10.5194/gmd-11-3883-2018, 2018. a, b
Benn, D. I., Fowler, A. C., Hewitt, I., and Sevestre, H.: A general theory of
glacier surges, J. Glaciol., 65, 701–716,
https://doi.org/10.1017/jog.2019.62, 2019. a
Bueler, E. and Brown, J.: Shallow shelf approximation as a “sliding law” in a
thermodynamically coupled ice sheet model, J. Geophys. Res., 114, F03008,
https://doi.org/10.1029/2008JF001179, 2009. a, b, c
Bueler, E. and van Pelt, W.: Mass-conserving subglacial hydrology in the Parallel Ice Sheet Model version 0.6, Geosci. Model Dev., 8, 1613–1635, https://doi.org/10.5194/gmd-8-1613-2015, 2015. a, b
Calov, R., Greve, R., Abe-Ouchi, A., Bueler, E., Huybrechts, P., Johnson,
J. V., Pattyn, F., Pollard, D., Ritz, C., Saito, F., and Tarasov, L.:
Results from the Ice-Sheet Model Intercomparison Project-Heinrich Event
INtercOmparison (ISMIP HEINO), J. Glaciol., 56, 371–383,
https://doi.org/10.3189/002214310792447789, 2010. a, b, c, d, e, f, g, h
Courant, R., Friedrichs, K., and Lewy, H.: Über die partiellen
Differenzengleichungen der mathematischen Physik, Math. Ann.,
100, 32–74, https://doi.org/10.1007/BF01448839, 1928. a
Cuffey, K. M. and Paterson, W. S. B.: The Physics of Glaciers,
Butterworth-Heinemann/Elsevier, Burlington, MA, 4th Edn., ISBN 9780123694614, 2010. a
Cuffey, K. M., Conway, H., Hallet, B., Gades, A. M., and Raymond, C. F.:
Interfacial water in polar glaciers and glacier sliding at -17 °C,
Geophys. Res. Lett., 26, 751–754, https://doi.org/10.1029/1999GL900096, 1999. a
Drew, M. and Tarasov, L.: Surging of a Hudson Strait Scale Ice Stream: Subglacial hydrology matters but the process details don't, The Cryosphere Discuss. [preprint], https://doi.org/10.5194/tc-2022-226, in review, 2022. a, b, c, d
Echelmeyer, K. and Zhongxiang, W.: Direct Observation of Basal Sliding and
Deformation of Basal Drift at Sub-Freezing Temperatures, J.
Glaciol., 33, 83–98, https://doi.org/10.3189/s0022143000005396, 1987. a
Flowers, G. E., Björnsson, H., and Pálsson, F.: New insights into
the subglacial and periglacial hydrology of Vatnajökull, Iceland, from
a distributed physical model, J. Glaciol., 49, 257–270,
https://doi.org/10.3189/172756503781830827, 2003. a
Fowler, A. C.: Sub-Temperate Basal Sliding, J. Glaciol., 32, 3–5,
https://doi.org/10.3189/S0022143000006808, 1986. a, b
Fowler, A. C. and Johnson, C.: Hydraulic run-away: a mechanism for thermally
regulated surges of ice sheets, J. Glaciol., 41, 554–561,
https://doi.org/10.3189/S002214300003478X, 1995. a
Fowler, A. C. and Schiavi, E.: A theory of ice-sheet surges, J.
Glaciol., 44, 104–118, https://doi.org/10.3189/s0022143000002409, 1998. a, b
Gandy, N., Gregoire, L. J., Ely, J. C., Cornford, S. L., Clark, C. D., and
Hodgson, D. M.: Exploring the ingredients required to successfully model the
placement, generation, and evolution of ice streams in the British-Irish Ice
Sheet, Quaternary Sci. Rev., 223, 105915,
https://doi.org/10.1016/j.quascirev.2019.105915, 2019. a
Hemming, S.: Heinrich events: Massive late Pleistocene detritus layers of the
North Atlantic and their global climate imprint, Rev. Geophys., 42, RG1005, https://doi.org/10.1029/2003RG000128, 2004. a
Joughin, I., Smith, B. E., Howat, I. M., Floricioiu, D., Alley, R. B., Truffer,
M., and Fahnestock, M.: Seasonal to decadal scale variations in the surface
velocity of Jakobshavn Isbrae, Greenland: Observation and model-based
analysis, J. Geophys. Res.-Earth, 117, 1–20,
https://doi.org/10.1029/2011JF002110, 2012. a
Joughin, I., Smith, B. E., Shean, D. E., and Floricioiu, D.: Brief Communication: Further summer speedup of Jakobshavn Isbræ, The Cryosphere, 8, 209–214, https://doi.org/10.5194/tc-8-209-2014, 2014. a
Khrulev, C., Bueler, E., Aschwanden, A., damaxwell, Albrecht, T., Brown, J.,
Seguinot, J., Mengel, M., Hinck, S., EnricoDeg, Ziemen, F., Blum, K.,
ronjareese, tkleiner, Bot, D., sschoell, and Kreuzer, M.: pism/pism: v2.0.2
bug fix release, Zenodo [code], https://doi.org/10.5281/zenodo.6001196, 2022. a
Kyrke-Smith, T. M., Katz, R. F., and Fowler, A. C.: Subglacial hydrology and
the formation of ice streams, P. Roy. Soc. A, 470, 20130494,
https://doi.org/10.1098/rspa.2013.0494, 2014. a
MacAyeal, D. R.: Binge/purge oscillations of the Laurentide Ice Sheet as a
cause of the North Atlantic's Heinrich events, Paleoceanography, 8, 775–784,
https://doi.org/10.1029/93PA02200, 1993. a
Mantelli, E., Bertagni, M. B., and Ridolfi, L.: Stochastic ice stream dynamics,
P. Natl. Acad. Sci. USA, 113, E4594–E4600,
https://doi.org/10.1073/pnas.1600362113, 2016. a, b
McCarthy, C., Savage, H., and Nettles, M.: Temperature dependence of
ice-on-rock friction at realistic glacier conditions, Philos.
T. Roy. Soc. A, 375, 20150348, https://doi.org/10.1098/rsta.2015.0348, 2017. a, b
Payne, A. J. and Dongelmans, P. W.: Self-organization in the thermomechanical
flow of ice sheets, J. Geophys. Res.-Sol. Ea., 102,
12219–12233, https://doi.org/10.1029/97jb00513, 1997. a
Payne, A. J., Huybrechts, P., Abe-Ouchi, A., Calov, R., Fastook, J. L., Greve,
R., Marshall, S. J., Marsiat, I., Ritz, C., Tarasov, L., and Thomassen,
M. P.: Results from the EISMINT model intercomparison: The effects of
thermomechanical coupling, J. Glaciol., 46, 227–238,
https://doi.org/10.3189/172756500781832891, 2000. a
PISM 2.0.6 documentation: PETSc options for PISM users,
https://www.pism.io/docs/manual/practical-usage/petsc-options.html (last access: 23 August 2023),
2023. a
Pollard, D. and DeConto, R. M.: A Coupled Ice-Sheet/Ice-Shelf/Sediment Model
Applied to a Marine-Margin Flowline: Forced and Unforced Variations, Glacial
Sedimentary Processes and Products, edited by: Montanez, I., Hambrey, M. J., Christoffersen, P., Glasser, N. F., and Hubbard, B., 37–52,
https://doi.org/10.1002/9781444304435.ch4, 2007. a
Pollard, D. and DeConto, R. M.: Description of a hybrid ice sheet-shelf model, and application to Antarctica, Geosci. Model Dev., 5, 1273–1295, https://doi.org/10.5194/gmd-5-1273-2012, 2012. a, b
Robel, A. A., Degiuli, E., Schoof, C., and Tziperman, E.: Dynamics of ice
stream temporal variability: Modes, scales, and hysteresis, J.
Geophys. Res.-Earth, 118, 925–936, https://doi.org/10.1002/jgrf.20072,
2013. a
Roberts, W. H. G., Payne, A. J., and Valdes, P. J.: The role of basal hydrology in the surging of the Laurentide Ice Sheet, Clim. Past, 12, 1601–1617, https://doi.org/10.5194/cp-12-1601-2016, 2016. a, b
Sayag, R. and Tziperman, E.: Interaction and variability of ice streams under
a triple-valued sliding law and non-Newtonian rheology, J.
Geophys. Res.-Earth, 116, F01009, https://doi.org/10.1029/2010JF001839, 2011. a
Schannwell, C., Mikolajewicz, U., Ziemen, F., and Kapsch, M.-L.: Sensitivity of Heinrich-type ice-sheet surge characteristics to boundary forcing perturbations, Clim. Past, 19, 179–198, https://doi.org/10.5194/cp-19-179-2023, 2023. a
Shreve, R. L.: Glacier sliding at subfreezing temperatures, J.
Glaciol., 30, 341–347, https://doi.org/10.1017/S0022143000006195, 1984. a
Takahama, R.: Heinrich Event Intercomparison with the ice-sheet model
SICOPOLIS, Master's thesis,
http://hdl.handle.net/2115/28749 (last access: 9 August 2022), 2006. a
Tarasov, L. and Peltier, W. R.: A high-resolution model of the 100 ka ice-age
cycle, Ann. Glaciol., 25, 58–65, https://doi.org/10.3189/s026030550001380x,
1997. a
Tarasov, L., Dyke, A. S., Neal, R. M., and Peltier, W. R.: A data-calibrated
distribution of deglacial chronologies for the North American ice complex
from glaciological modeling, Earth Planet. Sc. Lett., 315–316,
30–40, https://doi.org/10.1016/j.epsl.2011.09.010, 2012. a
Tarasov, L., Hank, K., and Lecavalier, B. S.: GSMv01.31.2023 code archive for
LISsq experiments, Zenodo [code], https://doi.org/10.5281/zenodo.7668472, 2023. a, b, c
Tulaczyk, S., Kamb, W. B., and Engelhardt, H. F.: Basal mechanics of Ice
Stream B, west Antarctica: 2. Undrained plastic bed model, J.
Geophys. Res.-Sol. Ea., 105, 483–494,
https://doi.org/10.1029/1999JB900328, 2000a. a
Tulaczyk, S., Kamb, W. B., and Engelhardt, H. F.: Basal mechanics of Ice
Stream B, West Antarctica 1. Till mechanics, J.
Geophys. Res.-Sol. Ea., 105, 463–481, https://doi.org/10.1029/1999jb900329,
2000b.
a
Van Pelt, W. J. and Oerlemans, J.: Numerical simulations of cyclic behaviour
in the Parallel Ice Sheet Model (PISM, J. Glaciol., 58, 347–360,
https://doi.org/10.3189/2012JoG11J217, 2012. a, b
Werder, M. A., Hewitt, I. J., Schoof, C. G., and Flowers, G. E.: Modeling
channelized and distributed subglacial drainage in two dimensions, J.
Geophys. Res.-Earth, 118, 2140–2158,
https://doi.org/10.1002/jgrf.20146, 2013. a
Winkelmann, R., Martin, M. A., Haseloff, M., Albrecht, T., Bueler, E., Khroulev, C., and Levermann, A.: The Potsdam Parallel Ice Sheet Model (PISM-PIK) – Part 1: Model description, The Cryosphere, 5, 715–726, https://doi.org/10.5194/tc-5-715-2011, 2011. a, b
Winsborrow, M. C., Clark, C. D., and Stokes, C. R.: What controls the location
of ice streams?, Earth-Sci. Rev., 103, 45–59,
https://doi.org/10.1016/j.earscirev.2010.07.003, 2010. a
Ziemen, F. A., Rodehacke, C. B., and Mikolajewicz, U.: Coupled ice sheet–climate modeling under glacial and pre-industrial boundary conditions, Clim. Past, 10, 1817–1836, https://doi.org/10.5194/cp-10-1817-2014, 2014. a
Ziemen, F. A., Kapsch, M.-L., Klockmann, M., and Mikolajewicz, U.: Heinrich events show two-stage climate response in transient glacial simulations, Clim. Past, 15, 153–168, https://doi.org/10.5194/cp-15-153-2019, 2019. a, b, c
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.
Physically meaningful modeling of geophysical system instabilities is numerically challenging,...
