Consistent approximations and boundary conditions for ice-sheet dynamics from a principle of least action

Abstract

ABSTRACT. The formulation of a physical problem in terms of a variational (or action) principle conveys significant advantages for the analytical formulation and numerical solution of that problem. One such problem is ice-sheet dynamics as described by non-Newtonian Stokes flow, for which the variational principle can be interpreted as stating that a measure of heat dissipation, due to internal deformation and boundary friction, plus the rate of loss of total potential energy is minimized under the constraint of incompressible flow. By carrying out low-aspect-ratio approximations to the Stokes flow problem within this variational principle, we obtain approximate dynamical equations and boundary conditions that are internally consistent and preserve the analytical structure of the full Stokes system. This also allows us to define an action principle for the popular first-order or 'Blatter-Pattyn' shallow-ice approximation that is distinct from the action principle for the Stokes problem yet preserves its most important properties and elucidates various details about this approximation. Further approximations within this new action functional yield the standard zero-order shallow-ice and shallow-shelf approximations, with their own action principles and boundary conditions. We emphasize the specification of boundary conditions, which are problematic to derive and implement consistently in approximate models but whose formulation is greatly simplified in a variational setting.

Committed sea-level rise for the next century from Greenland ice sheet dynamics during the past decade

“…Our three-dimensional, dynamic ice flow model (25) obeys conservation of linear momentum according to the first-order approximation to the nonlinear Stokes flow equations (24,(32)(33)(34). The model participated in the ice sheet model intercomparison project for higher-order models (35) and model output for tests A-E fall within one standard deviation of the mean defined by all participating models of the same type.…”

<i>Albany/FELIX</i>: a parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis

“…We consider a power-law viscous, incompressible fluid in a low Reynolds number flow, described by the first-order approximation to the nonlinear Stokes flow equations for glaciers and ice sheets (Dukowicz et al, 2010;Schoof and Hindmarsh, 2010). The first-order (FO) approximation, also referred to as the Blatter-Pattyn model (Pattyn, 2003;Blatter, 1995), follows from assumptions of a small geometric aspect ratio, δ = H /L (where H and L are characteristic length scales for the vertical and horizontal dimensions, respectively, and H L), and the assumption that the normal vectors to the ice sheet's upper and lower surfaces, n ∈ R 3 , are nearly vertical:…”

Incorporating modelled subglacial hydrology into inversions for basal drag

“…This formulation is derived from the Stokes equations using variational principles (Dukowicz et al, 2010;Goldberg, 2011) and is a hybrid of the shallow ice approximation (Cuffey and Paterson, 2010) and the shallow shelf approximation (MacAyeal, 1989;Morland, 1987).…”