Abstract
Numerical simulation of ice sheet dynamics is crucial for improved assessment and prediction of glacier mass balance and grounding line migration. For the quantitative simulation of continental-scale ice sheets in Antarctica and Greenland, ice flow models have been described using the Shallow Shelf Approximation (SSA), Shallow Ice Approximation (SIA), and combinations thereof. Here, we developed the COM-SSA, a finite element SSA module implemented in COMSOL Multiphysics® via the Weak Form module. To verify the applicability of the COM-SSA to ice sheet dynamics, five independent cases were modeled and compared with the Ice Sheet System Model (ISSM), which is a widely used ice sheet dynamics package. In Cases 1 and 2, which considered the nonlinearity problems of Glen’s viscosity and Weertman-type friction laws, we confirmed that the COM-SSA successfully handled the convergence of viscosity and basal friction with an increasing number of iterations. In Case 3, we tested the ice thickness evolution computed by solving the equation of mass conservation. Grounding line migration in response to ice thickness changes was investigated in Case 4. The COM-SSA showed a similar pattern of glacier geometry evolution compared to ISSM in Cases 3 and 4, indicating that the COM-SSA could precisely calculate ice sheet dynamics in a time-dependent manner. We demonstrated the capabilities of the COM-SSA in Case 5, to which the Pine Island Glacier in Antarctica was subjected. Similar to the ISSM, the COM-SSA predicted the velocity field, ice thickness, and grounding line for the Pine Island Glacier after 300 years of ablation.
Similar content being viewed by others
Data availability
COM-SSA module application files associated with this research are available and can be obtained by contacting the corresponding author.
References
Adhikari S, Ivins ER, Larour E, Seroussi H, Morlighem M, Nowicki S (2014) Future antarctic bed topography and its implications for ice sheet dynamics. Solid Earth 5:569–584. https://doi.org/10.5194/se-5-569-2014
Agrusta R, van Hunen J, Goes S (2018) Strong plates enhance mantle mixing in early Earth. Nat Commun 9:2708. https://doi.org/10.1038/s41467-018-05194-5
Ahrens J, Geveci B, law C (2005) Paraview: an end-user tool for large data visualization. The visualization handbook, New York, p 717
Amante C, Eakins BW (2009) ETOPO1 1 arc-minute global relief model: procedures, data sources and analysis, NOAA technical memorandum NESDIS NGDC-24. Natl Geophys Data Center, NOAA. https://doi.org/10.7289/V5C8276M
Bindschadler RA, Nowicki S, Abe-OUCHI A, Aschwanden A, Choi H, Fastook J, Granzow G, Greve R, Gutowski G, Herzfeld U, Jackson C, Johnson J, Khroulev C, Levermann A, Lipscomb WH, Martin MA, Morlighem M, Parizek BR, Pollard D, Wang WL (2013) Ice-sheet model sensitivities to environmental forcing and their use in projecting future sea level (the SeaRISE project). J Glaciol 59:195–224. https://doi.org/10.3189/2013JoG12J125
Blatter H (1995) 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
Bondzio JH, Seroussi H, Morlighem M, Kleiner T, Rückamp M, Humbert A, Larour EY (2016) Modelling calving front dynamics using a level-set method: application to Jakobshavn Isbræ, West Greenland. Cryosphere 10:497–510. https://doi.org/10.5194/tc-10-497-2016
Brondex J, Gagliardini O, Gillet-Chaulet F, Durand G (2017) 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
Bueler E, Brown J (2009) Shallow shelf approximation as a “sliding law” in a thermomechanically coupled ice sheet model. J Geophys Res 114:F03008. https://doi.org/10.1029/2008JF001179
Cole DM (2020) On the physical basis for the creep of ice: the high temperature regime. J Glaciol 66:401–414. https://doi.org/10.1017/jog.2020.15
Cook AJ, Copland L, Noël BPY, Stokes CR, Bentley MJ, Sharp MJ, Bingham RG, van den Broeke MR (2019) Atmospheric forcing of rapid marine-terminating glacier retreat in the Canadian Arctic Archipelago. Sci Adv 5:8507. https://doi.org/10.1126/sciadv.aau8507
dos Santos TD, Morlighem M, Seroussi H (2021) Assessment of numerical schemes for transient, finite-element ice flow models using ISSM v4.18. Geosci Mod Dev 14:2545–2573. https://doi.org/10.5194/gmd-14-2545-2021
Favier L, Gagliardini O, Durand G, Zwinger T (2012) A three-dimensional full Stokes model of the grounding line dynamics: effect of a pinning point beneath the ice shelf. Cryosphere 6:101–112. https://doi.org/10.5194/tc-6-101-2012
Favier L, Durand G, Cornford SL, Gudmundsson GH, Gagliardini O, Gillet-Chaulet F, Zwinger T, Payne AJ, Le Brocq AM (2014) Retreat of Pine Island Glacier controlled by marine ice-sheet instability. Nat Clim Change 4:117–121. https://doi.org/10.1038/nclimate2094
Fischler Y, Rückamp M, Bischof C, Aizinger V, Morlighem M, Humbert A (2022) A scalability study of the Ice-sheet and Sea-level System Model (ISSM, version 4.18). Geoscientific Model Dev 15:3753–3771. https://doi.org/10.5194/gmd-15-3753-2022
Fujita K, Ageta Y (2000) Effect of summer accumulation on glacier mass balance on the Tibetan Plateau revealed by mass-balance model. J Glaciol 46:244–252. https://doi.org/10.3189/172756500781832945
Fyke J, Sergienko O, Löfverström M, Price S, Lenaerts JTM (2018) An overview of interactions and feedbacks between ice sheets and the earth system. Rev Geophys 56:361–408. https://doi.org/10.1029/2018RG000600
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, Thies J (2013) Capabilities and performance of Elmer/Ice, a new-generation ice sheet model. Geosci Mod Dev 6:1299–1318. https://doi.org/10.5194/gmd-6-1299-2013
Glen JW, PRSL A (1955) The creep of polycrystalline ice. P R Soc London A 228:519–538. https://doi.org/10.1098/rspa.1955.0066
Hutter K (1982) Dynamics of glaciers and large ice masses. Ann Rev Fluid Mech 14:87–130
IPCC (2007) Climate change 2007: the physical science basis. Intergovernmental Panel on Climate Change, Cambridge, p 996
Jamieson SSR, Vieli A, Livingstone SJ, Cofaigh CÓ, Stokes C, Hillenbrand C-D, Dowdeswell JA (2012) Ice-stream stability on a reverse bed slope. Nat Geosci 5:799–802. https://doi.org/10.1038/ngeo1600
Jenkins A, Dutrieux P, Jacobs SS, McPhail SD, Perrett JR, Webb AT, White D (2010) Observations beneath Pine Island Glacier in West Antarctica and implications for its retreat. Nat Geosci 3:468–472. https://doi.org/10.1038/ngeo890
Khan SA, Choi Y, Morlighem M, Rignot E, Helm V, Humbert A, Mouginot J, Millan R, Kjær KH, Bjørk AA (2022) Extensive inland thinning and speed-up of Northeast Greenland Ice Stream. Nature 611:727–732. https://doi.org/10.1038/s41586-022-05301-z
King MD, Howat IM, Candela SG, Noh MJ, Jeong S, Noël BPY, van den Broeke MR, Wouters B, Negrete A (2020) Dynamic ice loss from the Greenland Ice Sheet driven by sustained glacier retreat. Commun Earth Environ 1:1–7. https://doi.org/10.1038/s43247-020-0001-2
Krenk S, Hededal O (1995) A dual orthogonality procedure for non-linear finite element equations. Comput Methods Appl Mech Eng 123:95–107. https://doi.org/10.1016/0045-7825(95)00779-Z
Kuslits LB, Farkas MP, Galsa A (2014) Effect of temperature-dependent viscosity on mantle convection. Acta Geod Geophys 49:249–263. https://doi.org/10.1007/s40328-014-0055-7
Larour E, Morlighem M, Seroussi H, Schiermeier J, Rignot E (2012a) Ice flow sensitivity to geothermal heat flux of Pine Island Glacier, Antarctica. J Geophys Res-Earth 117:1–12. https://doi.org/10.1029/2012JF002371
Larour E, Schiermeier J, Rignot E, Seroussi H, Morlighem M, Paden J (2012b) Sensitivity analysis of Pine Island Glacier ice flow using ISSM and DAKOTA. J Geophys Res-Earth 117:1–16. https://doi.org/10.1029/2011JF002146
MacAyeal DR (1989) Large-scale ice flow over a viscous basal sediment: theory and application to ice stream B, Antarctica. J Geophys Res-Sol Ea 94:4071–4087. https://doi.org/10.1029/JB094iB04p04071
MacAyeal DR (1992) The basal stress distribution of Ice Stream E, Antarctica, inferred by control methods. J Geophys Res 97:595. https://doi.org/10.1029/91JB02454
MacAyeal DR, Sergienko OV, Banwell AF, Macdonald GJ, Willis IC, Stevens LA (2021) Treatment of ice-shelf evolution combining flow and flexure. J Glaciol 67:885–902. https://doi.org/10.1017/jog.2021.39
Machguth H, Thomsen HH, Weidick A, Ahlstrøm AP, Abermann J, Andersen ML, Andersen SB, Bjørk AA, Box JE, Braithwaite RJ, Bøggild CE, Citterio M, Clement P, Colgan W, Fausto RS, Gleie K, Gubler S, Hasholt B, Hynek B, Van De Wal RSW (2016) Greenland surface mass-balance observations from the ice-sheet ablation area and local glaciers. J Glaciol 62:861–887. https://doi.org/10.1017/jog.2016.75
Maier N, Gimbert F, Gillet-Chaulet F, Gilbert A (2021) Basal traction mainly dictated by hard-bed physics over grounded regions of Greenland. Cryosphere 15:1435–1451. https://doi.org/10.5194/tc-15-1435-2021
Marzeion B, Hock R, Anderson B, Bliss A, Champollion N, Fujita K, Huss M, Immerzeel WW, Kraaijenbrink P, Malles J, Maussion F, Radić V, Rounce DR, Sakai A, Shannon S, van de Wal R, Zekollari H (2020) Partitioning the uncertainty of ensemble projections of global glacier mass change. Earth Fut 8:417–426. https://doi.org/10.1029/2019EF001470
Millstein JD, Minchew BM, Pegler SS (2022) Ice viscosity is more sensitive to stress than commonly assumed. Commun Ea Environ 3:1–7. https://doi.org/10.1038/s43247-022-00385-x
Morlighem M, Wood M, Seroussi H, Choi Y, Rignot E (2019) Modeling the response of northwest Greenland to enhanced ocean thermal forcing and subglacial discharge. Cryosphere 13:723–734. https://doi.org/10.5194/tc-13-723-2019
Morlighem M, Seroussi H, Larour É, Schlegel N, Borstad C, de Fleurian B (2015) Ice sheet system model 2015 (4.9) user guide. http://www.ccpo.odu.edu/∼klinck/Reprints/PDF/ISSMguide2015.pdf Accessed 15 Jan 2023
Muszynski I, Birchfield GE (1987) A coupled Marine Ice-Stream – Ice-shelf model. J Glaciol 33:3–15. https://doi.org/10.3189/s0022143000005281
Pattyn F (2003) 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
Pattyn F, Schoof C, Perichon L, Hindmarsh RCA, Bueler E, de Fleurian B, Durand G, Gagliardini O, Gladstone R, Goldberg D, Gudmundsson GH, Huybrechts P, Lee V, Nick FM, Payne AJ, Pollard D, Rybak O, Saito F, Vieli A (2012) Results of the marine ice sheet model intercomparison project, MISMIP. Cryosphere 6:573–588. https://doi.org/10.5194/tc-6-573-2012
Pelle T, Morlighem M, Nakayama Y, Seroussi H (2021) Widespread grounding line retreat of totten glacier, East Antarctica, over the 21st century. Geophys Res Lett 48:1–11. https://doi.org/10.1029/2021GL093213
Pollack HN, Chapman DS (1977) Mantle heat flow. Earth Planet Sci Lett 34:174–184. https://doi.org/10.1016/0012-821X(77)90002-4
Pollard D, Deconto RM (2012) Description of a hybrid ice sheet-shelf model, and application to Antarctica. Geoscientific Model Dev 5:1273–1295. https://doi.org/10.5194/gmd-5-1273-2012
Qi C, Goldsby DL (2021) An experimental investigation of the effect of grain size on “Dislocation Creep” of ice. J Geophys Res-Sol Ea 126:e2021JB0218. https://doi.org/10.1029/2021JB021824
Rignot E, Mouginot J, Scheuchl B (2011) Ice flow of the Antarctic Ice sheet. Science 333:1427. https://doi.org/10.1126/science.1208336
Rückamp M, Greve R, Humbert A (2019) Comparative simulations of the evolution of the Greenland ice sheet under simplified Paris agreement scenarios with the models SICOPOLIS and ISSM. Polar Sci 21:14–25. https://doi.org/10.1016/j.polar.2018.12.003
Schlegel NJ, Seroussi H, Schodlok MP, Larour EY, Boening C, Limonadi D, Watkins MM, Morlighem M, van den Broeke MR (2018) Exploration of Antarctic Ice Sheet 100 year contribution to sea level rise and associated model uncertainties using the ISSM framework. Cryosphere 12:3511–3534. https://doi.org/10.5194/tc-12-3511-2018
Seddik H, Greve R, Sakakibara D, Tsutaki S, Minowa M, Sugiyama S (2019) Response of the flow dynamics of Bowdoin Glacier, northwestern Greenland, to basal lubrication and tidal forcing. J Glaciol 65:225–238. https://doi.org/10.1017/jog.2018.106
Seroussi H, Morlighem M, Larour E, Rignot E, Khazendar A (2014) Hydrostatic grounding line parameterization in ice sheet models. Cryosphere 8:2075–2087. https://doi.org/10.5194/tc-8-2075-2014
Seroussi H, Nakayama Y, Larour E, Menemenlis D, Morlighem M, Rignot E, Khazendar A (2017) Continued retreat of Thwaites Glacier, West Antarctica, controlled by bed topography and ocean circulation. Geophys Res Lett 44:6191–6199. https://doi.org/10.1002/2017GL072910
Stansell ND, Licciardi JM, Rodbell DT, Mark BG (2017) Tropical ocean-atmospheric forcing of Late Glacial and Holocene glacier fluctuations in the Cordillera Blanca, Peru. Geophy Res Lett 44:4176–4185. https://doi.org/10.1002/2016GL072408
Tackley P (2000) Mantle convection and plate tectonics: toward an integrated physical and chemical theory. Science 288:2002–2007. https://doi.org/10.1126/science.288.5473.2002
Trim SJ, Butler SL, Spiteri RJ (2021) Benchmarking multiphysics software for mantle convection. Comput Geosci 154:104797. https://doi.org/10.1016/j.cageo.2021.104797
van Dongen ECH, Kirchner N, Van Gijzen MB, Van De Wal RSW, Zwinger T, Cheng G, Lötstedt P, Von Sydow L (2018) 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
Weertman J (1957) On the sliding of glaciers. J Glaciol 3:33–38. https://doi.org/10.3189/S0022143000024709
Weis M, Greve R, Hutter K (1999) Theory of shallow ice shelves. Contin Mech Thermodyn 11:15–50. https://doi.org/10.1007/s001610050102
Woidt WD (1978) Finite element calculations applied to salt dome analysis. Tectonophysics 50:369–386. https://doi.org/10.1016/0040-1951(78)90143-9
Zemp M, Huss M, Thibert E, Eckert N, McNabb R, Huber J, Barandun M, Machguth H, Nussbaumer SU, Gärtner-Roer I, Thomson L, Paul F, Maussion F, Kutuzov S, Cogley JG (2019) Global glacier mass changes and their contributions to sea-level rise from 1961 to 2016. Nature 568:382–386. https://doi.org/10.1038/s41586-019-1071-0
Acknowledgements
This research was supported by the National Research Foundation of Korea (NRF-2022R1A4A3027001, NRF-2022R1A2C1009742, and NRF-2019R1A6A1A03033167) and the Korea Institute of Marine Science & Technology Promotion (KIMST) funded by the Ministry of Oceans and Fisheries (RS-2023-00256330, and RS-2023-00259686), which was awarded to B-D S. Y-J B. and S-J L. was partly supported by the Ministry of the Interior and Security as a human resources development project in the field of disaster management.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Baek, YJ., Lim, SJ. & So, BD. Development and Benchmarking of the Shallow Shelf Approximation Ice Sheet Dynamics Module. Ocean Sci. J. 58, 27 (2023). https://doi.org/10.1007/s12601-023-00120-3
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s12601-023-00120-3