We present a method for tracking radio-tagged pebbles and cobbles through subglacial meltwater channels under shallow temperate glaciers. Natural particles tagged with active radio transmitters were injected directly into a large subglacial channel 300 m up-glacier from the terminus of the Glacier d'Otemma, Switzerland. A roving antenna was developed to localise tagged particles planimetrically in subglacial and proglacial channel reaches (350 and 150 m long, respectively) using a probabilistic technique, delivering records of the change in particle location and transport distance over time with uncertainty. The roving antenna had a ±5−15 m planimetric precision, a 75% particle localisation rate and operated at a maximum ice depth of 47 m. Additionally, stationary supraglacial and proglacial antennas continuously monitored the passage of tagged particles through consecutive reaches of the channel, constraining the timing of particle transport events. The proglacial antenna system had a 98.1% detection rate and was operational to 0.89 m water depth during testing. Roving and stationary antenna records were combined to create a transport distance model for each particle, which may be used in conjunction with hydraulic data to investigate the kinematics of particle motion. When applied at scale in future studies, this method may be used to reveal the mechanisms and timescales of coarse sediment export from Alpine glaciers.

This study presents the first light detection and ranging (LiDAR)-based morphometric description of a water blister from a past ice-sheet bed caused by rapid supraglacial drainage. The blister formed during the rapid early Holocene deglaciation of the Fennoscandian Ice Sheet (FIS). It is located in southwest Finland within a subglacial meltwater route interpreted to represent the transition from a distributed to a channelized drainage system. A LiDAR digital elevation model was supplemented with sedimentological and ground-penetrating radar data on blister outflow channels and sedimentology of downflow polymorphous mounds and ridges (PMRs). Unlike the water blisters recorded from the rapid drainage of supraglacial lakes on the Greenland Ice Sheet, the smaller blister size here was either due to crevasse or moulin drainage, or was a supraglacial lake drainage that tapped into a pre-existing, relatively efficient drainage system and related semi-sorted sediments, promoting rapid drainage and reworking of PMRs along the meltwater route. The preservation potential or exposure probability of blister marks is presumably low but they can provide important information about evolution of subglacial drainage systems that is of value to modern interpretations of glacial hydrology.

Distributed Acoustic Sensing (DAS) is increasingly recognised as a valuable tool for glaciological seismic applications, although analysing the large data volumes generated in acquisitions poses computational challenges. We show the potential of active-source DAS to image and characterise subglacial sediment beneath a fast-flowing Greenlandic outlet glacier, estimating the thickness of sediment layers to be 20–30 m. However, the lack of subglacial velocity constraint limits the accuracy of this estimate. Constraint could be provided by analysing cryoseismic events in a counterpart 3-day record of passive seismicity through, for example, seismic tomography, but locating them within the 9 TB data volume is computationally inefficient. We describe experiments with data compression using the frequency-wavenumber (f-k) transform ahead of training a convolutional neural network, that provides a ~300-fold improvement in efficiency. In combining active and passive-source and our machine learning framework, the potential of large DAS datasets could be unlocked for a range of future applications.

We use a mathematical model to investigate the effect of basal topography and ice surface slope on transport and deposition of sediment within a water-filled subglacial channel. In our model, three zones of different behaviour occur. In the zone furthest upstream, variations in basal topography lead to sediment deposition under a wide range of conditions. In this first zone, even very small and gradually varying basal undulations (~5 m amplitude) can lead to the deposition of sediment within a modelled channel. Deposition is concentrated on the downstream gradient of subglacial ridges, and on the upstream gradient of subglacial troughs. The thickness and steepness of the ice sheet has a substantial impact on deposition rates, with shallow ice profiles strongly promoting both the magnitude and extent of sediment deposition. In a second zone, all sediment is transported downstream. Finally, a third zone close to the ice margin is characterised by high rates of sediment deposition. The existence of these zones has implications for esker formation and the dynamics of the subglacial environment.

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.

Rapid Access Ice Drill is a new drilling technology capable of quickly accessing the glacial bed of Antarctic ice sheets, retrieving ice core and rock core samples, and providing boreholes for downhole logging of physical properties. Scientific goals include searching for old ice near the glacial bed and sampling subglacial bedrock. During field trials near McMurdo Station on a piedmont glacier at Minna Bluff in the 2019–20 austral summer, we successfully completed a ‘top-to-bottom’ operational sequence in three boreholes by (1) augering through firn, (2) creating a borehole packer seal in non-porous ice, (3) establishing fluid circulation, (4) quickly drilling a borehole in ice at penetration rates up to 1.2 m min−1, (5) acquiring a short ice core at depth, (6) penetrating the glacial bed at a depth of ~677 m, (7) recovering a 3.2 m core of ice, basal till and subglacial bedrock, (8) optically logging the borehole on wireline, (9) testing hydrofracture potential by overpressuring the borehole fluid and (10) operating in an environmentally benign yet rapid field mode. Minna Bluff testing, therefore, demonstrates the effectiveness of this integrated system to drill rapidly through thick ice and penetrate across the glacial bed to take cores of bedrock.

During the 2018/19 Antarctic field season, the British Antarctic Survey (BAS) Basal conditions on Rutford Ice Stream: BEd Access, Monitoring and Ice Sheet History’ (BEAMISH) project drilled three holes through the Rutford Ice Stream, West Antarctica. At up to 2154 m, these are the deepest hot water drilled subglacial access holes yet created, enabling the recovery of sediment from the subglacial environment, and instrumenting the ice stream and its bed. The BEAMISH hot-water drill system was built on extensive experience with the BAS ice shelf hot-water drill and utilises many identical components. With up to 1 MW of heating power available, the hot water drill produces 140 L min−1 of water at 85°C to create a 300 mm diameter access hole to the base of the ice stream. New systems and processes were developed for BEAMISH to aid critical aspects of deep access drilling, most notably the creation of cavities interlinking boreholes at 230 m below the surface and enabling water recirculation throughout the deep drilling operations. The modular design of the BEAMISH drill offers many benefits in its adaptability, redundancy, and minimal logistical footprint. These design features can easily accommodate the modifications needed for future deep, clean access hole creation in the exploration of subglacial environments.

Three holes were drilled to the bed of Rutford Ice Stream, through ice up to 2154 m thick, to investigate the basal processes and conditions associated with fast ice flow and the glacial history of the West Antarctic Ice Sheet. A narrative of the drilling, measuring and sampling activities, as well as some preliminary results and initial interpretations of subglacial conditions, is given. These were the deepest subglacial access holes ever drilled using the hot-water drilling method. Samples of bed and englacial sediments were recovered, and a number of instruments were installed in the ice column and the bed. The ice–bed interface was found to be unfrozen, with an existing, well-developed subglacial hydrological system at high pressure, within ~1% of the ice overburden. The bed itself comprises soft, water-saturated sediments, consistent with previous geophysical interpretations. Englacial sediment quantity varies significantly between two locations ~2 km apart, and possibly over even shorter (~20 m) distances. Difficulties and unusual observations encountered while connecting to the subglacial hydrological system in one hole possibly resulted from the presence of a large clast embedded in the bottom of the ice.

Subglacial sediments have the potential to reveal information about the controls on glacier flow, changes in ice-sheet history and characterise life in those environments. Retrieving sediments from beneath the ice, through hot water drilled access holes at remote field locations, present many challenges. Motivated by the need to minimise weight, corer diameter and simplify assembly and operation, British Antarctic Survey, in collaboration with UWITEC, developed a simple mechanical percussion corer. At depths over 1000 m however, manual operation of the percussion hammer is compromised by the lack of clear operator feedback at the surface. To address this, we present a new auto-release-recovery percussion hammer mechanism that makes coring operations depth independent and improves hammer efficiency. Using a single rope tether for both the corer and hammer operation, this modified percussion corer is relatively simple to operate, easy to maintain, and has successfully operated at a depth of >2130 m.

The Winkie Drill is an agile, commercially available rock coring system. The U.S. Ice Drilling Program has modified a Winkie Drill for subglacial rock and ice/rock interface coring, as well as drilling and coring access holes through ice. The original gasoline engine was replaced with an electric motor though the two-speed gear reducer and Unipress hand feed system were maintained. Using standard aluminum AW34 drill rod (for 33.5 mm diameter core), the system has a depth capability of 120 m. The drill uses forward fluid circulation in a closed loop system. The drilling fluid is Isopar K, selected for favorable properties in polar environment. When firn or snow is present at the drill site, casing with an inflatable packer can be deployed to contain the drill fluid. The Winkie Drill will operate from sea level to high altitudes and operation results in minimal environmental impact. The drill can be easily and quickly assembled and disassembled in the field by two people. All components can be transported by Twin Otter or helicopter to the field site.

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.

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.

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.

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.

Fast ice flow is associated with the deformation of subglacial sediment. Seismic shear velocities, Vs, increase with the rigidity of material and hence can be used to distinguish soft sediment from hard bedrock substrates. Depth profiles of Vs can be obtained from inversions of Rayleigh wave dispersion curves, from passive or active-sources, but these can be highly ambiguous and lack depth sensitivity. Our novel Bayesian transdimensional algorithm, MuLTI, circumvents these issues by adding independent depth constraints to the inversion, also allowing comprehensive uncertainty analysis. We apply MuLTI to the inversion of a Rayleigh wave dataset, acquired using active-source (Multichannel Analysis of Surface Waves) techniques, to characterise sediment distribution beneath the frontal margin of Midtdalsbreen, an outlet of Norway's Hardangerjøkulen ice cap. Ice thickness (0–20 m) is constrained using co-located GPR data. Outputs from MuLTI suggest that partly-frozen sediment (Vs 500–1000 m s−1), overlying bedrock (Vs 2000–2500 m s−1), is present in patches with a thickness of ~4 m, although this approaches the resolvable limit of our Rayleigh wave frequencies (14–100 Hz). Uncertainties immediately beneath the glacier bed are <280 m s−1, implying that MuLTI cannot only distinguish bedrock and sediment substrates but does so with an accuracy sufficient for resolving variations in sediment properties.

Subglacial tills play an important role in glacier dynamics but are difficult to characterize in situ. Amplitude Variation with Angle (AVA) analysis of seismic reflection data can distinguish between stiff tills and deformable tills. However, AVA analysis in mountain glacier environments can be problematic: reflections can be obscured by Rayleigh wave energy scattered from crevasses, and complex basal topography can impede the location of reflection points in 2-D acquisitions. We use a forward model to produce challenging synthetic seismic records in order to test the efficacy of AVA in crevassed and geometrically complex environments. We find that we can distinguish subglacial till types in moderately crevassed environments, where ‘moderate’ depends on crevasse spacing and orientation. The forward model serves as a planning tool, as it can predict AVA success or failure based on characteristics of the study glacier. Applying lessons from the forward model, we perform AVA on a seismic dataset collected from Taku Glacier in Southeast Alaska in March 2016. Taku Glacier is a valley glacier thought to overlay thick sediment deposits. A near-offset polarity reversal confirms that the tills are deformable.

Stable isotopes of water (δ18O and δ2H) were measured in the debris-laden ice underlying an Antarctic blue ice moraine, and in adjoining Law Glacier in the central Transantarctic Mountains. Air bubble content and morphology were assessed in shallow ice core samples. Stable isotope measurements plot either on the meteoric waterline or are enriched from it. The data cluster in two groups: the ice underlying the moraine has a δ2H:δ18O slope of 5.35 ± 0.92; ice from adjoining portions of Law Glacier has a slope of 6.69 ± 1.39. This enrichment pattern suggests the moraine's underlying blue ice entrained sediment through refreezing processes acting in an open system. Glaciological conditions favorable to warm-based sediment entrainment occur 30–50 km upstream. Basal melting and refreezing are further evidenced by abundant vapor figures formed from internal melting of the ice crystals. Both the moraine ice and Law Glacier are sufficiently depleted of heavy isotopes that their ice cannot be sourced locally, but instead must be derived from far-field interior regions of the higher polar plateau. Modeled ice flow speeds suggest the ice must be at least 80 ka old, with Law Glacier's ice possibly dating to OIS 5 and moraine ice older still.

Flow-frictional resistance at the base of glaciers and ice sheets is strongly linked to subglacial water pressure. Understanding the physical mechanisms that govern meltwater fluxes in subglacial channels is hence critical for constraining variations in ice flow. Previous mathematical descriptions of soft-bed subglacial channels assume a viscous till rheology, which is inconsistent with laboratory data and the majority of field studies. Here, we use a grain-scale numerical formulation coupled to pore-water dynamics to analyze the structural stability of channels carved into soft beds. Contrary to the soft-bed channel models assuming viscous till rheology, we show that the flanks of till channels can support substantial ice loads without creep closure of the channel, because the sediment has finite frictional strength. Increased normal stress on the channel flanks causes plastic failure of the sediment, and the channel rapidly shrinks to increase the ice-bed contact area. We derive a new parameterization for subglacial channelized flow on soft beds and show that channel dynamics are dominated by fluvial erosion and deposition processes with thresholds linked to the plastic rheology of subglacial tills. We infer that the described limits to channel size may cause subglacial drainage to arrange in networks of multiple closely spaced channels.