Journal topic
The Cryosphere, 12, 1767–1778, 2018
https://doi.org/10.5194/tc-12-1767-2018
The Cryosphere, 12, 1767–1778, 2018
https://doi.org/10.5194/tc-12-1767-2018

Research article 25 May 2018

Research article | 25 May 2018

# Seasonal variations of the backscattering coefficient measured by radar altimeters over the Antarctic Ice Sheet

Seasonal variations of the backscattering coefficient measured by radar altimeters over the Antarctic Ice Sheet
Fifi Ibrahime Adodo1,2, Frédérique Remy1, and Ghislain Picard2 Fifi Ibrahime Adodo et al.
• 1Laboratoire d'Etudes en Géophysique et Oceanographie Spatiale (LEGOS), Centre National de la Recherche Scientifique (CNRS), Toulouse, 31400, France
• 2Institut des Géosciences de l'Environnement (IGE), Grenoble, 38402, Saint-Martin-d'Hères CEDEX, France

Abstract

Spaceborne radar altimeters are a valuable tool for observing the Antarctic Ice Sheet. The radar wave interaction with the snow provides information on both the surface and the subsurface of the snowpack due to its dependence on the snow properties. However, the penetration of the radar wave within the snowpack also induces a negative bias on the estimated surface elevation. Empirical corrections of this space- and time-varying bias are usually based on the backscattering coefficient variability. We investigate the spatial and seasonal variations of the backscattering coefficient at the S (3.2 GHz  9.4 cm), Ku (13.6 GHz  2.3 cm) and Ka (37 GHz  0.8 cm) bands. We identified that the backscattering coefficient at Ku band reaches a maximum in winter in part of the continent (Region 1) and in the summer in the remaining (Region 2), while the evolution at other frequencies is relatively uniform over the whole continent. To explain this contrasting behavior between frequencies and between regions, we studied the sensitivity of the backscattering coefficient at three frequencies to several parameters (surface snow density, snow temperature and snow grain size) using an electromagnetic model. The results show that the seasonal cycle of the backscattering coefficient at Ka frequency is dominated by the volume echo and is mainly driven by snow temperature evolution everywhere. In contrast, at S band, the cycle is dominated by the surface echo. At Ku band, the seasonal cycle is dominated by the volume echo in Region 1 and by the surface echo in Region 2. This investigation provides new information on the seasonal dynamics of the Antarctic Ice Sheet surface and provides new clues to build more accurate corrections of the radar altimeter surface elevation signal in the future.

1 Introduction

Radar altimeters are the most widely used sensors for measuring the surface elevation of the polar ice sheets (Remy et al., 1999; Allison et al., 2009). It is a valuable tool for monitoring and quantifying the volume change of the Antarctic Ice Sheet (AIS) (Zwally et al., 2005; Wingham et al., 2006; Flament and Rémy, 2012; Helm et al., 2014). However, altimetric observations are affected by several errors: error due to atmospheric or ionospheric propagations, slope error and error due to the radar wave penetration into the cold and dry snow (Ridley and Partington, 1988). The first two errors are usually corrected with good accuracy (Remy et al., 2012; Nilsson et al., 2016), while the last one is the most critical and the most challenging problem to tackle (Remy et al., 2012) as it results in an overestimation of the observed distance between the satellite and the target, leading to a negative bias in the surface elevation estimation. The magnitude of the penetration error on the estimated surface elevation is between a few tens of centimeters and few meters (Remy and Parouty, 2009). For instance, Michel et al. (2014) have found a surface elevation difference of 0.5 m between ENVIronment SATellite (ENVISAT) and ICESat crossover points over Antarctica. Authors relate this negative bias to the difference in the penetration depth between the radar altimeter wave that penetrates within the snowpack and the laser altimeter beam that does not penetrate within the snowpack. The temporal variation in the penetration error is therefore critical for accurate interpretation of ice sheet volume changes (Remy et al., 2012).

The radar wave interaction with snow provides information on the snowpack surface and subsurface properties, but it complicates the altimetric signal interpretation because the latter would be sensitive to many more snow parameters than if the signal only came from the surface. To clarify the impacts of snow parameters on the backscattering coefficient, this paper investigates the spatial and seasonal variations of the radar backscattering coefficient at the S, Ku and Ka bands. To this end, electromagnetic models are used to assess the backscattering coefficient sensitivity to snow properties at the three frequencies. The aim of this paper is to determine the prevailing snow parameters that drive the seasonal cycle of the observed backscattering coefficient at different radar frequencies and locations over the AIS.

The ENVISAT carries two radar altimeter sensors (RA-2) that operate at 13.6 GHz (Ku band  2.3 cm) and 3.2 GHz (S band  9.4 cm). The S band was originally intended for ionospheric corrections while the Ku band provides more accurate surface elevation due to the lower penetration depth. Comparison of the altimetric waveform characteristics between the Ku and S bands revealed different seasonal variations over the AIS (Lacroix et al., 2008b). The dual-frequency information can therefore be useful for retrieving information on snowpack properties. The launch in March 2013 of the radar altimeter Satellite for ARgos and ALtiKa (SARAL)/AltiKa that operates at the Ka band (37 GHz  0.8 cm) and had the same 35-day phased orbit as ENVISAT until March 2016 allowed comparisons with much higher frequencies for the first time. Temporal variations of the estimated surface elevation with respect to the backscattering coefficient are 6 times lower at the Ka band than that of the Ku band, which implies that the volume echo at the Ka band comes from the near subsurface (< 1 m) and is mostly controlled by ice grain size and temperature (Remy et al., 2015).

2 Data and methods

## 2.1 Altimetric observations

Radar altimeters data were acquired by ENVISAT launched on March 2002 by the European Space Agency (ESA). Acquisitions are simultaneous at the S and Ku bands, every 330 m along track on a 35-day repeat cycle orbit from September 2002 to October 2010 (the end of its repeat cycle orbit). The S band sensor failed after 5 years of measurements. The satellite footprint has typically a 5 km radius and no data were acquired above 81.5 S due to its orbit maximum inclination. The range gate resolution is about 94 and 47 cm at the S and Ku bands, respectively.

To ensure a long and homogeneous time series of post-ENVISAT missions and to complement the Ocean Surface Topography Mission (OSTM)/Jason (Steunou et al., 2015), the SARAL/AltiKa was launched on 25 February 2013 by a joint CNES-ISRO (Centre National d'Etudes Spatiales – Indian Space Research Organisation) mission, on the same 35-day repeat cycle orbit as ENVISAT. On March 2016, SARAL/AltiKa orbit was shifted onto a new orbit. Unlike classical Ku band radar altimeter, the SARAL/AltiKa altimeter operates at the Ka band (37 GHz  0.8 cm) and has a range gate resolution of 30 cm. The ICE-2 retracking process was applied to the Ka, Ku and S band waveforms, allowing estimation of the range, the backscattering coefficient (σ0), the leading edge width and the trailing edge slope. The difference between the Ka and Ku bands and between the Ka and S bands is up to a factor of 2.7 and 11.6, respectively, which results in different sensitivity to the surface and the subsurface characteristics.

The ENVISAT and AltiKa datasets used in this study were averaged at a 1 km scale on the ENVISAT nominal orbit. We processed 84 cycles of the backscattering coefficient from October 2002 until September 2010 for the Ku band and 55 cycles from October 2002 until December 2007 for the S band. Moreover, we consider 3 years of AltiKa altimeter data from March 2013 to March 2016, i.e., a total of 32 cycles of the backscattering coefficient over the whole Antarctic continent.

Figure 1Time series of the backscattering coefficient at the S (blue), Ku (black) and Ka (red) bands at location (69.468 S, 134.28 E) between October 2002 and December 2007 at S band, October 2002 and September 2010 at Ku band and March 2013 and March 2016 at Ka band. The dashed lines represent the best fits to the time series (see Eq. 1). The observations show the seasonal cycle with a 1-year period at the different frequencies.

## 2.2 Amplitude and date of maximum backscattering coefficient in the seasonal cycle

The amplitude and the date at which the backscattering coefficient (σ0) reaches a maximum within a seasonal cycle were calculated at the S, Ku and Ka bands for the entire Antarctic continent. Figure 1 shows an example of the temporal evolution of σ0 at a location (69.46 S, 134.28 E) at the three frequencies. The time series of σ0 exhibit a clear and well-marked cycle with a 1-year period (called seasonal cycle hereafter). The amplitude and the phase of the seasonal cycle of σ0 were computed by fitting the observations with the following model Eq. (1):

$\begin{array}{}\text{(1)}& {\mathit{\sigma }}_{i}^{\mathrm{0}}\left(t\right)={\mathit{\alpha }}_{i}\mathrm{sin}\left(\mathrm{2}\mathit{\pi }\frac{t}{T}\right)+{\mathit{\beta }}_{i}\mathrm{cos}\left(\mathrm{2}\mathit{\pi }\frac{t}{T}\right)+{C}_{i},\end{array}$

with ${A}_{i}=\sqrt{{\mathit{\alpha }}_{i}^{\mathrm{2}}+{\mathit{\beta }}_{i}^{\mathrm{2}}}$ and Φi = arctan (βiαi), where Ai and Φi is the amplitude and the phase of the seasonal cycle of ${\mathit{\sigma }}_{i}^{\mathrm{0}}$, respectively, deduced from constants αi and βi returned by the model; T= 365 days; t ranges from 0 to 5 years for the S band, from 0 to 8 years for the Ku band and from 0 to 3 years for the Ka band, with steps of 35 days; and i represents the data point over the continent. The fit was done with the ordinary least squares method and all data points with time series length less than 11 cycles (about a year) were discarded. The date at which σ0 reaches a maximum within a seasonal cycle is obtained by converting the seasonal phase Φi to fraction of a year (assuming a year counts for 360 days). We have found that along-track analysis of the seasonal parameters of σ0 showed no dependence to anisotropic effects. In the following, both ascending and descending measurements are mixed to keep a high density of observations and cover most AIS ( 1.9 million data points). For visualization needs, seasonal parameters are interpolated on a map of 5 km × 5 km grid by averaging with Gaussian weights. We considered all data points within a 25 km radius and weighted with a decorrelation radius of 10 km.

## 2.3 Backscattering coefficient modeling

To explore the snowpack properties that drive the seasonal cycle of σ0, we investigated its sensitivity to the snowpack surface and subsurface properties using an altimetric echo model of snow. This model accounts for the surface and the volume echoes. The surface echo results from the interactions of the radar wave with the snow surface (air–snow interactions) while the volume echo results from the interactions of the radar wave with the scatterers within the snowpack (snow–snow interactions). The physics involved in both surface and volume echoes have been previously studied by Lacroix et al. (2008b).

### 2.3.1 Surface echo modeling

The snow surface can be modeled as a randomly rough surface because most naturally occurring surfaces are irregular. The surface scattering coefficient from a rough surface is thus controlled by the effective dielectric constant of the medium and the surface roughness characteristics (Ulaby et al., 1982; Fung, 1994). The effective dielectric constant of the snow is a function of the snow density and the ice dielectric constant, while the roughness is usually modeled by two parameters: the surface correlation length (l) and the standard deviation of the surface elevation (σh) (Ulaby et al., 1982). In the case of large standard deviations of the surface elevation (σh) (compared to the radar wavelength), the backscattering coefficient from a rough surface ${\mathit{\sigma }}_{\mathrm{sur}}^{\mathrm{0}}$ can be estimated assuming a Gaussian autocorrelation function (Ulaby et al., 1982):

$\begin{array}{}\text{(2)}& {\mathit{\sigma }}_{\mathrm{sur}}^{\mathrm{0}}=\frac{{\left|R\left(\mathrm{0}\right)\right|}^{\mathrm{2}}}{\mathrm{2}{S}^{\mathrm{2}}},\end{array}$

where R(0) is the Fresnel reflection coefficient at the normal incident angle and $S=l/{\mathit{\sigma }}_{\mathrm{h}}$ the root mean square (RMS) of the surface slope at the radar wavelength scale. Equation (2) is almost independent of the radar wave frequency and ${\mathit{\sigma }}_{\mathrm{sur}}^{\mathrm{0}}$ increases with increasing surface snow density and decreasing surface slope RMS. Surface snow density variations from 300 to 400 kg m−3 induce a variation of ±2.17 dB in the surface echo.

### 2.3.2 Volume echo modeling

The volume echo is mainly controlled by the scattering coefficient (Ks), depending on the size of the scatterers and the radar frequency. The power extinction in the snowpack is the sum of the scattering coefficient (Ks) and the absorption (Kab) coefficient. The latter depends on snow temperature and radar frequency. In the following, the scatterers are assumed to be spherical. Ks and Kab are given by Mätzler (1998):

$\begin{array}{}\text{(3)}& & {K}_{\mathrm{s}}=\frac{\mathrm{3}}{\mathrm{32}}{p}_{\mathrm{c}}^{\mathrm{3}}{K}_{\mathrm{0}}^{\mathrm{4}}\mathit{\nu }\left(\mathrm{1}-\mathit{\nu }\right){\left({\mathit{ϵ}}_{i}^{\prime }-\mathrm{1}\right)}^{\mathrm{2}}{K}_{\mathrm{d}}^{\mathrm{2}},\text{(4)}& & {K}_{\mathrm{ab}}={K}_{\mathrm{0}}\mathit{\nu }{\mathit{ϵ}}_{i}^{\prime \prime }{K}_{\mathrm{d}}^{\mathrm{2}},\end{array}$

where ${K}_{\mathrm{0}}=\mathrm{2}\mathit{\pi }/\mathit{\lambda }$ is the wave number and λ the wavelength, ν is the fractional volume of the scatterers, ${\mathit{ϵ}}_{i}^{\prime }$ and ${\mathit{ϵ}}_{i}^{\prime \prime }$ are the real and imaginary parts of the effective dielectric constant of pure ice, ${p}_{\mathrm{c}}=\left(\mathrm{4}{r}_{\mathrm{g}}\right)/\mathrm{3}$ (Mätzler, 1998) is the correlation length (used here as the effective size parameter) with rg the scatterers radius and ${K}_{\mathrm{d}}^{\mathrm{2}}={\left|\mathrm{2}{\mathit{ϵ}}^{\prime }+\mathrm{1}\right|}^{\mathrm{2}}/{\left|\mathrm{2}{\mathit{ϵ}}^{\prime }+{\mathit{ϵ}}_{i}^{\prime }\right|}^{\mathrm{2}}$ with ϵ the real part of the effective dielectric constant of snow (Tiuri et al., 1984).

For snow grain radius increasing from 0.3 to 0.5 mm, Ks increases from 1.05 to 4.85 m−1 at the Ka band, from 0.02 to 0.08 m−1 at the Ku band and from 0.58 × 10−4 to 2.7 × 10−4 m−1 at the S band. As snow temperature varies from 220 to 250 K, Kab increases from 0.194 to 0.287 m−1 at the Ka band, from 0.026 to 0.039 m−1 at the Ku band and from 0.002 to 0.003 m−1 at the S band. The extinction coefficient at the Ka band is dominated by the scattering coefficient. In contrast, the losses by absorption dominate the extinction at the S band while at the Ku band, both coefficients are of the same order of magnitude. Volume scattering mainly affects the Ka and Ku bands. Finally, the losses by absorption increase with snow temperature while the scattering coefficient is mainly driven by snow grain size. Both the losses by absorption and scattering coefficient increase with increasing radar frequency.

## Snow property profiles

For the simulations, we considered the same vertical density profile as Lacroix et al. (2008b) with a variation only in the top 10 m given by

$\begin{array}{}\text{(5)}& \mathit{\rho }\left(z\right)={\mathit{\rho }}_{\mathrm{0}}+pz+{c}_{\mathrm{2}}{z}^{\mathrm{2}}+{c}_{\mathrm{3}}{z}^{\mathrm{3}},\end{array}$

where c2 and c3 are constant values taken from the Talos Dome density profile: 1.35 × 10−4 and 5.86 × 10−7, respectively (Frezzotti et al., 2004). ρ0 is the mean surface density and p= 1.40 × 10−2 is calculated as a function of ρ so that the density at the depth below the surface z=10 m is the density measured at the Talos Dome (72.78 S, 159.06 E). The choice of the vertical density profile has a negligible effect on the results of the sensitivity test. Snow temperature is computed using the solution of the thermal diffusion equation (e.g Bingham and Drinkwater, 2000; Surdyk, 2002), assuming a sinusoidal seasonal surface temperature and constant snow thermal diffusivity κ. The temperature at depth z is of the form

$\begin{array}{}\text{(6)}& T\left(z,t\right)={A}_{\mathrm{m}}\text{exp}\left(\frac{-z}{l}\right)\mathrm{cos}\left(\mathit{\omega }t-\frac{z}{l}\right)+{T}_{\mathrm{m}},\end{array}$

where Am and Tm are the seasonal amplitude and mean temperatures, respectively, ω is the angular frequency, t is the time, z is the depth and $l=\sqrt{\mathrm{2}\mathit{\kappa }/\mathit{\omega }}$. κ is the ratio of the thermal conductivity (κd) to the heat capacity and the snow density (ρ). We used the quadratic relationship of the thermal conductivity derived by Sturm et al. (1997): ${\mathit{\kappa }}_{\mathrm{d}}=\mathrm{0.138}-\mathrm{1.01}\mathit{\rho }+\mathrm{3.233}{\mathit{\rho }}^{\mathrm{2}}$. In the computing of κd and κ, snow density, ρ, is assumed equal to an average of the density profile of Eq. (5) (Bingham and Drinkwater, 2000). The temperature wave propagating in the snowpack has decreasing amplitude with respect to depth. The snow grain growth rate is mainly dependent on snow temperature (Brucker et al., 2010) and the snow grain profile with depth (Bingham and Drinkwater, 2000) is expressed by

$\begin{array}{}\text{(7)}& {r}_{\mathrm{g}}{\left(z\right)}^{\mathrm{2}}={r}_{\mathrm{0}}^{\mathrm{2}}+{K}_{\mathrm{g}}z/\mathit{\pi }D,\end{array}$

where Kg=0.00042 mm2 yr−1 is the typical snow grain growth rate, D is the mean annual snow accumulation (mm yr−1), z is the depth and r0 is the spherical scatterer mean radius at the surface. Tests of variation of D show no significant effect on the volume echo trend, and we therefore set D to 50 mm yr−1 (Bingham and Drinkwater, 2000).

Figure 2Spatial distribution of the seasonal date of maximum backscattering coefficient at the S (a), Ku (b) and Ka (c) bands. Black contour lines delineate regions where the backscattering coefficient at the Ku band peaks before April. Blue defines a maximum in the winter while the magenta a maximum in the summer. The cross mark represents the location of the time series shown in Fig. 1. White areas indicate regions where no observations are available (latitudinal orbit limit of 81.5 S). Color bar is cyclic and defines Julian days.

Figure 3Histogram of the seasonal date of maximum backscattering coefficient at the S (blue), Ku (black) and Ka (red) bands. The gray bars represent periods referred to as summer (January to April) and winter (June to September).

Figure 4Difference of the seasonal date of maximum backscattering coefficient between the Ku and Ka bands. Blue defines a maximum in the Ka band before the Ku band while the magenta the inverse. Black contour lines delineate regions where the backscattering coefficient at the Ku band peaks before April. The cross mark represents the location of the time series shown in Fig. 1. White areas indicate regions where no observations are available (latitudinal orbit limit of 81.5 S). The color bar is cyclic and defines the Julian days.

Figure 5Spatial distribution of the seasonal amplitude of the backscattering coefficient at the S (a), Ku (b) and Ka (c) bands. Black contour lines delineate regions where the backscattering coefficient at the Ku band peaks before April. The cross mark represents the location of the time series shown in Fig. 1. White areas indicate regions where no observations are available (latitudinal orbit limit of 81.5 S). Values are expressed in dB.

3 Results

## 3.1 Spatial patterns of the amplitude and date of maximum backscattering coefficient

The spatial distribution and the histogram of the seasonal date of maximum σ0 at the S, Ku and Ka bands are shown in Figs. 2 and 3, respectively. Among the three bands, the Ku band presents the most contrasted geographical patterns. In the zone that appears in magenta, the seasonal cycle of σ0 reaches a maximum early in the year (summer peak zone, SP hereafter). This zone covers the eastern-central part of the AIS, which encompasses the domes and high-altitude regions (> 3000 m a.s.l.). It extends from Wilkes Land to Dronning Maud Land and is characterized by a decrease in σ0 from late fall to early spring followed by an increase at the end of the summer. The zone appearing in blue (hereafter winter peak zone, WP), encompasses the lower regions (< 3000 m a.s.l.) including coastal steeply sloped regions. It is characterized by an increase in σ0 from late fall to early spring. In contrast to the Ku band, the seasonal cycles of σ0 over the AIS are generally maximum in the summer at the S band but maximum in the winter at the Ka band. In Fig. 3, the Ku band date of maximum σ0 histogram is clearly bimodal with peaks between Julian days 1 and 100 (1 January to mid-April) and between Julian days 175 and 275 (June to September). In the following, these two periods are referred to as summer and winter, respectively. With these definitions, the WP and SP represent 42 and 45 % of the observed area, respectively. The histogram of the date of maximum σ0 at the S and Ka bands are unimodal with a peak in summer for a lower frequency (WP: 11 %, SP: 66 %, using the summer and winter periods previously defined) and a peak in winter for a higher frequency (WP: 50 %, SP: 14 %, using the summer and winter periods previously defined). The difference of the seasonal date of maximum σ0 between the Ku and Ka bands (Fig. 4), over the AIS, shows a geographical pattern similar to that observed in Fig. 2b. Negative values indicate that σ0 is maximum at the Ku band before the Ka band while positive values indicate the opposite. Negative values account for about 36 % of the observations and coincide with the SP where σ0 is maximum in summer at the Ku band. Positive values, the zone appearing in blue, cover 48 % of the AIS and coincide with the WP. Hence, we note a positive lag of the date of maximum σ0 between the Ku and Ka bands only in the zone where σ0 is maximum in the winter in both frequencies and a negative lag in the other zones. The spatial distribution of the seasonal amplitude of σ0 at the Ka band (Fig. 5c) shows an obvious geographical pattern close to that of the seasonal date of maximum σ0 at the Ku band. The Ka band seasonal amplitude of σ0 is the highest in the WP (1.02 ± 0.56 dB) and weakest in the SP (0.53 ± 0.41 dB) as shown in Fig. 6. By contrast, the seasonal amplitude of σ0 at the S band (Fig. 5a) appears anticorrelated with that at the Ka band, exhibiting a large seasonal amplitude in the SP (0.79 ± 0.40 dB) and a weak amplitude in the WP (0.42 ± 0.28 dB). The seasonal amplitude of σ0 in the SP is almost twice as large as that of the WP at the S band and the inverse is true at the Ka band. The seasonal amplitude of σ0 at the Ku band shows no evident regional patterns and is almost of the same magnitude in both zones (Fig. 5b), except in the interior of Wilkes Land, Princess Elisabeth Land and the Ronne Ice Shelf, which showed greater amplitudes.

Figure 6Mean seasonal amplitude with respect to the date of maximum backscattering coefficient at the S (blue), Ku (black) and Ka (red) bands. The gray bars represent periods referred to as summer (January to April) and winter (June to September).

## 3.2 Temporal variations of the surface elevation with respect to the backscattering coefficient

Figure 7 shows the spatial distribution of temporal variations of the estimated surface elevation residuals with respect to σ0 residuals at the Ku band, hereafter denoted dh dσ0. The surface elevation was indirectly estimated from the retracked range (computed with the ICE-2 retracker) at each data point and was corrected for atmospheric errors. dh and dσ0 were derived by subtracting the mean value from the time series of the elevation and backscatter, respectively. dh dσ0 represents the correlation gradient or the slope at each data points over the AIS. Negative values of dh dσ0 indicate that surface elevation decreases when σ0 increases, implying that temporal variations in σ0 are due to changes in the deep snowpack properties, i.e., in the volume echo. In fact, the inverse relationship between surface elevation and σ0 is related to a greater backscatter from depth that shifts more power to greater delay times in the received waveform, thus increasing the retracked range and decreasing the estimated elevation (Armitage et al., 2014). In contrast, positive values of dh dσ0 indicate that the surface elevation increases with σ0. In this case, the temporal variations of σ0 are related to changes in the surface echo. The map in Fig. 7 shows that near-zero and negative values of dh dσ0 (in blue) are found in the WP. This means that the WP undergoes large variations of volume echo.

Figure 7Temporal variations of the surface elevation residuals with respect to the backscattering coefficient residuals at the Ku band (denoted hereafter, dh dσ0). Black contour lines delineate regions where the backscattering coefficient at the Ku band peaks before April. The cross mark represents the location of the time series shown in Fig. 1. White areas indicate regions where no observations are available (latitudinal orbit limit of 81.5 S). Values are expressed in m dB−1.

## 3.3 Sensitivity test

Since there are few, if any, studies on the seasonal cycle of snow surface roughness, it is poorly known. The sensitivity study of the surface echo is thus limited by the lack of information on snow surface roughness, in particular over the AIS. Consequently, we have focused on the modeling of the seasonal cycle of the volume echo. In this subsection, the sensitivity test of the volume echo at the S, Ku and Ka bands to snow properties is explored considering three parameters snow temperature, snow grain size and snow density in the analysis of the seasonal cycle of σ0.

The model shows an increase in the volume echo with snow density at the three frequencies (Fig. 8a). Snow density controls the thermal conductivity of the medium. Increasing surface snow density increases thermal diffusivity, which attenuates the propagation of the temperature wave in the snowpack. Figure 8b and c show that the volume echo at the S band is not sensitive to snow temperature and grain size variations, while the volume echo at the Ku and Ka bands is affected by both parameters. Snow density, temperature and grain size impacts on the volume echo are more significant at the Ka band than at the Ku and S band levels. The volume echo increases with the snow density at the three frequencies, and at the S band the volume echo is less significant.

Figure 8Sensitivity tests of the volume echo with respect to the surface snow density (a), snow temperature (b) and snow grain size (c) at the S (blue), Ku (black) and Ka (red) bands.

Figure 9Distribution of the date of maximum backscattering coefficient at the Ku band superimposed on the RADARSAT mosaic (RAMP). Blue contour lines show the boundaries between the WP and the SP over the Antarctica Ice Sheet. SPs are regions where the backscattering coefficient reaches a maximum in summer (inset of the contours), and WPs are regions where the backscattering coefficient reaches a maximum in winter (where snow surface features are apparent). No observations are available beyond 81.5 S (black circle).

4 Discussion

The sensitivity of the volume echo to snow temperature shown in Fig. 8b implies that the volume echo is maximum in winter at the Ku and Ka bands and constant at the S band. This sensitivity is explained by the fact that increasing snow temperature increases absorption, resulting in a decrease of the radar wave penetration in the medium and thus limiting the volume echo. Also, increasing snow grain size increases the scattering coefficient, which in turn increases the radar wave extinction in the medium. This results in a decrease of the radar wave penetration and therefore may limit the volume echo. Moreover, the positive lag observed between the Ku and Ka bands in the WP in Fig. 4 can be explained by the difference of the radar wave penetration depth between the Ku ( 10 m) and Ka (> 1 m) bands in the snowpack. This lag is related to the propagation of the temperature gradient from the surface into the snowpack. As the temperature controls the snow grain metamorphism and the radar wave penetration depth, the variation in the volume echo would be predominantly driven by the seasonal variations of snow temperature.

Snow density is involved in both the surface and volume echoes. The magnitude of these echoes increase with increasing surface snow density, and thus similar seasonal cycle of σ0 would be expected at any frequency if snow density were the main driver. This is in contrast to the observations (Fig. 3). Therefore, the seasonal cycle of σ0 cannot be explained solely by snow density. Being insensitive to snow temperature and grain size (Fig. 8b, c), the observed seasonal cycle of σ0 at the S band cannot be explained by the volume echo. This implies that snow surface properties (surface snow density and roughness) are the main factors driving the seasonal cycle of σ0 at the S band.

Figure 10Seasonal wind speed amplitude (a) and average (b). Data are extracted from ERA-Interim reanalysis provided by ECMWF on 25×25 km2 grid cells, from the period 2002 to 2010 corresponding to the ENVISAT lifetime. Black contour lines delineate regions where the backscattering coefficient at the Ku band peaks before April. The white star shows the location of the time series plotted in Fig. 1. No observations are available beyond 81.5 S (black dotted circle).

From the S to Ka band, the radar wavelength decreases by a factor of 12 from 9.4 to 0.8 cm corresponding to a scale change from centimeter to millimeter. The scale at which the surface roughness plays a role in radar backscattering coefficient depends on the radar wavelength (Ulaby et al., 1982). On a rough surface, the surface scattering consists of two components: the coherent and incoherent scattering (Ulaby et al., 1982). The former is the scattered component in the specular direction while the latter is the scattered component in all directions. As the radar wavelength is shortened to less than a centimeter, the surface appears rougher and the surface coherent component vanishes (Ulaby et al., 1982). The surface incoherent component magnitude is small and thus is concealed by the volume scattering, which consists of only incoherent scattering. The backscattering coefficient at a smaller wavelength or on a rougher surface consists of only incoherent components and therefore appears as a volume-scattering medium. Simulations in Fig. 8 emphasize this contention, showing a greater amplitude of the volume echo at higher frequencies. We can therefore argue that the seasonal cycle of the observed σ0 at the Ka band is governed by the volume echo. This explains the peak of the observed σ0 in the winter at the Ka band over the AIS.

Several observations show that sastrugi (10 cm to 1 m height) are the main contributors to surface roughness (Kotlyakov, 1966; Inoue, 1989; Lacroix et al., 2007). Since the biggest features (hectometer to kilometer scales) change little over time, it is likely that the most influential roughness scale in the seasonal cycle of the surface echo is the sastrugi (Lacroix et al., 2008a). Despite the increase in magnitude of the surface and volume echoes with surface snow density, evidence from Fig. 3 suggests that the seasonal cycle of σ0 cannot be explained by the seasonal cycle in surface snow density. Therefore, it is likely that the seasonal cycle of the observed σ0 at the S band, predominantly driven by the surface echo, stems from the seasonal cycle of the snow surface roughness. There is no field observation that confirms this fact, but our findings suggest that such information would help to understand the altimetric signal in the future. However, in this study it is difficult to differentiate with certainty between the surface snow density and the snow surface roughness, which drives the seasonal cycle of the surface echo. There are three main reasons for this: (i) The snow surface roughness is poorly known, in particular its seasonal variability; (ii) surface snow properties evolve rapidly with the wind; and (iii) the relation between the surface snow roughness and density is complex because both variables are interdependent. The denser the snow surface, the larger the effect of surface roughness. This amplification is due to the increase of the effective dielectric discontinuity with density (Fung, 1994).

Considering that σ0 at the Ku band shows two opposing seasonal cycle patterns over the AIS and its wavelength is between that of the S and Ka bands, we suggest that σ0 at the Ku band is dominated by the seasonal cycle of the surface echo, similar to the S band in the SP, and by the seasonal cycle of the volume echo, similar to the Ka band in the WP. We support this hypothesis with ancillary data and by modeling. By overlaying the Antarctica RADARSAT mosaic with the SP boundaries (Fig. 9), we find that the WP matches regions of large heterogeneous backscatter from RADARSAT, where megadunes (Frezzotti et al., 2002) and wind-glazed surfaces (Scambos et al., 2012) have been observed. The seasonal cycle of σ0 at the Ku band is maximum in the winter in heterogeneous RADARSAT backscatter regions while it is maximum in the summer in the other regions. In fact, areas of megadunes are characterized by slightly steeper regional slope and the presence of highly persistent katabatic winds (Frezzotti et al., 2002) and wind-glazed surfaces have been formed by persistent katabatic winds in areas of megadunes (Scambos et al., 2012). There exists therefore a relationship between the wind and the seasonal cycle of σ0. To further investigate this point, we used ERA-Interim reanalysis wind speed data supplied by ECMWF (European Centre For Medium-Range Weather Forecasts) on the period 2002 to 2010, corresponding to that of the Ku band. Equation (1) is used to compute the seasonal characteristics of the wind speed by replacing σ0 with the wind speed. A visual inspection shows a high spatial coherence of the seasonal amplitude of the wind speed (Fig. 10a) patterns with the date of maximum σ0 over the seasonal cycle at the Ku band (Fig. 2b). Wind speed average (8.2 ± 1.6 m s−1) and seasonal amplitude (1.7 ± 0.4 m s−1) are higher in the WP than in the SP (6.6 ± 1.58 m s−1 and 1.0 ± 0.3 m s−1, respectively).

The striking similarity in the spatial distribution of the seasonal amplitude of σ0 at the Ka band (Fig. 5c) and the seasonal date of maximum σ0 at the Ku band (Fig. 2b), which is itself correlated to the seasonal amplitude of the wind speed (Fig. 10a), suggests that the wind plays a significant role in the spatial distribution of the seasonal amplitude of σ0 at the Ka band. Although the wind effects on the snowpack are numerous and complex, we retained two for which we simulated the impacts on the volume echo (Fig. 8):

• a.

Wind may smash snow grains so that the surface snow density increases with wind speed (Male, 1980); this leads to an enhancement of the volume echo at the three frequencies as shown in Fig. 8a. Surface snow density is a good candidate for explaining the spatial distribution of the seasonal amplitude of σ0 at the Ka band because snow compaction can occur at different times of the year depending on the snow accumulation rate and the temperature gradient (Li and Zwally, 2002, 2004).

• b.

Increasing wind speed leads to an increase in snow erosion and transport that removes all or almost all the precipitated or wind deposited snow that may temporarily accumulate (Scambos et al., 2012; Lenaerts et al., 2012). This implies that there is no significant change in the surface mass balance over an annual cycle, i.e., near-zero net accumulation (Scambos et al., 2012), allowing snow surface to be almost constant and smooth. This corroborates our contention that the seasonal variation of the observed σ0 at the Ku band in the WP emanates exclusively from the volume echo (i.e., a greater backscatter from depth). Thus, it is presumed that these variations are due to depth hoar formation during winter in the WP. Indeed, the wind speed is on average maximum between Julian days 170 and 230 (June to August), when air temperature is colder than the snow temperature. Cold and persistent winds may unusually accelerate the cooling of the surface snow temperature (Remy and Minster, 1991). This causes an important temperature gradient, which determines the rate of metamorphism of snow grains within the snowpack. This specific increase of the temperature gradient would promote the formation of depth hoar in winter (Champollion et al., 2013), which creates coarse cup-shaped ice crystals, acts as more effective volume scatterers and hence increases the volume echo magnitude as predicted in Fig. 8c. For instance, Brucker et al. (2010) have found the highest vertical gradient in grain size, obtained over a multiyear average from 1987 to 2002, in the regions of the WP.

Finally, the combined effects of wind speed and temperature may explain the observed difference between the seasonal cycle of σ0 at the Ka and Ku bands. Similarly, the spatial distribution of the seasonal amplitude of σ0 at the Ka band is ascribed to the wind effects mentioned above on the snowpack.

5 Conclusions

This study, using 35-day repeat radar altimetry data, carries out spatial and temporal comparative analysis of the seasonal amplitude and date of maximum σ0 at the S, Ku and Ka bands. We used an 8-year time series of σ0 for the Ku band, a 5-year time series of σ0 for the S band and a 3-year time series of σ0 for the Ka band, covering the time period of 2002 to 2010 for ENVISAT sensors and 2013 to 2016 for the SARAL/AltiKa sensor. The backscattering coefficient shows seasonal variations with varying amplitude and phase over the AIS and with a marked dependence on radar frequency. In general, it is maximum in winter at the Ka band and maximum in summer at the S band. At the Ku band, both behaviors are found with a maximum in the winter in the so-called WP and a maximum in the summer in the SP.

We investigated snow properties that dominate the seasonal changes in the volume echo with electromagnetic models of the backscattering coefficient. As a result, we showed that variations in snow properties, such as temperature and grain size, cannot explain the seasonal cycle of σ0 observed at the S band due to its small sensitivity to those parameters. In contrast, the temperature cycle reasonably explain the seasonal cycle of the observed σ0 at the Ka band. We explain that the contrasted seasonal cycle of the observed σ0 at the Ku band is due to its high sensitivity to the volume echo in the WP and to the surface echo in the SP. The geographical patterns of the WP and SP are related to the seasonal amplitude of the wind speed. This is a result of the presence or lack of wind-glazed surfaces, induced by strong and persistent winds in the megadune areas.

This investigation provides new information on the Antarctic Ice Sheet surface seasonal dynamics and provides new clues to build a robust correction of the altimetric surface elevation signal. Multi-frequency sensors are the key for improving the understanding of the physics of radar altimeter measurements over the AIS. An important limitation of this study is the lack of information on the seasonal variability of the snow surface roughness in Antarctica, which will be the topic of future work.

Data availability
Data availability.

Data can be access on request on the CTOH and AVISO websites. ENVISAT datasets at (http://ctoh.legos.obs-mip.fr/products, last access: 23 May 2018) and SARAL/AltiKa datasets are downloadable from AVISO website (https://aviso-data-center.cnes.fr/\#altika_LG.htm, last access: 23 May 2018).

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

This work is a contribution to the ASUMA (improving the Accuracy of the SUrface Mass balance of Antarctica) project funded by the Agence Nationale de la Recherche, contract ANR-14-CE01-0001-01. ENVISAT and AltiKa data were provided by the Center for Topographic studies of the Oceans and Hydrosphere (CTOH) at LEGOS and are available at http://ctoh.legos.obs-mip.fr/, last access: 23 May 2018. The authors would like to thank Etienne Berthier and Jessica Klar from LEGOS for their helpful comments and suggestions. We are grateful to the anonymous reviewers and the editor, whose comments significantly improved the manuscript.

Edited by: Robert Arthern
Reviewed by: two anonymous referees

References

Allison, I., Alley, R. B., Fricker, H. A., Thomas, R. H., and Warner, R. C.: Ice sheet mass balance and sea level, Antarct. Sci., 21, 413–426, 2009.

Armitage, T. W. K., Wingham, D. J., and Ridout, A. L.: Meteorological Origin of the Static Crossover Pattern Present in Low-Resolution-Mode CryoSat-2 Data Over Central Antarctica, IEEE Geosci. Remote Sens. Lett., 11, 1295–1299, https://doi.org/10.1109/LGRS.2013.2292821, 2014.

Arthern, R. J., Wingham, D. J., and Ridout, A. L.: Controls on ERS altimeter measurements over ice sheets: Footprint-scale topography, backscatter fluctuations, and the dependence of microwave penetration depth on satellite orientation, J. Geophys. Res.-Atmos., 106, 33471–33484, 2001.

Bingham, A. W. and Drinkwater, M. R.: Recent changes in the microwave scattering properties of the Antarctic ice sheet, IEEE Trans. Geosci. Remote Sens., 38, 1810–1820, https://doi.org/10.1109/36.851765, 2000.

Brenner, A. C., DiMarzio, J. P., and Zwally, H. J.: Precision and accuracy of satellite radar and laser altimeter data over the continental ice sheets, IEEE Trans. Geosci. Remote Sens., 45, 321–331, 2007.

Brown, G.: The average impulse response of a rough surface and its applications, IEEE Trans. Antennas Propag., 25, 67–74, 1977.

Brucker, L., Picard, G., and Fily, M.: Snow grain-size profiles deduced from microwave snow emissivities in Antarctica, J. Glaciol., 56, 514–526, 2010.

Champollion, N., Picard, G., Arnaud, L., Lefebvre, E., and Fily, M.: Hoar crystal development and disappearance at Dome C, Antarctica: observation by near-infrared photography and passive microwave satellite, The Cryosphere, 7, 1247–1262, https://doi.org/10.5194/tc-7-1247-2013, 2013.

Flament, T. and Rémy, F.: Dynamic thinning of Antarctic glaciers from along-track repeat radar altimetry, J. Glaciol., 58, 830–840, https://doi.org/10.3189/2012JoG11J118, 2012.

Frezzotti, M., Gandolfi, S., and Urbini, S.: Snow megadunes in Antarctica: Sedimentary structure and genesis, J. Geophys. Res.-Atmos., 107, 4334, https://doi.org/10.1029/2001JD000673, 2002.

Frezzotti, M., Bitelli, G., De Michelis, P., Deponti, A., Forieri, A., Gandolfi, S., Maggi, V., Mancini, F., Remy, F., Tabacco, I. E., Urbini, S., Vittuari, L., and Zirizzottl, A.: Geophysical survey at Talos Dome, East Antarctica: the search for a new deep-drilling site, Ann. Glaciol., 39, 423–432, 2004.

Fung, A. K.: Microwave Scattering and Emission Models and their Applications, Artech house, 1994.

Helm, V., Humbert, A., and Miller, H.: Elevation and elevation change of Greenland and Antarctica derived from CryoSat-2, The Cryosphere, 8, 1539–1559, https://doi.org/10.5194/tc-8-1539-2014, 2014.

Inoue, J.: Surface drag over the snow surface of the Antarctic Plateau: 1. Factors controlling surface drag over the katabatic wind region, J. Geophys. Res.-Atmos., 94, 2207–2217, 1989.

Kotlyakov, V. M.: The Snow Cover of the Antarctic and its role in the Present-Day Glaciation of the Continent (Snezhni pokrov antarktidy i ego rol'v somvremennom oledenenii materika), Transl. Russ., 1966.

Lacroix, P., Legresy, B., Coleman, R., Dechambre, M., and Remy, F.: Dual-frequency altimeter signal from Envisat on the Amery ice-shelf, Remote Sens. Environ., 109, 285–294, 2007.

Lacroix, P., Legresy, B., Langley, K., Hamran, S. E., Kohler, J., Roques, S., Remy, F., and Dechambre, M.: In situ measurements of snow surface roughness using a laser profiler, J. Glaciol., 54, 753–762, 2008a.

Lacroix, P., Dechambre, M., Legresy, B., Blarel, F., and Remy, F.: On the use of the dual-frequency ENVISAT altimeter to determine snowpack properties of the Antarctic ice sheet, Remote Sens. Environ., 112, 1712–1729, 2008b.

Legresy, B. and Remy, F.: Using the temporal variability of satellite radar altimetric observations to map surface properties of the Antarctic ice sheet, J. Glaciol., 44, 197–206, 1998.

Legresy, B., Papa, F., Remy, F., Vinay, G., van den Bosch, M., and Zanife, O. Z.: ENVISAT radar altimeter measurements over continental surfaces and ice caps using the ICE-2 retracking algorithm, Remote Sens. Environ., 95, 150–163, 2005.

Lenaerts, J. T. M., Den Broeke, M. R., Berg, W. J., van Meijgaard, E., and Kuipers Munneke, P.: A new, high-resolution surface mass balance map of Antarctica (1979–2010) based on regional atmospheric climate modeling, Geophys. Res. Lett., 39, L04501, https://doi.org/10.1029/2011GL05071, 2012.

Li, J. and Zwally, H. J.: Modeled seasonal variations of firn density induced by steady-state surface air-temperature cycle, Ann. Glaciol., 34, 299–302, 2002.

Li, J. and Zwally, H. J.: Modeling the density variation in the shallow firn layer, Ann. Glaciol., 38, 309–313, 2004.

Li, J. and Zwally, H. J.: Modeling of firn compaction for estimating ice-sheet mass change from observed ice-sheet elevation change, Ann. Glaciol., 52, 1–7, 2011.

Male, D. H.: 6 – THE SEASONAL SNOWCOVER, in: Dynamics of Snow and Ice Masses, edited by: Colbeck, S. C., 305–395, Academic Press., 1980.

Mätzler, C.: Improved Born approximation for scattering of radiation in a granular medium, J. Appl. Phys., 83, 6111–6117, 1998.

Michel, A., Flament, T., and Remy, F.: Study of the Penetration Bias of ENVISAT Altimeter Observations over Antarctica in Comparison to ICESat Observations, Remote Sens., 6, 9412–9434, https://doi.org/10.3390/rs6109412, 2014.

Nilsson, J., Gardner, A., Sandberg Sørensen, L., and Forsberg, R.: Improved retrieval of land ice topography from CryoSat-2 data and its impact for volume-change estimation of the Greenland Ice Sheet, The Cryosphere, 10, 2953–2969, https://doi.org/10.5194/tc-10-2953-2016, 2016.

Remy, F. and Minster, J. F.: A Comparison between Active and Passive Microwave Measurements of the Antarctic Ice-Sheet and Their Association with the Surface Katabatic Winds, J. Glaciol., 37, 3–10, 1991.

Remy, F. and Parouty, S.: Antarctic ice sheet and radar altimetry: A review, Remote Sens., 1, 1212–1239, https://doi.org/10.3390/rs1041212, 2009.

Remy, F., Shaeffer, P., and Legresy, B.: Ice flow physical processes derived from the ERS-1 high-resolution map of the Antarctica and Greenland ice sheets, Geophys. J. Int., 139, 645–656, 1999.

Remy, F., Flament, T., Blarel, F., and Benveniste, J.: Radar altimetry measurements over antarctic ice sheet: A focus on antenna polarization and change in backscatter problems, Adv. Space Res., 50, 998–1006, 2012.

Remy, F., Flament, T., Michel, A., and Blumstein, D.: Envisat and SARAL/AltiKa Observations of the Antarctic Ice Sheet: A Comparison Between the Ku-band and Ka-band, Mar. Geod., 38, 510–521, https://doi.org/10.1080/01490419.2014.985347, 2015.

Ridley, J. and Partington, K.: A Model of Satellite Radar Altimeter Return from Ice Sheets, Int. J. Remote Sens., 9, 601–624, 1988.

Scambos, T. A., Frezzotti, M., Haran, T., Bohlander, J., Lenaerts, J. T. M., Van Den Broeke, M. R., Jezek, K., Long, D., Urbini, S., Farness, K., Neumann, T., Winther, J.-G., and Albert, M.: Extent of low-accumulation'wind glaze'areas on the East Antarctic plateau: implications for continental ice mass balance, J. Glaciol., 58, 633–647, https://doi.org/10.3189/2012JoG11J232, 2012.

Steunou, N., Desjonqueres, J. D., Picot, N., Sengenes, P., Noubel, J., and Poisson, J. C.: AltiKa Altimeter: Instrument Description and In Flight Performance, Mar. Geod., 38, 22–42, https://doi.org/10.1080/01490419.2014.988835, 2015.

Sturm, M., Holmgren, J., König, M., and Morris, K.: The thermal conductivity of seasonal snow, J. Glaciol., 43, 26–41, 1997.

Surdyk, S.: Using microwave brightness temperature to detect short-term surface air temperature changes in Antarctica: An analytical approach, Remote Sens. Environ., 80, 256–271, 2002.

Tiuri, M., Sihvola, A., Nyfors, E., and Hallikaiken, M.: The complex dielectric constant of snow at microwave frequencies, IEEE J. Ocean. Eng., 9, 377–382, https://doi.org/10.1109/JOE.1984.1145645, 1984.

Ulaby, F. T., Moore, R. K., and Fung, A. K.: Microwave Remote Sensing Active and Passive-Volume II: Radar Remote Sensing and Surface Scattering and Emission Theory, Artech house, Boston, London, available at: https://infoscience.epfl.ch/record/51982 (last access: 23 May 2018), 1982.

Wingham, D. J., Ridout, A. J., Scharroo, R., Arthern, R. J., and Shum, C. K.: Antarctic Elevation Change from 1992 to 1996, Science, 282, 456–458, https://doi.org/10.1126/science.282.5388.456, 1998.

Wingham, D. J., Shepherd, A., Muir, A., and Marshall, G. J.: Mass balance of the Antarctic ice sheet, Philos. Trans. R. Soc. Lond. Math. Phys. Eng. Sci., 364, 1627–1635, 2006.

Zwally, H. J., Giovinetto, M. B., Li, J., Cornejo, H. G., Beckley, M. A., Brenner, A. C., Saba, J. L., and Yi, D.: Mass changes of the Greenland and Antarctic ice sheets and shelves and contributions to sea-level rise: 1992–2002, J. Glaciol., 51, 509–527, 2005.