Articles | Volume 18, issue 12
https://doi.org/10.5194/gmd-18-3895-2025
© Author(s) 2025. 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-18-3895-2025
© Author(s) 2025. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Model description paper
| 
01 Jul 2025
Model description paper |
| 01 Jul 2025
| 01 Jul 2025
Description and validation of the ice-sheet model Nix v1.0
Daniel Moreno-Parada
CORRESPONDING AUTHOR
Laboratoire de Glaciologie, Faculty of Sciences, Université Libre de Bruxelles, Av. F. D. Roosevelt 50, 1050 Brussels, Belgium
Alexander Robinson
Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research, Potsdam, Germany
Marisa Montoya
Departamento de Física de la Tierra y Astrofísica, Facultad de Ciencias Físicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
Instituto de Geociencias, Consejo Superior de Investigaciones Cientifícas, Universidad Complutense de Madrid, 28040 Madrid, Spain
Jorge Alvarez-Solas
CORRESPONDING AUTHOR
Departamento de Física de la Tierra y Astrofísica, Facultad de Ciencias Físicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
Instituto de Geociencias, Consejo Superior de Investigaciones Cientifícas, Universidad Complutense de Madrid, 28040 Madrid, Spain
Related authors
Lucía Gutiérrez-González, Alexander Robinson, Jorge Alvarez-Solas, Ilaria Tabone, Jan Swierczek-Jereczek, Daniel Moreno-Parada, and Marisa Montoya
EGUsphere, https://doi.org/10.5194/egusphere-2025-2616, https://doi.org/10.5194/egusphere-2025-2616, 2025
This preprint is open for discussion and under review for The Cryosphere (TC).
Short summary
Short summary
The Greenland ice sheet is considered a tipping element: if temperatures exceed its threshold, it would transition to a virtually ice-free state. We analyze its stability across the full range of glacial-interglacial temperatures, as well as those expected in the coming centuries. We find a future critical threshold between +1.5 and +2 K of global warming, another under colder climates, and persistent hysteresis across the full range of study.
Sergio Pérez-Montero, Jorge Alvarez-Solas, Jan Swierczek-Jereczek, Daniel Moreno-Parada, Alexander Robinson, and Marisa Montoya
Earth Syst. Dynam., 16, 915–937, https://doi.org/10.5194/esd-16-915-2025, https://doi.org/10.5194/esd-16-915-2025, 2025
Short summary
Short summary
The climate of the last 3 Myr has varied between cold and warm periods. Numerous independent mechanisms have been proposed to explain this; however, no effort has been made to study their competing effects. Here we present a simple but physically motivated model that includes these mechanisms in a modular way. We identify ice-sheet dynamics and lithosphere displacement as main triggers, but reproducing the climate records additionally requires the natural darkening of ice.
Sergio Pérez-Montero, Jorge Alvarez-Solas, Jan Swierczek-Jereczek, Daniel Moreno-Parada, Alexander Robinson, and Marisa Montoya
EGUsphere, https://doi.org/10.5194/egusphere-2025-2467, https://doi.org/10.5194/egusphere-2025-2467, 2025
This preprint is open for discussion and under review for Climate of the Past (CP).
Short summary
Short summary
Almost 3 million years ago, the planet began to experience a succession of cold and warm periods every 40,000 years. However, about 1 million years ago, they began to occur every 100,000 years. In this paper we explore how the change in the basal velocity of the ice sheets could have produced this change in behavior. On the other hand, we also see that in our model, decreasing in time the sensitivity of snowfall to temperature is also an effective mechanism with which to reproduce the records.
Daniel Moreno-Parada, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas
The Cryosphere, 18, 4215–4232, https://doi.org/10.5194/tc-18-4215-2024, https://doi.org/10.5194/tc-18-4215-2024, 2024
Short summary
Short summary
Our study tries to understand how the ice temperature evolves in a large mass as in the case of Antarctica. We found a relation that tells us the ice temperature at any point. These results are important because they also determine how the ice moves. In general, ice moves due to slow deformation (as if pouring honey from a jar). Nevertheless, in some regions the ice base warms enough and melts. The liquid water then serves as lubricant and the ice slides and its velocity increases rapidly.
Javier Blasco, Ilaria Tabone, Daniel Moreno-Parada, Alexander Robinson, Jorge Alvarez-Solas, Frank Pattyn, and Marisa Montoya
Clim. Past, 20, 1919–1938, https://doi.org/10.5194/cp-20-1919-2024, https://doi.org/10.5194/cp-20-1919-2024, 2024
Short summary
Short summary
In this study, we assess Antarctic tipping points which may had been crossed during the mid-Pliocene Warm Period. For this, we use data from the PlioMIP2 ensemble. Additionally, we investigate various sources of uncertainty, like ice dynamics and bedrock configuration. Our research significantly enhances our comprehension of Antarctica's response to a warming climate, shedding light on potential future tipping points that may be surpassed.
Daniel Moreno-Parada, Jorge Alvarez-Solas, Javier Blasco, Marisa Montoya, and Alexander Robinson
The Cryosphere, 17, 2139–2156, https://doi.org/10.5194/tc-17-2139-2023, https://doi.org/10.5194/tc-17-2139-2023, 2023
Short summary
Short summary
We have reconstructed the Laurentide Ice Sheet, located in North America during the Last Glacial Maximum (21 000 years ago). The absence of direct measurements raises a number of uncertainties. Here we study the impact of different physical laws that describe the friction as the ice slides over its base. We found that the Laurentide Ice Sheet is closest to prior reconstructions when the basal friction takes into account whether the base is frozen or thawed during its motion.
Lucía Gutiérrez-González, Alexander Robinson, Jorge Alvarez-Solas, Ilaria Tabone, Jan Swierczek-Jereczek, Daniel Moreno-Parada, and Marisa Montoya
EGUsphere, https://doi.org/10.5194/egusphere-2025-2616, https://doi.org/10.5194/egusphere-2025-2616, 2025
This preprint is open for discussion and under review for The Cryosphere (TC).
Short summary
Short summary
The Greenland ice sheet is considered a tipping element: if temperatures exceed its threshold, it would transition to a virtually ice-free state. We analyze its stability across the full range of glacial-interglacial temperatures, as well as those expected in the coming centuries. We find a future critical threshold between +1.5 and +2 K of global warming, another under colder climates, and persistent hysteresis across the full range of study.
Sergio Pérez-Montero, Jorge Alvarez-Solas, Jan Swierczek-Jereczek, Daniel Moreno-Parada, Alexander Robinson, and Marisa Montoya
Earth Syst. Dynam., 16, 915–937, https://doi.org/10.5194/esd-16-915-2025, https://doi.org/10.5194/esd-16-915-2025, 2025
Short summary
Short summary
The climate of the last 3 Myr has varied between cold and warm periods. Numerous independent mechanisms have been proposed to explain this; however, no effort has been made to study their competing effects. Here we present a simple but physically motivated model that includes these mechanisms in a modular way. We identify ice-sheet dynamics and lithosphere displacement as main triggers, but reproducing the climate records additionally requires the natural darkening of ice.
Ricarda Winkelmann, Donovan P. Dennis, Jonathan F. Donges, Sina Loriani, Ann Kristin Klose, Jesse F. Abrams, Jorge Alvarez-Solas, Torsten Albrecht, David Armstrong McKay, Sebastian Bathiany, Javier Blasco Navarro, Victor Brovkin, Eleanor Burke, Gokhan Danabasoglu, Reik V. Donner, Markus Drüke, Goran Georgievski, Heiko Goelzer, Anna B. Harper, Gabriele Hegerl, Marina Hirota, Aixue Hu, Laura C. Jackson, Colin Jones, Hyungjun Kim, Torben Koenigk, Peter Lawrence, Timothy M. Lenton, Hannah Liddy, José Licón-Saláiz, Maxence Menthon, Marisa Montoya, Jan Nitzbon, Sophie Nowicki, Bette Otto-Bliesner, Francesco Pausata, Stefan Rahmstorf, Karoline Ramin, Alexander Robinson, Johan Rockström, Anastasia Romanou, Boris Sakschewski, Christina Schädel, Steven Sherwood, Robin S. Smith, Norman J. Steinert, Didier Swingedouw, Matteo Willeit, Wilbert Weijer, Richard Wood, Klaus Wyser, and Shuting Yang
EGUsphere, https://doi.org/10.5194/egusphere-2025-1899, https://doi.org/10.5194/egusphere-2025-1899, 2025
This preprint is open for discussion and under review for Earth System Dynamics (ESD).
Short summary
Short summary
The Tipping Points Modelling Intercomparison Project (TIPMIP) is an international collaborative effort to systematically assess tipping point risks in the Earth system using state-of-the-art coupled and stand-alone domain models. TIPMIP will provide a first global atlas of potential tipping dynamics, respective critical thresholds and key uncertainties, generating an important building block towards a comprehensive scientific basis for policy- and decision-making.
Sergio Pérez-Montero, Jorge Alvarez-Solas, Jan Swierczek-Jereczek, Daniel Moreno-Parada, Alexander Robinson, and Marisa Montoya
EGUsphere, https://doi.org/10.5194/egusphere-2025-2467, https://doi.org/10.5194/egusphere-2025-2467, 2025
This preprint is open for discussion and under review for Climate of the Past (CP).
Short summary
Short summary
Almost 3 million years ago, the planet began to experience a succession of cold and warm periods every 40,000 years. However, about 1 million years ago, they began to occur every 100,000 years. In this paper we explore how the change in the basal velocity of the ice sheets could have produced this change in behavior. On the other hand, we also see that in our model, decreasing in time the sensitivity of snowfall to temperature is also an effective mechanism with which to reproduce the records.
Antonio Juarez-Martinez, Javier Blasco, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas
The Cryosphere, 18, 4257–4283, https://doi.org/10.5194/tc-18-4257-2024, https://doi.org/10.5194/tc-18-4257-2024, 2024
Short summary
Short summary
We present sea level projections for Antarctica in the context of ISMIP6-2300 with several forcings but extend the simulations to 2500, showing that more than 3 m of sea level contribution could be reached. We also test the sensitivity on a basal melting parameter and determine the timing of the loss of ice in the west region. All the simulations were carried out with the ice sheet model Yelmo.
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.
Daniel Moreno-Parada, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas
The Cryosphere, 18, 4215–4232, https://doi.org/10.5194/tc-18-4215-2024, https://doi.org/10.5194/tc-18-4215-2024, 2024
Short summary
Short summary
Our study tries to understand how the ice temperature evolves in a large mass as in the case of Antarctica. We found a relation that tells us the ice temperature at any point. These results are important because they also determine how the ice moves. In general, ice moves due to slow deformation (as if pouring honey from a jar). Nevertheless, in some regions the ice base warms enough and melts. The liquid water then serves as lubricant and the ice slides and its velocity increases rapidly.
Javier Blasco, Ilaria Tabone, Daniel Moreno-Parada, Alexander Robinson, Jorge Alvarez-Solas, Frank Pattyn, and Marisa Montoya
Clim. Past, 20, 1919–1938, https://doi.org/10.5194/cp-20-1919-2024, https://doi.org/10.5194/cp-20-1919-2024, 2024
Short summary
Short summary
In this study, we assess Antarctic tipping points which may had been crossed during the mid-Pliocene Warm Period. For this, we use data from the PlioMIP2 ensemble. Additionally, we investigate various sources of uncertainty, like ice dynamics and bedrock configuration. Our research significantly enhances our comprehension of Antarctica's response to a warming climate, shedding light on potential future tipping points that may be surpassed.
Jan Swierczek-Jereczek, Marisa Montoya, Konstantin Latychev, Alexander Robinson, Jorge Alvarez-Solas, and Jerry Mitrovica
Geosci. Model Dev., 17, 5263–5290, https://doi.org/10.5194/gmd-17-5263-2024, https://doi.org/10.5194/gmd-17-5263-2024, 2024
Short summary
Short summary
Ice sheets present a thickness of a few kilometres, leading to a vertical deformation of the crust of up to a kilometre. This process depends on properties of the solid Earth, which can be regionally very different. We propose a model that accounts for this often-ignored heterogeneity and run 100 000 simulation years in minutes. Thus, the evolution of ice sheets is modeled with better accuracy, which is critical for a good mitigation of climate change and, in particular, sea-level rise.
Daniel Moreno-Parada, Jorge Alvarez-Solas, Javier Blasco, Marisa Montoya, and Alexander Robinson
The Cryosphere, 17, 2139–2156, https://doi.org/10.5194/tc-17-2139-2023, https://doi.org/10.5194/tc-17-2139-2023, 2023
Short summary
Short summary
We have reconstructed the Laurentide Ice Sheet, located in North America during the Last Glacial Maximum (21 000 years ago). The absence of direct measurements raises a number of uncertainties. Here we study the impact of different physical laws that describe the friction as the ice slides over its base. We found that the Laurentide Ice Sheet is closest to prior reconstructions when the basal friction takes into account whether the base is frozen or thawed during its motion.
Matteo Willeit, Andrey Ganopolski, Alexander Robinson, and Neil R. Edwards
Geosci. Model Dev., 15, 5905–5948, https://doi.org/10.5194/gmd-15-5905-2022, https://doi.org/10.5194/gmd-15-5905-2022, 2022
Short summary
Short summary
In this paper we present the climate component of the newly developed fast Earth system model CLIMBER-X. It has a horizontal resolution of 5°x5° and is designed to simulate the evolution of the Earth system on temporal scales ranging from decades to >100 000 years. CLIMBER-X is available as open-source code and is expected to be a useful tool for studying past climate changes and for the investigation of the long-term future evolution of the climate.
Alexander Robinson, Daniel Goldberg, and William H. Lipscomb
The Cryosphere, 16, 689–709, https://doi.org/10.5194/tc-16-689-2022, https://doi.org/10.5194/tc-16-689-2022, 2022
Short summary
Short summary
Here we investigate the numerical stability of several commonly used methods in order to determine which of them are capable of resolving the complex physics of the ice flow and are also computationally efficient. We find that the so-called DIVA solver outperforms the others. Its representation of the physics is consistent with more complex methods, while it remains computationally efficient at high resolution.
Andreas Born and Alexander Robinson
The Cryosphere, 15, 4539–4556, https://doi.org/10.5194/tc-15-4539-2021, https://doi.org/10.5194/tc-15-4539-2021, 2021
Short summary
Short summary
Ice penetrating radar reflections from the Greenland ice sheet are the best available record of past accumulation and how these layers have been deformed over time by the flow of ice. Direct simulations of this archive hold great promise for improving our models and for uncovering details of ice sheet dynamics that neither models nor data could achieve alone. We present the first three-dimensional ice sheet model that explicitly simulates individual layers of accumulation and how they deform.
Javier Blasco, Jorge Alvarez-Solas, Alexander Robinson, and Marisa Montoya
The Cryosphere, 15, 215–231, https://doi.org/10.5194/tc-15-215-2021, https://doi.org/10.5194/tc-15-215-2021, 2021
Short summary
Short summary
During the Last Glacial Maximum the Antarctic Ice Sheet was larger and more extended than at present. However, neither its exact position nor the total ice volume are well constrained. Here we investigate how the different climatic boundary conditions, as well as basal friction configurations, affect the size and extent of the Antarctic Ice Sheet and discuss its potential implications.
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.
Jorge Alvarez-Solas, Marisa Montoya, and Alexander Robinson
Clim. Past Discuss., https://doi.org/10.5194/cp-2019-96, https://doi.org/10.5194/cp-2019-96, 2019
Publication in CP not foreseen
Short summary
Short summary
Modelling the past abrupt climate changes often resorts to the use of freshwater flux (FWF) in the North Atlantic as an effective method to cause reorganizations of the Atlantic Meridional Overturning Circulation. This procedure has allowed to reproduce the timing of the events. However, the required FWF is inconsistent with reconstructions. Conversely, using a forcing derived from the sea-level record results in a poor fit with the data, highlighting the need of exploring other mechanisms.
Ilaria Tabone, Alexander Robinson, Jorge Alvarez-Solas, and Marisa Montoya
The Cryosphere, 13, 1911–1923, https://doi.org/10.5194/tc-13-1911-2019, https://doi.org/10.5194/tc-13-1911-2019, 2019
Short summary
Short summary
Recent reconstructions show that the North East Greenland Ice Stream (NEGIS) retreated away from its present-day position by 20–40 km during MIS-3. Atmospheric and external forcings were proposed as potential causes of this retreat, but the role of the ocean was not considered. Here, using a 3-D ice-sheet model, we suggest that oceanic warming is sufficient to induce a retreat of the NEGIS margin of many tens of kilometres during MIS-3, helping to explain this conundrum.
Jorge Alvarez-Solas, Rubén Banderas, Alexander Robinson, and Marisa Montoya
Clim. Past, 15, 957–979, https://doi.org/10.5194/cp-15-957-2019, https://doi.org/10.5194/cp-15-957-2019, 2019
Short summary
Short summary
The last glacial period was marked by the existence of of abrupt climatic changes; it is generally accepted that the presence of ice sheets played an important role in their occurrence. While an important effort has been made to investigate the dynamics and evolution of the Laurentide ice sheet during this period, the Eurasian ice sheet (EIS) has not received much attention. Here we investigate the response of the EIS to millennial-scale climate variability using a hybrid 3-D ice-sheet model.
Ilaria Tabone, Alexander Robinson, Jorge Alvarez-Solas, and Marisa Montoya
Clim. Past, 15, 593–609, https://doi.org/10.5194/cp-15-593-2019, https://doi.org/10.5194/cp-15-593-2019, 2019
Short summary
Short summary
By using a 3-D hybrid ice-sheet–shelf model, we investigate the impact of millennial-scale oceanic variability on the Greenland Ice Sheet (GrIS) evolution during the last glacial period (LGP). We show that the GrIS may have strongly reacted to oceanic temperature fluctuations associated with Dansgaard–Oeschger cycles, contributing to sea-level variations of more than 1 m. Our results open the chance for a non-negligible role of the GrIS in millennial-scale oceanic reorganisations during the LGP.
Javier Blasco, Ilaria Tabone, Jorge Alvarez-Solas, Alexander Robinson, and Marisa Montoya
Clim. Past, 15, 121–133, https://doi.org/10.5194/cp-15-121-2019, https://doi.org/10.5194/cp-15-121-2019, 2019
Short summary
Short summary
The LGP is a period punctuated by the presence of several abrupt climate events and sea-level variations of up to 20 m at millennial timescales. The origin of those fluctuations is attributed to NH paleo ice sheets, but a contribution from the AIS cannot be excluded. Here, for the first time, we investigate the response of the AIS to millennial climate variability using an ice sheet–shelf model. We shows that the AIS produces substantial sea-level rises and grounding line migrations.
Rubén Banderas, Jorge Alvarez-Solas, Alexander Robinson, and Marisa Montoya
Geosci. Model Dev., 11, 2299–2314, https://doi.org/10.5194/gmd-11-2299-2018, https://doi.org/10.5194/gmd-11-2299-2018, 2018
Short summary
Short summary
Here we present a new approach to force ice-sheet models offline, which accounts for a more realistic treatment of millennial-scale climate variability as compared to the existing methods. Our results reveal that an incorrect representation of the characteristic pattern of millennial-scale climate variability within the climate forcing not only affects NH ice-volume variations at millennial timescales but has consequences for glacial–interglacial ice-volume changes too.
Ilaria Tabone, Javier Blasco, Alexander Robinson, Jorge Alvarez-Solas, and Marisa Montoya
Clim. Past, 14, 455–472, https://doi.org/10.5194/cp-14-455-2018, https://doi.org/10.5194/cp-14-455-2018, 2018
Short summary
Short summary
The response of the Greenland Ice Sheet (GrIS) to palaeo-oceanic changes on a glacial–interglacial timescale is studied from a modelling perspective. A 3-D hybrid ice-sheet–shelf model which includes a parameterization of the basal melting rate at the GrIS marine margins is used. The results show that the oceanic forcing plays a key role in the GrIS evolution, not only by controlling the ice retreat during the deglaciation but also by driving the ice expansion in glacial periods.
Jorge Alvarez-Solas, Rubén Banderas, Alexander Robinson, and Marisa Montoya
Clim. Past Discuss., https://doi.org/10.5194/cp-2017-143, https://doi.org/10.5194/cp-2017-143, 2017
Revised manuscript not accepted
Short summary
Short summary
The last glacial period was marked by the existence of of abrupt climatic changes. It is generally accepted that the presence of ice sheets played an
important role in their occurrence. While an important effort has been made to investigate the dynamics and evolution of the Laurentide Ice Sheet during this period, the Eurasian Ice Sheet (EIS) has not received much attention. Here we investigate the response of the EIS to millennial-scale climate variability. We use a hybrid 3D ice-sheet model.
Mario Krapp, Alexander Robinson, and Andrey Ganopolski
The Cryosphere, 11, 1519–1535, https://doi.org/10.5194/tc-11-1519-2017, https://doi.org/10.5194/tc-11-1519-2017, 2017
Short summary
Short summary
We present the snowpack model SEMIC. It calculates snow height, surface temperature, surface albedo, and the surface mass balance of snow- and ice-covered surfaces while using meteorological data as input. In this paper we describe how SEMIC works and how well it compares with snowpack data of a more sophisticated regional climate model applied to the Greenland ice sheet. Because of its simplicity and efficiency, SEMIC can be used as a coupling interface between atmospheric and ice sheet models.
A. Robinson and M. Perrette
Geosci. Model Dev., 8, 1877–1883, https://doi.org/10.5194/gmd-8-1877-2015, https://doi.org/10.5194/gmd-8-1877-2015, 2015
Short summary
Short summary
Here we present a concise interface to the NetCDF library designed to simplify reading and writing tasks of up to 6-D arrays in Fortran programs.
R. Calov, A. Robinson, M. Perrette, and A. Ganopolski
The Cryosphere, 9, 179–196, https://doi.org/10.5194/tc-9-179-2015, https://doi.org/10.5194/tc-9-179-2015, 2015
Short summary
Short summary
Ice discharge into the ocean from outlet glaciers is an important
component of mass loss of the Greenland ice sheet. Here, we present a
simple parameterization of ice discharge for coarse resolution ice
sheet models, suitable for large ensembles or long-term palaeo
simulations. This parameterization reproduces in a good approximation
the present-day ice discharge compared with estimates, and the
simulation of the present-day ice sheet elevation is considerably
improved.
A. Robinson and H. Goelzer
The Cryosphere, 8, 1419–1428, https://doi.org/10.5194/tc-8-1419-2014, https://doi.org/10.5194/tc-8-1419-2014, 2014
Related subject area
Cryosphere
Computationally efficient subglacial drainage modelling using Gaussian process emulators: GlaDS-GP v1.0
Anisotropic metric-based mesh adaptation for ice flow modelling in Firedrake
The Utrecht Finite Volume Ice-Sheet Model (UFEMISM) version 2.0 – Part 1: Description and idealised experiments
A Flexible Snow Model (FSM 2.1.1) including a forest canopy
CMIP6 models overestimate sea ice melt, growth and conduction relative to ice mass balance buoy estimates
Coupling framework (1.0) for the Úa (2023b) ice sheet model and the FESOM-1.4 z-coordinate ocean model in an Antarctic domain
A gradient-boosted tree framework to model the ice thickness of the world's glaciers (IceBoost v1.1)
Comparison of calibration methods of a PICO basal ice shelf melt module implemented in the GRISLI v2.0 ice sheet model
A Python library for solving ice sheet modeling problems using Physics Informed Neural Networks, PINNICLE v1.0
Towards deep-learning solutions for classification of automated snow height measurements (CleanSnow v1.0.2)
Constraining CMIP6 sea ice simulations with ICESat-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
Implementation of an intermediate complexity snow-physics scheme (ISBA-Explicit Snow) into a sea-ice model (SI3): 1D thermodynamic coupling and validation
A framework for three-dimensional dynamic modeling of mountain glaciers in the Community Ice Sheet Model (CISM v2.2)
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
Design and performance of ELSA v2.0: an isochronal model for ice-sheet layer tracing
Southern Ocean Ice Prediction System version 1.0 (SOIPS v1.0): description of the system and evaluation of synoptic-scale sea ice forecasts
Lagrangian tracking of sea ice in Community Ice CodE (CICE; version 5)
openAMUNDSEN v1.0: an open-source snow-hydrological model for mountain regions
OpenFOAM-avalanche 2312: depth-integrated models beyond dense-flow avalanches
Refactoring the elastic–viscous–plastic solver from the sea ice model CICE v6.5.1 for improved performance
A new 3D full-Stokes calving algorithm within Elmer/Ice (v9.0)
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)
Tim Hill, Derek Bingham, Gwenn E. Flowers, and Matthew J. Hoffman
Geosci. Model Dev., 18, 4045–4074, https://doi.org/10.5194/gmd-18-4045-2025, https://doi.org/10.5194/gmd-18-4045-2025, 2025
Short summary
Short summary
Subglacial drainage models represent water flow beneath glaciers and ice sheets. Here, we train fast statistical models called Gaussian process (GP) 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.
Davor Dundovic, Joseph G. Wallwork, Stephan C. Kramer, Fabien Gillet-Chaulet, Regine Hock, and Matthew D. Piggott
Geosci. Model Dev., 18, 4023–4044, https://doi.org/10.5194/gmd-18-4023-2025, https://doi.org/10.5194/gmd-18-4023-2025, 2025
Short summary
Short summary
Accurate numerical studies of glaciers often require high-resolution simulations, which often prove too demanding even for modern computers. In this paper we develop a method that identifies whether different parts of a glacier require high or low resolution based on its physical features, such as its thickness and velocity. We show that by doing so we can achieve a more optimal simulation accuracy for the available computing resources compared to uniform-resolution simulations.
Constantijn J. Berends, Victor Azizi, Jorge A. Bernales, and Roderik S. W. van de Wal
Geosci. Model Dev., 18, 3635–3659, https://doi.org/10.5194/gmd-18-3635-2025, https://doi.org/10.5194/gmd-18-3635-2025, 2025
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 that 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.
Richard Essery, Giulia Mazzotti, Sarah Barr, Tobias Jonas, Tristan Quaife, and Nick Rutter
Geosci. Model Dev., 18, 3583–3605, https://doi.org/10.5194/gmd-18-3583-2025, https://doi.org/10.5194/gmd-18-3583-2025, 2025
Short summary
Short summary
How forests influence accumulation and melt of snow on the ground is of long-standing interest, but uncertainty remains in how best to model forest snow processes. We developed the Flexible Snow Model version 2 to quantify these uncertainties. In a first model demonstration, how unloading of intercepted snow from the forest canopy is represented is responsible for the largest uncertainty. Global mapping of forest distribution is also likely to be a large source of uncertainty in existing models.
Alex E. West and Edward W. Blockley
Geosci. Model Dev., 18, 3041–3064, https://doi.org/10.5194/gmd-18-3041-2025, https://doi.org/10.5194/gmd-18-3041-2025, 2025
Short summary
Short summary
This study uses ice mass balance buoys – temperature- and height-measuring devices frozen into sea ice – to find how well climate models simulate (1) melt and growth of Arctic sea ice and (2) conduction of heat through Arctic sea ice. This may help understand why models produce varying amounts of sea ice in the present day. We find that models tend to show more melt, growth or conduction for a given ice thickness than the buoys, although the difference is smaller for models with more physically realistic thermodynamics.
Ole Richter, Ralph Timmermann, G. Hilmar Gudmundsson, and Jan De Rydt
Geosci. Model Dev., 18, 2945–2960, https://doi.org/10.5194/gmd-18-2945-2025, https://doi.org/10.5194/gmd-18-2945-2025, 2025
Short summary
Short summary
The new coupled ice sheet–ocean model addresses challenges related to horizontal resolution through advanced mesh flexibility, enabled by the use of unstructured grids. We describe the new model, verify its functioning in an idealised setting and demonstrate its advantages in a global-ocean–Antarctic ice sheet domain. The results of this study comprise an important step towards improving predictions of the Antarctic contribution to sea level rise over centennial timescales.
Niccolò Maffezzoli, Eric Rignot, Carlo Barbante, Troels Petersen, and Sebastiano Vascon
Geosci. Model Dev., 18, 2545–2568, https://doi.org/10.5194/gmd-18-2545-2025, https://doi.org/10.5194/gmd-18-2545-2025, 2025
Short summary
Short summary
In this work we introduce IceBoost, a machine learning framework to model the ice thickness distribution of all the world's glaciers with greater accuracy than state-of-the-art methods. The model is trained on 3.7 million measurements globally available and provides skilful estimates across all regions. This advancement will help in better assessing future sea level changes and freshwater resources, with significance for both the scientific community and society at large.
Maxence Menthon, Pepijn Bakker, Aurélien Quiquet, Didier M. Roche, and Ronja Reese
EGUsphere, https://doi.org/10.5194/egusphere-2025-777, https://doi.org/10.5194/egusphere-2025-777, 2025
Short summary
Short summary
The ice-ocean interaction is a large source of uncertainty in future projections of the Antarctic ice sheet. Here we implement a basal ice shelf melt module (PICO) in a ice sheet model (GRISLI) and test six simple statistical methods to calibrate this module. We show that calculating the mean absolute error of bins best fits the observational datasets under multiple conditions. We also assess the impact of the module implementation and calibration choice on future projections until 2300.
Gong Cheng, Mansa Krishna, and Mathieu Morlighem
EGUsphere, https://doi.org/10.5194/egusphere-2025-1188, https://doi.org/10.5194/egusphere-2025-1188, 2025
Short summary
Short summary
Predicting ice sheet contributions to sea level rise is challenging due to limited data and uncertainties in key processes. Traditional models require complex methods that lack flexibility. We developed PINNICLE, an open-source Python library that integrates machine learning with physical laws to improve ice sheet modeling. By combining data and physics, PINNICLE enhances predictions and adaptability, providing a powerful tool for climate research and sea level rise projections.
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.
Alek Aaron Petty, Christopher Cardinale, and Madison Smith
EGUsphere, https://doi.org/10.5194/egusphere-2025-766, https://doi.org/10.5194/egusphere-2025-766, 2025
Short summary
Short summary
We put global climate models to the test against NASA’s ICESat-2 satellite to see how well they simulate global sea ice cover. By adding fancy laser data from ICESat-2, we can better assess how well the models are performing compared to the standard assessments of sea ice area. Overall the models do a good job but there’s room for improvement, especially across the Southern Ocean. We should think a bit more about sea ice density if we want more reliable freeboard comparisons.
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.
Théo Brivoal, Virginie Guemas, Martin Vancoppenolle, Clément Rousset, and Bertrand Decharme
EGUsphere, https://doi.org/10.5194/egusphere-2024-3220, https://doi.org/10.5194/egusphere-2024-3220, 2025
Short summary
Short summary
Snow in polar regions is key to sea ice formation and the Earth's climate, but current climate models simplify snow cover on sea ice. This study integrates an intermediate complexity snow-physics scheme into a sea-ice model designed for climate applications. We show that modelling the temporal changes in properties such as the density and thermal conductivity of the snow layers leads to a more accurate representation of heat transfer between the underlying sea ice and the atmosphere.
Samar Minallah, William Lipscomb, Gunter Leguy, and Harry Zekollari
EGUsphere, https://doi.org/10.5194/egusphere-2024-4152, https://doi.org/10.5194/egusphere-2024-4152, 2025
Short summary
Short summary
We implemented a new modeling framework within an Earth system model to study the evolution of mountain glaciers under different climate scenarios and applied it to the European Alps. Alpine glaciers will lose a large volume fraction under current temperatures, with near complete ice loss under warmer scenarios. This is the first use of a 3D, higher-order ice flow model for regional-scale glacier simulations that will enable assessments of coupled land ice and Earth system processes.
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.
Therese Rieckh, Andreas Born, Alexander Robinson, Robert Law, and Gerrit Gülle
Geosci. Model Dev., 17, 6987–7000, https://doi.org/10.5194/gmd-17-6987-2024, https://doi.org/10.5194/gmd-17-6987-2024, 2024
Short summary
Short summary
We present the open-source model ELSA, which simulates the internal age structure of large ice sheets. It creates layers of snow accumulation at fixed times during the simulation, which are used to model the internal stratification of the ice sheet. Together with reconstructed isochrones from radiostratigraphy data, ELSA can be used to assess ice sheet models and to improve their parameterization. ELSA can be used coupled to an ice sheet model or forced with its output.
Fu Zhao, Xi Liang, Zhongxiang Tian, Ming Li, Na Liu, and Chengyan Liu
Geosci. Model Dev., 17, 6867–6886, https://doi.org/10.5194/gmd-17-6867-2024, https://doi.org/10.5194/gmd-17-6867-2024, 2024
Short summary
Short summary
In this work, we introduce a newly developed Antarctic sea ice forecasting system, namely the Southern Ocean Ice Prediction System (SOIPS). The system is based on a regional sea ice‒ocean‒ice shelf coupled model and can assimilate sea ice concentration observations. By assessing the system's performance in sea ice forecasts, we find that the system can provide reliable Antarctic sea ice forecasts for the next 7 d and has the potential to guide ship navigation in the Antarctic sea ice zone.
Chenhui Ning, Shiming Xu, Yan Zhang, Xuantong Wang, Zhihao Fan, and Jiping Liu
Geosci. Model Dev., 17, 6847–6866, https://doi.org/10.5194/gmd-17-6847-2024, https://doi.org/10.5194/gmd-17-6847-2024, 2024
Short summary
Short summary
Sea ice models are mainly based on non-moving structured grids, which is different from buoy measurements that follow the ice drift. To facilitate Lagrangian analysis, we introduce online tracking of sea ice in Community Ice CodE (CICE). We validate the sea ice tracking with buoys and evaluate the sea ice deformation in high-resolution simulations, which show multi-fractal characteristics. The source code is openly available and can be used in various scientific and operational applications.
Ulrich Strasser, Michael Warscher, Erwin Rottler, and Florian Hanzer
Geosci. Model Dev., 17, 6775–6797, https://doi.org/10.5194/gmd-17-6775-2024, https://doi.org/10.5194/gmd-17-6775-2024, 2024
Short summary
Short summary
openAMUNDSEN is a fully distributed open-source snow-hydrological model for mountain catchments. It includes process representations of an empirical, semi-empirical, and physical nature. It uses temperature, precipitation, humidity, radiation, and wind speed as forcing data and is computationally efficient, of a modular nature, and easily extendible. The Python code is available on GitHub (https://github.com/openamundsen/openamundsen), including documentation (https://doc.openamundsen.org).
Matthias Rauter and Julia Kowalski
Geosci. Model Dev., 17, 6545–6569, https://doi.org/10.5194/gmd-17-6545-2024, https://doi.org/10.5194/gmd-17-6545-2024, 2024
Short summary
Short summary
Snow avalanches can form large powder clouds that substantially exceed the velocity and reach of the dense core. Only a few complex models exist to simulate this phenomenon, and the respective hazard is hard to predict. This work provides a novel flow model that focuses on simple relations while still encapsulating the significant behaviour. The model is applied to reconstruct two catastrophic powder snow avalanche events in Austria.
Till Andreas Soya Rasmussen, Jacob Poulsen, Mads Hvid Ribergaard, Ruchira Sasanka, Anthony P. Craig, Elizabeth C. Hunke, and Stefan Rethmeier
Geosci. Model Dev., 17, 6529–6544, https://doi.org/10.5194/gmd-17-6529-2024, https://doi.org/10.5194/gmd-17-6529-2024, 2024
Short summary
Short summary
Earth system models (ESMs) today strive for better quality based on improved resolutions and improved physics. A limiting factor is the supercomputers at hand and how best to utilize them. This study focuses on the refactorization of one part of a sea ice model (CICE), namely the dynamics. It shows that the performance can be significantly improved, which means that one can either run the same simulations much cheaper or advance the system according to what is needed.
Iain Wheel, Douglas I. Benn, Anna J. Crawford, Joe Todd, and Thomas Zwinger
Geosci. Model Dev., 17, 5759–5777, https://doi.org/10.5194/gmd-17-5759-2024, https://doi.org/10.5194/gmd-17-5759-2024, 2024
Short summary
Short summary
Calving, the detachment of large icebergs from glaciers, is one of the largest uncertainties in future sea level rise projections. This process is poorly understood, and there is an absence of detailed models capable of simulating calving. A new 3D calving model has been developed to better understand calving at glaciers where detailed modelling was previously limited. Importantly, the new model is very flexible. By allowing for unrestricted calving geometries, it can be applied at any location.
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.
Cited articles
Aboelaze, M., Chrisochoides, N., and Houstis, E.: The Parallelization of Level 2 and 3 BLAS Operations on Distributed Memory Machines, Tech. Rep. CSD-TR-91-007, Purdue University, West Lafayette, IN, 1991. a
Alley, R. B., Blankenship, D. D., Bentley, C. R., and Rooney, S. T.: Deformation of till beneath ice stream B, West Antarctica, Nature, 322, 57–59, https://doi.org/10.1038/322057a0, 1986. a
Arakawa, A. and Lamb, V. R.: Computational Design of the Basic Dynamical Processes of the UCLA General Circulation Model, in: 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
Arthern, R. J. and Williams, C. R.: The sensitivity of West Antarctica to the submarine melting feedback, Geophys. Res. Lett., 44, 2352–2359, https://doi.org/10.1002/2017gl072514, 2017. a
Arthern, R. J., Hindmarsh, R. C. A., and Williams, C. R.: Flow speed within the Antarctic ice sheet and its controls inferred from satellite observations, J. Geophys. Res.-Earth, 120, 1171–1188, https://doi.org/10.1002/2014jf003239, 2015. a
Asay-Davis, X. S., Jourdain, N. C., and Nakayama, Y.: Developments in Simulating and Parameterizing Interactions Between the Southern Ocean and the Antarctic Ice Sheet, Current Climate Change Reports, 3, 316–329, https://doi.org/10.1007/s40641-017-0071-0, 2017. a
Aschwanden, A., Aðalgeirsdóttir, G., and Khroulev, C.: Hindcasting to measure ice sheet model sensitivity to initial states, The Cryosphere, 7, 1083–1093, https://doi.org/10.5194/tc-7-1083-2013, 2013. a
Bamber, J. L., Vaughan, D. G., and Joughin, I.: Widespread Complex Flow in the Interior of the Antarctic Ice Sheet, Science, 287, 1248–1250, https://doi.org/10.1126/science.287.5456.1248, 2000. a
Bamber, J. L., Riva, R. E. M., Vermeersen, B. L. A., and LeBrocq, A. M.: Reassessment of the Potential Sea-Level Rise from a Collapse of the West Antarctic Ice Sheet, Science, 324, 901–903, https://doi.org/10.1126/science.1169335, 2009. a, b
Bassis, J. and Jacobs, S.: Diverse calving patterns linked to glacier geometry, Nat. Geosci., 6, 833–836, https://doi.org/10.1038/ngeo1887, 2013. a
Bassis, J. N. and Walker, C. C.: Upper and lower limits on the stability of calving glaciers from the yield strength envelope of ice, P. R. Soc. A, 468, 913–931, https://doi.org/10.1098/rspa.2011.0422, 2012. a
Bassis, J. N., Petersen, S. V., and Cathles, L. M.: Heinrich events triggered by ocean forcing and modulated by isostatic adjustment, Nature, 542, 332–334, https://doi.org/10.1038/nature21069, 2017. a, b, c
Bentley, C. R.: Rapid sea-level rise from a West Antarctic ice-sheet collapse: a short-term perspective, J. Glaciol., 44, 157–163, https://doi.org/10.3189/s0022143000002446, 1998. a, b
Blankenship, D. D., Bentley, C. R., Rooney, S. T., and Alley, R. B.: Seismic measurements reveal a saturated porous layer beneath an active Antarctic ice stream, Nature, 322, 54–57, https://doi.org/10.1038/322054a0, 1986. a
Blatter, H.: Velocity and stress fields in grounded glaciers: a simple algorithm for including deviatoric stress gradients, J. Glaciol., 41, 333–344, https://doi.org/10.3189/s002214300001621x, 1995. a, b
Bougamont, M. and Tulaczyk, S.: Glacial erosion beneath ice streams and ice-stream tributaries: constraints on temporal and spatial distribution of erosion from numerical simulations of a West Antarctic ice stream, Boreas, 32, 178–190, https://doi.org/10.1111/j.1502-3885.2003.tb01436.x, 2003. a
Bougamont, M., Christoffersen, P., and Tulaczyk, S.: Response of subglacial sediments to basal freeze-on 1. Theory and comparison to observations from beneath the West Antarctic Ice Sheet, J. Geophys. Res.-Sol. Ea., 108, 2222, https://doi.org/10.1029/2002jb001935, 2003. a
Brown, S. A., Folk, M., Goucher, G., Rew, R., and Dubois, P. F.: Software for Portable Scientific Data Management, Comput. Phys., 7, 304–308, https://doi.org/10.1063/1.4823180, 1993. a
Bueler, E. and Brown, J.: Shallow shelf approximation as a “sliding law” in a thermomechanically coupled ice sheet model, J. Geophys. Res., 114, F03008, https://doi.org/10.1029/2008JF001179, 2009. a, b, c
Catalyurek, U. and Aykanat, C.: Hypergraph-partitioning-based decomposition for parallel sparse-matrix vector multiplication, IEEE T. Parall. Distr., 10, 673–693, https://doi.org/10.1109/71.780863, 1999. a
Christian, J. E., Robel, A. A., and Catania, G.: A probabilistic framework for quantifying the role of anthropogenic climate change in marine-terminating glacier retreats, The Cryosphere, 16, 2725–2743, https://doi.org/10.5194/tc-16-2725-2022, 2022. a
Chugunov, V. A. and Wilchinsky, A. V.: Modelling of a marine glacier and ice-sheet-ice-shelf transition zone based on asymptotic analysis, Ann. Glaciol., 23, 59–67, https://doi.org/10.3189/s0260305500013264, 1996. a
DeConto, R. M., Pollard, D., Alley, R. B., Velicogna, I., Gasson, E., Gomez, N., Sadai, S., Condron, A., Gilford, D. M., Ashe, E. L., Kopp, R. E., Li, D., and Dutton, A.: The Paris Climate Agreement and future sea-level rise from Antarctica, Nature, 593, 83–89, https://doi.org/10.1038/s41586-021-03427-0, 2021. a
Donat-Magnin, M., Jourdain, N. C., Spence, P., Le Sommer, J., Gallée, H., and Durand, G.: Ice-shelf melt response to changing winds and glacier dynamics in the Amundsen Sea Sector, Antarctica, J. Geophys. Res.-Oceans, 122, 10206–10224, https://doi.org/10.1002/2017JC013059, 2017. a
Dukowicz, J. K., Price, S. F., and Lipscomb, W. H.: Consistent approximations and boundary conditions for ice-sheet dynamics from a principle of least action, J. Glaciol., 56, 480–496, https://doi.org/10.3189/002214310792447851, 2010. a
Dupont, T. K. and Alley, R. B.: Assessment of the importance of ice-shelf buttressing to ice-sheet flow, Geophys. Res. Lett., 32, L04503, https://doi.org/10.1029/2004gl022024, 2005. a
Engelhardt, H., Humphrey, N., Kamb, B., and Fahnestock, M.: Physical Conditions at the Base of a Fast Moving Antarctic Ice Stream, Science, 248, 57–59, https://doi.org/10.1126/science.248.4951.57, 1990. a
Favier, L., Jourdain, N. C., Jenkins, A., Merino, N., Durand, G., Gagliardini, O., Gillet-Chaulet, F., and Mathiot, P.: Assessment of sub-shelf melting parameterisations using the ocean–ice-sheet coupled model NEMO(v3.6)–Elmer/Ice(v8.3), Geosci. Model Dev., 12, 2255–2283, https://doi.org/10.5194/gmd-12-2255-2019, 2019. a, b, c, d, e, f, g, h
Feldmann, J. and Levermann, A.: Collapse of the West Antarctic Ice Sheet after local destabilization of the Amundsen Basin, P. Natl. Acad. Sci. USA, 112, 14191–14196, https://doi.org/10.1073/pnas.1512482112, 2015. a, b
Fischler, Y., Rückamp, M., Bischof, C., Aizinger, V., Morlighem, M., and Humbert, A.: A scalability study of the Ice-sheet and Sea-level System Model (ISSM, version 4.18), Geosci. Model Dev., 15, 3753–3771, https://doi.org/10.5194/gmd-15-3753-2022, 2022. 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
Gagliardini, O., Cohen, D., Råback, P., and Zwinger, T.: Finite-element modeling of subglacial cavities and related friction law, J. Geophys. Res., 112, F02027, https://doi.org/10.1029/2006jf000576, 2007. a
Gagliardini, O., Zwinger, T., Gillet-Chaulet, F., Durand, G., Favier, L., de Fleurian, B., Greve, R., Malinen, M., Martín, C., Råback, P., Ruokolainen, J., Sacchettini, M., Schäfer, M., Seddik, H., and Thies, J.: Capabilities and performance of Elmer/Ice, a new-generation ice sheet model, Geosci. Model Dev., 6, 1299–1318, https://doi.org/10.5194/gmd-6-1299-2013, 2013. a
Garbe, J., Albrecht, T., Levermann, A., Donges, J. F., and Winkelmann, R.: The hysteresis of the Antarctic Ice Sheet, Nature, 585, 538–544, https://doi.org/10.1038/s41586-020-2727-5, 2020. a, b
Glen, J. W.: The creep of polycrystalline ice, P. R. Soc. Lond. A, 228, 519–538, https://doi.org/10.1098/rspa.1955.0066, 1995. a
Goldberg, D. N.: A variationally derived, depth-integrated approximation to a higher-order glaciological flow model, J. Glaciol., 57, 157–170, https://doi.org/10.3189/002214311795306763, 2011. a, b
Graham, R. L., Woodall, T. S., and Squyres, J. M.: Open MPI: A Flexible High Performance MPI, Springer, Berlin, Heidelberg, 228–239, https://doi.org/10.1007/11752578_29, 2006. a
Greve, R. and Blatter, H.: Dynamics of Ice Sheets and Glaciers, Springer Berlin Heidelberg, https://doi.org/10.1007/978-3-642-03415-2, 2009. a, b, c
Grosfeld, K., Gerdes, R., and Determann, J.: Thermohaline circulation and interaction between ice shelf cavities and the adjacent open ocean, J. Geophys. Res.-Oceans, 102, 15595–15610, https://doi.org/10.1029/97jc00891, 1997. a
Haseloff, M. and Sergienko, O. V.: The effect of buttressing on grounding line dynamics, J. Glaciol., 64, 417–431, https://doi.org/10.1017/jog.2018.30, 2018. a
Hill, E. A., Urruty, B., Reese, R., Garbe, J., Gagliardini, O., Durand, G., Gillet-Chaulet, F., Gudmundsson, G. H., Winkelmann, R., Chekki, M., Chandler, D., and Langebroek, P. M.: The stability of present-day Antarctic grounding lines – Part 1: No indication of marine ice sheet instability in the current geometry, The Cryosphere, 17, 3739–3759, https://doi.org/10.5194/tc-17-3739-2023, 2023. a, b
Hindmarsh, R. C.: The role of membrane-like stresses in determining the stability and sensitivity of the Antarctic ice sheets: back pressure and grounding line motion, Philos. T. R. Soc. A, 364, 1733–1767, https://doi.org/10.1098/rsta.2006.1797, 2006. a
Hindmarsh, R. C. and le Meur, E.: Dynamical processes involved in the retreat of marine ice sheets, J. Glaciol., 47, 271–282, https://doi.org/10.3189/172756501781832269, 2001. a
Hindmarsh, R. C. A.: Qualitative Dynamics of Marine Ice Sheets, in: Ice in the Climate System, Springer Berlin Heidelberg, 67–99, https://doi.org/10.1007/978-3-642-85016-5_5, 1993. a, b, c
Hindmarsh, R. C. A.: Stability of ice rises and uncoupled marine ice sheets, Ann. Glaciol., 23, 105–115, https://doi.org/10.3189/s0260305500013318, 1996. 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
Holland, P. R., Jenkins, A., and Holland, D. M.: The Response of Ice Shelf Basal Melting to Variations in Ocean Temperature, J. Climate, 21, 2558–2572, https://doi.org/10.1175/2007JCLI1909.1, 2008. a, b
Jacob, B. and Guennebaud, G.: Eigen v3, http://eigen.tuxfamily.org (last access: 27 March 2025), 2010. a
Jamieson, S. S. R., Vieli, A., Livingstone, S. J., Cofaigh, C. O., Stokes, C., Hillenbrand, C.-D., and Dowdeswell, J. A.: Ice-stream stability on a reverse bed slope, Nat. Geosci., 5, 799–802, https://doi.org/10.1038/ngeo1600, 2012. a, b
Jenkins, A., Shoosmith, D., Dutrieux, P., Jacobs, S., Kim, T. W., Lee, S. H., Ha, H. K., and Stammerjohn, S.: West Antarctic Ice Sheet retreat in the Amundsen Sea driven by decadal oceanic variability, Nat. Geosci., 11, 733–738, https://doi.org/10.1038/s41561-018-0207-4, 2018. a
Joughin, I., MacAyeal, D. R., and Tulaczyk, S.: Basal shear stress of the Ross ice streams from control method inversions, J. Geophys. Res.-Sol. Ea., 109, B09405, https://doi.org/10.1029/2003jb002960, 2004. a
Joughin, I., Smith, B. E., and Schoof, C. G.: Regularized Coulomb Friction Laws for Ice Sheet Sliding: Application to Pine Island Glacier, Antarctica, Geophys. Res. Lett., 46, 4764–4771, https://doi.org/10.1029/2019gl082526, 2019. a, b
Joughin, I., Shapero, D., Dutrieux, P., and Smith, B.: Ocean-induced melt volume directly paces ice loss from Pine Island Glacier, Science Advances, 7, eabi5738, https://doi.org/10.1126/sciadv.abi5738, 2021. a, b
Lipscomb, W. H., Price, S. F., Hoffman, M. J., Leguy, G. R., Bennett, A. R., Bradley, S. L., Evans, K. J., Fyke, J. G., Kennedy, J. H., Perego, M., Ranken, D. M., Sacks, W. J., Salinger, A. G., Vargo, L. J., and Worley, P. H.: Description and evaluation of the Community Ice Sheet Model (CISM) v2.1, Geosci. Model Dev., 12, 387–424, https://doi.org/10.5194/gmd-12-387-2019, 2019. a, b
Ma, Y., Gagliardini, O., Ritz, C., Gillet-Chaulet, F., Durand, G., and Montagnat, M.: Enhancement factors for grounded ice and ice shelves inferred from an anisotropic ice-flow model, J. Glaciol., 56, 805–812, https://doi.org/10.3189/002214310794457209, 2010. a
MacAyeal, D. R. and Barcilon, V.: Ice-shelf Response to Ice-stream Discharge Fluctuations: I. Unconfined Ice Tongues, J. Glaciol., 34, 121–127, https://doi.org/10.3189/s002214300000914x, 1988. a
Martin, D. F., Cornford, S. L., and Payne, A. J.: Millennial-Scale Vulnerability of the Antarctic Ice Sheet to Regional Ice Shelf Collapse, Geophys. Res. Lett., 46, 1467–1475, https://doi.org/10.1029/2018gl081229, 2019. a
Meier, M. F. and Post, A.: Fast tidewater glaciers, J. Geophys. Res., 92, 9051–9058, https://doi.org/10.1029/JB092iB09p09051, 1987. a
Minchew, B. M., Meyer, C. R., Robel, A. A., Gudmundsson, G. H., and Simons, M.: Processes controlling the downstream evolution of ice rheology in glacier shear margins: case study on Rutford Ice Stream, West Antarctica, J. Glaciol., 64, 583–594, https://doi.org/10.1017/jog.2018.47, 2018. a
Moreno-Parada, D.: Nix, Zenodo [code], https://doi.org/10.5281/zenodo.14943834, 2025a. a
Moreno-Parada, D.: Nix simulation under oceanic forcing, Zenodo [video], https://doi.org/10.5281/zenodo.15082468, 2025b. a
Moreno-Parada, D.: Description and validation of the ice sheet model Nix v1.0, Zenodo [data set], https://doi.org/10.5281/zenodo.15721846, 2025c. a
Moreno-Parada, D., Robinson, A., Montoya, M., and Alvarez-Solas, J.: Nix ice sheet model v1.0.0, https://doi.org/10.5281/ZENODO.10228874, 2023. a
Moreno-Parada, D., Robinson, A., Montoya, M., and Alvarez-Solas, J.: Analytical solutions for the advective–diffusive ice column in the presence of strain heating, The Cryosphere, 18, 4215–4232, https://doi.org/10.5194/tc-18-4215-2024, 2024. a
Mouginot, J., Rignot, E., and Scheuchl, B.: Sustained increase in ice discharge from the Amundsen Sea Embayment, West Antarctica, from 1973 to 2013, Geophys. Res. Lett., 41, 1576–1584, https://doi.org/10.1002/2013gl059069, 2014. a
Nick, F., Vieli, A., Howat, I., and Joughin, I.: Large-scale changes in Greenland outlet glacier dynamics triggered at the terminus, Nat. Geosci., 2, 110–114, https://doi.org/10.1038/ngeo394, 2009. a
Nick, F., van der Veen, C., Vieli, A., and Benn, D.: A physically based calving model applied to marine outlet glaciers and implications for the glacier dynamics, J. Glaciol., 56, 781–794, https://doi.org/10.3189/002214310794457344, 2010. a
Nickolls, J., Buck, I., Garland, M., and Skadron, K.: Scalable Parallel Programming with CUDA: Is CUDA the parallel programming model that application developers have been waiting for?, Queue, 6, 40–53, https://doi.org/10.1145/1365490.1365500, 2008. a
Nowicki, S. and Wingham, D.: Conditions for a steady ice sheet–ice shelf junction, Earth Planet. Sc. Lett., 265, 246–255, https://doi.org/10.1016/j.epsl.2007.10.018, 2008. a
Nye, J. F.: The distribution of stress and velocity in glaciers and ice-sheets, P. R. Soc. Lond. A, 239, 113–133, https://doi.org/10.1098/rspa.1957.0026, 1957. a
Paolo, F. S., Fricker, H. A., and Padman, L.: Volume loss from Antarctic ice shelves is accelerating, Science, 348, 327–331, https://doi.org/10.1126/science.aaa0940, 2015. a
Pattyn, F.: A new three-dimensional higher-order thermomechanical ice sheet model: Basic sensitivity, ice stream development, and ice flow across subglacial lakes, J. Geophys. Res., 108, 2382, https://doi.org/10.1029/2002jb002329, 2003. a, b, c
Pattyn, F. and Morlighem, M.: The uncertain future of the Antarctic Ice Sheet, Science, 367, 1331–1335, https://doi.org/10.1126/science.aaz5487, 2020. a, b
Pattyn, F., Huyghe, A., Brabander, S. D., and Smedt, B. D.: Role of transition zones in marine ice sheet dynamics, J. Geophys. Res., 111, F02004, https://doi.org/10.1029/2005jf000394, 2006. a, b, c
Pattyn, F., Perichon, L., Aschwanden, A., Breuer, B., de Smedt, B., Gagliardini, O., Gudmundsson, G. H., Hindmarsh, R. C. A., Hubbard, A., Johnson, J. V., Kleiner, T., Konovalov, Y., Martin, C., Payne, A. J., Pollard, D., Price, S., Rückamp, M., Saito, F., Souček, O., Sugiyama, S., and Zwinger, T.: Benchmark experiments for higher-order and full-Stokes ice sheet models (ISMIP–HOM), The Cryosphere, 2, 95–108, https://doi.org/10.5194/tc-2-95-2008, 2008. a
Pattyn, F., Schoof, C., Perichon, L., Hindmarsh, R. C. A., Bueler, E., de Fleurian, B., Durand, G., Gagliardini, O., Gladstone, R., Goldberg, D., Gudmundsson, G. H., Huybrechts, P., Lee, V., Nick, F. M., Payne, A. J., Pollard, D., Rybak, O., Saito, F., and Vieli, A.: Results of the Marine Ice Sheet Model Intercomparison Project, MISMIP, The Cryosphere, 6, 573–588, https://doi.org/10.5194/tc-6-573-2012, 2012. a, b, c, d, e, f, g, h, i, j, k, l, m
Payne, A. J.: Limit cycles in the basal thermal regime of ice sheets, J. Geophys. Res.-Sol. Ea., 100, 4249–4263, https://doi.org/10.1029/94jb02778, 1995. 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. A.: Results from the EISMINT model intercomparison: the effects of thermomechanical coupling, J. Glaciol., 46, 227–238, https://doi.org/10.3189/172756500781832891, 2000. a
Petitet, A. and Dongarra, J.: Algorithmic redistribution methods for block-cyclic decompositions, IEEE T. Parall. Distr., 10, 1201–1216, https://doi.org/10.1109/71.819944, 1999. a
Quiquet, A., Dumas, C., Ritz, C., Peyaud, V., and Roche, D. M.: The GRISLI ice sheet model (version 2.0): calibration and validation for multi-millennial changes of the Antarctic ice sheet, Geosci. Model Dev., 11, 5003–5025, https://doi.org/10.5194/gmd-11-5003-2018, 2018. a
Rew, R. and Davis, G.: NetCDF: an interface for scientific data access, IEEE Comput. Graph., 10, 76–82, https://doi.org/10.1109/38.56302, 1990. a
Rignot, E., Mouginot, J., Scheuchl, B., van den Broeke, M., van Wessem, M. J., and Morlighem, M.: Four decades of Antarctic Ice Sheet mass balance from 1979–2017, P. Natl. Acad. Sci. USA, 116, 1095–1103, https://doi.org/10.1073/pnas.1812883116, 2019. 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
Robel, A. A., Schoof, C., and Tziperman, E.: Rapid grounding line migration induced by internal ice stream variability, J. Geophys. Res.-Earth, 119, 2430–2447, https://doi.org/10.1002/2014jf003251, 2014. a
Robel, A. A., Seroussi, H., and Roe, G. H.: Marine ice sheet instability amplifies and skews uncertainty in projections of future sea-level rise, P. Natl. Acad. Sci. USA, 116, 14887–14892, https://doi.org/10.1073/pnas.1904822116, 2019. a, b, c
Robinson, A., Goldberg, D., and Lipscomb, W. H.: A comparison of the stability and performance of depth-integrated ice-dynamics solvers, The Cryosphere, 16, 689–709, https://doi.org/10.5194/tc-16-689-2022, 2022. a
Schoof, C.: The effect of cavitation on glacier sliding, P. Roy. Soc. A-Math. Phy., 461, 609–627, https://doi.org/10.1098/rspa.2004.1350, 2005. a, b
Schoof, C.: A variational approach to ice stream flow, J. Fluid Mech., 556, 227, https://doi.org/10.1017/s0022112006009591, 2006a. a
Schoof, C.: Marine ice-sheet dynamics. Part 1. The case of rapid sliding, J. Fluid Mech., 573, 27–55, https://doi.org/10.1017/s0022112006003570, 2007b. a, b
Schoof, C.: Ice-sheet acceleration driven by melt supply variability, Nature, 468, 803–806, https://doi.org/10.1038/nature09618, 2010. a
Schoof, C.: Marine ice sheet dynamics. Part 2. A Stokes flow contact problem, J. Fluid Mech., 679, 122–155, https://doi.org/10.1017/jfm.2011.129, 2011. a, b
Schoof, C. and Hindmarsh, R. C. A.: Thin-Film Flows with Wall Slip: An Asymptotic Analysis of Higher Order Glacier Flow Models, Q. J. Mech. Appl. Math., 63, 73–114, https://doi.org/10.1093/qjmam/hbp025, 2010. a
Sergienko, O. V., Goldberg, D. N., and Little, C. M.: Alternative ice shelf equilibria determined by ocean environment, J. Geophys. Res.-Earth, 118, 970–981, https://doi.org/10.1002/jgrf.20054, 2013. a
Shepherd, A. and Wingham, D.: Recent Sea-Level Contributions of the Antarctic and Greenland Ice Sheets, Science, 315, 1529–1532, https://doi.org/10.1126/science.1136776, 2007. a
Shepherd, A., Fricker, H. A., and Farrell, S. L.: Trends and connections across the Antarctic cryosphere, Nature, 558, 223–232, https://doi.org/10.1038/s41586-018-0171-6, 2018. a, b
Shewchuk, J. R.: An Introduction to the Conjugate Gradient Method Without the Agonizing Pain, https://api.semanticscholar.org/CorpusID:6491967 (last access: 27 March 2025), 1994. a
Smedt, B. D., Pattyn, F., and Groen, P. D.: Using the unstable manifold correction in a Picard iteration to solve the velocity field in higher-order ice-flow models, J. Glaciol., 56, 257–261, https://doi.org/10.3189/002214310791968395, 2010. a
Stearns, L. A. and van der Veen, C. J.: Friction at the bed does not control fast glacier flow, Science, 361, 273–277, https://doi.org/10.1126/science.aat2217, 2018. a
Stone, P. H. and Carlson, J. H.: Atmospheric Lapse Rate Regimes and Their Parameterization, Journal of Atmospheric Sciences, 36, 415–423, https://doi.org/10.1175/1520-0469(1979)036<0415:ALRRAT>2.0.CO;2, 1979. a
Teschl, G.: Ordinary Differential Equations and Dynamical Systems, Graduate studies in mathematics, American Mathematical Society, https://books.google.es/books?id=FZ0CAQAAQBAJ (last access: 27 March 2025), 2012. a
Thomas, R. H. and Bentley, C. R.: A Model for Holocene Retreat of the West Antarctic Ice Sheet, Quaternary Res., 10, 150–170, https://doi.org/10.1016/0033-5894(78)90098-4, 1978. a, b
Truffer, M. and Fahnestock, M.: Rethinking Ice Sheet Time Scales, Science, 315, 1508–1510, https://doi.org/10.1126/science.1140469, 2007. a
Tulaczyk, S.: Ice sliding over weak, fine-grained tills: Dependence of ice-till interactions on till granulometry, in: Glacial Processes Past and Present, Geological Society of America, https://doi.org/10.1130/0-8137-2337-x.159, 1999. a
Tulaczyk, S., Kamb, B., Scherer, R. P., and Engelhardt, H. F.: Sedimentary processes at the base of a West Antarctic ice stream; constraints from textural and compositional properties of subglacial debris, J. Sediment. Res., 68, 487–496, https://doi.org/10.2110/jsr.68.487, 1998. 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, 2000. a
Van der Veen, C. J.: Tidewater calving, J. Glaciol., 42, 375–385, 1996. a
Van der Veen, C. J.: Calving glaciers, Prog. Phys. Geogr., 26, 96–122, https://doi.org/10.1191/0309133302pp327ra, 2002. a
van der Vorst, H. A.: Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems, SIAM J. Sci. Stat. Comp., 13, 631–644, https://doi.org/10.1137/0913035, 1992. a
van der Wel, N., Christoffersen, P., and Bougamont, M.: The influence of subglacial hydrology on the flow of Kamb Ice Stream, West Antarctica, J. Geophys. Res.-Earth, 118, 97–110, https://doi.org/10.1029/2012jf002570, 2013. a
Vaughan, D. G. and Arthern, R.: Why Is It Hard to Predict the Future of Ice Sheets?, Science, 315, 1503–1504, https://doi.org/10.1126/science.1141111, 2007. a
Veen, C. V. D. and Whillans, I.: Force Budget: I. Theory and Numerical Methods, J. Glaciol., 35, 53–60, https://doi.org/10.3189/002214389793701581, 1989. a
Weertman, J.: On the Sliding of Glaciers, J. Glaciol., 3, 33–38, https://doi.org/10.3189/s0022143000024709, 1957. a
Weertman, J.: Stability of the Junction of an Ice Sheet and an Ice Shelf, J. Glaciol., 13, 3–11, https://doi.org/10.3189/s0022143000023327, 1974. a, b, c, d
Whillans, I. M. and van der Veen, C. J.: The role of lateral drag in the dynamics of Ice Stream B, Antarctica, J. Glaciol., 43, 231–237, https://doi.org/10.3189/s0022143000003178, 1997. a
Wolf, M. M., Hendrickson, B. A., and Boman, E. G.: Optimizing parallel sparse matrix-vector multiplication by partitioning, Sandia National Laboratories (SNL), Albuquerque, NM, and Livermore, CA (United States), https://www.osti.gov/biblio/949000 (last access: 27 March 2025), 2008. a
Zoet, L. K. and Iverson, N. R.: A slip law for glaciers on deformable beds, Science, 368, 76–78, https://doi.org/10.1126/science.aaz1183, 2020. a, b
Short summary
We introduce Nix, an ice-sheet model designed for understanding how large masses of ice behave. Nix is a computer programme that simulates the movement and temperature evolution in ice sheets. It helps us study how ice sheets respond to changes in the atmosphere and ocean. We found that ice temperatures play an essential role in determining the motion and stability of ice sheets. Nix is a useful tool for learning how climate change affects polar ice sheets.
We introduce Nix, an ice-sheet model designed for understanding how large masses of ice behave....
