Articles | Volume 18, issue 13
https://doi.org/10.5194/gmd-18-4023-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-4023-2025
© Author(s) 2025. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Development and technical paper
| 
02 Jul 2025
Development and technical paper |
| 02 Jul 2025
| 02 Jul 2025
Anisotropic metric-based mesh adaptation for ice flow modelling in Firedrake
Department of Geosciences, University of Oslo, Oslo, Norway
Joseph G. Wallwork
Institute of Computing for Climate Science, University of Cambridge, Cambridge, UK
Stephan C. Kramer
Department of Earth Science and Engineering, Imperial College London, London, UK
Fabien Gillet-Chaulet
Univ. Grenoble Alpes, CNRS, INRAE, IRD, Grenoble INP, IGE, 38000 Grenoble, France
Regine Hock
Department of Geosciences, University of Oslo, Oslo, Norway
Geophysical Institute, University of Alaska Fairbanks, Fairbanks, USA
Matthew D. Piggott
Department of Earth Science and Engineering, Imperial College London, London, UK
Related authors
No articles found.
Justine Caillet, Nicolas C. Jourdain, Pierre Mathiot, Fabien Gillet-Chaulet, Benoit Urruty, Clara Burgard, Charles Amory, Mondher Chekki, and Christoph Kittel
Earth Syst. Dynam., 16, 293–315, https://doi.org/10.5194/esd-16-293-2025, https://doi.org/10.5194/esd-16-293-2025, 2025
Short summary
Short summary
Internal climate variability, resulting from processes intrinsic to the climate system, modulates the Antarctic response to climate change by delaying or offsetting its effects. Using climate and ice-sheet models, we highlight that irreducible internal climate variability significantly enlarges the likely range of Antarctic contribution to sea-level rise until 2100. Thus, we recommend considering internal climate variability as a source of uncertainty for future ice-sheet projections.
Eliot Jager, Fabien Gillet-Chaulet, Nicolas Champollion, Romain Millan, Heiko Goelzer, and Jérémie Mouginot
The Cryosphere, 18, 5519–5550, https://doi.org/10.5194/tc-18-5519-2024, https://doi.org/10.5194/tc-18-5519-2024, 2024
Short summary
Short summary
Inspired by a previous intercomparison framework, our study better constrains uncertainties in glacier evolution using an innovative method to validate Bayesian calibration. Upernavik Isstrøm, one of Greenland's largest glaciers, has lost significant mass since 1985. By integrating observational data, climate models, human emissions, and internal model parameters, we project its evolution until 2100. We show that future human emissions are the main source of uncertainty in 2100, making up half.
Sia Ghelichkhan, Angus Gibson, D. Rhodri Davies, Stephan C. Kramer, and David A. Ham
Geosci. Model Dev., 17, 5057–5086, https://doi.org/10.5194/gmd-17-5057-2024, https://doi.org/10.5194/gmd-17-5057-2024, 2024
Short summary
Short summary
We introduce the Geoscientific ADjoint Optimisation PlaTform (G-ADOPT), designed for inverse modelling of Earth system processes, with an initial focus on mantle dynamics. G-ADOPT is built upon Firedrake, Dolfin-Adjoint and the Rapid Optimisation Library, which work together to optimise models using an adjoint method, aligning them with seismic and geologic datasets. We demonstrate G-ADOPT's ability to reconstruct mantle evolution and thus be a powerful tool in geosciences.
Eric Petersen, Regine Hock, and Michael G. Loso
Earth Surf. Dynam., 12, 727–745, https://doi.org/10.5194/esurf-12-727-2024, https://doi.org/10.5194/esurf-12-727-2024, 2024
Short summary
Short summary
Ice cliffs are melt hot spots that increase melt rates on debris-covered glaciers which otherwise see a reduction in melt rates. In this study, we show how surface runoff streams contribute to the generation, evolution, and survival of ice cliffs by carving into the glacier and transporting rocky debris. On Kennicott Glacier, Alaska, 33 % of ice cliffs are actively influenced by streams, while nearly half are within 10 m of streams.
Shuaib Rasheed, Simon C. Warder, Yves Plancherel, and Matthew D. Piggott
Nat. Hazards Earth Syst. Sci., 24, 737–755, https://doi.org/10.5194/nhess-24-737-2024, https://doi.org/10.5194/nhess-24-737-2024, 2024
Short summary
Short summary
Here we use a high-resolution bathymetry dataset of the Maldives archipelago, as well as corresponding high numerical model resolution, to carry out a scenario-based tsunami hazard assessment for the entire Maldives archipelago to investigate the potential impact of plausible far-field tsunamis across the Indian Ocean at the island scale. The results indicate that several factors contribute to mitigating and amplifying tsunami waves at the island scale.
Baptiste Vandecrux, Robert S. Fausto, Jason E. Box, Federico Covi, Regine Hock, Åsa K. Rennermalm, Achim Heilig, Jakob Abermann, Dirk van As, Elisa Bjerre, Xavier Fettweis, Paul C. J. P. Smeets, Peter Kuipers Munneke, Michiel R. van den Broeke, Max Brils, Peter L. Langen, Ruth Mottram, and Andreas P. Ahlstrøm
The Cryosphere, 18, 609–631, https://doi.org/10.5194/tc-18-609-2024, https://doi.org/10.5194/tc-18-609-2024, 2024
Short summary
Short summary
How fast is the Greenland ice sheet warming? In this study, we compiled 4500+ temperature measurements at 10 m below the ice sheet surface (T10m) from 1912 to 2022. We trained a machine learning model on these data and reconstructed T10m for the ice sheet during 1950–2022. After a slight cooling during 1950–1985, the ice sheet warmed at a rate of 0.7 °C per decade until 2022. Climate models showed mixed results compared to our observations and underestimated the warming in key regions.
Emily A. Hill, Benoît Urruty, Ronja Reese, Julius Garbe, Olivier Gagliardini, Gaël Durand, Fabien Gillet-Chaulet, G. Hilmar Gudmundsson, Ricarda Winkelmann, Mondher Chekki, David Chandler, and Petra M. Langebroek
The Cryosphere, 17, 3739–3759, https://doi.org/10.5194/tc-17-3739-2023, https://doi.org/10.5194/tc-17-3739-2023, 2023
Short summary
Short summary
The grounding lines of the Antarctic Ice Sheet could enter phases of irreversible retreat or advance. We use three ice sheet models to show that the present-day locations of Antarctic grounding lines are reversible with respect to a small perturbation away from their current position. This indicates that present-day retreat of the grounding lines is not yet irreversible or self-enhancing.
Ronja Reese, Julius Garbe, Emily A. Hill, Benoît Urruty, Kaitlin A. Naughten, Olivier Gagliardini, Gaël Durand, Fabien Gillet-Chaulet, G. Hilmar Gudmundsson, David Chandler, Petra M. Langebroek, and Ricarda Winkelmann
The Cryosphere, 17, 3761–3783, https://doi.org/10.5194/tc-17-3761-2023, https://doi.org/10.5194/tc-17-3761-2023, 2023
Short summary
Short summary
We use an ice sheet model to test where current climate conditions in Antarctica might lead. We find that present-day ocean and atmosphere conditions might commit an irreversible collapse of parts of West Antarctica which evolves over centuries to millennia. Importantly, this collapse is not irreversible yet.
Jennifer R. Shadrick, Dylan H. Rood, Martin D. Hurst, Matthew D. Piggott, Klaus M. Wilcken, and Alexander J. Seal
Earth Surf. Dynam., 11, 429–450, https://doi.org/10.5194/esurf-11-429-2023, https://doi.org/10.5194/esurf-11-429-2023, 2023
Short summary
Short summary
This study uses a coastal evolution model to interpret cosmogenic beryllium-10 concentrations and topographic data and, in turn, quantify long-term cliff retreat rates for four chalk sites on the south coast of England. By using a process-based model, clear distinctions between intertidal weathering rates have been recognised between chalk and sandstone rock coast sites, advocating the use of process-based models to interpret the long-term behaviour of rock coasts.
Mariana C. A. Clare, Tim W. B. Leijnse, Robert T. McCall, Ferdinand L. M. Diermanse, Colin J. Cotter, and Matthew D. Piggott
Nat. Hazards Earth Syst. Sci., 22, 2491–2515, https://doi.org/10.5194/nhess-22-2491-2022, https://doi.org/10.5194/nhess-22-2491-2022, 2022
Short summary
Short summary
Assessing uncertainty is computationally expensive because it requires multiple runs of expensive models. We take the novel approach of assessing uncertainty from coastal flooding using a multilevel multifidelity (MLMF) method which combines the efficiency of less accurate models with the accuracy of more expensive models at different resolutions. This significantly reduces the computational cost but maintains accuracy, making previously unfeasible real-world studies possible.
D. Rhodri Davies, Stephan C. Kramer, Sia Ghelichkhan, and Angus Gibson
Geosci. Model Dev., 15, 5127–5166, https://doi.org/10.5194/gmd-15-5127-2022, https://doi.org/10.5194/gmd-15-5127-2022, 2022
Short summary
Short summary
Firedrake is a state-of-the-art system that automatically generates highly optimised code for simulating finite-element (FE) problems in geophysical fluid dynamics. It creates a separation of concerns between employing the FE method and implementing it. Here, we demonstrate the applicability and benefits of Firedrake for simulating geodynamical flows, with a focus on the slow creeping motion of Earth's mantle over geological timescales, which is ultimately the engine driving our dynamic Earth.
Anna Derkacheva, Fabien Gillet-Chaulet, Jeremie Mouginot, Eliot Jager, Nathan Maier, and Samuel Cook
The Cryosphere, 15, 5675–5704, https://doi.org/10.5194/tc-15-5675-2021, https://doi.org/10.5194/tc-15-5675-2021, 2021
Short summary
Short summary
Along the edges of the Greenland Ice Sheet surface melt lubricates the bed and causes large seasonal fluctuations in ice speeds during summer. Accurately understanding how these ice speed changes occur is difficult due to the inaccessibility of the glacier bed. We show that by using surface velocity maps with high temporal resolution and numerical modelling we can infer the basal conditions that control seasonal fluctuations in ice speed and gain insight into seasonal dynamics over large areas.
Jennifer R. Shadrick, Martin D. Hurst, Matthew D. Piggott, Bethany G. Hebditch, Alexander J. Seal, Klaus M. Wilcken, and Dylan H. Rood
Earth Surf. Dynam., 9, 1505–1529, https://doi.org/10.5194/esurf-9-1505-2021, https://doi.org/10.5194/esurf-9-1505-2021, 2021
Short summary
Short summary
Here we use topographic and 10Be concentration data to optimise a coastal evolution model. Cliff retreat rates are calculated for two UK sites for the past 8000 years and, for the first time, highlight a strong link between the rate of sea level rise and long-term cliff retreat rates. This method enables us to study past cliff response to sea level rise and so to greatly improve forecasts of future responses to accelerations in sea level rise that will result from climate change.
Stephan C. Kramer, D. Rhodri Davies, and Cian R. Wilson
Geosci. Model Dev., 14, 1899–1919, https://doi.org/10.5194/gmd-14-1899-2021, https://doi.org/10.5194/gmd-14-1899-2021, 2021
Short summary
Short summary
Computational models of Earth's mantle require rigorous verification and validation. Analytical solutions of the underlying Stokes equations provide a method to verify that these equations are accurately solved for. However, their derivation in spherical and cylindrical shell domains with physically relevant boundary conditions is involved. This paper provides a number of solutions. They are provided in a Python package (Assess) and their use is demonstrated in a convergence study with Fluidity.
Nathan Maier, Florent Gimbert, Fabien Gillet-Chaulet, and Adrien Gilbert
The Cryosphere, 15, 1435–1451, https://doi.org/10.5194/tc-15-1435-2021, https://doi.org/10.5194/tc-15-1435-2021, 2021
Short summary
Short summary
In Greenland, ice motion and the surface geometry depend on the friction at the bed. We use satellite measurements and modeling to determine how ice speeds and friction are related across the ice sheet. The relationships indicate that ice flowing over bed bumps sets the friction across most of the ice sheet's on-land regions. This result helps simplify and improve our understanding of how ice motion will change in the future.
Shuaib Rasheed, Simon C. Warder, Yves Plancherel, and Matthew D. Piggott
Ocean Sci., 17, 319–334, https://doi.org/10.5194/os-17-319-2021, https://doi.org/10.5194/os-17-319-2021, 2021
Short summary
Short summary
Environmental issues arising due to coastal modification and future sea level scenarios are a major environmental hazard facing the Maldives today. Here, we carry out high-resolution tidal modelling of a Maldivian atoll for the first time and show that coastal modification in the island scale is capable of driving large-scale change in the wider atoll basin in a short time, comparable to that of long-term sea level rise scenarios and on par with observations.
Vincent Peyaud, Coline Bouchayer, Olivier Gagliardini, Christian Vincent, Fabien Gillet-Chaulet, Delphine Six, and Olivier Laarman
The Cryosphere, 14, 3979–3994, https://doi.org/10.5194/tc-14-3979-2020, https://doi.org/10.5194/tc-14-3979-2020, 2020
Short summary
Short summary
Alpine glaciers are retreating at an accelerating rate in a warming climate. Numerical models allow us to study and anticipate these changes, but the performance of a model is difficult to evaluate. So we compared an ice flow model with the long dataset of observations obtained between 1979 and 2015 on Mer de Glace (Mont Blanc area). The model accurately reconstructs the past evolution of the glacier. We simulate the future evolution of Mer de Glace; it could retreat by 2 to 6 km by 2050.
Fabien Gillet-Chaulet
The Cryosphere, 14, 811–832, https://doi.org/10.5194/tc-14-811-2020, https://doi.org/10.5194/tc-14-811-2020, 2020
Short summary
Short summary
Marine-based sectors of the Antarctic Ice Sheet are increasingly contributing to sea-level rise. The basal conditions exert an important control on the ice dynamics. For obvious reasons of inaccessibility, they are an important source of uncertainties in numerical ice flow models used for sea-level projections. Here we assess the performance of an ensemble Kalman filter for the assimilation of transient observations of surface elevation and velocities in a marine ice sheet model.
Andrew Bliss, Regine Hock, Gabriel Wolken, Erin Whorton, Caroline Aubry-Wake, Juliana Braun, Alessio Gusmeroli, Will Harrison, Andrew Hoffman, Anna Liljedahl, and Jing Zhang
Earth Syst. Sci. Data, 12, 403–427, https://doi.org/10.5194/essd-12-403-2020, https://doi.org/10.5194/essd-12-403-2020, 2020
Short summary
Short summary
Extensive field observations were conducted in the Upper Susitna basin in central Alaska in 2012–2014. This paper describes the weather, glacier mass balance, snow cover, and soils of the basin. We found that temperatures over the glacier are cooler than over land at the same elevation. The glaciers have been losing mass faster in recent years than in the 1980s. Measurements of glacier mass change with traditional methods closely matched radar measurements.
Clemens Schannwell, Reinhard Drews, Todd A. Ehlers, Olaf Eisen, Christoph Mayer, and Fabien Gillet-Chaulet
The Cryosphere, 13, 2673–2691, https://doi.org/10.5194/tc-13-2673-2019, https://doi.org/10.5194/tc-13-2673-2019, 2019
Short summary
Short summary
Ice rises are important ice-sheet features that archive the ice sheet's history in their internal structure. Here we use a 3-D numerical ice-sheet model to simulate mechanisms that lead to changes in the geometry of the internal structure. We find that changes in snowfall result in much larger and faster changes than similar changes in ice-shelf geometry. This result is integral to fully unlocking the potential of ice rises as ice-dynamic archives and potential ice-core drilling sites.
Lionel Favier, Nicolas C. Jourdain, Adrian Jenkins, Nacho Merino, Gaël Durand, Olivier Gagliardini, Fabien Gillet-Chaulet, and Pierre Mathiot
Geosci. Model Dev., 12, 2255–2283, https://doi.org/10.5194/gmd-12-2255-2019, https://doi.org/10.5194/gmd-12-2255-2019, 2019
Short summary
Short summary
The melting at the base of floating ice shelves is the main driver of the Antarctic ice sheet current retreat. Here, we use an ideal set-up to assess a wide range of melting parameterisations depending on oceanic properties with regard to a new ocean–ice-sheet coupled model, published here for the first time. A parameterisation that depends quadratically on thermal forcing in both a local and a non-local way yields the best results and needs to be further assessed with more realistic set-ups.
Hélène Seroussi, Sophie Nowicki, Erika Simon, Ayako Abe-Ouchi, Torsten Albrecht, Julien Brondex, Stephen Cornford, Christophe Dumas, Fabien Gillet-Chaulet, Heiko Goelzer, Nicholas R. Golledge, Jonathan M. Gregory, Ralf Greve, Matthew J. Hoffman, Angelika Humbert, Philippe Huybrechts, Thomas Kleiner, Eric Larour, Gunter Leguy, William H. Lipscomb, Daniel Lowry, Matthias Mengel, Mathieu Morlighem, Frank Pattyn, Anthony J. Payne, David Pollard, Stephen F. Price, Aurélien Quiquet, Thomas J. Reerink, Ronja Reese, Christian B. Rodehacke, Nicole-Jeanne Schlegel, Andrew Shepherd, Sainan Sun, Johannes Sutter, Jonas Van Breedam, Roderik S. W. van de Wal, Ricarda Winkelmann, and Tong Zhang
The Cryosphere, 13, 1441–1471, https://doi.org/10.5194/tc-13-1441-2019, https://doi.org/10.5194/tc-13-1441-2019, 2019
Short summary
Short summary
We compare a wide range of Antarctic ice sheet simulations with varying initialization techniques and model parameters to understand the role they play on the projected evolution of this ice sheet under simple scenarios. Results are improved compared to previous assessments and show that continued improvements in the representation of the floating ice around Antarctica are critical to reduce the uncertainty in the future ice sheet contribution to sea level rise.
Ian Allison, Charles Fierz, Regine Hock, Andrew Mackintosh, Georg Kaser, and Samuel U. Nussbaumer
Hist. Geo Space. Sci., 10, 97–107, https://doi.org/10.5194/hgss-10-97-2019, https://doi.org/10.5194/hgss-10-97-2019, 2019
Short summary
Short summary
The International Association of Cryospheric Sciences (IACS) became the eighth and most recent association of IUGG in July 2007. IACS was launched in recognition of the importance of the cryosphere, particularly at a time of significant global change. The forbears of IACS, however, start with the 1894 Commission Internationale des Glaciers (CIG). This paper traces the transition from CIG to IACS; scientific objectives that drove activities and changes, and key events and individuals involved.
Julien Brondex, Fabien Gillet-Chaulet, and Olivier Gagliardini
The Cryosphere, 13, 177–195, https://doi.org/10.5194/tc-13-177-2019, https://doi.org/10.5194/tc-13-177-2019, 2019
Short summary
Short summary
Here, we apply a synthetic perturbation to the most active drainage basin of Antarctica and show that centennial mass loss projections obtained through ice flow models depend strongly on the implemented friction law, i.e. the mathematical relationship between basal drag and sliding velocities. In particular, the commonly used Weertman law considerably underestimates the sea-level contribution of this basin in comparison to two water pressure-dependent laws which rely on stronger physical bases.
Tuomas Kärnä, Stephan C. Kramer, Lawrence Mitchell, David A. Ham, Matthew D. Piggott, and António M. Baptista
Geosci. Model Dev., 11, 4359–4382, https://doi.org/10.5194/gmd-11-4359-2018, https://doi.org/10.5194/gmd-11-4359-2018, 2018
Short summary
Short summary
Unstructured meshes are attractive for coastal ocean modeling, as they allow more accurate representation of complex coastal topography. Unstructured mesh models are, however, often perceived as slow and inaccurate. We present a novel discontinuous Galerkin ocean model: Thetis. We demonstrate that the model is able to simulate baroclinic ocean flows with high accuracy on a triangular prismatic mesh. This work paves the way for highly accurate and efficient three-dimensional coastal ocean models.
Denis Cohen, Fabien Gillet-Chaulet, Wilfried Haeberli, Horst Machguth, and Urs H. Fischer
The Cryosphere, 12, 2515–2544, https://doi.org/10.5194/tc-12-2515-2018, https://doi.org/10.5194/tc-12-2515-2018, 2018
Short summary
Short summary
As part of an integrative study about the safety of repositories for radioactive waste under ice age conditions in Switzerland, we modeled the flow of ice of the Rhine glacier at the Last Glacial Maximum to determine conditions at the ice–bed interface. Results indicate that portions of the ice lobes were at the melting temperature and ice was sliding, two conditions necessary for erosion by glacier. Conditions at the bed of the ice lobes were affected by climate and also by topography.
Olivier Passalacqua, Marie Cavitte, Olivier Gagliardini, Fabien Gillet-Chaulet, Frédéric Parrenin, Catherine Ritz, and Duncan Young
The Cryosphere, 12, 2167–2174, https://doi.org/10.5194/tc-12-2167-2018, https://doi.org/10.5194/tc-12-2167-2018, 2018
Short summary
Short summary
Locating a suitable drill site is a key step in the Antarctic oldest-ice challenge. Here we have conducted a 3-D ice flow simulation in the region of Dome C using a refined bedrock description. Five selection criteria are computed that together provide an objective overview on the local ice flow conditions. We delineate kilometer-scale favorable areas that overlap with the ones recently proposed by another group. We propose a few drill sites that should be surveyed during the next field seasons.
Heiko Goelzer, Sophie Nowicki, Tamsin Edwards, Matthew Beckley, Ayako Abe-Ouchi, Andy Aschwanden, Reinhard Calov, Olivier Gagliardini, Fabien Gillet-Chaulet, Nicholas R. Golledge, Jonathan Gregory, Ralf Greve, Angelika Humbert, Philippe Huybrechts, Joseph H. Kennedy, Eric Larour, William H. Lipscomb, Sébastien Le clec'h, Victoria Lee, Mathieu Morlighem, Frank Pattyn, Antony J. Payne, Christian Rodehacke, Martin Rückamp, Fuyuki Saito, Nicole Schlegel, Helene Seroussi, Andrew Shepherd, Sainan Sun, Roderik van de Wal, and Florian A. Ziemen
The Cryosphere, 12, 1433–1460, https://doi.org/10.5194/tc-12-1433-2018, https://doi.org/10.5194/tc-12-1433-2018, 2018
Short summary
Short summary
We have compared a wide spectrum of different initialisation techniques used in the ice sheet modelling community to define the modelled present-day Greenland ice sheet state as a starting point for physically based future-sea-level-change projections. Compared to earlier community-wide comparisons, we find better agreement across different models, which implies overall improvement of our understanding of what is needed to produce such initial states.
Johannes Jakob Fürst, Fabien Gillet-Chaulet, Toby J. Benham, Julian A. Dowdeswell, Mariusz Grabiec, Francisco Navarro, Rickard Pettersson, Geir Moholdt, Christopher Nuth, Björn Sass, Kjetil Aas, Xavier Fettweis, Charlotte Lang, Thorsten Seehaus, and Matthias Braun
The Cryosphere, 11, 2003–2032, https://doi.org/10.5194/tc-11-2003-2017, https://doi.org/10.5194/tc-11-2003-2017, 2017
Short summary
Short summary
For the large majority of glaciers and ice caps, there is no information on the thickness of the ice cover. Any attempt to predict glacier demise under climatic warming and to estimate the future contribution to sea-level rise is limited as long as the glacier thickness is not well constrained. Here, we present a two-step mass-conservation approach for mapping ice thickness. Measurements are naturally reproduced. The reliability is readily assessible from a complementary map of error estimates.
Samuel D. Parkinson, Simon W. Funke, Jon Hill, Matthew D. Piggott, and Peter A. Allison
Geosci. Model Dev., 10, 1051–1068, https://doi.org/10.5194/gmd-10-1051-2017, https://doi.org/10.5194/gmd-10-1051-2017, 2017
Torbjørn Ims Østby, Thomas Vikhamar Schuler, Jon Ove Hagen, Regine Hock, Jack Kohler, and Carleen H. Reijmer
The Cryosphere, 11, 191–215, https://doi.org/10.5194/tc-11-191-2017, https://doi.org/10.5194/tc-11-191-2017, 2017
Short summary
Short summary
We present modelled climatic mass balance for all glaciers in Svalbard for the period 1957–2014 at 1 km resolution using a coupled surface energy balance and snowpack model, thereby closing temporal and spatial gaps in direct and geodetic mass balance estimates.
Supporting previous studies, our results indicate increased mass loss over the period.
A detailed analysis of the involved energy fluxes reveals that increased mass loss is caused by atmospheric warming further amplified by feedbacks.
O. Gagliardini, J. Brondex, F. Gillet-Chaulet, L. Tavard, V. Peyaud, and G. Durand
The Cryosphere, 10, 307–312, https://doi.org/10.5194/tc-10-307-2016, https://doi.org/10.5194/tc-10-307-2016, 2016
Short summary
Short summary
In this paper it is shown that the sensitivity to the mesh resolution is not
improved for a vanishing friction at the grounding line (GL). For a discontinuous friction at the GL, we further show that the results are moreover very sensitive to the way the friction is interpolated in the close vicinity of the GL. In the light of these new insights, new results for the MISMIP3d experiments obtained for higher resolutions than previously published are made available for future comparisons.
C. T. Jacobs and M. D. Piggott
Geosci. Model Dev., 8, 533–547, https://doi.org/10.5194/gmd-8-533-2015, https://doi.org/10.5194/gmd-8-533-2015, 2015
A. S. Candy, A. Avdis, J. Hill, G. J. Gorman, and M. D. Piggott
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmdd-7-5993-2014, https://doi.org/10.5194/gmdd-7-5993-2014, 2014
Revised manuscript has not been submitted
S. D. Parkinson, J. Hill, M. D. Piggott, and P. A. Allison
Geosci. Model Dev., 7, 1945–1960, https://doi.org/10.5194/gmd-7-1945-2014, https://doi.org/10.5194/gmd-7-1945-2014, 2014
J. Hill, E. E. Popova, D. A. Ham, M. D. Piggott, and M. Srokosz
Ocean Sci., 10, 323–343, https://doi.org/10.5194/os-10-323-2014, https://doi.org/10.5194/os-10-323-2014, 2014
T. L. Edwards, X. Fettweis, O. Gagliardini, F. Gillet-Chaulet, H. Goelzer, J. M. Gregory, M. Hoffman, P. Huybrechts, A. J. Payne, M. Perego, S. Price, A. Quiquet, and C. Ritz
The Cryosphere, 8, 181–194, https://doi.org/10.5194/tc-8-181-2014, https://doi.org/10.5194/tc-8-181-2014, 2014
T. L. Edwards, X. Fettweis, O. Gagliardini, F. Gillet-Chaulet, H. Goelzer, J. M. Gregory, M. Hoffman, P. Huybrechts, A. J. Payne, M. Perego, S. Price, A. Quiquet, and C. Ritz
The Cryosphere, 8, 195–208, https://doi.org/10.5194/tc-8-195-2014, https://doi.org/10.5194/tc-8-195-2014, 2014
F. Gillet-Chaulet, O. Gagliardini, H. Seddik, M. Nodet, G. Durand, C. Ritz, T. Zwinger, R. Greve, and D. G. Vaughan
The Cryosphere, 6, 1561–1576, https://doi.org/10.5194/tc-6-1561-2012, https://doi.org/10.5194/tc-6-1561-2012, 2012
Related subject area
Cryosphere
Description and validation of the ice-sheet model Nix v1.0
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)
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
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
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
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)
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
Daniel Moreno-Parada, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas
Geosci. Model Dev., 18, 3895–3919, https://doi.org/10.5194/gmd-18-3895-2025, https://doi.org/10.5194/gmd-18-3895-2025, 2025
Short summary
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.
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.
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.
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.
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.
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.
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.
Cited articles
Alauzet, F. and Loseille, A.: A decade of progress on anisotropic mesh adaptation for computational fluid dynamics, Comput.-Aid. Design, 72, 13–39, 2016. a
Alauzet, F. and Olivier, G.: Extension of metric-based anisotropic mesh adaptation to time-dependent problems involving moving geometries, in: 49th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, Orlando, Florida, 4–7 January 2011, p. 896, https://doi.org/10.2514/6.2011-896, 2011. a
Alnæs, M. S., Logg, A., Ølgaard, K. B., Rognes, M. E., and Wells, G. N.: Unified form language: A domain-specific language for weak formulations of partial differential equations, ACM Trans. Math. Softw., 40, 1–37, https://doi.org/10.1145/2566630, 2014. a
Amestoy, P. R., Duff, I. S., L'Excellent, J.-Y., and Koster, J.: A fully asynchronous multifrontal solver using distributed dynamic scheduling, SIAM J. Matr. Anal. Appl., 23, 15–41, 2001. a
Amestoy, P. R., Guermouche, A., L'Excellent, J.-Y., and Pralet, S.: Hybrid scheduling for the parallel solution of linear systems, Parallel Comput., 32, 136–156, 2006. a
Asay-Davis, X. S., Cornford, S. L., Durand, G., Galton-Fenzi, B. K., Gladstone, R. M., Gudmundsson, G. H., Hattermann, T., Holland, D. M., Holland, D., Holland, P. R., Martin, D. F., Mathiot, P., Pattyn, F., and Seroussi, H.: Experimental design for three interrelated marine ice sheet and ocean model intercomparison projects: MISMIP v.3 (MISMIP+), ISOMIP v.2 (ISOMIP+) and MISOMIP v.1 (MISOMIP1), Geosci. Model Dev., 9, 2471–2497, https://doi.org/10.5194/gmd-9-2471-2016, 2016. a, b
Aschwanden, A., Bartholomaus, T. C., Brinkerhoff, D. J., and Truffer, M.: Brief communication: A roadmap towards credible projections of ice sheet contribution to sea level, The Cryosphere, 15, 5705–5715, https://doi.org/10.5194/tc-15-5705-2021, 2021. a
Balarac, G., Basile, F., Bénard, P., Bordeu, F., Chapelier, J.-B., Cirrottola, L., Caumon, G., Dapogny, C., Frey, P., Froehly, A., Ghigliotti, G., Laraufie, R., Lartigue, G., Legentil, C., Mercier, R., Moureau, V., Nardoni, C., Pertant, S., and Zakari, M.: Tetrahedral remeshing in the context of large-scale numerical simulation and high performance computing, Math. Action, 11, 129–164, https://doi.org/10.5802/msia.22, 2022. a
Balay, S., Gropp, W. D., McInnes, L. C., and Smith, B. F.: Efficient Management of Parallelism in Object Oriented Numerical Software Libraries, in: Modern Software Tools in Scientific Computing, edited by: Arge, E., Bruaset, A. M., and Langtangen, H. P., Birkhäuser Press, 163–202, https://doi.org/10.1007/978-1-4612-1986-6_8, 1997. a
Balay, S., Abhyankar, S., Adams, M. F., Benson, S., Brown, J., Brune, P., Buschelman, K., Constantinescu, E., Dalcin, L., Dener, A., Eijkhout, V., Faibussowitsch, J., Gropp, W. D., Hapla, V., Isaac, T., Jolivet, P., Karpeev, D., Kaushik, D., Knepley, M. G., Kong, F., Kruger, S., May, D. A., McInnes, L. C., Mills, R. T., Mitchell, L., Munson, T., Roman, J. E., Rupp, K., Sanan, P., Sarich, J., Smith, B. F., Zampini, S., Zhang, H., Zhang, H., and Zhang, J.: PETSc/TAO Users Manual, Tech. Rep. ANL-21/39 – Revision 3.20, Argonne National Laboratory, https://doi.org/10.2172/1968587, 2023. a
Borouchaki, H., Hecht, F., and Frey, P. J.: Mesh gradation control, Int. J. Numer. Meth. Eng., 43, 1143–1165, 1998. a
Ciarlet, P.: Basic error estimates for elliptic problems, in: Finite Element Methods (Part 1), vol. 2 of Handbook of Numerical Analysis, Elsevier, 17–351, https://doi.org/10.1016/S1570-8659(05)80039-0, 1991. a, b
Ciarlet, P. G.: The finite element method for elliptic problems, Classics in Applied Mathematics, SIAM, https://doi.org/10.1137/1.9780898719208, 2002. a
Cornford, S. L., Seroussi, H., Asay-Davis, X. S., Gudmundsson, G. H., Arthern, R., Borstad, C., Christmann, J., Dias dos Santos, T., Feldmann, J., Goldberg, D., Hoffman, M. J., Humbert, A., Kleiner, T., Leguy, G., Lipscomb, W. H., Merino, N., Durand, G., Morlighem, M., Pollard, D., Rückamp, M., Williams, C. R., and Yu, H.: Results of the third Marine Ice Sheet Model Intercomparison Project (MISMIP+), The Cryosphere, 14, 2283–2301, https://doi.org/10.5194/tc-14-2283-2020, 2020. a, b, c, d
Crameri, F.: Scientific colour maps, Zenodo [data set], https://doi.org/10.5281/zenodo.8409685, 2023. a
Cuzzone, J. K., Morlighem, M., Larour, E., Schlegel, N., and Seroussi, H.: Implementation of higher-order vertical finite elements in ISSM v4.13 for improved ice sheet flow modeling over paleoclimate timescales, Geosci. Model Dev., 11, 1683–1694, https://doi.org/10.5194/gmd-11-1683-2018, 2018. a
Dapogny, C., Dobrzynski, C., and Frey, P.: Three-dimensional adaptive domain remeshing, implicit domain meshing, and applications to free and moving boundary problems, J. Comput. Phys., 262, 358–378, https://doi.org/10.1016/j.jcp.2014.01.005, 2014a. a
Dapogny, C., Dobrzynski, C., and Frey, P.: Three-dimensional adaptive domain remeshing, implicit domain meshing, and applications to free and moving boundary problems, J. Comput. Phys., 262, 358–378, 2014b. a
Davies, D. R., Wilson, C. R., and Kramer, S. C.: Fluidity: A fully unstructured anisotropic adaptive mesh computational modeling framework for geodynamics, Geochem. Geophy. Geosy., 12, Q06001, https://doi.org/10.1029/2011GC003551, 2011. a
de Diego, G. G., Farrell, P. E., and Hewitt, I. J.: Numerical approximation of viscous contact problems applied to glacial sliding, J. Fluid Mech., 938, A21, https://doi.org/10.1017/jfm.2022.178, 2022. a
Dobrzynski, C. and Frey, P.: Anisotropic Delaunay Mesh Adaptation for Unsteady Simulations, in: Proceedings of the 17th International Meshing Roundtable, edited by: Garimella, R. V., Springer, Berlin, Heidelberg, 177–194, ISBN 978-3-540-87921-3, 2008. a
dos Santos, T. D., Morlighem, M., Seroussi, H., Devloo, P. R. B., and Simões, J. C.: Implementation and performance of adaptive mesh refinement in the Ice Sheet System Model (ISSM v4.14), Geosci. Model Dev., 12, 215–232, https://doi.org/10.5194/gmd-12-215-2019, 2019. a, b, c
Dundovic, D.: ddundo/mismip-adapt-paper: Code for Anisotropic metric-based mesh adaptation for ice flow modelling in Firedrake paper, Zenodo [code], https://doi.org/10.5281/zenodo.14913405, 2025. a
Durand, G. and Pattyn, F.: Reducing uncertainties in projections of Antarctic ice mass loss, The Cryosphere, 9, 2043–2055, https://doi.org/10.5194/tc-9-2043-2015, 2015. a
Farrell, P. and Maddison, J.: Conservative interpolation between volume meshes by local Galerkin projection, Comput. Meth. Appl. Mech. Eng., 200, 89–100, 2011. a
Farrell, P. E.: Galerkin projection of discrete fields via supermesh construction, PhD thesis, Imperial College, London, https://doi.org/10.25560/5515, 2010. a
Gagliardini, O., Brondex, J., Gillet-Chaulet, F., Tavard, L., Peyaud, V., and Durand, G.: Brief communication: Impact of mesh resolution for MISMIP and MISMIP3d experiments using Elmer/Ice, The Cryosphere, 10, 307–312, https://doi.org/10.5194/tc-10-307-2016, 2016. a
George, P. L., Hecht, F., and Vallet, M. G.: Creation of internal points in Voronoi's type method. Control adaptation, Adv. Eng. Softw. Workstations, 13, 303–312, 1991. a
Geuzaine, C. and Remacle, J.-F.: Gmsh: A 3-D finite element mesh generator with built-in pre- and post-processing facilities, Int. J. Numer. Meth. Eng., 79, 1309–1331, 2009. a
Gillet-Chaulet, F., Gagliardini, O., Seddik, H., Nodet, M., Durand, G., Ritz, C., Zwinger, T., Greve, R., and Vaughan, D. G.: Greenland ice sheet contribution to sea-level rise from a new-generation ice-sheet model, The Cryosphere, 6, 1561–1576, https://doi.org/10.5194/tc-6-1561-2012, 2012. a
Gladstone, R., Lee, V., Vieli, A., and Payne, A.: Grounding line migration in an adaptive mesh ice sheet model, J. Geophys. Res.-Earth, 115, F04014, https://doi.org/10.1029/2009JF001615, 2010. a
Gladstone, R. M., Payne, A. J., and Cornford, S. L.: Resolution requirements for grounding-line modelling: sensitivity to basal drag and ice-shelf buttressing, Ann. Glaciol., 53, 97–105, 2012. a
Goldberg, D., Holland, D., and Schoof, C.: Grounding line movement and ice shelf buttressing in marine ice sheets, J. Geophys. Res.-Earth, 114, F04026, https://doi.org/10.1029/2008JF001227, 2009. a, b, c, d
Greve, R. and Blatter, H.: Dynamics of ice sheets and glaciers, Springer Science & Business Media, https://doi.org/10.1007/978-3-642-03415-2, 2009. a, b
Gudmundsson, G. H., Krug, J., Durand, G., Favier, L., and Gagliardini, O.: The stability of grounding lines on retrograde slopes, The Cryosphere, 6, 1497–1505, https://doi.org/10.5194/tc-6-1497-2012, 2012. a, b
Ham, D. A., Kelly, P. H. J., Mitchell, L., Cotter, C. J., Kirby, R. C., Sagiyama, K., Bouziani, N., Vorderwuelbecke, S., Gregory, T. J., Betteridge, J., Shapero, D. R., Nixon-Hill, R. W., Ward, C. J., Farrell, P. E., Brubeck, P. D., Marsden, I., Gibson, T. H., Homolya, M., Sun, T., McRae, A. T. T., Luporini, F., Gregory, A., Lange, M., Funke, S. W., Rathgeber, F., Bercea, G.-T., and Markall, G. R.: Firedrake User Manual, in: 1st Edn., Imperial College London and University of Oxford and Baylor University and University of Washington, https://doi.org/10.25561/104839, 2023. a
Hunter, J. D.: Matplotlib: A 2D graphics environment, Comput. Sci. Eng., 9, 90–95, https://doi.org/10.1109/MCSE.2007.55, 2007. a
Kirby, R .C. and Shapero, D.: High-order bounds-satisfying approximation of partial differential equations via finite element variational inequalities, Numer. Math. 156, 927–947https://doi.org/10.1007/s00211-024-01405-y, 2024. a
Lange, M., Mitchell, L., Knepley, M. G., and Gorman, G. J.: Efficient mesh management in Firedrake using PETSc-DMPlex, SIAM J. Sci. Comput., 38, S143–S155, https://doi.org/10.1137/15M1026092, 2016. a
Leguy, G. R., Asay-Davis, X. S., and Lipscomb, W. H.: Parameterization of basal friction near grounding lines in a one-dimensional ice sheet model, The Cryosphere, 8, 1239–1259, https://doi.org/10.5194/tc-8-1239-2014, 2014. a, b
Löhner, R. and Baum, J. D.: Adaptive h-refinement on 3D unstructured grids for transient problems, Int. J. Numer. Meth. Fluids, 14, 1407–1419, 1992. a
Loseille, A. and Alauzet, F.: Continuous mesh model and well-posed continuous interpolation error estimation, PhD thesis, INRIA, https://inria.hal.science/inria-00370235v1 (last access: 25 June 2025), 2009. a
McRae, A. T. T., Cotter, C. J., and Budd, C. J.: Optimal-Transport-Based Mesh Adaptivity on the Plane and Sphere Using Finite Elements, SIAM J. Sci. Comput., 40, A1121–A1148, https://doi.org/10.1137/16M1109515, 2018. a
Moreno-Parada, D., Robinson, A., Montoya, M., and Alvarez-Solas, J.: Description and validation of the ice sheet model Nix v1.0, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2023-2690, 2023. a
Morlighem, M., Rignot, E., Seroussi, H., Larour, E., Ben Dhia, H., and Aubry, D.: Spatial patterns of basal drag inferred using control methods from a full-Stokes and simpler models for Pine Island Glacier, West Antarctica, Geophys. Res. Lett., 37, L14502, https://doi.org/10.1029/2010GL043853, 2010. a
Pain, C., Umpleby, A., de Oliveira, C., and Goddard, A.: Tetrahedral mesh optimisation and adaptivity for steady-state and transient finite element calculations, Comput. Meth. Appl. Mech. Eng., 190, 3771–3796, https://doi.org/10.1016/S0045-7825(00)00294-2, 2001. a
Pattyn, F., Huyghe, A., De Brabander, S., and De Smedt, B.: Role of transition zones in marine ice sheet dynamics, J. Geophys. Res.-Earth, 111, F02004, https://doi.org/10.1029/2005JF000394, 2006. 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
attyn, F., Perichon, L., Durand, G., Favier, L., Gagliardini, O., Hindmarsh, R. C., Zwinger, T., Albrecht, T., Cornford, S., Docquier, D., Fürst, J. J., Goldberg, D., Gudmundsson, G. H., Humbert, A., Hütten, M., Huybrechts, P., Jouvet, G., Kleiner, T., Larour, E., Martin, D., Morlighem, M., Payne, A. J., Pollard, D., Rückamp, M., Rybak, O., Seroussi, H., Thoma, M., and Wilkens, N.: Grounding-line migration in plan-view marine ice-sheet models: results of the ice2sea MISMIP3d intercomparison, J. Glaciol., 59, 410–422, 2013. a
Peraire, J., Vahdati, M., Morgan, K., and Zienkiewicz, O. C.: Adaptive remeshing for compressible flow computations, J. Comput. Phys., 72, 449–466, 1987. a
Piggott, M., Farrell, P., Wilson, C., Gorman, G., and Pain, C.: Anisotropic mesh adaptivity for multi-scale ocean modelling, Philos. T. Roy. Soc. A, 367, 4591–4611, 2009a. a
Piggott, M. D., Farrell, P. E., Wilson, C. R., Gorman, G. J., and Pain, C. C.: Anisotropic mesh adaptivity for multi-scale ocean modelling, Philos. T. Roy. Soc. A, 367, 4591–4611, https://doi.org/10.1098/rsta.2009.0155, 2009b. a
Roache, P. J.: Code Verification by the Method of Manufactured Solutions, J. Fluids Eng., 124, 4–10, https://doi.org/10.1115/1.1436090, 2001. a
Schöberl, J.: NETGEN An advancing front 2D/3D-mesh generator based on abstract rules, Comput. Visual. Sci., 1, 41–52, 1997. a
Schoof, C.: Ice sheet grounding line dynamics: Steady states, stability, and hysteresis, J. Geophys. Res.-Earth, 112, F03S28, https://doi.org/10.1029/2006JF000664, 2007. a, b
Seddik, H., Greve, R., Zwinger, T., Gillet-Chaulet, F., and Gagliardini, O.: Simulations of the Greenland ice sheet 100 years into the future with the full Stokes model Elmer/Ice, J. Glaciol., 58, 427–440, 2012. a
Seroussi, H., Morlighem, M., Larour, E., Rignot, E., and Khazendar, A.: Hydrostatic grounding line parameterization in ice sheet models, The Cryosphere, 8, 2075–2087, https://doi.org/10.5194/tc-8-2075-2014, 2014. a
Smith, R. C., Hill, J., Collins, G. S., Piggott, M. D., Kramer, S. C., Parkinson, S. D., and Wilson, C.: Comparing approaches for numerical modelling of tsunami generation by deformable submarine slides, Ocean Model., 100, 125–140, 2016. a
Stokes, C. R.: Geomorphology under ice streams: moving from form to process, Earth Surf. Proc. Land., 43, 85–123, 2018. a
Sun, S., Pattyn, F., Simon, E. G., Albrecht, T., Cornford, S., Calov, R., Dumas, C., Gillet-Chaulet, F., Goelzer, H., Golledge, N. R., Greve, R., Hoffman, M. J., Humbert, A., Kazmierczak, E., Kleiner, T., Leguy, G. R., Lipscomb, W. H., Martin, D., Morlighem, M., Nowicki, S., Pollard, D., Price, S., Quiquet, A., Seroussi, H., Schlemm, T., Sutter, J., van de Wal, R. S. W., Winkelmann, R., and Zhang, T.: Antarctic ice sheet response to sudden and sustained ice-shelf collapse (ABUMIP), J. Glaciol., 66, 891–904, 2020. a
Vieli, A. and Payne, A. J.: Assessing the ability of numerical ice sheet models to simulate grounding line migration, J. Geophys. Res.-Earth, 110, F01003, https://doi.org/10.1029/2004JF000202, 2005. a, b, c
Wallwork, J. G., Knepley, M. G., Barral, N., and Piggott, M. D.: Parallel Metric-Based Mesh Adaptation in PETSc using ParMmg, 30th International Meshing Roundtable, arXiv [preprint], https://doi.org/10.48550/arXiv.2201.02806, 2022. a
Wallwork, J. G., Dundovic, D., Johnson, E., and Kramer, S. C.: Animate, Zenodo [code], https://doi.org/10.5281/zenodo.11537707, 2024a. a, b
Wallwork, J. G., Dundovic, D., Johnson, E., and Kramer, S. C.: Goalie, Zenodo [code], https://doi.org/10.5281/zenodo.13353191, 2024b. a, b
Wilson, C.: Modelling multiple-material flows on adaptive unstructured meshes, PhD thesis, Imperial College, London, https://doi.org/10.25560/5526, 2010. a
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.
Accurate numerical studies of glaciers often require high-resolution simulations, which often...
