A matched-filter-based reverse-time migration algorithm for ground-penetrating radar data

Ground-penetrating radar (GPR) is a remote sensing technique used to obtain information on subsurface features from data collected over the surface. The process of collecting data may be viewed as mapping from the object space to an image space. Since most GPRs use broad beamwidth antennas, the energy reflected from a buried structure is recorded over a large lateral aperture in the image space. Migration algorithms are used to reconstruct an accurate scattering map by refocusing the recorded scattering events to their true spatial locations through a backpropagation process. The goal of this paper is to present a pair of finite-difference time-domain (FDTD) reverse-time migration algorithms for CPR data processing. Linear inverse scattering theory is used to develop a matched-filter response for the GPR problem. The reverse-time migration algorithms, developed for both bistatic and monostatic antenna configurations, are implemented via FDTD in the object space. Several examples are presented

Validation of Airborne FMCW Radar Measurements of Snow Thickness Over Sea Ice in Antarctica

Antarctic sea ice and its snow cover are integral components of the global climate system, yet many aspects of their vertical dimensions are poorly understood, making their representation in global climate models poor. Remote sensing is the key to monitoring the dynamic nature of sea ice and its snow cover. Reliable and accurate snow thickness data are currently a highly sought after data product. Remotely sensed snow thickness measurements can provide an indication of precipitation levels, predicted to increase with effects of climate change in the polar regions. Airborne techniques provide a means for regional-scale estimation of snow depth and distribution. Accurate regional-scale snow thickness data will also facilitate an increase in the accuracy of sea ice thickness retrieval from satellite altimeter freeboard estimates. The airborne data sets are easier to validate with in situ measurements and are better suited to validating satellite algorithms when compared with in situ techniques. This is primarily due to two factors: better chance of getting coincident in situ and airborne data sets and the tractability of comparison between an in situ data set and the airborne data set averaged over the footprint of the antennas. A 28-GHz frequency modulated continuous wave (FMCW) radar loaned by the Center for Remote Sensing of Ice Sheets to the Australian Antarctic Division is used to measure snow thickness over sea ice in East Antarctica. Provided with the radar design parameters, the expected performance parameters of the radar are summarized. The necessary conditions for unambiguous identification of the airsnow and snowice layers for the radar are presented. Roughnesses of the snow and ice surfaces are found to be dominant determinants in the effectiveness of layer identification for this radar. Finally, this paper presents the first in situ validated snow thickness estimates over sea ice in Antarctica derived from an FMCW radar on a helicopterborne platform

Radar mapping of Isunnguata Sermia, Greenland

This is the published version. Copyright 2013 International Glaciological SocietyIce thickness estimates using advanced nadir sounding and tomographic radar processing techniques are compared and combined in a study of Isunnguata Sermia glacier, Greenland. Using an ensemble of Operation IceBridge flight lines spaced at 500 m intervals and running approximately along the flow direction, we find there is a statistically excellent comparison between subglacial terrains derived from two-dimensional tomography and gridded nadir sounding. Analysis shows that tomographic data better capture short wavelength (1–2 km) patterns in basal terrain, but interpolated nadir sounding data yield more spatially extensive and continuous coverage across the glacier, especially in deep subglacial troughs. Using derived surface and basal topography maps, we find that driving stress and measured and modeled surface velocity comparisons indicate that basal sliding is an important component of the glacier motion, but is also only weakly coupled to the detailed bed topography save for the deepest troughs. As might be expected for this land-terminating, relatively slow-moving glacier, the subglacial and proglacial topography is similar, suggesting the erosional processes acting on the modern glacier bed once helped sculpt the now exposed land

CReSIS airborne radars and platforms for ice and snow sounding

This work is licensed under a Creative Commons Attribution 4.0 International License.This paper provides an update and overview of the Center for Remote Sensing of Ice Sheets (CReSIS) radars and platforms, including representative results from these systems. CReSIS radar systems operate over a frequency range of 14–38 GHz. Each radar system's specific frequency band is driven by the required depth of signal penetration, measurement resolution, allocated frequency spectra, and antenna operating frequencies (often influenced by aircraft integration). We also highlight recent system advancements and future work, including (1) increasing system bandwidth; (2) miniaturizing radar hardware; and (3) increasing sensitivity. For platform development, we are developing smaller, easier to operate and less expensive unmanned aerial systems. Next-generation platforms will further expand accessibility to scientists with vertical takeoff and landing capabilities

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

An Ultra-Wideband, Microwave Radar for Measuring Snow Thickness on Sea Ice and Mapping Near-Surface Internal Layers in Polar Firn

Sea ice is generally covered with snow, which can vary in thickness from a few centimeters to >1 m. Snow cover acts as a thermal insulator modulating the heat exchange between the ocean and the atmosphere, and it impacts sea-ice growth rates and overall thickness, a key indicator of climate change in polar regions. Snow depth is required to estimate sea-ice thickness using freeboard measurements made with satellite altimeters. The snow cover also acts as a mechanical load that depresses ice freeboard (snow and ice above sea level). Freeboard depression can result in flooding of the snow/ice interface and the formation of a thick slush layer, particularly in the Antarctic sea-ice cover. The Center for Remote Sensing of Ice Sheets (CReSIS) has developed an ultra-wideband, microwave radar capable of operation on long-endurance aircraft to characterize the thickness of snow over sea ice. The low-power, 100mW signal is swept from 2 to 8GHz allowing the air/snow and snow/ ice interfaces to be mapped with 5 c range resolution in snow; this is an improvement over the original system that worked from 2 to 6.5 GHz. From 2009 to 2012, CReSIS successfully operated the radar on the NASA P-3B and DC-8 aircraft to collect data on snow-covered sea ice in the Arctic and Antarctic for NASA Operation IceBridge. The radar was found capable of snow depth retrievals ranging from 10cm to >1 m. We also demonstrated that this radar can be used to map near-surface internal layers in polar firn with fine range resolution. Here we describe the instrument design, characteristics and performance of the radar

The trough-system algorithm and its application to spatial modeling of Greenland subglacial topography

This is the published version. Copyright 2014 Herzfeld et al.Dynamic ice-sheet models are used to assess the contribution of mass loss from the Greenland ice sheet to sea-level rise. Mass transfer from ice sheet to ocean is in a large part through outlet glaciers. Bed topography plays an important role in ice dynamics, since the acceleration from the slow-moving inland ice to an ice stream is in many cases caused by the existence of a subglacial trough or trough system. Problems are that most subglacial troughs are features of a scale not resolved in most ice-sheet models and that radar measurements of subglacial topography do not always reach the bottoms of narrow troughs. The trough-system algorithm introduced here employs mathematical morphology and algebraic topology to correctly represent subscale features in a topographic generalization, so the effects of troughs on ice flow are retained in ice-dynamic models. The algorithm is applied to derive a spatial elevation model of Greenland subglacial topography, integrating recently collected radar measurements (CReSIS data) of the Jakobshavn Isbræ, Helheim, Kangerdlussuaq and Petermann glacier regions. The resultant JakHelKanPet digital elevation model has been applied in dynamic ice-sheet modeling and sea-level-rise assessment

An ultra-wideband, microwave radar for measuring snow thickness on sea ice and mapping near-surface internal layers in polar firn

This is the published version. Copyright 2013 International Glaciological SocietySea ice is generally covered with snow, which can vary in thickness from a few centimeters to >1 m. Snow cover acts as a thermal insulator modulating the heat exchange between the ocean and the atmosphere, and it impacts sea-ice growth rates and overall thickness, a key indicator of climate change in polar regions. Snow depth is required to estimate sea-ice thickness using freeboard measurements made with satellite altimeters. The snow cover also acts as a mechanical load that depresses ice freeboard (snow and ice above sea level). Freeboard depression can result in flooding of the snow/ice interface and the formation of a thick slush layer, particularly in the Antarctic sea-ice cover. The Center for Remote Sensing of Ice Sheets (CReSIS) has developed an ultra-wideband, microwave radar capable of operation on long-endurance aircraft to characterize the thickness of snow over sea ice. The low-power, 100mW signal is swept from 2 to 8GHz allowing the air/snow and snow/ ice interfaces to be mapped with 5 cm range resolution in snow; this is an improvement over the original system that worked from 2 to 6.5 GHz. From 2009 to 2012, CReSIS successfully operated the radar on the NASA P-3B and DC-8 aircraft to collect data on snow-covered sea ice in the Arctic and Antarctic for NASA Operation IceBridge. The radar was found capable of snow depth retrievals ranging from 10cm to >1 m. We also demonstrated that this radar can be used to map near-surface internal layers in polar firn with fine range resolution. Here we describe the instrument design, characteristics and performance of the radar

Measuring Height Change Around the Periphery of the Greenland Ice Sheet With Radar Altimetry

This work is licensed under a Creative Commons Attribution 4.0 International License.Ice loss measurements around the periphery of the Greenland Ice Sheet can provide key information on the response to climate change. Here we use the excellent spatial and temporal coverage provided by the European Space Agency (ESA) CryoSat satellite, together with NASA airborne Operation IceBridge and automatic weather station data, to study the influence of changing conditions on the bias between the height estimated by the satellite radar altimeter and the ice sheet surface. Surface and near-surface conditions on the ice sheet periphery change with season and geographic position in a way that affects the returned altimeter waveform and can therefore affect the estimate of the surface height derived from the waveform. Notwithstanding the possibility of a varying bias between the derived and real surface, for the lower accumulation regions in the western and northern ice sheet periphery (<∼1 m snow accumulation yearly) we show that the CryoSat altimeter can measure height change throughout the year, including that associated with ice dynamics, summer melt and winter accumulation. Further, over the 9-year CryoSat lifetime it is also possible to relate height change to change in speed of large outlet glaciers, for example, there is significant height loss upstream of two branches of the Upernavik glacier in NW Greenland that increased in speed during this time, but much less height loss over a third branch that slowed in the same time period. In contrast to the west and north, winter snow accumulation in the south-east periphery can be 2–3 m and the average altimeter height for this area can decrease by up to 2 m during the fall and winter when the change in the surface elevation is much smaller. We show that vertical downward movement of the dense layer from the last summer melt, coupled with overlying dry snow, is responsible for the anomalous altimeter height change. However, it is still possible to estimate year-to-year height change measurements in this area by using data from the late-summer to early fall when surface returns dominate the altimeter signal