Glacier-bed characteristics that are poorly known and modeled are important in projected sea-level rise from ice-sheet changes under strong warming, especially in the Thwaites Glacier drainage of West Antarctica. Ocean warming may induce ice-shelf thinning or loss, or thinning of ice in estuarine zones, reducing backstress on grounded ice. Models indicate that, in response, more-nearly-plastic beds favor faster ice loss by causing larger flow acceleration, but more-nearly-viscous beds favor localized near-coastal thinning that could speed grounding-zone retreat into interior basins where marine-ice-sheet instability or cliff instability could develop and cause very rapid ice loss. Interpretation of available data indicates that the bed is spatially mosaicked, with both viscous and plastic regions. Flow against bedrock topography removes plastic lubricating tills, exposing bedrock that is eroded on up-glacier sides of obstacles to form moats with exposed bedrock tails extending downglacier adjacent to lee-side soft-till bedforms. Flow against topography also generates high-ice-pressure zones that prevent inflow of lubricating water over distances that scale with the obstacle size. Extending existing observations to sufficiently large regions, and developing models assimilating such data at the appropriate scale, present large, important research challenges that must be met to reliably project future forced sea-level rise.

Sea-level projections depend sensitively on the parameterization used for basal slip in glacier flow models. During slip over rock-beds, ice-bed separation increases with slip velocity and basal water pressure. We present a method for using these variables and measured bed topography to estimate the average bed slope in contact with ice, ${\bar m}$. Three-dimensional numerical modeling of slip over small areas of former beds has shown that changes in ${\bar m}$ with increasing slip velocity and water pressure mimic changes in basal drag. Computed values of ${\bar m}$ can thus provide the form of the slip law that relates drag to velocity and water pressure, avoiding computationally expensive numerical modeling. The method is applied to 618 sections from four former glacier beds. Results generally show an increase in ${\bar m}$, and hence inferred basal drag, with slip velocity up to a limiting value, consistent with a regularized Coulomb slip law.

Glacier motion responds dynamically to changing meltwater inputs, but the multi-decadal response of basal sliding to climate remains poorly constrained due to its sensitivity across multiple timescales. Observational records of glacier motion provide critical benchmarks to decode processes influencing glacier dynamics, but multi-decadal records that precede satellite observation and modern warming are rare. Here we present a record of motion in the ablation zone of Saskatchewan Glacier that spans seven decades. We combine in situ and remote-sensing observations to inform a first-order glacier flow model used to estimate the relative contributions of sliding and internal deformation on dynamics. We find a significant increase in basal sliding rates between melt-seasons in the 1950s and those in the 1990s and 2010s and explore three process-based explanations for this anomalous behavior: (i) the glacier surface steepened over seven decades, maintaining flow-driving stresses despite sustained thinning; (ii) the formation of a proglacial lake after 1955 may support elevated basal water pressures; and (iii) subglacial topography may cause dynamic responses specific to Saskatchewan Glacier. Although further constraints are necessary to ascertain which processes are of greatest importance for Saskatchewan Glacier's dynamic evolution, this record provides a benchmark for studies of multi-decadal glacier dynamics.

Glacier slip is usually described using steady-state sliding laws that relate drag, slip velocity and effective pressure, but where subglacial conditions vary rapidly transient effects may influence slip dynamics. Here we use results from a set of laboratory experiments to examine the transient response of glacier slip over a hard bed to velocity perturbations. The drag and cavity evolution from lab experiments are used to parameterize a rate-and-state drag model that is applied to observations of surface velocity and ice-bed separation from the Greenland ice sheet. The drag model successfully predicts observed lags between changes in ice-bed separation and sliding speed. These lags result from the time (or displacement) required for cavities to evolve from one steady-state condition to another. In comparing drag estimates resulting from applying rate-and-state and steady-state slip laws to transient data, we find the peaks in drag are out of phase. This suggests that in locations where subglacial conditions vary on timescales shorter than those needed for cavity adjustment transient slip processes control basal drag.

Drumlins form at the ice/bed interface through subglacial processes that are not directly observable. The internal stratigraphy of drumlins provides insight into how they developed and associated subglacial processes, but traditional stratigraphic logging techniques are limited to natural exposures and excavations. Using ground-penetrating radar, we imaged the internal stratigraphy of seven drumlins from a recently exposed drumlin field in the forefield of Múlajökull, Iceland. Data were collected with 100 and 200 MHz antennas with maximum resolvable depths of 8 and 4 m, respectively. Longitudinal echograms contained coherent down-ice dipping reflectors over the lengths of the drumlins. Near the drumlin heads (i.e., stoss sides), down-glacier dipping beds lie at high angles to the surface, whereas on the lee sides, the down-glacier dipping beds lie at low angles, or conform, to drumlin surfaces. Transverse echograms exhibited unconformities along the flanks of drumlin heads and conformable bedding across the lee side widths of the drumlins. These observations were ground-truthed with stratigraphic logs from a subset of drumlins and good agreement was found. The stratigraphic patterns support previous conclusions that drumlins at Múlajökull formed on a deformable bed through both depositional and erosional processes which may alternate between its surge and quiescent phases.

Seismic-reflection surveys of the Isle Royale sub-basin, central Lake Superior, reveal two large end moraines and associated glacial sediments deposited during the last cycle of the Laurentide Ice Sheet in the basin. The Isle Royale moraines directly overlie bedrock and are cored with dense, acoustically massive till intercalated down-ice with acoustically stratified outwash. Till and outwash are overlain by glacial varves, a lower red unit and an upper gray unit.

The maximum extent of late Younger Dryas-age readvance into the western Lake Superior basin is uncertain, but it was probably controlled by both ice dynamics and climate. Our data indicate that during retreat from the maximum, the ice paused just long enough to construct the outer of the two moraines, >100 m high, and then retreated to the inner moraine, during which time most of the lower glacial-lacustrine sequence (red varves) was deposited. Retreat from the inner moraine coincided with a marked flux of icebergs at the calving margin and a change to gray varves. Rapid retreat may be related to both an influx of meltwater from Glacial Lake Agassiz about 10,500 cal yr BP and retreat of the calving margin down an adverse slope into the Isle Royale sub-basin.

Seismic studies of glaciers yield insights into spatio-temporal processes within and beneath glaciers on scales relevant to flow and deformation of the ice. These methods enable direct monitoring of the bed in ways that complement other geophysical techniques, such as geodetic or ground penetrating radar observations. In this work, we report on the analysis of passive seismic data collected from the interior of the North East Greenland Ice Stream, the Greenland ice sheet's largest outlet glacier. We record thousands of basal earthquakes, many of which repeat with nearly identical waveforms. We also record many long-duration glacial tremor episodes that migrate across the seismic network with slow velocities (e.g. ~4–12 m s−1). Analysis of the basal earthquakes indicates a transition between times of individual event activity and times of tremor activity. We suggest that both processes are produced by shear slip at localized asperities along the bed. The transition between discrete and quasi-continuous slipping modes may be driven by pore-water pressure transients or heterogeneous strain accumulation in the ice due to strength contrasts of the underlying till.

Theory and experiments indicate that ice–bed separation during glacier slip over 2-D hard beds causes basal shear stress to reach a maximum at a particular slip velocity and decrease at higher velocities. We use the sliding theory of Lliboutry (1968) to explore how friction between debris particles in sliding ice and a rock bed affects this relationship between shear stress and slip velocity. Particle–bed contact forces and associated debris friction increase with increasing slip velocity, owing to increased rates of ice convergence with up-glacier facing surfaces. However, debris friction on diminished areas of the bed counteracts this effect as cavities grow. Thus, friction from debris alone increases only slightly with slip velocity, and for sediment particles larger than ~60 mm in diameter, debris friction peaks and decreases with increasing slip velocity. The effect on the sliding relationship is to steepen its rising limb and shift its shear stress peak to a slightly higher velocity. These results, which exclude the effect of debris friction on cavity size and debris concentrations above ~15%, indicate that the effect of debris in ice is to increase basal shear stress but not significantly change the form of the sliding relationship.

Subglacial rock friction is an important control on the sliding dynamics and erosive potential of hard-bedded glaciers, yet it remains largely unconstrained. To explore the relative influence of basal melt rate, effective stress and ice temperature on frictional resistance, we conducted abrasion experiments in which limestone beds were slid beneath a fixed slab of ice laden with granitic rock fragments. Shear stress scales linearly with melt rate and cryostatic stress, confirming that both viscous drag and effective stress are first-order controls on the contact force in drained conditions. Furthermore, temperature gradients in the ice increase the contribution of viscous drag on basal shear stress. In all experiments, the relationship between melt rate and shear stress is best explained by a model that accounts for the effects of regelation and viscous creep on the bed-normal drag force. We interpret this to mean fluid flow around entrained clasts contributed to basal drag even at subfreezing temperatures. Incorporating premelting dynamics into the Watts/Hallet model for subglacial rock friction, we find that the predicted debris-bed drag decreases by approximately an order of magnitude, with a corresponding ~3.5 × increase in the transition radius. This is lower than we observe for ice slightly below the pressure melting point.

The evolution of glaciers and ice sheets depends on processes in the subglacial environment. Shear seismicity along the ice–bed interface provides a window into these processes. Such seismicity requires a rapid loss of strength that is typically ascribed to rate-weakening friction, i.e., decreasing friction with sliding or sliding rate. Many friction experiments have investigated glacial materials at the temperate conditions typical of fast flowing glacier beds. To our knowledge, however, these studies have all found rate-strengthening friction. Here, we investigate the possibility that rate-weakening rock-on-rock friction between sediments frozen to the bottom of the glacier and the underlying water-saturated sediments or bedrock may be responsible for subglacial shear seismicity along temperate glacier beds. We test this ‘entrainment-seismicity hypothesis’ using targeted laboratory experiments and simple models of glacier sliding, seismicity and sediment entrainment. These models suggest that sediment entrainment may be a necessary but not sufficient condition for the occurrence of basal shear seismicity. We propose that stagnation at the Whillans Ice Stream, West Antarctica may be caused by the growth of a frozen fringe of entrained sediment in the ice stream margins. Our results suggest that basal shear seismicity may indicate geomorphic activity.

The contribution of glaciers to sea-level rise and their effects on landscape evolution depend on the poorly known relationship between sliding speed and drag at the ice/bed interface. Results from experiments with a new rotary laboratory device demonstrate empirically for the first time a double-valued drag relationship like that suggested by some sliding theories: steady drag on a rigid, sinusoidal bed increases, peaks and declines at progressively higher sliding speeds due to growth of cavities in the lee sides of bed undulations. Drag decreases with increased sliding speed if cavities extend beyond the inflection points of up-glacier facing surfaces, so that adverse bed slopes in contact with ice diminish with further cavity growth. These results indicate that shear tractions on glacier beds can potentially decrease due to increases in sliding speed driven by weather or climate variability, promoting even more rapid glacier motion by requiring greater strain rates to produce resistive stresses. Although a double-valued drag relationship has not yet been demonstrated for the complicated geometries of real glacier beds, both its potential major implications and the characteristically convex stoss surfaces of bumps on real glacier beds provide stimulus for exploring the effects of this relationship in ice-sheet models.