Ice-penetrating Radar Data Research Articles

A newly digitized ice-penetrating radar data set acquired over the Greenland ice sheet in 1971–1979

Abstract. We present an ice-penetrating radar data set acquired over the Greenland ice sheet by aircraft during the years 1971, 1972, 1974, 1978, and 1979. The data set comprises over 177 000 km of flight lines and contains a wealth of information on the state of the Greenland ice sheet, including information on ice thickness and englacial properties. During data collection in the 1970s, the data were recorded on optical film rolls, and in this paper, we document the digitization of these film rolls and their associated geographical information. Our data digitization enables interaction with and analysis of the data and facilitates comparison with modern-day radar observations. The complete data set in full resolution is available in the Stanford Digital Repository (https://doi.org/10.25740/wm135gp2721; Karlsson et al., 2023), with the associated technical reports. Part of the data set is available as low-resolution JPG files at the Technical University of Denmark's data repository (https://doi.org/10.6084/m9.figshare.c.7235299.v1; Karlsson et al., 2024), with associated technical reports and digitized geographical information . The Stanford Digital Repository serves as long-term storage, providing archival historic preservation in perpetuity, and is not intended as a primary data access point. The Technical University of Denmark (DTU) data repository serves as a primary entry point for data access, with files organized according to acquisition year and flight line in a simple folder structure. Here, we release the full data sets to enable the larger community to access and interact with the data (Karlsson et al., 2023).

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.

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.

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.

Investigating the spatial representativeness of East Antarctic ice cores: a comparison of ice core and radar-derived surface mass balance over coastal ice rises and Dome Fuji

Abstract. Surface mass balance (SMB) of the Antarctic Ice Sheet must be better understood to document the current Antarctic contribution to sea-level rise. In situ point data using snow stakes and ice cores are often used to evaluate the state of the ice sheet's mass balance, as well as to assess SMB derived from regional climate models, which are then used to produce future climate projections. However, spatial representativeness of individual point data remains largely unknown, particularly in the coastal regions of Antarctica with highly variable terrain. Here, we compare ice core data collected at the summit of eight ice rises along the coast of Dronning Maud Land, as well as at the Dome Fuji site, and shallow ice-penetrating radar data over these regions. Shallow radar data have the advantage of being spatially extensive, with a temporal resolution that varies between a yearly and multi-year resolution, from which we can derive a SMB record over the entire radar survey. This comparison therefore allows us to evaluate the spatial variability of SMB and the spatial representativeness of ice-core-derived SMB. We found that ice core mean SMB is very local, and the difference with radar-derived SMB increases in a logarithmic fashion as the surface covered by the radar data increases, with a plateau ∼ 1–2 km away from the ice crest for most ice rises, where there are strong wind–topography interactions, and ∼ 10 km where the ice shelves begin. The relative uncertainty in measuring SMB also increases rapidly as we move away from the ice core sites. Five of our ice rise sites show a strong spatial representativeness in terms of temporal variability, while the other three sites show that it is limited to a surface area between 20–120 km2. The Dome Fuji site, on the other hand, shows a small difference between pointwise and area mean SMB, as well as a strong spatial representativeness in terms of temporal variability. We found no simple parameterization that could represent the spatial variability observed at all the sites. However, these data clearly indicate that local spatial SMB variability must be considered when assessing mass balance, as well as comparing modeled SMB values to point field data, and therefore must be included in the estimate of the uncertainty of the observations.

Linear analysis of ice-shelf topography response to basal melting and freezing

Floating ice shelves in Antarctica and Greenland limit land-ice contributions to sea-level rise by resisting the flow of grounded ice. Melting at the surface and base of ice shelves can lead to destabilization by promoting thinning and fracturing. Basal melting often results in channelized features that manifest as surface topography due to buoyancy. The assumption of hydrostatic flotation commonly underlies estimates of basal melting rates. However, numerical simulations and ice-penetrating radar data have shown that narrow topographic features do not necessarily satisfy the local flotation condition. Here, we introduce a linearized model for ice-shelf topographic response to basal melting perturbations to quantify deviations from hydrostatic flotation and the stability of topography. While hydrostatic flotation is the dominant behaviour at wavelengths greater than the ice thickness, ice elevation can deviate from the perfect flotation condition at smaller wavelengths. The linearized analysis shows that channelized features can be stable when the time scale of extensional thinning is small relative to the time scale of viscous flow towards the channel. When extension is non-negligible, channels can break through the ice column. We validate the linearized analysis by comparing numerical solutions to a nonlinear ice-flow model with steady-state solutions obtained via a Green’s function.

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.

Quantifying Basal Roughness and Internal Layer Continuity Index of Ice Sheets by an Integrated Means with Radar Data and Deep Learning

Understanding englacial and subglacial structures is a fundamental method of inferring ice sheets’ historical evolution and surface mass balance. The internal layer continuity index and the basal roughness are key parameters and indicators for the speculation of the relationship between the ice sheet’s internal structure or bottom and ice flow. Several methods have been proposed in the past two decades to quantitatively calculate the continuity index of ice layer geometry and the roughness of the ice–bedrock interface based on radar echo signals. These methods are mainly based on the average of the absolute value of the vertical gradient of the echo signal amplitude and the standard deviation of the horizontal fluctuation of the bedrock interface. However, these methods are limited by the amount and quality of unprocessed radar datasets and have not been widely used, which also hinders further research, such as the analysis of the englacial reflectivity, the subglacial conditions, and the history of the ice sheets. In this paper, based on geophysical processing methods for radar image denoising and deep learning for ice layer and bedrock interface extraction, we propose a new method for calculating the layer continuity index and basal roughness. Using this method, we demonstrate the ice-penetrating radar data processing and compare the imaging and calculation of the radar profiles from Dome A to Zhongshan Station, East Antarctica. We removed the noise from the processed radar data, extracted ice layer continuity features, and used other techniques to verify the calculation. The potential application of this method in the future is illustrated by several examples. We believe that this method can become an effective approach for future Antarctic geophysical and glaciological research and for obtaining more information about the history and dynamics of ice sheets from their radar-extracted internal structure.

Paths forward in radioglaciology

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

From ice core to ground-penetrating radar: representativeness of SMB at three ice rises along the Princess Ragnhild Coast, East Antarctica

AbstractThe future contributions of the Antarctic Ice Sheet to sea level rise will depend on the evolution of its surface mass balance (SMB), which could amplify/dampen mass losses increasingly observed at the ice sheet's edge. In situ constraints of SMB over annual-to-decadal timescales consist mostly of firn/ice cores that have a surface footprint $\sim$cm$^{2}$. SMB constraints also come from climate models, which have a higher temporal resolution but a larger surface footprint of several km$^{2}$. We use ice-penetrating radar data to obtain an intermediate spatial and temporal resolution SMB record over three ice rises along the Princess Ragnhild Coast. The co-located ice cores allow us to obtain absolute radar-derived SMB rates at a multi-annual-to-decadal temporal resolution. By comparing the ice core SMB measurements and the radar-derived SMB records, we determine that pointwise measurements of SMB are representative of a small surface area, $\sim 200-500$ m radius extending from the ice core drill site for the ice rises studied here, and that the pointwise measurements are systematically 7–15 cm w.e. a$^{-1}$ lower than the mean SMB value calculated for the whole ice rises. However, ice core records are representative of an entire ice rise's temporal variability at the temporal resolution examined.

Ice Thickness From Deep Learning and Conditional Random Fields: Application to Ice-Penetrating Radar Data With Radiometric Validation

Identifying the location of the ice-bedrock interface of glaciers and ice sheets is crucial for a wide range of geophysical applications, such as searching for liquid water in basal regions and computing ice thickness to quantify ice sheet and glacier mass balance. Simple, record-by-record, approaches to detecting the bottom of the ice echo may be affected by spurious off-nadir noise that requires significant manual interaction to correct. In this paper, we propose a deep learning model based on convolutional neural networks (CNNs) and continuous conditional random fields (CCRFs) to automate ice bed identification and better capture fine-grained basal detail. We deploy this approach on HiCARS radargrams, the first time deep learning methods have been applied to this dataset. Intuitively, our CNN captures the global geometry of the ice bed, while the CCRF adjusts initial CNN outputs to better incorporate fine-scale spatial information into the final prediction. We also develop a coherent geophysical framework using three echo characters (along-track continuity, relative delay, and signal coherency) to compare our model’s outputs with those of a manually targeted approach. Our analysis suggests that our CNN + CCRF model is as suitable as the manual approach for radiometric applications, and it outperforms the manual technique in identifying the first continuous return, which is most often the near-nadir reflection. Thus, our approach is more universal than the current manual labeling methodology in the range of geophysical applications where it can be used, and it provides better confidence regarding the source of the basal return detected.

Radiometric analysis of digitized Z-scope records in archival radar sounding film

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

A detailed radiostratigraphic data set for the central East Antarctic Plateau spanning from the Holocene to the mid-Pleistocene

Abstract. We present an ice-penetrating radar data set which consists of 26 internal reflecting horizons (IRHs) that cover the entire Dome C area of the East Antarctic plateau, the most extensive to date in the region. This data set uses radar surveys collected over the space of 10 years, starting with an airborne international collaboration in 2008 to explore the region, up to the detailed ground-based surveys in support of the Beyond EPICA – Oldest Ice (BE-OI) European Consortium. Through direct correlation with the EPICA-DC ice core, we date 19 IRHs that span the past four glacial cycles, from 10 ka, beginning of the Holocene, to over 350 ka, ranging from 10 % to 83 % of the ice thickness at the EPICA-DC ice core site. We indirectly date and provide stratigraphic information for seven older IRHs using a 1D ice flow inverse model, going back to an estimated 700 ka. Depth and age uncertainties are quantified for all IRHs and provided as part of the data set. The IRH data set presented in this study is available at the US Antarctic Program Data Center (USAP-DC) (https://doi.org/10.15784/601411, Cavitte et al., 2020) and represents a contribution to the SCAR AntArchitecture action group (AntArchitecture, 2017).

Radio-echo sounding and waveform modeling reveal abundant marine ice in former rifts and basal crevasses within Crary Ice Rise, Antarctica

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

An imbalancing act: the delayed dynamic response of the Kaskawulsh Glacier to sustained mass loss

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

Stochastic modeling of subglacial topography exposes uncertainty in water routing at Jakobshavn Glacier

AbstractSubglacial topography is an important feature in numerous ice-sheet analyses and can drive the routing of water at the bed. Bed topography is primarily measured with ice-penetrating radar. Significant gaps, however, remain in data coverage that require interpolation. Topographic interpolations are typically made with kriging, as well as with mass conservation, where ice flow dynamics are used to constrain bed geometry. However, these techniques generate bed topography that is unrealistically smooth at small scales, which biases subglacial water flowpath models and makes it difficult to rigorously quantify uncertainty in subglacial drainage patterns. To address this challenge, we adapt a geostatistical simulation method with probabilistic modeling to stochastically simulate bed topography such that the interpolated topography retains the spatial statistics of the ice-penetrating radar data. We use this method to simulate subglacial topography using mass conservation topography as a secondary constraint. We apply a water routing model to each of these realizations. Our results show that many of the flowpaths significantly change with each topographic realization, demonstrating that geostatistical simulation can be useful for assessing confidence in subglacial flowpaths.

A framework for attenuation method selection evaluated with ice-penetrating radar data at South Pole Lake

AbstractAll radar power interpretations require a correction for attenuative losses. Moreover, radar attenuation is a proxy for ice-column properties, such as temperature and chemistry. Prior studies use either paired thermodynamic and conductivity models or the radar data themselves to calculate attenuation, but there is no standard method to do so; and, before now, there has been no robust methodological comparison. Here, we develop a framework meant to guide the implementation of empirical attenuation methods based on survey design and regional glaciological conditions. We divide the methods into the three main groups: (1) those that infer attenuation from a single reflector across many traces; (2) those that infer attenuation from multiple reflectors within one trace; and (3) those that infer attenuation by contrasting the measured power from primary and secondary reflections. To assess our framework, we introduce a new ground-based radar survey from South Pole Lake, comparing selected empirical methods to the expected attenuation from a temperature- and chemistry-dependent Arrhenius model. Based on the small surveyed area, lack of a sufficient calibration surface and low reflector relief, the attenuation methods that use multiple reflectors are most suitable at South Pole Lake.

Interferometric discrimination of cross-track bed clutter in ice-penetrating radar sounding data

AbstractThe interpretations of relevant interfaces (i.e. the surface and bed) in radar sounding datasets over glaciers and ice sheets are primary boundary conditions in a variety of climate studies and particularly subglacial water routing models. It is therefore necessary to ensure these interpretations are consistent and not affected by cross-track clutter. For the surface interface, interferometry and a family of methods relying on digital elevation models have been used to successfully discriminate cross-track surface clutter. Here we present how interferometry can be applied to the problem of basal clutter from cross-track bed topography. Our approach is based on a comparison of the differential phases of ambiguous reflectors that may represent bed clutter and the differential phase of a reflector in an adjacent area that appears unaffected by basal clutter. The reflector yielding the smallest interferometric phase difference relative to the unambiguous bed reflector is considered to represent its consistent continuation. We successfully demonstrate our approach using 60 MHz center frequency MARFA data collected over Devon Ice Cap in the Canadian Arctic. Finally, we investigate the effects of clutter-affected and interferometry-corrected bed interpretations on ice layer thickness estimates, basal hydraulic head gradients and the potential extent of inferred subglacial water bodies.

Multidecadal Basal Melt Rates and Structure of the Ross Ice Shelf, Antarctica, Using Airborne Ice Penetrating Radar

AbstractBasal melting of ice shelves is a major source of mass loss from the Antarctic Ice Sheet. In situ measurements of ice shelf basal melt rates are sparse, while the more extensive estimates from satellite altimetry require precise information about firn density and characteristics of near‐surface layers. We describe a novel method for estimating multidecadal basal melt rates using airborne ice penetrating radar data acquired during a 3‐year survey of the Ross Ice Shelf. These data revealed an ice column with distinct upper and lower units whose thicknesses change as ice flows from the grounding line toward the ice front. We interpret the lower unit as continental meteoric ice that has flowed across the grounding line and the upper unit as ice formed from snowfall onto the relatively flat ice shelf. We used the ice thickness difference and strain‐induced thickness change of the lower unit between the survey lines, combined with ice velocities, to derive basal melt rates averaged over one to six decades. Our results are similar to satellite laser altimetry estimates for the period 2003–2009, suggesting that the Ross Ice Shelf melt rates have been fairly stable for several decades. We identify five sites of elevated basal melt rates, in the range 0.5–2 m a−1, near the ice shelf front. These hot spots indicate pathways into the sub‐ice‐shelf ocean cavity for warm seawater, likely a combination of summer‐warmed Antarctic Surface Water and modified Circumpolar Deep Water, and are potential areas of ice shelf weakening if the ocean warms.

Seawater softening of suture zones inhibits fracture propagation in Antarctic ice shelves

Suture zones are abundant on Antarctic ice shelves and widely observed to impede fracture propagation, greatly enhancing ice-shelf stability. Using seismic and radar observations on the Larsen C Ice Shelf of the Antarctic Peninsula, we confirm that such zones are highly heterogeneous, consisting of multiple meteoric and marine ice bodies of diverse provenance fused together. Here we demonstrate that fracture detainment is predominantly controlled by enhanced seawater content in suture zones, rather than by enhanced temperature as previously thought. We show that interstitial seawater can reduce fracture-driving stress by orders of magnitude, promoting both viscous relaxation and the development of micro cracks, the incidence of which scales inversely with stress intensity. We show how simple analysis of viscous buckles in ice-penetrating radar data can quantify the seawater content of suture zones and their modification of the ice-shelf’s stress regime. By limiting fracture, enhancing stability and restraining continental ice discharge into the ocean, suture zones act as vital regulators of Antarctic mass balance.