Ionospheric correction of interferometric SAR data with application to the cryospheric sciences

Thesis (Ph.D.) University of Alaska Fairbanks, 2018The ionosphere has been identified as an important error source for spaceborne Synthetic Aperture Radar (SAR) data and SAR Interferometry (InSAR), especially for low frequency SAR missions, operating, e.g., at L-band or P-band. Developing effective algorithms for the correction of ionospheric effects is still a developing and active topic of remote sensing research. The focus of this thesis is to develop robust and accurate techniques for ionospheric correction of SAR and InSAR data and evaluate the benefit of these techniques for cryospheric research fields such as glacier ice velocity tracking and permafrost deformation monitoring. As both topics are mostly concerned with high latitude areas where the ionosphere is often active and characterized by turbulence, ionospheric correction is particularly relevant for these applications. After an introduction to the research topic in Chapter 1, Chapter 2 will discuss open issues in ionospheric correction including processing issues related to baseline-induced spectrum shifts. The effect of large baseline on split spectrum InSAR technique has been thoroughly evaluated and effective solutions for compensating this effect are proposed. In addition, a multiple sub-band approach is proposed for increasing the algorithm robustness and accuracy. Selected case studies are shown with the purpose of demonstrating the performance of the developed algorithm. In Chapter 3, the developed ionospheric correction technology is applied to optimize InSAR-based ice velocity measurements over the big ice sheets in Greenland and the Antarctic. Selected case studies are presented to demonstrate and validate the effectiveness of the proposed correction algorithms for ice velocity applications. It is shown that the ionosphere signal can be larger than the actual glacier motion signal in the interior of Greenland and Antarctic, emphasizing the necessity for operational ionospheric correction. The case studies also show that the accuracy of ice velocity estimates was significantly improved once the developed ionospheric correction techniques were integrated into the data processing flow. We demonstrate that the proposed ionosphere correction outperforms the traditionally-used approaches such as the averaging of multi-temporal data and the removal of obviously affected data sets. For instance, it is shown that about one hundred multi-temporal ice velocity estimates would need to be averaged to achieve the estimation accuracy of a single ionosphere-corrected measurement. In Chapter 4, we evaluate the necessity and benefit of ionospheric-correction for L-band InSAR-based permafrost research. In permafrost zones, InSAR-based surface deformation measurements are used together with geophysical models to estimate permafrost parameters such as active layer thickness, soil ice content, and permafrost degradation. Accurate error correction is needed to avoid biases in the estimated parameters and their co-variance properties. Through statistical analyses of a large number of L-band InSAR data sets over Alaska, we show that ionospheric signal distortions, at different levels of magnitude, are present in almost every InSAR dataset acquired in permafrost-affected regions. We analyze the ionospheric correction performance that can be achieved in permafrost zones by statistically analyzing correction results for large number of InSAR data. We also investigate the impact of ionospheric correction on the performance of the two main InSAR approaches that are used in permafrost zones: (1) we show the importance of ionospheric correction for permafrost deformation estimation from discrete InSAR observations; (2) we demonstrate that ionospheric correction leads to significant improvements in the accuracy of time-series InSAR-based permafrost products. Chapter 5 summarizes the work conducted in this dissertation and proposes next steps in this field of research

Advancing the Monitoring Capabilities of Mountain Snowpack Fluctuations at Various Spatial and Temporal Scales

Snow is a critical water resource for the western US and many regions across the globe. However, our ability to accurately monitor changes in snow mass from satellite remote sensing, specifically its water equivalent, remains a challenge in mountain regions. No single sensor currently has the ability to directly measure snow water equivalent (SWE) from space at a spatial scale suitable for water supply forecasting in mountain environments. This knowledge gap calls for the innovative use of remote sensing techniques, computational tools, and data science methods to advance our ability to estimate mountain snowpacks across a range of spatial and temporal scales. The goal of this dissertation is to advance our capabilities for understanding snowpack across watershed-relevant spatial and temporal scales. Two research approaches were used to accomplish this goal: quantifying the physiographic controls and sensitivities of hydrologically important snow metrics and progressing our ability to use L-band interferometric synthetic aperture radar (InSAR) to measure SWE changes. First, we quantify the physiographic controls and various snowpack metrics in the Sierra Nevada using a novel gridded SWE reanalysis dataset. Such work demonstrates the complexity of snowpack processes and the need for fine-resolution snowpack information. Next, using L-band Interferometric Synthetic Aperture Radar (InSAR) from the NASA SnowEx campaign, both snow ablation and accumulation are estimated in the Jemez Mountains, NM. The radar-derived retrievals are evaluated utilizing a combination of optical snow-cover data, snow pits, meteorological station data, in situ snow depth sensors, and ground-penetrating radar (GPR). Lastly, we compare multisensor optical-radar approaches for SWE retrievals and find that moderate-resolution legacy satellite products provide sufficient results. The results of this work show that L-band InSAR is a suitable technique for global SWE monitoring when used synergistically with optical SCA data and snowpack modeling. While two distinctive methods are present in this research, they both work towards advancing our ability to understand the dynamics of mountain snowpack

Autonomous Extraction of Millimeter-scale Deformation in InSAR Time Series Using Deep Learning

Systematic characterization of slip behaviours on active faults is key to unraveling the physics of tectonic faulting and the interplay between slow and fast earthquakes. Interferometric Synthetic Aperture Radar (InSAR), by enabling measurement of ground deformation at a global scale every few days, may hold the key to those interactions. However, atmospheric propagation delays often exceed ground deformation of interest despite state-of-the art processing, and thus InSAR analysis requires expert interpretation and a priori knowledge of fault systems, precluding global investigations of deformation dynamics. Here we show that a deep auto-encoder architecture tailored to untangle ground deformation from noise in InSAR time series autonomously extracts deformation signals, without prior knowledge of a fault's location or slip behaviour. Applied to InSAR data over the North Anatolian Fault, our method reaches 2 mm detection, revealing a slow earthquake twice as extensive as previously recognized. We further explore the generalization of our approach to inflation/deflation-induced deformation, applying the same methodology to the geothermal field of Coso, California

Statistical comparison of InSAR tropospheric correction techniques

Correcting for tropospheric delays is one of the largest challenges facing the interferometric synthetic aperture radar (InSAR) community. Spatial and temporal variations in temperature, pressure, and relative humidity create tropospheric signals in InSAR data, masking smaller surface displacements due to tectonic or volcanic deformation. Correction methods using weather model data, GNSS and/or spectrometer data have been applied in the past, but are often limited by the spatial and temporal resolution of the auxiliary data. Alternatively a correction can be estimated from the interferometric phase by assuming a linear or a power-law relationship between the phase and topography. Typically the challenge lies in separating deformation from tropospheric phase signals. In this study we performed a statistical comparison of the state-of-the-art tropospheric corrections estimated from the MERIS and MODIS spectrometers, a low and high spatial-resolution weather model (ERA-I and WRF), and both the conventional linear and new power-law empirical methods. Our test-regions include Southern Mexico, Italy, and El Hierro. We find spectrometers give the largest reduction in tropospheric signal, but are limited to cloud-free and daylight acquisitions. We find a ~ 10–20% RMSE increase with increasing cloud cover consistent across methods. None of the other tropospheric correction methods consistently reduced tropospheric signals over different regions and times. We have released a new software package called TRAIN (Toolbox for Reducing Atmospheric InSAR Noise), which includes all these state-of-the-art correction methods. We recommend future developments should aim towards combining the different correction methods in an optimal manner

Estimating Snow Accumulation and Ablation with L-Band Interferometric Synthetic Aperture Radar (InSAR)

Snow is a critical water resource for the western United States and many regions across the globe. However, our ability to accurately measure and monitor changes in snow mass from satellite remote sensing, specifically its water equivalent, remains a challenge. To confront these challenges, NASA initiated the SnowEx program, a multiyear effort to address knowledge gaps in snow remote sensing. During SnowEx 2020, the Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) team acquired an L-band interferometric synthetic aperture radar (InSAR) data time series to evaluate the capabilities and limitations of repeat-pass L-band InSAR for tracking changes in snow water equivalent (SWE). The goal was to develop a more comprehensive understanding of where and when L-band InSAR can provide SWE change estimates, allowing the snow community to leverage the upcoming NASA–ISRO (NASA–Indian Space Research Organization) SAR (NISAR) mission. Our study analyzed three InSAR image pairs from the Jemez Mountains, NM, between 12 and 26 February 2020. We developed a snow-focused multi-sensor method that uses UAVSAR InSAR data synergistically with optical fractional snow-covered area (fSCA) information. Combining these two remote sensing datasets allows for atmospheric correction and delineation of snow-covered pixels within the radar swath. For all InSAR pairs, we converted phase change values to SWE change estimates between the three acquisition dates. We then evaluated InSAR-derived retrievals using a combination of fSCA, snow pits, meteorological station data, in situ snow depth sensors, and ground-penetrating radar (GPR). The results of this study show that repeat-pass L-band InSAR is effective for estimating both snow accumulation and ablation with the proper measurement timing, reference phase, and snowpack conditions

Characterizing slope instability kinematics by integrating multi-sensor satellite remote sensing observations

Over the past few decades, the occurrence and intensity of geological hazards, such as landslides, have substantially risen due to various factors, including global climate change, seismic events, rapid urbanization and other anthropogenic activities. Landslide disasters pose a significant risk in both urban and rural areas, resulting in fatalities, infrastructure damages, and economic losses. Nevertheless, conventional ground-based monitoring techniques are often costly, time-consuming, and require considerable resources. Moreover, some landslide incidents occur in remote or hazardous locations, making ground-based observation and field investigation challenging or even impossible. Fortunately, the advancements in spaceborne remote sensing technology have led to the availability of large-scale and high-quality imagery, which can be utilized for various landslide-related applications, including identification, monitoring, analysis, and prediction. This efficient and cost-effective technology allows for remote monitoring and assessment of landslide risks and can significantly contribute to disaster management and mitigation efforts. Consequently, spaceborne remote sensing techniques have become vital for geohazard management in many countries, benefiting society by providing reliable downstream services. However, substantial effort is required to ensure that such benefits are provided. For establishing long-term data archives and reliable analyses, it is essential to maintain consistent and continued use of multi-sensor spaceborne remote sensing techniques. This will enable a more thorough understanding of the physical mechanisms responsible for slope instabilities, leading to better decision-making and development of effective mitigation strategies. Ultimately, this can reduce the impact of landslide hazards on the general public. The present dissertation contributes to this effort from the following perspectives: 1. To obtain a comprehensive understanding of spaceborne remote sensing techniques for landslide monitoring, we integrated multi-sensor methods to monitor the entire life cycle of landslide dynamics. We aimed to comprehend the landslide evolution under complex cascading events by utilizing various spaceborne remote sensing techniques, e.g., the precursory deformation before catastrophic failure, co-failure procedures, and post-failure evolution of slope instability. 2. To address the discrepancies between spaceborne optical and radar imagery, we present a methodology that models four-dimensional (4D) post-failure landslide kinematics using a decaying mathematical model. This approach enables us to represent the stress relaxation for the landslide body dynamics after failure. By employing this methodology, we can overcome the weaknesses of the individual sensor in spaceborne optical and radar imaging. 3. We assessed the effectiveness of a newly designed small dihedral corner reflector for landslide monitoring. The reflector is compatible with both ascending and descending satellite orbits, while it is also suitable for applications with both high-resolution and medium-resolution satellite imagery. Furthermore, although its echoes are not as strong as those of conventional reflectors, the cost of the newly designed reflectors is reduced, with more manageable installation and maintenance. To overcome this limitation, we propose a specific selection strategy based on a probability model to identify the reflectors in satellite images

The SAR Handbook: Comprehensive Methodologies for Forest Monitoring and Biomass Estimation

This Synthetic Aperture Radar (SAR) handbook of applied methods for forest monitoring and biomass estimation has been developed by SERVIR in collaboration with SilvaCarbon to address pressing needs in the development of operational forest monitoring services. Despite the existence of SAR technology with all-weather capability for over 30 years, the applied use of this technology for operational purposes has proven difficult. This handbook seeks to provide understandable, easy-to-assimilate technical material to remote sensing specialists that may not have expertise on SAR but are interested in leveraging SAR technology in the forestry sector

Potential and limits of InSAR to characterize interseismic deformation independently of GPS data: Application to the southern San Andreas Fault system

The evaluation of long-wavelength deformation associated with interseismic strain accumulation traditionally relies on spatially sparse GPS measurements, or on high spatial-resolution InSAR velocity fields aligned to a GPS-based model. In this approach the InSAR contributes only short-wavelength deformation and the two data sets are dependent, thereby challenging the evaluation of the InSAR uncertainties and the justification of atmospheric corrections. Here we present an analysis using 7 years of Envisat InSAR data to characterize interseismic deformation along the southern San Andreas Fault (SAF) and the San Jacinto Fault (SJF) in southern California, where the SAF bifurcates onto the Mission Creek (MCF) and the Banning (BF) fault strands. We outline the processing steps for using InSAR alone to characterize both the short- and long-wavelength deformation, and evaluate the velocity field uncertainties with independent continuous GPS data. InSAR line-of-sight (LOS) and continuous GPS velocities agree within ∼1–2 mm/yr in the study area, suggesting that multiyear InSAR time series can be used to characterize interseismic deformation with a higher spatial resolution than GPS. We investigate with dislocation models the ability of this mean LOS velocity field to constrain fault slip rates and show that a single viewing geometry can help distinguish between different slip-rate scenarios on the SAF and SJF (∼35 km apart) but multiple viewing geometries are needed to differentiate slip on the MCF and BF (<12 km apart). Our results demonstrate that interseismic models of strain accumulation used for seismic hazards assessment would benefit from the consideration of InSAR mean velocity maps

The BIOMASS level 2 prototype processor : design and experimental results of above-ground biomass estimation

BIOMASS is ESA’s seventh Earth Explorer mission, scheduled for launch in 2022. The satellite will be the first P-band SAR sensor in space and will be operated in fully polarimetric interferometric and tomographic modes. The mission aim is to map forest above-ground biomass (AGB), forest height (FH) and severe forest disturbance (FD) globally with a particular focus on tropical forests. This paper presents the algorithms developed to estimate these biophysical parameters from the BIOMASS level 1 SAR measurements and their implementation in the BIOMASS level 2 prototype processor with a focus on the AGB product. The AGB product retrieval uses a physically-based inversion model, using ground-canceled level 1 data as input. The FH product retrieval applies a classical PolInSAR inversion, based on the Random Volume over Ground Model (RVOG). The FD product will provide an indication of where significant changes occurred within the forest, based on the statistical properties of SAR data. We test the AGB retrieval using modified airborne P-Band data from the AfriSAR and TropiSAR campaigns together with reference data from LiDAR-based AGB maps and plot-based ground measurements. For AGB estimation based on data from a single heading, comparison with reference data yields relative Root Mean Square Difference (RMSD) values mostly between 20% and 30%. Combining different headings in the estimation process significantly improves the AGB retrieval to slightly less than 20%. The experimental results indicate that the implemented retrieval scheme provides robust results that are within mission requirements