Kirsty Langley

Use of C-Band Ground Penetrating Radar to Determine Backscatter Sources Within Glaciers

The question of penetration of synthetic aperture radar (SAR) signals at C-band frequency into polar glaciers is addressed by comparing ground penetrating radar (GPR) and SAR backscatter signatures. Profiles of the Kongsvegen glacier, Svalbard, were obtained with a C-band GPR. The received signal is converted to the equivalent radar cross section using the standard radar equation, thus mapping the effective scattering sources within the glacier at this frequency. The depth of the observed scattering sources is greatest in the superimposed ice where layers are clearly seen to a depth of approximately 14 m. The very high scattering properties of the upper firn layers preclude layers deeper than approximately 6 m from being imaged. Integrating the radar cross sections over the depth gives a single backscatter value that we compare with the backscatter coefficient (scattering cross section per unit area) of the processed SAR data for the same profile. The comparison indicates that for coincidentally acquired GPR and SAR data, the radar cross section measured by the GPR does represent the features that contribute to the SAR signal

From Glacier Facies to SAR Backscatter Zones via GPR

We present a comparison between data acquired with frequency-modulated ground-penetrating radar (GPR) and satellite synthetic aperture radar (SAR). Both radars are polarimetric and operate at a center frequency of 5.3 GHz. The field site is the polythermal glacier Kongsvegen, Svalbard. Along glacier GPR profiles cover the ablation area and the accumulation area, where the latter consists of superimposed ice (SI) and firn. The glacier facies are clearly identifiable on the GPR profiles, although we show that the copolarized response is better for distinguishing different ice zones, whereas the SI-firn boundary is most obvious in the cross-polarized response. A calibrated backscatter coefficient is calculated for the GPR data and compared with the SAR backscatter coefficient. The SAR zones are in very good agreement with the GPR-derived glacier facies. We show that, in the ablation area, the SAR response is dominated by backscatter from the previous summer surface. In the SI and firn areas, it is dominated by sources below the previous summer surface.

Complex network of channels beneath an Antarctic ice shelf

Ice shelves play an important role in stabilizing the interior grounded ice of the large ice sheets. The thinning of major ice shelves observed in recent years, possibly in connection to warmer ocean waters coming into contact with the ice-shelf base, has focused attention on the ice-ocean interface. Here we reveal a complex network of sub ice-shelf channels under the Fimbul Ice Shelf, Antarctica, mapped using ground-penetrating radar over a 100 km2 grid. The channels are 300–500 m wide and 50 m high, among the narrowest of any reported. Observing narrow channels beneath an ice shelf that is mainly surrounded by cold ocean waters, with temperatures close to the surface freezing point, shows that channelized basal melting is not restricted to rapidly melting ice shelves, indicating that spatial melt patterns around Antarctica are likely to vary on scales that are not yet incorporated in ice-ocean models.

Observation-Based Estimates of Global Glacier Mass Change and Its Contribution to Sea-Level Change

Glaciers have strongly contributed to sea-level rise during the past century and will continue to be an important part of the sea-level budget during the twenty-first century. Here, we review the progress in estimating global glacier mass change from in situ measurements of mass and length changes, remote sensing methods, and mass balance modeling driven by climate observations. For the period before the onset of satellite observations, different strategies to overcome the uncertainty associated with monitoring only a small sample of the world’s glaciers have been developed. These methods now yield estimates generally reconcilable with each other within their respective uncertainty margins. Whereas this is also the case for the recent decades, the greatly increased number of estimates obtained from remote sensing reveals that gravimetry-based methods typically arrive at lower mass loss estimates than the other methods. We suggest that strategies for better interconnecting the different methods are needed to ensure progress and to increase the temporal and spatial detail of reliable glacier mass change estimates.

Sources of backscatter at 5.3 GHz from a superimposed ice and firn area revealed by multi-frequency GPR and cores

We investigate the major sources of backscatter at 5.3 GHz, within the superimposed ice and firn areas of a polythermal glacier. Two ground-penetrating radar systems, an 800 MHz impulse system and a polarimetric 5.3 GHz frequency-modulated continuous-wave system, are used to acquire along-glacier profiles in the accumulation area of Kongsvegen, Svalbard. The 800 MHz response is used to map reflection horizons in the glacier. Using cores from the superimposed ice and firn areas, the causes of these reflection horizons, in terms of snow, firn and ice layers, are investigated. Superimposing the reflection horizons on the co-polarized and cross-polarized 5.3 GHz profile, we are able to determine how the 5.3 GHz frequency responds to the different media. Scattering at rough interfaces and volume scattering occur in the superimposed ice area and are apparently caused by air-bubble number, size and distribution. In the firn the strongest return originates from below the previous summer surface, consistent with previous findings. At approximately the same depth, strong incoherent scattering begins. The rapid decrease in coherent reflections indicates the significance of scattering in the firn.

In situ measurements of snow surface roughness using a laser profiler

The snow surface roughness at centimetre and millimetre scales is an important parameter related to wind transport, snowdrifts, snowfall, snowmelt and snow grain size. Knowledge of the snow surface roughness is also of high interest for analyzing the signal from radar sensors such as SAR, altimeters and scatterometers. Unfortunately, this parameter has seldom been measured over snow surfaces. The techniques used to measure the roughness of other surfaces, such as agricultural or sand soils, are difficult to implement in polar regions because of the harsh climatic conditions. In this paper we develop a device based on a laser profiler coupled with a GPS receiver on board a snowmobile. This instrumentation was tested successfully in midre Lovenbreen, Svalbard, in April 2006. It allowed us to generate profiles of 3 km sections of the snow-covered glacier surface. Because of the motion of the snowmobile, the roughness signal is mixed with the snowmobile signal. We use a distance/frequency analysis (the empirical mode decomposition) to filter the signal. This method allows us to recover the snow surface structures of wavelengths between 4 and 50 cm with amplitudes of >1 mm. Finally, the roughness parameters of snow surfaces are retrieved. The snow surface roughness is found to be dependent on the scales of the observations. The retrieved RMS of the height distribution is found to vary between 0.5 and 9.2 mm, and the correlation length is found to be between 0.6 and 46 cm. This range of measurements is particularly well adapted to the analysis of GHz radar response on snow surfaces.

Low melt rates with seasonal variability at the base of Fimbul Ice Shelf, East Antarctica, revealed by in situ interferometric radar measurements

Basal melt is a major cause of ice shelf thinning affecting the stability of the ice shelf and reducing its buttressing effect on the inland ice. The Fimbul ice shelf (FIS) in Dronning Maud Land (DML), East Antarctica, is fed by the fast-flowing Jutulstraumen glacier, responsible for 10% of ice discharge from the DML sector of the ice sheet. Current estimates of the basal melt rates of the FIS come from regional ocean models, autosub measurements, and satellite observations, which vary considerably. This discrepancy hampers evaluation of the stability of the Jutulstraumen catchment. Here, we present estimates of basal melt rates of the FIS using ground-based interferometric radar. We find a low average basal melt rate on the order of 1 m/yr, with the highest rates located at the ice shelf front, which extends beyond the continental shelf break. Furthermore, our results provide evidence for a significant seasonal variability.

Analysis and Classification of high Arctic Glaciers with ASAR Data

We apply a polarimetric classification scheme to Envisat Alternating Polarisation mode ASAR images and compare to ground truth data. By analysing images with a range of acquisition conditions and comparing the classification accuracy with the ground truth data, we investigate the influence of the acquisition conditions for glacier facies discrimination. We find that the main influence is the image swath angle, which affects the image pixel geometry. There was no obvious preference for the different polarisation channels, although the dual polarisation classifications were consistently better than single polarisation classifications.

Validation of CryoSat-2 SARIn Data over Austfonna Ice Cap Using Airborne Laser Scanner Measurements

The study presented here is focused on the assessment of surface elevations derived from CryoSat-2 SARIn level 1b data over the Austfonna ice cap, Svalbard, in 2016. The processing chain that must be applied to the CryoSat-2 waveforms to derive heights is non-trivial, and consists of multiple steps, all requiring subjective choices of methods such as the choice of retracker, geo-relocation, and outlier rejection. Here, we compare six CryoSat-2 level-2 type data sets of surface elevations derived using different SARIn processing chains. These data sets are validated against surface elevation data collected from an airborne laser scanner, during a dedicated CryoSat validation experiment field campaign carried out in April 2016. The flight pattern of the airborne campaign was designed so that elevations were measured in a grid pattern rather than along single lines, as has previously been the standard procedure. The flight grid pattern was chosen to optimize the comparison with the CryoSat-2 SARIn elevation data, the location of which can deviate from nadir by several kilometers due to topography within the satellite footprint. The processing chains behind the six data sets include different outlier/error rejection approaches, and do not produce the same number of data points in our region of interest. To make a consistent analysis, we provide statistics from the validation of both the full data sets from each processing chain, and on only those data that all the six data sets provide a geo-located elevation estimate for. We find that the CryoSat-2 data sets that agree best with the validation data are those derived from dedicated land ice processing schemes. This study may serve as a benchmark for future CryoSat-2 retracker developments, and the evaluation software and data set are made publicly available.