Modeling of Subsurface Scattering from Ice Sheets for Pol-InSAR Applications

Remote sensing is a fundamental tool to measure the dynamics of ice sheets and provides valuable information for ice sheet projections under a changing climate. There is, however, the potential to further reduce the uncertainties in these projections by developing innovative remote sensing methods. One of these remote sensing techniques, the polarimetric synthetic aperture radar interferometry (Pol-InSAR), is known since decades to have the potential to assess the geophysical properties below the surface of ice sheets, because of the penetration of microwave signals into dry snow, firn, and ice. Despite this, only very few studies have addressed this topic and the development of robust Pol-InSAR applications is at an early stage. Two potential Pol-InSAR applications are identified as the motivation for this thesis. First, the estimation and compensation of the penetration bias in digital elevation models derived with SAR interferometry. This bias can lead to errors of several meters or even tens of meters in surface elevation measurements. Second, the estimation of geophysical properties of the subsurface of glaciers and ice sheets using Pol-InSAR techniques. There is indeed potential to derive information about melt-refreeze processes within the firn, which are related to density and affect the mass balance. Such Pol-InSAR applications can be a valuable information source with the potential for monthly ice sheet wide coverage and high spatial resolution provided by the next generation of SAR satellites. However, the required models to link the Pol-InSAR measurements to the subsurface properties are not yet established. The aim of this thesis is to improve the modeling of the vertical backscattering distribution in the subsurface of ice sheets and its effect on polarimetric interferometric SAR measurements at different frequencies. In order to achieve this, polarimetric interferometric multi-baseline SAR data at different frequencies and from two different test sites on the Greenland ice sheet are investigated. This thesis contributes with three concepts to a better understanding and to a more accurate modeling of the vertical backscattering distribution in the subsurface of ice sheets. First, the integration of scattering from distinct subsurface layers. These are formed by refrozen melt water in the upper percolation zone and cause an interesting coherence undulation pattern, which cannot be explained with previously existing models. This represents a first link between Pol-InSAR data and geophysical subsurface properties. The second step is the improved modeling of the general vertical backscattering distribution of the subsurface volume. The advantages of more flexible volume models are demonstrated, but interestingly, the simple modification of a previously existing model with a vertical shift parameter lead to the best agreement between model and data. The third contribution is the model based compensation of the penetration bias, which is experimentally validated. At the investigated test sites, it becomes evident that the model based estimates of the surface elevations are more accurate than the interferometric phase center locations, which are conventionally used to derive surface elevations of ice sheets. This thesis therefore improves the state of the art of subsurface scattering modeling for Pol-InSAR applications, demonstrates the model-based penetration bias compensation, and makes a further research step towards the retrieval of geophysical subsurface information with Pol-InSAR

Cryosphere Applications

Synthetic aperture radar (SAR) provides large coverage and high resolution, and it has been proven to be sensitive to both surface and near-surface features related to accumulation, ablation, and metamorphism of snow and firn. Exploiting this sensitivity, SAR polarimetry and polarimetric interferometry found application to land ice for instance for the estimation of wave extinction (which relates to sub surface ice volume structure) and for the estimation of snow water equivalent (which relates to snow density and depth). After presenting these applications, the Chapter proceeds by reviewing applications of SAR polarimetry to sea ice for the classification of different ice types, the estimation of thickness, and the characterisation of its surface. Finally, an application to the characterisation of permafrost regions is considered. For each application, the used (model-based) decomposition and polarimetric parameters are critically described, and real data results from relevant airborne campaigns and space borne acquisitions are reported

Investigation of Sea Ice Using Multiple Synthetic Aperture Radar Acquisitions

The papers of this thesis are not available in Munin. Paper I: Yitayew, T. G., Ferro-Famil, L., Eltoft, T. & Tebaldini, S. (2017). Tomographic imaging of fjord ice using a very high resolution ground-based SAR system. Available in IEEE Transactions on Geoscience and Remote Sensing, 55 (2):698-714. Paper II: Yitayew, T. G., Ferro-Famil, L., Eltoft, T. & Tebaldini, S. (2017). Lake and fjord ice imaging using a multifrequency ground-based tomographic SAR system. Available in IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 10(10):4457-4468. Paper III: Yitayew, T. G., Divine, D. V., Dierking, W., Eltoft, T., Ferro-Famil, L., Rosel, A. & Negrel, J. Validation of Sea ice Topographic Heights Derived from TanDEMX Interferometric SAR Data with Results from Laser Profiler and Photogrammetry. (Manuscript).The thesis investigates imaging in the vertical direction of different types of ice in the arctic using synthetic aperture radar (SAR) tomography and SAR interferometry. In the first part, the magnitude and the positions of the dominant scattering contributions within snow covered fjord and lake ice layers are effectively identified by using a very high resolution ground-based tomographic SAR system. Datasets collected at multiple frequencies and polarizations over two test sites in Tromsø area, northern Norway, are used for characterizing the three-dimensional response of snow and ice. The presented experimental results helped to improve our understanding of the interaction between radar waves and snow and ice layers. The reconstructed radar responses are also used for estimating the refractive indices and the vertical positions of the different sub-layers of snow and ice. The second part of the thesis deals with the retrieval of the surface topography of multi-year sea ice using SAR interferometry. Satellite acquisitions from TanDEM-X over the Svalbard area are used for analysis. The retrieved surface height is validated by using overlapping helicopter-based stereo camera and laser profiler measurements, and a very good agreement has been found. The work contributes to an improved understanding regarding the potential of SAR tomography for imaging the vertical scattering distribution of snow and ice layers, and for studying the influence of both sensor parameters such as its frequency and polarization and scene properties such as layer stratification, air bubbles and small-scale roughness of the interfaces on snow and ice backscattered signal. Moreover, the presented results reveal the potential of SAR interferometry for retrieving the surface topography of sea ice

Multi-Baseline Pol-InSAR Inversion of the Subsurface Scattering Structure of Ice Sheets

The influence of the subsurface properties of ice sheets on polarimetric synthetic aperture radar interferometry (Pol-InSAR) measurements is well known. In order to invert this relationship for the extraction of geophysical parameters from Pol-InSAR data, models of the subsurface scattering structure are required. One potential application is the estimation of the penetration bias in interferometric surface elevation measurements of ice sheets, which was demonstrated based on single-baseline data. However, the model complexity and performance are constrained by the limited observation space. This study, therefore, investigates the inversion of subsurface scattering structures with multi-baseline fully polarimetric Pol-InSAR data, which allows accounting for more realistic scattering scenarios. Preliminary results indicate a more robust inversion of the penetration bias compared to the single-baseline case

The European Space Agency BIOMASS mission: Measuring forest above-ground biomass from space

The primary objective of the European Space Agency's 7th Earth Explorer mission, BIOMASS, is to determine the worldwide distribution of forest above-ground biomass (AGB) in order to reduce the major uncertainties in calculations of carbon stocks and fluxes associated with the terrestrial biosphere, including carbon fluxes associated with Land Use Change, forest degradation and forest regrowth. To meet this objective it will carry, for the first time in space, a fully polarimetric P-band synthetic aperture radar (SAR). Three main products will be provided: global maps of both AGB and forest height, with a spatial resolution of 200 m, and maps of severe forest disturbance at 50 m resolution (where “global” is to be understood as subject to Space Object tracking radar restrictions). After launch in 2022, there will be a 3-month commissioning phase, followed by a 14-month phase during which there will be global coverage by SAR tomography. In the succeeding interferometric phase, global polarimetric interferometry Pol-InSAR coverage will be achieved every 7 months up to the end of the 5-year mission. Both Pol-InSAR and TomoSAR will be used to eliminate scattering from the ground (both direct and double bounce backscatter) in forests. In dense tropical forests AGB can then be estimated from the remaining volume scattering using non-linear inversion of a backscattering model. Airborne campaigns in the tropics also indicate that AGB is highly correlated with the backscatter from around 30 m above the ground, as measured by tomography. In contrast, double bounce scattering appears to carry important information about the AGB of boreal forests, so ground cancellation may not be appropriate and the best approach for such forests remains to be finalized. Several methods to exploit these new data in carbon cycle calculations have already been demonstrated. In addition, major mutual gains will be made by combining BIOMASS data with data from other missions that will measure forest biomass, structure, height and change, including the NASA Global Ecosystem Dynamics Investigation lidar deployed on the International Space Station after its launch in December 2018, and the NASA-ISRO NISAR L- and S-band SAR, due for launch in 2022. More generally, space-based measurements of biomass are a core component of a carbon cycle observation and modelling strategy developed by the Group on Earth Observations. Secondary objectives of the mission include imaging of sub-surface geological structures in arid environments, generation of a true Digital Terrain Model without biases caused by forest cover, and measurement of glacier and icesheet velocities. In addition, the operations needed for ionospheric correction of the data will allow very sensitive estimates of ionospheric Total Electron Content and its changes along the dawn-dusk orbit of the mission

Polarimetric Synthetic Aperture Radar

This open access book focuses on the practical application of electromagnetic polarimetry principles in Earth remote sensing with an educational purpose. In the last decade, the operations from fully polarimetric synthetic aperture radar such as the Japanese ALOS/PalSAR, the Canadian Radarsat-2 and the German TerraSAR-X and their easy data access for scientific use have developed further the research and data applications at L,C and X band. As a consequence, the wider distribution of polarimetric data sets across the remote sensing community boosted activity and development in polarimetric SAR applications, also in view of future missions. Numerous experiments with real data from spaceborne platforms are shown, with the aim of giving an up-to-date and complete treatment of the unique benefits of fully polarimetric synthetic aperture radar data in five different domains: forest, agriculture, cryosphere, urban and oceans

Polarimetric Synthetic Aperture Radar, Principles and Application

Demonstrates the benefits of the usage of fully polarimetric synthetic aperture radar data in applications of Earth remote sensing, with educational and development purposes. Includes numerous up-to-date examples with real data from spaceborne platforms and possibility to use a software to support lecture practicals. Reviews theoretical principles in an intuitive way for each application topic. Covers in depth five application domains (forests, agriculture, cryosphere, urban, and oceans), with reference also to hazard monitorin

Processing of optic and radar images.Application in satellite remote sensing of snow, ice and glaciers

Ce document présente une synthèse de mes activités de recherche depuis la soutenance de ma thèse en 1999. L'activité rapportée ici est celle d'un ingénieur de recherche, et donc s'est déroulée en parallèle d'une activité ``technique'' comprenant des taches d'instrumentation en laboratoire, d'instrumentation de plateformes en montagne, de raids scientifiques sur les calottes polaires, d'élaboration de projets scientifiques, d'organisation d'équipes ou d'ordre administratif. Je suis Ingénieur de recherche CNRS depuis 2004 affecté au laboratoire Gipsa-lab, une unité mixte de recherche du CNRS, de Grenoble-INP, de l'université Joseph Fourier et de l'université Stendhal. Ce laboratoire (d'environ 400 personnes), conventionné avec l'INRIA, l'Observatoire de Grenoble et l'université Pierre Mendès France, est pluridisciplinaire et développe des recherches fondamentales et finalisées sur les signaux et les systèmes complexes.}Lors de la préparation de ma thèse (mi-temps 1995-99) au LGGE, je me suis intéressé au traitement des images de microstructures de la neige, du névé et de la glace. C'est assez naturellement que j'ai rejoint le laboratoire LIS devenu Gipsa-lab pour y développer des activités de traitement des images Radar à Synthèse d'Ouverture (RSO) appliqué aux milieux naturels neige, glace et glaciers. Etant le premier à générer un interférogramme différentiel des glaciers des Alpes, j'ai continué à travailler sur la phase interférométrique pour extraire des informations de déplacement et valider ces méthodes sur le glacier d'Argentière (massif du Mont-Blanc) qui présente l'énorme avantage de se déplacer de quelques centimètres par jour. Ces activités m'ont amené à développer, en collaboration avec les laboratoires LISTIC, LTCI et IETR, des méthodes plus générales pour extraire des informations dans les images RSO.Ma formation initiale en électronique, puis de doctorat en physique m'ont amené à mettre à profit mes connaissances en traitement d'images et des signaux, en électromagnétisme, en calcul numérique, en informatique et en physique de la neige et de la glace pour étudier les problèmes de traitement des images RSO appliqués à la glace, aux glaciers et à la neige