An array of information about the Antarctic ice sheet can be extracted from ice-sheet internal architecture imaged by airborne ice-penetrating radar surveys. We identify, trace and date three key internal reflection horizons (IRHs) across multiple radar surveys from South Pole to Dome A, East Antarctica. Ages of ~38 ± 2.2, ~90 ± 3.6 and ~162 ± 6.7 ka are assigned to the three IRHs, with verification of the upper IRH age from the South Pole ice core. The resultant englacial stratigraphy is used to identify the locations of the oldest ice, specifically in the upper Byrd Glacier catchment and the Gamburtsev Subglacial Mountains. The distinct glaciological conditions of the Gamburtsev Mountains, including slower ice flow, low geothermal heat flux and frozen base, make it the more likely to host the oldest ice. We also observe a distinct drawdown of IRH geometry around South Pole, indicative of melting from enhanced geothermal heat flux or the removal of deeper, older ice under a previous faster ice flow regime. Our traced IRHs underpin the wider objective to develop a continental-scale database of IRHs which will constrain and validate future ice-sheet modelling and the history of the Antarctic ice sheet.

One of the largest contributors to uncertainty in predictions of sea-level rise from ice-sheet models is a lack of knowledge about the bed topography beneath ice sheets. Bed topography maps are normally made by interpolating between linear radar surveys using methods that include kriging, mass conservation and flowline diffusion, all of which may miss influential mesoscale (2–30 km) bedforms. Previous works have explored an Ice-Flow Perturbation Analysis (IFPA) approach for estimating bed topography using the surface expression of these mesoscale bedforms. Using regions of Pine Island Glacier that have been intensively surveyed by ice-penetrating radar as test sites, and a refined IFPA methodology, we find that IFPA detects bedforms capable of influencing ice flow which are not represented in Bedmachine Antarctica and other interpolated bed products. We further explore the ability of IFPA to estimate relative bed slipperiness, finding higher slipperiness in the main trunk and tributaries. Alongside other methods which estimate ice thickness, bed topography maps from IFPA have the potential to constrain projections of future sea-level rise, especially where radar data are sparse.

Supraglacial debris cover regulates the melt rates of many glaciers in mountainous regions around the world, thereby modifying the availability and quality of downstream water resources. However, the influence of supraglacial debris is often poorly represented within glaciological models, due to the absence of a technique to provide high-precision, spatially continuous measurements of debris thickness. Here, we use high-resolution UAV-derived thermal imagery, in conjunction with local meteorological data, visible UAV imagery and vertically profiled debris temperature time series, to model the spatially distributed debris thickness across a portion of Llaca Glacier in the Cordillera Blanca of Peru. Based on our results, we simulate daily sub-debris melt rates over a 3-month period during 2019. We demonstrate that, by effectively calibrating the radiometric thermal imagery and accounting for temporal and spatial variations in meteorological variables during UAV surveys, thermal UAV data can be used to more precisely represent the highly heterogeneous patterns of debris thickness and sub-debris melt on debris-covered glaciers. Additionally, our results indicate a mean sub-debris melt rate nearly three times greater than the mean melt rate simulated from satellite-derived debris thicknesses, emphasising the importance of acquiring further high-precision debris thickness data for the purposes of investigating glacier-scale melt processes, calibrating regional melt models and improving the accuracy of runoff predictions.

In this study, we compare equilibrium-line altitudes (ELAs) calculated using the area–altitude balance ratio (AABR) and the accumulation–area ratio (AAR) methods, with measured ELAs derived from direct field observations. We utilise a GIS toolbox to calculate the ELA for 64 extant glaciers by applying the AABR and AAR methods to DEMs and polygons of their geometry. The calculated ELAs (c-ELAs) are then compared to measured zero-net balance ELAs (znb-ELAs) obtained from mass-balance time series held by the WGMS for the same glaciers. The correlation between znb-ELAs and AABR (1.56)/AAR (0.58) c-ELAs is very strong, with an r2 = 0.99. The smallest median difference between znb-ELAs and c-ELAs (i.e. 65.5 m) is obtained when a globally representative AABR of 1.56 is used. When applied to palaeoglacier-climate applications, this difference translates to ~0.42°C, well within the uncertainty of palaeotemperature proxies used to determine mean summer temperature at the ELA. The more widely used mean AABR of 1.75 is shown to be statistically invalid due to the skewness of the dataset. On this basis, when calculating glacier ELAs, we recommend the use of a global AABR value of 1.56.

Radar surveys across ice sheets typically measure numerous englacial layers that can often be regarded as isochrones. Such layers are valuable for extrapolating age–depth relationships away from ice-core locations, reconstructing palaeoaccumulation variability, and investigating past ice-sheet dynamics. However, the use of englacial layers in Antarctica has been hampered by underdeveloped techniques for characterising layer continuity and geometry over large distances, with techniques developed independently and little opportunity for inter-comparison of results. In this paper, we present a methodology to assess the performance of automated layer-tracking and layer-dip-estimation algorithms through their ability to propagate a correct age–depth model. We use this to assess isochrone-tracking techniques applied to two test case datasets, selected from CreSIS MCoRDS data over Antarctica from a range of environments including low-dip, continuous layers and layers with terminations. We find that dip-estimation techniques are generally successful in tracking englacial dip but break down in the upper and lower regions of the ice sheet. The results of testing two previously published layer-tracking algorithms show that further development is required to attain a good constraint of age–depth relationships away from dated ice cores. We recommend that auto-tracking techniques focus on improved linking of picked stratigraphy across signal disruptions to enable accurate determination of the Antarctic-wide age–depth structure.

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.

The basic theory of crevasse formation suggests that crevasses initiate at or near the surface. However, due to variations in stress with depth, it has been suggested that it is possible for crevasses to initiate at depths of 10–30m. From December 2006 to January 2007, hot-water drilling on Pine Island Glacier, West Antarctica, was found to trigger crevasses. Satellite imagery and field investigations in 2008, including ice cores, radar and GPS, revealed that these formed a new band of arcuate (curvilinear) crevasses around 70 km long and 100 m deep. This new band is located 10 km upstream from the previous limit of the arcuate crevasse zone. The crevasses were triggered on drilling through an exceptional ice layer at >20m depth. Ice layers within the firn will change both the strength and stress intensity. As the firn changes spatially and temporally (e.g. with the burial of an ice layer), it is possible for the position of crevasse initiation to change whilst the along-stream strain-rate profile remains constant. However, the main cause of an upstream migration of the arcuate crevasse zone on Pine Island Glacier is still likely to be an increase in strain rate.

Measurements of surface dynamics on polythermal John Evans Glacier, Nunavut, Canada, over two winter periods and every 7–10 days throughout two melt seasons (June–July 2000, 2001) provide new insight into spatio-temporal patterns of High Arctic glacier dynamics. In the lower ablation zone, mean annual surface velocities are 10–21 m a–1, but peak velocities up to 50% higher are attained during late June/early July. In the upper ablation zone and lower accumulation zone, mean annual surface velocities are typically 10–18 m a–1, and peak velocities up to 40% higher occur during late July. In the upper accumulation zone, mean annual surface velocities are 2–9 m a–1, and motion in mid- to late July exceeds this by up to 10%. Rapid drainage of ponded supraglacial water in the upper ablation zone to an initially distributed subglacial drainage system in mid-June may force excess surface motion in the warm-based lower glacier. The data indicate that the duration of the velocity response may be related to the rate of channelization of the basal drainage, and the velocity response may be transmitted up-glacier by longitudinal coupling. An increase in surface velocities in the middle glacier in late July occurs in conjunction with the opening of two further moulins in the accumulation zone.

The longitudinal pattern of surface velocity of a large, predominantly cold, polythermal glacier (John Evans Glacier, Ellesmere Island, Canada) was measured over summer and winter periods. In the accumulation and upper ablation areas, where ice is predominantly cold-based, summer velocities were slightly higher than overwinter velocities. Predicted velocities due to ice deformation alone in these areas closely matched these observations in the winter, with limited basal motion likely in the summer. In the lower ablation area, where ice is likely warm-based, measured summer velocities were up to double overwinter velocities. Predicted ice deformation could not account for all of these measured velocities in either summer or winter. This suggests that basal motion occurs throughout the year over at least part of the lower ablation area. This finding is supported by radio-echo sounding, subglacial drainage reconstructions and analyses of early-summer meltwater chemistry, which suggest that subglacial water is present throughout the year in this region. In summer, basal motion may account for up to 75% of the total surface velocity throughout the lower ablation area. The inferred rate of basal motion increases sharply directly below a set of moulins by which most surface meltwater reaches the glacier bed.

Dye-tracer experiments undertaken over two summer melt seasons at polythermal John Evans Glacier, Ellesmere Island, Canada, were designed to investigate the character of the subglacial drainage system and its evolution over a melt season. In both summers, dye injections were conducted at several moulins and traced to a single subglacial outflow. Tracer breakthrough curves suggest that supraglacial meltwater initially encounters a distributed subglacial drainage system in late June. The subsequent development and maintenance of a channelled subglacial network are dependent upon sustained high rates of surface melting maintaining high supraglacial inputs. In a consistently warm summer (2000), subglacial drainage became rapidly and persistently channelled. In a cooler summer (2001), distributed subglacial drainage predominated. These observations confirm that supraglacial meltwater can access the bed of a High Arctic glacier in summer, and induce significant structural evolution of the subglacial drainage system. They do not support the view that subglacial drainage systems beneath polythermal glaciers are always poorly developed. They do suggest that the effects on ice flow of surface water penetration to the bed of predominantly cold glaciers may be short-lived.

Previous studies using Fourier transformation (FT) methods to analyze subglacial roughness have shown promise for distinguishing between different types of subglacial landscape from raw subglacial elevation data. We derive a two-parameter FT roughness index {ξ, η}, where • is based on the FT of elevation (as previously considered in isolation), and η is based on both the FT of elevation and the FT of bed-slope profile. In this way, we take account of both vertical and horizontal irregularities in subglacial surfaces. We demonstrate the statistical veracity of using {ξ, η} to consider roughness in terms of obstacle amplitudes and spacing, and consider the use of {ξ, η} in studies of ice dynamics and subglacial geomorphological interpretation. We show that {ξ, η} can be linked to basal sliding rates on the metre scale, and can be used to differentiate further than single-parameter roughness indices between different classes of subglacial landscape, in particular between erosional and depositional settings.

Isolated areas of high basal drag, or ‘sticky spots’, are important and poorly understood features in the force balance and dynamics of West Antarctic ice streams. Characterizing sticky spots formed by thin or drying subglacial till using ice-penetrating radar is theoretically possible, as high radar bed-returned power (BRP) is commonly related to an abundance of free water at the ice/bed interface, provided losses from englacial attenuation can be estimated. In this study we use airborne radar data collected over Evans Ice Stream to extract BRP profiles and test the sensitivity of BRP to the adopted englacial attenuation correction. We analyse 11 ~ 2 0 km profiles in four fast-flow areas where sticky spots have been inferred to exist on the basis of model and surface data inversions. In the majority of profiles we note that the increase in basal drag is accompanied by a decrease in BRP and suggest that this is evidence both for the presence of a sticky spot in those locations and that local variations in subglacial hydrology are responsible for their existence. A comparison is made between empirical and numerical modelling approaches for deriving englacial attenuation, and our findings generally support previous studies that advocate a modelling approach.

The potential for future dynamical instability of Pine Island Glacier, West Antarctica, has been addressed in a number of studies, but information on its past remains limited. In this study we use airborne radio-echo sounding (RES) data acquired over Pine Island Glacier to investigate past variations in accumulation pattern. In the dataset a distinctive pattern of layers was identified in the central part of the glacier basin. We use these layers as chronological identifiers in order to construct elevation maps of the internal stratigraphy. The observed internal layer stratigraphy is then compared to calculated stratigraphy from a three-dimensional ice-flow model that has been forced with different accumulation scenarios. The model results indicate that the accumulation pattern is likely to have changed at least twice since the deposition of the deepest identified layer. Additional RES data linked to the Byrd ice core provide an approximate timescale. This timescale suggests that the layers were deposited at the beginning of or during the Holocene period. Thus the widespread changes occurring in the coastal extent of the West Antarctic ice sheet at the end of the last glacial period could have been accompanied by changes in accumulation pattern.

In 1976, David Sugden and Brian John developed a classification for Antarctic landscapes of glacial erosion based upon exposed and eroded coastal topography, providing insight into the past glacial dynamics of the Antarctic ice sheets. We extend this classification to cover the continental interior of Antarctica by analysing the hypsometry of the subglacial landscape using a recently released dataset of bed topography (BEDMAP2). We used the existing classification as a basis for first developing a low-resolution description of landscape evolution under the ice sheet before building a more detailed classification of patterns of glacial erosion. Our key finding is that a more widespread distribution of ancient, preserved alpine landscapes may survive beneath the Antarctic ice sheets than has been previously recognized. Furthermore, the findings suggest that landscapes of selective erosion exist further inland than might be expected, and may reflect the presence of thinner, less extensive ice in the past. Much of the selective nature of erosion may be controlled by pre-glacial topography, and especially by the large-scale tectonic structure and fluvial valley network. The hypotheses of landscape evolution presented here can be tested by future surveys of the Antarctic ice sheet bed.

Flood-carved landforms across the deglaciated terrain of Victoria Land, East Antarctica, provide convincing geomorphological evidence for the existence of subglacial drainage networks beneath the Antarctic ice sheet, and motivate research into the inaccessible environment beneath the contemporary ice sheet. Through this research, our understanding of Antarctic subglacial hydrology is steadily building, and this paper presents an overview of the current state of knowledge. The conceptualization of subglacial hydrological behaviour was developed at temperate and Arctic glaciers, and is thus less mature in the Antarctic. Geophysical and remote sensing observations have demonstrated that many subglacial lakes form part of a highly dynamic network of subglacial drainage beneath the Antarctic ice sheet. Recent research into subglacial water flows, other than those directly concerned with lakes, has discovered potentially significant impacts on ice stream dynamics, ice sheet mass balance, and supplies of water to the ocean potentially affecting circulation and nutrient productivity. Despite considerable advances in understanding there remain a number of grand challenges that must be overcome in order to improve our knowledge of these subglacial hydrological processes.