Airborne Radar Sounding Research Articles

Characterizing sub-glacial hydrology using radar simulations

Abstract. The structure and distribution of sub-glacial water directly influences Antarctic ice mass loss by reducing or enhancing basal shear stress and accelerating grounding line retreat. A common technique for detecting sub-glacial water involves analyzing the spatial variation in reflectivity from an airborne radar echo sounding (RES) survey. Basic RES analysis exploits the high dielectric contrast between water and most other substrate materials, where a reflectivity increase ≥ 15 dB is frequently correlated with the presence of sub-glacial water. There are surprisingly few additional tools to further characterize the size, shape, or extent of hydrological systems beneath large ice masses. We adapted an existing radar backscattering simulator to model RES reflections from sub-glacial water structures using the University of Texas Institute for Geophysics (UTIG) Multifrequency Airborne Radar Sounder with Full-phase Assessment (MARFA) instrument. Our series of hypothetical simulation cases modeled water structures from 5 to 50 m wide, surrounded by bed materials of varying roughness. We compared the relative reflectivity from rounded Röthlisberger channels and specular flat canals, showing both types of channels exhibit a positive correlation between size and reflectivity. Large (> 20 m), flat canals can increase reflectivity by more than 20 dB, while equivalent Röthlisberger channels show only modest reflectivity gains of 8–13 dB. Changes in substrate roughness may also alter observed reflectivity by 3–6 dB. All of these results indicate that a sophisticated approach to RES interpretation can be useful in constraining the size and shape of sub-glacial water features. However, a highly nuanced treatment of the geometric context is necessary. Finally, we compared simulated outputs to actual reflectivity from a single RES flight line collected over Thwaites Glacier in 2022. The flight line crosses a previously proposed Röthlisberger channel route, with an obvious bright bed reflection in the radargram. Through multiple simulations comparing various water system geometries, such as canals and sub-glacial lakes, we demonstrated the important role that topography and water geometry can play in observed RES reflectivity. From the scenarios that we tested, we concluded the bright reflector from our RES flight line cannot be a Röthlisberger channel but could be consistent with a series of flat canals or a sub-glacial lake. However, we note our simulations were not exhaustive of all possible sub-glacial water configurations. The approach outlined here has broad applicability for studying the basal environment of large glaciers. We expect to apply this technique when constraining the geometry and extent of many sub-glacial hydrologic structures in the future. Further research may also include comprehensive investigations of the impact of sub-glacial roughness, substrate heterogeneity, and computational efficiencies enabling more complex and complete simulations.

Comprehensive Radar Mapping of Malaspina Glacier (Sít' Tlein), Alaska—The World's Largest Piedmont Glacier—Reveals Potential for Instability

AbstractMalaspina Glacier, located on the coast of southern Alaska, is the world's largest piedmont glacier. A narrow ice‐cored foreland zone undergoing rapid thermokarst erosion separates the glacier from the relatively warm waters of the Gulf of Alaska. Glacier‐wide thinning rates for Malaspina are greater than 1 m/yr, and previous geophysical investigations indicated that bed elevation exceeds 300 m below sea level in some places. These observations together give rise to the question of glacial stability. To address this question, glacier evolution models are dependent upon detailed observations of Malaspina's subglacial topography. Here, we map 2,000 line‐km of the glacier's bed using airborne radar sounding data collected by NASA's Operation IceBridge. When compared to gridded radar measurements, we find that glaciological models overestimate Malaspina's volume by more than 30%. While we report a mean bed elevation 100 m greater than previous models, we find that Malaspina inhabits a broad basin largely grounded below sea level. Several subglacial channels dissect the glacier's bed: the most prominent of these channels extends at least 35 km up‐glacier from the terminus toward the throat of Seward Glacier. Provided continued foreland erosion, an ice‐ocean connection may promote rapid retreat along these overdeepened subglacial channels, with a global sea‐level rise potential of 1.4 mm.

GBaTSv2: a revised synthesis of the likely basal thermal state of the Greenland Ice Sheet

Abstract. The basal thermal state (frozen or thawed) of the Greenland Ice Sheet is under-constrained due to few direct measurements, yet knowledge of this state is becoming increasingly important to interpret modern changes in ice flow. The first synthesis of this state relied on inferences from widespread airborne and satellite observations and numerical models, for which most of the underlying datasets have since been updated. Further, new and independent constraints on the basal thermal state have been developed from analysis of basal and englacial reflections observed by airborne radar sounding. Here we synthesize constraints on the Greenland Ice Sheet's basal thermal state from boreholes, thermomechanical ice-flow models that participated in the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6; Coupled Model Intercomparison Project Phase 6), IceBridge BedMachine Greenland v4 bed topography, Making Earth Science Data Records for Use in Research Environments (MEaSUREs) Multi-Year Greenland Ice Sheet Velocity Mosaic v1 and multiple inferences of a thawed bed from airborne radar sounding. Most constraints can only identify where the bed is likely thawed rather than where it is frozen. This revised synthesis of the Greenland likely Basal Thermal State version 2 (GBaTSv2) indicates that 33 % of the ice sheet's bed is likely thawed, 40 % is likely frozen and the remainder (28 %) is too uncertain to specify. The spatial pattern of GBaTSv2 is broadly similar to the previous synthesis, including a scalloped frozen core and thawed outlet-glacier systems. Although the likely basal thermal state of nearly half (46 %) of the ice sheet changed designation, the assigned state changed from likely frozen to likely thawed (or vice versa) for less than 6 % of the ice sheet. This revised synthesis suggests that more of northern Greenland is likely thawed at its bed and conversely that more of southern Greenland is likely frozen, both of which influence interpretation of the ice sheet's present subglacial hydrology and models of its future evolution. The GBaTSv2 dataset, including both code that performed the analysis and the resulting datasets, is freely available at https://doi.org/10.5281/zenodo.6759384 (MacGregor, 2022).

Noise Removal and Feature Extraction in Airborne Radar Sounding Data of Ice Sheets

The airborne ice-penetrating radar (IPR) is an effective method used for ice sheet exploration and is widely applied for detecting the internal structures of ice sheets and for understanding the mechanism of ice flow and the characteristics of the bottom of ice sheets. However, because of the ambient influence and the limitations of the instruments, IPR data are frequently overlaid with noise and interference, which further impedes the extraction of layer features and the interpretation of the physical characteristics of the ice sheet. In this paper, we first applied conventional filtering methods to remove the feature noise and interference in IPR data. Furthermore, machine learning methods were introduced in IPR data processing for noise removal and feature extraction. Inspired by a comparison of the filtering methods and machine learning methods, we propose a fusion method combining both filtering methods and machine-learning-based methods to optimize the feature extraction in IPR data. Field data tests indicated that, under different conditions of IPR data, the application of different methods and strategies can improve the layer feature extraction.

Polarimetric Airborne Radar Sounding as an Approach to Characterizing Subglacial Röthlisberger Channels

Conventional airborne radar sounding techniques are well suited to the detection and characterization of flat-lying, specular subglacial water bodies. However, topographically positive, diffusively scattering Rthlisberger (R-)channels are more difficult to image, while also exerting substantial control on basal friction and ice dynamics. As subglacial R-channels share geometrical similarities with targets of interest in polarimetric ground-penetrating radar studies (i.e., cylindrical pipes), in this study we investigate whether similar concepts can be adapted to detect and characterize R-channels. While closed form analytical solutions exist for the scattering widths of perfect electrically conducting (PEC) and dielectric circular cylinders, the insight they provide for the polarimetric response of half-cylinder R-channels is limited. As such, a series of modeling experiments have been performed to characterize the scattering widths of half cylinders. Our results demonstrate that scattering from subglacial R-channels depends on numerous factors including the polarization and frequency of the incident radar wave, the size of the R-channel, and the relative orientation of the R-channel to a pair of orthogonally oriented, linearly polarized radar antennas. The results imply that patterns in the like-polarized echo powers across the signal bandwidth may be useful in inferring R-channel existence and possibly estimate R-channel size. However, as differences in scattering width are small and decrease with increasing misalignment between the acquisition system antennas and R-channel strike, achieving a high degree of radiometric resolution and stability should be a driving factor in polarimetric radar system design.

Multisystem Synthesis of Radar Sounding Observations of the Amundsen Sea Sector From the 2004-2005 Field Season.

The Amundsen Sea Embayment of the West Antarctic Ice Sheet contains Thwaites and Pine Island Glaciers, two of the most rapidly changing glaciers in Antarctica. To date, Pine Island and Thwaites Glaciers have only been observed by independent airborne radar sounding surveys, but a combined cross‐basin analysis that investigates the basal conditions across the Pine Island‐Thwaites Glaciers boundary has not been performed. Here, we combine two radar surveys and correct for their differences in system parameters to produce unified englacial attenuation and basal relative reflectivity maps spanning both Pine Island and Thwaites Glaciers. Relative reflectivities range from −24.8 to +37.4 dB with the highest values beneath fast‐flowing ice at the ice sheet margin. By comparing our reflectivity results with previously derived radar specularity and trailing bed echoes at Thwaites Glacier, we find a highly diverse subglacial landscape and hydrologic conditions that evolve along‐flow. Together, these findings highlight the potential for joint airborne radar analysis with ground‐based seismic and geomorphological observations to understand variations in the bed properties and cross‐catchment interactions of ice streams and outlet glaciers.

Brief communication: An empirical relation between center frequency and measured thickness for radar sounding of temperate glaciers

Abstract. Radar sounding of the thickness of temperate glaciers is challenged by substantial volume scattering, surface scattering and high attenuation rates. Lower-frequency radar sounders are often deployed to mitigate these effects, but the lack of a global synthesis of their success limits progress in system and survey design. Here we extend a recent global compilation of glacier thickness measurements (GlaThiDa) with the center frequency for radar-sounding surveys. From a maximum reported thickness of ∼ 1500 m near 1 MHz, the maximum thickness sounded decreases by ∼ 500 m per frequency decade. Between 25–100 MHz, newer airborne radar sounders generally outperform older, ground-based ones. Based on globally modeled glacier thicknesses, we conclude that a multi-element, ≤30 MHz airborne radar sounder could survey most temperate glaciers more efficiently.

The regional-scale surface mass balance of Pine Island Glacier, West Antarctica, over the period 2005–2014, derived from airborne radar soundings and neutron probe measurements

Abstract. We derive recent surface mass balance (SMB) estimates from airborne radar observations along the iSTAR traverse (2013, 2014) at Pine Island Glacier (PIG), West Antarctica. Ground-based neutron probe measurements provide information of snow and firn density with depth at 22 locations and were used to date internal annual reflection layers. The 2005 layer was traced for a total distance of 2367 km to determine annual mean SMB for the period 2005–2014. Using complementary SMB estimates from two regional climate models, RACMO2.3p2 and MAR, and a geostatistical kriging scheme, we determine a regional-scale SMB distribution with similar main characteristics to that determined for the period 1985–2009 in previous studies. Local departures exist for the northern PIG slopes, where the orographic precipitation shadow effect appears to be more pronounced in our observations, and the southward interior, where the SMB gradient is more pronounced in previous studies. We derive total mass inputs of 79.9±19.2 and 82.1±19.2 Gt yr−1 to the PIG basin based on complementary ASIRAS–RACMO and ASIRAS–MAR SMB estimates, respectively. These are not significantly different to the value of 78.3±6.8 Gt yr−1 for the period 1985–2009. Thus, there is no evidence of a secular trend at decadal scales in total mass input to the PIG basin. We note, however, that our estimated uncertainty is more than twice the uncertainty for the 1985–2009 estimate on total mass input. Our error analysis indicates that uncertainty estimates on total mass input are highly sensitive to the selected krige methodology and assumptions made on the interpolation error, which we identify as the main cause for the increased uncertainty range compared to the 1985–2009 estimates.

Five decades of radioglaciology

AbstractRadar sounding is a powerful geophysical approach for characterizing the subsurface conditions of terrestrial and planetary ice masses at local to global scales. As a result, a wide array of orbital, airborne, ground-based, and in situ instruments, platforms and data analysis approaches for radioglaciology have been developed, applied or proposed. Terrestrially, airborne radar sounding has been used in glaciology to observe ice thickness, basal topography and englacial layers for five decades. More recently, radar sounding data have also been exploited to estimate the extent and configuration of subglacial water, the geometry of subglacial bedforms and the subglacial and englacial thermal states of ice sheets. Planetary radar sounders have observed, or are planned to observe, the subsurfaces and near-surfaces of Mars, Earth's Moon, comets and the icy moons of Jupiter. In this review paper, and the thematic issue of theAnnals of Glaciologyon ‘Five decades of radioglaciology’ to which it belongs, we present recent advances in the fields of radar systems, missions, signal processing, data analysis, modeling and scientific interpretation. Our review presents progress in these fields since the last radio-glaciologicalAnnals of Glaciologyissue of 2014, the context of their history and future prospects.

Multidecadal observations of the Antarctic ice sheet from restored analog radar records

Airborne radar sounding can measure conditions within and beneath polar ice sheets. In Antarctica, most digital radar-sounding data have been collected in the last 2 decades, limiting our ability to understand processes that govern longer-term ice-sheet behavior. Here, we demonstrate how analog radar data collected over 40 y ago in Antarctica can be combined with modern records to quantify multidecadal changes. Specifically, we digitize over 400,000 line kilometers of exploratory Antarctic radar data originally recorded on 35-mm optical film between 1971 and 1979. We leverage the increased geometric and radiometric resolution of our digitization process to show how these data can be used to identify and investigate hydrologic, geologic, and topographic features beneath and within the ice sheet. To highlight their scientific potential, we compare the digitized data with contemporary radar measurements to reveal that the remnant eastern ice shelf of Thwaites Glacier in West Antarctica had thinned between 10 and 33% between 1978 and 2009. We also release the collection of scanned radargrams in their entirety in a persistent public archive along with updated geolocation data for a subset of the data that reduces the mean positioning error from 5 to 2.5 km. Together, these data represent a unique and renewed extensive, multidecadal historical baseline, critical for observing and modeling ice-sheet change on societally relevant timescales.

Rapid submarine ice melting in the grounding zones of ice shelves in West Antarctica.

Enhanced submarine ice-shelf melting strongly controls ice loss in the Amundsen Sea embayment (ASE) of West Antarctica, but its magnitude is not well known in the critical grounding zones of the ASE's major glaciers. Here we directly quantify bottom ice losses along tens of kilometres with airborne radar sounding of the Dotson and Crosson ice shelves, which buttress the rapidly changing Smith, Pope and Kohler glaciers. Melting in the grounding zones is found to be much higher than steady-state levels, removing 300–490 m of solid ice between 2002 and 2009 beneath the retreating Smith Glacier. The vigorous, unbalanced melting supports the hypothesis that a significant increase in ocean heat influx into ASE sub-ice-shelf cavities took place in the mid-2000s. The synchronous but diverse evolutions of these glaciers illustrate how combinations of oceanography and topography modulate rapid submarine melting to hasten mass loss and glacier retreat from West Antarctica.

Estimating Subglacial Water Geometry Using Radar Bed Echo Specularity: Application to Thwaites Glacier, West Antarctica

Airborne radar sounding is an established tool for observing the bed conditions and subglacial hydrology of ice sheets and glaciers. The specularity content of radar bed echoes has also been used to detect the hydrologic transition of a subglacial water system from a network of distributed canals to a network of concentrated channels beneath the Thwaites Glacier. However, the physical dimensions of the distributed water bodies in these networks have not been constrained by observations. In this letter, we use a variety of simple radar scattering, attenuation, and cross-sectional models to provide a first estimate of the subglacial water body geometries capable of producing the observed anisotropic specularity of the Thwaites Glacier catchment. This approach leads to estimates of ice/water interface root mean square roughnesses less than about 15 cm, thicknesses of more than about 5 cm, lengths of more than about 15 m, and widths between about 0.5 and 5 m.

Two-Dimensional Imaging of Ice Sheets of Airborne Radar Sounder via a Combined Modified Range Migration Algorithm Based on ISFT and Beamforming Using Curvelets

We develop and demonstrate a new method to address the surface clutter problem in imaging of ice sheets by airborne radar sounder. This new method utilizes the modified range migration algorithm based on inverse scaled Fourier transform (ISFT) to reduce along-track clutter and the beamforming using curvelets to reduce cross-track clutter. First, theory analysis has been carried on to the proposed method. Then, we give the particular realizing steps to implement this method. At last, we apply this method to the simulation point targets and MCoRDS data to prove its validity of imaging of ice sheets. Furthermore, we show the advantages of the proposed method through comparison with two methods in four major aspects—the power of along-track clutter reduction, the power of cross-track clutter reduction, the equivalent number of looks (ENLs), and the calculating time.

A Model-Based Technique for the Automatic Detection of Earth Continental Ice Subsurface Targets in Radar Sounder Data

The continuous melting of the ice at the Earth continental polar caps highlights the importance of an exhaustive study of the properties of the ice subsurface targets in order to provide a reliable analysis of their past and future evolution. Such study can be efficiently performed by automatically analyzing radargrams of the ice cross section acquired by radar sounder (RS) instruments. In this letter, we propose an automatic technique for a large-scale detection of the ice subsurface targets and the estimation of their properties (e.g., layered area thickness and bedrock scattering area) from radargrams acquired by RS operated at the Earth continental polar caps. This is done by using the parameters of the RS acquisition system combined with the output of an automatic image segmentation algorithm. The segmentation operation is applied to the radargrams after a preliminary processing phase aimed to emphasize the relevant subsurface targets. The segmentation criterion considers the radar signal backscattering properties and a model of the spatial distribution of the investigated targets that takes into account the effects of the wave propagation though the subsurface. Experimental results obtained on real radargrams acquired by an airborne RS in Antarctica confirm the effectiveness of the proposed technique.

Surface slope control on firn density at Thwaites Glacier, West Antarctica: Results from airborne radar sounding

AbstractWe derive the surface density variations over Thwaites Glacier from a statistical analysis of airborne High‐Capability Radar Sounder surface returns acquired in the 2004–2005 summer. We produce a 5 km gridded map with an estimated ±12.5 kg m−3 accuracy. The background pattern of densities decreases inland from ~ 480 kg m−3 to ~ 350 kg m−3. A remarkable ~ 30 km wide vein‐shaped anomaly of up to 570 kg m−3 is located ~ 100 km from the coastline. Density anomalies correspond with surface slope breaks but are not necessarily coincident with smaller slopes. They could result from complex wind‐driven snow redistribution and/or refreezing of former snowmelt. This inversion technique can significantly improve surface mass balance calculations to understanding of glacier dynamics at regional scales and is valuable to verify and improve Antarctic climate models. It is also a promising approach for future surface analysis of icy moons by planetary radars.

Planetary landing-zone reconnaissance using ice-penetrating radar data: Concept validation in Antarctica

The potential for a nadir-looking radar sounder to retrieve significant surface roughness/permittivity information valuable for planetary landing site selection is demonstrated using data from an airborne survey of the Thwaites Glacier Catchment, West Antarctica using the High Capability Airborne Radar Sounder (HiCARS). The statistical method introduced by Grima et al. (2012. Icarus 220, 84–99. http://dx.doi.org/10.1007/s11214-012-9916-y) for surface characterization is applied systematically along the survey flights. The coherent and incoherent components of the surface signal, along with an internally generated confidence factor, are extracted and mapped in order to show how a radar sounder can be used as both a reflectometer and a scatterometer to identify regions of low surface roughness compatible with a planetary lander. These signal components are used with a backscattering model to produce a landing risk assessment map by considering the following surface properties: Root mean square (RMS) heights, RMS slopes, roughness homogeneity/stationarity over the landing ellipse, and soil porosity. Comparing these radar-derived surface properties with simultaneously acquired nadir-looking imagery and laser-altimetry validates this method. The ability to assess all of these parameters with an ice penetrating radar expands the demonstrated capability of a principle instrument in icy planet satellite science to include statistical reconnaissance of the surface roughness to identify suitable sites for a follow-on lander mission.

Low-frequency radar sounding of ice in East Antarctica and southern Greenland

AbstractWe discuss a decameter-wavelength airborne radar sounder, the Warm Ice Sounding Explorer (WISE), that provides ice thickness in areas where radar signal penetration at higher frequencies is expected to be limited. Here we report results for three campaigns conducted in Greenland (2008, 2009, 2010) and two in Antarctica (2009, 2010). Comparisons with higher-frequency radar data indicate an accuracy of ±55 m for ice-thickness measurements in Greenland and ±25 m in Antarctica. We also estimate ice thickness of the Qassimiut lobe in southwest Greenland, where few ice-thickness measurements have been made, demonstrating that WISE penetrates in strongly scattering environments.

A mass conservation approach for mapping glacier ice thickness

The traditional method for interpolating ice thickness data from airborne radar sounding surveys onto regular grids is to employ geostatistical techniques such as kriging. While this approach provides continuous maps of ice thickness, it generates products that are not consistent with ice flow dynamics and are impractical for high resolution ice flow simulations. Here, we present a novel approach that combines sparse ice thickness data collected by airborne radar sounding profilers with high resolution swath mapping of ice velocity derived from satellite synthetic-aperture interferometry to obtain a high resolution map of ice thickness that conserves mass and minimizes the departure from observations. We apply this approach to the case of Nioghalvfjerdsfjorden (79North) Glacier, a major outlet glacier in northeast Greenland that has been relatively stable in recent decades. The results show that our mass conserving method removes the anomalies in mass flux divergence, yields interpolated data that are within about 5% of the original data, and produces thickness maps that are directly usable in high spatial-resolution, high-order ice flow models. We discuss the application of this method to the broad and detailed radar surveys of ice sheet and glacier thickness. Copyright 2011 by the American Geophysical Union.

Ice flux divergence anomalies on 79north Glacier, Greenland

The ice flux divergence of a glacier is an important quantity to examine because it determines the rate of temporal change of its thickness. Here, we combine high‐resolution ice surface velocity observations of Nioghalvfjerdsfjorden (79north) Glacier, a major outlet glacier in north Greenland, with a dense grid of ice thickness data collected with an airborne radar sounder in 1998, to examine its ice flux divergence. We detect large variations, up to 100 m/yr, in flux divergence on grounded ice that are incompatible with what we know of the glacier surface mass balance, basal mass balance and thinning rate. We examine the hypothesis that these anomalies are due to the three‐dimensional flow of ice around and atop bumps and hollows in basal topography by comparing the flux divergence of three‐dimensional numerical models with its surface equivalent. We find that three‐dimensional effects have only a small contribution to the observed anomalies. On the other hand, if we degrade the spatial resolution of the data to 10 km the anomalies disappear. Further analysis shows that the source of the anomalies is not the ice velocity data but the interpolation of multiple tracks of ice thickness data onto a regular grid using a scheme (here block kriging) that does not conserve mass or ice flux. This problem is not unique to 79north Glacier but is common to all conventional ice thickness surveys of glaciers and ice sheets; and fundamentally limits the application of ice thickness grids to high‐resolution numerical modeling of glacier flow.

Marine ice in Larsen Ice Shelf

It is argued that Larsen Ice Shelf contains marine ice formed by oceanic freezing and other mechanisms. Missing basal returns in airborne radar soundings and observations of a smooth and healed surface coincide downstream of regions where an ocean model predicts freezing. Visible imagery suggests that marine ice currently stabilizes Larsen C Ice Shelf and implicates failure of marine flow bands in the 2002 Larsen B Ice Shelf collapse. Ocean modeling indicates that any regime change towards the incursion of warmer Modified Weddell Deep Water into the Larsen C cavity could curtail basal freezing and its stabilizing influence.