Airborne fine-resolution UHF radar: an approach to the study of englacial reflections, firn compaction and ice attenuation rates

This is the published version. Copyright 2015 International Glaciological SocietyWe have built and operated an ultra-wideband UHF pulsed-chirp radar for measuring firn stratigraphy from airborne platforms over the ice sheets of Greenland and West Antarctica. Our analysis found a wide range of capabilities, including imaging of post firn–ice transition horizons and sounding of shallow glaciers and ice shelves. Imaging of horizons to depths exceeding 600 m was possible in the colder interior regions of the ice sheet, where scattering from the ice surface and inclusions was minimal. The radar's high sensitivity and large dynamic range point to loss tangent variations as the dominant mechanism for these englacial reflective horizons. The radar is capable of mapping interfaces with reflection coefficients as low as –80 dB near the firn–ice transition and as low as –64 dB at depths of 600 m. We found that firn horizon reflectivity strongly mirrored density variance, a result of the near-unity interfacial transmission coefficients. Zones with differing compaction mechanisms were also apparent in the data. We were able to sound many ice shelves and areas of shallow ice. We estimated ice attenuation rates for a few locations, and our attenuation estimates for the Ross Ice Shelf, West Antarctica, appear to agree well with earlier reported results

Radar Systems for Glaciology

This chapter deals with radar systems, measurements and instrumentation employed to study the internal core and bedrock of ice sheets in glaciology. The Earth's ice sheets are in Greenland and Antarctica. They cover about 10% of the land surface of the planet. The total accumulated ice comprises 90% of the global fresh water reserve. These ice sheets, associated with the ocean environment, provide a major heat sink which significantly modulates climate. Glaciology studies aim to understand the various process involved in the flow (dynamics), thermodynamics, and long-term behaviour of ice sheets. Studies of large ice masses are conducted in adverse environmental conditions (extreme cold, long periods of darkness). The development of remote sensing techniques have played an important role in obtaining useful results. The most widely used techniques are radar systems, employed since the 1950s in response to a need to provide a rapid and accurate method of measuring ice thickness. Year by year, polar research has become increasingly important because of global warming. Moreover, the discovery of numerous subglacial lake areas (water entrapped beneath the ice sheets) has attracted scientific interest in the possible existence of water circulation between lakes or beneath the ice (Kapitsa et al., 2006; Wingham et al., 2006; Bell et al., 2007). Recent studies in radar signal shape and amplitude could provide evidence of water circulation below the ice (Carter 2007, Oswald and Gogineni 2008). In this chapter the radar systems employed in glaciology, radio echo sounding (RES), are briefly described with some interesting results. RES are active remote sensing systems that utilize electromagnetic waves that penetrate the ice. They are used to obtain information about the electromagnetic properties of different interfaces (for example rock-ice, ice-water, seawater-ice) that reflect the incoming signal back to the radar. RES systems are characterized by a high energy (peak power from 10 W to 10 KW) variable transmitted pulse width (about from 0.5 ns to several microseconds) in order to investigate bedrock characteristics even in the thickest zones of the ice sheets (4755 m is the deepest ice thickness measured in Antarctica using a RES system). Changing the pulse length or the transmitted signal frequencies it is possible to investigate particular ice sheet details with different resolution. Long pulses allows transmission of higher power than short pulses, penetrating the thickest parts of the ice sheets but, as a consequence, resolution decreases. For example, the GPR system, commonly used in geophysics for rock, soil, ice, fresh water, pavement and structure characterization, employs a very short transmitted pulse (0.5 ns to 10 ns) that allow detailing of the shallow parts of an ice sheet (100-200 m in depth) (Reynolds 1997). Consequently, in recent years, GPR systems are also employed by explorers to find hidden crevasses on glaciers for safety. RES surveys have been widely employed in Antarctic ice sheet exploration and they are still an indispensable tool for mapping bedrock morphologies and properties of the last unexplored continent on Earth. The advantage of using these remote sensing techniques is that they allow large areas to be covered, in good detail and in short times using platforms like aeroplanes and surface vehicles

A wideband radar for high-resolution mapping of near-surface internal layers in glacial ice

©2004 IEEE. Personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional purposes or for creating new collective works for resale or redistribution to servers or lists, or to reuse any copyrighted component of this work in other works must be obtained from the IEEE.Snow accumulation rate is an important parameter in determining the mass balance of polar ice sheets. Accumulation rate is currently determined by analyzing ice cores and snow pits. Inadequate sampling of the spatial variations in the ice sheet accumulation has resulted in accumulation rate uncertainties as large as 24%. We designed and developed a 600-900-MHz airborne radar system for high-resolution mapping of the near-surface internal layers for estimating the accumulation rate of polar ice sheets. Our radar system can provide improved spatial and temporal coverage by mapping a continuous profile of the isochronous layers in the ice sheet. During the 2002 field season in Greenland, we successfully mapped the near-surface layers to a depth of 200 m in the dry-snow zone, 120 m in the percolation zone, and 20 m in the melt zone. We determined the water equivalent accumulation rate at the NASA-U_1 site to be 34.9 +/- 5.1 cm/year from 1964 to 1992. This is in close agreement with the ice-core derived accumulation rate of 34.6 cm/year for the same period

Origin of englacial stratigraphy at three deep ice core sites of the Greenland Ice Sheet by synthetic radar modelling

During the past 20 years, multi-channel radar emerged as a key tool for deciphering an ice sheet's internal architecture. To assign ages to radar reflections and connect them over large areas in the ice sheet, the layer genesis has to be understood on a microphysical scale. Synthetic radar trace modelling based on the dielectric profile of ice cores allows for the assignation of observed physical properties’ variations on the decimetre scale to radar reflectors extending from the coring site to a regional or even whole-ice-sheet scale. In this paper we rely on the available dielectric profiling data of the northern Greenland deep ice cores: NGRIP, NEEM and EGRIP. The three records are well suited for assigning an age model to the stratigraphic radar-mapped layers, and linking up the reflector properties to observations in the cores. Our modelling results show that the internal reflections are mainly due to conductivity changes. Furthermore, we deduce fabric characteristics at the EGRIP drill site from two-way-travel-time differences of along and across-flow polarized radarwave reflections of selected horizons (below 980 m). These indicate in deeper parts of the ice column an across-flow concentrated c-axis fabric

The precision of radar-derived subglacial bed topography: a case study from Pine Island Glacier, Antarctica

Recent advances in the measurement of bedforms beneath active ice streams have been made using ground-based grid profiling using impulse radar systems operating with centre frequencies in the 3–5 MHz range. Surveys of Rutford Ice Stream and Pine Island Glacier have shown that features such as mega-scale glacial lineations with topographic relief of as little as 3 m can be traced for many kilometres downstream under more than 2 km of fast-moving ice. In the discussion of these data, it is often asked ‘How is it possible to map such fine-scale topography with such a low-frequency radar’. In answering that question, the key point is the distinction between the precision of a radar range measurement to a single, isolated reflective interface and the ability to resolve the presence of two closely-spaced interfaces of similar reflectivity (commonly referred to as the vertical resolution). This paper will discuss and illustrate this distinction and use the case study of data acquired over Pine Island Glacier to examine the limits of precision of the radar range measurement

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

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

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

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

Five decades of radioglaciology

Radar 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 the Annals of Glaciology on ‘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-glaciological Annals of Glaciology issue of 2014, the context of their history and future prospects

Bed topography of Jakobshavn Isbrae, Greenland, and Byrd Glacier, Antarctica

This is the published version. Copyright 2015 International Glaciological SocietyThis paper presents the bed topography of Jakobshavn Isbrae, Greenland, and Byrd Glacier, Antarctica, derived from sounding these glaciers with high-sensitivity radars. To understand the processes causing the speed-up and retreat of outlet glaciers, and to enable the development of next-generation ice-sheet models, we need information on bed topography and basal conditions. To this end, we performed measurements with the progressively improved Multichannel Coherent Radar Depth Sounder/Imager (MCoRDS/I). We processed the data from each antenna-array element using synthetic aperture radar algorithms to improve radar sensitivity and reduce along-track surface clutter. We then applied array and image-processing algorithms to extract the weak bed echoes buried in off-vertical scatter (cross-track surface clutter). At Jakobshavn Isbrae, we observed 2.7 km thick ice ∼30 km upstream of the calving front and ∼850 m thick ice at the calving front. We also observed echoes from multiple interfaces near the bed. We applied the MUSIC algorithm to the data to derive the direction of arrival of the signals. This analysis revealed that clutter is dominated by the ice surface at Jakobshavn Isbrae. At Byrd Glacier, we found ∼3.62 km thick ice, as well as a subglacial trench ∼3.05 km below sea level. We used ice thickness information derived from radar data in conjunction with surface elevation data to generate bed maps for these two critical glaciers. The performance of current radars must be improved further by ∼15 dB to fully sound the deepest part of Byrd Glacier. Unmanned aerial systems equipped with radars that can be flown over lines spaced as close as 5 m apart in the cross-track direction to synthesize a two-dimensional aperture would be ideal for collecting fine-resolution data over glaciers like Jakobshavn near their grounding lines