Penetration of interferometric radar signals in Antarctic snow

Abstract

Abstract. Synthetic aperture radar interferometry (InSAR) is an efficient technique for mapping the surface elevation and its temporal change over glaciers and ice sheets. However, due to the penetration of the SAR signal into snow and ice, the apparent elevation in uncorrected InSAR digital elevation models (DEMs) is displaced versus the actual surface. We studied relations between interferometric radar signals and physical snow properties and tested procedures for correcting the elevation bias. The work is based on satellite and in situ data over Union Glacier in the Ellsworth Mountains, West Antarctica, including interferometric data of the TanDEM-X mission, topographic data from optical satellite sensors and field measurements on snow structure, and stratigraphy undertaken in December 2016. The study area comprises ice-free surfaces, bare ice, dry snow and firn with a variety of structural features related to local differences in wind exposure and snow accumulation. Time series of laser measurements of NASA's Ice, Cloud and land Elevation Satellite (ICESat) and ICESat-2 show steady-state surface topography. For area-wide elevation reference we use the Reference Elevation Model of Antarctica (REMA). The different elevation data are vertically co-registered on a blue ice area that is not affected by radar signal penetration. Backscatter simulations with a multilayer radiative transfer model show large variations for scattering of individual snow layers, but the vertical backscatter distribution can be approximated by an exponential function representing uniform absorption and scattering properties. We obtain estimates of the elevation bias by inverting the interferometric volume correlation coefficient (coherence), applying a uniform volume model for describing the vertical loss function. Whereas the mean values of the computed elevation bias and the elevation difference between the TanDEM-X DEMs and the REMA show good agreement, a trend towards overestimation of penetration is evident for heavily wind-exposed areas with low accumulation and towards underestimation for areas with higher accumulation rates. In both cases deviations from the uniform volume structure are the main reason. In the first case the dense sequence of horizontal structures related to internal wind crust, ice layers and density stratification causes increased scattering in near-surface layers. In the second case the small grain size of the top snow layers causes a downward shift in the scattering phase centre.