-
The Cryosphere, 12, 2969–2979,
2018https://doi.org/10.5194/tc-12-2969-2018© Author(s) 2018. This
work is distributed underthe Creative Commons Attribution 4.0
License.
Coherent large beamwidth processing of radio-echo sounding
dataAnton Heister and Rolf ScheiberGerman Aerospace Center (DLR),
Microwaves and Radar Institute, Wessling, Germany
Correspondence: Anton Heister ([email protected])
Received: 23 March 2018 – Discussion started: 18 April
2018Revised: 19 August 2018 – Accepted: 22 August 2018 – Published:
19 September 2018
Abstract. Coherent processing of radio-echo sounding dataof
polar ice sheets is known to provide an indication ofbedrock
properties and detection of internal layers. We in-vestigate the
benefits of coherent processing of a large az-imuth beamwidth to
retrieve and characterize the orienta-tion and angular
backscattering properties of the surface andsubsurface features.
MCRDS data acquired over two distincttest areas in Greenland are
used to demonstrate the specu-lar backscattering properties of the
ice surface and of the in-ternal layers, as well as the much wider
angular response ofthe bedrock. The coupling of internal layers’
orientation withthe bed topography is shown to increase with depth.
Spec-tral filtering can be used to increase the SNR of the
internallayers and mitigate the surface multiple. The variance of
thebed backscattering can be used to characterize the bed
returnspecularity. The use of the SAR-focused RES data ensuresthe
correct azimuth positioning of the internal layers for
thesubsequent slope estimation.
1 Introduction
Radio-echo sounding (RES) is a well-established techniquefor
remotely measuring the thickness of ice sheets. The use ofsynthetic
aperture radar (SAR) focusing improves gain andazimuth resolution
of the echograms. Overall, state-of-the-artSAR processing offers
information about the spatial proper-ties of the ice sheet and the
strength of the response, which isused to determine ice thickness,
internal layers’ orientationand bedrock conditions, i.e., presence
or absence of water.There exist several SAR algorithms for focusing
RES data,among them 1-D matched filtering (Legarsky et al.,
2001),the ω−k migration (Leuschen et al., 2000), 2-D matched
fil-tering (Heliere et al., 2007; Peters et al., 2007), and
multilooktime-domain back-projection (Mishra et al., 2016).
Addition-
ally, Holschuh et al. (2014) offer a method for improvingSAR
focusing of internal layers by introducing a correctionof
attenuation, migration and radial spreading for the returnsfrom
tilted internal layers.
Previous studies of angular backscattering properties ofthe ice
sheet and bed were performed in Jezek et al. (2009);Smith (2014);
Schroeder et al. (2013); MacGregor et al.(2015). Jezek et al.
(2009) offer a technique for studyingthe backscattering properties
of the ice sheet and bed usinga special subaperture SAR approach.
The authors study thedependency of the surface and bed return power
on the in-cidence angle and the effect of the surface slope on the
sur-face return power. They show that the response of the inter-nal
layers is specular and propose incoherent summation ofsubapertures
to improve the signal-to-noise ratio (SNR) ofinternal layers. Smith
(2014) estimates an optimal value forthe SAR beamwidth based on the
bedrock SNR. Schroederet al. (2013) offer an approach for detecting
the presence ofsubglacial water at the bed based on its angular
backscat-tering characteristics. The authors estimate the
specularityof the bed returns by comparing power contributions in
twoHiCARS 60 MHz (Peters et al., 2007) SAR echograms withsynthetic
apertures of 700 and 2000 m. MacGregor et al.(2015) introduce two
new methods for estimating the slopeof internal layers, among them
the Doppler centroid method,which leverages the fact that internal
layers’ returns arehighly specular. The authors use azimuth Fourier
transformof short overlapping range-compressed RES data blocks,
andderive the slope of internal layers from the wavenumber ofthe
corresponding Doppler centroids.
Mishra et al. (2016) introduce a novel approach for
SARprocessing of RES data, where the processing chain gener-ates a
number of SAR echograms, each corresponding toa particular
incidence angle in along-track. A subset of the
Published by Copernicus Publications on behalf of the European
Geosciences Union.
-
2970 A. Heister and R. Scheiber: Processing of radio-echo
sounding data
echograms with the highest SNR is then selected for
furtherprocessing.
In this paper we introduce a new flexible technique to ana-lyze
the angular backscattering properties of the ice sheet andbed,
which can be applied to previously conventionally SAR-focused
complex-valued echograms. A better understandingof those properties
allows us to offer novel strategies for im-proving internal layer
and possibly bed SNR, to mitigate thesurface multiple return and to
train sparsifying dictionariesfor model-based cross-track focusing
methods such as Wuet al. (2011); Heister and Scheiber (2016).
This paper begins with a description of the employedSAR-focusing
algorithm for RES data in Sect. 2. After thatwe introduce the
technique for analyzing angular backscat-tering properties of the
ice sheet and bed in Sect. 3. In Sect. 4we analyze the processing
results for two RES surveyscollected by the Center for Remote
Control of Ice Sheets(CReSIS), Kansas, USA using their
Multi-Channel RadarDepth Sounder (MCRDS) (Lohoefener, 2006;
Marathe,2008) during the Greenland campaign in 2008 (CReSIS,2012).
Based on the results of Sect. 4, we discuss and demon-strate
approaches for improving internal-layer visibility andfor
mitigating the surface multiple in Sect. 5. Potential im-pacts for
the scientific evaluation of SAR-focused RES datawith large
beamwidth are discussed in Sect. 6. Finally, a sum-mary and
conclusions are given.
2 SAR focusing
We perform SAR focusing of RES data using a modificationof the
range-Doppler algorithm. The processing is done inoverlapping 8000
m long azimuth blocks, with each blockprocessed as described in
Algorithm 1. For each block weassume the platform to fly with a
constant velocity v, the icesurface to have a constant along-track
slope ψ , and the icesheet to have a constant refractive index nice
= 1.78 with anequivalent real part of relative permittivity εice =
3.17. Wealso assume that the electromagnetic wave propagation
obeysSnell’s law for a two-layer air-ice model shown in Fig. 1.The
number of azimuth samples in each block is selectedto satisfy at
least twice the desired SAR beamwidth of 1θ =30◦. We additionally
assume that the azimuth antenna patternis broad enough so that its
variation for incidence angles inthe interval θ = [−15◦,15◦] can be
safely ignored.
We now describe the algorithm inputs using the notationwhere τ
denotes range time, fτ denotes range frequency,η denotes azimuth
time, and fη denotes azimuth frequency(Algorithm 1).
Range compression, which is a signal-processing tech-nique for
improving the radar range resolution, is imple-mented using a
matched filter HRC equal to a complex con-jugate of the Fourier
transform of the transmitted signalweighted by Hamming or Blackman
windows for the side-lobe suppression.
Ricenice=1.78
nair=1
Rair
Point target
Radar
d
s
R0
Azimuthx
ψ
θ Range
Figure 1. Along-track geometry.
The range-Doppler algorithm assumes a linear motion tra-jectory
of the platform; therefore the motion compensation,a procedure that
corrects the platform’s trajectory deviationfrom a linear reference
trajectory, is needed. We implement itusing a filter HMOCO, which
only corrects for a vertical com-ponent of the platform’s deviation
from a reference track inthe range frequency domain.
Algorithm 1 SAR focusing
Require: raw data DATA, filters HRC, HMOCO, and HREF,amount of
RCM 1RRCM.
Ensure: SAR-focused echogram DATASAR1: DATA := FFTrange(DATA)2:
DATA := DATA ·HRC ·HMOCO3: DATA := IFFTrange(DATA)4: DATA :=
FFTazimuth(DATA)5: for fη ∈ [−Baz/2,Baz/2] do6: DATA[:,fη] :=
interp(DATA[:,fη],1RRCM[:,fη])7: end for8: DATA := DATA ·HREF9:
DATASAR := IFFTazimuth(DATA)
10: return DATASAR
As the platform moves in azimuth, the response from atarget
spreads across multiple range bins. Range cell migra-tion
correction (RCMC) is an operation that removes thisrange variation,
bringing the target response to a fixed rangebin at every azimuth
position. RCMC is performed in therange-Doppler domain. During RCMC
every range line isshifted by the time corresponding to the amount
of rangecell migration 1RRCM. The spatially variant shift in
rangeis implemented using a sinc interpolator with Lanczos win-dow
(a = 2) and the length exceeding 3 times the maximal1RRCM. Finally,
azimuth compression is done by applyingthe HREF filter. We now
derive equations for 1RRCMC andHREF.
The Cryosphere, 12, 2969–2979, 2018
www.the-cryosphere.net/12/2969/2018/
-
A. Heister and R. Scheiber: Processing of radio-echo sounding
data 2971
From Fig. 1 the optical path length R from the radar atazimuth x
= η · v to a point target at depth d is
R(d,η)= Rair+ niceRice, (1)
where the geometric lengths that the electromagnetic wavetravels
in air and ice are
Rair =
√(R20 + s · tanψ)
2+ (x− s)2, (2)
Rice =
√s2+ (d − s · tanψ)2. (3)
Both Eqs. (2) and (3) depend on an unknown location ofthe
refraction point s, which is a function of time η and depthd . The
refraction point s can be found by solving a fourth-order
polynomial equation (Heliere et al., 2007; Scheiberet al., 2008)
or, more efficiently, by using Newton’s optimiza-tion method, which
iteratively finds s that minimizes (Eq. 1)with the following update
rule at (i+ 1)th iteration:
si+1 = si −R′(s)/R′′(s), (4)
where we initialize the refraction point with s0 = 0.Knowing s,
we calculate the phase shift of the received
signal with respect to the time η = 0, when the platformcrosses
the origin of the x axis
1φ(η)=4πλ01R(η)=
4πλ0(R(η)−R(0)), (5)
where λ0 is the wavelength of the transmitted wave in the
air.The Doppler frequency shift of the received signal is pro-
portional to the derivative of Eq. (5) in time:
1f (η)=1
2π∂φ(η)
∂η. (6)
Knowing Eq. (6) for each depth d and azimuth posi-tion x, we
compute the amount of range cell migrationin the range-Doppler
domain 1RRCM(d,fη) by interpolat-ing its time domain counterpart
1RRCM(d,η)= τ · c0/2−1R(d,η) onto a regularly sampled azimuth
frequency gridfη ∈ [−Baz/2,Baz/2], where fη =±Baz/2 corresponds
toincidence angles θ =±15◦.
Finally, we compute 1φ(fη) by interpolating Eq. (5) ontofη and
calculate the matched filter for SAR-focusing HREFas
HREF(τ,fη)= exp(−j1φ(fη)). (7)
We note that more precise and less restrictive SAR-focusing
algorithms for ice-sounder data exist, such as time-domain
back-projection (Mishra et al., 2016); our choice of aparticular
approach described above is based on simplicity ofimplementation
and its sufficiency for the subsequent analy-sis of the ice sheet
and bed angular backscattering properties.
3 Multiple subbands processing
In order to analyze the dependency of the
backscatteringproperties of the ice sheet and bed on the incidence
angle,we divide the azimuth spectrum of an echogram intoN
over-lapping subbands of beamwidth 1θsub = 2◦ and an overlapbetween
two adjacent subbands of 1◦, with each subbandweighted by a
rectangular window. The central frequency ofthe n ∈ (1,N) subband,
f0(θn), corresponds to the incidenceangle of interest θn ∈
[−14◦,14◦]:
f0(θn)=2v sinθnλ0
. (8)
Each subband is then accordingly zero-padded in azimuthso that
all N subbands have the same size. An inverse az-imuth Fourier
transform is then applied to each subband toget a set of N
echograms In, each containing returns comingpredominantly from the
corresponding incidence angle θn.
The positions of ice-sheet features of interest, such as
sur-face, internal layers and bed, are then manually selected
froman echogram Iincoh, calculated as the incoherent sum
Iincoh =
N∑n=1|In|. (9)
We note that the ice sheet features can be
trackedsemi-automatically or automatically; however, for the
smallamount of data we analyze in the paper, manual selection
isfeasible.
4 Greenland MCRDS data
We apply the approach presented in Sect. 3 to RES data
col-lected by CReSIS (CReSIS, 2012) using their MCRDS sys-tem
(Lohoefener, 2006; Marathe, 2008). The main parame-ters of the
radar and the acquisitions are summarized in Ta-ble 1. Two chirps
with different durations were transmittedalternately on a
pulse-to-pulse basis, with a 3 µs chirp in-tended to capture the
surface and the shallow internal-layerreturns (shallow mode), and a
10 µs chirp intended to cap-ture deeper internal layers and bed
returns (deep mode). Weemploy the availability of multiple
cross-track channels ofMCRDS to increase the SNR of nadir returns
by combiningthe SAR echograms of the cross-track channels using
con-ventional delay and sum beamforming.
We select two tracks, both flown over Greenland in sum-mer 2008
and approximately 70 km long. The track flownfrom the inland
towards Jakobshavn Glacier is referred toas track 1, the track
flown over south-eastern Greenlandin a north-easterly direction is
referred to as track 2. Theregions of interest, their topography
and flight trajectoriesare shown in Fig. 2. These particular
surveys are chosento demonstrate how different bed topography
affects the re-flective properties of the internal layers; the bed
in track 1
www.the-cryosphere.net/12/2969/2018/ The Cryosphere, 12,
2969–2979, 2018
-
2972 A. Heister and R. Scheiber: Processing of radio-echo
sounding data
Table 1. Parameters of MCRDS acquisitions.
Parameters Track 1 Track 2
Central frequency 150 MHzChirp bandwidth 20 MHzChirp duration
3/10 µsSampling frequency 120 MHz
Effective PRF 78 Hz 156 HzNumber of cross-track channels 16
6Effective cross-track aperture 14.34 m 4.79 mAcquisition date 20
June 2008 1 August 2008Acquisition start UTC 18:32:29
16:49:49Acquisition end UTC 18:47:33 17:07:23Average height over
surface 160 m 800 mAverage velocity 78 m s−1 65 ms−1
Figure 2. MCRDS surveys on a map. The map of Greenland isplotted
using a stereographic projection with a central meridian of41◦W and
a central parallel of 72◦ N. Isolines on maps correspondto a
surface elevation change of 250 m.
has depth varying in the interval dbed ∈ [2170m,3030m] andslopes
varying in the interval ψbed ∈ [−35◦,33◦]. The corre-sponding
intervals for track 2 are dbed ∈ [640m,1970m] andψbed ∈ [−62◦,65◦].
Here we calculate the slopes of the bedand the internal layers
as
ψbed/layer(x)= tan−1(∂dbed/layer(x)
∂x· nice
), (10)
where nice scales the geometric slope to correspond to
theincidence angle observed by the radar.
The full bandwidth echogram of track 1 is shown inFig. 3a. To
produce the figure we combine echograms ofthe shallow and deep
modes, average the intensities of eacheight adjacent azimuth
samples, downsample the result in az-imuth by a factor of 8, and
add a depth-dependent amplitude
Figure 3. Backscattering characteristics of the ice sheet and
bed fortrack 1.
The Cryosphere, 12, 2969–2979, 2018
www.the-cryosphere.net/12/2969/2018/
-
A. Heister and R. Scheiber: Processing of radio-echo sounding
data 2973
Figure 4. Backscattering characteristics of the ice sheet and
bed fortrack 2.
ramp of 20 dB km−1 to improve the visibility of the deep
in-ternal layers and the bed. The internal layers are visible tod ≈
2 km. The gaps in internal-layer visibility occur at az-imuth
positions where the bed slope is the steepest.
First, we investigate reflective properties of the ice
surface.Figure 3b shows the normalized reflectivity power of the
sur-face as a function of incidence angle. The surface response
isspecular, with the incidence angle corresponding to the max-imum
intensity θmax(I ), varying slowly in azimuth.
To study backscattering properties of internal layers weselect a
single internal layer tracked with a solid red linein Fig. 3a. A
deep layer is selected in order to avoid unde-sired contributions
of the off-nadir surface returns. Figure 3cshows the internal
layer’s normalized power together with itsslope, computed from Eq.
(10) and drawn as a white line. Weuse bicubic interpolation to plot
the figure. The internal-layerresponse is specular, with θmax(I )
proportional to the layer’sgeometric slope.
A further insight into the behavior of the internal
layer’sresponse is given in Fig. 3d, where for each pixel of Iincoh
wecolor coded the incidence angle corresponding to the maxi-mum
intensity θmax(I ). Prior to plotting, we additionally ap-plied a
median filter of size (5,5) and bicubic interpolation.The black
lines on the figure correspond to the surface andbed return
positions. The figure shows a correlation betweenθmax(I ) and the
bed slope, with the blue and the red color ap-pearing at azimuth
positions with negative and positive bedslope correspondingly.
Moreover, for a given azimuth posi-tion x0 the absolute value of
θmax(I ) increases with depth;therefore, according to Fig. 3d, the
absolute value of inter-nal layers’ slope also increases with
depth. This implies thatthe deeper the internal layer is located,
the more its shaperesembles the shape of the bed.
Figure 3e shows the normalized power of the bed response,where,
prior to the normalization, we additionally compen-sate for the
two-way propagation power loss of 20 dB km−1.The incidence angle
θmax(I ) of the bed response varies in az-imuth, and overall the
response is wide, meaning the bed is arough surface for a radar
with λ0 = 2 m.
Figure 5 shows the dependency of the return power of
thepreviously selected internal layer and bed at four fixed
az-imuth positions. Those particular positions are selected
todemonstrate the variety of shapes of reflective signatures forthe
bed and the persistent signature shape for the internallayer.
The full bandwidth echogram for track 2 is shown inFig. 4a,
where we add a depth-dependent amplitude ramp of20 dB km−1. Here
the bed topography varies more comparedto track 1, the internal
layers are visible close to the bed withgaps appearing at azimuth
positions where the absolute valueof bed slope is the highest, and
the surface multiple is alsopresent in the echogram.
Figure 4b shows the normalized reflectivity power of thesurface.
The surface response is similar to the one for track 1,
www.the-cryosphere.net/12/2969/2018/ The Cryosphere, 12,
2969–2979, 2018
-
2974 A. Heister and R. Scheiber: Processing of radio-echo
sounding data
Figure 5. Internal layer (a) and bed (b) returns for track 1 at
az-imuth positions x = (3,9,46,59) km, marked with vertical red
linesin Fig. 3a. Quadratic interpolation was applied to smooth the
signa-tures.
with higher variation of θmax(I ) starting from azimuth x =65
km.
Reflective properties of a single internal layer tracked witha
solid red line in Fig. 4a are shown in Fig. 4c. Here we selecta
shallow layer because θmax(I ) for deeper layers would lieoutside
the interval θn ∈ [−14◦,14◦] previously selected inSect. 3. The
incidence angle θmax(I ) in Fig. 4c varies morestrongly and
frequently compared to that in Fig. 3c.
Figure 4d is plotted similarly to Fig. 3d. As expected,
weobserve larger color gradients for internal layers for track
2,whereas incidence angles of the surface multiple lie aroundθn ≈
0◦ in white, corresponding to the ice surface.
The normalized power of the bed response for track 2 isshown in
Fig. 4e.
In Fig. 5 we compare the responses of the previously se-lected
internal layer and bed at four fixed azimuth positions.We again
observe specular reflections from the internal layerand wider
reflections from the bed, with θmax(I ) for the bedand the internal
layer positively correlated for each selectedposition.
For both tracks the mean value of the surface and the
layerbeamwidth at the −6 dB level is 2.2◦. The shape of the
bedresponse varies and does not necessarily have a prominentsingle
peak; therefore we do not calculate its beamwidth.We suggest using
the variance of the bed angular responseto quantify its spread when
the bed roughness or the pres-ence of basal water is of interest.
For more details we referthe readers to Sect. 6.
Figure 6. Internal layer (a) and bed (b) returns for track 2 at
azimuthpositions x = (3,25,32,54) km, marked with vertical red
lines inFig. 4a. Quadratic interpolation was applied to smooth the
signa-tures.
5 Enhancement of SAR echograms
In this section we offer two straightforward applications ofthe
results provided in Sect. 4 to improve the SAR-focusedRES data.
First, the fact that an internal-layer response is narrowmeans
that for a given depth and azimuth it contributes onlyto a limited
azimuth frequency range of a SAR echogramspectrum. For the azimuth
spectrum of a small block of aSAR echogram we see that at each
depth internal-layer con-tribution is clustered around the
frequency corresponding tothe layer slope, which is demonstrated in
Fig. 7. This allowsus to use spectral filtering to improve the SNR
of the in-ternal layers. We perform the filtering using 250 m
azimuthblocks with a 70 % overlap. We apply a Fourier transformin
azimuth to each block, after which for each depth we se-lect a
frequency corresponding to the maximum spectrum in-tensity fmax(I
)(d) and fit it to a piecewise linear regressionwith npieces = 3 to
make the estimate f̂layer(d) of the internallayers’ frequency
flayer(d) more robust against outliers. De-pending on the shape of
flayer(d), other values for npieces aswell as other types of
regression (e.g., polynomial regression)might be used to calculate
f̂layer(d). After that we nullify thepart of the spectrum at
frequencies lying outside the intervalf̂layer(d)± 0.05Baz, whereas
the part of the spectrum lyinginside the interval is kept intact.
Finally we apply an inverseFourier transform in azimuth to each
block and re-assemblethe overlapping blocks. We apply this method
to the track 1SAR echogram. The results are shown in Fig. 8.
The procedure results in a 21.8 % sharpness improvementin terms
of intensity squared metric (Fienup and Miller,2003) (Eq. 11) with
the mean intensity of the echograms nor-
The Cryosphere, 12, 2969–2979, 2018
www.the-cryosphere.net/12/2969/2018/
-
A. Heister and R. Scheiber: Processing of radio-echo sounding
data 2975
Figure 7. Azimuth spectrum of a 250 m long SAR block in track1,
azimuth position x = 59 km. The frequency interval
f̂layer(d)±0.05Baz is marked with white lines.
Figure 8. Improvement of internal-layer visibility for track 1.
Thesubsets of the SAR echogram before and after the processing
areshown at the top and the bottom correspondingly.
malized prior to the comparison.
s(I )=∑i,j
Ii,j · I∗
i,j (11)
Figure 9 depicts the antenna power pattern, the layer
powerspectrum centered around flayer, and the power spectrumof the
noise. The noise power is proportional to the inte-gral of its
power spectrum. Prior to the spectral filtering,this frequency
range contains the entire processed bandwidthBaz, whereas after the
filtering the interval is fη ∈ [flayer−0.05Baz,flayer+ 0.05Baz].
Assuming white Gaussian noise,and that the signal energy is not
affected by the filtering, theSNR improvement is 10 dB.
P
f layerPsignal+Pnoise
Psignal+P'noise0.1 Baz
fη
Baz
Antenna pattern
Signal
Noise
Figure 9. Noise reduction for internal-layer SNR
improvement.
Second, according to the Fig. 4d, the contribution of thesurface
multiple return, which is a multipath reflection fromthe ice
surface as well as from the upper internal layers andthe bottom of
the aircraft fuselage and wings, can be miti-gated by identifying
and filtering out its contribution in theazimuth frequency domain,
therefore revealing previouslymasked internal layers. This
approach, however, only worksin areas where the θmax(I ) for
internal layers and the surfacemultiple differ.
To demonstrate the approach we process a subset of a SARechogram
of track 2. The processing is done in blocks withsimilar parameters
as for internal-layer visibility improve-ment. In each block we
locate the strongest surface multiplecontribution at a depth dmult
corresponding to the doubledheight of the aircraft over the ice
surface, and at azimuth fre-quency fmult which corresponds to the
frequency of strongestsurface return. After that in each block we
apply a 2-D notchfilter located at depth d ∈ [dmult,dmult+ 100m]
and at fre-quency fη ∈ [fmult− 0.05Baz,fmult+ 0.05Baz]. The
resultsare shown in Fig. 10.
6 Scientific utility of large beamwidth SAR processing
In this section we discuss the practical benefit of
largebeamwidth SAR processing of RES data for two
scientificapplications, namely basal water detection and estimation
ofinternal-layer slope.
Presence of water bodies can be detected in RES data us-ing
either amplitude or angular information. The amplitudedetection is
based on the fact that the basal water producesstronger reflection
compared to the grounded bed. Howeverthis method is prone to errors
(Matsuoka, 2011) as the radarattenuation depends on the chemical
composition of the iceas well as its temperature (MacGregor et al.,
2007). The anal-ysis of the angular backscattering distribution of
the bed re-turn, on the other hand, is free of the aforementioned
limita-
www.the-cryosphere.net/12/2969/2018/ The Cryosphere, 12,
2969–2979, 2018
-
2976 A. Heister and R. Scheiber: Processing of radio-echo
sounding data
Figure 10. Mitigation of the surface multiple for track 2. The
sub-sets of the SAR echogram before and after the mitigation are
shownat the top and the bottom correspondingly.
tions, and a specular bed response indicates the presence
ofbasal water.
In order to quantify specularity of the bed return,Schroeder et
al. (2013) introduce the specularity content,a measure which is
calculated as a ratio of energies ofthe bed return in two SAR
echograms, I1 and I2, focusedwith synthetic apertures L1 = 700 m
and L2 = 2000 m cor-respondingly. The echogram I1 contains specular
returns,whereas the echogram I2 contains both specular and
non-specular returns. For a typical height above the ice surfaceR0
= 500 m and ice thickness d = 2km used in the survey(Schroeder et
al., 2014), L1 and L2 respectively correspondto beamwidths 1θ1 ≈
10◦ and 1θ2 ≈ 28◦.
We compare the specularity content used so far with thevariance
of angular backscattering of the bed, which is ameasure of the bed
angular distribution spread, introduced inSect. 4. We calculate the
specularity content using1θ1 = 10◦
and 1θ2 = 30◦. To calculate the variance we normalize theenergy
of the angular backscattering in each azimuth posi-tion. In Fig. 11
the specularity content is low at azimuth po-sitions A and D,
failing to detect specular bed reflections lo-cated outside of 1θ1;
the specularity content is high at az-imuth positions B and C, but
decays rapidly as the bed re-flection moves outside of 1θ1. The
variance, on the otherhand, is insensitive to an angular shift of
the reflected energy,therefore making it possible to detect
additional specular re-flections. Using the variance instead of the
specularity con-tent can potentially lead to better detection and
mapping ofsubglacial water bodies, especially in areas where the
bed istilted in azimuth.
Internal layers, frequently observed in RES echograms,are widely
attributed to the changes of electrical conductiv-ity within the
ice sheet. Owing to the fact that the internallayers are considered
isochrones (Hempel et al., 2000), whentracked in RES data, they
provide information about changes
Figure 11. Specularity content and variance as measures of
spec-ularity of the bed return for track 2. For comparison we take
twoparts of track 2 each 2.3 km long. The power of the angular
energyreturns (a) and (c) is normalized with respect to the highest
value ofthe distribution at each azimuth position. The area inside
the dashedhorizontal lines in (a) and (c) corresponds to 1θ1.
in the ice flow (Siegert, 2004) and snow accumulation rate inthe
past (Fahnestock et al., 2001). The availability of a largeamount
of RES data and the fact that the semi-automaticlayer tracking is
prohibitively expensive (Sime et al., 2011)motivate the development
of automatic layer-tracking algo-rithms.
Some of the tracking algorithms use reflection slope topredict
the internal layering, which in turn simplifies sub-sequent
tracing. MacGregor et al. (2015) introduce two newmethods for the
slope estimation, namely horizontal phase
The Cryosphere, 12, 2969–2979, 2018
www.the-cryosphere.net/12/2969/2018/
-
A. Heister and R. Scheiber: Processing of radio-echo sounding
data 2977
RCSAR
(a)
(b)
RCSAR
Δx
Depth
Azimuth
Figure 12. Azimuth displacement of a tilted internal layer in
range-compressed data.
gradient and Doppler centroid methods. Both methods usecoherent
RES data (phase preserved). As stated in MacGre-gor et al. (2015)
currently the data have been range com-pressed but without SAR
focusing. The use of range com-pressed as opposed to SAR-focused
RES data might lead toerroneous estimation of the slope due to the
displacementof internal layers in azimuth. The displacement is due
to thefact that the internal-layer reflection is specular, which in
turnmeans that prior to SAR focusing the return from an
internallayer appears at the azimuth position, where the incident
en-ergy is normal to the layer’s surface. SAR focusing registersthe
return at its zero Doppler position, which corresponds tothe nadir
direction.
Figure 12 schematically illustrates this effect with two
ex-amples. For the sake of simplicity we ignore the ray bendingdue
to the difference in refractive properties of the air and ice.The
azimuth positions at which the layer returns are regis-tered in
range compressed and SAR-focused data are markedwith circles and
squares, respectively. A convex internal layerin Fig. 12a appears
stretched in azimuth in range-compresseddata, whereas a concave
internal layer in Fig. 12b appearsshrunk (in extreme cases
reflections from the left and theright would be overlaid).
The amount of the displacement1x depends on the layer’sgeometric
slope α, its depth d, platform height above the sur-face R0, the
refractive index of the ice nice, and the surfaceslope ψ . When ψ =
0◦, the displacement is calculated as
1x = R0 tan(sin−1(nice sinα))+ niced tanα. (12)
Figure 13. Internal-layer displacement in azimuth.
Figure 14. Stretching and shrinkage of internal layers in
range-compressed data of track 2. The range-compressed echogram
isshown in (a), the SAR-focused echogram is shown in (b).
Figure 13 shows the dependency of 1x on the layer ge-ometric
slope α and depth d when the height over the sur-face is R0 = 800
m, the ice refractive index is nice = 1.78,and ψ = 0◦.
Real data examples for both cases are demonstrated inFig. 14,
where we compare range-compressed and SAR-focused subsets of RES
data for track 2. For the range-compressed data, the internal
layers appear stretched in az-imuth in areas A and C, whereas for
the area B the layersshrink and even overlay at the lower
depth.
7 Summary and conclusions
In this paper we offered a new approach to study
scatteringcharacteristics of ice sheets, which is based on the
division ofa conventionally focused large beamwidth ice sounder
SARechogram into a set of subbands, each of which correspondto a
particular incidence angle in an along-track direction.We estimated
and compared scattering characteristics of theice surface, internal
layers and bed for two surveys in Green-land. For those surveys,
the surface and internal layers have
www.the-cryosphere.net/12/2969/2018/ The Cryosphere, 12,
2969–2979, 2018
-
2978 A. Heister and R. Scheiber: Processing of radio-echo
sounding data
narrow responses, which correspond to a smooth specularsurface,
while the bed response is wide, which correspondsto a rough
surface. The scattering properties carry informa-tion which can be
used to estimate characteristics of the bedroughness (Fung and Eom,
1983), with the specular bed re-sponse indicating the likely
presence of subglacial water atthe bed (Schroeder et al.,
2013).
Based on the scattering characteristics of internal layers,we
offered a post-processing technique to improve their vis-ibility.
By taking a small azimuth block of a SAR echogram,within which the
orientation of internal layers varies slightlyin along-track, we
observe that the internal layer’s contribu-tion to the block’s
azimuth spectrum is sparse and is clusteredaround the frequency
corresponding to the internal layer’sslope. This observation
directly suggests a way to improveinternal layers’ SNR by keeping
only those spectral com-ponents where the internal layers’
contributions are present.This post-processing technique can
improve spatial trackingand interpretation of both horizontal and
tilted internal lay-ers. As a subject for further studies, we
suggest that denois-ing all ice sheet features in a SAR echogram is
possible byfinding a sparse representation of the echogram given a
spar-sifying dictionary learned on patches with high SNR.
We also demonstrated a way to reduce the undesired con-tribution
of the surface multiple return, which masks inter-nal layers at
corresponding depths. The reduction is possiblewhenever the surface
multiple and the masked layer contri-butions come from different
incidence angles, in which casethey are separable in the azimuth
frequency domain.
Finally, we discussed the potential benefit offered by
theanalysis of RES data focused using large beamwidth SARwith
regards to the bed specularity characterization and forthe correct
azimuth positioning of the tilted internal layers.
Data availability. The raw MCRDS data used in the paper
areavailable from CReSIS upon request.
Competing interests. The authors declare that they have no
conflictof interest.
Acknowledgements. We acknowledge the use of data fromCReSIS
generated with support from the University of Kansas,NASA Operation
IceBridge grant NNX16AH54G and NSF grantACI-1443054. The authors
would also like to thank John Paden ofCReSIS for answering the
sensor and data-related questions.
The article processing charges for this open-accesspublication
were covered by a ResearchCentre of the Helmholtz Association.
Edited by: Joseph MacGregorReviewed by: John Paden and one
anonymous referee
References
CReSIS: Radar Depth Sounder Data, Lawrence, Kansas, USA,
Dig-ital Media, available at: https://data.cresis.ku.edu/ (last
access:1 June 2018), 2012.
Fahnestock, M., Abdalati, W., Luo, S., and Gogineni, S.:
Inter-nal layer tracing and age-depth-accumulation relationships
forthe northern Greenland ice sheet, J. Geophys. Res.-Atmos.,
106,33789–33797, https://doi.org/10.1029/2001jd900200, 2001.
Fienup, J. R. and Miller, J. J.: Aberration correction by
maximizinggeneralized sharpness metrics, J. Opt. Soc. Am., 20,
609–620,https://doi.org/10.1364/JOSAA.20.000609, 2003.
Fung, A. and Eom, H.: Coherent scattering of a spherical
wavefrom an irregular surface, IEEE T. Antenn. Propag., 31,
68–72,https://doi.org/10.1109/tap.1983.1142979, 1983.
Heister, A. and Scheiber, R.: First Analysis of Sparse Signal
Recon-struction for Radar Imaging of Ice Sheets, in: Proceedings of
the11th European Conference on Synthetic Aperture Radar,
VDE,Berlin, Germany, 6–9 June 2016, 788–791, 2016.
Heliere, F., Lin, C.-C., Corr, H., and Vaughan, D.: Ra-dio echo
sounding of Pine Island Glacier, West Antarc-tica: aperture
synthesis processing and analysis of feasibil-ity from space, IEEE
T. Geosci. Remote, 45,
2573–2582,https://doi.org/10.1109/tgrs.2007.897433, 2007.
Hempel, L., Thyssen, F., Gundestrup, N., Clausen, H. B.,
andMiller, H.: A comparison of radio-echo sounding data and
elec-trical conductivity of the GRIP ice core, J. Glaciol., 46,
369–374,https://doi.org/10.3189/172756500781833070, 2000.
Holschuh, N., Christianson, K., and Anandakrishnan, S.:Power
loss in dipping internal reflectors, imaged us-ing ice-penetrating
radar, Ann. Glaciol., 55,
49–56,https://doi.org/10.3189/2014AoG67A005, 2014.
Jezek, K., Gogineni, S., Rodriguez, E., Wu, X., Sonntag,
J.,Rodriguez, F., Freeman, A., and Curlander, J.: GlobalIce Sheet
Mapping observatory: final report, Tech. rep.,Earth Science
Technology Office, NASA, Greenbelt, MD,USA, available at:
https://pdfs.semanticscholar.org/d361/3bc68c51ebcffbaa0b9aa6056c1cf8108adc.pdf
(last access:20 March 2018), 2009.
Legarsky, J., Gogineni, S., and Akins, T.: Focused synthetic
aper-ture radar processing of ice-sounder data collected over
theGreenland ice sheet, IEEE T. Geosci. Remote, 39,
2109–2117,https://doi.org/10.1109/36.957274, 2001.
Leuschen, C., Gogineni, S., and Tammana, D.: SAR processing
ofradar echo sounder data, in: IEEE 2000 International
Geoscienceand Remote Sensing Symposium, Honolulu, HI, USA,
24–28July 2000, https://doi.org/10.1109/igarss.2000.859643,
2000.
Lohoefener, A.: Design and Development of a Multi Channel
RadarDepth Sounder, Master’s thesis, University of Kansas,
2006.
MacGregor, J. A., Winebrenner, D. P., Conway, H., Matsuoka,
K.,Mayewski, P. A., and Clow, G. D.: Modeling englacial
radarattenuation at Siple Dome, West Antarctica, using ice
chem-istry and temperature data, J. Geophys. Res., 112,
F03008,https://doi.org/10.1029/2006jf000717, 2007.
MacGregor, J. A., Fahnestock, M. A., Catania, G. A., Paden, J.
D.,Prasad Gogineni, S., Young, S. K., Rybarski, S. C., Mabrey,A.
N., Wagman, B. M., and Morlighem, M.: Radiostratigraphyand age
structure of the Greenland ice sheet, J. Geophys. Res.-Earth, 120,
212–241, https://doi.org/10.1002/2014JF003215,2015.
The Cryosphere, 12, 2969–2979, 2018
www.the-cryosphere.net/12/2969/2018/
https://data.cresis.ku.edu/https://doi.org/10.1029/2001jd900200https://doi.org/10.1364/JOSAA.20.000609https://doi.org/10.1109/tap.1983.1142979https://doi.org/10.1109/tgrs.2007.897433https://doi.org/10.3189/172756500781833070https://doi.org/10.3189/2014AoG67A005https://pdfs.semanticscholar.org/d361/3bc68c51ebcffbaa0b9aa6056c1cf8108adc.pdfhttps://pdfs.semanticscholar.org/d361/3bc68c51ebcffbaa0b9aa6056c1cf8108adc.pdfhttps://doi.org/10.1109/36.957274https://doi.org/10.1109/igarss.2000.859643https://doi.org/10.1029/2006jf000717https://doi.org/10.1002/2014JF003215
-
A. Heister and R. Scheiber: Processing of radio-echo sounding
data 2979
Marathe, K. C.: Dual-band multi-channel airborne radar for
map-ping the internal and basal layers of polar ice sheets,
Master’sthesis, University of Kansas, available at:
https://kuscholarworks.ku.edu/handle/1808/4036 (last access: 1 June
2018), 2008.
Matsuoka, K.: Pitfalls in radar diagnosis of ice-sheet bed
conditions:Lessons from englacial attenuation models, Geophys. Res.
Lett.,38, 1–5, https://doi.org/10.1029/2010gl046205, 2011.
Mishra, A., Yan, J.-B., Leuschen, C. J., and Gogineni,
S.:Multilook SAR processing and array optimization appliedto radio
echo ice sounding data, in: 2016 IEEE Interna-tional Symposium on
Phased Array Systems and Technol-ogy (PAST), Waltham, MA, USA,
18–21 October 2016,https://doi.org/10.1109/array.2016.7832586,
2016.
Peters, M., Blankenship, D., Carter, S., Kempf, S., Young, D.,
andHolt, J.: Along-track focusing of airborne radar sounding
datafrom West Antarctica for improving basal reflection analysisand
layer detection, IEEE T. Geosci. Remote, 45,
2725–2736,https://doi.org/10.1109/tgrs.2007.897416, 2007.
Scheiber, R., Prats, P., and Heliere, F.: Surface clutter
suppressiontechniques for ice sounding radars: analysis of airborne
data, in:Proceedings of the 7th European Conference on Synthetic
Aper-ture Radar, Friedrichshafen, Germany, 2–5 June 2008, 1–4,
2008.
Schroeder, D. M., Blankenship, D. D., and Young, D. A.:
Evi-dence for a water system transition beneath Thwaites
Glacier,West Antarctica, P. Natl. Acad. Sci., 110,
12225–12228,https://doi.org/10.1073/pnas.1302828110, 2013.
Schroeder, D. M., Blankenship, D. D., Young, D. A., Witus, A.
E.,and Anderson, J. B.: Airborne radar sounding evidence for
de-formable sediments and outcropping bedrock beneath
ThwaitesGlacier, West Antarctica, Geophys. Res. Lett., 41,
7200–7208,https://doi.org/10.1002/2014gl061645, 2014.
Siegert, M. J.: Ice flow direction change in inte-rior West
Antarctica, Science, 305,
1948–1951,https://doi.org/10.1126/science.1101072, 2004.
Sime, L. C., Hindmarsh, R. C., and Corr, H.: Auto-mated
processing to derive dip angles of englacialradar reflectors in ice
sheets, J. Glaciol., 57,
260–266,https://doi.org/10.3189/002214311796405870, 2011.
Smith, L.: Validation of CReSIS synthetic aperture radar
processorand optimal processing parameters, Master’s thesis,
Universityof Kansas, available at:
https://kuscholarworks.ku.edu/handle/1808/18381 (last access: 1
June 2018), 2014.
Wu, X., Jezek, K. C., Rodriguez, E., Gogineni, S.,
Rodriguez-Morales, F., and Freeman, A.: Ice sheet bed mapping with
air-borne SAR tomography, IEEE T. Geosci. Remote, 49, 3791–3802,
https://doi.org/10.1109/tgrs.2011.2132802, 2011.
www.the-cryosphere.net/12/2969/2018/ The Cryosphere, 12,
2969–2979, 2018
https://kuscholarworks.ku.edu/handle/1808/4036https://kuscholarworks.ku.edu/handle/1808/4036https://doi.org/10.1029/2010gl046205https://doi.org/10.1109/array.2016.7832586https://doi.org/10.1109/tgrs.2007.897416https://doi.org/10.1073/pnas.1302828110https://doi.org/10.1002/2014gl061645https://doi.org/10.1126/science.1101072https://doi.org/10.3189/002214311796405870https://kuscholarworks.ku.edu/handle/1808/18381https://kuscholarworks.ku.edu/handle/1808/18381https://doi.org/10.1109/tgrs.2011.2132802
AbstractIntroductionSAR focusingMultiple subbands
processingGreenland MCRDS dataEnhancement of SAR
echogramsScientific utility of large beamwidth SAR
processingSummary and conclusionsData availabilityCompeting
interestsAcknowledgementsReferences