scholarly journals Interferometric discrimination of cross-track bed clutter in ice-penetrating radar sounding data

2020 ◽  
Vol 61 (81) ◽  
pp. 68-73 ◽  
Author(s):  
Kirk M. Scanlan ◽  
Anja Rutishauser ◽  
Duncan A. Young ◽  
Donald D. Blankenship

AbstractThe interpretations of relevant interfaces (i.e. the surface and bed) in radar sounding datasets over glaciers and ice sheets are primary boundary conditions in a variety of climate studies and particularly subglacial water routing models. It is therefore necessary to ensure these interpretations are consistent and not affected by cross-track clutter. For the surface interface, interferometry and a family of methods relying on digital elevation models have been used to successfully discriminate cross-track surface clutter. Here we present how interferometry can be applied to the problem of basal clutter from cross-track bed topography. Our approach is based on a comparison of the differential phases of ambiguous reflectors that may represent bed clutter and the differential phase of a reflector in an adjacent area that appears unaffected by basal clutter. The reflector yielding the smallest interferometric phase difference relative to the unambiguous bed reflector is considered to represent its consistent continuation. We successfully demonstrate our approach using 60 MHz center frequency MARFA data collected over Devon Ice Cap in the Canadian Arctic. Finally, we investigate the effects of clutter-affected and interferometry-corrected bed interpretations on ice layer thickness estimates, basal hydraulic head gradients and the potential extent of inferred subglacial water bodies.

2013 ◽  
Vol 7 (6) ◽  
pp. 1857-1867 ◽  
Author(s):  
L. Gray ◽  
D. Burgess ◽  
L. Copland ◽  
R. Cullen ◽  
N. Galin ◽  
...  

Abstract. We have derived digital elevation models (DEMs) over the western part of the Devon Ice Cap in Nunavut, Canada, using "swath processing" of interferometric data collected by Cryosat between February 2011 and January 2012. With the standard ESA (European Space Agency) SARIn (synthetic aperture radar interferometry) level 2 (L2) data product, the interferometric mode is used to map the cross-track position and elevation of the "point-of-closest-approach" (POCA) in sloping glacial terrain. However, in this work we explore the extent to which the phase of the returns in the intermediate L1b product can also be used to map the heights of time-delayed footprints beyond the POCA. We show that there is a range of average cross-track slopes (~ 0.5 to ~ 2°) for which the returns will be dominated by those beneath the satellite in the main beam of the antenna so that the resulting interferometric phase allows mapping of heights in the delayed range window beyond the POCA. In this way a swath of elevation data is mapped, allowing the creation of DEMs from a sequence of L1b SARIn Cryosat data takes. Comparison of the Devon results with airborne scanning laser data showed a mean difference of order 1 m with a standard deviation of about 1 m. The limitations of swath processing, which generates almost 2 orders of magnitude more data than traditional radar altimetry, are explored through simulation, and the strengths and weaknesses of the technique are discussed.


2013 ◽  
Vol 7 (3) ◽  
pp. 3133-3162 ◽  
Author(s):  
L. Gray ◽  
D. Burgess ◽  
L. Copland ◽  
R. Cullen ◽  
N. Galin ◽  
...  

Abstract. We have derived Digital Elevation Models (DEMs) over the western part of the Devon Ice Cap in Nunavut Canada using "swath processing" of interferometric data collected by Cryosat 2 between February 2011 and January 2012. In the standard ESA "SARIn" level 2 (L2) data product the interferometric mode is used to map the cross-track position and elevation of the "point-of-closest-approach" (POCA) in sloping glacial terrain. However, in this work we explore the extent to which the phase of the returns in the intermediate L1b product can also be used to map the heights of time delayed footprints beyond the POCA. We show that there is a range of average cross-track slopes (~0.5° to ~2°) for which the returns will be dominated by those beneath the satellite in the main beam of the antenna and that the resulting interferometric phase allows mapping of heights in the delayed range window beyond the POCA. In this way a swath of elevation data is mapped allowing the creation of Digital Elevation Models (DEMs) from a sequence of L1b SARIn Cryosat-2 data takes. Comparison of the Devon results with airborne scanning laser data showed a mean difference of order a meter with a standard deviation < 1 m. The limitations of swath processing, which generates almost 2 orders of magnitude more data than traditional radar altimetry, are explored through simulation, and the strengths and weaknesses of the technique are discussed.


2017 ◽  
Vol 63 (239) ◽  
pp. 436-449 ◽  
Author(s):  
WESLEY VAN WYCHEN ◽  
JAMIE DAVIS ◽  
LUKE COPLAND ◽  
DAVID O. BURGESS ◽  
LAURENCE GRAY ◽  
...  

ABSTRACTFeature tracking of approximately annually separated Landsat-7 ETM+ imagery acquired from 1999 to 2010 and speckle tracking of 24-day separated RADARSAT-2 imagery acquired from 2009 to 2015 reveal that motion of the major tidewater glaciers of Devon Ice Cap is more variable than previously described. The flow of almost half (six of 14) of the outlet glaciers slowed over the observation period, while that of the terminus regions of three of 14 of the glaciers sped up in the most recent years of observation. The North Croker Bay Glacier of southern Devon Ice Cap showed the greatest variability in motion, oscillating between multi-year (three or more) periods of slower and faster flow and exhibited a pattern of velocity variability that is different from that of the rest of the ice cap's outlet glaciers. Comparisons between areas of dynamic variability and glacier bed topography indicate that velocity variability is largely restricted to regions where the glacier bed is grounded below sea level. Derived velocities are combined with measurements of ice thickness at the fronts of tidewater glacier to determine a mean annual (2009; 2011–15) dynamic ice discharge of 0.41 ± 0.11 Gt a−1for Devon Ice Cap. The Belcher Glacier is becoming a larger source of mass loss via ice discharge.


2014 ◽  
Vol 60 (223) ◽  
pp. 813-833 ◽  
Author(s):  
S. Gogineni ◽  
J.-B. Yan ◽  
J. Paden ◽  
C. Leuschen ◽  
J. Li ◽  
...  

AbstractThis paper presents the bed topography of Jakobshavn Isbræ, 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 Isbræ, 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 Isbræ. 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.


2016 ◽  
Author(s):  
Laurence Gray ◽  
David Burgess ◽  
Luke Copland ◽  
Thorben Dunse ◽  
Kirsty Langley ◽  
...  

Abstract. We compare geocoded heights derived from the interferometric mode (SARIn) of CryoSat to surface heights from calibration-validation sites on Devon Ice Cap and West Greenland. Comparisons are included for both the heights derived from the first return (the "point-of-closest-approach" or POCA) as well as heights derived from delayed waveform returns ("swath" processing). While swath processed heights are normally less accurate than edited POCA heights, of order 1–5 m instead of order 1–2 m, the increased coverage possible with swath data complements the POCA data and provides useful information for both system calibration and improving digital elevation models (DEMs). We show that the pre-launch interferometric baseline coupled with an additional roll correction (~ 0.0075°), or equivalent phase correction (~ 0.0435 radians), provides an improved calibration of the interferometric SARIn mode. We extend the potential use of SARIn data by showing the influence of surface conditions, especially melt, on the return waveforms, and that it is possible to detect and measure the height of summer supraglacial lakes in West Greenland. A supraglacial lake can provide a strong radar target in the waveform, stronger than the initial POCA return, if viewed at near normal incidence. This provides an ideal situation for swath processing and we demonstrate height accuracies of ~ 0.5 m for two lake sites, one in the accumulation zone and one in the ablation zone, which were measured every year from 2010 or 2011 to 2016. Each year the lake in the ablation zone was viewed in June by ascending passes and then 5.5 days later by descending passes which allows an approximate estimate of the filling rate. The results suggest that CryoSat waveform data and measurements of supraglacial lake height change could complement the use of optical satellite and be helpful as proxy indicators for surface melt around Greenland.


2005 ◽  
Vol 42 ◽  
pp. 83-89 ◽  
Author(s):  
Donghui Yi ◽  
H. Jay Zwally ◽  
Xiaoli Sun

AbstractThe Ice, Cloud and land Elevation Satellite (ICESat) in its 8 day repeat orbit mode provided data not only on the along-track surface slope, but also on the cross-track surface slope from adjacent repeat ground tracks. During the first 36 days of operation, four to five such repeat orbits occurred within 1 km in the cross-track direction. This provided an opportunity to use ICESat data to measure surface slope in the cross-track direction at 1 km scale. An algorithm was developed to calculate the cross-track surface slope. Combining the slopes in the cross-track and along-track directions gives a three-dimensional surface slope at 1 km scale. The along-track surface slope and surface roughness at 10km scale are also calculated. A comparison between ICESat surface elevation and a European Remote-sensing Satellite (ERS-1) 5 km digital elevation model shows a difference of 1–2 m in central Greenland where the surface slope is small, and >20m at the edge of Greenland where the surface slope is large. The large elevation difference at the edge is most likely due to the slope-induced error in radar altimeter measurement. Accurate surface slope data from ICESat will help to correct the slope-induced error of radar altimeter missions such as Geosat, ERS-1 and ERS-2.


2017 ◽  
Vol 11 (3) ◽  
pp. 1041-1058 ◽  
Author(s):  
Laurence Gray ◽  
David Burgess ◽  
Luke Copland ◽  
Thorben Dunse ◽  
Kirsty Langley ◽  
...  

Abstract. We compare geocoded heights derived from the interferometric mode (SARIn) of CryoSat to surface heights from calibration–validation sites on Devon Ice Cap and western Greenland. Comparisons are included for both the heights derived from the first return (the point-of-closest-approach or POCA) and heights derived from delayed waveform returns (swath processing). While swath-processed heights are normally less precise than edited POCA heights, e.g. standard deviations of  ∼  3 and  ∼  1.5 m respectively for the western Greenland site, the increased coverage possible with swath data complements the POCA data and provides useful information for both system calibration and improving digital elevation models (DEMs). We show that the pre-launch interferometric baseline coupled with an additional roll correction ( ∼  0.0075° ± 0.0025°), or equivalent phase correction ( ∼  0.0435 ± 0.0145 radians), provides an improved calibration of the interferometric SARIn mode. We extend the potential use of SARIn data by showing the influence of surface conditions, especially melt, on the return waveforms and that it is possible to detect and measure the height of summer supraglacial lakes in western Greenland. A supraglacial lake can provide a strong radar target in the waveform, stronger than the initial POCA return, if viewed at near-normal incidence. This provides an ideal situation for swath processing and we demonstrate a height precision of  ∼  0.5 m for two lake sites, one in the accumulation zone and one in the ablation zone, which were measured every year from 2010 or 2011 to 2016. Each year the lake in the ablation zone was viewed in June by ascending passes and then 5.5 days later by descending passes, which allows an approximate estimate of the filling rate. The results suggest that CryoSat waveform data and measurements of supraglacial lake height change could complement the use of optical satellite imagery and be helpful as proxy indicators for surface melt around Greenland.


2013 ◽  
Vol 59 (217) ◽  
pp. 874-882 ◽  
Author(s):  
Tyler Sylvestre ◽  
Luke Copland ◽  
Michael N. Demuth ◽  
Martin Sharp

AbstractGround-penetrating radar (GPR) surveys at a center frequency of 500 MHz were used to determine winter (2007/08) and net annual (2005–07) snow water equivalent (SWE) patterns across the upper parts of Belcher Glacier, Devon Ice Cap, Nunavut, Canada. The GPR measurements were validated against snow depths determined from avalanche probe measurements, and converted to SWE values using densities measured with a down-borehole neutron density probe and in shallow snow pits. Distinct internal reflection horizons (IRHs) in the GPR record were formed during warm summers in 2007 and 2005, and a large rain event in summer 2006 which caused ice to accumulate above the 2005 melt surface. Elevation provides the dominant control on winter SWE distribution across the basin, with surface topography (e.g. gullies) also being locally important. Based on the location where IRHs intersected the ice-cap surface, the basin-wide firn line occurred at an altitude of 1260–1300 m over the period 2005–08. Net mass balance across the accumulation area of Belcher Glacier averaged 0.24 m w.e. a−1 over the period 2005–07, mainly dependent on altitude. This is a little higher than most previous estimates for the period since the 1960s, although the differences lie within error limits.


Sensors ◽  
2018 ◽  
Vol 18 (7) ◽  
pp. 2265 ◽  
Author(s):  
Qingqing Feng ◽  
Huaping Xu ◽  
Zhefeng Wu ◽  
Wei Liu

Deceptive jamming against synthetic aperture radar (SAR) can create false targets or deceptive scenes in the image effectively. Based on the difference in interferometric phase between the target and deceptive jamming signals, a novel method for detecting deceptive jamming using cross-track interferometry is proposed, where the echoes with deceptive jamming are received by two SAR antennas simultaneously and the false targets are identified through SAR interferometry. Since the derived false phase is close to a constant in interferogram, it is extracted through phase filtering and frequency detection. Finally, the false targets in the SAR image are obtained according to the detected false part in the interferogram. The effectiveness of the proposed method is validated by simulation results based on the TanDEM-X system.


2015 ◽  
Vol 9 (3) ◽  
pp. 2821-2865 ◽  
Author(s):  
L. Gray ◽  
D. Burgess ◽  
L. Copland ◽  
M. N. Demuth ◽  
T. Dunse ◽  
...  

Abstract. We show that the CryoSat-2 radar altimeter can provide useful estimates of surface elevation change on a variety of Arctic ice caps, on both monthly and yearly time scales. Changing conditions, however, can lead to a varying bias between the elevation estimated from the radar altimeter and the physical surface due to changes in the contribution of subsurface to surface backscatter. Under melting conditions the radar returns are predominantly from the surface so that if surface melt is extensive across the ice cap estimates of summer elevation loss can be made with the frequent coverage provided by CryoSat-2. For example, the average summer elevation decreases on the Barnes Ice Cap, Baffin Island, Canada were 2.05 ± 0.36 m (2011), 2.55 ± 0.32 m (2012), 1.38 ± 0.40 m (2013) and 1.44 ± 0.37 m (2014), losses which were not balanced by the winter snow accumulation. As winter-to-winter conditions were similar, the net elevation losses were 1.0 ± 0.2 m (winter 2010/2011 to winter 2011/2012), 1.39 ± 0.2 m (2011/2012 to 2012/2013) and 0.36 ± 0.2 m (2012/2013 to 2013/2014); for a total surface elevation loss of 2.75 ± 0.2 m over this 3 year period. In contrast, the uncertainty in height change results from Devon Ice Cap, Canada, and Austfonna, Svalbard, can be up to twice as large because of the presence of firn and the possibility of a varying bias between the true surface and the detected elevation due to changing year-to-year conditions. Nevertheless, the surface elevation change estimates from CryoSat for both ice caps are consistent with field and meteorological measurements. For example, the average 3 year elevation difference for footprints within 100 m of a repeated surface GPS track on Austfonna differed from the GPS change by 0.18 m.


Sign in / Sign up

Export Citation Format

Share Document