Brief communication: Time step dependence (and fixes) in Stokes simulations of calving ice shelves

• Published in: The Cryosphere
• Main Authors: B. Berg, J. Bassis
• Format: Article in Journal/Newspaper
• Language: English
• Published: Copernicus Publications 2020
• Subjects: Environmental sciences; GE1-350; Geology; QE1-996.5; Ice Shelves; The Cryosphere
• Online Access: https://doi.org/10.5194/tc-14-3209-2020; https://doaj.org/article/baacce07730e4b3e9b01e8528b53f892

The buoyancy boundary condition applied to floating portions of ice sheets and glaciers in Stokes models requires special consideration when the glacier rapidly departs from hydrostatic equilibrium. This boundary condition can manifest in velocity fields that are unphysically (and strongly) dependent on time step size, thereby contaminating diagnostic stress fields. This can be especially problematic for models of calving glaciers, where rapid changes in geometry cause configurations that suddenly depart from hydrostatic equilibrium and lead to inaccurate estimates of the stress field. Here we show that the unphysical behavior can be cured with minimal computational cost by reintroducing a regularization that corresponds to the acceleration term in the stress balance. This regularization provides consistent velocity solutions for all time step sizes.