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

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

Spaceborne synthetic aperture radar: Current status and future directions. A report to the Committee on Earth Sciences, Space Studies Board, National Research Council

This report provides a context in which questions put forth by NASA's Office of Mission to Planet Earth (OMPTE) regarding the next steps in spaceborne synthetic aperture radar (SAR) science and technology can be addressed. It summarizes the state-of-the-art in theory, experimental design, technology, data analysis, and utilization of SAR data for studies of the Earth, and describes potential new applications. The report is divided into five science chapters and a technology assessment. The chapters summarize the value of existing SAR data and currently planned SAR systems, and identify gaps in observational capabilities needing to be filled to address the scientific questions. Cases where SAR provides complementary data to other (non-SAR) measurement techniques are also described. The chapter on technology assessment outlines SAR technology development which is critical not only to NASA's providing societally relevant geophysical parameters but to maintaining competitiveness in SAR technology, and promoting economic development

Monitoring Snow Cover and Snowmelt Dynamics and Assessing their Influences on Inland Water Resources

Snow is one of the most vital cryospheric components owing to its wide coverage as well as its unique physical characteristics. It not only affects the balance of numerous natural systems but also influences various socio-economic activities of human beings. Notably, the importance of snowmelt water to global water resources is outstanding, as millions of populations rely on snowmelt water for daily consumption and agricultural use. Nevertheless, due to the unprecedented temperature rise resulting from the deterioration of climate change, global snow cover extent (SCE) has been shrinking significantly, which endangers the sustainability and availability of inland water resources. Therefore, in order to understand cryo-hydrosphere interactions under a warming climate, (1) monitoring SCE dynamics and snowmelt conditions, (2) tracking the dynamics of snowmelt-influenced waterbodies, and (3) assessing the causal effect of snowmelt conditions on inland water resources are indispensable. However, for each point, there exist many research questions that need to be answered. Consequently, in this thesis, five objectives are proposed accordingly. Objective 1: Reviewing the characteristics of SAR and its interactions with snow, and exploring the trends, difficulties, and opportunities of existing SAR-based SCE mapping studies; Objective 2: Proposing a novel total and wet SCE mapping strategy based on freely accessible SAR imagery with all land cover classes applicability and global transferability; Objective 3: Enhancing total SCE mapping accuracy by fusing SAR- and multi-spectral sensor-based information, and providing total SCE mapping reliability map information; Objective 4: Proposing a cloud-free and illumination-independent inland waterbody dynamics tracking strategy using freely accessible datasets and services; Objective 5: Assessing the influence of snowmelt conditions on inland water resources

SCAR Report on Antarctic Climate Change and the Environment

The first comprehensive review of the state of Antarctica’s climate and its relationship to the global climate by the Scientific Committee on Antarctic Research (SCAR). The review - Antarctic Climate Change and the Environment – presents the latest research from the icy continent, identifies areas for future scientific research, and addresses the urgent questions that policy makers have about Antarctic melting, sea-level rise and biodiversity

Modeling Multifrequency Pol-InSAR Data From the Percolation Zone of the Greenland Ice Sheet

The analysis of data from an airborne synthetic aperture radar (SAR) campaign in the percolation zone of Greenland revealed an interferometric coherence undulation behavior with respect to vertical wavenumber, which cannot be explained with existing models. We propose a model extension that accounts for scattering from distinct layers below the surface. Simulations show that the periodicity of the coherence undulation is mainly driven by the vertical distance between dominant subsurface layers, while the amplitude of the undulation is determined by the ratio between scattering from distinct layers and scattering from the firn volume. We use the model to interpret quad-pol SAR data at X-, C-, S-, L-, and P-bands. The inferred layer depths match layer detections in ground-based radar data and in situ measurements. We conclude that in the percolation zone, scattering from subsurface layers has to be taken into account to correctly interpret SAR data and demonstrate the potential to retrieve geophysical information about the vertical subsurface structure

Towards a Pol-InSAR Firn Density Retrieval

Polarimetric and (multi-baseline) interferometric SAR techniques are promising tools to investigate the subsurface properties of glaciers and ice sheets, due to the signal penetration of up to several tens of meters into dry snow, firn, and ice. Two different lines of research were addressed in recent years. The first is based on PolSAR, which provides not only information about the scattering mechanisms, but also has the uniqueness of being sensitive to anisotropic signal propagation in non-scattering layers of snow and firn. The second line of research is related to the use of Pol-InSAR and TomoSAR to retrieve the 3D location of scatterers within the subsurface. So far, the different SAR techniques were mainly assessed separately. In the field of PolSAR modeling efforts have been dedicated to establish a link between the co-polarization HH-VV phase difference (CPD) and the structural properties of firn [1]. CPDs have then been interpreted as the result of propagation effects due to the dielectric anisotropy of the firn volume. This modeling approach establishes a relationship between the measured CPD and firn density, firn anisotropy and the vertical backscattering distribution in the subsurface of the glacier or ice sheet. By assuming bulk values for density and anisotropy and employing a constant signal extinction for the vertical backscattering function, i.e. a uniform volume, a first attempt to retrieve firn thickness from PolSAR data was presented [1]. In the fields of Pol-InSAR and TomoSAR for the investigation of the subsurface scattering structure of glaciers and ice sheets, recent studies were concerned with the estimation of the vertical backscatter distribution, either model-based or through tomographic imaging techniques. The complexity of (Pol-)InSAR models for the retrieval of subsurface structure information is mainly limited by the available observation space. Thus, constant signal extinction volumes [2], with additional Dirac deltas to represent refrozen melt layers [3] and variable extinction volumes [4] have been modelled and used to retrieve information about the subsurface. With TomoSAR, the imaging of subsurface features in glaciers [5], and ice sheets [4][6][7] was demonstrated and the effect of subsurface layers, different ice types, firn bodies, crevasses, and bed rock was recognized in the tomograms. Those studies on PolSAR, Pol-InSAR and TomoSAR techniques made promising steps towards a subsurface structure information retrieval on glaciers and ice sheets, but the direct relationship to geophysical parameters, which are useful for the glaciological community, is only limited. The most promising way to achieve this goal is the combination of PolSAR and Pol-InSAR/TomoSAR techniques in order to exploit the synergies between the individual methods and by integrating the models and algorithms into one common framework. A first experiment in this direction was the integration of TomoSAR vertical scattering profiles into the PolSAR firn anisotropy model [8]. This allowed the inversion of the firn density from experimental airborne F-SAR data over Greenland, which would be a promising geophysical information product. However, this experiment enabled only the inversion of a single bulk density value for the entire penetration depth and was based on hundreds of samples of CPD measurements and TomoSAR profiles across a wide range of incidence angles, which is of limited applicability. Therefore, this study will continue this investigation, assessing to which degree an inversion is possible on reduced and more feasible observation spaces. Open questions are, which number and range of incidence angles as well as which and how many baselines are required, and whether depth-varying densities could be retrieved. The results should give an indication if such an approach might be feasible with future spaceborne SAR systems. [1] G. Parrella, I. Hajnsek, and K. P. Papathanassiou, “Retrieval of Firn Thickness by Means of Polarisation Phase Differences in L-Band SAR Data,” Remote Sensing, vol. 13, no. 21, p. 4448, Nov. 2021, doi: 10.3390/rs13214448. [2] E. W. Hoen and H. Zebker, “Penetration depths inferred from interferometric volume decorrelation observed over the Greenland ice sheet,” IEEE Transactions on Geoscience and Remote Sensing, vol. 38, no. 6, pp. 2571–2583, 2000. [3] G. Fischer, K. P. Papathanassiou and I. Hajnsek, "Modeling Multifrequency Pol-InSAR Data from the Percolation Zone of the Greenland Ice Sheet," IEEE Transactions on Geoscience and Remote Sensing, vol. 57, no. 4, pp. 1963-1976, 2019. [4] G. Fischer, M. Jäger, K. P. Papathanassiou and I. Hajnsek, "Modeling the Vertical Backscattering Distribution in the Percolation Zone of the Greenland Ice Sheet with SAR Tomography," IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, vol. 12, no. 11, pp. 4389-4405, 2019. [5] S. Tebaldini, T. Nagler, H. Rott, and A. Heilig, “Imaging the Internal Structure of an Alpine Glacier via L-Band Airborne SAR Tomography,” IEEE Transactions on Geoscience and Remote Sensing, vol. 54, no. 12, pp. 7197–7209, 2016. [6] F. Banda, J. Dall, and S. Tebaldini, “Single and Multipolarimetric P-Band SAR Tomography of Subsurface Ice Structure,” IEEE Transactions on Geoscience and Remote Sensing, vol. 54, no. 5, pp. 2832–2845, 2016. [7] M. Pardini, G. Parrella, G. Fischer, and K. Papathanassiou, “A Multi-Frequency SAR Tomographic Characterization of Sub-Surface Ice Volumes,” in Proceedings of EUSAR, Hamburg, Germany, 2016. [8] G. Fischer, K. Papathanassiou, I. Hajnsek, and G. Parrella, “Combining PolSAR, Pol-InSAR and TomoSAR for Snow and Ice Subsurface Characterization,” in Proceedings of the ESA POLinSAR Workshop, Online, Apr. 2021