Glacier-erosion rates range across orders of magnitude, and much of this variation cannot be attributed to basal sliding rates. Subglacial till acts as lubricating ‘fault gouge’ or ‘sawdust’, and must be removed for rapid subglacial bedrock erosion. Such erosion occurs especially where and when moulin-fed streams access the bed and are unconstrained by supercooling or other processes. Streams also may directly erode bedrock, likely with strong time-evolution. Erosion is primarily by quarrying, aided by strong fluctuations in the water system driven by variable surface melt and by subglacial earthquakes. Debris-bed friction significantly affects abrasion, quarrying and general glacier flow. Frost heave drives cirque headwall erosion as winter cold air enters bergschrunds, creating temperature gradients to drive water flow along premelted films to growing ice lenses that fracture rock, and the glacier removes the resulting blocks. Recent subglacial bedrock erosion and sediment flux are in many cases much higher than long-term averages. Over glacial cycles, evolution of glacial-valley form feeds back strongly on erosion and deposition. Most of this is poorly quantified, with parts open to argument. Glacial erosion and interactions are important to tectonic and volcanic processes as well as climate and biogeochemical fluxes, motivating vigorous research.

Late Holocene sediment deposits in Pine Island Bay, West Antarctica, are hypothesized to be linked to intensive meltwater drainage during the retreat of the paleo-Pine Island Ice Stream after the Last Glacial Maximum. The uppermost sediment units show an abrupt transition from ice-proximal debris to a draped silt during the late Holocene, which is interpreted to coincide with rapid deglaciation. The small scale and fine sorting of the upper unit could be attributed to origins in subglacial meltwater; however the thickness and deposition rate for this unit imply punctuated- rather than continuous-deposition. This, combined with the deposit's location seaward of large, bedrock basins, has led to the interpretation of this unit as the result of subglacial lake outbursts in these basins. However, the fine-scale sorting of the silt unit is problematic for this energetic interpretation, which should mobilize and deposit a wider range of sediment sizes. To resolve this discrepancy, we present an alternative mechanism in which the silt was sorted by a distributed subglacial water system, stored in bedrock basins far inland of the grounding line, and subsequently eroded at higher flow speeds during retreat. We demonstrate that this mechanism is physically plausible given the subglacial conditions during the late Holocene. We hypothesize that similar silt units observed elsewhere in Antarctica downstream of bedrock basins could be the result of the same mechanism.

Johannes Weertman published his first glaciological paper in 1957 only 5 years after getting his DSc in metallurgy from the Carnegie Institute of Technology. The paper presented the very first sliding law developed quantitatively from first principles, and involved the unconventional idealization of bed roughness using cubic ‘tombstones’ of rock. Since 1957, there has been a great deal of progress in understanding glacier sliding, but few studies, if any, have preserved the original tombstone geometry that was a hallmark of this first theory. The current study presents a partial reanalysis of the sliding process over a bed with tombstone obstacles using modern numerical methods. The result confirms the enduring applicability of Weertman's model as a pedagogical tool and motivates new questions about (1) folding flow near bedrock obstacles that invert normal ice stratigraphy, (2) the presence and role of stress singularities on sharp edges of bedrock, and (3) the validity of a presumption that regelation flow can be plug-like.

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.

Glacial erosion produces vast quantities of fine-grained sediment that has a far-reaching impact on Earth surface processes. To gain a better understanding of the production of glacial silt and clay, we use automated mineralogy to quantify the microstructure and mineralogy of rock and sediment samples from 20 basins in the St. Elias Mountains, Yukon, Canada. Sediments were collected from proglacial streams, while rock samples were collected from ice marginal outcrops and fragmented using electrical pulse disaggregation. For both rock fragments and sediments, we observe a log-normal distribution of grain sizes and a sub-micrometer terminal grain size. We find that the abrasion of silt and clay results in both rounding and the exploitation of through-going fractures. The abundance of inter- versus intragranular fractures depends on mineralogy and size. Unlike the relatively larger grains, where crushing and abrasion are thought to exploit and produce discrete populations of grain sizes, the comminution of fines leads to a grain size, composition and rounding that is continuously distributed across size, and highly dependent on source-rock properties.

Quarrying and abrasion are the two principal processes responsible for glacial erosion of bedrock. The morphologies of glacier hard beds depend on the relative effectiveness of these two processes, as abrasion tends to smooth bedrock surfaces and quarrying tends to roughen them. Here we analyze concentrations of bedrock discontinuities in the Tsanfleuron forefield, Switzerland, to help determine the geologic conditions that favor glacial quarrying over abrasion. Aerial discontinuity concentrations are measured from scaled drone-based photos where fractures and bedding planes in the bedrock are manually mapped. A Tukey honest significant difference test indicates that aerial concentration of bed-normal bedrock discontinuities is not significantly different between quarried and non-quarried areas of the forefield. Thus, an alternative explanation is needed to account for the spatial variability of quarried areas. To investigate the role that bed-parallel discontinuities might play in quarrying, we use a finite element model to simulate bed-normal fracture propagation within a stepped bed with different step heights. Results indicate that higher steps (larger spacing of bed-parallel discontinuities) propagate bed-normal fractures more readily than smaller steps. Thus, the spacing of bed-parallel discontinuities could exert strong control on quarrying by determining the rate that blocks can be loosened from the host rock.

We quantify the magnitude of millennial-scale glacial erosion at Engabreen, a temperate glacier in coastal northern Norway, using the in situ cosmogenic nuclides carbon-14 (14C) and beryllium-10 (10Be) in bedrock exposed recently by glacial retreat. Nuclide concentrations show no dependence on distance down or across the valley. As such, resulting Holocene erosion depths along two transects perpendicular to glacier flow are highly variable with no systematic distribution, ranging from 0.10 to 2.95 m. We observed 14C–10Be ratios elevated above the production ratio in samples of abraded bedrock, which is counter to the expectation for surfaces covered during the Holocene and exposed only recently. Muon reactions produce nuclides at greater depths than do spallation reactions and 14C at production rates at higher than those of 10Be, resulting in 14C–10Be ratios that increase with depth. Therefore, elevated 14C–10Be ratios indicate that sampled sites were deeply plucked during recent cover, the Little Ice Age in this case, and then rapidly abraded prior to retreat. Our results suggest that, while glacial erosion can generate a u-shaped valley profile over long periods of time (e.g., 105–107 years), the discontinuous nature of glacial plucking produces spatially variable patterns of erosion over shorter millennial timescales.

Recent seismic measurements from upper Thwaites Glacier indicate that the bed-type variability is closely related to the along-flow basal topography. In high-relief subglacial highlands, stoss sides of topographic highs have a relatively higher acoustic impedance (‘hard’ bed) with lower acoustic impedance (‘soft’ till) on lee sides. This pattern is similar to observations of many deglaciated terrains. Subglacial hydraulic-potential gradient and its divergence show a tendency for water to diverge over the stoss sides and converge into the lee sides. Convergence favors a thicker or more widespread water system, which can more efficiently decouple ice from the underlying till. Under such circumstances, till deformation does occur but, fluxes are relatively small. Till carried from the lee sides onto stoss sides of downstream bumps should couple to the ice more efficiently, increasing the ability for transport by till deformation. In turn, this suggests that steady-state till transport can be achieved if the stoss-side till layer is thin or discontinuous. In addition, the large basal shear stress generated in the highlands seems too high for a bed lubricated by a continuous although thin deforming till, suggesting till discontinuity, which would allow debris-laden ice to erode bedrock on stoss sides, supplying additional till for transport.

We describe elongate, wet, subglacial bedforms in the shear margins of the NE Greenland Ice Stream and place some constraints on their formation. Lateral shear margin moraines have been observed across the previously glaciated landscape, but little is known about the ice-flow conditions necessary to form these bedforms. Here we describe in situ sediment bedforms under the NE Greenland Ice Stream shear margins that are observed in active-source seismic and ground-penetrating radar surveys. We find bedforms in the shear margins that are ~500 m wide, ~50 m tall, and elongated nearly parallel to ice-flow, including what we believe to be the first subglacial observation of a shear margin moraine. Acoustic impedance analysis of the bedforms shows that they are composed of unconsolidated, deformable, water-saturated till. We use these geophysical observations to place constraints on the possible formation mechanism of these subglacial features.

Late Wisconsinan glacial sediments, exposed on Whidbey Island and Camano Island, Puget Sound (Washington State, USA), were deposited in a proglacial shallow marine/outwash environment during northward retreat of the Puget Lobe of the Cordilleran ice sheet. Sediments mainly comprise massive and cross-bedded sand and gravels, and rhythmically-bedded clay and silt/fine sand couplets, interbedded with diamictons that were deposited by a range of mass flows of different viscosities. Although sediment stratigraphy and ice advance–retreat patterns are well established for the Puget Lobe, brittle and ductile deformation structures within, and separating, these sediment units are less well understood. These structures record the nature of ice–bed interactions taking place in subglacial and proglacial environments. This study examines evidence for these processes and environments. Key deformation structures identified include open to overturned folds, normal and reverse faults, clastic dikes and hydrofractures and passive-loading structures. Evidence for coeval development of ductile and brittle deformation structures shows the close relationship between porewater changes, sediment rheology and sediment system responses to changes in strain caused by ice–bed interactions.

High-resolution digital elevation models of Finland and Sweden based on LiDAR (Light Detection and Ranging) reveal subglacial landforms in great detail. We describe the ice-sheet scale distribution and morphometric characteristics of a glacial landform that is distinctive in morphology and occurs commonly in the central parts of the former Scandinavian Ice Sheet, especially up-ice of the Younger Dryas end moraine zone. We refer to these triangular or V-shaped landforms as murtoos (singular, ‘murtoo’). Murtoos are typically 30–200 m in length and 30–200 m in width with a relief of commonly <5 m. Murtoos have straight and steep edges, a triangular tip oriented parallel to ice-flow direction, and an asymmetric longitudinal profile with a shorter, but steeper down-ice slope. The spatial distribution of murtoos and their geomorphic relation to other landforms indicate that they formed subglacially during times of climate warming and rapid retreat of the Scandinavian Ice Sheet when large amounts of meltwater were delivered to the bed. Murtoos are formed under warm-based ice and may be associated with a non-channelized subglacial hydraulic system that evacuated large discharges of subglacial water.

How landscapes respond to, and evolve from, large jökulhlaups (glacial outburst floods) is poorly constrained due to limited observations and detailed monitoring. We investigate how melt of glacier ice transported and deposited by multiple jökulhlaups during the 2010 eruption of Eyjafjallajökull, Iceland, modified the volume and surface elevation of jökulhlaup deposits. Jökulhlaups generated by the eruption deposited large volumes of sediment and ice, causing significant geomorphic change in the Gígjökull proglacial basin over a 4-week period. Observation of these events enabled robust constraints on the physical properties of the floods which informs our understanding of the deposits. Using ground-based LiDAR, GPS observations and the satellite-image-derived ArcticDEMs, we quantify the post-depositional response of the 60 m-thick Gígjökull sediment package to the meltout of buried ice and other geomorphic processes. Between 2010 and 2016, total deposit volume reduced by −0.95 × 106 m3 a−1, with significant surface lowering of up to 1.88 m a−1. Surface lowering and volumetric loss of the deposits is attributed to three factors: (i) meltout of ice deposited by the jökulhlaups; (ii) rapid melting of the buried Gígjökull glacier snout; and (iii) incision of the proglacial meltwater system into the jökulhlaup deposits.

Ramparted depressions (doughnut-shaped debris-cored ridges with peat- and/or sediment-filled central basins) are commonly perceived to represent the relict collapsed forms of permafrost ground-ice mounds (i.e. pingos or lithalsas). In Wales, UK, ramparted depressions of Late Pleistocene age have been widely attributed to permafrost-related processes. However, a variety of alternative glacial origins for these enigmatic landforms are also consistent with the available geological and geomorphological evidence, although previous studies have barely considered such alternative processes of formation. From detailed geophysical, sedimentological and remote-sensing studies at two field sites, we demonstrate that: (i) the wastage of stagnating glacier ice is a viable alternative explanation for the formation of ramparted depressions in Wales; (ii) the glacial geomorphology and geology of these landforms is analogous to supraglacial and subglacial landforms from the last Laurentide and Fennoscandian ice sheets; (iii) these landforms have significant potential for characterising the nature of deglaciation around the margins of the Irish Sea during the last glacial cycle, and may record evidence for the overextension and stagnation of the south-eastern margin of the Irish Sea Ice Stream; and (iv) investigations of ramparted depressions within formerly glaciated terrains must consider both glacial and periglacial mechanisms of formation.

Glacier thinning and retreat drives initial acceleration of glacier sliding and erosion, de-buttressing of steep valley walls, and destabilization of ice-marginal deposits and bedrock, which can lead to massive rock avalanching and accelerated incision of tributary watersheds. A compelling example of these changes occurred in Taan Fjord in SE Alaska due to the rapid thinning and retreat of Tyndall Glacier over the past half century. Increased glacier sliding speeds led to both increased rates of subglacial erosion and the evacuation of subglacially stored sediments into the proglacial basins. The shrinking glacier also exposed proglacial tributary watersheds to rapid incision and denudation driven by >350 m of baselevel fall in a few decades. Moreover, in October 2015 a large tsunamigenic landslide occurred at the terminus of Tyndall Glacier, largely due to thinning exposing oversteepened, unstable slopes. Sediment yields from the glacier, the landslide and the tributary watersheds, measured from surveys of the sediments in the fjord collected in 1999 and 2016, are compared to ongoing changes in glacier and fjord geometry to investigate the magnitude of glacial and paraglacial denudation in Taan Fjord during retreat. In the last 50 years, sediment yields from the glacier and non-glacial tributaries kept pace with the rapid rate of retreat, and were on par with each other. Notably, basin-averaged erosion rates from the paraglacial landscape were twice that from the glacier, averaging 58 ± 9 and 26 ± 5 mm a−1, respectively. The sharp increases in sediment yields during retreat observed from both the glacier and the adjacent watersheds, including the landslide, highlight the rapid evolution of landscapes undergoing glacier shrinkage.