Radio-echo Sounding Research Articles

Thickness of Ruth Glacier, Alaska, and depth of its Great Gorge from ice-penetrating radar and mass conservation

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

Entrained Water in Basal Ice Suppresses Radar Bed‐Echo Power at Active Subglacial Lakes

AbstractSubglacial lakes have been mapped across Antarctica with two methods, radio‐echo sounding (RES) and ice‐surface deformation. At sites where both are coincident, these methods typically provide conflicting interpretations about the ice‐bed interface. With a single exception, active subglacial lakes identified by surface deformation do not display the expected flat, bright, and specular bed reflection in RES data, characteristic of non‐active lakes. This observational conundrum suggests that our understanding of Antarctic subglacial hydrology, especially beneath important fast‐moving ice streams, remains incomplete. Here, we use an airborne RES campaign that surveyed a well‐characterized group of active subglacial lakes on lower Mercer and Whillans ice streams, West Antarctica, to explore inconsistency between the two observational techniques. We test hypotheses of increased scattering and attenuation due to the presence of an active subglacial lake system that could suppress reflected bed‐echo power for RES observations in these locations, finding that entrained water is most plausible.

Sentinel-1 detection of ice slabs on the Greenland Ice Sheet

Abstract. Ice slabs are multi-meter-thick layers of refrozen ice that limit meltwater storage in firn, leading to enhanced surface runoff and ice sheet mass loss. To date, ice slabs have primarily been mapped using airborne ice-penetrating radar, which has limited spatial and temporal coverage. This makes it difficult to fully assess the current extent and continuity of ice slabs or to validate predictive models of ice slab evolution that are key to understanding their impact on Greenland's surface mass balance. Here, for the first time, we map the extent of ice slabs and superimposed ice facies across the entire Greenland Ice Sheet at 500 m resolution using dual-polarization Sentinel-1 (S-1) synthetic-aperture radar (SAR) data collected in winter 2016–2017. We do this by selecting empirical thresholds for the cross-polarized backscatter ratio and HV backscattered power that jointly optimize the agreement between airborne ice-penetrating radar data detections of ice slabs and the S-1 estimates of ice slab extent. Our results show that there is a strong correlation between C-band backscatter and the ice content of the upper ∼ 7 m of the firn column that enables ice slab mapping with S-1. Our new mapping shows that ice slabs are nearly continuous around the entire margin of the ice sheet. This includes regions in southwest Greenland where ice slabs have not been previously identified by ice-penetrating radar but where the S-1-inferred ice slab extent shows strong agreement with the extent of visible runoff mapped from optical imagery. The algorithm developed here lays the groundwork for the long-term monitoring of ice slab expansion with current and future C-band satellite systems and highlights the potential added value of future L-band missions for near-surface studies in Greenland.

Alpine topography of the Gamburtsev Subglacial Mountains, Antarctica, mapped from ice sheet surface morphology

Abstract. Landscapes buried beneath the Antarctic Ice Sheet preserve information about the geologic and geomorphic evolution of the continent both before and during the wide-scale glaciation that began roughly 34×106 years ago. Since the inception of this ice sheet, some areas have remained cold-based and non-erosive, preserving ancient landscapes remarkably intact. The Gamburtsev Subglacial Mountains in central East Antarctica are one such landscape, maintaining evidence of tectonic, fluvial and glacial controls on their distinctly alpine morphology. The central Gamburtsev Mountains have previously been surveyed using airborne ice-penetrating radar; however, many questions remain as to their evolution and their influence on the East Antarctic Ice Sheet, including where in the region to drill for a 1.5×106 year-long “oldest-ice” core. Here, we derive new maps of the planform geometry of the Gamburtsev Subglacial Mountains from satellite remote sensing datasets of the ice sheet surface, based on the relationship between bed roughness and ice surface morphology. Automated and manual approaches to mapping were tested and validated against existing radar data and elevation models. Manual mapping was more effective than automated approaches at reproducing bed features observed in radar data, but a hybrid approach is suggested for future work. The maps produced here show the detail of mountain ridges and valleys on wavelengths significantly smaller than the spacing of existing radar flightlines, and mapping has extended well beyond the confines of existing radar surveys. Morphometric analysis of the mapped landscape reveals that it constitutes a preserved (>34 Ma) dendritic valley network, with some evidence for modification by topographically confined glaciation prior to ice sheet inception. The planform geometry of the landscape is a significant control on locations of basal melting, subglacial hydrological flows and the stability of the ice sheet over time, so the maps presented here may help to guide decisions about where to search for oldest ice.

Evaluating and locating a suitable bedrock drilling site near zhongshan station with airborne and ground-based observations

Understanding the interaction between the Antarctic biogeochemical cycle and global environmental change requires a comprehensive investigation of the ice sheet and underlying bedrock. An effective approach to investigating subglacial characteristics is through bedrock drilling, making the identification of a suitable drilling site a crucial prerequisite. In this study, data from airborne ice-penetrating radar and an Uncrewed Aerial Vehicle (UAV) remote sensing were used in combination to obtain real-time and high-resolution information about the surface, ice thickness, and subglacial conditions in a local area near Zhongshan Station, and further to identify a drilling site that is both safe and optimal for the successful acquisition of bedrock samples beneath the East Antarctic Ice Sheet. By analyzing the airborne ice-penetrating radar data from the 35th Chinese National Antarctic Research Expedition (CHINARE) (2018–2019), a potential drilling site was selected based on ice thickness (≤200 m) and a dry basal condition with no channels or basal crevasses, to avoid drilling fluid leakage. UAV data confirmed stable surface topography, ice flow, and the absence of surface crevasses. The optimal site was identified as located close to the western flank of the upstream part of the Dålk Glacier. In 2019, CHINARE successfully retrieved a valuable bedrock core sample from the borehole at this location. This study provides a feasible and effective methodology for locating a bedrock drilling site by using data from airborne ice-penetrating radar and UAV remote sensing.

Deep clustering in subglacial radar reflectance reveals subglacial lakes

Abstract. Ice-penetrating radar (IPR) imaging is a valuable tool for observing the internal structure and bottom of ice sheets. Subglacial water bodies, also known as subglacial lakes, generally appear as distinct, bright, flat, and continuous reflections in IPR images. In this study, we use available IPR images from the Gamburtsev Subglacial Mountains to extract one-dimensional reflector waveform features of the ice–bedrock interface. We apply a deep-learning method to reduce the dimension of the reflector features. An unsupervised clustering method is then used to separate different types of reflector features, including a reflector type corresponding to subglacial lakes. The derived clustering labels are then used to detect features of subglacial lakes in IPR images. Using this method, we compare the new detections with a known-lakes inventory. The results indicate that this new method identified additional subglacial lakes that were not previously detected, and some previously known lakes are found to correspond to other reflector clusters. This method can offer automatic detections of subglacial lakes and provide new insight for subglacial studies.

Ice Thickness Assessment of Non-Freshwater Lakes in the Qinghai–Tibetan Plateau Based on Unmanned Aerial Vehicle-Borne Ice-Penetrating Radar: A Case Study of Qinghai Lake and Gahai Lake

Ice thickness has a significant effect on the physical and biogeochemical processes of a lake, and it is an integral focus of research in the field of ice engineering. The Qinghai–Tibetan Plateau, known as the Third Pole of the world, contains numerous lakes. Compared with some information, such as the area, water level, and ice phenology of its lakes, the ice thickness of these lakes remains poorly understood. In this study, we used an unmanned aerial vehicle (UAV) with a 400/900 MHz ice-penetrating radar to detect the ice thickness of Qinghai Lake and Gahai Lake. Two observation fields were established on the western side of Qinghai Lake and Gahai Lake in January 2019 and January 2021, respectively. Based on the in situ ice thickness and the propagation time of the radar, the accuracy of the ice thickness measurements of these two non-freshwater lakes was comprehensively assessed. The results indicate that pre-processed echo images from the UAV-borne ice-penetrating radar identified non-freshwater lake ice, and we were thus able to accurately calculate the propagation time of radar waves through the ice. The average dielectric constants of Qinghai Lake and Gahai Lake were 4.3 and 4.6, respectively. This means that the speed of the radar waves that propagated through the ice of the non-freshwater lake was lower than that of the radio waves that propagated through the freshwater lake. The antenna frequency of the radar also had an impact on the accuracy of ice thickness modeling. The RMSEs were 0.034 m using the 400 MHz radar and 0.010 m using the 900 MHz radar. The radar with a higher antenna frequency was shown to provide greater accuracy in ice thickness monitoring, but the control of the UAV’s altitude and speed should be addressed.

Removal of Diagonal Strip Noise in Airborne Ice Penetrating Radar Data

Airborne ice-penetrating radar (IPR) is a powerful geophysical method for detecting ice thickness, subglacial topography and internal layers in the vast and frigid Antarctic ice sheet. Data processing is an essential step in improving the signal-to-noise ratio (SNR). In order to eliminate the diagonal strip noise mainly present in airborne IPR data of the Amery ice shelf and Grove Mountain regions, traditional ground penetrating radar processing methods such as F-K (frequency-wavenumber domain) filtering and F-K migration are applied. The validity of the data processing methods is proved by the quantitative index of image quality evaluation. It provides a reference for subsequent airborne radar data processing and application.

Ice and Snow Thickness of the IGAN Glacier in the Polar Urals from Ground-Based Radio-Echo Sounding in 2019 and 2021

Ice and Snow Thickness of the IGAN Glacier in the Polar Urals from Ground-Based Radio-Echo Sounding in 2019 and 2021

Airborne Radio-Echo Sounding Data Denoising Using Particle Swarm Optimization and Multivariate Variational Mode Decomposition

Radio-echo sounding (RES) is widely used for polar ice sheet detection due to its wide coverage and high efficiency. The multivariate variational mode decomposition (MVMD) algorithm for the processing of RES data is an improvement to the variational mode decomposition (VMD) algorithm. It processes data encompassing multiple channels. Determining the most effective component combination of the penalty parameter (α) and the number of intrinsic mode functions (IMFs) (K) is fundamental and affects the decomposition results. α and K in traditional MVMD are provided by subjective experience. We integrated the particle swarm optimization (PSO) algorithm to iteratively optimize these parameters—specifically, α and K—with high precision. This was then combined with the four quantitative parameters: energy entropy, signal-to-noise ratio (SNR), peak signal-to-noise ratio (PSNR), and root-mean-square error (RMSE). The RES signal decomposition results were judged, and the most effective component combination for noise suppression was selected. We processed the airborne RES data from the East Antarctic ice sheet using the combined PSO–MVMD method. The results confirmed the quality of the proposed method in attenuating the RES signal noise, enhancing the weak signal of the ice base, and improving the SNR. This combined PSO–MVMD method may help to enhance weak signals in deeper parts of ice sheets and may be an effective tool for RES data interpretation.

Stagnant ice and age modelling in the Dome C region, Antarctica

Abstract. The European Beyond EPICA project aims to extract a continuous ice core of up to 1.5 Ma, with a maximum age density of 20 kyr m−1 at Little Dome C (LDC). We present a 1D numerical model which calculates the age of the ice around Dome C. The model inverts for basal conditions and accounts either for melting or for a layer of stagnant ice above the bedrock. It is constrained by internal reflecting horizons traced in radargrams and dated using the EPICA Dome C (EDC) ice core age profile. We used three different radar datasets ranging from a 10 000 km2 airborne survey down to 5 km long ground-based radar transects over LDC. We find that stagnant ice exists in many places, including above the LDC relief where the new Beyond EPICA drill site (BELDC) is located. The modelled thickness of this layer of stagnant ice roughly corresponds to the thickness of the basal unit observed in one of the radar surveys and in the autonomous phase-sensitive radio-echo sounder (ApRES) dataset. At BELDC, the modelled stagnant ice thickness is 198±44 m and the modelled oldest age of ice is 1.45±0.16 Ma at a depth of 2494±30 m. This is very similar to all sites situated on the LDC relief, including that of the Million Year Ice Core project being conducted by the Australian Antarctic Division. The model was also applied to radar data in the area 10–15 km north of EDC (North Patch), where we find either a thin layer of stagnant ice (generally <60 m) or a negligible melt rate (<0.1 mm yr−1). The modelled maximum age at North Patch is over 2 Ma in most places, with ice at 1.5 Ma having a resolution of 9–12 kyr m−1, making it an exciting prospect for a future Oldest Ice drill site.

A comparison of contemporaneous airborne altimetry and ice-thickness measurements of Antarctic ice shelves

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

Spatial characterization of near-surface structure and meltwater runoff conditions across the Devon Ice Cap from dual-frequency radar reflectivity

Abstract. Melting and refreezing processes in the firn of the Devon Ice Cap control meltwater infiltration and runoff across the ice cap, but their full spatial extent and effect on near-surface structure is difficult to measure with surface-based traverses or existing satellite remote sensing. Here, we derive the coherent component of the near-surface return from airborne ice-penetrating radar surveys over the Devon Ice Cap, Canadian Arctic, to characterize firn containing centimeter- to meter-thick ice layers (i.e., ice slabs) formed from refrozen meltwater in firn. We assess the use of dual-frequency airborne ice-penetrating radar to characterize the spatial and vertical near-surface structure of the Devon Ice Cap by leveraging differences in range resolution of the radar systems. Comparison with reflectivities using a thin layer reflectivity model, informed by surface-based radar and firn core measurements, indicates that the coherent component is sensitive to the near-surface firn structure composed of quasi-specular ice and firn layers, limited by the bandwidth-constrained radar range resolution. Our results suggest that average ice slab thickness throughout the Devon Ice Cap percolation zone ranges from 4.2 to 5.6 m. This implies conditions that can enable lateral meltwater runoff and potentially contribute to the total surface runoff routed through supraglacial rivers down glacier. Together with the incoherent component of the surface return previously studied, our dual-frequency approach provides an alternative method for characterizing bulk firn properties, particularly where high-resolution radar data are not available.

Basal Melt Patterns around the Deep Ice Core Drilling Site in the Dome A Region from Ice-Penetrating Radar Measurements

Basal melt in the Dome A region will influence the deep-ice-core drilling at Kunlun Station. The melting point (wet bedrock) has a higher reflectivity than the surrounding area, which can be assessed using radar echoes from the bedrock. This paper uses a linear absorption model to determine wet and dry ice–bedrock interfaces around the Kunlun drilling site. In the determination process, an artificial intelligence model was applied to extract the ice–bedrock interface for inferring the ice thickness. Additionally, the various depth-averaged attenuation rates were used to identify the maximal range of basal melting. We mapped the patterns of the wet points on the bottom of the ice sheet and the modeled basal temperature to verify the results of the wet bed conditions. According to these maps of wet bed conditions, the areas with basal melting around the drilling site primarily appear in valley walls with low basal temperatures and are linked with hydraulic potential and bedrock elevation. Identifying the ice–bedrock interface is challenging in the valley bottom area, where the melting points are less than at the valley walls. Additionally, the melting proportions of 11.8% and 3.62% were calculated from two ice-penetrating radar data in this research. The mapped melting points around the site of Kunlun ice core drilling suggest complex ice flow effects and can be used to better interpret archives of old ice for paleoclimate research.

Changes in the Structure of the Snow Cover of Hansbreen (S Spitsbergen) Derived from Repeated High-Frequency Radio-Echo Sounding

This paper explores the potential of ground-penetrating radar (GPR) monitoring for an advanced understanding of snow cover processes and structure. For this purpose, the study uses the Hansbreen (SW Spitsbergen) records that are among the longest and the most comprehensive snow-cover GPR monitoring records available on Svalbard. While snow depth (HS) is frequently the only feature derived from high-frequency radio-echo sounding (RES), this study also offers an analysis of the physical characteristics (grain shape, size, hardness, and density) of the snow cover structure. We demonstrate that, based on GPR data (800 MHz) and a single snow pit, it is possible to extrapolate the detailed features of snow cover to the accumulation area. Field studies (snow pits and RES) were conducted at the end of selected accumulation seasons in the period 2008–2019, under dry snow conditions and HS close to the maximum. The paper shows that although the snow cover structure varies in space and from season to season, a single snow pit site can represent the entire center line of the accumulation zone. Numerous hard layers (HLs) (up to 30% of the snow column) were observed that reflect progressive climate change, but there is no trend in quantity, thickness, or percentage contribution in total snow depth in the study period. HLs with strong crystal bonds create a “framework” in the snowpack, which reduces compaction and, consequently, the ice formation layers slow down the rate of snowpack metamorphosis. The extrapolation of snow pit data through radar profiling is a novel solution that can improve spatial recognition of snow cover characteristics and the accuracy of calculation of snow water equivalent (SWE).

Comparison of Hydrothermal Structure of Two Glaciers in Spitsbergen and Tian Shan Based on Radio-Echo Sounding Data

Comparison of Hydrothermal Structure of Two Glaciers in Spitsbergen and Tian Shan Based on Radio-Echo Sounding Data

An overview of GPR subsurface exploration of planets and moons

Geophysical techniques were first tested beyond Earth during the Apollo program. Of those examined, radio-wave propagation methods appeared to be the most suitable for the moon and other solar system bodies. This was due to the electromagnetic characteristics of planetary subsurfaces and the possibility to remotely perform measurements on board spacecrafts and rovers. After the first successful experiment on the moon, more than 20 years passed before ground-penetrating radar (GPR) was included in the payload of a planetary mission. Technological advancements in GPR design and successful results of radio echo sounding measurements for the detection of basal water below terrestrial ice sheets paved the way for the application of similar techniques to search for liquid water in the Martian subsurface. Since deployment of the first two subsurface radar sounders above Mars, the number of proposed planetary missions relying on GPR for surveying the subsurface of planets, moons, and other objects has grown progressively. Six orbiting radar sounders and five GPRs mounted on rovers/landers have been employed so far to explore the moon, Mars, and comet 67P/GC. Some of these are in full operation and some are just starting to operate. Planned missions to the icy moons of Jupiter will also depend heavily on radar sounders to detect evidence of an internal ocean on Europa and to understand the habitability conditions on Europa, Ganymede, and Callisto. Finally, planetary missions to Earth's twin, the planet Venus, could take advantage of GPR to understand the cause of its drastic change in climatic conditions and the geologic phenomena that contributed to changing a watery and hospitable surface into a hot and asphyxiating inhabitable planet.

A site for deep ice coring at West Hercules Dome: results from ground-based geophysics and modeling

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

Towards a common terminology in radioglaciology

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

Mapping of the Subglacial Topography of Folgefonna Ice Cap in Western Norway—Consequences for Ice Retreat Patterns and Hydrological Changes

Folgefonna consists of three ice caps which are rapidly retreating in response to warmer temperatures. The melting of Folgefonna has implications for meltwater drainage and hydropower production, as well as the potential for geohazards and impacts to tourism, the communities and infrastructures surrounding the glacier. To support future adaptation strategies, we need to know the subglacial topography of the ice caps to identify water divides and possible areas for geohazards. Therefore, we mapped the subglacial topography at Sørfonna, the largest of the Folgefonna ice caps, using an ice-penetrating radar (2.5 MHz antennas; 1,000 × 500 m grid). The results show a highly irregular subglacial landscape, with deep valleys and high mountain peaks. The maximum ice thickness is 570 m and the mean ice thickness is 190 m. We examined the retreat pattern of Sørfonna using the subglacial topography map in combination with a simple ice flow model and simulated the ice retreat 150 years into the future. We used two climate scenarios (one with a 1.5°C warming and a 3% increase in precipitation, and a second with a 3.5°C warming together with 15% increase in precipitation) and focused on how the glacial retreat will cause hydrological changes in the catchments surrounding the glacier. The main drainage pattern shifts during glacial retreat, with a larger proportion of southward drainage compared to the present day. The ice flow modelling also reveals that the southern part of Sørfonna is more durable during climate change whereas the thinner part of the ice cap, in the north, melts faster. We suggest that increased winter precipitation in a future warmer climate makes the southern part of Sørfonna more resilient than many other glaciers in southern Norway. The subglacial topography map and the retreat pattern also reveal areas that may accumulate water and could potentially generate a future glacial outburst flood. Sediments from distal glacier-fed lakes around Sørfonna have been used to constrain the thresholds identified on the subglacial topography map. Combining sedimentological evidence from distal glacier-fed lakes with the new subglacial topography map confirms that the retreat of specific outlet glaciers, such as Bondhusbreen, Buerbreen, and Møsevassbreen, will have a large impact on meltwater routing, as they are situated behind bedrock thresholds in the upper part of the glacier’s catchment area.