skip to main content


Title: Layer optimized SAR processing and slope estimation in radar sounder data
Abstract Englacial layers in Antarctica and Greenland are indicators of the dynamic, rheological and subglacial configuration of the ice sheets. Airborne radar sounder data is the primary remote sensing solution for directly observing englacial layers and structures at the glacier-catchment to ice-sheet scale. However, when traditional along-track synthetic aperture radar (SAR) processing is applied, steep layers can disappear, limiting the detectability and interpretability of englacial layer geometry. This study provides a reconstruction algorithm to address the problem of destructive phase interference during the radargram formation. We develop and apply a novel SAR processor optimized for layer detection that enhances the Signal-to-Noise ratio (SNR) of specular reflectors. The algorithm also enables the automatic estimation of layer slope. We demonstrate the algorithm using data acquired at the Institute Ice Stream, West Antarctica.  more » « less
Award ID(s):
1745137
NSF-PAR ID:
10143101
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of Glaciology
Volume:
65
Issue:
254
ISSN:
0022-1430
Page Range / eLocation ID:
983 to 988
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract Radar surveys across ice sheets typically measure numerous englacial layers that can often be regarded as isochrones. Such layers are valuable for extrapolating age–depth relationships away from ice-core locations, reconstructing palaeoaccumulation variability, and investigating past ice-sheet dynamics. However, the use of englacial layers in Antarctica has been hampered by underdeveloped techniques for characterising layer continuity and geometry over large distances, with techniques developed independently and little opportunity for inter-comparison of results. In this paper, we present a methodology to assess the performance of automated layer-tracking and layer-dip-estimation algorithms through their ability to propagate a correct age–depth model. We use this to assess isochrone-tracking techniques applied to two test case datasets, selected from CreSIS MCoRDS data over Antarctica from a range of environments including low-dip, continuous layers and layers with terminations. We find that dip-estimation techniques are generally successful in tracking englacial dip but break down in the upper and lower regions of the ice sheet. The results of testing two previously published layer-tracking algorithms show that further development is required to attain a good constraint of age–depth relationships away from dated ice cores. We recommend that auto-tracking techniques focus on improved linking of picked stratigraphy across signal disruptions to enable accurate determination of the Antarctic-wide age–depth structure. 
    more » « less
  2. null (Ed.)
    Abstract Basal units – visibly distinct englacial structures near the ice-bed interface – warrant investigation for a number of reasons. Many are of unknown composition and origin, characteristics that could provide substantial insight into subglacial processes and ice-sheet history. Their significance, moreover, is not limited to near-bed depths; these units appear to dramatically influence the flow of surrounding ice. In order to enable improved characterization of these features, we develop and apply an algorithm that allows for the automatic detection of basal units. We use a tunable layer-optimized SAR processor to distinguish these structures from the bed, isochronous englacial layers and the ice-sheet surface, presenting a conceptual framework for the use of radio-echo character in the identification of ice-sheet features. We also outline a method by which our processor could be used to place observational constraints on basal units’ configuration, composition and provenance. 
    more » « less
  3. Abstract 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. 
    more » « less
  4. Abstract

    The Amundsen Sea Embayment of the West Antarctic Ice Sheet contains Thwaites and Pine Island Glaciers, two of the most rapidly changing glaciers in Antarctica. To date, Pine Island and Thwaites Glaciers have only been observed by independent airborne radar sounding surveys, but a combined cross‐basin analysis that investigates the basal conditions across the Pine Island‐Thwaites Glaciers boundary has not been performed. Here, we combine two radar surveys and correct for their differences in system parameters to produce unified englacial attenuation and basal relative reflectivity maps spanning both Pine Island and Thwaites Glaciers. Relative reflectivities range from −24.8 to +37.4 dB with the highest values beneath fast‐flowing ice at the ice sheet margin. By comparing our reflectivity results with previously derived radar specularity and trailing bed echoes at Thwaites Glacier, we find a highly diverse subglacial landscape and hydrologic conditions that evolve along‐flow. Together, these findings highlight the potential for joint airborne radar analysis with ground‐based seismic and geomorphological observations to understand variations in the bed properties and cross‐catchment interactions of ice streams and outlet glaciers.

     
    more » « less
  5. IODP Expedition 379 deep-sea drilling in 2019 (Gohl et al. 2021, doi:10.14379/iodp.proc.379.2021), offered an opportunity to obtain chronostratigraphic control for seismic reflection data for Amundsen Sea shelf and slope deposits that record Miocene to Present fluctuations in volume of the West Antarctic ice sheet. Here we report the age and interpret the provenance of a volcanic ash bed recovered at/near the Plio-Pleistocene boundary at 31.51 meters below sea level in Hole U1533B and 33.94 mbsf in Hole U1533D. With distinctive geochemistry and inferred wide regional distribution, the bed may serve as a reliable age marker. In Hole 1533B, the fresh tephra forms a discrete layer interstratified within uniform brown marine mud. The layer has a sharp base and upper boundary that is gradational over 5 cm into overlying mud. Color reflectance and density data aided identification of the tephra horizon (diffuse) in Hole 1533D, ~1000m away. A possible on-land source for ash is the Miocene to Pleistocene Marie Byrd Land volcanic province, comprising 18 large alkaline volcanoes dominated by effusive lavas. Products of pyroclastic eruptions are uncommon, mainly occurring as distal englacial, and probably marine, tephra. We undertook an offshore-onshore comparison by first characterizing samples of Site U1533 tephra from a petrographic and geochemical standpoint, using thin section observations, EMPA-WDS glass compositions, and 40Ar/39Ar dating. We then identified onshore exposures with similar characteristics. The offshore tephra are composed of coarse (50-300µm) cuspate glass shards with elongated vesicles. The glass composition is rhyolite, with 75-79wt.% SiO2, ~4wt.% FeO and 0.0wt.% MgO. Single-crystal feldspar 40Ar/39Ar dates are 2.55±0.12 and 2.92±0.02 Ma for U1533B and 2.87 ±0.45 Ma for U1533D. The geochemistry, shard morphology, discrete bed expression, and lateral continuity between Holes U1533B-U1533D indicate that the rhyolite tephra formed as airfall settled to the deep seabed. The ca. 2.55 Ma age based on youngest feldspar grains differs slightly from the 2.1 to 2.2 Ma result obtained from in-progress core bio-magnetostratigraphy. Rare exposures of rhyolite are found in the Chang Peak/Mt. Waesche centers, 1080 km from Site U1533. We obtained pumice sample MB.7.3 (prior-published age of 1.6±0.2 Ma), which displays elevated FeO and F content, and MB.8.1, a specimen of porphyritic cryptocrystalline lava. Single-crystal sanidine 40Ar/39Ar dates are 1.315±0.007 Ma (MB.7.3) and 1.385±0.003 Ma (MB.8.1). Site U1533 samples share a geochemical affinity with these on-land rhyolites, expressed as similar SiO2, CaO, TiO2, MgO and FeO content, suggesting an origin for Site U1533 tephra in the Chang-Waesche volcanoes. A possible explanation for the distinctly greater age, and observed contrasts in Al2O3, Na2O and F percentages, is that Site U1533 tephra are older and erupted from a source entirely concealed beneath subsequent eruptions and the ice sheet. Our results suggest that rhyolite volcanism initiated earlier, was of longer duration than previously known (2.92 to 1.315 Ma), and dispersed tephra far offshore. The finding is significant because ash and aerosols produced by large eruptions may influence regional climate. Antarctica cooled significantly and ice sheets expanded in latest Pliocene time (McKay et al. 2012, doi:10.1073/pnas.1112248109). 
    more » « less