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

Accurate knowledge of the ice thickness distribution and glacier bed topography is essential for predicting dynamic glacier changes and the future developments of downstream hydrology, which are impacting the energy sector, tourism industry and natural hazard management. Using AIR-ETH, a new helicopter-borne ground-penetrating radar (GPR) platform, we measured the ice thickness of all large and most medium-sized glaciers in the Swiss Alps during the years 2016–20. Most of these had either never or only partially been surveyed before. With this new dataset, 251 glaciers – making up 81% of the glacierized area – are now covered by GPR surveys. For obtaining a comprehensive estimate of the overall glacier ice volume, ice thickness distribution and glacier bed topography, we combined this large amount of data with two independent modeling algorithms. This resulted in new maps of the glacier bed topography with unprecedented accuracy. The total glacier volume in the Swiss Alps was determined to be 58.7 ± 2.5 km3 in the year 2016. By projecting these results based on mass-balance data, we estimated a total ice volume of 52.9 ± 2.7 km3 for the year 2020. Data and modeling results are accessible in the form of the SwissGlacierThickness-R2020 data package.

Basal motion of ice sheets depends in part on the roughness and material properties of the subglacial bed and the occurrence of water. To date, basal motion represents one of the largest uncertainties in ice-flow models. It is that component of the total flow velocity that can change most rapidly and can, therefore, facilitate rapid variations in dynamic behaviour. In this study, we investigate the subglacial properties of the East Antarctic Ice Sheet by statistically analysing the roughness of the bed topography, inferred from radio-echo sounding measurements. We analyse two sets of roughness parameters, one derived in the spatial and the other in the spectral domain, with two roughness parameters each. This enables us to compare the suitability of the four roughness parameters to classify the subglacial landscapes below the ice sheet. We further investigate the relationship of the roughness parameters with observed surface flow velocity and modelled basal temperatures of the ice sheet. We find that one of the roughness parameters, the Hurst exponent derived in the spatial domain, coincides with the thermal condition at the base of the ice sheet for slow flow velocities and varies with flow velocity.

We present sea-ice surface roughness estimates, i.e. the standard deviation of relative surface elevation, in the Arctic regions of Fram Strait and the Nansen Basin north of Svalbard acquired by an airborne laser scanner and a single-beam laser altimeter in 2010. We compare the scanner to the altimeter and compare the differences between the two survey regions. We estimate and correct sensor roll from the scanner data using the hyperbolic response of the scanner over a flat surface. Measurement surveys had to be longer than 5 km north of Svalbard and longer than 15 km in Fram Strait before the statistical distribution in surface roughness from the scanner and altimeter became similar. The shape of the surface roughness probability distributions agrees with those of airborne electromagnetic induction measurements of ice thickness. The ice in Fram Strait had a greater mean surface roughness, 0.16 m vs 0.09 m, and a wider distribution in roughness values than the ice in the Nansen Basin. An increase in surface roughness with increasing ice thickness was observed over fast ice found in Fram Strait near the coast of Greenland but not for the drift ice.

Data from the Seasonal Ice Zone Observing Network (SIZONet) acquired near Barrow, Alaska, during the 2009/10 ice season allow novel comparisons between measurements of ice thickness and velocity. An airborne electromagnetic survey that passed over a moored Ice Profiling Sonar (IPS) provided coincident independent measurements of total ice and snow thickness and ice draft at a scale of 10 km. Once differences in sampling footprint size are accounted for, we reconcile the respective probability distributions and estimate the thickness of level sea ice at 1.48 ± 0.1 m, with a snow depth of 0.12 ± 0.07 m. We also complete what we believe is the first independent validation of radar-derived ice velocities by comparing measurements from a coastal radar with those from an under-ice acoustic Doppler current profiler (ADCP). After applying a median filter to reduce high-frequency scatter in the radar-derived data, we find good agreement with the ADCP bottom-tracked ice velocities. With increasing regulatory and operational needs for sea-ice data, including the number and thickness of pressure ridges, coordinated observing networks such as SIZONet can provide the means of reducing uncertainties inherent in individual datasets.

During spring 2013, we performed 500 MHz, helicopter-borne impulsive ground-penetrating radar surveys of several glaciers and glacier forelands in south-central Alaska, USA. These surveys were designed to obtain spatially distributed measurements of snow accumulation spanning a broad range of continental and maritime climatic zones. Visual assessment of radar images shows that data quality varied with the terrains and was optimal for snow that covered smooth glacier ice and firn, smooth debris-covered areas and moraines, freshwater lake and river ice, tundra, and taiga. Conversely, returns from the base of the snowpack were unrecognizable over rough debris-covered glacier termini, icefalls and some high-altitude accumulation basins. Optimal flying speed was 15-20ms–1 (30–40kt). At these speeds, which are two to three times faster than previously reported for such surveys, we could still identify snow-depth data with confidence, at a point spacing of ~1.5-2.0m. Data quality on glaciers decreased with increased air speed, though useful echoes from the base of the snowpack were still obtained at 40-45 ms–1 (87 kt; data point spacing of 6-8 m). Similar high-speed surveys over non-glacial terrains were unsuccessful, as basal reflections were no longer recognizable.

Isolated areas of high basal drag, or ‘sticky spots’, are important and poorly understood features in the force balance and dynamics of West Antarctic ice streams. Characterizing sticky spots formed by thin or drying subglacial till using ice-penetrating radar is theoretically possible, as high radar bed-returned power (BRP) is commonly related to an abundance of free water at the ice/bed interface, provided losses from englacial attenuation can be estimated. In this study we use airborne radar data collected over Evans Ice Stream to extract BRP profiles and test the sensitivity of BRP to the adopted englacial attenuation correction. We analyse 11 ~ 2 0 km profiles in four fast-flow areas where sticky spots have been inferred to exist on the basis of model and surface data inversions. In the majority of profiles we note that the increase in basal drag is accompanied by a decrease in BRP and suggest that this is evidence both for the presence of a sticky spot in those locations and that local variations in subglacial hydrology are responsible for their existence. A comparison is made between empirical and numerical modelling approaches for deriving englacial attenuation, and our findings generally support previous studies that advocate a modelling approach.

Despite their importance as freshwater reservoirs for downstream river systems, few glaciers in central Chile have been comprehensively surveyed. This study presents ground-penetrating radar (GPR) and field-based observations for characterizing the englacial and basal conditions of Glaciar Olivares Alfa (33°110 S, 70°130 W), central Chilean Andes. Using a 50 MHz radar mounted onto a helicopter platform, data were collected covering large portions of the glacier accumulation and ablation zones. The radar data revealed boundaries of a temperate-ice layer at the base of the eastern body of Glaciar Olivares Alfa which appears to be covered by colder ice that extends throughout large parts of the glacier. The thickness of the temperate ice layer is highly variable across the glacier, being on average 40% of the total ice thickness. Radar data analyses reveal regions of cold ice at the bottom/base of the glacier and also patterns of highly saturated sediments beneath the glacier. Using GPR data, this study represents the most exhaustive analysis of glacier ice structure performed in the central Chilean Andes. The results will enable improved estimations of the glacier’s mass balance and ice dynamics, helping us to understand its further development and its impact on water availability.

We analyze the internal stratigraphy in radio-echo sounding data of the northeast Greenland ice stream to infer past and present ice dynamics. In the upper reaches of the ice stream, we propose that shear-margin steady-state folds in internal reflecting horizons (IRHs) form due to the influence of ice flow over spatially varying basal lubrication. IRHs are generally lower in the ice stream than outside, likely because of greater basal melting in the ice stream from enhanced geothermal flux and heat of sliding. Strain-rate modeling of IRHs deposited during the Holocene indicates no recent major changes in ice-stream vigor or extent in this region. Downstream of our survey, IRHs are disrupted as the ice flows into a prominent overdeepening. When combined with additional data from other studies, these data suggest that upstream portions of the ice stream are controlled by variations in basal lubrication whereas downstream portions are confined by basal topography.

Satellite altimetric time series allow high-precision monitoring of ice-sheet mass balance. Understanding elevation changes in these regions is important because outlet glaciers along ice-sheet margins are critical in controlling flow of inland ice. Here we discuss a new airborne altimetry dataset collected as part of the ICECAP (International Collaborative Exploration of the Cryosphere by Airborne Profiling) project over East Antarctica. Using the ALAMO (Airborne Laser Altimeter with Mapping Optics) system of a scanning photon-counting lidar combined with a laser altimeter, we extend the 2003–09 surface elevation record of NASA’s ICESat satellite, by determining cross-track slope and thus independently correcting for ICESat’s cross-track pointing errors. In areas of high slope, cross-track errors result in measured elevation change that combines surface slope and the actual Δz/Δt signal. Slope corrections are particularly important in coastal ice streams, which often exhibit both rapidly changing elevations and high surface slopes. As a test case (assuming that surface slopes do not change significantly) we observe a lack of ice dynamic change at Cook Ice Shelf, while significant thinning occurred at Totten and Denman Glaciers during 2003–09.