One fifth of Earth's volcanoes are covered by snow or ice and many have active geothermal systems that interact with the overlying ice. These glaciovolcanic interactions can melt voids into glaciers, and are subject to controls exerted by ice dynamics and geothermal heat output. Glaciovolcanic voids have been observed to form prior to volcanic eruptions, which raised concerns when such features were discovered within Job Glacier on Qw̓elqw̓elústen (Mount Meager Volcanic Complex), British Columbia, Canada. In this study we model the formation, evolution, and steady-state morphology of glaciovolcanic voids using analytical and numerical models. Analytical steady-state void geometries show cave height limited to one quarter of the ice thickness, while numerical model results suggest the void height h scales with ice thickness H and geothermal heat flux $\dot {Q}$ as $h/H = a H^b \dot {Q}^c$, with exponents b = −n/2 and c = 1/2 where n is the creep exponent. Applying this scaling to the glaciovolcanic voids within Job Glacier suggests the potential for total geothermal heat flux in excess of 10 MW. Our results show that relative changes in ice thickness are more influential in glaciovolcanic void formation and evolution than relative changes in geothermal heat flux.

We investigate unusual discontinuous glacier motion on Thompson Glacier, Umingmat Nunaat, Arctic Canada, using synthetic aperture radar (SAR) images and ice-flow modeling. A novel intensity-rescaling scheme is developed to reduce errors in high-resolution speckle tracking, resulting in a ~25% improvement in accuracy. Interferometric SAR (InSAR) and speckle tracking using high resolution RADARSAT-2 data indicate velocity discontinuities of up to 1 cm d−1 across deep and longitudinally extensive supraglacial channels on Thompson Glacier. We use a cross-sectional finite-element ice-flow model to determine the conditions under which velocity discontinuities of the observed magnitude and signature are possible. The modeling suggests that discontinuous motion across (long and straight) supraglacial channels can occur without ice fracture and under a wide variety of glacier thermal structures, including in fully temperate glaciers. Despite the wide range of conditions conducive to discontinuous motion, the form we observe requires that the associated channels be deep, longitudinally extensive and located in regions of lateral shearing. We speculate that these combined conditions are rare except on polythermal glaciers, where drainage features such as moulins are comparatively scarce and lower deformation rates allow channels to incise consistently and persist over many years.

Theory and observation show that glacier-flow regimes characterized by high basal slip enhance the projection of topographic detail to the surface, motivating this investigation into the efficacy of using glacier surges to improve bed estimation. Here we adapt a Bayesian inversion scheme and apply it to real and synthetic data as a proof of concept. Synthetic tests show a reduction in mean RMSE between true and inferred beds by more than half, and an increase in the mean correlation coefficient of ~0.5, when data from slip- versus deformation-dominated regimes are used. Multi-epoch inversions, which partition slip- and deformation-dominated regimes, are shown to outperform inversions that average over these flow regimes thereby squandering information. Tests with real data from a surging glacier in Yukon, Canada, corroborate these results, while highlighting the challenges of limited or inconsistent data. With the growing torrent of satellite-based observations, fast-flow events such as glacier surges offer potential to improve bed estimation for some of the world's most dynamic glaciers.

The Kaskawulsh Glacier is an iconic outlet draining the icefields of the St. Elias Mountains in Yukon, Canada. We determine and attempt to interpret its catchment-wide mass budget since 2007. Using SPOT5/6/7 data we estimate a 2007–18 geodetic balance of −0.46 ± 0.17 m w.e. a−1. We then compute balance fluxes and observed ice fluxes at nine flux gates to examine the discrepancy between the climatic mass balance and internal mass redistribution by glacier flow. Balance fluxes are computed using a fully distributed mass-balance model driven by downscaled and bias-corrected climate-reanalysis data. Observed fluxes are calculated using NASA ITS_LIVE surface velocities and glacier cross-sectional areas derived from ice-penetrating radar data. We find the glacier is still in the early stages of dynamic adjustment to its mass imbalance. We estimate a committed terminus retreat of ~23 km under the 2007–18 climate and a lower bound of 46 km3 of committed ice loss, equivalent to ~15% of the total glacier volume.

Radar sounding is a powerful geophysical approach for characterizing the subsurface conditions of terrestrial and planetary ice masses at local to global scales. As a result, a wide array of orbital, airborne, ground-based, and in situ instruments, platforms and data analysis approaches for radioglaciology have been developed, applied or proposed. Terrestrially, airborne radar sounding has been used in glaciology to observe ice thickness, basal topography and englacial layers for five decades. More recently, radar sounding data have also been exploited to estimate the extent and configuration of subglacial water, the geometry of subglacial bedforms and the subglacial and englacial thermal states of ice sheets. Planetary radar sounders have observed, or are planned to observe, the subsurfaces and near-surfaces of Mars, Earth's Moon, comets and the icy moons of Jupiter. In this review paper, and the thematic issue of the Annals of Glaciology on ‘Five decades of radioglaciology’ to which it belongs, we present recent advances in the fields of radar systems, missions, signal processing, data analysis, modeling and scientific interpretation. Our review presents progress in these fields since the last radio-glaciological Annals of Glaciology issue of 2014, the context of their history and future prospects.

Stationary ice-penetrating radar (sIPR) systems can be used to monitor temporal changes in electromagnetically sensitive properties of glaciers and ice sheets. We describe a system intended for autonomous operation in remote glacial environments, and document its performance during deployments in cold and temperate settings. The design is patterned after an existing impulse radar system, with the addition of a fibre-optic link and timing module to control transmitter pulses, a micro-UPS (uninterruptable power supply) to prevent uncontrolled system shutdown and a customized satellite telemetry scheme. Various implementations of the sIPR were deployed on the Kaskawulsh Glacier near an ice-marginal lake in Yukon, Canada, for 44–77 days in summers 2014, 2015 and 2017. Pronounced perturbations to englacial radiostratigraphy were observed commensurate with lake filling and drainage, and are interpreted as changes in englacial water storage. Another sIPR was deployed in 2015–2016 on ice island PII-A-1-f, which originated from the Petermann Glacier in northwest Greenland. This system operated autonomously for almost a year during which changes in thickness of the ice column were clearly detected.

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.

Accurately estimating winter surface mass balance on glaciers is central to assessing glacier health and predicting glacier run-off. However, measuring and modelling snow distribution is inherently difficult in mountainous terrain. Here, we explore rigorous statistical methods of estimating winter balance and its uncertainty from multiscale measurements of snow depth and density. In May 2016, we collected over 9000 manual measurements of snow depth across three glaciers in the St. Elias Mountains, Yukon, Canada. Linear regression, combined with cross-validation and Bayesian model averaging, as well as ordinary kriging are used to interpolate point-scale values to glacier-wide estimates of winter balance. Elevation and a wind-redistribution parameter exhibit the highest correlations with winter balance, but the relationship varies considerably between glaciers. A Monte Carlo analysis reveals that the interpolation itself introduces more uncertainty than the assignment of snow density or the representation of grid-scale variability. For our study glaciers, the winter balance uncertainty from all assessed sources ranges from 0.03 to 0.15 m w.e. (5–39%). Despite the challenges associated with estimating winter balance, our results are consistent with a regional-scale winter-balance gradient.

Surges are common in all the major ice caps in Iceland, and historical reports of surge occurrence go back several centuries. Data collection and regular observation over the last several decades have permitted a detailed description of several surges, from which it is possible to generalize on the nature of surging in Icelandic glaciers. Combining the historical records of glacier-front variations and recent field research, we summarize the geographic distribution of surging glaciers, their subglacial topography and geology, the frequency and duration of surges, changes in glacier surface geometry during the surge cycle, and measured velocity changes compared to calculated balance velocities. We note the indicators of surge onset and describe changes in ice, water and sediment fluxes during a surge. Surges accomplish a significant fraction of the total mass transport through the main outlet glaciers of ice caps in Iceland and have important implications for their hydrology. Our analysis of the data suggests that surge-type glaciers in Iceland are characterized by gently sloping surfaces and that they move too slowly to remain in balance given their accumulation rate. Surge frequency is neither regular nor clearly related to glacier size or mass balance. Steeply sloping glaciers, whether hard- or soft-bedded, seem to move sufficiently rapidly to keep in balance with the annual accumulation.

Outbursts of subglacial water from numerous alpine glaciers have been observed and documented. Such events tend to occur in spring and are thus attributed to an inability of the winter subglacial drainage system (characterized by high water pressure and low capacity) to accommodate a sudden and profuse influx of surface meltwater. Prior to a release event, bursts of glacier motion are common, and the release then precipitates the restoration of summer plumbing that damps or terminates surface acceleration. The events bear witness to the importance of interactions between surface melt, runoff, en-glacial water storage and internal routing, in addition to subglacial drainage morphology. Using a distributed numerical model to simultaneously solve surficial, englacial and subglacial water-transport equations, we investigate the role of these components in a hydro-mechanical event observed at Trapridge Glacier, YukonTerritory, Canada, in July 1990.

Digital elevation models of the bed and surface of Unteraargletscher, Switzerland, are used to reconstruct the theoretical pattern of basal water drainage for the years 1927, 1947, 1961 and 1997, during which period the glacier was thinning and receding. The theoretical drainage pattern for 1997 compares well, in a broad sense, with the locations of active moulins and the hydraulic connection status of boreholes drilled to the glacier bed. Changes in the basal water-flow pattern over the period 1927–97 that are revealed by the theoretical reconstructions of the subglacial drainage system structure are likely to have resulted from changes in glacier geometry. Concurrent with the retreat and thinning of the glacier, the height of medial moraines increased, probably due to the insulating effect of the debris cover reducing the melt of the underlying ice. This increase of moraine heights has led to the formation of hydraulic barriers at the glacier bed such that water flow has become channelized beneath the ice along drainage axes that parallel the course of the medial moraines on the glacier surface.

We describe a portable low-frequency impulse radar system intended for ground-based surveys that employs off-the-shelf hardware integrated with custom-designed software. The hardware comprises a 1–200 MHz transmitter, digitizer, computer and GPS receiver, which together weigh ∼1.5 kg. The entire system, including waterproof enclosures and batteries suited for >8 hours of continuous operation, weighs <10 kg plus the weight of the antenna housing. The system design is flexible, permitting hardware components such as the digitizer or navigation device to be exchanged. The software includes acquisition parameter control, real-time visual ice-depth rendering and data management capabilities using a hierarchical data format. The system described here has been successfully used to sound polythermal ice up to ∼220 m thick in ski-based surveys in the Yukon, Canada, and temperate ice up to ∼550 m thick in machine-based surveys in Iceland.

The association between basal hydrology and glacier sliding has become nearly synonymous with the early work of Almut Iken and colleagues. Their research published in the Journal of Glaciology from 1981 to 1986 made an indelible impact on the study of glacier hydromechanics by documenting strong correlations between basal water pressure and short-term ice-flow variations. With a passion for elucidating the physics of glacier-bed processes, Iken herself made fundamental contributions to our theoretical and empirical understanding of the sliding process. From the theoretical bound on basal shear stress, to the inferences drawn from detailed horizontal and vertical velocity measurements, the work of Iken and colleagues continues to inform the interpretation of data from alpine glaciers and has found increasing relevance to observations from the ice sheets.

We apply a time-dependent distributed glaciohydraulic model to Vatnajökull ice cap, Iceland, aiming to determine the large-scale subglacial drainage structure, the importance of basally derived meltwater, the influence of a permeable glacier bed and Vatnajökull’s discharge contribution to major rivers in Iceland. The model comprises two coupled layers that represent the subglacial horizon perched on a subsurface aquifer in the western sector and bedrock in the eastern sector. To initialize and drive the simulations, we use digital elevation models of the ice surface and bed, the 1999/2000 measured mass balance and an estimate of subglacial geothermal heat fluxes. The modelled subglacial flow field differs substantially from that derived by hydraulic-potential calculations, and the corresponding distribution of basal effective pressure shows a strong correlation between low effective pressure and surge-prone areas in northeastern and southern sectors of Vatnajökull. Simulations suggest that geothermally derived basal melt may account for up to ∼5% of the annual glacial discharge, and buried aquifers may evacuate up to ∼30% of subglacialwater.Time-dependent tests yield estimates of the glacial discharge component in various outlet rivers and suggest a possible seasonal migration of subglacial hydraulic divides. This study of present-day Vatnajökull hydrology forms the starting point for investigations of its future evolution.

We apply a point-scale energy-balance model to a small polythermal glacier in the St Elias Mountains of Canada in order to investigate the applicability and limitations of different treatments of the glacier surface temperature and subsurface heat flux. These treatments range in complexity from a multilayer subsurface model that simulates snowpack evolution, to the assumption of a constant glacier surface temperature equal to 0°C. The most sophisticated model includes dry densification of the snowpack, penetration of shortwave radiation into the subsurface, internal melting, refreezing of percolating meltwater and generation of slush layers. Measurements of subsurface temperature and surface lowering are used for model validation, and highlight the importance of including subsurface penetration of shortwave radiation in the model. Using an iterative scheme to solve for the subsurface heat flux as the residual of the energy-balance equation results in an overestimation of total ablation by 18%, while the multilayer subsurface model underestimates ablation by 6%. By comparison, the 0°C surface assumption leads to an overestimation of ablation of 29% in this study where the mean annual air temperature is about −8°C.

We investigate the dynamics of a small surge-type valley glacier as part of a study to characterize glacier response to climate in the Donjek Range, southwest Yukon, Canada. Pole displacements were measured using kinematic GPS techniques during three consecutive summer field seasons. Measured surface velocities range from <10 m a−1 over the lower 1500 m of the 5 km long glacier to a maximum of ∼25–35 m a−1 over the upper 3500 m. Basal velocities along an approximate flowline are reconstructed from the measured surface velocities using inverse methods. Control tests are used to validate the inversion scheme, and sensitivity tests are performed to evaluate the influence of the flow-law coefficient, shape factor and longitudinal averaging length. Inversion of the real data shows that basal motion accounts for 50–100% of the total surface motion along the flowline. Based on these results, and several other lines of evidence, we suggest this glacier may be undergoing a slow surge.

Significant intra-annual variability in flow rates of tidewater-terminating Arctic glaciers has been observed in recent years. These changes may result from oceanic and/or atmospheric forcing through (1) perturbations at the terminus, such as enhanced submarine melt and changes in sea-ice buttressing, or (2) increased surface melt, in response to atmospheric warming, reaching the bed and promoting glacier slip. We examine the influence of these processes on Belcher Glacier, a large fast-flowing tidewater outlet of the Devon Island ice cap in the Canadian Arctic. A hydrologically-coupled higher-order ice flow model is used to estimate changes in glacier flow speed as a result of changes in sea-ice buttressing and hydrologically-driven melt-season dynamics. Daily run-off from five sub-catchments over the 2008 and 2009 melt seasons provides meltwater forcing for the model simulations. Model results are compared with remotely-sensed and in situ ice-surface velocity measurements. Sea-ice effects are found to have a minor influence on glacier flow speed relative to that of meltwater drainage, which is clearly implicated in short-term velocity variations during the melt season. We find that threshold drainage is essential in determining the timing of these short-lived accelerations.

Bed-reflection power (BRP) from ice-penetrating radar has been used to make inferences about subglacial conditions and processes, yet is subject to confounding influences, including englacial attenuation and bed geometry. We use radar data collected in 2008–11 from a polythermal glacier to compute BRP with the aim of relating BRP to basal conditions. We examine the relationship between raw BRP and ice thickness, apparent bed slope and thickness of the englacial scattering layer as a proxy for internal reflection power. We then analyze a corrected form of the BRP with a graph-segmentation algorithm to delineate areas of high and low reflection power. Low corrected BRP values are found near the glacier terminus where the bed is most likely to be cold, while high corrected BRP is found in the region thought to be undergoing a slow surge. We find a spatial correlation between high BRP and high values of subglacial hydraulic upstream area, suggestive of a hydrological control on BRP. Whereas in dominantly cold glaciers BRP seems to distinguish cold from temperate regions of the bed, BRP in a polythermal glacier with a substantial volume of temperate ice may be a more complex product of thermal and hydrological conditions.

The subglacial chemical weathering environment is largely controlled by low temperatures and the presence of freshly comminuted minerals with a high surface area. These characteristics are believed to promote dissolution processes that give rise to low silica and high Ca2+ fluxes emanating from glacierized basins. We test an alternative hypothesis, that mineral precipitation reactions in the subglacial environment play an equally important role in controlling the water chemistry in glacierized basins. We analyze borehole and proglacial water chemistry from a subarctic polythermal glacier, complemented by mineral XRD analysis of suspended sediment, till and bedrock samples. In conjunction with a thermodynamic analysis of the water and mineral chemistry, we use reaction-path modelling to study the chemical enrichment of water through the glacier system. We find that the high pH of the subglacial environment is conducive to secondary mineral precipitation, and that it is not possible to balance the water chemistry using dissolution reactions alone. We show that low silica can be explained by standard weathering reactions without having to invoke mineral-leaching reactions. Our results suggest that subglacial weathering intensity may be significantly underestimated if the production of secondary minerals is not considered.