Spatial variability in bed topography, characterized as bed roughness, impacts ice-sheet flow and organization and can be used to infer subglacial conditions and processes, yet is difficult to quantify due to sparse observations. Paleo-subglacial beds of formerly expanded glaciers found across the Antarctic continental shelf are well preserved, have relatively limited post-glacial sediment cover and contain glacial landforms that can be resolved at sub-meter vertical scales. We analyze high-resolution bathymetry offshore of Pine Island and Thwaites glaciers in the Amundsen Sea to explore spatial variability of bed roughness where streamlined subglacial landforms allow for the determination of ice-flow direction. We quantify bed roughness using std dev. and Fast Fourier Transform methods, each employed at local (100 km) and regional (101–2 km) scales and in along- and across-flow orientations to determine roughness expressions across spatial scales. We find that the magnitude of roughness is impacted by the parameters selected – which are often not sufficiently reported in studies – to quantify roughness. Important spatial patterns can be discerned from high-resolution bathymetry, highlighting both its usefulness in identifying patterns of streaming ice flow and underscores the need for a standardized way of characterizing topographic variability.

Bed topography is a critical parameter for determining the modern-day and future dynamics of ice sheets and their outlet glaciers. This is because the topography controls the state of stress for glaciers. At glacier termini, topography can influence the timing of terminus retreat by controlling access to warm ocean waters and/or by influencing the ability of a glacier terminus to retreat over bed bumps (moraines). Inland from the terminus, the topography can also influence where glacier retreat and thinning can stabilize. In part, this is because of knickpoints in bed topography created through glacial erosion that may influence the extent to which thinning can diffuse inland for an individual glacier and thus, the timing and magnitude of long-term mass loss. Here we provide a review of the current literature on these topics. While much of the reviewed literature assumes that topography is stable on relevant timescales to humans, new research suggests that topography may change much faster than previously thought and this further complicates our ability to project future outlet glacier change.

Glacial ripping involves glaciotectonic disintegration of rock hills and extensive removal of rock at the ice-sheet bed, triggered by hydraulic jacking caused by fluctuating water pressures. Evidence from eastern Sweden shows that glacial ripping caused significant subglacial erosion during the final deglaciation of the Fennoscandian ice sheet, distinct from abrasion and plucking (quarrying). Here we analyse the ice drag forces exerted onto rock obstacles at the base of an ice sheet, and the resisting forces of such rock obstacles: glaciotectonic disintegration requires that ice drag forces exceed the resisting forces of the rock obstacle. We consider rock obstacles of different sizes, shapes and fracture patterns, informed by natural examples from eastern Sweden. Our analysis shows that limited overpressure events, unfavourable fracture patterns, low-transmissivity fractures, slow ice and streamlined rock hamper rock hill disintegration. Conversely, under fast ice flow and fluctuating water pressures, disintegration is possible if the rock hill contains subhorizontal, transmissive fractures. Rock steps on previously smooth, abraded surfaces, caused by hydraulic jacking, also enhance drag forces and can cause disintegration of a rock hill. Glacial ripping is a physically plausible erosion mechanism, under realistic glaciological conditions prevalent near ice margins.

The recently-advancing Taku Glacier is excavating subglacial sediments at high rates over multi-decadal timescales. However, sediment redistribution over shorter timescales remains unquantified. We use a variety of methods to study subglacial and proglacial sediment redistribution on decadal, seasonal, and daily timescales to gain insight into sub- and proglacial landscape formation. Both excavation and deposition were observed from 2003 to 2015 (2.8 ± 0.9 m a−1 to +2.9 ± 0.9 m a−1). The observed patterns imply that a subglacial conduit has occupied the same site over the past decade. Outwash fans on the subaerial end moraine experience fluvial sediment reworking almost year-round, with net sediment gain in winter and net sediment loss in summer, and an overall mass gain between 2005 and 2015. We estimate that tens of meters of sediment still underlie the glacier terminus, sediments which can be remobilized during future activity. However, imminent retreat from the proglacial moraine will limit its sediment supply, leaving the moraine vulnerable to erosion by bordering rivers. Retreat into an over-deepened basin will leave the glacier vulnerable to increased frontal ablation and accelerating retreat.

The Northeast Greenland Ice Stream (NEGIS) is an important dynamic component for the total mass balance of the Greenland ice sheet, as it reaches up to the central divide and drains 12% of the ice sheet. The geometric boundary conditions and in particular the nature of the subglacial bed of the NEGIS are essential to understand its ice flow dynamics. We present a record of more than 8000 km of radar survey lines of multi-channel, ultra-wideband radio echo sounding data covering an area of 24 000 km2, centered on the drill site for the East Greenland Ice-core Project (EGRIP), in the upper part of the NEGIS catchment. Our data yield a new detailed model of ice-thickness distribution and basal topography in the region. The enhanced resolution of our bed topography model shows features which we interpret to be caused by erosional activity, potentially over several glacial–interglacial cycles. Off-nadir reflections from the ice–bed interface in the center of the ice stream indicate a streamlined bed with elongated subglacial landforms. Our new bed topography model will help to improve the basal boundary conditions of NEGIS prescribed for ice flow models and thus foster an improved understanding of the ice-dynamic setting.

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.

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.

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.

Proglacial environments are ideal for studying the development of soils through the changes of rocks exposed by glacier retreat to weathering and microbial processes. Carbon (C) and nitrogen (N) contents as well as soil pH and soil elemental compositions are thought to be dominant factors structuring the bacterial, archaeal and fungal communities in the early stages of soil ecosystem formation. However, the functional linkages between C and N contents, soil composition and microbial community structures remain poorly understood. Here, we describe a multivariate analysis of geochemical properties and associated microbial community structures between a moraine and a glaciofluvial outwash in the proglacial area of a High Arctic glacier (Longyearbreen, Svalbard). Our results reveal distinct differences in developmental stages and heterogeneity between the moraine and the glaciofluvial outwash. We observed significant relationships between C and N contents, δ13Corg and δ15N isotopic ratios, weathering and microbial abundance and community structures. We suggest that the observed differences in microbial and geochemical parameters between the moraine and the glaciofluvial outwash are primarily a result of geomorphological variations of the proglacial terrain.

Surface digital elevation models (DEMs) and slope-related estimates of glacier thickness enable modelling of glacier-bed topographies over large ice-covered areas. Due to the erosive power of glaciers, such bed topographies can contain numerous overdeepenings, which when exposed following glacier retreat may fill with water and form new lakes. In this study, the bed overdeepenings for ~28 000 glaciers (40 775 km2) of the Himalaya-Karakoram region are modelled using GlabTop2 (Glacier Bed Topography model version 2), in which ice thickness is inferred from surface slope by parameterizing basal shear stress as a function of elevation range for each glacier. The modelled ice thicknesses are uncertain (±30%), but spatial patterns of ice thickness and bed elevation primarily depend on surface slopes as derived from the DEM and, hence, are more robust. About 16 000 overdeepenings larger than 104m2 were detected in the modelled glacier beds, covering an area of ~2200 km2 and having a volume of ~120km3 (3-4% of present-day glacier volume). About 5000 of these overdeepenings (1800 km2) have a volume larger than 106m3. The results presented here are useful for anticipating landscape evolution and potential future lake formation with associated opportunities (tourism, hydropower) and risks (lake outbursts).