Ice Till Research Articles

An experimental baseline for ice-till strain indicators

Subglacial till can deform when overriding ice exerts shear traction at the ice–till interface. This deformation leaves a strain signature in the till, aligning grains in the direction of ice flow and producing a range of diagnostic microstructures. Constraining the conditions that produce these kinematic indicators is key to interpreting the myriad of features found in basal till deposits. Here, we used a cryogenic ring shear device with transparent sample chamber walls to slip a ring of temperate ice over a till bed from which we examined the strain signature in the till. We used cameras mounted to the side of the ring shear and bead strings inserted in the till to estimate the strain distribution within the till layer. Following the completion of the experiment, we extracted and analyzed anisotropy of magnetic susceptibility (AMS) samples and created thin sections of the till bed for microstructure analysis. We then compared the AMS and microstructures with the observed strain history to examine the relationship between kinematic indicators and strain in a setting where shear traction is supplied by ice. We found that AMS fabrics show a high degree of clustering in regions of high strain near the ice–till interface. In the uppermost zone of till, k1 eigenvector azimuths are generally aligned with ice flow, and S1 eigenvalues are high. However, S1 eigenvalues and the alignment of the k1 eigenvector with ice flow decrease nonlinearly with distance from the ice–till interface. There is a high occurrence of microshears in the zone of increased deformation.

Shear dilation of subglacial till results in time-dependent sliding laws

The dynamics of glacial sliding over water-saturated tills are poorly constrained and difficult to capture realistically in large-scale models. Experiments characterize till as a plastic material with a pressure-dependent yield stress, but the subglacial water pressure may fluctuate on annual to daily timescales, leading to transient adjustment of the till. We construct a continuum two-phase model of coupled fluid and solid deformation, describing the movement of water through the pore space of a till that is itself dilating and deforming. By forcing the model with time-dependent effective pressure at the ice–till interface, we infer the resulting relationships between basal traction, solid fraction and rate of deformation. We find that shear dilation introduces internal pressure variations and transient dilatant strengthening emerges, leading to hysteretic behaviour in low-permeability materials. The result is a time-dependent effective sliding law, with permeability-dependent lag between changes in effective pressure and the slidingspeed. This deviation from traditional steady-state sliding laws may play an important role in a wide range of transient ice-sheet phenomena, from glacier surges to the tidal response of ice streams.

Centennial-scale variability of sea-ice cover in the Chukchi Sea since AD 1850 based on biomarker reconstruction

Paleo-climate proxy records documenting sea-ice extent are important sources of information to assess the time of emergence and magnitude of on-going changes in the Arctic Ocean and better predict future climate and environmental evolution in that region. In this study, a suite of geochemical tracers including total organic carbon (TOC), total nitrogen (TN), carbon/nitrogen ratio (C/N), stable isotope composition of organic carbon and nitrogen (δ 13C, δ 15N), and phytoplankton biomarkers (highly branched isoprenoids (HBIs) and sterols) were measured in a marine sediment core to document the sea-ice variability in the Chukchi Sea since the beginning of the Industrial Era. The downcore profile of the sea-ice proxy HBIs suggests a transition from extensive sea ice in the late 19th century to Marginal Ice Zone (MIZ) in AD 1930–1990s and then moderate sea-ice cover since 1990s. Rising of all HBI abundances between AD 1865–1875 indicate a transient retreat of summer ice edge off the shelf followed by a return to near-perennial sea ice till 1920–1930 as revealed by the absence of HBIs and brassicasterol. Sea ice retreat occurred again in AD 1920–1930 and followed by colder decades in 1940s–1960s before a sustained decline since the 1990s. The downcore profile of C/N, δ 13C of organic matter and sterols indicates a gradual increase of terrigenous inputs accelerating during the most recent decades likely due to enhanced fluvial run-off and sediment transport by sediment-laden sea ice. Concomitantly, increasing δ 15N values suggest limited nutrient utilization due to enhanced stratification of the surface ocean caused by increased freshening. The role of the Arctic oscillation (AO), the Pacific decadal oscillation (PDO) and Atlantic Multidecadal Oscillation (AMO) are discussed to explore potential drivers of the observed sea-ice changes.

Experimental Evidence of Velocity-Weakening Friction during Ice Slip over Frozen Till: Implications for Basal Seismicity in Fast Moving, Soft-Bed Glaciers and Ice Streams

Abstract Observations of glacier slip over till beds, across a range of spatial and temporal scales, show abundant seismicity ranging from Mw∼−2 microearthquakes and tremor (submeter asperities and millisecond duration) to Mw∼7 slow-slip events (∼50 km rupture lengths and ∼30 min durations). A complete understanding of the mechanisms capable of producing seismic signals in these environments represents a strong constraint on bed conditions. In particular, there is a lack of experimental confirmation of velocity-weakening behavior of ice slipping on till, where friction decreases with increasing velocity—a necessity for nucleating seismic slip. To measure the frictional strength and stability of ice sliding against till, we performed a series of double-direct-shear experiments at controlled temperatures slightly above and below the ice melting point. Our results confirm velocity-strengthening ice–till slip at melting temperatures, as has been found in the few previous studies. We provide best-fit rate-and-state friction parameters and their standard deviations from averaging 13 experiments at equivalent conditions. We find evidence of similar velocity-strengthening behavior with 50% by volume debris-laden ice slid against till under the same conditions. In contrast, velocity-weakening and linear time-dependent healing of ice–till slip is present at temperatures slightly below the melting point, providing an experimentally supported mechanism for subglacial seismicity on soft-beds. The stability parameter (a−b) decreases with slip velocity, and evolution occurs over large (mm scale) displacements, suggesting that shear heating and melt buildup is responsible for the weakening. These measurements provide insight into subglacial stiffness in which seismicity of this type might be expected. We discuss glaciological circumstances pointing to potential field targets in which to test this frozen seismic asperity hypothesis.

The influence of a poroelastic till on rapid subglacial flooding and cavity formation

We develop a model of the rapid propagation of water at the contact between elastic glacial ice and a poroelastic subglacial till, motivated by observations of the rapid drainage of supraglacial lakes in Greenland. By treating the ice as an elastic bending beam, the fluid dynamics of contact with the subglacial hydrological network, which is modelled as a saturated poroelastic till, can be examined in detail. The model describes the formation and dynamics of an axisymmetric subglacial cavity, and the spread of pore pressure, in response to injection of fluid. A combination of numerical simulation and asymptotic analysis is used to describe these dynamics for both a rigid and a deformable porous till, and for both laminar and turbulent fluid flow. For constant injection rates and laminar flow, the cavity is isostatic and its spread is controlled by bending of the ice and suction of pore water in the vicinity of the ice–till contact. For a deformable till, this control can be modified: generically, a flexural wave that is initially trapped in advance of the contact point relaxes over time by diffusion of pore pressure ahead of the cavity. While the dynamics are found to be relatively insensitive to the properties of the subglacial till during injection with a constant flux, significant dependence on the till properties is manifest during the subsequent spread of a constant volume. A simple hybrid turbulent–laminar model is presented to account for fast injection rates of water: in this case, self-similar turbulent propagation can initially control the spread of the cavity, but there is a transition to laminar control in the vicinity of the ice–till contact point as the flow slows. Finally, the model results are compared with recent geophysical observations of the rapid drainage of supraglacial lakes in Greenland; the comparison provides qualitative agreement and raises suggestions for future quantitative comparison.

Formation of the snow isotopic composition on the Arctic glaciers (case study of Spitsbergen and Severnaya Zemlya)

Stable isotope studies of snow cover along several itineraries through Spitsbergen were carried out (from Gronfjorden coast till ice divide of Gronfjord – Fridtjov glacier and from Barentsburg settlement till Olav mountain). Results obtained were compared with ones fixed for snow on Severnaya Zemlya archipelago. Uncertainty of condensation processes were shown for precipitation falling on Arctic polar ice domes (sometimes we fixed altitudinal isotope gradient and in other cases the gradient was absent). In periglacial zone of arctic islands isotope values are scattered randomly, covering the entire annual precipitation range due to strong snow-storms.

Importance of basal processes in simulations of a surging Svalbard outlet glacier

Abstract. The outlet glacier of Basin 3 (B3) of Austfonna ice cap, Svalbard, is one of the fastest outlet glaciers in Svalbard, and shows dramatic changes since 1995. In addition to previously observed seasonal summer speed-up associated with the melt season, the winter speed of B3 has accelerated approximately fivefold since 1995. We use the Elmer/Ice full-Stokes model for ice dynamics to infer spatial distributions of basal drag for the winter seasons of 1995, 2008 and 2011. This "inverse" method is based on minimising discrepancy between modelled and observed surface velocities, using satellite remotely sensed velocity fields. We generate steady-state temperature distributions for 1995 and 2011. Frictional heating caused by basal sliding contributes significantly to basal temperatures of the B3 outlet glacier, with heat advection (a longer-timescale process than frictional heating) also being important in the steady state. We present a sensitivity experiment consisting of transient simulations under present-day forcing to demonstrate that using a temporally fixed basal drag field obtained through inversion can lead to thickness change errors of the order of 2 m year−1. Hence it is essential to incorporate the evolution of basal processes in future projections of the evolution of B3. Informed by a combination of our inverse method results and previous studies, we hypothesise a system of processes and feedbacks involving till deformation and basal hydrology to explain both the seasonal accelerations (short residence time pooling of meltwater at the ice–till interface) and the ongoing interannual speed-up (gradual penetration of water into the till, reducing till strength).

An assessment of fluting and “till esker” formation on the foreland of Sandfellsjökull, Iceland

Intensive sedimentological analysis of tills, parallel-sided flutings, and a conduit fill (till esker) on the recently deglaciated foreland of Sandfellsjökull (a temperate, actively receding piedmont lobe draining the eastern margin of the Mýrdalsjökull ice cap in southern Iceland) reveals that subglacial sediment was emplaced predominantly as a dilatant A horizon after an initial phase of glacitectonite production from local outwash. Clast macrofabric strengths weakened through time, presumably in response to increasingly wetter substrate conditions, until the emplacement of the fluted till when clast alignment was more tightly constrained. The fluting phase was initiated only after boulders were plucked from an adjacent bedrock cliff and lodged at the ice–till interface. This suggests that the parallel-sided flutings constrained shearing at the ice–bed interface in a similar fashion to a fault gouge, thereby strengthening clast alignments but allowing the till to mould around larger obstacles to produce herringbone fabric patterns. Strong clast macrofabric strengths in the till esker are not parallel to the ice flow direction, but instead align parallel with the local axis of the meandering ridge. This suggests that the till esker was likely produced by the squeezing of a dilatant till into an elongated cavity or R-channel after meltwater evacuation, followed by localized deformation, constrained by the channel walls, as the channel moved downvalley because of ice flow. The short period over which the deformation took place was enough to impart weaker fabrics in parts of the conduit that lie transverse to the ice flow direction and stronger fabrics where the conduit trends toward parallelism with the ice flow.

Using multi-beam RADARSAT-1 imagery to augment mapping surficial geology in northwest Alberta, Canada

This paper evaluates the applicability of multi-beam RADARSAT-1 satellite imagery for mapping surficial sediments in northwest Alberta, Canada, using simple arithmetic algorithms. The goal is to make the RADARSAT-1 imagery acquired by the Alberta Geological Survey (AGS) usable to common users who have no access to complex image classification software and (or) expertise. Unlike the Canadian Shield, northern Alberta is mantled with much thicker glacial deposits that are covered by the Northern Boreal Forest. This limits the use of conventional image classification techniques, as there is not a one-to-one relationship between vegetation and the underlying surficial sediments. The RADARSAT-1 data utilized consists of RADARSAT-1 imagery acquired using Standard Beam 1 and 7 in both ascending and descending passes. Moisture and surface roughness information can be extracted by applying algorithms that calculate the differences between images acquired from different incidence angles and look directions (ascending versus descending looks). The image created using different incidence angles highlights flat surfaces, including wetlands, agricultural areas, and forest fire scars. The image generated using opposing look directions enhances geomorphologic features, including ice flow features, meltwater channels, till hummocks, and sand dune ridges. In addition, the enhanced RADARSAT-1 imagery (includes integration of topographic information provided by a digital elevation model (DEM)) provides an additional level of confidence to the air photo interpretation conducted in regions with little or no ground access and also provides a regional perspective on large-scale subtle surface morphological elements (e.g., ice flow patterns) that might be overlooked in individual aerial photographs.

Modeling the subglacial hydrology of the James Lobe of the Laurentide Ice Sheet

To investigate the drainage conditions that might be expected to develop beneath soft-bedded ice sheets, we modeled the subglacial hydrology of the James Lobe of the Laurentide Ice Sheet from Hudson Bay to the Missouri River. Simulations suggest the James Lobe had little effect on regional groundwater flow because the poorly conductive Upper-Cretaceous shale that occupies the upper layer of the bedrock would have functioned as a regional aquitard. This implies that general northward groundwater flow out of the Williston Basin has likely persisted throughout the Quaternary. Moreover, the simulations indicate that the regional aquifer system could not have drained even the minimum amount of basal meltwater that might have been produced from at the glacier bed. Therefore, excess drainage must have occurred by some sort of channelized drainage network at the ice–till interface. Using a regional groundwater model to determine the hydraulic conductivity for an equivalent porous medium in a 1-m thick zone between the ice and underlying sediment, and assuming conduit dimensions from previous theoretical work, we use a theoretical karst aquifer analog as a heuristic approach to estimate the spacing of subglacial conduits that would have been required at the ice–till interface to evacuate the minimum water flux. Results suggest that for conduits assumed to be on the order of a tenth of a meter deep and up to a meter wide, inter-conduit spacing must be on the order of tens–hundreds of meters apart to maintain basal water pressures below the ice overburden pressure while evacuating the hypothesized minimum meltwater flux.

In situ measurements of till deformation and water pressure

AbstractA newly developed hammer was used to insert two autonomous probes 0.8 m and 2.1 m into clast-rich subglacial till under Black Rapids Glacier, Alaska, USA. Both probes were instrumented with a dual-axis tilt sensor and a pore-water pressure transducer. The data are compared to a 75 day record of surface velocities. Till deformation at depth was found to be highly seasonal: it is significant during an early-season speed-up event, but during long periods thereafter measured till deformation rates are negligible. Both tilt records show rotation around the probe axis, which indicates a change in tilt direction of about 30°. The tilt records are very similar, suggesting spatial homogeneity on the scale of the probe separation (4 m horizontal and 3.3 m vertical). There is evidence that during much of the year sliding of ice over till or deformation of a thin till layer (<20 cm) accounts for at least two-thirds of total basal motion. Basal motion accounts for 50–70% of the total surface motion. The inferred amount of ice–till sliding is larger than that found at the same location in a previous study, when surface velocities were about 10% lower. We suggest that variations in ice–till coupling account for the observed variations in mean annual speed.

Probing the till beneath Black Rapids Glacier, Alaska, USA

AbstractA heavy down-hole hammer actuated from the surface by a light composition rope was used to place instrumented probes into the active, 7m thick, clast-rich till underlying a site on Black Rapids Glacier, Alaska, USA, where the ice is 500m thick. A till penetration of about 2.5m was obtained, and greater depths seem possible. The probes measured pore-water pressure and two axes of tilt, which they broadcasted, without wires, to a receiver just above the ice–till interface.

Mechanical and hydraulic properties of till beneath Bakaninbreen, Svalbard

AbstractBakaninbreen, Svalbard, was in the late phase of a prolonged surge during 1995–96. In order to investigate the nature of the basal mechanical and hydraulic regime, sets of tilt sensors equipped with pore-pressure sensors were installed into the basal till layer. Water-pressure transducers were also installed in the boreholes above the ice–till interface. During 1995, the measured mean strain rate in the till was 0.44 a−1 up-glacier of the surge front, while down-glacier the rate was 0.04 a−1. Our interpretation indicates that deformation was the dominant component of basal motion. Strain rates fluctuated strongly, and periods of both positive (down-glacier) and negative (up-glacier) strain occurred. No consistent relationship was apparent between effective pressure and strain rate within the basal till layer. The viscosity of the basal sediments ranged from 1.3 × 1012 to 3.4 × 1012 Pa s. Such high values suggest that the till layer beneath Bakaninbreen may be partially frozen. Analysis of simultaneous measurements of pore-water and borehole water pressure suggests the hydraulic diffusivity of the till layer ranges from 1.2 to 5.0 × 10−6 m2 s−1, and the hydraulic conductivity from 3.2 × 10−7 to 7.7 × 10−8 m s−1.

Coupled ice–till deformation near subglacial channels and cavities

AbstractPrevious models of ice–till deformation near subglacial channels or cavities neglect the fact that the motions of the two materials are coupled, and thus the interface between ice and till may not remain stationary. Here, we analyze in succession two models which address the effect of such coupling via specification of appropriate continuity conditions for stress and velocity across the interface. The modelled scenario is that of a shallow channel–cavity, with its long axis parallel to the principal ice-flow direction, overlying actively deforming till sediments. By applying asymptotic techniques, we investigate how the pattern and velocity of the creep flow depend generally on the ratio between the ice and till viscosities, and on the deforming-till thickness. A more sophisticated, non-linear rheology for till sediments is then introduced. It reveals that the two-way interaction between water percolation and deformation in the till will enhance the localization of sediment flow near the channel margins. The length scale over which transition of effective stress in the till takes place — from its relatively high, far-field value to the low, channel value — is found to depend critically on a dimensionless permeability parameter (Λ). In any case, coupled deformation causes sediment (and ice) flow towards the channel, subsidence of the ice–till interface just outside the channel, and extension of the area over which the ice is in contact with till. Apart from having direct implications for subglacial sediment transport, these results indicate that coupled deformation can contribute significantly to the spatial evolution of stress distribution under a glacier, and thus its incorporation into future sliding and drainage theories for a soft bed should be considered essential.

Sliding versus till deformation in the fast motion of an ice stream over a viscous till

AbstractThe partitioning between till deformation and sliding in the fast flow of ice streams with active basal melting is examined assuming no adhesion of the till to the ice base and incompressible viscous fluid behavior for the till. For deforming-till thickness of 10 m or less the predicted contribution to basal motion by sliding is larger than shearing in the till unless there is short-scale roughness with wavelengths less than order 0.1 m on the ice sole. At such short scales strain heating within the till and focused melting on the ice sole would quickly eliminate the roughness. Thus, fast flow over a till bed would be expected to be mostly by sliding over the subglacial till. More realistic continuum behavior of the till including non-linear and compressible deformation strengthens the conclusion. If sliding is not dominant, then there must be adhesion of the till to the ice base, some mechanism that continuously generates short-scale roughness on the ice–till interface, or very weak internal slip boundaries within the till.

Glacier motion dominated by processes deep in underlying till

AbstractBlack Rapids Glacier is a 40 km long surge-type glacier in the central Alaska Range. In spring 1997 a wireline drill rig was set up at a location where the measured surface velocities are high and seasonal and annual velocity variations are large. The drilling revealed a layer of subglacial “till”, up to 7 m thick, that is believed to be water-saturated. At one location a string of instruments, containing three dual-axis tiltmeters and one piezometer, was successfully introduced into the till. The tiltmeters monitored the inclination of the borehole at the ice–till interface and at 1 and 2 m into the till, for 410 days. They showed that no significant deformation occurred in the upper 2 m of the till layer, and no significant amount of the basal motion was due to sliding of the ice over the till. The measured surface velocity at the drill site is about 60 m a−1, of which 20–30 m a can be accounted for by ice deformation. Almost the entire amount of basal motion, 30–40 m a−1, was taken up at a depth of > 2 m in the till, possibly in discrete shear layers, or as sliding of till over the underlying bedrock. We propose that the large-scale mobilization of such till layers is a key factor in initiating glacier surges.

Coupled ice–till dynamics and the seeding of drumlins and bedrock forms

AbstractThe geomorphological effects of ice sliding over till, internal deformation of till and till sliding over bedrock are considered. Two questions are examined: (1) is the till-sheet flow unstable, i.e. is a layer of uniform thickness maintained or not, and (2) does the slip of till over bedrock cause amplification of relief of the bedrock? Such instabilities seem to be necessary to explain such features as drumlins and whaleback forms.It is found that the answer to (1) and (2) depends on the position of the system in a parameter space, defined by the till rheology, and applied shear stress, the effective pressure at the ice—till interface, the thickness of ice and till and the wavelength of the instability. Two configurations are considered: one where the wavelength of the perturbation is much less than the the ice-thickness, which is related to the classical Nye—Kamb solution for flow over bumps; and one where the wavelength is much greater than the ice thickness, where the mechanics are described by the shallow-ice approximation. In both cases, substantial areas of parameter space, where till-sheet and bedrock modes are unstable, are found. The conceptually related Smalley—Unwin bifurcation is re-examined.The physical mechanisms by which ice and till flows couple are examined. At very short wavelengths (~10 m), the ice is so rigid that it forces till waves to move at the ice velocity; while at long wavelengths (~1000 m), the flows become essentially uncoupled and till waves move at the kinematic velocity. At intermediate wavelengths (~100 m), high growth rates occur ; this is postulated to be the scale of drumlin seeding.

Measurement and control of ice adhesion to aluminum 6061 alloy

A new experimental strategy for measuring the tensile strength of ice coatings to structural surfaces is presented. In this experiment, a laser-induced compressive stress pulse travels through a 1 mm-thick substrate disc that has a layer of ice grown on its front surface. The compressive stress pulse reflects into a tensile wave from the free surface of the ice and pulls the ice\\interface apart, given a sufficient amplitude. The interface strength was calculated by recording the free surface velocity of an Al substrate using a Doppler interferometer and calculating the stress at the interface using a finite-difference elastic wave mechanics simulation with the free surface velocity as an input. The test procedure was used to study ice adhesion on 6061 aluminum alloy sheets. It was found that the adhesion strength of ice to unpolished aluminum substrates was 274 MPa at -10°C. This value decreased with temperature, down to 179 MPa at -40°C. Interestingly, this decrement in the tensile strength could be directly related to the existence of a liquid-like layer that is known to exist on the surface of solid ice till -30°C. The interface strength was also shown to decrease by polishing the Al substrate surface or by adding thin polymer coatings on the unpolished Al substrate. The sensitivity of the technique to such microstructural changes in the interfacial region is indicative of the experiments ability to provide basic adhesion data, which in turn, can be used to solve the deicing problem from a fundamental standpoint. 1998 Elsevier Science Ltd. All rights reserved.

Ice-stream surface texture, sticky spots, waves and breathers: the coupled flow of ice, till and water

AbstractThis paper addresses the coupling of flows of ice, till and water, and the issue of whether such coupling provides mechanisms for meso-scale (kilometres to tens of kilometres) variability in ice-sheet flow and texture. The question of whether effective pressures at the ice-bed interface are statically or hydraulically controlled is examined in this paper. The answer is scale dependent, and has a significant effect on the relationship between ice surface and basal topography.The consequences of these considerations on till flow, coupled ice–till flow and coupled ice, till and water flow are examined. An analysis of till-flow kinematics and shock formation is carried out. The linear stability of coupled long-wavelength ice-till flow is analysed, and regions in parameter space where this flow is unstable, with rather small rate constants are found. Upstream-moving ice surface waves are predicted. The linear stablity of coupled ice–till–water flow is examined, where water flow is modelled using a basal flow system with effective-pressure-dependent properties. Again, regions in parameter space where the system is linearly unstable are found, this time with relatively rapid rate constants. The water pressure exhibits “breather” modes.These analyses assume that there is a substantial basal traction. A problem with models of ice streams wholly restrained at the side is identified: they seem to predict erosion rates which are unfeasibly large.There appears to be sufficient variability in the ice–till– water system to potentially explain texture in ice-stream surfaces, variations in ice-stream thickness of tens of metres not directly relatable to topography, and waves moving upstream or downstream. Most importantly, the ice-stream–bed system is shown to exhibit meso-scale variability simply by coupling ice flow according to the shallow-ice approximation, till flow according to the hydrostatic thin-till approximation and water flow according to an effective-pressure-dependent hydraulics.

Ice-stream surface texture, sticky spots, waves and breathers: the coupled flow of ice, till and water

AbstractThis paper addresses the coupling of flows of ice, till and water, and the issue of whether such coupling provides mechanisms for meso-scale (kilometres to tens of kilometres) variability in ice-sheet flow and texture. The question of whether effective pressures at the ice-bed interface are statically or hydraulically controlled is examined in this paper. The answer is scale dependent, and has a significant effect on the relationship between ice surface and basal topography.The consequences of these considerations on till flow, coupled ice–till flow and coupled ice, till and water flow are examined. An analysis of till-flow kinematics and shock formation is carried out. The linear stability of coupled long-wavelength ice-till flow is analysed, and regions in parameter space where this flow is unstable, with rather small rate constants are found. Upstream-moving ice surface waves are predicted. The linear stablity of coupled ice–till–water flow is examined, where water flow is modelled using a basal flow system with effective-pressure-dependent properties. Again, regions in parameter space where the system is linearly unstable are found, this time with relatively rapid rate constants. The water pressure exhibits “breather” modes.These analyses assume that there is a substantial basal traction. A problem with models of ice streams wholly restrained at the side is identified: they seem to predict erosion rates which are unfeasibly large.There appears to be sufficient variability in the ice–till– water system to potentially explain texture in ice-stream surfaces, variations in ice-stream thickness of tens of metres not directly relatable to topography, and waves moving upstream or downstream. Most importantly, the ice-stream–bed system is shown to exhibit meso-scale variability simply by coupling ice flow according to the shallow-ice approximation, till flow according to the hydrostatic thin-till approximation and water flow according to an effective-pressure-dependent hydraulics.