Dynamically coupling the non-linear Stokes equations with the shallow ice approximation in glaciology: Description and first applications of the ISCAL method
Number of Authors: 4
2016 (English)In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 308, 1-19 p.Article in journal (Refereed) Published
We propose and implement a new method, called the Ice Sheet Coupled Approximation Levels (ISCAL) method, for simulation of ice sheet flow in large domains during long time-intervals. The method couples the full Stokes (FS) equations with the Shallow Ice Approximation (SIA). The part of the domain where SIA is applied is determined automatically and dynamically based on estimates of the modeling error. For a three dimensional model problem, ISCAL computes the solution substantially faster with a low reduction in accuracy compared to a monolithic FS. Furthermore, ISCAL is shown to be able to detect rapid dynamic changes in the flow. Three different error estimations are applied and compared. Finally, ISCAL is applied to the Greenland Ice Sheet on a quasi-uniform grid, proving ISCAL to be a potential valuable tool for the ice sheet modeling community.
Place, publisher, year, edition, pages
2016. Vol. 308, 1-19 p.
Full Stokes equations, Shallow Ice Approximation, Ice sheet modeling, Error estimation, Finite element method
Physical Geography Computer Science Computational Mathematics
IdentifiersURN: urn:nbn:se:su:diva-127342DOI: 10.1016/j.jcp.2015.12.025ISI: 000369086700001OAI: oai:DiVA.org:su-127342DiVA: diva2:911324