The aim of this paper is to characterise the internal structures and ice-flow history of representative valley glaciers in Svalbard and infer from them dynamic changes over centennial timescales. Three polythermal and one cold valley glacier are investigated using field- and laboratory-based techniques and remote sensing. Structures along flow-unit boundaries indicate that ice-flow configuration in three of the glaciers has remained stable spanning the residence time of the ice. Deformation of a flow-unit boundary in the fourth reveals an ice-flow instability, albeit one that has been maintained since its most recent advance. Macro-crystallographic, sedimentological and isotopic analyses indicate that basal ice is elevated to the glacier surface, as shown by entrained sediments and enrichment in heavy isotopes. In narrow zones of enhanced cumulative strain, new ice facies are generated through dynamic recrystallisation. The surface density of longitudinal foliation is shown to represent the relative magnitude of cumulative strain. Geometric similarities between flow-unit boundaries in Svalbard valley glaciers and larger scale longitudinal surface structures in ice sheets suggest that deformation mechanisms are common to both.

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.

The density structure of firn has implications for hydrological and climate modelling, and ice-shelf stability. The structure of firn can be evaluated from depth models of seismic velocity, widely obtained with Herglotz–Wiechert inversion (HWI), an approach that considers the slowness of refracted seismic arrivals. However, HWI is strictly appropriate only for steady-state firn profiles and the inversion accuracy can be compromised where firn contains ice layers. In these cases, full waveform inversion (FWI) may yield more success than HWI. FWI extends HWI capabilities by considering the full seismic waveform and incorporates reflected arrivals. Using synthetic firn density profiles, assuming both steady- and non-steady-state accumulation, we show that FWI outperforms HWI for detecting ice slab boundaries (5–80 m thick, 5–80 m deep) and velocity anomalies within firn. FWI can detect slabs thicker than one wavelength (here, 20 m, assuming a maximum frequency of 60 Hz) but requires the starting velocity model to be accurate to ±2.5%. We recommend for field practice that the shallowest layers of velocity models are constrained with ground-truth data. Nonetheless, FWI shows advantages over established methods, and should be considered when the characterisation of firn ice slabs is the goal of the seismic survey.

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.

The debris that covers the ablation areas of high-elevation debris-covered glaciers contributes to the distinctive features and processes occurring both on and within such glaciers. Despite recent advances, knowledge of the subsurface environments of high-elevation debris-covered glaciers is still extremely limited. In particular, targeted field-based data are needed to parameterise and refine the projections of these glaciers in numerical models. Here, we outline the current understanding of the internal properties of high-elevation debris-covered glaciers based on direct field-based methods and suggest future research directions for field-based studies.

The density structure of firn has implications for hydrological and climate modelling and for ice shelf stability. The firn structure can be evaluated from depth models of seismic velocity, widely obtained with Herglotz-Wiechert inversion (HWI), an approach that considers the slowness of refracted seismic arrivals. However, HWI is appropriate only for steady-state firn profiles and the inversion accuracy can be compromised where firn contains ice layers. In these cases, Full Waveform Inversion (FWI) can be more successful than HWI. FWI extends HWI capabilities by considering the full seismic waveform and incorporates reflected arrivals, thus offering a more accurate estimate of a velocity profile. We show the FWI characterisation of the velocity model has an error of only 1.7% for regions (vs. 4.2% with HWI) with an ice slab (20 m thick, 40 m deep) in an otherwise steady-state firn profile.

Characterising the structures within glaciers can give unique insight into ice motion processes. On debris-covered glaciers, traditional structural glaciological mapping is challenging because the lower glacier is hidden by the supraglacial debris layer. Here, we use high-resolution optical televiewer (OPTV) image logs from four boreholes drilled into Khumbu Glacier, Nepal, to overcome this limitation and investigate englacial structural features within a Himalayan debris-covered glacier. The OPTV logs show structural features that are up to an order of magnitude thinner than those observed at the glacier surface and reveal five structural units: (I) primary stratification of ice; (II) debris-rich planes that conform with the primary stratification; (III) water-healed crevasse traces; (IV) healed crevasse traces; and (V) steeply dipping planes of basally derived fine sediment near the glacier terminus. The OPTV logs also reveal that the primary stratification both decreases in dip with depth (by up to 56° over 20 m) and rotates with depth (by up to 100° over 20 m) towards parallelism with the proximal lateral moraine. This transformation and the presence of relict layers of basally derived sediment raised into an englacial position – possibly involving thrusting – near the glacier's now stagnant terminus reveal a previously more dynamic glacier regime.

Subglacial hydrology modulates basal motion but remains poorly constrained, particularly for soft-bedded Greenlandic outlet glaciers. Here, we report detailed measurements of the response of subglacial water pressure to the connection and drainage of adjacent water-filled boreholes drilled through kilometre-thick ice on Sermeq Kujalleq (Store Glacier). These measurements provide evidence for gap opening at the ice-sediment interface, Darcian flow through the sediment layer, and the forcing of water pressure in hydraulically-isolated cavities by stress transfer. We observed a small pressure drop followed by a large pressure rise in response to the connection of an adjacent borehole, consistent with the propagation of a flexural wave within the ice and underlying deformable sediment. We interpret the delayed pressure rise as evidence of no pre-existing conduit and the progressive decrease in hydraulic transmissivity as the closure of a narrow (< 1.5 mm) gap opened at the ice-sediment interface, and a reversion to Darcian flow through the sediment layer with a hydraulic conductivity of ≤ 10−6 m s−1. We suggest that gap opening at the ice-sediment interface deserves further attention as it will occur naturally in response to the rapid pressurisation of water at the bed.

Subglacial hydrological systems require innovative technological solutions to access and observe. Wireless sensor platforms can be used to collect and return data, but their performance in deep and fast-moving ice requires quantification. We report experimental results from Cryoegg: a spherical probe that can be deployed into a borehole or moulin and transit through the subglacial hydrological system. The probe measures temperature, pressure and electrical conductivity in situ and returns all data wirelessly via a radio link. We demonstrate Cryoegg's utility in studying englacial channels and moulins, including in situ salt dilution gauging. Cryoegg uses VHF radio to transmit data to a surface receiving array. We demonstrate transmission through up to 1.3 km of cold ice – a significant improvement on the previous design. The wireless transmission uses Wireless M-Bus on 169 MHz; we present a simple radio link budget model for its performance in cold ice and experimentally confirm its validity. Cryoegg has also been tested successfully in temperate ice. The battery capacity should allow measurements to be made every 2 h for more than a year. Future iterations of the radio system will enable Cryoegg to transmit data through up to 2.5 km of ice.

Direct measurements of spatially distributed vertical strain within ice masses are scientifically valuable but challenging to acquire. We use manual marker tracking and automatic cross correlation between two repeat optical televiewer (OPTV) images of an ~100 m-long borehole at Derwael Ice Rise (DIR), Antarctica, to reconstruct discretised, vertical strain rate and velocity at millimetre resolution. The resulting profiles decay with depth, from −0.07 a−1 at the surface to ~−0.002 a−1 towards the base in strain and from −1.3 m a−1 at the surface to ~−0.5 m a−1 towards the base in velocity. Both profiles also show substantial local variability. Three coffee-can markers installed at different depths into adjacent boreholes record consistent strain rates and velocities, although averaged over longer depth ranges and subject to greater uncertainty. Measured strain-rate profiles generally compare closely with output from a 2-D ice-flow model, while the former additionally reveal substantial high-resolution variability. We conclude that repeat OPTV borehole logging represents an effective means of measuring distributed vertical strain at millimetre scale, revealing high-resolution variability along the uppermost ~100 m of DIR, Antarctica.

While hot-water drilling is a well-established technique used to access the subsurface of ice masses, drilling into high-elevation (≳ 4000 m a.s.l.) debris-covered glaciers faces specific challenges. First, restricted transport capacity limits individual equipment items to a volume and mass that can be slung by small helicopters. Second, low atmospheric oxygen and pressure reduces the effectiveness of combustion, limiting a system's ability to pump and heat water. Third, thick supraglacial debris, which is both highly uneven and unstable, inhibits direct access to the ice surface, hinders the manoeuvring of equipment and limits secure sites for equipment placement. Fourth, englacial debris can slow the drilling rate such that continued drilling becomes impracticable and/or boreholes deviate substantially from vertical. Because of these challenges, field-based englacial and subglacial data required to calibrate numerical models of high-elevation debris-covered glaciers are scarce or absent. Here, we summarise our experiences of hot-water drilling over two field seasons (2017–2018) at the debris-covered Khumbu Glacier, Nepal, where we melted 27 boreholes up to 192 m length, at elevations between 4900 and 5200 m a.s.l. We describe the drilling equipment and operation, evaluate the effectiveness of our approach and suggest equipment and methodological adaptations for future use.

The phase-sensitive radio-echo sounder (pRES) is a powerful new instrument that can measure the depth of internal layers and the glacier bed to millimetre accuracy. We use a stationary 16-antenna pRES array on Store Glacier in West Greenland to measure the three-dimensional orientation of dipping internal reflectors, extending the capabilities of pRES beyond conventional depth sounding. This novel technique portrays the effectiveness of pRES in deriving the orientation of dipping internal layers that may complement profiles obtained through other geophysical surveying methods. Deriving ice vertical strain rates from changes in layer depth as measured by a sequence of pRES observations assumes that the internal reflections come from vertically beneath the antenna. By revealing the orientation of internal reflectors and the potential deviation from nadir of their associated reflections, the use of an antenna array can correct this assumption. While the array configuration was able to resolve the geometry of englacial layers, the same configuration could not be used to accurately image the glacier bed. Here, we use simulations of the performance of different array geometries to identify configurations that can be tailored to study different types of basal geometry for future deployments.

We have drilled 13 boreholes within and around a through-cutting rift on the (unofficially named) Roi Baudouin Ice Shelf, East Antarctica. Logging by optical televiewer (OPTV) combined with core inspection has resulted in the identification and characterization of several material facies. Outside the rift, OPTV-imaged annual layering indicates ~150 years of accumulation over the 66 m length of one of the boreholes. Luminosity analysis of this image also reveals the presence of numerous dark melt layers as well as a systematic decrease in background luminosity, interpreted in terms of a progressive increase in light transmission during firnification. We identify four material facies within the rift: snow, granular ice, marine ice and unconsolidated platelets. We interpret the granular ice facies as snow that has been saturated by percolating sea water, and the underlying marine ice as compacted buoyant platelets that have adhered to the rift base. Core sections reveal the presence of tubular channels within the marine ice, indicating that it is macroporous and permeable to sea water. The lower boundary of this facies merges into a mushy layer of unconsolidated platelets that were successfully imaged by OPTV, revealing irregular sub-horizontal layering similar to that reported previously on the basis of (directional) borehole video.

Ice-penetrating radar and modelling data are presented suggesting the presence of a zone of temperate ice, water ponding or saturated sediment beneath the tongue of Taylor Glacier, Dry Valleys, Antarctica. The proposed subglacial zone lies 3–6 km up-glacier of the terminus and is 400– 1000m across. The zone coincides with an extensive topographic overdeepening to 80m below sea level. High values of residual bed reflective power across this zone compared to other regions and the margins of the glacier require a high dielectric contrast between the ice and the bed and are strongly indicative of the presence of basal water or saturated sediment. Analysis of the hydraulic equipotential surface also indicates strong convergence into this zone of subglacial water flow paths. However, thermodynamic modelling reveals that basal temperatures in this region could not exceed –7˚C relative to the pressure-melting point. Such a result is at odds with the radar observations unless the subglacial water is a hypersaline brine.

The physical properties of eight ice cores recovered from along a flowline at Glacier de Tsanfleuron, Switzerland, have led to the identification of three distinctive internal zones. We use variations in the bulk ionic chemistry of these zones to approximate their relative liquid-water concentrations and ice viscosities. Results suggest that relative bulk water concentration and ice softness vary by over an order of magnitude between the zones. Implications of this variability for predictions of the glacier’s response to climate change are evaluated by incorporating these relative softnesses into a multi-layered (two-dimensional) model of ice flow. Model output is compared with that from an identical model constrained with a spatially uniform ice viscosity under advance and retreat modelling scenarios. The former scenario is used to tune viscosity by growing a glacier to its present long-section geometry, resulting in best-fit ice hardness values of 1.2 a1 bar–3 for the englacial ice in the multi-layered model and 7.0 a–1 bar–3 for all of the ice in the single-layered model. Both result in close approximations to the current long profile, yielding rms deviationsbetween measured and modelled ice thicknesses that are 5 5 m. In contrast, a single-layered model constrained with a hardness of 1.2 a–1 bar–3 overestimates the current measured long-section area by 31%, having a rms ice-thickness error of 15.0 m. Under the retreat modelling scenario, which gauges the response of the glacier to an imposed 75 m rise in equilibrium-line altitude (ELA), the multi-layered model predicts a long-section area reduction of 78%. This contrasts with a reduction of 64% for the single-layered model (hardness = 7.0 a–1 bar–3) and 85% for the single-layered model (hardness = 1.2 a–1 bar–3).

Records of surface motion, englacial tilt and repeat inclinometry are used to determine patterns of surface, internal and basal motion across the tongue of Haut Glacier d’Arolla, Switzerland, over temporal scales ranging from days to months. Findings are interpreted with reference to contemporaneous measurements of subglacial water pressures, and prior knowledge of the glacier’s subglacial drainage-system structure. Long-term inclinometry results show pronounced extrusion flow over a subglacial drainage axis, with basal velocities up to twice those measured at the glacier surface. Deformation profiles are more conventional away from the drainage axis, with basal velocities ∼60–70% of surface velocities. Comparison of long-term tilt rates from repeat inclinometry and englacial tiltmeters shows close correspondence. Englacial tiltmeter data are used to reconstruct internal velocity profiles and to split surface velocities into internal deformation and basal motion contributions over spring, summer and autumn/winter periods. Although, spatial patterns of surface movement are similar between periods, patterns of internal and basal motion are not. Results are interpreted in terms of the location of sticky and slippery spots, with temporally changing patterns of basal drag reflecting changing distributions of water pressure.

Channels incised into bedrock, or Nye channels, often form an important component of subglacial drainage at temperate glaciers, and their structure exerts control over patterns and rates of (a) channel erosion, (b) water flow-velocity and (c) water pressure. The latter, in turn, exerts a strong control over basal traction and, thus, ice dynamics. In order to investigate these controls, it is necessary to quantify detailed flow processes in subglacial Nye channels. However, it is effectively impossible to acquire such measurements from fully pressurized, subglacial channels. To solve this problem, we here apply a three-dimensional, finite-volume solution of the Reynolds averaged Navier– Stokes (RANS) equations with a one-equation mixing-length turbulence closure to simulate flow in a 3 m long section of an active Nye channel located in the immediate foreground of Glacier de Tsanfleuron, Switzerland. Numerical model output permits high-resolution visualization of water flow through the channel reach, and enables evaluation of the experimental manipulation of the pressure field adopted across the overlying ice lid. This yields an increased theoretical understanding of the hydraulic behaviour of Nye channels, and, in the future, of their effect on glacier drainage, geomorphology and ice dynamics.

Fourteen subglacial debris samples have been recovered from the margins of, or beneath, Haut Glacier d’Arolla, Switzerland. The grain-size distributions of these sediments are presented and compared with each other as bivariate plots of percentage weight against (sieve-defined) particle size and log number of particles against log particle diameter. All of the samples recovered are composed of a broad range of clast sizes and approach self-similarity over the four orders of magnitude of grain-sizes analysed. Fractal dimensions range from 2.47 to 2.77. Sample intercomparison reveals the operation of at least two processes of textural evolution: the production of fines by in-situ weathering, interpreted in terms of abrasion associated with subglacial sediment deformation, and the loss of fines, interpreted in terms of eluviation by percolating subglacial meltwaters. These interpretations are supported and refined through comparison of the grain-size fractions gained (in the case of deformation) and lost (in the case of eluviation) with those fractions respectively generated in a laboratory-based simulation of sediment deformation and exiting the glacier suspended in the proglacial meltwater stream. While sediment deformation has the effect of increasing the fine fraction between 0 and 10ϕ and of raising the fractal dimension of undeformed sediments from 2.47 to 2.77, eluviation removes particles between 2 and 100, driving the fractal dimension of deformed sediments down from 2.77 to 2.54. These fractal dimensions are generally lower than those recorded at other comparable glaciers, consistent with the relatively low rates of sediment deformation inferred from other studies at Haut Glacier d’Arolla.

Glen’s flow law is a well-established general law for steady-state glacier ice deformation, and many laboratory tests and field measurements have been undertaken which have shown the generality of the law to be correct. In Nature, ice deformation is the response of the glacier/ice sheet to the applied self-weight stress of the ice mass (i.e. ice thickness, gravity and ice density) which produces a stress gradient within the ice column. Detailed experimental analyses of ice samples in the laboratory have until now only been undertaken using uniform stress fields in uniaxial or triaxial tests. Obviously the best method for investigating ice in the laboratory would be if stress gradients similar to those found in Nature could be replicated. In the following paper we describe the physical modelling of two (laboratory-prepared) isotropic, polycrystalline ice models (0.75 × 0.25 × 0.18 m) at enhanced gravity levels (80g) in a geotechnical beam centrifuge. Steel plate was placed on top of the ice model, replicating an overburden of approximately 36 m of ice (at 80g). Thus we were able to model the deformation of the lower 14 m of an ice mass approximately 50 m thick. Models are confined laterally by the Perspex strongbox walls, preventing lateral extension within the sample during testing. Models are unconfined on their downslope ends, rendering longitudinal stresses negligible. Deformation can therefore be treated as simple shear. Samples are instrumented with displacement markers and thermocouples. Values for A and n in the flow law derived from the experiments are reasonable and indicate the potential of this method for ice-deformation studies.

The mass balance of polythermal ice masses is critically dependent on the proportion of surface-generated meltwater that subsequently refreezes in the snowpack and firn. In order to quantify this effect and to characterize its spatial variability, we measured near-surface (<10 m) snow and firn densities at an elevation of ~1945ma.s.l. in the percolation zone of the Greenland ice sheet in spring and autumn 2004. Results indicate that local snowpack depth above the previous end-of-summer 2003 melt surface increased by ±5% (7.6 cm) from spring to autumn while, over the same period, snowpack density increased by >26%, resulting in a 32% increase in net accumulation. This ‘seasonal densification’ increased at lower elevations, rising to 47% 10 km closer to the ice-sheet margin at 1860ma.s.l. Density/depth profiles from nine sites within 1 km2 at ~1945ma.s.l. reveal complex stratigraphies that change over short spatial scales and seasonally. We conclude that estimates of mass-balance change cannot be calculated solely from observed changes in surface elevation, but that near-surface densification must also be considered. However, predicting spatial and temporal variations in densification may not be straightforward. Further, the development of complex firn-density profiles both masks discernible annual layers in the near-surface firn and ice stratigraphy and is likely to introduce error into radar-derived estimates of surface elevation.