Brief communication: Time step dependence (and fixes) in Stokes simulations of calving ice shelves
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) dependen...
| Published in: | The Cryosphere |
|---|---|
| Main Authors: | , |
| Format: | Article in Journal/Newspaper |
| Language: | English |
| Published: |
Copernicus Publications
2020
|
| Subjects: | |
| Online Access: | https://doi.org/10.5194/tc-14-3209-2020 https://doaj.org/article/baacce07730e4b3e9b01e8528b53f892 |
| Summary: | 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. |
|---|