There is a growing collection of ice-penetrating radar systems being used in a variety of glaciological settings to collect a body of data that is expanding in its scale, diversity and information richness. While there are nuanced differences among particular instruments, there are four main classes of radar systems used to observe the subsurface of ice sheets and glaciers: monotone pulse airborne sounders, chirped airborne sounders, ground-traversed impulsive systems and stationary continuous-wave systems. While these classes are often treated as fundamentally different instruments with distinct communities and literature, we seek to reframe them as choices on a continuum of potential waveforms. In particular, we focus on the tradeoff between achieving large coherent processing gains and mitigating interference from the matched-filter response of strong echoes at shorter delays as the key limiting factors for subsurface target detection. We also discuss the glaciological settings in which different waveform choices perform best, as well as potential avenues this raises for new instrument and experiment designs.

This paper introduces a methodology for detecting and classifying railway infrastructure elements using ground-penetrating radar (GPR). Seven trips, covering $57.4\,$km across a German track network were recorded. The network includes ballasted tracks, switches, bridges with dominant metal elements, and overpasses. GPR data from a multi-channel sensor array is transformed through a preprocessing pipeline to generate spatial B-scans. Statistical and structural features, including energy, depth, and skewness, are extracted via an overlapping sliding window. The detection algorithm operates as a streaming process, combining correlation-based pattern matching with feature-based thresholds in a state-machine architecture. The evaluation results demonstrate overall good performance for infrastructure detection and classification ($98$–$100\%$), with only minor misclassifications occurring for switch orientation estimation in scenarios involving closely spaced or physically interconnected switches ($82$–$100\%$).

Meltwater drainage through glaciers strongly influences ice dynamics and related hazards, yet detailed observations of active englacial and subglacial networks remain scarce due to challenges in direct observation. We present a novel high-density, UAV-based 3D ground-penetrating radar (GPR) survey of Switzerland’s Otemma Glacier, integrated with dye tracing experiments, photogrammetry and hydraulic potential modeling, to map internal water pathways with unprecedented spatial resolution and coverage. Advanced 3D imaging techniques, adapted from seismic diffraction processing, enhance the detection of small-scale conduits within the ice. We identify two primary subglacial channels having distinct hydraulic efficiencies and transport behaviors. Englacial drainage is found to be structurally complex, comprising lateral conduits, branching networks and localized pooling zones. These observations provide critical constraints on the geometry and dynamics of glacier drainage, supporting and extending existing models, and are essential for forecasting glacier behavior under future warming. UAV-based 3D GPR, in combination with novel processing strategies, emerges as an innovative approach for large-scale, repeatable surveys of glacier hydrology.

An alpine glacier below Sunlight Peak in northwest Wyoming was first photographically documented in 1893, near the end of the Little Ice Age and during the time of industrialization. Since then, evolving technologies have been applied to observe this glacier and nearby discontinuous permafrost for studies spanning Earth, environmental, and planetary sciences. Surveys in the 21st century indicate negative mass balance coinciding with rising average air temperature. This paper reviews the geological and geophysical data on record for the Sunlight Glacier system, presents new results from a 2023 fieldwork campaign combined with remote sensing analysis and comments on likely scenarios of future evolution for this individual body of ice within a broader alpine cryosphere feeding the watersheds of western North America.

The results of a ground-penetrating radar survey and multiproxy studies of the sediment cores collected from two lakes in the Valdai Highlands (East European Plain) provide new insights into the late glacial and Holocene environmental history of the region situated in the marginal zone of the last Scandinavian ice sheet. The cores were analyzed for organic carbon and nitrogen content, as well as for pollen and diatoms. The chronology of the cores is based on radiocarbon dates and pollen-based stratigraphy. The studied records document that vast dead ice masses and associated ice-dammed lakes existed in the Valdai Highlands area until ∼14 cal ka BP. Open tundra-steppe communities dominated the study area during the Oldest Dryas, Bølling, and Older Dryas (between ca. 17 and 14 cal ka BP), but dwarf birch (Betula nana), shrub alder (Alnus fruticosa), and willow (Salix) were also common. Scots pine forest (Pinus sylvestris) became common for a short interval during the Bølling warming (ca. 14.9 and 14.4 cal ka BP). The appearance of spruce (Picea) forest in the landscape occurred in the beginning of the Allerød warming (∼14 cal ka BP), but the open steppe-like plant communities remained common until the onset of the Holocene. The modern lake systems emerged at ∼10 cal ka BP, marked by an onset of organic-type sedimentation and the appearance of modern-type forests. The Mid-Holocene (∼8–4 cal ka BP) was the warmest time, as documented by the maximal distribution of temperate and broadleaved taxa in the region. The onset of agricultural land use and simultaneous trend of increasing lake trophic state and increasing paludification in the area is recorded at ∼2.5 cal ka BP.

Recent ice cores from the Allan Hills, a blue ice area in Antarctica, are nearly 3 million years old. These cores extend ice core chronologies, enabling new insight into key climate periods such as the Mid-Pleistocene Transition. The interpretation of these climate records is complex because of the disturbed stratigraphy in this ice. Here, we present a new three-dimensional multitrack electrical conductivity measurement method (3D ECM) to resolve layer structure. We demonstrate this technique on a cumulative 60 m of two large-diameter (241 mm) ice cores, ALHIC2201 and ALHIC2302. Measurements were taken on the upper section of both cores due to better ice core quality in this shallow ice. We find well-defined and dipping layering in both cores, averaging 29° in ALHIC2201 and 69° in ALHIC2302 from horizontal. We observe a slight decrease in dip with depth in both cores, although it only achieves statistical significance in ALHIC2302. We discuss how this new method can be applied to enable accurate, high-resolution multi-proxy record development even in ice cores with steeply dipping layers. 3D ECM improves interpretation of blue ice area cores by providing accurate, non-destructive constraints on stratigraphy.

Greenland’s peripheral glaciers and ice caps contribute disproportionately to sea-level rise relative to their small area. Winter snow accumulation directly influences glacier mass balance and downstream hydrology, but spatially extensive observations of this important mass balance component remain sparse. In this study, we present a unique multi-year (2008–2024) dataset of winter snow accumulation over A.P. Olsen Ice Cap, Northeast Greenland, from ground-penetrating radar surveys covering an average of 47 km per survey year. Our results reveal strong spatial heterogeneity that is likely influenced by wind redistribution and local topography, especially in the ablation zone. We compare our findings with automatic weather station data from three sites and outputs from the Copernicus Arctic Regional Reanalysis (CARRA). Governed by the high spatial variability, the automatic weather station point-based observations significantly underestimate regional accumulation by 40–45%. Despite the high spatial variability, the CARRA accumulated precipitation variable provides a reasonable overall mean winter snow accumulation (RMSE of 0.07 m w.e.); however, it fails to reproduce the complex non-linear relationship between snow depth and elevation observed in the radar data. Our findings emphasize the need for high-resolution, spatially extensive measurements to better understand snow accumulation on ice caps and glaciers and improve reanalysis assessments.

Quantifying snow water equivalent (SWE) with ground-penetrating radar (GPR) in a warming climate is complicated by the incidence and variability of liquid water in snow. Snow surveys conducted during the melt season serve as a valuable analog to conditions under future warming. Here, we determine the variability of wet snowpack properties (relative permittivity and density) to quantify their impact on SWE estimates using GPR. We collected spatially continuous snowpack measurements with 400 MHz GPR in 2012 and 2021 across repeat transects (∼150 km each year) along with spring and summer snow depth and density measurements from snow pits and snow cores. Snow relative permittivity values ranged from 2.06 to 2.62 in 2012 and 2.11 to 5.11 in 2021, resulting in calculated volumetric liquid water content (LWC) between 1.7% and 5.7% in 2012 and 2.1% and 16% in 2021. This variability in snow relative permittivity results in SWE uncertainties of 8% —33%, with more extreme cases reaching 13%—45%. We attribute this uncertainty to spatial and temporal variability in LWC when using GPR to estimate SWE. As snowpacks become wetter with rising atmospheric temperatures, GPR surveys should include in situ relative permittivity measurements to reduce depth and SWE interpretation uncertainties.

Peru hosts a significant portion of the world’s tropical glaciers, which are undergoing rapid mass loss due to climate change. Knowledge of the ice volume and bedrock topography of these glaciers is important for predicting changes in glacier dynamics, runoff, and interpreting ice-core records. This study presents results from glaciological and geophysical surveys conducted during a 2019 expedition to Nevado Huascarán, Peru’s highest mountain when four ice cores were extracted from the col and summit. Ground-penetrating radar measurements provided detailed ice thickness and snow accumulation data, highlighting complex internal glacier structure and indicated that the climatic records obtained from ice cores recovered in 2019 were continuous and extended past the Holocene. Ice flow modeling enabled investigation of glacier dynamics. It was shown that the upstream effect on ice-core record is minimal. Comparison with ice thickness modeling data for Huascarán from various sources revealed significant discrepancies with measured ice thicknesses, suggesting that the inversion methods underestimate ice thickness for the accumulation zones of mountain glaciers. This research contributes data for understanding glacier behavior in the context of climate change and for modeling efforts for better assessments of water resources, potential geohazards and paleoclimatic interpretations.

Recent observations have shown a fast decrease in thickness and area of Pyrenean glaciers in some cases leading to a stagnation of ice flow. However, their transition to a new paraglacial stage is not well understood. Through the combination of uncrewed aerial vehicles imagery, airborne LiDAR, ground-penetrating radar and ground temperature observations, we characterized the recent evolution of Infiernos Glacier. In 2021, this glacier had small sectors thicker than 25 m, but most of area did not exceed 10 m. The thickness losses from 2011 to 2023 reached 9 m in average, of which 5 m occurring during the period 2020–23. This trend demonstrates the significant ice melt under current climatic conditions. In the last years, the glacier has also shown a remarkable increase of debris cover extent. In these areas, the ice loss was reduced by half when compared to the thickness decrease in the entire glacier. Sub-freezing ground temperatures evidence the highly probable presence of permafrost or buried ice in the surroundings of the glacier. The clear signs of ice stagnation and the magnitude of area and thickness decrease support the main hypothesis of this work: After 2023, the Infiernos Glacier can no longer be considered a glacier and has become an ice patch.

Digital elevation models (DEMs) from the spaceborne interferometric radar mission TanDEM-X hold a large potential for glacier change assessments. However, a bias is potentially introduced through the penetration of the X-band signal into snow and firn. To improve our understanding of radar penetration on glaciers, we compare DEMs derived from the almost synchronous acquisition of TanDEM-X and Pléiades optical stereo-images of Grosser Aletschgletscher in March 2021. We found that the elevation bias – averaged per elevation bin – can reach up to 4–8 m in the accumulation area, depending on post co-registration corrections. Concurrent in situ measurements (ground-penetrating radar, snow cores, snow pits) reveal that the signal is not obstructed by the last summer horizon but reaches into perennial firn. Because of volume scattering, the TanDEM-X surface is determined by the scattering phase centre and does not coincide with a specific firn layer. We show that the bias corresponds to more than half of the decadal ice loss rate. To minimize the radar penetration bias, we recommend to select DEMs from the same time of the year and over long observation periods. A correction of the radar penetration bias is recommended, especially when combining optical and TanDEM-X DEMs.

We present a quantitative ground-penetrating radar (GPR) analysis of the Lower Calderone Glacieret to highlight the recent evolution of one of the southernmost glacial systems in Europe. The Upper and Lower Calderone Glacierets are the last two perennial ice bodies in the Apennine Mountains (Italy), and their continuous monitoring is important for glaciology, hydrology and climate science. We applied a proprietary auto-picking algorithm to track reflections accurately and objectively within three pseudo-3-D GPR data sets that were acquired over the Lower Glacieret in different years. After the time-to-depth conversion, the basal reflections were projected onto the normal versors of a 3-D topographic model of the survey area, at the different GPR trace positions. We then applied an Adjusted Inverse-Distance Weighted Spatial Interpolation method to extrapolate the ice-bedrock interface within the areas not directly covered by the GPR profiles and compare it with the topographic surface to recover the glacieret volume. In this paper, we critically examine the accuracy of the reconstructed models, including possible issues related to the challenging survey areas, such as local artifacts in the interpolated interface caused by irregular GPR coverage. We further discuss the various advantages of the implemented procedure with respect to more traditional glacier monitoring techniques.

This article describes a series of steps to integrate multiple modes of archaeological mapping in arid and agricultural settings. We use the coastal region of Peru as a case study and share our recent field experience at Cerro San Isidro, a multicomponent hill site located in the agriculture-intensive and mid-elevation (about 500 m asl) Moro region of the Nepeña Valley. In June and July 2022, we spent eight weeks deploying a combination of drone aerial imagery, pedestrian GPS reconnaissance, and GPR survey to map the surface and subsurface features at the site and in the adjacent agricultural fields. Our efforts suggest that the ancient settlement extended over an area of at least 50 ha, well beyond the visible surface architecture. Using a multimodal approach to confirming the partial destruction of archaeological vestiges by modern agricultural encroachment is both time-effective and noninvasive. The article offers insights from our experience, including the sequence of field operations, technical troubleshooting, and the collection and integration of datasets. We discuss the methodological potential and implications of this combination of multimodal mapping and its deployment in coastal Peru, a region that, like many others in the world, is increasingly subject to rapid agricultural expansion and other anthropogenic developments.

We review successes and challenges from five recent subglacial bedrock drilling campaigns intended to find evidence for Antarctic Ice Sheet retreat during warm periods in the geologic past. Insights into times when the polar ice sheets were smaller than present serve as guiding information for modeling efforts that aim to predict the rate and magnitude of future sea level rise that would accompany major retreat of the Antarctic Ice Sheet. One method to provide direct evidence for the timing of deglaciations and minimum extent of prior ice sheets is to extract subglacial bedrock cores for cosmogenic nuclide analysis from beneath the modern ice sheet surface. Here we summarize the lessons learned from five field seasons tasked with obtaining bedrock cores from shallow depths (<120 m beneath ice surface) across West Antarctica since 2016. We focus our findings on drilling efforts and technology and geophysical surveys with ground-penetrating radar. Shallow subglacial drilling provides a high risk, high reward means to test for past instabilities of the Antarctic Ice Sheet, and we highlight key challenges and solutions to increase the likelihood of success for future subglacial drilling efforts in polar regions.

Ice shelves, which regulate ice flow from the Antarctic ice sheet towards the ocean, are shaped by spatiotemporal patterns of surface accumulation, surface/basal melt and ice dynamics. Therefore, an ice dynamic and accumulation history are imprinted in the internal ice stratigraphy, which can be imaged by radar in the form of internal reflection horizons (IRHs). Here, IRHs were derived from radar data combined across radar platforms (airborne and ground-based) in coastal eastern Dronning Maud Land (East Antarctica), comprising three ice rises and adjacent two ice shelves. To facilitate interpretation of dominant spatiotemporal patterns of processes shaping the local IRH geometry, traced IRHs are classified into three different types (laterally continuous, discontinuous or absent/IRH-free). Near-surface laterally continuous IRHs reveal local accumulation patterns, reflecting the mean easterly wind direction, and correlate with surface slopes. Areas of current and past increased ice flow and internal deformation are marked by discontinuous or IRH-free zones, and can inform about paleo ice-stream dynamics. The established IRH datasets extend continent-wide mapping efforts of IRHs to an important and climatically sensitive ice marginal region of Antarctica and are ready for integration into ice-flow models to improve predictions of Antarctic ice drainage.

We analyze the internal structure of two polythermal glaciers, Hurd and Johnsons, located on Livingston Island, Antarctica, using 200 and 750 MHz GPR data collected in 2003/04, 2008/09 and 2016/17 field campaigns. Based on the different permittivities of snow and ice, we determined the thickness distribution of the end-of winter snow cover and of the cold ice layer. Their knowledge is fundamental for mass balance and glacier dynamics studies due to the different densities and rheological properties of such media. The average measured thicknesses for the snow and cold ice layers (the latter including the snow layer) were of 1.44 ± 0.09 and 29.1 ± 1.5 m, and their corresponding maxima were of 2.45 ± 0.21 and 80.8 ± 2.5 m. GPR snow profiling allowed for extension of the coverage of the snow thickness survey, but added little information to that supplied by snow pits, stake readings and manual snow probing, because of the multiplicity of reflections within the seasonal snowpack caused by internal ice layers and lenses. The polythermal structure determined for Hurd Glacier fits into the so-called Scandinavian type, seldom reported for the Antarctic region.

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.

Ground penetrating radar (GPR) has been extensively used in glaciology to infer glacier's ice thickness, liquid water content, water drainage pathways, and other properties. The interpretation of such GPR data is not always straightforward and for temperate glaciers, the signal is often affected by strong scattering and attenuation. It has often been suggested that such effects originate from englacial water inclusions, since water and ice have a large contrast in their di-electric permittivity. To investigate such effects quantitatively, we perform an extensive numerical modeling study of GPR signals. By exploring how different liquid water contents (LWC) and water-inclusions size affect the GPR signal, we show that their effects are much larger than the potential presence of a wet snowpack or a heterogeneous distribution of ice permittivity. In particularly, we show that the presence of such water inclusions is a necessary and sufficient condition for reproducing the typical characteristics of GPR data acquired in the field. Further, we find that for 25 MHz GPR antennas, a bulk LWC $\gtrsim$ 0.2%, associated with decimeters-scale water inclusions already limits bedrock detectability for ice thicknesses $\gtrsim 100$ m. Since these values are typical for Alpine glaciers, they clarify why the quality of GPR data is often poor in such environments.

Measuring surface mass-balance in the accumulation areas of glaciers is challenging because of the high spatial variability of snow accumulation and the difficulty of conducting annual field glaciological measurements. Here, we propose a method that can solve both these problems for many locations. Ground-penetrating radar measurements and firn cores extracted from a site in the French Alps were first used to reconstruct the topography of a buried end-of-summer snow horizon from a past year. Using these data and surface elevation observations from LiDAR and Global Navigation Satellite System instruments, we calculated the submergence velocities over the period between the buried horizon and more recent surface elevation observations. The differences between the changes in surface elevation and the submergence velocities were then used to calculate the annual surface mass-balances with an accuracy of ±0.34 m w.e. Assuming that the submergence velocities remain stable over several years, the surface mass-balance can be reconstructed for subsequent years from the differences in surface elevation alone. As opposed to the glaciological method that requires substantial fieldwork year after year to provide only point observations, this method, once submergence velocities have been calculated, requires only remote-sensing data to provide spatially distributed annual mass-balances in accumulation areas.

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)).