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

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

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

Radar sounder evidence of thick, porous sediments in Meridiani Planum and implications for ice‐filled deposits on Mars

Meridiani Planum is one of the most intensely studied regions on Mars, yet little is known about the physical properties of the deposits below those examined by the Opportunity rover. We report the detection of subsurface echoes within the Meridiani Planum deposits from data obtained by the Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) instrument. The delay time between the surface and subsurface returns is indicative of materials with a real dielectric constant of 3.6 ± 0.6. The real dielectric constant is strongly modulated by bulk density. Newly derived compaction relationships for Mars indicate that the relatively low dielectric constant of the Meridiani Planum deposits is consistent with a thick layer of ice‐free, porous, basaltic sand. The unique physiographic and hydrologic setting of Meridiani Planum may have provided an ideal sediment trap for eolian sands. The relatively low gravity and the cold, dry climate that has dominated Mars for billions of years may have allowed thick eolian sand deposits to remain porous and only weakly indurated. Minimally compacted sedimentary deposits may offer a possible explanation for other nonpolar region units with low apparent bulk dielectric constants

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

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

A Comparison of Snow Depth on Sea Ice Retrievals Using Airborne Altimeters and an AMSR-E Simulator

A comparison of snow depths on sea ice was made using airborne altimeters and an Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E) simulator. The data were collected during the March 2006 National Aeronautics and Space Administration (NASA) Arctic field campaign utilizing the NASA P-3B aircraft. The campaign consisted of an initial series of coordinated surface and aircraft measurements over Elson Lagoon, Alaska and adjacent seas followed by a series of large-scale (100 km ? 50 km) coordinated aircraft and AMSR-E snow depth measurements over portions of the Chukchi and Beaufort seas. This paper focuses on the latter part of the campaign. The P-3B aircraft carried the University of Colorado Polarimetric Scanning Radiometer (PSR-A), the NASA Wallops Airborne Topographic Mapper (ATM) lidar altimeter, and the University of Kansas Delay-Doppler (D2P) radar altimeter. The PSR-A was used as an AMSR-E simulator, whereas the ATM and D2P altimeters were used in combination to provide an independent estimate of snow depth. Results of a comparison between the altimeter-derived snow depths and the equivalent AMSR-E snow depths using PSR-A brightness temperatures calibrated relative to AMSR-E are presented. Data collected over a frozen coastal polynya were used to intercalibrate the ATM and D2P altimeters before estimating an altimeter snow depth. Results show that the mean difference between the PSR and altimeter snow depths is -2.4 cm (PSR minus altimeter) with a standard deviation of 7.7 cm. The RMS difference is 8.0 cm. The overall correlation between the two snow depth data sets is 0.59

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

MARSIS Subsurface Radar Sounding of Medusae Fossae Formation Deposits on Mars: Ice Rich or Ice Poor, That is the Question

Newly acquired MARSIS radar sounder data shows evidence of layering in all major units of the MFF deposits, establishing another key similarity to PLD