Increasing numerical stability of mountain valley glacier simulations: implementation and testing of free-surface stabilization in Elmer/Ice

This paper concerns a numerical stabilization method for free-surface ice flow called the free-surface stabilization algorithm (FSSA). In the current study, the FSSA is implemented into the numerical ice-flow software Elmer/Ice and tested on synthetic two-dimensional (2D) glaciers, as well as on the real-world glacier of Midtre Lovénbreen, Svalbard. For the synthetic 2D cases it is found that the FSSA method increases the largest stable time-step size at least by a factor of 5 for the case of a gently sloping ice surface (∼ 3°) and by at least a factor of 2 for cases of moderately to steeply inclined surfaces (∼ 6° to 12°) on a fine mesh. Compared with other means of stabilization, the FSSA is the only one in this study that increases largest stable time-step sizes when used alone. Furthermore, the FSSA method increases the overall accuracy for all surface slopes. The largest stable time-step size is found to be smallest for the case of a low sloping surface, despite having overall smaller velocities. For an Arctic-type glacier, Midtre Lovénbreen, the FSSA method doubles the largest stable time-step size; however, the accuracy is in this case slightly lowered in the deeper parts of the glacier, while it increases near edges. The implication is that the non-FSSA method might be more accurate at predicting glacier thinning, while the FSSA method is more suitable for predicting future glacier extent. A possible application of the larger time-step sizes allowed for by the FSSA is for spin-up simulations, where relatively fast-changing climate data can be incorporated on short timescales, while the slow-changing velocity field is updated over larger timescales.</p

Accuracy of the zeroth- and second-order shallow-ice approximation – numerical and theoretical results

In ice sheet modelling, the shallow-ice approximation (SIA) and second-order shallow-ice approximation (SOSIA) schemes are approaches to approximate the solution of the full Stokes equations governing ice sheet dynamics. This is done by writing the solution to the full Stokes equations as an asymptotic expansion in the aspect ratio &epsilon;, i.e. the quotient between a characteristic height and a characteristic length of the ice sheet. SIA retains the zeroth-order terms and SOSIA the zeroth-, first-, and second-order terms in the expansion. Here, we evaluate the order of accuracy of SIA and SOSIA by numerically solving a two-dimensional model problem for different values of &epsilon;, and comparing the solutions with afinite element solution to the full Stokes equations obtained from Elmer/Ice. The SIA and SOSIA solutions are also derived analytically for the model problem. For decreasing &epsilon;, the computed errors in SIA and SOSIA decrease, but not always in the expected way. Moreover, they depend critically on a parameter introduced to avoid singularities in Glen's flow law in the ice model. This is because the assumptions behind the SIA and SOSIA neglect a thick, high-viscosity boundary layer near the ice surface. The sensitivity to the parameter is explained by the analytical solutions. As a verification of the comparison technique, the SIA and SOSIA solutions for a fluid with Newtonian rheology are compared to the solutions by Elmer/Ice, with results agreeing very well with theory

Shallow flows of generalised Newtonian fluids on an inclined plane

We derive a general evolution equation for a shallow layer of a generalised Newtonian fluid undergoing two-dimensional gravity-driven flow on an inclined plane. The flux term appearing in this equation is expressed in terms of an integral involving the prescribed constitutive relation and, crucially, does not require explicit knowledge of the velocity profile of the flow; this allows the equation to be formulated for any generalised Newtonian fluid. In particular, we develop general solutions for travelling waves on a mild slope and for kinematic waves on a moderately steep slope; these results provide simple and accessible models of, for example, the propagation of non-Newtonian mud and debris flows