Antarctic ice streams are associated with pressurized subglacial meltwater but the role this water plays in the dynamics of the streams is not known. To address this, we present a model of subglacial water flow below ice sheets, and particularly below ice streams. The base-level flow is fed by subglacial melting and is presumed to take the form of a rough-bedded film, in which the ice is supported by larger clasts, but there is a millimetric water film which submerges the smaller particles. A model for the film is given by two coupled partial differential equations, representing mass conservation of water and ice closure. We assume that there is no sediment transport and solve for water film depth and effective pressure. This is coupled to a vertically integrated, higher order model for ice-sheet dynamics. If there is a sufficiently small amount of meltwater produced (e.g. if ice flux is low), the distributed film and ice sheet are stable, whereas for larger amounts of melt the ice–water system can become unstable, and ice streams form spontaneously as a consequence. We show that this can be explained in terms of a multi-valued sliding law, which arises from a simplified, one-dimensional analysis of the coupled model.

The existence of both water and sediment at the bed of ice streams is well documented, but there is a lack of fundamental understanding about the mechanisms of ice, water and sediment interaction. We pose a model to describe subglacial water flow below ice sheets, in the presence of a deformable sediment layer. Water flows in a rough-bedded film; the ice is supported by larger clasts, but there is a millimetric water layer submerging the smaller particles. Partial differential equations describing the water film are derived from a description of the dynamics of ice, water and mobile sediment. We assume that sediment transport is possible, either as fluvial bedload, but more significantly by ice-driven shearing and by internal squeezing. This provides an instability mechanism for rivulet formation; in the model, downstream sediment transport is compensated by lateral squeezing of till towards the incipient streams. We show that the model predicts the formation of shallow, swamp-like streams, with a typical depth of the order of centimetres. The swamps are stable features, typically with a width of the order of tens to hundreds of metres.

No abstract

ABSTRACT. One challenge in improving our understanding of ice-stream dynamics is to develop models of the spatial and temporal transition from ice-sheet to ice-stream flow. We address this with a new, vertically integrated, higher-order formulation for ice-sheet dynamics that captures the leading-order physics of low aspect ratio, viscous fluid flow, regardless of the amount of slip at the bed. The theory introduces a parameter, , which approximates the ratio of the basal stress to the shear stress scale, providing a measure of the relative importance of sliding and internal deformation. Our model is able to simultaneously describe the dynamics of both a slow-moving sheet and rapidly flowing ice streams. To test the formulation, we apply a triple-valued sliding law as the basal boundary condition and obtain numerical solutions that can be compared with previous work. We investigate the sensitivity of flow regimes and shear margin width to parameter variation, deriving a scaling for the latter. We also consider a double-valued sliding law, which enforces a constant, low basal stress beneath the ice stream. Comparisons of the resultant stress fields illustrate the different stress balances that can maintain ice-stream flow.