[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,

Predictions of marine ice-sheet behaviour require models that are able to robustly simulate grounding line migration. We present results of an intercomparison exercise for marine ice-sheet models. Verification is effected by comparison with approximate analytical solutions for flux across the grounding line using simplified geometrical configurations (no lateral variations, no effects of lateral buttressing). Unique steady state grounding line positions exist for ice sheets on a downward sloping bed, while hysteresis occurs across an overdeepened bed, and stable steady state grounding line positions only occur on the downward-sloping sections. Models based on the shallow ice approximation, which does not resolve extensional stresses, do not reproduce the approximate analytical results unless appropriate parameterizations for ice flux are imposed at the grounding line. For extensional-stress resolving "shelfy stream" models, differences between model results were mainly due to the choice of spatial discretization. Moving grid methods were found to be the most accurate at capturing grounding line evolution, since they track the grounding line explicitly. Adaptive mesh refinement can further improve accuracy, including fixed grid models that generally perform poorly at coarse resolution. Fixed grid models, with nested grid representations of the grounding line, are able to generate accurate steady state positions, but can be inaccurate over transients. Only one full-Stokes model was included in the intercomparison, and consequently the accuracy of shelfy stream models as approximations of full-Stokes models remains to be determined in detail, especially during transients

Polythermal conditions are ubiquitous among glaciers, from small valley glaciers to ice sheets. Conventional temperature-based ‘cold-ice’ models of such ice masses cannot account for that portion of the internal energy which is latent heat of liquid water within temperate ice, so such schemes are not energy-conserving when temperate ice is present. Temperature and liquid water fraction are, however, functions of a single enthalpy variable: a small enthalpy change in cold ice is a change in temperature, while a small enthalpy change in temperate ice is a change in liquid water fraction. The unified enthalpy formulation described here models the mass and energy balance for the threedimensional ice fluid, for the surface runoff layer and for the subglacial hydrology layer, together in a single energy-conserving theoretical framework. It is implemented in the Parallel Ice Sheet Model. Results for the Greenland ice sheet are compared with those from a cold-ice scheme. This paper is intended to be an accessible foundation for enthalpy formulations in glaciology.

Abstract. We present the Potsdam Parallel

We present a dynamic equilibrium simulation of the ice sheet-shelf system on Antarctica with the Potsdam Parallel Ice Sheet Model (PISM-PIK). The simulation is initialized with present-day conditions for bed topography and ice thickness and then run to steady state with constant present-day surface mass balance. Surface temperature and sub-shelf basal melt distribution are parameterized. Grounding lines and calving fronts are free to evolve, and their modeled equilibrium state is compared to observational data. A physically-motivated calving law based on horizontal spreading rates allows for realistic calving fronts for various types of shelves. Steady-state dynamics including surface velocity and ice flux are analyzed for whole Antarctica and the Ronne-Filchner and Ross ice shelf areas in particular. The results show that the different flow regimes in sheet and shelves, and the transition zone between them, are captured reasonably well, supporting the approach of superposition of SIA and SSA for the representation of fast motion of grounded ice. This approach also leads to a natural emergence of sliding-dominated flow in stream-like features in this new 3-D marine ice sheet model.