Work

Weak form shallow ice approximation models with an improved time-step restriction

Year: 2025

Type: article

Abstract: The shallow ice approximation (SIA) model in strong form is commonly used for inferring the flow dynamics of grounded ice sheets. The solution to the SIA model is a closed-form expression for the velocity field. When that velocity field is used to advance the ice surface in time, the time steps have to take small values due to quadratic scaling in the horizontal mesh size. In this paper, we write the SIA model in weak form and add in the free-surface stabilization algorithm (FSSA) terms. We find numerically that the time-step restriction scaling is improved from quadratic to linear, but only for large horizontal mesh sizes. We then modify the weak form by adding the initially neglected normal stress terms. This allows for a linear time-step restriction across the whole range of horizontal mesh sizes, leading to improve efficiency. Theoretical analysis demonstrates that the inclusion of FSSA stabilization terms transitions the explicit time-stepping treatment of second derivative surface terms to an implicit approach. Moreover, a computational cost analysis, combined with numerical results on stability and accuracy, advocates for preferring the SIA models written in weak form over the standard SIA model.

Source: Journal of Glaciology

DOI: 10.1017/jog.2024.105

Authors: Igor Tominec, Josefin Ahlkrona

Institutions: Stockholm University, Bolin Centre for Climate Research, Swedish e-Science Research Centre

Language: English

Topic: Cryospheric studies and observations

Subfield: Atmospheric Science

Field: Earth and Planetary Sciences

Domain: Physical Sciences

Open Access status: gold

Estimated APC paid: $2,600

Funders: Swedish e-Science Research Centre, Indiana Family and Social Services Administration, Vetenskapsrådet, Stockholms Universitet