Shallow shelf approximation as a “sliding law” in a thermomechanically coupled ice sheet model

Abstract

[1] The shallow shelf approximation, a balance of membrane stresses for ice flow, is an effective ''sliding law'' for ice sheet modeling. Our use of it as a sliding law becomes a standard model for ice stream flow when the sliding velocity is large (100 m a À1 and faster). Following Schoof (2006a), we describe the basal resistance as plastic till for which the yield stress is given by a Mohr-Coulomb formula. Pore water pressure is related to basal melt rate. The velocity field used in the mass continuity and conservation of energy equations is an average of velocities from the shallow shelf approximation and the nonsliding shallow ice approximation. Using this scheme, our model has realistic, time-dependent ice streams which exhibit the range of surface velocities seen in actual ice streams. We demonstrate the model at high spatial resolution (5 km grid) over multiple millenia using its implementation in the Parallel Ice Sheet Model. Numerical experiments show that the entire scheme is stable with respect to many parameter changes. Some experiments reveal significant ice stream variability in a hypothetical steady climate, with characteristic cycles on the order of 1000 years. We believe this is the first practical whole ice sheet model with a unified treatment of vertical shear stresses and membrane stresses. It is capable of high-resolution, thermomechanically coupled, multimillenia simulations of ice sheets containing ice streams.Citation: Bueler, E., and J. Brown (2009), Shallow shelf approximation as a ''sliding law'' in a thermomechanically coupled ice sheet model,

Understanding and Modelling Rapid Dynamic Changes of Tidewater Outlet Glaciers: Issues and Implications

“…Such models are becoming readily available now (Pattyn et al 2008) but are computationally far more expensive and therefore unlikely to be used for longer-term simulations on an full ice sheet scale in the near future. Thus, it will be crucial to rigorously test and assess the validity of different approximations of flow physics such as the shallow-shelf approximation (MacAyeal (1989); Bueler and Brown (2009); and as used here), a depth-integrated hybrid approximation (Bassis 2010;Schoof and Hindmarsh 2010;Goldberg 2011) or the first-order Blatter/Pattyn type approximation (Blatter 1995;Pattyn 2003) against the full-Stokes case for a marine terminating boundary and in view of computational expense. Importantly, this requires also validation against real world data.…”

A fast Newton–Krylov solver for seasonally varying global ocean biogeochemistry models

“…These models rely on different iceflow approximations, parameterizations of physical processes and initialization procedures. Most ice-flow models rely on simplified first-order approximations of the momentum balance equations, such as the shallow-ice approximation (SIA; Hutter, 1982), the shelfy-stream approximation (SSA;MacAyeal, 1989;Greve and others, 2011;Aschwanden and others, 2012a) or a combination of SIA and SSA (Bueler and Brown, 2009;Pollard and DeConto, 2009). Initialization of these numerical models generally consists of running paleoclimate spin-ups over at least the last glacial cycle, in order to obtain a suitable present-day configuration (Greve and others, 2011;Aschwanden and others, 2012a).…”

On the reconstruction of palaeo-ice sheets: Recent advances and future challenges

“…Earlier recommendations to employ second order SIA and SSA models had to be rectified (Kirchner et al, 2011) and have led to the development of new adaptive, error-based couplings schemes (Ahlkrona et al, submitted), while so-called 'hybrid' models (Pollard andDeConto, 2012, Bueler andBrown, 2009) make use of an heuristic approach of coupling an ice sheet and ice shelf. Hybrid models can also be used for simulations of ice sheet-shelf systems, including grounding line migration, at comparatively low computational cost (e.g.…”