Long-term monitoring of geodynamic surface deformation using SAR interferometry

Thesis (Ph.D.) University of Alaska Fairbanks, 2014Synthetic Aperture Radar Interferometry (InSAR) is a powerful tool to measure surface deformation and is well suited for surveying active volcanoes using historical and existing satellites. However, the value and applicability of InSAR for geodynamic monitoring problems is limited by the influence of temporal decorrelation and electromagnetic path delay variations in the atmosphere, both of which reduce the sensitivity and accuracy of the technique. The aim of this PhD thesis research is: how to optimize the quantity and quality of deformation signals extracted from InSAR stacks that contain only a low number of images in order to facilitate volcano monitoring and the study of their geophysical signatures. In particular, the focus is on methods of mitigating atmospheric artifacts in interferograms by combining time-series InSAR techniques and external atmospheric delay maps derived by Numerical Weather Prediction (NWP) models. In the first chapter of the thesis, the potential of the NWP Weather Research & Forecasting (WRF) model for InSAR data correction has been studied extensively. Forecasted atmospheric delays derived from operational High Resolution Rapid Refresh for the Alaska region (HRRRAK) products have been compared to radiosonding measurements in the first chapter. The result suggests that the HRRR-AK operational products are a good data source for correcting atmospheric delays in spaceborne geodetic radar observations, if the geophysical signal to be observed is larger than 20 mm. In the second chapter, an advanced method for integrating NWP products into the time series InSAR workflow is developed. The efficiency of the algorithm is tested via simulated data experiments, which demonstrate the method outperforms other more conventional methods. In Chapter 3, a geophysical case study is performed by applying the developed algorithm to the active volcanoes of Unimak Island Alaska (Westdahl, Fisher and Shishaldin) for long term volcano deformation monitoring. The volcano source location at Westdahl is determined to be approx. 7 km below sea level and approx. 3.5 km north of the Westdahl peak. This study demonstrates that Fisher caldera has had continuous subsidence over more than 10 years and there is no evident deformation signal around Shishaldin peak.Chapter 1. Performance of the High Resolution Atmospheric Model HRRR-AK for Correcting Geodetic Observations from Spaceborne Radars -- Chapter 2. Robust atmospheric filtering of InSAR data based on numerical weather prediction models -- Chapter 3. Subtle motion long term monitoring of Unimak Island from 2003 to 2010 by advanced time series SAR interferometry -- Chapter 4. Conclusion and future work

Advanced satellite radar interferometry for small-scale surface deformation detection

Synthetic aperture radar interferometry (InSAR) is a technique that enables generation of Digital Elevation Models (DEMs) and detection of surface motion at the centimetre level using radar signals transmitted from a satellite or an aeroplane. Deformation observations can be performed due to the fact that surface motion, caused by natural and human activities, generates a local phase shift in the resultant interferogram. The magnitude of surface deformation can be estimated directly as a fraction of the wavelength of the transmitted signal. Moreover, differential InSAR (DInSAR) eliminates the phase signal caused by relief to yield a differential interferogram in which the signature of surface deformation can be seen. Although InSAR applications are well established, the improvement of the interferometry technique and the quality of its products is highly desirable to further enhance its capabilities. The application of InSAR encounters problems due to noise in the interferometric phase measurement, caused by a number of decorrelation factors. In addition, the interferogram contains biases owing to satellite orbit errors and atmospheric heterogeneity These factors dramatically reduce the stlectiveness of radar interferometry in many applications, and, in particular, compromise detection and analysis of small-scale spatial deformations. The research presented in this thesis aim to apply radar interferometry processing to detect small-scale surface deformations, improve the quality of the interferometry products, determine the minimum and maximum detectable deformation gradient and enhance the analysis of the interferometric phase image. The quality of DEM and displacement maps can be improved by various methods at different processing levels. One of the methods is filtering of the interferometric phase.However, while filtering reduces noise in the interferogram, it does not necessarily enhance or recover the signal. Furthermore, the impact of the filter can significantly change the structure of the interferogram. A new adaptive radar interferogram filter has been developed and is presented herein. The filter is based on a modification to the Goldstein radar interferogram filter making the filter parameter dependent on coherence so that incoherent areas are filtered more than coherent areas. This modification minimises the loss of signal while still reducing the level of noise. A methodology leading to the creation of a functional model for determining minimum and maximum detectable deformation gradient, in terms of the coherence value, has been developed. The sets of representative deformation models have been simulated and the associated phase from these models has been introduced to real SAR data acquired by ERS-1/2 satellites. A number of cases of surface motion with varying magnitudes and spatial extent have been simulated. In each case, the resultant surface deformation has been compared with the 'true' surface deformation as defined by the deformation model. Based on those observations, the functional model has been developed. Finally, the extended analysis of the interferometric phase image using a wavelet approach is presented. The ability of a continuous wavelet transform to reveal the content of the wrapped phase interferogram, such as (i) discontinuities, (ii) extent of the deformation signal, and (iii) the magnitude of the deformation signal is examined. The results presented represent a preliminary study revealing the wavelet method as a promising technique for interferometric phase image analysis

Geodetic monitoring of complex shaped infrastructures using Ground-Based InSAR

In the context of climate change, alternatives to fossil energies need to be used as much as possible to produce electricity. Hydroelectric power generation through the utilisation of dams stands out as an exemplar of highly effective methodologies in this endeavour. Various monitoring sensors can be installed with different characteristics w.r.t. spatial resolution, temporal resolution and accuracy to assess their safe usage. Among the array of techniques available, it is noteworthy that ground-based synthetic aperture radar (GB-SAR) has not yet been widely adopted for this purpose. Despite its remarkable equilibrium between the aforementioned attributes, its sensitivity to atmospheric disruptions, specific acquisition geometry, and the requisite for phase unwrapping collectively contribute to constraining its usage. Several processing strategies are developed in this thesis to capitalise on all the opportunities of GB-SAR systems, such as continuous, flexible and autonomous observation combined with high resolutions and accuracy. The first challenge that needs to be solved is to accurately localise and estimate the azimuth of the GB-SAR to improve the geocoding of the image in the subsequent step. A ray tracing algorithm and tomographic techniques are used to recover these external parameters of the sensors. The introduction of corner reflectors for validation purposes confirms a significant error reduction. However, for the subsequent geocoding, challenges persist in scenarios involving vertical structures due to foreshortening and layover, which notably compromise the geocoding quality of the observed points. These issues arise when multiple points at varying elevations are encapsulated within a singular resolution cell, posing difficulties in pinpointing the precise location of the scattering point responsible for signal return. To surmount these hurdles, a Bayesian approach grounded in intensity models is formulated, offering a tool to enhance the accuracy of the geocoding process. The validation is assessed on a dam in the black forest in Germany, characterised by a very specific structure. The second part of this thesis is focused on the feasibility of using GB-SAR systems for long-term geodetic monitoring of large structures. A first assessment is made by testing large temporal baselines between acquisitions for epoch-wise monitoring. Due to large displacements, the phase unwrapping can not recover all the information. An improvement is made by adapting the geometry of the signal processing with the principal component analysis. The main case study consists of several campaigns from different stations at Enguri Dam in Georgia. The consistency of the estimated displacement map is assessed by comparing it to a numerical model calibrated on the plumblines data. It exhibits a strong agreement between the two results and comforts the usage of GB-SAR for epoch-wise monitoring, as it can measure several thousand points on the dam. It also exhibits the possibility of detecting local anomalies in the numerical model. Finally, the instrument has been installed for continuous monitoring for over two years at Enguri Dam. An adequate flowchart is developed to eliminate the drift happening with classical interferometric algorithms to achieve the accuracy required for geodetic monitoring. The analysis of the obtained time series confirms a very plausible result with classical parametric models of dam deformations. Moreover, the results of this processing strategy are also confronted with the numerical model and demonstrate a high consistency. The final comforting result is the comparison of the GB-SAR time series with the output from four GNSS stations installed on the dam crest. The developed algorithms and methods increase the capabilities of the GB-SAR for dam monitoring in different configurations. It can be a valuable and precious supplement to other classical sensors for long-term geodetic observation purposes as well as short-term monitoring in cases of particular dam operations

Phi-Net: Deep Residual Learning for InSAR Parameters Estimation

Nowadays, deep learning (DL) finds application in a large number of scientific fields, among which the estimation and the enhancement of signals disrupted by the noise of different natures. In this article, we address the problem of the estimation of the interferometric parameters from synthetic aperture radar (SAR) data. In particular, we combine convolutional neural networks together with the concept of residual learning to define a novel architecture, named Phi-Net, for the joint estimation of the interferometric phase and coherence. Phi-Net is trained using synthetic data obtained by an innovative strategy based on the theoretical modeling of the physics behind the SAR acquisition principle. This strategy allows the network to generalize the estimation problem with respect to: 1) different noise levels; 2) the nature of the imaged target on the ground; and 3) the acquisition geometry. We then analyze the Phi-Net performance on an independent data set of synthesized interferometric data, as well as on real InSAR data from the TanDEM-X and Sentinel-1 missions. The proposed architecture provides better results with respect to state-of-the-art InSAR algorithms on both synthetic and real test data. Finally, we perform an application-oriented study on the retrieval of the topographic information, which shows that Phi-Net is a strong candidate for the generation of high-quality digital elevation models at a resolution close to the one of the original single-look complex data

Displacement Estimation by Maximum Likelihood Texture Tracking

International audienceThis paper presents a novel method to estimate displacement by maximum-likelihood (ML) texture tracking. The observed polarimetric synthetic aperture radar (PolSAR) data-set is composed by two terms: the scalar texture parameter and the speckle component. Based on the Spherically Invariant Random Vectors (SIRV) theory, the ML estimator of the texture is computed. A generalization of the ML texture tracking based on the Fisher probability density function (pdf) modeling is introduced. For random variables with Fisher distributions, the ratio distribution is established. The proposed method is tested with both simulated PolSAR data and spaceborne PolSAR images provided by the TerraSAR-X (TSX) and the RADARSAT-2 (RS-2) sensors

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

Non-local methods for InSAR parameters estimation

In the thesis work the nonlocal paradigm has been investigated in the framework of Multitemporal SAR Interferometry, e.g. Differential Interferometry, Tomography, etc., and single InSAR pair, e.g. DEM generation. In the former, Adaptive Multi-Looking methods have been developed for the generation of interferometric data-stacks. Following the nonlocal approach, the proposed methods rely only on similar pixels according to a suitable similarity measure that exploits the stack's temporal information. An hybrid approach that jointly uses the nonlocal paradigm and transform domain filtering has been investigated for InSAR pair phase estimation. On the track of the BM3D and SARBM3D algorithms, different approaches to the filtering in the transform domain are investigated. Furthermore, a novel approach to the similarity computation and filtering, based on a relative-topography content of the interferometric phase rather than its absolute value, is proposed

Tracking hidden crisis in India's capital from space: implications of unsustainable groundwater use.

Funder: Helmholtz-Zentrum Potsdam Deutsches GeoForschungsZentrum - GFZNational Capital Region (NCR, Delhi) in India is one of the fastest-growing metropolitan cities which is facing a severe water crisis due to increasing water demand. The over-extraction of groundwater, particularly from its unconsolidated alluvial deposits makes the region prone to subsidence. In this study, we investigated the effects of plummeting groundwater levels on land surface elevations in Delhi NCR using Sentinel-1 datasets acquired during the years 2014-2020. Our analysis reveals two distinct subsidence features in the study area with rates exceeding 11 cm/year in Kapashera-an urban village near IGI airport Delhi, and 3 cm/year in Faridabad throughout the study period. The subsidence in these two areas are accelerating and follows the depleting groundwater trend. The third region, Dwarka shows a shift from subsidence to uplift during the years which can be attributed to the strict government policies to regulate groundwater use and incentivizing rainwater harvesting. Further analysis using a classified risk map based on hazard risk and vulnerability approach highlights an approximate area of 100 square kilometers to be subjected to the highest risk level of ground movement, demanding urgent attention. The findings of this study are highly relevant for government agencies to formulate new policies against the over-exploitation of groundwater and to facilitate a sustainable and resilient groundwater management system in Delhi NCR