Advanced SBAS-DInSAR technique for controlling large civil infrastructures: an application to the Genzano di Lucania dam

Monitoring surface deformation on dams is commonly carried out by in situ geodetic surveying, which is time consuming and characterized by some limitations in space coverage and frequency. More recently microwave satellite-based technologies, such as advanced-DInSAR (Differential Synthetic Aperture Radar Interferometry), have allowed the integration and improvement of the observation capabilities of ground-based methods thanks to their effectiveness in collecting displacement measurements on many non-destructive control points, corresponding to radar reflecting targets. The availability of such a large number of points of measurement, which are distributed along the whole structure and are characterized by millimetric accuracy on displacement rates, can be profitably adopted for the calibration of numerical models. These models are implemented to simulate the structural behaviour of a dam under conditions of stress thus improving the ability to maintain safety standards. In this work, after having analysed how advanced DInSAR can effectively enhance the results from traditional monitoring systems that provide comparable accuracy measurements on a limited number of points, an FEM model of the Genzano di Lucania earth dam is developed and calibrated. This work is concentrated on the advanced DInSAR technique referred to as Small BAseline Subset (SBAS) approach, benefiting from its capability to generate deformation time series at full spatial resolution and from multi-sensor SAR data, to measure the vertical consolidation displacement of the Genzano di Lucania earth dam

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

Magma and fluid migration at Yellowstone Caldera in the last three decades inferred from InSAR, leveling and gravity measurements

We studied the Yellowstone caldera geological unrest between 1977 and 2010 by investigating temporal changes in differential Interferometric Synthetic Aperture Radar (InSAR), precise spirit leveling and gravity measurements. The analysis of the 1992–2010 displacement time series, retrieved by applying the SBAS InSAR technique, allowed the identification of three areas of deformation: (i) the Mallard Lake (ML) and Sour Creek (SC) resurgent domes, (ii) a region close to the Northern Caldera Rim (NCR), and (iii) the eastern Snake River Plain (SRP). While the eastern SRP shows a signal related to tectonic deformation, the other two regions are influenced by the caldera unrest. We removed the tectonic signal from the InSAR displacements, and we modeled the InSAR, leveling, and gravity measurements to retrieve the best fitting source parameters. Our findings confirmed the existence of different distinct sources, beneath the brittle-ductile transition zone, which have been intermittently active during the last three decades. Moreover, we interpreted our results in the light of existing seismic tomography studies. Concerning the SC dome, we highlighted the role of hydrothermal fluids as the driving force behind the 1977–1983 uplift; since 1983–1993 the deformation source transformed into a deeper one with a higher magmatic component. Furthermore, our results support the magmatic nature of the deformation source beneath ML dome for the overall investigated period. Finally, the uplift at NCR is interpreted as magma accumulation, while its subsidence could either be the result of fluids migration outside the caldera or the gravitational adjustment of the source from a spherical to a sill-like geometr

Comparison of Small Baseline Interferometric SAR Processors for Estimating Ground Deformation

The small Baseline Synthetic Aperture Radar (SAR) Interferometry (SBI) technique has been widely and successfully applied in various ground deformation monitoring applications. Over the last decade, a variety of SBI algorithms have been developed based on the same fundamental concepts. Recently developed SBI toolboxes provide an open environment for researchers to apply different SBI methods for various purposes. However, there has been no thorough discussion that compares the particular characteristics of different SBI methods and their corresponding performance in ground deformation reconstruction. Thus, two SBI toolboxes that implement a total of four SBI algorithms were selected for comparison. This study discusses and summarizes the main differences, pros and cons of these four SBI implementations, which could help users to choose a suitable SBI method for their specific application. The study focuses on exploring the suitability of each SBI module under various data set conditions, including small/large number of interferograms, the presence or absence of larger time gaps, urban/vegetation ground coverage, and temporally regular/irregular ground displacement with multiple spatial scales. Within this paper we discuss the corresponding theoretical background of each SBI method. We present a performance analysis of these SBI modules based on two real data sets characterized by different environmental and surface deformation conditions. The study shows that all four SBI processors are capable of generating similar ground deformation results when the data set has sufficient temporal sampling and a stable ground backscatter mechanism like urban area. Strengths and limitations of different SBI processors were analyzed based on data set configuration and environmental conditions and are summarized in this paper to guide future users of SBI techniques

Land subsidence over oilfields in the Yellow River Delta

Subsidence in river deltas is a complex process that has both natural and human causes. Increasing human activities like aquaculture and petroleum extraction are affecting the Yellow River delta, and one consequence is subsidence. The purpose of this study is to measure the surface displacements in the Yellow River delta region and to investigate the corresponding subsidence source. In this paper, the Stanford Method for Persistent Scatterers (StaMPS) package was employed to process Envisat ASAR images collected between 2007 and 2010. Consistent results between two descending tracks show subsidence with a mean rate up to 30 mm/yr in the radar line of sight direction in Gudao Town (oilfield), Gudong oilfield and Xianhe Town of the delta, each of which is within the delta, and also show that subsidence is not uniform across the delta. Field investigation shows a connection between areas of non-uniform subsidence and of petroleum extraction. In a 9 km2 area of the Gudao Oilfield, a poroelastic disk reservoir model is used to model the InSAR derived displacements. In general, good fits between InSAR observations and modeled displacements are seen. The subsidence observed in the vicinity of the oilfield is thus suggested to be caused by fluid extraction

Coherency Matrix Decomposition-Based Polarimetric Persistent Scatterer Interferometry

© 2019 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes,creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.The rationale of polarimetric optimization techniques is to enhance the phase quality of the interferograms by combining adequately the different polarization channels available to produce an improved one. Different approaches have been proposed for polarimetric persistent scatterer interferometry (PolPSI). They range from the simple and computationally efficient BEST, where, for each pixel, the polarimetric channel with the best response in terms of phase quality is selected, to those with high-computational burden like the equal scattering mechanism (ESM) and the suboptimum scattering mechanism (SOM). BEST is fast and simple, but it does not fully exploit the potentials of polarimetry. On the other side, ESM explores all the space of solutions and finds the optimal one but with a very high-computational burden. A new PolPSI algorithm, named coherency matrix decomposition-based PolPSI (CMD-PolPSI), is proposed to achieve a compromise between phase optimization and computational cost. Its core idea is utilizing the polarimetric synthetic aperture radar (PolSAR) coherency matrix decomposition to determine the optimal polarization channel for each pixel. Three different PolSAR image sets of both full- (Barcelona) and dual-polarization (Murcia and Mexico City) are used to evaluate the performance of CMD-PolPSI. The results show that CMD-PolPSI presents better optimization results than the BEST method by using either $D_{\mathrm{ A}}$ or temporal mean coherence as phase quality metrics. Compared with the ESM algorithm, CMD-PolPSI is 255 times faster but its performance is not optimal. The influence of the number of available polarization channels and pixel's resolutions on the CMD-PolPSI performance is also discussed.Peer ReviewedPostprint (author's final draft

Performance Analysis of GPS Aided Geo Augmented Navigation (GAGAN) Over Sri Lanka

Satellite-Based Augmentation Systems (SBAS) are being developed worldwide due to their unique advantage of wide area coverage. GPS Aided Geo Augmented Navigation (GAGAN) is an Indian implementation of SBAS, with three (03) geo stationary satellites in space covering a huge area even beyond Indian Territory. This study focused on analyzing the improvement in position solution with GAGAN corrections over Sri Lanka. In order to test its performances, several dual and single frequency GNSS receivers have been used for this experiment and one receiver is configured as SBAS receiver and other two were kept as GPS stand-alone receivers. For investigate its coverage over Sri Lanka; observations were carried out over seven (07) known control stations of 6 different districts. At each of the tested stations the GAGAN active L1 receiver has always shown a significant accuracy improvement over L1 uncorrected observations. Further, five out of the 7 observation locations the calculated average 3D positional errors are lower than 1m. Almost 79% of observations (out of 24 hours of observations) have shown acceptable 3D positional accuracy, of less than 1m, for many spatial data collection applications. However, the local DGPS correction has shown higher reliability than GAGAN corrections with almost 85% of observations with less tan 1m, 3D positional error.Satellite-Based Augmentation Systems (SBAS) are being developed worldwide due to their unique advantage of wide area coverage. GPS Aided Geo Augmented Navigation (GAGAN) is an Indian implementation of SBAS, with three (03) geo stationary satellites in space covering a huge area even beyond Indian Territory. This study focused on analyzing the improvement in position solution with GAGAN corrections over Sri Lanka. In order to test its performances, several dual and single frequency GNSS receivers were used in this experiment, one receiver was configured as SBAS receiver and other two were kept as GPS stand-alone receivers. Observations were carried out over seven (07) known control stations of six (6) different districts to investigate its coverage over Sri Lanka. At each of the tested stations the GAGAN active L1 receiver has always shown a significant accuracy improvement over L1 uncorrected observations. Further, five out of the seven (7) observation locations the calculated average 3D positional errors were lower than 1 m. Almost 79% of observations (out of 24 hours of observations) have shown acceptable 3D positional accuracy, of less than 1m, for many spatial data collection applications. However, the local DGPS correction has shown higher reliability than GAGAN corrections with almost 85% of observations with less than 1 m, 3D positional error

A user-oriented methodology for DInSAR time series analysis and interpretation: landslides and subsidence case studies

Recent advances in multi-temporal Differential Synthetic Aperture Radar (SAR) Interferometry (DInSAR) have greatly improved our capability to monitor geological processes. Ground motion studies using DInSAR require both the availability of good quality input data and rigorous approaches to exploit the retrieved Time Series (TS) at their full potential. In this work we present a methodology for DInSAR TS analysis, with particular focus on landslides and subsidence phenomena. The proposed methodology consists of three main steps: (1) pre-processing, i.e., assessment of a SAR Dataset Quality Index (SDQI) (2) post-processing, i.e., application of empirical/stochastic methods to improve the TS quality, and (3) trend analysis, i.e., comparative implementation of methodologies for automatic TS analysis. Tests were carried out on TS datasets retrieved from processing of SAR imagery acquired by different radar sensors (i.e., ERS-1/2 SAR, RADARSAT-1, ENVISAT ASAR, ALOS PALSAR, TerraSAR-X, COSMO-SkyMed) using advanced DInSAR techniques (i.e., SqueeSAR™, PSInSAR™, SPN and SBAS). The obtained values of SDQI are discussed against the technical parameters of each data stack (e.g., radar band, number of SAR scenes, temporal coverage, revisiting time), the retrieved coverage of the DInSAR results, and the constraints related to the characterization of the investigated geological processes. Empirical and stochastic approaches were used to demonstrate how the quality of the TS can be improved after the SAR processing, and examples are discussed to mitigate phase unwrapping errors, and remove regional trends, noise and anomalies. Performance assessment of recently developed methods of trend analysis (i.e., PS-Time, Deviation Index and velocity TS) was conducted on two selected study areas in Northern Italy affected by land subsidence and landslides. Results show that the automatic detection of motion trends enhances the interpretation of DInSAR data, since it provides an objective picture of the deformation behaviour recorded through TS and therefore contributes to the understanding of the on-going geological processes

An integrated approach for monitoring soil settlements at the VIrgo site

The Virgo detector, currently in its 2nd generation configuration Advanced Virgo (AdV), is a Michelson interferometer aimed at the gravitational waves research and at opening a new window on the study of the Universe. It is made of two orthogonal arms being each 3 kilometers long and is located at the site of the European Gravitational Observatory (EGO), in the countryside near Pisa, Italy. After the construction of the Virgo facilities completed in 2002, over the years a steady subsidence process has been observed as a consequence of the building and embankment overloads. In consideration of the subsoil characteristics, whose surface portion is mainly formed by a 25÷60 m thickness layer of clay with limited thin layers of sands, the evolution of settlements was expected and properly considered for the design of the civil engineering infrastructures, so that the vacuum tubes can be readjusted to keep the original alignment. However, along 15 years of time life, the initial estimates of the expected displacements were continuously compared with the observed effects. The measured settlements have been regularly monitored and adopted for implementing the necessary realignment activities. This paper reports the monitoring activities conducted over the years, mainly consisting of regular high accuracy levelling surveys, periodically integrated by GPS and classical theodolite measurements. These sets of measurement were adopted to perform the Virgo realignment procedure needed to keep the interferometer rigidly tied in a 3x3km plane. In order to improve the knowledge on the trend of the settlements affecting the Virgo infrastructures, an analysis based on differential interferometry using satellite Synthetic Aperture Radar (SAR) data has been performed and compared with the outcome from in-situ data

Surface deformation of active volcanic areas retrieved with the SBAS-DInSAR technique: an overview

This paper presents a comprehensive overview of the surface deformation retrieval capability of the Differential Synthetic Aperture Radar Interferometry (DInSAR) algorithm, referred to as Small BAseline Subset (SBAS) technique, in the context of active volcanic areas. In particular, after a brief description of the algorithm some experiments relevant to three selected case-study areas are presented. First, we concentrate on the application of the SBAS algorithm to a single-orbit scenario, thus considering a set of SAR data composed by images acquired on descending orbits by the European Remote Sensing (ERS) radar sensors and relevant to the Long Valley caldera (eastern California) area. Subsequently, we address the capability of the SBAS technique in a multipleorbit context by referring to Mt. Etna volcano (southern Italy) test site, with respect to which two different ERS data set, composed by images acquired both on ascending and descending orbits, are available. Finally, we take advantage of the capability of the algorithm to work in a multi-platform scenario by jointly exploiting two different sets of SAR images collected by the ERS and the Environment Satellite (ENVISAT) radar sensors in the Campi Flegrei caldera (southern Italy) area. The presented results demonstrate the effectiveness of the algorithm to investigate the deformation field in active volcanic areas and the potential of the DInSAR methodologies within routine surveillance scenario