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.

The increase in iceberg discharge into the polar oceans highlights the importance of understanding how quickly icebergs are deteriorating and where the resulting freshwater injection is occurring. Recent advances in quantifying iceberg deterioration through combinations of modeling, remote sensing and direct in situ measurements have successfully calculated overall ablation rates, and surface and sidewall ablation; however, in situ measurements of basal melt rates have been difficult to obtain. Radar has successfully measured iceberg thickness, but repeat measurements, which would capture a change in iceberg thickness with time, have not yet been collected. Here we test the applicability of using an on-iceberg autonomous phase-sensitive radar (ApRES) to quantify basal ablation rates of a large (~800 m long) non-tabular Arctic iceberg during an intensive 2019 summer field campaign in Sermilik Fjord, southeast Greenland. We find that ApRES can be used to measure basal ablation even over a short deployment period (10 d), and also provide a lower bound on sidewall melt. This study fills a critical gap in iceberg research and pushes the limits of field instrumentation.

Ruth Glacier is situated in the Central Alaska Range, with the Don Sheldon Amphitheater comprising much of its broad accumulation area, directly adjacent to North America's tallest mountain, Denali. From there it funnels through the ‘Great Gorge,’ flanked by steep valley walls reaching over 1500 m. We combine airborne and ground-based radar measurements of ice thickness with satellite-derived surface velocities to constrain ice flux above and below the gorge, and employ a mass conservation approach to estimate the glacier's thickness within the gorge. We measure ice thickness in the amphitheater to reach 950 m, and estimate centerline thickness in the gorge to range from 610 to 960 m. Our estimates are up to two times greater than those suggested by global models, and allow us to confirm that the Great Gorge rivals Hells Canyon as the deepest gorge in North America. We found that the geometry of the gorge prevents radar measurements of ice thickness there since returns from the subglacial valley walls would precede and potentially occlude nadir bed returns. The same may be true of other unmapped mountain glaciers; however, thickness may be determined using appropriately located flux gates where radar sounding is feasible, combined with mass conservation methods.

Sub-glacial canyon features up to 580 m deep between flat terraces were identified beneath Devon Ice Cap during a 2023 radar echo sounding (RES) survey. The largest canyon connects a hypothesized brine network near the Devon Ice Cap summit with the marine-terminating Sverdrup outlet glacier. This canyon represents a probable drainage route for the hypothesized water system. Radar bed reflectivity is consistently 30 dB lower along the canyon floor than on the terraces, contradicting the signature expected for sub-glacial water. We compare these data with backscattering simulations to demonstrate that the reflectivity pattern may be topographically induced. Our simulated results indicated a 10 m wide canal-like water feature is unlikely along the canyon floor, but smaller features may be difficult to detect via RES. We calculated basal temperature profiles using a 2D finite difference method and found the floor may be up to 18°C warmer than the terraces. However, temperatures remain below the pressure melting point, and there is limited evidence that the canyon floor supports a connected drainage system between the DIC summit and Sverdrup Glacier. The terrain beneath Devon Ice Cap demonstrates limitations for RES. Future studies should evaluate additional correction methods near complex terrain, such as RES simulation as we demonstrate here.

Ice sheet models use observations to infer basal shear stress, but the variety of methods and datasets available has resulted in a wide range of estimates. Radar-based metrics such as reflectivity and specularity content have been used to characterize subglacial hydrologic conditions that are linked to spatial variations in basal shear stress. We explore whether radar metrics can be used to inform models about basal shear stress. At Thwaites Glacier, West Antarctica, we sample basal shear stress inversions across a wide range of ice sheet models to see how the basal shear stress distribution changes in regions of varying relative reflectivity and specularity content. Our results reveal three key findings: (1) Regions of high specularity content exhibit lower mean basal shear stresses (2) Wet and bumpy regions, as characterized by high relative reflectivity and low specularity content, exhibit higher mean basal shear stresses (3) Models disagree about what basal shear stress should be at the onset of rapid ice flow and high basal melt where relative reflectivity and specularity content are low.

Observations remain sparse for peripheral glaciers and ice caps in Greenland. Here, we present the results of a multi-frequency radar survey of Lyngmarksbræen Ice Cap in West Greenland conducted in April 2017. Radar measurements show thick ice of up to ~120 m in subglacial valleys associated with the largest outlet glaciers, while relatively thin ice cover the upper plateau ice divides, suggesting future vulnerability to ice cap fragmentation. At the time of the radar survey, Lyngmarksbræen Ice Cap had a total volume of 0.82 ± 0.1 km3. Measurements show a 1.5–2 m thick end-of-winter snowpack, and that firn is largely absent, signifying a prolonged period of negative mass balance for most of the ice cap. The thermal regime of Lyngmarksbræen Ice Cap is investigated through analysis of scattering observed along radar profiles. Results show that the ice cap is largely below the pressure melting point, but that temperate ice exists both in deep basal pockets and in shallow zones that some places extend from ~15 m depth and to the ice base. The distribution of shallow temperate ice appears unrelated to variations in ice thickness; instead we find a strong correlation to the presence of nearby surface crevasses.

We recently reported on the radio-frequency attenuation length of cold polar ice at Summit Station, Greenland, based on bi-static radar measurements of radio-frequency bedrock echo strengths taken during the summer of 2021. Those data also allow studies of (a) the relative contributions of coherent (such as discrete internal conducting layers with sub-centimeter transverse scale) vs incoherent (e.g. bulk volumetric) scattering, (b) the magnitude of internal layer reflection coefficients, (c) limits on signal propagation velocity asymmetries (‘birefringence’) and (d) limits on signal dispersion in-ice over a bandwidth of ~100 MHz. We find that (1) attenuation lengths approach 1 km in our band, (2) after averaging 10 000 echo triggers, reflected signals observable over the thermal floor (to depths of ~1500 m) are consistent with being entirely coherent, (3) internal layer reflectivities are ≈–60$\to$–70 dB, (4) birefringent effects for vertically propagating signals are smaller by an order of magnitude relative to South Pole and (5) within our experimental limits, glacial ice is non-dispersive over the frequency band relevant for neutrino detection experiments.

Estimates of ice shelf mass loss are typically based on surface height measurements, assuming hydrostatic equilibrium and estimated firn thickness. Recent investigations, however, challenge the assumption that ice shelves are freely floating, particularly in proximity to narrow structures such as basal channels and shear margins. We compare contemporaneous measurements of Antarctic ice shelf thickness, from ice-penetrating radar, to freeboard height, from laser altimetry, acquired during multiple airborne surveys. On average, the hydrostatic thickness differs from observed thickness by at least ~17 ± 98 m, but this difference varies well beyond the propagated error within and among ice shelves, and depends on the corrections applied. We find that uncertainty in firn thickness can account for most, but not all, of the imbalance. Overall, errors in hydrostatic thickness do not significantly impact estimated basal melt rates. Our results indicate that localized approaches to estimating ice shelf thickness and rates of change are not applicable at large scales, and vice versa, and point to the need for more abundant and accurate firn and ice thickness measurements to improve estimates and predictions of ice shelf mass loss.

Over the past 70 years, many different components of the cryosphere have been imaged with a variety of radar systems using increasingly sophisticated processing techniques. These systems use various pulse lengths, signal frequencies and, in some cases, modulated signals. The increasing diversity of radar systems has created the potential for confusion due to the use of non-consistent terminology. Here we provide an overview of state-of-the-science radar technologies and suggest a simplified and unified terminology for use by the cryosphere community. We recommend a terminology that is target independent but specifies the characteristics of the signal. Following this recommendation, commercial impulse systems that penetrate the subsurface should be referred to as ground-penetrating radar (GPR), and pulse radars as radio-echo sounding (RES). Continuous-wave (CW) radar systems should be referred to as ground-penetrating CW radars. We further suggest any additional characterisation of the system be expressed using descriptors that specify the platform it is mounted on (e.g. airborne) or the frequency range (e.g. HF (high frequency)) or modulation (e.g. FM (frequency modulated)).

Ice-penetrating radar sounding is a powerful geophysical tool for studying terrestrial and planetary ice with a rich glaciological heritage reaching back over half a century. Recent years have also seen rapid growth in both the radioglaciological community itself and in the scope and sophistication of its analysis of ice-penetrating radar data. This has been spurred by a combination of growing datasets and improvements in computational resources as well as advances in radar sounding instrumentation and platforms. Together, these developments are transforming the field and highlight exciting paths forward for future innovation and investigation.

We use the apparent horizontal shift of an englacial tephra layer outcrop to calculate local glacier mass balance on Sléttjökull, a lobe of Mýrdalsjökull in Southern Iceland. For this approach, the dipping angle of the englacial tephra layer in the glacier upstream of the outcrop and the flow velocity of the ice need to be known. An earlier investigation was expanded by the application of ground-penetrating radar, detecting the depth of the tephra along tracks with a total length of 10 km. Interpolation between the tracks enables us to derive the dipping angle of the layer along several flow lines. Together with glacier surface velocities, determined from feature tracking, we are able to estimate the local surface mass balance from the horizontal displacement of the tephra outcrop using freely available satellite imagery without additional fieldwork. The earlier local balance series was extended to the period 2014/15 to 2019/20. Although the results for the individual profiles differ slightly from each other, they show the same temporal pattern and clear variations from year to year. The results are compared to traditional mass-balance data from Hofsjökull. The two series show a good agreement in their interannual variability.

Ice streams are warmed by shear strain, both vertical shear near the bed and lateral shear at the margins. Warm ice deforms more easily, establishing a positive feedback loop in an ice stream where fast flow leads to warm ice and then to even faster flow. Here, we use radar attenuation measurements to show that the Siple Coast ice streams are colder than previously thought, which we hypothesize is due to along-flow advection of cold ice from upstream. We interpret the attenuation results within the context of previous ice-temperature measurements from nearby sites where hot-water boreholes were drilled. These in-situ temperatures are notably colder than model predictions, both in the ice streams and in an ice-stream shear margin. We then model ice temperature using a 1.5-dimensional numerical model which includes a parameterization for along-flow advection. Compared to analytical solutions, we find depth-averaged temperatures that are colder by 0.7°C in the Bindschadler Ice Stream, 2.7°C in the Kamb Ice Stream and 6.2–8.2°C in the Dragon Shear Margin of Whillans Ice Stream, closer to the borehole measurements at all locations. Modelled cooling corresponds to shear-margin thermal strengthening by 3–3.5 times compared to the warm-ice case, which must be compensated by some other weakening mechanism such as material damage or ice-crystal fabric anisotropy.

Hercules Dome, Antarctica, has long been identified as a prospective deep ice core site due to the undisturbed internal layering, climatic setting and potential to obtain proxy records from the Last Interglacial (LIG) period when the West Antarctic ice sheet may have collapsed. We performed a geophysical survey using multiple ice-penetrating radar systems to identify potential locations for a deep ice core at Hercules Dome. The surface topography, as revealed with recent satellite observations, is more complex than previously recognized. The most prominent dome, which we term ‘West Dome’, is the most promising region for a deep ice core for the following reasons: (1) bed-conformal radar reflections indicate minimal layer disturbance and extend to within tens of meters of the ice bottom; (2) the bed is likely frozen, as evidenced by both the shape of the measured vertical ice velocity profiles beneath the divide and modeled ice temperature using three remotely sensed estimates of geothermal flux and (3) models of layer thinning have 132 ka old ice at 45–90 m above the bed with an annual layer thickness of ~1 mm, satisfying the resolution and preservation needed for detailed analysis of the LIG period.

Over the last 25 years, radiowave detection of neutrino-generated signals, using cold polar ice as the neutrino target, has emerged as perhaps the most promising technique for detection of extragalactic ultra-high energy neutrinos (corresponding to neutrino energies in excess of 0.01 Joules, or 1017 electron volts). During the summer of 2021 and in tandem with the initial deployment of the Radio Neutrino Observatory in Greenland (RNO-G), we conducted radioglaciological measurements at Summit Station, Greenland to refine our understanding of the ice target. We report the result of one such measurement, the radio-frequency electric field attenuation length $L_\alpha$. We find an approximately linear dependence of $L_\alpha$ on frequency with the best fit of the average field attenuation for the upper 1500 m of ice: $\langle L_\alpha \rangle = ( ( 1154 \pm 121) - ( 0.81 \pm 0.14) \, ( \nu /{\rm MHz}) ) \,{\rm m}$ for frequencies ν ∈ [145 − 350] MHz.

During the past 20 years, multi-channel radar emerged as a key tool for deciphering an ice sheet's internal architecture. To assign ages to radar reflections and connect them over large areas in the ice sheet, the layer genesis has to be understood on a microphysical scale. Synthetic radar trace modelling based on the dielectric profile of ice cores allows for the assignation of observed physical properties’ variations on the decimetre scale to radar reflectors extending from the coring site to a regional or even whole-ice-sheet scale. In this paper we rely on the available dielectric profiling data of the northern Greenland deep ice cores: NGRIP, NEEM and EGRIP. The three records are well suited for assigning an age model to the stratigraphic radar-mapped layers, and linking up the reflector properties to observations in the cores. Our modelling results show that the internal reflections are mainly due to conductivity changes. Furthermore, we deduce fabric characteristics at the EGRIP drill site from two-way-travel-time differences of along and across-flow polarized radarwave reflections of selected horizons (below 980 m). These indicate in deeper parts of the ice column an across-flow concentrated c-axis fabric.

The earliest airborne geophysical campaigns over Antarctica and Greenland in the 1960s and 1970s collected ice penetrating radar data on 35 mm optical film. Early subglacial topographic and englacial stratigraphic analyses of these data were foundational to the field of radioglaciology. Recent efforts to digitize and release these data have resulted in geometric and ice-thickness analysis that constrain subsurface change over multiple decades but stop short of radiometric interpretation. The primary challenge for radiometric analysis is the poorly-characterized compression applied to Z-scope records and the sparse sampling of A-scope records. Here, we demonstrate the information richness and radiometric interpretability of Z-scope records. Z-scope pixels have uncalibrated fast-time, slow-time, and intensity scales. We develop approaches for mapping each of these scales to physical units (microseconds, seconds, and signal to noise ratio). We then demonstrate the application of this calibration and analysis approach to a flight in the interior of East Antarctica with subglacial lakes and to a reflight of an East Antarctic ice shelf that was observed by both archival and modern radar. These results demonstrate the potential use of Z-scope signals to extend the baseline of radiometric observations of the subsurface by decades.

Firn compaction models inform mass-balance estimates and paleo-climate reconstructions, but current models introduce key uncertainties. For example, models disagree on the dependence of density and compaction on accumulation rate. Observations of compaction to test these models are rare, partly because in situ methods for measuring englacial strain are time-consuming and expensive. Moreover, shallow measurements may confound strain due to compaction with strain due to ice-sheet flow. We show that phase-sensitive radio-echo sounder (pRES) systems, typically deployed to measure sub-shelf melting or ice-sheet deformation, can be used to measure firn compaction and test firn models. We present two complementary methods for extracting compaction information from pRES data, along with a method for comparing compaction models to pRES observations. The methods make different assumptions about the density structure and vary in their need for independent density measurements. Compaction profiles computed from pRES data collected on three ice rises in West Antarctica are largely consistent with measured densities and a physics-based model. With their minimal logistic requirements, new pRES systems, such as autonomous pRES, could be inexpensively deployed to monitor firn compaction more widely. Existing phase-sensitive radar data may contain compaction information even when surveys targeted other processes.

The density–depth relationship of the Antarctic ice sheet is important for establishing a high-precision surface mass balance model and predicting future ice-sheet contributions to global sea levels. A new algorithm is used to reconstruct firn density and densification rate by inverting monostatic radio wave echoes from ground-operated frequency-modulated continuous wave radar data collected near four ice cores along the transect from Zhongshan Station to Dome A. The inverted density profile is consistent with the core data within 5.54% root mean square error. Due to snow redistribution, the densification rate within 88 km of ice core DT401 is correlated with the accumulation rate and varies greatly over horizontal distances of <5 km. That is, the depth at which a critical density of 830 kg m−3 is reached decreases and densification rate increases in high-accumulation regions but decreases in low-accumulation regions. This inversion technique can be used to analyse more Antarctic radar data and obtain the density distribution trend, which can improve the accuracy of mass-balance estimations.

Crary Ice Rise formed after the Ross Ice Shelf re-grounded ~1 kyr BP. We present new ice-penetrating radar data from two systems operating at center frequencies of 7 and 750 MHz that confirm the ice rise is composed of a former ice shelf buried by subsequent accumulation. Stacks of englacial diffraction hyperbolas are present almost everywhere across the central ice rise and extend up to ~350 m above the bed. In many cases, bed reflections beneath the diffraction hyperbolas are obscured for distances up to 1 km. Waveform modeling indicates that the diffraction hyperbolas are likely caused by marine ice deposits in former basal crevasses and rifts. The in-filling of rifts and basal crevasses may have strengthened the connection between the ice rise and the surrounding ice shelf, which could have influenced local and regional ice dynamics. Three internal reflection horizons mark the upper limit of disturbed ice and diffraction hyperbolas in different sections of the ice rise, indicating at least three stages of flow stabilization across the ice rise. A surface lineation visible in MODIS imagery corresponds spatially to deepening and strong deformation of these layers, consistent with the characteristics of former grounding lines observed elsewhere in Antarctica.

During the last few decades, bed-elevation profiles from radar sounders have been used to quantify bed roughness. Various methods have been employed, such as the ‘two-parameter’ technique that considers vertical and slope irregularities in topography, but they struggle to incorporate roughness at multiple spatial scales leading to a breakdown in their depiction of bed roughness where the relief is most complex. In this article, we describe a new algorithm, analogous to wavelet transformations, to quantify the bed roughness at multiple scales. The ‘Self-Adaptive Two-Parameter’ system calculates the roughness of a bed profile using a frequency-domain method, allowing the extraction of three characteristic factors: (1) slope, (2) skewness and (3) coefficient of variation. The multi-scale roughness is derived by weighted-summing of these frequency-related factors. We use idealized bed elevations to initially validate the algorithm, and then actual bed-elevation data are used to compare the new roughness index with other methods. We show the new technique is an effective tool for quantifying bed roughness from radar data, paving the way for improved continental-wide depictions of bed roughness and incorporation of this information into ice flow models.