Articles | Volume 18, issue 8
https://doi.org/10.5194/tc-18-3453-2024
© Author(s) 2024. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/tc-18-3453-2024
© Author(s) 2024. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Research article
| 
07 Aug 2024
Research article |
| 07 Aug 2024
| 07 Aug 2024
Increasing numerical stability of mountain valley glacier simulations: implementation and testing of free-surface stabilization in Elmer/Ice
Department of Mathematics, Stockholm University, Stockholm, Sweden
Thomas Zwinger
CSC – IT Center for Science Ltd., Espoo, Finland
Peter Råback
CSC – IT Center for Science Ltd., Espoo, Finland
Christian Helanow
Department of Mathematics, Stockholm University, Stockholm, Sweden
Swedish e-Science Research Centre, Stockholm, Sweden
Josefin Ahlkrona
Department of Mathematics, Stockholm University, Stockholm, Sweden
Swedish e-Science Research Centre, Stockholm, Sweden
Related authors
No articles found.
Yu Wang, Chen Zhao, Rupert Gladstone, Thomas Zwinger, Benjamin K. Galton-Fenzi, and Poul Christoffersen
The Cryosphere, 18, 5117–5137, https://doi.org/10.5194/tc-18-5117-2024, https://doi.org/10.5194/tc-18-5117-2024, 2024
Short summary
Short summary
Our research delves into the future evolution of Antarctica's Wilkes Subglacial Basin (WSB) and its potential contribution to sea level rise, focusing on how basal melt is implemented at the grounding line in ice flow models. Our findings suggest that these implementation methods can significantly impact the magnitude of future ice loss projections. Under a high-emission scenario, the WSB ice sheet could undergo massive and rapid retreat between 2200 and 2300.
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.
Yiliang Ma, Liyun Zhao, Rupert Gladstone, Thomas Zwinger, Michael Wolovick, and John C. Moore
EGUsphere, https://doi.org/10.5194/egusphere-2024-1102, https://doi.org/10.5194/egusphere-2024-1102, 2024
Short summary
Short summary
Totten Glacier in Antarctica holds a sea level potential of 3.85 m. Basal sliding and sub-shelf melt rate have important impact on ice sheet dynamics. We simulate the evolution of Totten Glacier using an ice flow model with different basal sliding parameterizations as well as sub-shelf melt rates to quantify their effect on the projections. We found the modelled glacier retreat and mass loss is sensitive to the choice of basal sliding parameterizations and maximal sub-shelf melt rate.
Hélène Seroussi, Vincent Verjans, Sophie Nowicki, Antony J. Payne, Heiko Goelzer, William H. Lipscomb, Ayako Abe-Ouchi, Cécile Agosta, Torsten Albrecht, Xylar Asay-Davis, Alice Barthel, Reinhard Calov, Richard Cullather, Christophe Dumas, Benjamin K. Galton-Fenzi, Rupert Gladstone, Nicholas R. Golledge, Jonathan M. Gregory, Ralf Greve, Tore Hattermann, Matthew J. Hoffman, Angelika Humbert, Philippe Huybrechts, Nicolas C. Jourdain, Thomas Kleiner, Eric Larour, Gunter R. Leguy, Daniel P. Lowry, Chistopher M. Little, Mathieu Morlighem, Frank Pattyn, Tyler Pelle, Stephen F. Price, Aurélien Quiquet, Ronja Reese, Nicole-Jeanne Schlegel, Andrew Shepherd, Erika Simon, Robin S. Smith, Fiammetta Straneo, Sainan Sun, Luke D. Trusel, Jonas Van Breedam, Peter Van Katwyk, Roderik S. W. van de Wal, Ricarda Winkelmann, Chen Zhao, Tong Zhang, and Thomas Zwinger
The Cryosphere, 17, 5197–5217, https://doi.org/10.5194/tc-17-5197-2023, https://doi.org/10.5194/tc-17-5197-2023, 2023
Short summary
Short summary
Mass loss from Antarctica is a key contributor to sea level rise over the 21st century, and the associated uncertainty dominates sea level projections. We highlight here the Antarctic glaciers showing the largest changes and quantify the main sources of uncertainty in their future evolution using an ensemble of ice flow models. We show that on top of Pine Island and Thwaites glaciers, Totten and Moscow University glaciers show rapid changes and a strong sensitivity to warmer ocean conditions.
Denis Cohen, Guillaume Jouvet, Thomas Zwinger, Angela Landgraf, and Urs H. Fischer
E&G Quaternary Sci. J., 72, 189–201, https://doi.org/10.5194/egqsj-72-189-2023, https://doi.org/10.5194/egqsj-72-189-2023, 2023
Short summary
Short summary
During glacial times in Switzerland, glaciers of the Alps excavated valleys in low-lying regions that were later filled with sediment or water. How glaciers eroded these valleys is not well understood because erosion occurred near ice margins where ice moved slowly and was present for short times. Erosion is linked to the speed of ice and to water flowing under it. Here we present a model that estimates the location of water channels beneath the ice and links these locations to zones of erosion.
Maryam Zarrinderakht, Christian Schoof, and Thomas Zwinger
EGUsphere, https://doi.org/10.5194/egusphere-2023-807, https://doi.org/10.5194/egusphere-2023-807, 2023
Preprint withdrawn
Short summary
Short summary
We used a model to study how crevasses propagate in ice shelves. Our model combines a viscous model and a fracture mechanics model. We studied periodic crevasses on an ice shelf being stretched. We show that existing models based only on stress cannot fully explain how crevasses grow and lead to iceberg calving. This model can be a useful tool to train a low-dimensional representation calving law for an ice sheet model.
Douglas I. Benn, Adrian Luckman, Jan A. Åström, Anna J. Crawford, Stephen L. Cornford, Suzanne L. Bevan, Thomas Zwinger, Rupert Gladstone, Karen Alley, Erin Pettit, and Jeremy Bassis
The Cryosphere, 16, 2545–2564, https://doi.org/10.5194/tc-16-2545-2022, https://doi.org/10.5194/tc-16-2545-2022, 2022
Short summary
Short summary
Thwaites Glacier (TG), in West Antarctica, is potentially unstable and may contribute significantly to sea-level rise as global warming continues. Using satellite data, we show that Thwaites Eastern Ice Shelf, the largest remaining floating extension of TG, has started to accelerate as it fragments along a shear zone. Computer modelling does not indicate that fragmentation will lead to imminent glacier collapse, but it is clear that major, rapid, and unpredictable changes are underway.
Rupert Gladstone, Benjamin Galton-Fenzi, David Gwyther, Qin Zhou, Tore Hattermann, Chen Zhao, Lenneke Jong, Yuwei Xia, Xiaoran Guo, Konstantinos Petrakopoulos, Thomas Zwinger, Daniel Shapero, and John Moore
Geosci. Model Dev., 14, 889–905, https://doi.org/10.5194/gmd-14-889-2021, https://doi.org/10.5194/gmd-14-889-2021, 2021
Short summary
Short summary
Retreat of the Antarctic ice sheet, and hence its contribution to sea level rise, is highly sensitive to melting of its floating ice shelves. This melt is caused by warm ocean currents coming into contact with the ice. Computer models used for future ice sheet projections are not able to realistically evolve these melt rates. We describe a new coupling framework to enable ice sheet and ocean computer models to interact, allowing projection of the evolution of melt and its impact on sea level.
Hélène Seroussi, Sophie Nowicki, Antony J. Payne, Heiko Goelzer, William H. Lipscomb, Ayako Abe-Ouchi, Cécile Agosta, Torsten Albrecht, Xylar Asay-Davis, Alice Barthel, Reinhard Calov, Richard Cullather, Christophe Dumas, Benjamin K. Galton-Fenzi, Rupert Gladstone, Nicholas R. Golledge, Jonathan M. Gregory, Ralf Greve, Tore Hattermann, Matthew J. Hoffman, Angelika Humbert, Philippe Huybrechts, Nicolas C. Jourdain, Thomas Kleiner, Eric Larour, Gunter R. Leguy, Daniel P. Lowry, Chistopher M. Little, Mathieu Morlighem, Frank Pattyn, Tyler Pelle, Stephen F. Price, Aurélien Quiquet, Ronja Reese, Nicole-Jeanne Schlegel, Andrew Shepherd, Erika Simon, Robin S. Smith, Fiammetta Straneo, Sainan Sun, Luke D. Trusel, Jonas Van Breedam, Roderik S. W. van de Wal, Ricarda Winkelmann, Chen Zhao, Tong Zhang, and Thomas Zwinger
The Cryosphere, 14, 3033–3070, https://doi.org/10.5194/tc-14-3033-2020, https://doi.org/10.5194/tc-14-3033-2020, 2020
Short summary
Short summary
The Antarctic ice sheet has been losing mass over at least the past 3 decades in response to changes in atmospheric and oceanic conditions. This study presents an ensemble of model simulations of the Antarctic evolution over the 2015–2100 period based on various ice sheet models, climate forcings and emission scenarios. Results suggest that the West Antarctic ice sheet will continue losing a large amount of ice, while the East Antarctic ice sheet could experience increased snow accumulation.
Thomas Zwinger, Grace A. Nield, Juha Ruokolainen, and Matt A. King
Geosci. Model Dev., 13, 1155–1164, https://doi.org/10.5194/gmd-13-1155-2020, https://doi.org/10.5194/gmd-13-1155-2020, 2020
Short summary
Short summary
We present a newly developed flat-earth model, Elmer/Earth, for viscoelastic treatment of solid earth deformation under ice loads. Unlike many previous approaches with proprietary software, this model is based on the open-source FEM code Elmer, with the advantage for scientists to apply and alter the model without license constraints. The new-generation full-stress ice-sheet model Elmer/Ice shares the same code base, enabling future coupled ice-sheet–glacial-isostatic-adjustment simulations.
Shahbaz Memon, Dorothée Vallot, Thomas Zwinger, Jan Åström, Helmut Neukirchen, Morris Riedel, and Matthias Book
Geosci. Model Dev., 12, 3001–3015, https://doi.org/10.5194/gmd-12-3001-2019, https://doi.org/10.5194/gmd-12-3001-2019, 2019
Short summary
Short summary
Scientific workflows enable complex scientific computational scenarios, which include data intensive scenarios, parametric executions, and interactive simulations. In this article, we applied the UNICORE workflow management system to automate a formerly hard-coded coupling of a glacier flow model and a calving model, which contain many tasks and dependencies, ranging from pre-processing and data management to repetitive executions on heterogeneous high-performance computing (HPC) resources.
Joe Todd, Poul Christoffersen, Thomas Zwinger, Peter Råback, and Douglas I. Benn
The Cryosphere, 13, 1681–1694, https://doi.org/10.5194/tc-13-1681-2019, https://doi.org/10.5194/tc-13-1681-2019, 2019
Short summary
Short summary
The Greenland Ice Sheet loses 30 %–60 % of its ice due to iceberg calving. Calving processes and their links to climate are not well understood or incorporated into numerical models of glaciers. Here we use a new 3-D calving model to investigate calving at Store Glacier, West Greenland, and test its sensitivity to increased submarine melting and reduced support from ice mélange (sea ice and icebergs). We find Store remains fairly stable despite these changes, but less so in the southern side.
Eef C. H. van Dongen, Nina Kirchner, Martin B. van Gijzen, Roderik S. W. van de Wal, Thomas Zwinger, Gong Cheng, Per Lötstedt, and Lina von Sydow
Geosci. Model Dev., 11, 4563–4576, https://doi.org/10.5194/gmd-11-4563-2018, https://doi.org/10.5194/gmd-11-4563-2018, 2018
Short summary
Short summary
Ice flow forced by gravity is governed by the full Stokes (FS) equations, which are computationally expensive to solve. Therefore, approximations to the FS equations are used, especially when modeling an ice sheet on long time spans. Here, we report a combination of an approximation with the FS equations that allows simulating the dynamics of ice sheets over long time spans without introducing artifacts caused by application of approximations in parts of the domain where they are not valid.
Chen Zhao, Rupert M. Gladstone, Roland C. Warner, Matt A. King, Thomas Zwinger, and Mathieu Morlighem
The Cryosphere, 12, 2637–2652, https://doi.org/10.5194/tc-12-2637-2018, https://doi.org/10.5194/tc-12-2637-2018, 2018
Short summary
Short summary
A combination of computer modelling and observational data were used to infer the resistance to ice flow at the bed of the Fleming Glacier on the Antarctic Peninsula. The model was also used to simulate the distribution of temperature within the ice, which governs the rate at which the ice can deform. This is especially important for glaciers like the Fleming Glacier, which has both regions of rapid deformation and regions of rapid sliding at the bed.
Chen Zhao, Rupert M. Gladstone, Roland C. Warner, Matt A. King, Thomas Zwinger, and Mathieu Morlighem
The Cryosphere, 12, 2653–2666, https://doi.org/10.5194/tc-12-2653-2018, https://doi.org/10.5194/tc-12-2653-2018, 2018
Short summary
Short summary
A combination of computer modelling and observational data were used to infer the resistance to ice flow at the bed of the Fleming Glacier on the Antarctic Peninsula in both 2008 and 2015. The comparison suggests the grounding line retreated by ~ 9 km from 2008 to 2015. The retreat may be enhanced by a positive feedback between friction, melting and sliding at the glacier bed.
Sue Cook, Jan Åström, Thomas Zwinger, Benjamin Keith Galton-Fenzi, Jamin Stevens Greenbaum, and Richard Coleman
The Cryosphere, 12, 2401–2411, https://doi.org/10.5194/tc-12-2401-2018, https://doi.org/10.5194/tc-12-2401-2018, 2018
Short summary
Short summary
The growth of fractures on Antarctic ice shelves is important because it controls the amount of ice lost as icebergs. We use a model constructed of multiple interconnected blocks to predict the locations where fractures will form on the Totten Ice Shelf in East Antarctica. The results show that iceberg calving is controlled not only by fractures forming near the front of the ice shelf but also by fractures which formed many kilometres upstream.
Yongmei Gong, Thomas Zwinger, Jan Åström, Bas Altena, Thomas Schellenberger, Rupert Gladstone, and John C. Moore
The Cryosphere, 12, 1563–1577, https://doi.org/10.5194/tc-12-1563-2018, https://doi.org/10.5194/tc-12-1563-2018, 2018
Short summary
Short summary
In this study we apply a discrete element model capable of simulating ice fracturing. A microscopic-scale discrete process is applied in addition to a continuum ice dynamics model to investigate the mechanisms facilitated by basal meltwater production, surface meltwater and ice crack opening, for the surge in Basin 3, Austfonna ice cap. The discrete element model is used to locate the ice cracks that can penetrate though the full thickness of the glacier and deliver surface water to the bed.
Dorothée Vallot, Jan Åström, Thomas Zwinger, Rickard Pettersson, Alistair Everett, Douglas I. Benn, Adrian Luckman, Ward J. J. van Pelt, Faezeh Nick, and Jack Kohler
The Cryosphere, 12, 609–625, https://doi.org/10.5194/tc-12-609-2018, https://doi.org/10.5194/tc-12-609-2018, 2018
Short summary
Short summary
This paper presents a new perspective on the role of ice dynamics and ocean interaction in glacier calving processes applied to Kronebreen, a tidewater glacier in Svalbard. A global modelling approach includes ice flow modelling, undercutting estimation by a combination of glacier energy balance and plume modelling as well as calving by a discrete particle model. We show that modelling undercutting is necessary and calving is influenced by basal friction velocity and geometry.
Hakime Seddik, Ralf Greve, Thomas Zwinger, and Shin Sugiyama
The Cryosphere, 11, 2213–2229, https://doi.org/10.5194/tc-11-2213-2017, https://doi.org/10.5194/tc-11-2213-2017, 2017
Short summary
Short summary
The Shirase Glacier in Antarctica is studied by means of a computer model. This model implements two physical approaches to represent the glacier flow dynamics. This study finds that it is important to use the more precise and sophisticated method in order to better understand and predict the evolution of fast flowing glaciers. This may be important to more accurately predict the sea level change due to global warming.
Rupert Michael Gladstone, Roland Charles Warner, Benjamin Keith Galton-Fenzi, Olivier Gagliardini, Thomas Zwinger, and Ralf Greve
The Cryosphere, 11, 319–329, https://doi.org/10.5194/tc-11-319-2017, https://doi.org/10.5194/tc-11-319-2017, 2017
Short summary
Short summary
Computer models are used to simulate the behaviour of glaciers and ice sheets. It has been found that such models are required to be run at very high resolution (which means high computational expense) in order to accurately represent the evolution of marine ice sheets (ice sheets resting on bedrock below sea level), in certain situations which depend on sub-glacial physical processes.
T. Zwinger, T. Malm, M. Schäfer, R. Stenberg, and J. C. Moore
The Cryosphere, 9, 1415–1426, https://doi.org/10.5194/tc-9-1415-2015, https://doi.org/10.5194/tc-9-1415-2015, 2015
Short summary
Short summary
By deploying a large-scale high-resolution turbulent CFD simulation using the present-day topography of the Scharffenbergbotnen (SBB) valley, we show how the surrounding topography redirects incoming easterly katabatic storm fronts to impact the blue ice areas (BIA) inside the valley, where the snow cover frequently is removed. A further simulation of a reconstructed topography at the Late Glacial Maximum further reveals that the BIA at SBB must have formed after this period.
M. Schäfer, F. Gillet-Chaulet, R. Gladstone, R. Pettersson, V. A. Pohjola, T. Strozzi, and T. Zwinger
The Cryosphere, 8, 1951–1973, https://doi.org/10.5194/tc-8-1951-2014, https://doi.org/10.5194/tc-8-1951-2014, 2014
R. Gladstone, M. Schäfer, T. Zwinger, Y. Gong, T. Strozzi, R. Mottram, F. Boberg, and J. C. Moore
The Cryosphere, 8, 1393–1405, https://doi.org/10.5194/tc-8-1393-2014, https://doi.org/10.5194/tc-8-1393-2014, 2014
B. Sun, J. C. Moore, T. Zwinger, L. Zhao, D. Steinhage, X. Tang, D. Zhang, X. Cui, and C. Martín
The Cryosphere, 8, 1121–1128, https://doi.org/10.5194/tc-8-1121-2014, https://doi.org/10.5194/tc-8-1121-2014, 2014
S. Cook, I. C. Rutt, T. Murray, A. Luckman, T. Zwinger, N. Selmes, A. Goldsack, and T. D. James
The Cryosphere, 8, 827–841, https://doi.org/10.5194/tc-8-827-2014, https://doi.org/10.5194/tc-8-827-2014, 2014
T. Zwinger, M. Schäfer, C. Martín, and J. C. Moore
The Cryosphere, 8, 607–621, https://doi.org/10.5194/tc-8-607-2014, https://doi.org/10.5194/tc-8-607-2014, 2014
T. Sato, T. Shiraiwa, R. Greve, H. Seddik, E. Edelmann, and T. Zwinger
Clim. Past, 10, 393–404, https://doi.org/10.5194/cp-10-393-2014, https://doi.org/10.5194/cp-10-393-2014, 2014
B. de Fleurian, O. Gagliardini, T. Zwinger, G. Durand, E. Le Meur, D. Mair, and P. Råback
The Cryosphere, 8, 137–153, https://doi.org/10.5194/tc-8-137-2014, https://doi.org/10.5194/tc-8-137-2014, 2014
J. A. Åström, T. I. Riikilä, T. Tallinen, T. Zwinger, D. Benn, J. C. Moore, and J. Timonen
The Cryosphere, 7, 1591–1602, https://doi.org/10.5194/tc-7-1591-2013, https://doi.org/10.5194/tc-7-1591-2013, 2013
O. Gagliardini, T. Zwinger, F. Gillet-Chaulet, G. Durand, L. Favier, B. de Fleurian, R. Greve, M. Malinen, C. Martín, P. Råback, J. Ruokolainen, M. Sacchettini, M. Schäfer, H. Seddik, and J. Thies
Geosci. Model Dev., 6, 1299–1318, https://doi.org/10.5194/gmd-6-1299-2013, https://doi.org/10.5194/gmd-6-1299-2013, 2013
A. S. Drouet, D. Docquier, G. Durand, R. Hindmarsh, F. Pattyn, O. Gagliardini, and T. Zwinger
The Cryosphere, 7, 395–406, https://doi.org/10.5194/tc-7-395-2013, https://doi.org/10.5194/tc-7-395-2013, 2013
L. Zhao, L. Tian, T. Zwinger, R. Ding, J. Zong, Q. Ye, and J. C. Moore
The Cryosphere Discuss., https://doi.org/10.5194/tcd-7-145-2013, https://doi.org/10.5194/tcd-7-145-2013, 2013
Revised manuscript not accepted
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
Discipline: Glaciers | Subject: Numerical Modelling
Inter-model differences in 21st century glacier runoff for the world's major river basins
A minimal machine-learning glacier mass balance model
Physically based modelling of glacier evolution under climate change in the tropical Andes
New glacier thickness and bed topography maps for Svalbard
Sea-level rise contribution from Ryder Glacier in Northern Greenland varies by an order of magnitude by 2300 depending on future emissions
Quantifying the buttressing contribution of landfast sea ice and melange to Crane Glacier, Antarctic Peninsula
Application of a regularised Coulomb sliding law to Jakobshavn Isbræ, western Greenland
The demise of the world's largest piedmont glacier: a probabilistic forecast
A 3D glacier dynamics–line plume model to estimate the frontal ablation of Hansbreen, Svalbard
Impact of the Nares Strait sea ice arches on the long-term stability of the Petermann Glacier ice shelf
Reconciling ice dynamics and bed topography with a versatile and fast ice thickness inversion
Exploring the ability of the variable-resolution Community Earth System Model to simulate cryospheric–hydrological variables in High Mountain Asia
Modelling the development and decay of cryoconite holes in northwestern Greenland
Thermal regime of the Grigoriev ice cap and the Sary-Tor glacier in the inner Tien Shan, Kyrgyzstan
Modelling supraglacial debris-cover evolution from the single-glacier to the regional scale: an application to High Mountain Asia
The 21st-century fate of the Mocho-Choshuenco ice cap in southern Chile
Modelling steady states and the transient response of debris-covered glaciers
Twentieth century global glacier mass change: an ensemble-based model reconstruction
Mapping the age of ice of Gauligletscher combining surface radionuclide contamination and ice flow modeling
Modelling the evolution of Djankuat Glacier, North Caucasus, from 1752 until 2100 CE
Brief communication: Time step dependence (and fixes) in Stokes simulations of calving ice shelves
Modelling regional glacier length changes over the last millennium using the Open Global Glacier Model
The contrasting response of outlet glaciers to interior and ocean forcing
Deep learning applied to glacier evolution modelling
Initialization of a global glacier model based on present-day glacier geometry and past climate information: an ensemble approach
Contrasting thinning patterns between lake- and land-terminating glaciers in the Bhutanese Himalaya
Impact of frontal ablation on the ice thickness estimation of marine-terminating glaciers in Alaska
Modeling the response of Greenland outlet glaciers to global warming using a coupled flow line–plume model
Buoyant forces promote tidewater glacier iceberg calving through large basal stress concentrations
Global glacier volume projections under high-end climate change scenarios
Finn Wimberly, Lizz Ultee, Lilian Schuster, Matthias Huss, David R. Rounce, Fabien Maussion, Sloan Coats, Jonathan Mackay, and Erik Holmgren
The Cryosphere, 19, 1491–1511, https://doi.org/10.5194/tc-19-1491-2025, https://doi.org/10.5194/tc-19-1491-2025, 2025
Short summary
Short summary
Glacier models have historically been used to understand glacier melt’s contribution to sea level rise. The capacity to project seasonal glacier runoff is a relatively recent development for these models. In this study we provide the first model intercomparison of runoff projections for the glacier evolution models capable of simulating future runoff globally. We compare model projections from 2000 to 2100 for all major river basins larger than 3000 km2 with over 30 km2 of initial glacier cover.
Marijn van der Meer, Harry Zekollari, Matthias Huss, Jordi Bolibar, Kamilla Hauknes Sjursen, and Daniel Farinotti
The Cryosphere, 19, 805–826, https://doi.org/10.5194/tc-19-805-2025, https://doi.org/10.5194/tc-19-805-2025, 2025
Short summary
Short summary
Glacier retreat poses big challenges, making understanding how climate affects glaciers vital. But glacier measurements worldwide are limited. We created a simple machine-learning model called miniML-MB, which estimates annual changes in glacier mass in the Swiss Alps. As input, miniML-MB uses two climate variables: average temperature (May–Aug) and total precipitation (Oct–Feb). Our model can accurately predict glacier mass from 1961 to 2021 but struggles for extreme years (2022 and 2023).
Jonathan D. Mackay, Nicholas E. Barrand, David M. Hannah, Emily Potter, Nilton Montoya, and Wouter Buytaert
The Cryosphere, 19, 685–712, https://doi.org/10.5194/tc-19-685-2025, https://doi.org/10.5194/tc-19-685-2025, 2025
Short summary
Short summary
We combine two globally capable glacier evolution models to include processes that are typically neglected but thought to control tropical glacier retreat (e.g. sublimation). We apply the model to Peru's Vilcanota-Urubamba Basin. The model captures observed glacier mass changes,but struggles with surface albedo dynamics. Projections show glacier mass shrinking to 17 % or 6 % of 2000 levels by 2100 under moderate- and high-emission scenarios, respectively.
Ward van Pelt and Thomas Frank
The Cryosphere, 19, 1–17, https://doi.org/10.5194/tc-19-1-2025, https://doi.org/10.5194/tc-19-1-2025, 2025
Short summary
Short summary
Accurate information on the ice thickness of Svalbard's glaciers is important for assessing the contribution to sea level rise in a present and a future climate. However, direct observations of the glacier bed are scarce. Here, we use an inverse approach and high-resolution surface observations to infer basal conditions. We present and analyse the new bed and thickness maps, quantify the ice volume (6800 km3), and compare these against radar data and previous studies.
Felicity Alice Holmes, Jamie Barnett, Henning Åkesson, Mathieu Morlighem, Johan Nilsson, Nina Kirchner, and Martin Jakobsson
EGUsphere, https://doi.org/10.5194/egusphere-2024-3839, https://doi.org/10.5194/egusphere-2024-3839, 2024
Short summary
Short summary
Northern Greenland contains some of the ice sheet's last remaining glaciers with floating ice tongues. One of these is Ryder Glacier, which has been relatively stable in recent decades in contrast to nearby glaciers. Here, we use a computer model to simulate Ryder Glacier until 2300 under both a low and high emissions scenario. Very high levels of surface melt under a high emissions future leads to a sea-level rise contribution an order of magnitude higher than under a low emissions future.
Richard Parsons, Sainan Sun, G. Hilmar Gudmundsson, Jan Wuite, and Thomas Nagler
The Cryosphere, 18, 5789–5801, https://doi.org/10.5194/tc-18-5789-2024, https://doi.org/10.5194/tc-18-5789-2024, 2024
Short summary
Short summary
In 2022, multi-year landfast sea ice in Antarctica's Larsen B embayment disintegrated, after which time an increase in the rate at which Crane Glacier discharged ice into the ocean was observed. As the fast ice was joined to the glacier terminus, it could provide resistance against the glacier's flow, slowing down the rate of ice discharge. We used numerical modelling to quantify this resistive stress and found that the fast ice provided significant support to Crane prior to its disintegration.
Matt Trevers, Antony J. Payne, and Stephen L. Cornford
The Cryosphere, 18, 5101–5115, https://doi.org/10.5194/tc-18-5101-2024, https://doi.org/10.5194/tc-18-5101-2024, 2024
Short summary
Short summary
The form of the friction law which determines the speed of ice sliding over the bedrock remains a major source of uncertainty in ice sheet model projections of future sea level rise. Jakobshavn Isbræ, the fastest-flowing glacier in Greenland, which has undergone significant changes in the last few decades, is an ideal case for testing sliding laws. We find that a regularised Coulomb friction law reproduces the large seasonal and inter-annual flow speed variations most accurately.
Douglas Brinkerhoff, Brandon Tober, Michael Daniel, Victor Devaux-Chupin, Michael Christoffersen, John W. Holt, Christopher F. Larsen, Mark Fahnestock, Michael G. Loso, Kristin M. F. Timm, Russell Mitchell, and Martin Truffer
EGUsphere, https://doi.org/10.5194/egusphere-2024-2354, https://doi.org/10.5194/egusphere-2024-2354, 2024
Short summary
Short summary
Sít' Tlein is one of the largest glaciers in the world outside of the polar regions, and we know that it has been rapidly thinning. To forecast how this glacier will change in the future, we combine a computer model of ice flow with measurements from many different sources. Our model tells us that with high probability, Sít' Tlein's lower reaches are going to disappear in the next century and a half, creating a new bay or lake along Alaska's coastline.
José M. Muñoz-Hermosilla, Jaime Otero, Eva De Andrés, Kaian Shahateet, Francisco Navarro, and Iván Pérez-Doña
The Cryosphere, 18, 1911–1924, https://doi.org/10.5194/tc-18-1911-2024, https://doi.org/10.5194/tc-18-1911-2024, 2024
Short summary
Short summary
A large fraction of the mass loss from marine-terminating glaciers is attributed to frontal ablation. In this study, we used a 3D ice flow model of a real glacier that includes the effects of calving and submarine melting. Over a 30-month simulation, we found that the model reproduced the seasonal cycle for this glacier. Besides, the front positions were in good agreement with observations in the central part of the front, with longitudinal differences, on average, below 15 m.
Abhay Prakash, Qin Zhou, Tore Hattermann, and Nina Kirchner
The Cryosphere, 17, 5255–5281, https://doi.org/10.5194/tc-17-5255-2023, https://doi.org/10.5194/tc-17-5255-2023, 2023
Short summary
Short summary
Sea ice arch formation in the Nares Strait has shielded the Petermann Glacier ice shelf from enhanced basal melting. However, with the sustained decline of the Arctic sea ice predicted to continue, the ice shelf is likely to be exposed to a year-round mobile and thin sea ice cover. In such a scenario, our modelled results show that elevated temperatures, and more importantly, a stronger ocean circulation in the ice shelf cavity, could result in up to two-thirds increase in basal melt.
Thomas Frank, Ward J. J. van Pelt, and Jack Kohler
The Cryosphere, 17, 4021–4045, https://doi.org/10.5194/tc-17-4021-2023, https://doi.org/10.5194/tc-17-4021-2023, 2023
Short summary
Short summary
Since the ice thickness of most glaciers worldwide is unknown, and since it is not feasible to visit every glacier and observe their thickness directly, inverse modelling techniques are needed that can calculate ice thickness from abundant surface observations. Here, we present a new method for doing that. Our methodology relies on modelling the rate of surface elevation change for a given glacier, compare this with observations of the same quantity and change the bed until the two are in line.
René R. Wijngaard, Adam R. Herrington, William H. Lipscomb, Gunter R. Leguy, and Soon-Il An
The Cryosphere, 17, 3803–3828, https://doi.org/10.5194/tc-17-3803-2023, https://doi.org/10.5194/tc-17-3803-2023, 2023
Short summary
Short summary
We evaluate the ability of the Community Earth System Model (CESM2) to simulate cryospheric–hydrological variables, such as glacier surface mass balance (SMB), over High Mountain Asia (HMA) by using a global grid (~111 km) with regional refinement (~7 km) over HMA. Evaluations of two different simulations show that climatological biases are reduced, and glacier SMB is improved (but still too negative) by modifying the snow and glacier model and using an updated glacier cover dataset.
Yukihiko Onuma, Koji Fujita, Nozomu Takeuchi, Masashi Niwano, and Teruo Aoki
The Cryosphere, 17, 3309–3328, https://doi.org/10.5194/tc-17-3309-2023, https://doi.org/10.5194/tc-17-3309-2023, 2023
Short summary
Short summary
We established a novel model that simulates the temporal changes in cryoconite hole (CH) depth using heat budgets calculated independently at the ice surface and CH bottom based on hole shape geometry. The simulations suggest that CH depth is governed by the balance between the intensity of the diffuse component of downward shortwave radiation and the wind speed. The meteorological conditions may be important factors contributing to the recent ice surface darkening via the redistribution of CHs.
Lander Van Tricht and Philippe Huybrechts
The Cryosphere, 16, 4513–4535, https://doi.org/10.5194/tc-16-4513-2022, https://doi.org/10.5194/tc-16-4513-2022, 2022
Short summary
Short summary
We examine the thermal regime of the Grigoriev ice cap and the Sary-Tor glacier, both located in the inner Tien Shan in Kyrgyzstan. Our findings are important as the ice dynamics can only be understood and modelled precisely if ice temperature is considered correctly in ice flow models. The calibrated parameters of this study can be used in applications with ice flow models for individual ice masses as well as to optimise more general models for large-scale regional simulations.
Loris Compagno, Matthias Huss, Evan Stewart Miles, Michael James McCarthy, Harry Zekollari, Amaury Dehecq, Francesca Pellicciotti, and Daniel Farinotti
The Cryosphere, 16, 1697–1718, https://doi.org/10.5194/tc-16-1697-2022, https://doi.org/10.5194/tc-16-1697-2022, 2022
Short summary
Short summary
We present a new approach for modelling debris area and thickness evolution. We implement the module into a combined mass-balance ice-flow model, and we apply it using different climate scenarios to project the future evolution of all glaciers in High Mountain Asia. We show that glacier geometry, volume, and flow velocity evolve differently when modelling explicitly debris cover compared to glacier evolution without the debris-cover module, demonstrating the importance of accounting for debris.
Matthias Scheiter, Marius Schaefer, Eduardo Flández, Deniz Bozkurt, and Ralf Greve
The Cryosphere, 15, 3637–3654, https://doi.org/10.5194/tc-15-3637-2021, https://doi.org/10.5194/tc-15-3637-2021, 2021
Short summary
Short summary
We simulate the current state and future evolution of the Mocho-Choshuenco ice cap in southern Chile (40°S, 72°W) with the ice-sheet model SICOPOLIS. Under different global warming scenarios, we project ice mass losses between 56 % and 97 % by the end of the 21st century. We quantify the uncertainties based on an ensemble of climate models and on the temperature dependence of the equilibrium line altitude. Our results suggest a considerable deglaciation in southern Chile in the next 80 years.
James C. Ferguson and Andreas Vieli
The Cryosphere, 15, 3377–3399, https://doi.org/10.5194/tc-15-3377-2021, https://doi.org/10.5194/tc-15-3377-2021, 2021
Short summary
Short summary
Debris-covered glaciers have a greater extent than their debris-free counterparts due to insulation from the debris cover. However, the transient response to climate change remains poorly understood. We use a numerical model that couples ice dynamics and debris transport and varies the climate signal. We find that debris cover delays the transient response, especially for the extent. However, adding cryokarst features near the terminus greatly enhances the response.
Jan-Hendrik Malles and Ben Marzeion
The Cryosphere, 15, 3135–3157, https://doi.org/10.5194/tc-15-3135-2021, https://doi.org/10.5194/tc-15-3135-2021, 2021
Short summary
Short summary
To better estimate the uncertainty in glacier mass change modeling during the 20th century we ran an established model with an ensemble of meteorological data sets. We find that the total ensemble uncertainty, especially in the early 20th century, when glaciological and meteorological observations at glacier locations were sparse, increases considerably compared to individual ensemble runs. This stems from regions with a lot of ice mass but few observations (e.g., Greenland periphery).
Guillaume Jouvet, Stefan Röllin, Hans Sahli, José Corcho, Lars Gnägi, Loris Compagno, Dominik Sidler, Margit Schwikowski, Andreas Bauder, and Martin Funk
The Cryosphere, 14, 4233–4251, https://doi.org/10.5194/tc-14-4233-2020, https://doi.org/10.5194/tc-14-4233-2020, 2020
Short summary
Short summary
We show that plutonium is an effective tracer to identify ice originating from the early 1960s at the surface of a mountain glacier after a long time within the ice flow, giving unique information on the long-term former ice motion. Combined with ice flow modelling, the dating can be extended to the entire glacier, and we show that an airplane which crash-landed on the Gauligletscher in 1946 will likely soon be released from the ice close to the place where pieces have emerged in recent years.
Yoni Verhaegen, Philippe Huybrechts, Oleg Rybak, and Victor V. Popovnin
The Cryosphere, 14, 4039–4061, https://doi.org/10.5194/tc-14-4039-2020, https://doi.org/10.5194/tc-14-4039-2020, 2020
Short summary
Short summary
We use a numerical flow model to simulate the behaviour of the Djankuat Glacier, a WGMS reference glacier situated in the North Caucasus (Republic of Kabardino-Balkaria, Russian Federation), in response to past, present and future climate conditions (1752–2100 CE). In particular, we adapt a more sophisticated and physically based debris model, which has not been previously applied in time-dependent numerical flow line models, to look at the impact of a debris cover on the glacier’s evolution.
Brandon Berg and Jeremy Bassis
The Cryosphere, 14, 3209–3213, https://doi.org/10.5194/tc-14-3209-2020, https://doi.org/10.5194/tc-14-3209-2020, 2020
Short summary
Short summary
Computer models of ice sheets and glaciers are an important component of projecting sea level rise due to climate change. For models that seek to simulate the full balance of forces within the ice, if portions of the glacier are allowed to quickly break off in a process called iceberg calving, a numerical issue arises that can cause inaccurate results. We examine the issue and propose a solution so that future models can more accurately predict the future behavior of ice sheets and glaciers.
David Parkes and Hugues Goosse
The Cryosphere, 14, 3135–3153, https://doi.org/10.5194/tc-14-3135-2020, https://doi.org/10.5194/tc-14-3135-2020, 2020
Short summary
Short summary
Direct records of glacier changes rarely go back more than the last 100 years and are few and far between. We used a sophisticated glacier model to simulate glacier length changes over the last 1000 years for those glaciers that we do have long-term records of, to determine whether the model can run in a stable, realistic way over a long timescale, reproducing recent observed trends. We find that post-industrial changes are larger than other changes in this time period driven by recent warming.
John Erich Christian, Alexander A. Robel, Cristian Proistosescu, Gerard Roe, Michelle Koutnik, and Knut Christianson
The Cryosphere, 14, 2515–2535, https://doi.org/10.5194/tc-14-2515-2020, https://doi.org/10.5194/tc-14-2515-2020, 2020
Short summary
Short summary
We use simple, physics-based models to compare how marine-terminating glaciers respond to changes at their marine margin vs. inland surface melt. Initial glacier retreat is more rapid for ocean changes than for inland changes, but in both cases, glaciers will continue responding for millennia. We analyze several implications of these differing pathways of change. In particular, natural ocean variability must be better understood to correctly identify the anthropogenic role in glacier retreat.
Jordi Bolibar, Antoine Rabatel, Isabelle Gouttevin, Clovis Galiez, Thomas Condom, and Eric Sauquet
The Cryosphere, 14, 565–584, https://doi.org/10.5194/tc-14-565-2020, https://doi.org/10.5194/tc-14-565-2020, 2020
Short summary
Short summary
We introduce a novel approach for simulating glacier mass balances using a deep artificial neural network (i.e. deep learning) from climate and topographical data. This has been added as a component of a new open-source parameterized glacier evolution model. Deep learning is found to outperform linear machine learning methods, mainly due to its nonlinearity. Potential applications range from regional mass balance reconstructions from observations to simulations for past and future climates.
Julia Eis, Fabien Maussion, and Ben Marzeion
The Cryosphere, 13, 3317–3335, https://doi.org/10.5194/tc-13-3317-2019, https://doi.org/10.5194/tc-13-3317-2019, 2019
Short summary
Short summary
To provide estimates of past glacier mass changes, an adequate initial state is required. However, information about past glacier states at regional or global scales is largely incomplete. Our study presents a new way to initialize the Open Global Glacier Model from past climate information and present-day geometries. We show that even with perfectly known but incomplete boundary conditions, the problem of model initialization leads to nonunique solutions, and we propose an ensemble approach.
Shun Tsutaki, Koji Fujita, Takayuki Nuimura, Akiko Sakai, Shin Sugiyama, Jiro Komori, and Phuntsho Tshering
The Cryosphere, 13, 2733–2750, https://doi.org/10.5194/tc-13-2733-2019, https://doi.org/10.5194/tc-13-2733-2019, 2019
Short summary
Short summary
We investigate thickness change of Bhutanese glaciers during 2004–2011 using repeat GPS surveys and satellite-based observations. The thinning rate of Lugge Glacier (LG) is > 3 times that of Thorthormi Glacier (TG). Numerical simulations of ice dynamics and surface mass balance (SMB) demonstrate that the rapid thinning of LG is driven by both negative SMB and dynamic thinning, while the thinning of TG is minimised by a longitudinally compressive flow regime.
Beatriz Recinos, Fabien Maussion, Timo Rothenpieler, and Ben Marzeion
The Cryosphere, 13, 2657–2672, https://doi.org/10.5194/tc-13-2657-2019, https://doi.org/10.5194/tc-13-2657-2019, 2019
Short summary
Short summary
We have implemented a frontal ablation parameterization into the Open Global Glacier Model and have shown that inversion methods based on mass conservation systematically underestimate the mass turnover (and therefore the thickness) of tidewater glaciers when neglecting frontal ablation. This underestimation can rise up to 19 % on a regional scale. Not accounting for frontal ablation will have an impact on the estimate of the glaciers’ potential contribution to sea level rise.
Johanna Beckmann, Mahé Perrette, Sebastian Beyer, Reinhard Calov, Matteo Willeit, and Andrey Ganopolski
The Cryosphere, 13, 2281–2301, https://doi.org/10.5194/tc-13-2281-2019, https://doi.org/10.5194/tc-13-2281-2019, 2019
Short summary
Short summary
Submarine melting (SM) has been discussed as potentially triggering the recently observed retreat at outlet glaciers in Greenland. How much it may contribute in terms of future sea level rise (SLR) has not been quantified yet. When accounting for SM in our experiments, SLR contribution of 12 outlet glaciers increases by over 3-fold until the year 2100 under RCP8.5. Scaling up from 12 to all of Greenland's outlet glaciers increases future SLR contribution of Greenland by 50 %.
Matt Trevers, Antony J. Payne, Stephen L. Cornford, and Twila Moon
The Cryosphere, 13, 1877–1887, https://doi.org/10.5194/tc-13-1877-2019, https://doi.org/10.5194/tc-13-1877-2019, 2019
Short summary
Short summary
Iceberg calving is a major factor in the retreat of outlet glaciers of the Greenland Ice Sheet. Massive block overturning calving events occur at major outlet glaciers. A major calving event in 2009 was triggered by the release of a smaller block of ice from above the waterline. Using a numerical model, we investigate the feasibility of this mechanism to drive large calving events. We find that relatively small perturbations induce forces large enough to open cracks in ice at the glacier bed.
Sarah Shannon, Robin Smith, Andy Wiltshire, Tony Payne, Matthias Huss, Richard Betts, John Caesar, Aris Koutroulis, Darren Jones, and Stephan Harrison
The Cryosphere, 13, 325–350, https://doi.org/10.5194/tc-13-325-2019, https://doi.org/10.5194/tc-13-325-2019, 2019
Short summary
Short summary
We present global glacier volume projections for the end of this century, under a range of high-end climate change scenarios, defined as exceeding 2 °C global average warming. The ice loss contribution to sea level rise for all glaciers excluding those on the peripheral of the Antarctic ice sheet is 215.2 ± 21.3 mm. Such large ice losses will have consequences for sea level rise and for water supply in glacier-fed river systems.
Cited articles
Andrés-Martínez, M., Morgan, J. P., Pérez-Gussinyé, M., and Rüpke, L.: A new free-surface stabilization algorithm for geodynamical modelling: Theory and numerical tests, Phys. Earth Planet. In., 246, 41–51, https://doi.org/10.1016/j.pepi.2015.07.003, 2015. a
Babuška, I.: Error-bounds for finite element method, Numer. Math., 16, 322–333, https://doi.org/10.1007/bf02165003, 1971. a
Baiocchi, C., Brezzi, F., and Franca, L. P.: Virtual bubbles and Galerkin-least-squares type methods (Ga.L.S), Comput. Method. Appl. M., 105, 125–141, https://doi.org/10.1016/0045-7825(93)90119-I, 1993. 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.1017/S002214300001621X, 1995. a
Blatter, H., Greve, R., and Abe-Ouchi, A.: A short history of the thermomechanical theory and modelling of glaciers and ice sheets, J. Glaciol., 56, 1087–1094, https://doi.org/10.3189/002214311796406059, 2010. a
Brezzi, F.: On the existence, uniqueness and approximation of saddle-point problems arising from Lagrangian multipliers, Anal. Numer., 8, 129–151, https://doi.org/10.1051/m2an/197408R201291, 1974. a
Bueler, E.: Performance analysis of high-resolution ice-sheet simulations, J. Glaciol., 69, 1–6, https://doi.org/10.1017/jog.2022.113, 2022. a, b
Bueler, E., Lingle, C. S., Kallen-Brown, J. A., Covey, D. N., and Bowman, L. N.: Exact solutions and verification of numerical models for isothermal ice sheets, J. Glaciol, 51, 291–306, https://doi.org/10.3189/172756505781829449, 2005. a
Cuffey, K. M. and Paterson, W. S. B.: The Physics of Glaciers, Butterworth-Heinemann, Amsterdam, 4th Edn., ISBN 978-0-12-369461-4, 2010. a
DeConto, R. M. and Pollard, D.: Contribution of Antarctica to past and future sea-level rise, Nature, 531, 591–597, https://doi.org/10.1038/nature17145, 2016. a
Dukowicz, J. K., Price, S. F., and Lipscomb, W. H.: Incorporating arbitrary basal topography in the variational formulation of ice-sheet models, J. Glaciol., 57, 461–467, https://doi.org/10.3189/002214311796905550, 2011. a
Duretz, T., May, D. A., Gerya, T. V., and Tackley, P. J.: Discretization errors and free surface stabilization in the finite difference and marker-in-cell method for applied geodynamics: A numerical study, Geochem. Geophy. Geosy., 12, Q07004, https://doi.org/10.1029/2011GC003567, 2011. a
Fox-Kemper, B., Hewitt, H., Xiao, C., Aðalgeirsdóttir, G., Drijfhout, S., Edwards, T., Golledge, N., Hemer, M., Kopp, R., Krinner, G., Mix, A., Notz, D., Nowicki, S., Nurhati, I., Ruiz, L., Sallée, J.-B., Slangen, A., and Yu, Y.: Ocean, Cryosphere and Sea Level Change, Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 1211–1362, https://doi.org/10.1017/9781009157896.011, 2021. a
Franca, L. P. and Frey, S.: Stabilized finite element methods: II, the incompressible Navier-Stokes equations, Comput. Method. Appl. M., 99, 209–233, https://doi.org/10.1016/0045-7825(93)90119-I, 1992. 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, b, c, d
Glen, J. W.: The Creep of Polycrystalline Ice, P. R. Soc. London, 228, 519–538, https://doi.org/10.1098/rspa.1955.0066, 1955. a, b
Glerum, A., Brune, S., Stamps, D. S., and Strecker, M. R.: Victoria continental microplate dynamics controlled by the lithospheric strenght distribution of the East African Rift, Nat. Commun., 11, 2881, https://doi.org/10.1038/s41467-020-16176-x, 2020. a
Hanna, E., Navarro, F. J., Pattyn, F., Domingues, C. M., Fettweis, X., Ivins, E. R., Nicholls, R. J., Ritz, C., Smith, B., Tulaczyk, S., Whitehouse, P. L., and Zwally, H. J.: Ice-sheet mass balance and climate change, Nature, 498, 51–59, https://doi.org/10.1038/nature12238, 2013. a
Hindmarsh, R. C. A. and Payne, A. J.: Time-step limits for stable solutions of the ice-sheet equation, Ann. Glaciol., 23, 74–85, https://doi.org/10.3189/S0260305500013288, 1996. a
Hock, R., Bliss, A., Marzeion, B., Giesen, R. H., Hirabayashi, Y., Huss, M., Radić, V., and Slangen, A. B. A.: GlacierMIP – A model intercomparison of global-scale glacier mass-balance models and projections, J. Glaciol., 65, 453–467, https://doi.org/10.1017/jog.2019.22, 2019. a
Hughes, T. J. R., Franca, L. P., and Balestra, M.: A new finite element formulation for computational fluid dynamics: V. Circumventing the Babuska–Brezzi condition: a stable Petrov–Galerkin formulation of the stokes problem accommodating equal-order interpolations, Comput. Meth. Applt. Mech. Eng., 59, 85–99, https://doi.org/10.1016/0045-7825(86)90025-3, 1986. a
Hutter, K.: Theoretical Glaciology: Material Science of Ice and the Mechanics of Glaciers and Ice Sheets, Mathematical Approaches to Geophysics 1, D. Reidel Publishing Company, Dordrecht, Holland and Terra Scientific Publishing Company, Tokyo, Japan, https://doi.org/10.1007/978-94-015-1167-4, 1983. a
James, T. D., Murray, T., Barrand, N. E., and Barr, S. L.: Extracting photogrammetric ground control from lidar DEMs for change detection, Photogramm. Rec., 21, 312–328, 2006. a
Kaus, B. J., Mühlhaus, H., and May, D. A.: A stabilization algorithm for geodynamic numerical simulations with a free surface, Phys. Earth Planet. In., 181, 12–20, https://doi.org/10.1016/j.pepi.2010.04.007, 2010. a, b, c
Kohler, J., James, T. D., Murray, T., Nuth, C., Brandt, O., Barrand, N. E., Aas, H. F., and Luckman, A.: Acceleration in thinning rate on western Svalbard glaciers, Geophys. Res. Lett., 34, L18502, https://doi.org/10.1029/2007GL030681, 2007. a
Kramer, S. C., Wilson, C. R., and Davies, D. R.: An implicit free surface algorithm for geodynamical simulations, Phys. Earth Planet. In., 194–195, 25–37, https://doi.org/10.1016/j.pepi.2012.01.001, 2012. a
Ladyzhenskaya, O. A.: The Mathematical Theory of Viscous Incompressible Flow, Mathematic and its applications, Gordon and Breach New York, 2nd Edn., ISBN-10 067720230X, ISBN-13 9780677202303, 1969. a
Leng, W., Ju, L., Gunzburger, M., Price, S., and Ringler, T.: A parallel high-order accurate finite element nonlinear Stokes ice sheet model and benchmark experiments, J. Geophys. Res., 117, F01001, https://doi.org/10.1029/2011JF001962, 2012. a
Löfgren, A.: Perlin glacier data, Zenodo [data set], https://doi.org/10.5281/zenodo.11643074, 2024. a
Meredith, M., Sommerkorn, M., Cassotta, S., Derksen, C., Ekaykin, A., Hollowed, A., Kofinas, G., Mackintosh, A., Melbourne-Thomas, J., Muelbert, M., Ottersen, G., Pritchard, H., and Schuur, E.: Polar Regions, in: IPCC Special Report on the Ocean and Cryosphere in a Changing Climate, Cambridge University Press, 203–320, https://doi.org/10.1017/9781009157964.005, 2019. a
Meur, E. L., Gagliardini, O., Zwinger, T., and Ruokolainen, J.: Glacier flow modelling: a comparison of the Shallow Ice Approximation and the full-Stokes solution, C. R. Phys., 5, 709–722, https://doi.org/10.1016/j.crhy.2004.10.001, 2004. a
Morland, L. W.: Thermomechanical balances of ice sheet flows, Geophys. Astro Fluid, 29, 237–266, https://doi.org/10.1080/03091928408248191, 1984. a
Musgrave, F. K., Kolb, C. E., and Mace, R. S.: The Synthesis and Rendering of Eroded Fractal Terrain, ACM SIGGRAPH Comput. Graph., 23, 41–50, https://doi.org/10.1145/74334.74337, 1989. a
Norwegian Polar Institute: Norwegian Polar Institute Map Data and Services, https://geodata.npolar.no, last access: 14 February 2024. a
Nye, J. F.: The distribution of stress and velocity in glaciers and ice-sheets, Proc. R. Soc. Lon. Ser.-A, 239, 113–133, https://doi.org/10.1098/rspa.1957.0026, 1957. a
Paterson, W. S. B.: The Physics of Glaciers, Elsevier Science Ltd, The Boulevard, Langford Lane, Kidlington, Oxford, 3rd Edn., ISBN 978-0-08-037944-9, https://doi.org/10.1016/C2009-0-14802-X, 1994. 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
Pattyn, F.: The paradigm shift in Antarctic ice sheet modelling, Nat. Commun., 9, 2728, https://doi.org/10.1038/s41467-018-05003-z, 2018. a
Perlin, K.: An Image Synthesizer, ACM SIGGRAPH Comput. Graph., 19, 287–296, https://doi.org/10.1145/325165.325247, 1985. a, b
Rippin, D., Willis, I., Arnold, N., Hodson, A., Moore, J., Kohler, J., and Björnsson, H.: Changes in geometry and subglacial drainage of Midre Lovénbreen, Svalbard, determined from digital elevation models, Earth Surf. Proc. Land., 28, 273–298, https://doi.org/10.1002/esp.485, 2003. a
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, b, c
Rose, I., Buffet, B., and Heister, T.: Stability and accuracy of free surface time integration in viscous flows, Phys. Earth Planet. In., 262, 90–100, https://doi.org/10.1016/j.pepi.2016.11.007, 2017. a
Shepherd, A. and Nowicki, S.: Improvements in ice-sheet sea-level projections, Nat. Clim. Change, 7, 672–674, https://doi.org/10.1038/nclimate3400, 2017. a
Sullivan, C. B. and Kaszynski, A.: PyVista: 3D plotting and mesh analysis through a streamlined interface for the Visualization Toolkit (VTK), J. Open Source Softw., 4, 1450, https://doi.org/10.21105/joss.01450, 2019. a, b, c
Taylor, C. and Hood, P.: Numerical solution of the Navier-Stokes equations using the finite element technique, Comput. Fluids, 1, 73–100, https://doi.org/10.1016/0045-7930(73)90027-3, 1973. a
UK NERC Airborne Research and Survey Facility: Midtre Lovénbreen lidar data collected by UK NERC Airborne Research and Survey Facility, https://data.ceda.ac.uk/neodc/arsf/2003/03_04, last access: 13 June 2023. a
Weertman, J.: On the sliding of glaciers, J. Glaciol., 3, 38–42, https://doi.org/10.3189/S0022143000024709, 1957. a
Wirbel, A. and Jarosch, A. H.: Inequality-constrained free-surface evolution in a full Stokes ice flow model (evolve_glacier v1.1), Geosci. Model Dev., 13, 6425–6445, https://doi.org/10.5194/gmd-13-6425-2020, 2020. a
Zwinger, T. and Moore, J. C.: Diagnostic and prognostic simulations with a full Stokes model accounting for superimposed ice of Midtre Lovénbreen, Svalbard, The Cryosphere, 3, 217–229, https://doi.org/10.5194/tc-3-217-2009, 2009. a, b
Zwinger, T., Greve, R., Gagliardini, O., Shiraiwa, T., and Lyly, M.: A full Stokes-flow thermo-mechanical model for firn and ice applied to the Gorshkov crater glacier, Kamchatka, Ann. Glaciol., 45, 29–37, https://doi.org/10.3189/172756407782282543, 2007. a
Short summary
This paper investigates a stabilization method for free-surface flows in the context of glacier simulations. Previous applications of the stabilization on ice flows have only considered simple ice-sheet benchmark problems; in particular the method had not been tested on real-world glacier domains. This work addresses this shortcoming by demonstrating that the stabilization works well also in this case and increases stability and robustness without negatively impacting computation times.
This paper investigates a stabilization method for free-surface flows in the context of glacier...
