Articles | Volume 19, issue 16
https://doi.org/10.5194/gmd-19-7817-2026
© Author(s) 2026. 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-19-7817-2026
© Author(s) 2026. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Development and technical paper
| 
21 Aug 2026
Development and technical paper |
| 21 Aug 2026
| 21 Aug 2026
Grounding-line dynamics in a Stokes ice-flow model (Elmer/Ice v9.0): improved numerical stability allows larger time steps
A. Clara J. Henry
CORRESPONDING AUTHOR
Department of Mathematics, Stockholm University, Stockholm, Sweden
Bolin Centre for Climate Research, Stockholm University, Stockholm, Sweden
now at: British Antarctic Survey, Cambridge, UK
Thomas Zwinger
CSC–IT Center for Science, Espoo, Finland
Josefin Ahlkrona
Department of Mathematics, Stockholm University, Stockholm, Sweden
Bolin Centre for Climate Research, Stockholm University, Stockholm, Sweden
Swedish e-Science Research Centre, Stockholm, Sweden
Related authors
No articles found.
James R. Jordan, Frank Pattyn, Daniel Abele, Torsten Albrecht, Jorge Alvarez-Solas, Jowan M. Barnes, Tijn Berends, Jorge A. Bernales, Javier Blasco, Gong Cheng, Youngmin Choi, Stephen L. Cornford, Cruz Garcia-Molina, Fabien Gillet-Chaulet, G. Hilmar Gudmundsson, Angelika Humbert, Gunter R. Leguy, William H. Lipscomb, Marisa Montoya, Daniel Moreno-Parada, Mathieu Morlighem, Tyler Pelle, Alexander Robinson, Martin Rückamp, Hélène Seroussi, Yanmei Tian, Luisa Wagner, Roderick S. W. van de Wal, Liyun Zhao, and Thomas Zwinger
EGUsphere, https://doi.org/10.5194/egusphere-2026-1962, https://doi.org/10.5194/egusphere-2026-1962, 2026
Short summary
Short summary
Ice calving is a recent inclusion in ice sheet models and there has not been a systematic test of how well it has been implemented. We show that models can accurately represent a given rate of ice calving, properties at the ice front evolve smoothly during calving, and that there are no consistent differences in ice behaviour between various approaches to representing calving in numerical ice models. This gives us confidence in their ability for use in sea level rise prediction simulations.
Josefin Ahlkrona, A. Clara J. Henry, and André Löfgren
Geosci. Model Dev., 19, 2333–2348, https://doi.org/10.5194/gmd-19-2333-2026, https://doi.org/10.5194/gmd-19-2333-2026, 2026
Short summary
Short summary
This paper leverages the Free Surface Stabilization Algorithm of Kaus et al. (2010) to construct a 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.
Yiliang Ma, Liyun Zhao, Rupert Gladstone, Thomas Zwinger, Michael Wolovick, Junshun Wang, and John C. Moore
The Cryosphere, 19, 6187–6205, https://doi.org/10.5194/tc-19-6187-2025, https://doi.org/10.5194/tc-19-6187-2025, 2025
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 an important impact on ice sheet dynamics. We simulate the evolution of Totten Glacier using an ice flow model with different basal sliding parameterizations and sub-shelf melt rates to quantify their effect on the projections. We found that the modelled glacier retreat and mass loss are sensitive to the choice of basal sliding parameterizations and maximal sub-shelf melt rate.
Johanna Beckmann, Ronja Reese, Felicity S. McCormack, Sue Cook, Lawrence Bird, Dawid Gwyther, Daniel Richards, Matthias Scheiter, Yu Wang, Hélène Seroussi, Ayako Abe‐Ouchi, Torsten Albrecht, Jorge Alvarez‐Solas, Xylar S. Asay‐Davis, Jean‐Baptiste Barre, Constantijn J. Berends, Jorge Bernales, Javier Blasco, Justine Caillet, David M. Chandler, Violaine Coulon, Richard Cullather, Christophe Dumas, Benjamin K. Galton‐Fenzi, Julius Garbe, Fabien Gillet‐Chaulet, Rupert Gladstone, Heiko Goelzer, Nicholas R. Golledge, Ralf Greve, G. Hilmar Gudmundsson, Holly Kyeore Han, Trevor R. Hillebrand, Matthew J. Hoffman, Philippe Huybrechts, Nicolas C. Jourdain, Ann Kristin Klose, Petra M. Langebroek, Gunter R. Leguy, William H. Lipscomb, Daniel P. Lowry, Pierre Mathiot, Marisa Montoya, Mathieu Morlighem, Sophie Nowicki, Frank Pattyn, Antony J. Payne, Tyler Pelle, Aurélien Quiquet, Alexander Robinson, Leopekka Saraste, Erika G. Simon, Sainan Sun, Jake P. Twarog, Luke D. Trusel, Benoit Urruty, Jonas Van Breedam, Roderik S. W. van de Wal, Chen Zhao, and Thomas Zwinger
EGUsphere, https://doi.org/10.5194/egusphere-2025-4069, https://doi.org/10.5194/egusphere-2025-4069, 2025
Short summary
Short summary
Antarctica holds enough ice to raise sea levels by many meters, but its future is uncertain. Warm ocean water melts ice shelves from below, letting inland ice flow faster into the sea. By 2300, Antarctica could add 0.6–4.4 m to sea levels. Our study identifies two key factors—how strongly shelves melt and how the ice responds. These explain much of the range, and refining them in models may improve future predictions.
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.
André Löfgren, Thomas Zwinger, Peter Råback, Christian Helanow, and Josefin Ahlkrona
The Cryosphere, 18, 3453–3470, https://doi.org/10.5194/tc-18-3453-2024, https://doi.org/10.5194/tc-18-3453-2024, 2024
Short summary
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.
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.
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.
Cited articles
Adusumilli, S., Fricker, H. A., Medley, B., Padman, L., and Siegfried, M. R.: Interannual variations in meltwater input to the Southern Ocean from Antarctic ice shelves, Nat. Geosci., 13, 616–620, https://doi.org/10.1038/s41561-020-0616-z, 2020. a, b
Ahlkrona, J., Lötstedt, P., Kirchner, N., and Zwinger, T.: Dynamically coupling the non-linear Stokes equations with the shallow ice approximation in glaciology: Description and first applications of the ISCAL method, J. Comput. Phys., 308, 1–19, https://doi.org/10.1016/j.jcp.2015.12.025, 2016. a
Ahlkrona, J., Henry, A. C. J., and Löfgren, A.: A fully implicit second order method for viscous free surface Stokes flow – application to glacier simulations, Geosci. Model Dev., 19, 2333–2348, https://doi.org/10.5194/gmd-19-2333-2026, 2026. a, b
Amestoy, P., Duff, I. S., Koster, J., and L'Excellent, J.-Y.: A fully asynchronous multifrontal solver using distributed dynamic scheduling, SIAM J. Matrix Anal. A., 23, 15–41, https://doi.org/10.1137/S0895479899358194, 2001. 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
Bett, D. T., Bradley, A. T., Williams, C. R., Holland, P. R., Arthern, R. J., and Goldberg, D. N.: Coupled ice–ocean interactions during future retreat of West Antarctic ice streams in the Amundsen Sea sector, The Cryosphere, 18, 2653–2675, https://doi.org/10.5194/tc-18-2653-2024, 2024. a
Bingham, R. G., Bodart, J. A., Cavitte, M. G. P., Chung, A., Sanderson, R. J., Sutter, J. C. R., Eisen, O., Karlsson, N. B., MacGregor, J. A., Ross, N., Young, D. A., Ashmore, D. W., Born, A., Chu, W., Cui, X., Drews, R., Franke, S., Goel, V., Goodge, J. W., Henry, A. C. J., Hermant, A., Hills, B. H., Holschuh, N., Koutnik, M. R., Leysinger Vieli, G. J.-M. C., MacKie, E. J., Mantelli, E., Martín, C., Ng, F. S. L., Oraschewski, F. M., Napoleoni, F., Parrenin, F., Popov, S. V., Rieckh, T., Schlegel, R., Schroeder, D. M., Siegert, M. J., Tang, X., Teisberg, T. O., Winter, K., Yan, S., Davis, H., Dow, C. F., Fudge, T. J., Jordan, T. A., Kulessa, B., Matsuoka, K., Nyqvist, C. J., Rahnemoonfar, M., Siegfried, M. R., Singh, S., Višnjević, V., Zamora, R., and Zuhr, A.: Review article: AntArchitecture – building an age–depth model from Antarctica's radiostratigraphy to explore ice-sheet evolution, The Cryosphere, 19, 4611–4655, https://doi.org/10.5194/tc-19-4611-2025, 2025. a
Brondex, J., Gagliardini, O., Gillet-Chaulet, F., and Durand, G.: Sensitivity of grounding line dynamics to the choice of the friction law, J. Glaciol., 63, 854–866, https://doi.org/10.1017/jog.2017.51, 2017. a
Brondex, J., Gillet-Chaulet, F., and Gagliardini, O.: Sensitivity of centennial mass loss projections of the Amundsen basin to the friction law, The Cryosphere, 13, 177–195, https://doi.org/10.5194/tc-13-177-2019, 2019. a
Bueler, E.: Surface elevation errors in finite element Stokes models for glacier evolution, arXiv [preprint], https://doi.org/10.48550/arXiv.2408.06470, 2024. a
Cheng, G., Lötstedt, P., and von Sydow, L.: Accurate and stable time stepping in ice sheet modeling, J. Comput. Phys., 329, 29–47, https://doi.org/10.1016/j.jcp.2016.10.060, 2017. a, b
Cheng, G., Lötstedt, P., and von Sydow, L.: A full Stokes subgrid scheme in two dimensions for simulation of grounding line migration in ice sheets using Elmer/ICE (v8.3), Geosci. Model Dev., 13, 2245–2258, https://doi.org/10.5194/gmd-13-2245-2020, 2020. a, b
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, b
Donat-Magnin, M., Jourdain, N. C., Kittel, C., Agosta, C., Amory, C., Gallée, H., Krinner, G., and Chekki, M.: Future surface mass balance and surface melt in the Amundsen sector of the West Antarctic Ice Sheet, The Cryosphere, 15, 571–593, https://doi.org/10.5194/tc-15-571-2021, 2021. 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
Durand, G., Gagliardini, O., Zwinger, T., Le Meur, E., and Hindmarsh, R. C.: Full Stokes modeling of marine ice sheets: influence of the grid size, Ann. Glaciol., 50, 109–114, https://doi.org/10.3189/172756409789624283, 2009b. a, b
Eisermann, H., Eagles, G., Ruppel, A., Smith, E. C., and Jokat, W.: Bathymetry beneath ice shelves of Western Dronning Maud Land, East Antarctica, and implications on ice shelf stability, Geophys. Res. Lett., 47, e2019GL086724, https://doi.org/10.1029/2019GL086724, 2020. a
Fowler, A.: A sliding law for glaciers of constant viscosity in the presence of subglacial cavitation, P. Roy. Soc. A-Math. Phy., 407, 147–170, https://doi.org/10.1098/rspa.1986.0090, 1986. a
Fox-Kemper, B., Hewitt, H. T., 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, in: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S. L., Péan, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M. I., Huang, M., Leitzell, K., Lonnoy, E., Matthews, J. B. R., Maycock, T. K., Waterfield, T., Yelekçi, O., Yu, R., and Zhou, B., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, pp. 1211–1362, https://doi.org/10.1017/9781009157896.011, 2021. a
Franca, L. P. and Frey, S. L.: Stabilized finite element methods: II. The incompressible Navier–Stokes equations, Comput Method. Appl. M., 99, 209–233, https://doi.org/10.1016/0045-7825(92)90041-H, 1992. a
Gagliardini, O., Cohen, D., Råback, P., and Zwinger, T.: Finite-element modeling of subglacial cavities and related friction law, J. Geophys. Res., 112, F02027, https://doi.org/10.1029/2006JF000576, 2007. a
Gagliardini, O., Zwinger, T., Gillet-Chaulet, F., Durand, G., Favier, L., de Fleurian, B., Greve, R., Malinen, M., Martín, C., Råback, P., Ruokolainen, J., Sacchettini, M., Schäfer, M., Seddik, H., and Thies, J.: Capabilities and performance of Elmer/Ice, a new-generation ice sheet model, Geosci. Model Dev., 6, 1299–1318, https://doi.org/10.5194/gmd-6-1299-2013, 2013. a, b, c, d
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, b, c, d
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
Gladstone, R., Lee, V., Vieli, A., and Payne, A.: Grounding line migration in an adaptive mesh ice sheet model, J. Geophys. Res.-Earth, 115, 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, https://doi.org/10.3189/2012AoG60A148, 2012. a
Goldberg, D., Holland, D. M., and Schoof, C.: Grounding line movement and ice shelf buttressing in marine ice sheets, J. Geophys. Res., 114, F04026, https://doi.org/10.1029/2008JF001227, 2009. a
Helanow, C. and Ahlkrona, J.: Stabilized equal low-order finite elements in ice sheet modeling – accuracy and robustness, Computat. Geosci., 22, 951–974, https://doi.org/10.1007/s10596-017-9713-5, 2018. a
Henry, A., Schannwell, C., Višnjević, V., Millstein, J., Bons, P., Eisen, O., and Drews, R.: Predicting the three-dimensional stratigraphy of an ice rise, J. Geophys. Res.-Earth, 130, e2024JF007924, https://doi.org/10.1029/2024JF007924, 2025. a
Henry, A. C. J.: Model output data for the publication “Grounding-line dynamics in a Stokes ice-flow model: Improved numerical stability allows larger time steps”, Zenodo [data set], https://doi.org/10.5281/zenodo.16909985, 2025. a
Henry, A. C. J.: Code for the publication “Grounding-line dynamics in a Stokes ice-flow model: Improved numerical stability allows larger time steps”, Zenodo [code], https://doi.org/10.5281/zenodo.21681869, 2026. a, b
John, V., Knobloch, P., and Novo, J.: Finite elements for scalar convection-dominated equations and incompressible flow problems: a never ending story?, Computing and Visualization in Science, 19, 47–63, https://doi.org/10.1007/s00791-018-0290-5, 2018. 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, d
Kramer, S. C., Wilson, C. R., and Davies, D. R.: An implicit free surface algorithm for geodynamical simulations, Phys. Earth Planet. In., 194, 25–37, https://doi.org/10.1016/j.pepi.2012.01.001, 2012. a
Larour, E., Seroussi, H., Morlighem, M., and Rignot, E.: Continental scale, high order, high spatial resolution, ice sheet modeling using the Ice Sheet System Model (ISSM), J. Geophys. Res.-Earth, 117, https://doi.org/10.1029/2011JF002140, 2012. a
Löfgren, A., Zwinger, T., Råback, P., Helanow, C., and Ahlkrona, J.: Increasing numerical stability of mountain valley glacier simulations: implementation and testing of free-surface stabilization in Elmer/Ice, The Cryosphere, 18, 3453–3470, https://doi.org/10.5194/tc-18-3453-2024, 2024. a, b, c, d
Morlighem, M., Williams, C. N., Rignot, E., An, L., Arndt, J. E., Bamber, J. L., Catania, G., Chauché, N., Dowdeswell, J. A., Dorschel, B., Fenty, I., Hogan, K., Howat, I., Hubbard, A., Jakobsson, M., Jordan, T. M., Kjeldsen, K. K., Millan, R., Mayer, L., Mouginot, J., Noël, B. P. Y., O'Cofaigh, C., Palmer, S., Rysgaard, S., Seroussi, H., Siegert, M. J., Slabon, P., Straneo, F., van den Broeke, M. R., Weinrebe, W., Wood, M., and Zinglersen, K. B.: BedMachine v3: Complete bed topography and ocean bathymetry mapping of Greenland from multibeam echo sounding combined with mass conservation, Geophys. Res. Lett., 44, 11–051, https://doi.org/10.1002/2017GL074954, 2017. a
Morlighem, M., Rignot, E., Binder, T., Blankenship, D., Drews, R., Eagles, G., Eisen, O., Ferraccioli, F., Forsberg, R., Fretwell, P., Goel, V., Greenbaum, J. S., Gudmundsson, H., Guo, J., Helm, V., Hofstede, C., Howat, I., Humbert, A., Jokat, W., Karlsson, N. B., Lee, W. S., Matsuoka, K., Millan, R., Mouginot, J., Paden, J., Pattyn, F., Roberts, J., Rosier, S., Ruppel, A., Seroussi, H., Smith, E. C., Steinhage, D., Sun, B., van den Broeke, M. R., van Ommen, T. D., van Wessem, M., and Young, D. A.: Deep glacial troughs and stabilizing ridges unveiled beneath the margins of the Antarctic ice sheet, Nat. Geosci., 13, 132–137, https://doi.org/10.1038/s41561-019-0510-8, 2020. a, b, c, d
Moss, G., Višnjević, V., Eisen, O., Oraschewski, F. M., Schröder, C., Macke, J. H., and Drews, R.: Simulation-based inference of surface accumulation and basal melt rates of an Antarctic ice shelf from isochronal layers, J. Glaciol., 71, e44, https://doi.org/10.1017/jog.2025.13, 2025. a
Neckel, N., Drews, R., Rack, W., and Steinhage, D.: Basal melting at the Ekström Ice Shelf, Antarctica, estimated from mass flux divergence, Ann. Glaciol,, 53, 294–302, https://doi.org/10.3189/2012AoG60A167, 2012. a
Oetting, A., Smith, E. C., Arndt, J. E., Dorschel, B., Drews, R., Ehlers, T. A., Gaedicke, C., Hofstede, C., Klages, J. P., Kuhn, G., Lambrecht, A., Läufer, A., Mayer, C., Tiedemann, R., Wilhelms, F., and Eisen, O.: Geomorphology and shallow sub-sea-floor structures underneath the Ekström Ice Shelf, Antarctica, The Cryosphere, 16, 2051–2066, https://doi.org/10.5194/tc-16-2051-2022, 2022. a
Pattyn, F., Schoof, C., Perichon, L., Hindmarsh, R. C. A., Bueler, E., de Fleurian, B., Durand, G., Gagliardini, O., Gladstone, R., Goldberg, D., Gudmundsson, G. H., Huybrechts, P., Lee, V., Nick, F. M., Payne, A. J., Pollard, D., Rybak, O., Saito, F., and Vieli, A.: Results of the Marine Ice Sheet Model Intercomparison Project, MISMIP, The Cryosphere, 6, 573–588, https://doi.org/10.5194/tc-6-573-2012, 2012. a, b, c, d, e, f, g, h, i
Rignot, E., Mouginot, J., and Scheuchl., B.: MEaSUREs InSAR-based Antarctica ice velocity map, Version 2, https://doi.org/10.5067/D7GK8F5J8M8R, 2017. a
Rückamp, M., Kleiner, T., and Humbert, A.: Comparison of ice dynamics using full-Stokes and Blatter–Pattyn approximation: application to the Northeast Greenland Ice Stream, The Cryosphere, 16, 1675–1696, https://doi.org/10.5194/tc-16-1675-2022, 2022. a
Ruokolainen, J., Malinen, M., Råback, P., Zwinger, T., Takala, E., Kataja, J., Gillet-Chaulet, F., Ilvonen, S., Gladstone, R., Byckling, M., Chekki, M., Gong, C., Ponomarev, P., van Dongen, E., Robertsen, F., Wheel, I., Cook, S., t7saeki, luzpaz, and Rich_B: ElmerCSC/elmerfem: Elmer 9.0, Zenodo [code], https://doi.org/10.5281/zenodo.7892181, 2023. a, b
Schannwell, C., Drews, R., Ehlers, T. A., Eisen, O., Mayer, C., and Gillet-Chaulet, F.: Kinematic response of ice-rise divides to changes in ocean and atmosphere forcing, The Cryosphere, 13, 2673–2691, https://doi.org/10.5194/tc-13-2673-2019, 2019. a
Schoof, C.: The effect of cavitation on glacier sliding, P. Roy. Soc. A-Math. Phy., 461, 609–627, https://doi.org/10.1098/rspa.2004.1350, 2005. a
Schoof, C.: Ice sheet grounding line dynamics: Steady states, stability, and hysteresis, J. Geophys. Res., 112, F03S28, https://doi.org/10.1029/2006JF000664, 2007. a, b, c
Seddik, H., Greve, R., Zwinger, T., and Sugiyama, S.: Regional modeling of the Shirase drainage basin, East Antarctica: full Stokes vs. shallow ice dynamics, The Cryosphere, 11, 2213–2229, https://doi.org/10.5194/tc-11-2213-2017, 2017. a
Sergienko, O. and Haseloff, M.: “Stable” and “unstable” are not useful descriptions of marine ice sheets in the Earth's climate system, J. Glaciol., 69, 1483–1499, https://doi.org/10.1017/jog.2023.40, 2023. a
Sergienko, O. V.: No general stability conditions for marine ice-sheet grounding lines in the presence of feedbacks, Nat. Commun., 13, 2265, https://doi.org/10.1038/s41467-022-29892-3, 2022. a
Seroussi, H., Ben Dhia, H., Morlighem, M., Larour, E., Rignot, E., and Aubry, D.: Coupling ice flow models of varying orders of complexity with the Tiling method, J. Glaciol., 58, 776–786, https://doi.org/10.3189/2012JoG11J195, 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
Seroussi, H., Nowicki, S., Simon, E., Abe-Ouchi, A., Albrecht, T., Brondex, J., Cornford, S., Dumas, C., Gillet-Chaulet, F., Goelzer, H., Golledge, N. R., Gregory, J. M., Greve, R., Hoffman, M. J., Humbert, A., Huybrechts, P., Kleiner, T., Larour, E., Leguy, G., Lipscomb, W. H., Lowry, D., Mengel, M., Morlighem, M., Pattyn, F., Payne, A. J., Pollard, D., Price, S. F., Quiquet, A., Reerink, T. J., Reese, R., Rodehacke, C. B., Schlegel, N.-J., Shepherd, A., Sun, S., Sutter, J., Van Breedam, J., van de Wal, R. S. W., Winkelmann, R., and Zhang, T.: initMIP-Antarctica: an ice sheet model initialization experiment of ISMIP6, The Cryosphere, 13, 1441–1471, https://doi.org/10.5194/tc-13-1441-2019, 2019. a
Smith, E. C., Hattermann, T., Kuhn, G., Gaedicke, C., Berger, S., Drews, R., Ehlers, T. A., Franke, D., Gromig, R., Hofstede, C., Lambrecht, A., Läufer, A., Mayer, C., Tiedemann, R., Wilhelms, F., and Eisen, O.: Detailed seismic bathymetry beneath Ekström Ice Shelf, Antarctica: Implications for glacial history and ice-ocean interaction, Geophys. Res. Lett., 47, e2019GL086187, https://doi.org/10.1029/2019GL086187, 2020. a
Stubblefield, A. G., Spiegelman, M., and Creyts, T. T.: Variational formulation of marine ice-sheet and subglacial-lake grounding-line dynamics, J. Fluid Mech., 919, A23, https://doi.org/10.1017/jfm.2021.394, 2021. a
van den Broeke, M. R.: RACMO2.3p1 annual surface mass balance Antarctica (1979–2014), https://doi.org/10.1594/PANGAEA.896940, 2019. a, b
van Dongen, E. C. H., Kirchner, N., van Gijzen, M. B., van de Wal, R. S. W., Zwinger, T., Cheng, G., Lötstedt, P., and von Sydow, L.: Dynamically coupling full Stokes and shallow shelf approximation for marine ice sheet flow using Elmer/Ice (v8.3), Geosci. Model Dev., 11, 4563–4576, https://doi.org/10.5194/gmd-11-4563-2018, 2018. a
Višnjević, V., Moss, G., Henry, A. C. J., Wild, C. T., Steinhage, D., and Drews, R.: Mapping the composition of Antarctic ice shelves as a metric for their susceptibility to future climate change, Geophys. Res. Lett., 52, e2024GL112585, https://doi.org/10.1029/2024GL112585, 2025. a
Westling Dolling, T., Henry, A. C. J., and Ahlkrona, J.: A numerical stabilization scheme for the shallow shelf approximation, J. Comput. Phys., 565, 115182, https://doi.org/10.1016/j.jcp.2026.115182, 2026. a
Williams, C. R., Thodoroff, P., Arthern, R. J., Byrne, J., Hosking, J. S., Kaiser, M., Lawrence, N. D., and Kazlauskaite, I.: Calculations of extreme sea level rise scenarios are strongly dependent on ice sheet model resolution, Communications Earth and Environment, 6, 60, https://doi.org/10.1038/s43247-025-02010-z, 2025. a
Short summary
To overcome time-step restrictions, we implement the Free-Surface Stabilisation Algorithm (FSSA) at the ice-ocean interface in Stokes ice-sheet simulations. In 2D experiments, a time step of 10 years is generally numerically stable and accurate, whereas a time step of 50 years is stable, but cannot fully capture grounding-line dynamics. Implementation at the ice-ocean interface increases the applicability of Stokes models and motivates future coupling with adaptive time-stepping schemes.
To overcome time-step restrictions, we implement the Free-Surface Stabilisation Algorithm (FSSA)...
