Subglacial seismicity presents the opportunity to monitor inaccessible glacial beds at the epicentral location and time. Glaciers can be underlain by rock or till, a first order control on bed mechanics. Velocity-weakening, necessary for unstable slip, has been shown for each bed type, but is much stronger and evolves over more than an order of magnitude longer distances for till beds. Utilizing a de-stiffened double direct shear apparatus, we found conditions for instability at freezing temperatures and high slip rates for both bed types. During stick–slip stress-drops, we recorded acoustic emissions with piezoelectric transducers frozen into the ice. The two populations of event waveforms appear visually similar and overlap in their statistical features. We implemented a suite of supervised machine learning algorithms to classify the bed type of recorded waveforms and spectra, with prediction accuracy between 65–80%. The Random Forest Classifier is interpretable, showing the importance of initial oscillation peaks and higher frequency energy. Till beds have generally higher friction and resulting stress-drops, with more impulsive first arrivals and more high frequency content compared to rock emissions, but rock beds can produce many till-like events. Seismic signatures could enhance interpretation of bed conditions and mechanics from subglacial seismicity.

A 198.8 m deep borehole was drilled through ice to subglacial bedrock in the northwestern marginal part of Princess Elizabeth Land, ~12 km south of Zhongshan Station, in January–February 2019. Three years later, in February 2022, the borehole temperature profile was measured, and the geothermal heat flow (GHF) was estimated using a 1-D time-dependent energy-balance equation. For a depth corresponding to the base of the ice sheet, the GHF was calculated as 72.6 ± 2.3 mW m−2 and temperature −4.53 ± 0.27°C. The regional averages estimated for this area based, generally, on tectonic setting vary from 55 to 66 mW m−2. A higher GHF is interpreted to originate mostly from the occurrence of metamorphic complexes intruded by heat-producing elements in the subglacial bedrock below the drill site.

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.

Knowledge of glacier volume is crucial for ice flow modelling and predicting the impacts of climate change on glaciers. Rugged terrain, harsh weather conditions and logistic costs limit field-based ice thickness observations in the Himalaya. Remote-sensing applications, together with mathematical models, provide alternative techniques for glacier ice thickness and volume estimation. The objective of the present research is to assess the application of artificial neural network (ANN) modelling coupled with remote-sensing techniques to estimate ice thickness on individual glaciers with direct field measurements. We have developed two ANN models and estimated the ice thickness of Chhota Shigri Glacier (western Himalaya) on ten transverse cross sections and two longitudinal sections. The ANN model estimates agree well with ice thickness measurements from a ground-penetrating radar, available for five transverse cross sections on Chhota Shigri Glacier. The overall root mean square errors of the two ANN models are 24 and 13 m and the mean bias errors are ±13 and ±6 m, respectively, which are significantly lower than for other available models. The estimated mean ice thickness and volume for Chhota Shigri Glacier are 109 ± 17 m and 1.69 ± 0.26 km3, respectively.

Drilling to the bedrock of ice sheets and glaciers offers unique opportunities for examining the processes occurring in the bed. Basal and subglacial materials contain important paleoclimatic and paleoenvironmental records and provide a unique habitat for life; they offer significant information regarding the sediment deformation beneath glaciers and its effects on the subglacial hydraulic system and geology. The newly developed and tested Antarctic subglacial drilling rig (ASDR) is designed to recover ice and bedrock core samples from depths of up to 1400 m. All of the drilling equipment is installed inside a movable, sledge-mounted, temperature-controlled and wind-protected drilling shelter and workshop. To facilitate helicopter unloading of the research vessel, the shelter and workshop can be disassembled, with individual parts weighing <2–3 tons. The entire ASDR system weighs ~55 tons, including transport packaging. The ASDR is designed to be transported to the chosen site via snow vehicles and would be ready for drilling operations within 2–3 d after arrival. The ASDR was tested during the 2018–2019 summer season near Zhongshan Station, East Antarctica. At the test site, 2-week drilling operations resulted in a borehole that reached bedrock at a depth of 198 m.

A new, modified version of the cable-suspended Ice and Bedrock Electromechanical Drill (IBED) was designed for drilling in firn, ice, debris-rich ice and rock. The upper part of the drill is almost the same for all drill variants and comprises four sections: cable termination, a slip-ring section, an antitorque system and an electronic pressure chamber. The lower part of the IBED comprises an auger core barrel, reamers, a core barrel for ice/debris-ice drilling and a conventional geological single-tube core barrel or custom-made double-tube core barrel. First, the short and full-scale field versions of the IBED were tested at an outdoor testing stand and a testing facility with a 12.5 m-deep ice well. Then, in the 2018–2019 summer season, the IBED was tested in the field at a site ~12 km south of Zhongshan Station, East Antarctica, and a ~6 cm bedrock core was recovered from a 198 m-deep borehole. A total of 18 d was required to penetrate the ice sheet. The retrieved core samples of blue ice, basal ice and bedrock provided valuable information regarding the Earth's paleo-environment.

Recent advances in the measurement of bedforms beneath active ice streams have been made using ground-based grid profiling using impulse radar systems operating with centre frequencies in the 3–5 MHz range. Surveys of Rutford Ice Stream and Pine Island Glacier have shown that features such as mega-scale glacial lineations with topographic relief of as little as 3 m can be traced for many kilometres downstream under more than 2 km of fast-moving ice. In the discussion of these data, it is often asked ‘How is it possible to map such fine-scale topography with such a low-frequency radar’. In answering that question, the key point is the distinction between the precision of a radar range measurement to a single, isolated reflective interface and the ability to resolve the presence of two closely-spaced interfaces of similar reflectivity (commonly referred to as the vertical resolution). This paper will discuss and illustrate this distinction and use the case study of data acquired over Pine Island Glacier to examine the limits of precision of the radar range measurement.

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.

Radio-echo sounding (RES) can be used to understand ice-sheet processes, englacial flow structures and bed properties, making it one of the most popular tools in glaciological exploration. However, RES data are often subject to ‘strip noise’, caused by internal instrument noise and interference, and/or external environmental interference, which can hamper measurement and interpretation. For example, strip noise can result in reduced power from the bed, affecting the quality of ice thickness measurements and the characterization of subglacial conditions. Here, we present a method for removing strip noise based on combined wavelet and two-dimensional (2-D) Fourier filtering. First, we implement discrete wavelet decomposition on RES data to obtain multi-scale wavelet components. Then, 2-D discrete Fourier transform (DFT) spectral analysis is performed on components containing the noise. In the Fourier domain, the 2-D DFT spectrum of strip noise keeps its linear features and can be removed with a ‘targeted masking’ operation. Finally, inverse wavelet transforms are performed on all wavelet components, including strip-removed components, to restore the data with enhanced fidelity. Model tests and field-data processing demonstrate the method removes strip noise well and, incidentally, can remove the strong first reflector from the ice surface, thus improving the general quality of radar data.

Subglacial hydrologic systems in Antarctica and Greenland play a fundamental role in ice-sheet dynamics, yet critical aspects of these systems remain poorly understood due to a lack of observations. Ground-based electromagnetic (EM) geophysical methods are established for mapping groundwater in many environments, but have never been applied to imaging lakes beneath ice sheets. Here, we study the feasibility of passive- and active-source EM imaging for quantifying the nature of subglacial water systems beneath ice streams, with an emphasis on the interfaces between ice and basal meltwater, as well as deeper groundwater in the underlying sediments. We describe a suite of model studies that exam the data sensitivity as a function of ice thickness, water conductivity and hydrologic system geometry for models representative of a subglacial lake and a grounding zone estuary. We show that EM data are directly sensitive to groundwater and can image its lateral and depth extent. By combining the conductivity obtained from EM data with ice thickness and geological structure from conventional geophysical techniques, such as ground-penetrating radar and active seismic surveying, EM data have the potential to provide new insights on the interaction between ice, rock and water at critical ice-sheet boundaries.

There is significant interest in sampling subglacial environments for geobiological studies, but they are difficult to access. Existing ice-drilling technologies make it cumbersome to maintain microbiologically clean access for sample acquisition and environmental stewardship of potentially fragile subglacial aquatic ecosystems. The IceMole is a maneuverable subsurface ice probe for clean in situ analysis and sampling of glacial ice and subglacial materials. The design is based on the novel concept of combining melting and mechanical propulsion. It can change melting direction by differential heating of the melting head and optional side-wall heaters. The first two prototypes were successfully tested between 2010 and 2012 on glaciers in Switzerland and Iceland. They demonstrated downward, horizontal and upward melting, as well as curve driving and dirt layer penetration. A more advanced probe is currently under development as part of the Enceladus Explorer (EnEx) project. It offers systems for obstacle avoidance, target detection, and navigation in ice. For the EnEx-IceMole, we will pay particular attention to clean protocols for the sampling of subglacial materials for biogeochemical analysis. We plan to use this probe for clean access into a unique subglacial aquatic environment at Blood Falls, Antarctica, with return of a subglacial brine sample.

In the early hours of 25 December 2012, an attempt to explore Subglacial Lake Ellsworth, West Antarctica, using a specially designed hot-water drill, was halted. This UK project, involving several universities, the British Antarctic Survey and the National Oceanography Centre, had been in planning for 10 years. The project developed a full blueprint for subglacial lakes research, involving access to the subglacial environment through deep drilling, direct measurement and sampling of water and sediment by the construction of a probe and sediment corer, and environmental protocols to ensure cleanliness in line with international agreements on stewardship and protection of subglacial systems. Drilling was ceased after the main borehole failed to link with a subsurface cavity of water, built up over ∽40 hours. Without this link, insufficient water was available to continue drilling downwards to the lake, ∽3000 m beneath the surface. On return to the UK, an external review of the programme was undertaken to formally assess the reasons for the fieldwork failure, and to make recommendations on the modifications necessary for success. From this review, the Lake Ellsworth programme formulated a pathway along which a second attempt to explore the lake can be developed. Here details of the Lake Ellsworth field experiment, the circumstances that led to its failure and the corrections required are presented. Hot-water drilling is still regarded as the only feasible scheme for assuring clean access to the subglacial environment. The lessons learned from the Lake Ellsworth experience are substantial, however, and demonstrate that considerable technological and methodological advances are necessary for successful future research on subglacial lakes beneath thick (>2 km) ice.

Glacier dynamics are inextricably linked to the basal conditions of glaciers. Seismic reflection methods can image the glacier bed under certain conditions. However, where a seismically thin layer of material is present at the bed, traditional analyses may fail to fully characterize bed properties. We use a targeted full-waveform inversion algorithm to quantify the basal-layer parameters of a mountain glacier: thickness (d), P-wave velocity (α) and density (ρ). We simultaneously invert for the seismic quality factor (Q) of the bulk glacier ice. The inversion seeks to minimize the difference between the data and a one-dimensional reflectivity algorithm using a gradient-based search with starting values initialized from a Monte Carlo scheme. We test the inversion algorithm on four basal layer synthetic data traces with 5% added Gaussian noise. The inversion retrieved thin-layer parameters within 10% of synthetic test parameters with the exception of seismic Q. For the seismic dataset from Bench Glacier, Alaska, USA, inversion results indicate a thin basal layer of debris-rich ice within the study area having mean velocity 4000 ± 700 m s–1, density 1900 ± 200 kg m–3 and thickness 6 ± 1.5 m.