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
Josefin Ahlkrona, A. Clara J. Henry, and André Löfgren
EGUsphere, https://doi.org/10.5194/egusphere-2025-4359, https://doi.org/10.5194/egusphere-2025-4359, 2025
This preprint is open for discussion and under review for Geoscientific Model Development (GMD).
Short summary
Short summary
This paper leverages the Free Surface Stabilization Algorithm of Kaus et al. (2010) to construct the first fully implicit discretization of viscous free surface flows. It also presents the first second order accurate time-stepping scheme applicable to ice sheet models. We test the new method on an idealized problem and on a 2D glacier simulation. The results indicates that the method has great potential to speedup ice sheet models.
Josefin Ahlkrona, A. Clara J. Henry, and André Löfgren
EGUsphere, https://doi.org/10.5194/egusphere-2025-4359, https://doi.org/10.5194/egusphere-2025-4359, 2025
This preprint is open for discussion and under review for Geoscientific Model Development (GMD).
Short summary
Short summary
This paper leverages the Free Surface Stabilization Algorithm of Kaus et al. (2010) to construct the first fully implicit discretization of viscous free surface flows. It also presents the first second order accurate time-stepping scheme applicable to ice sheet models. We test the new method on an idealized problem and on a 2D glacier simulation. The results indicates that the method has great potential to speedup ice sheet models.
Cyrille Mosbeux, Peter Råback, Adrien Gilbert, Julien Brondex, Fabien Gillet-Chaulet, Nicolas C. Jourdain, Mondher Chekki, Olivier Gagliardini, and Gaël Durand
EGUsphere, https://doi.org/10.5194/egusphere-2025-3039, https://doi.org/10.5194/egusphere-2025-3039, 2025
This preprint is open for discussion and under review for Geoscientific Model Development (GMD).
Short summary
Short summary
Transport processes like rocks carried by ice flow and damage evolution – a proxy for crevasses – are key in ice sheet modeling and should occur without diffusion. Yet, standard numerical methods often blur these features. We explore a particle-based Semi-Lagrangian approach, comparing it to a Discontinuous Galerkin method, and show it can accurately simulate such transport when run at high enough resolution.
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.
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...
