DInSAR deformation time series for monitoring urban areas: The impact of the second generation SAR systems

We investigate the capability improvement of the DInSAR techniques to map deformation phenomena affecting urban areas, by performing a comparative analysis of the deformation time series retrieved by applying the full resolution Small BAseline Subset (SBAS) DInSAR technique to selected sequences of SAR data acquired by the ENVISAT, RADARSAT-1 and COSMO-SkyMed (CSK) SAR data. The presented study, focused on the city of Napoli (Italy), allows us to quantify the dramatic increase of the DInSAR coherent pixel density achieved by exploiting the high resolution X-Band CSK SAR images with respect to the RADARSAT-1 and ENVISAT products, respectively; this permits us to analyze nearly all the structures located within the investigated urbanized area and, in many cases, also portions of a same building. © 2012 IEEE

Post-failure evolution analysis of a rainfall-triggered landslide by multi-temporal interferometry SAR approaches integrated with geotechnical analysis

Persistent Scatterers Interferometry (PSI) represents one of the most powerful techniques for Earth's surface deformation processes' monitoring, especially for long-term evolution phenomena. In this work, a dataset of 34 TerraSAR-X StripMap images (October 2013–October 2014) has been processed by two PSI techniques - Coherent Pixel Technique-Temporal Sublook Coherence (CPT-TSC) and Small Baseline Subset (SBAS) - in order to study the evolution of a slow-moving landslide which occurred on February 23, 2012 in the Papanice hamlet (Crotone municipality, southern Italy) and induced by a significant rainfall event (185 mm in three days). The mass movement caused structural damage (buildings' collapse), and destruction of utility lines (gas, water and electricity) and roads. The results showed analogous displacement rates (30–40 mm/yr along the Line of Sight – LOS-of the satellite) with respect to the pre-failure phase (2008–2010) analyzed in previous works. Both approaches allowed detect the landslide-affected area, however the higher density of targets identified by means of CPT-TSC enabled to analyze in detail the slope behavior in order to design possible mitigation interventions. For this aim, a slope stability analysis has been carried out, considering the comparison between groundwater oscillations and time-series of displacement. Hence, the crucial role of the interaction between rainfall and groundwater level has been inferred for the landslide triggering. In conclusion, we showed that the integration of geotechnical and remote sensing approaches can be seen as the best practice to support stakeholders to design remedial works.Peer ReviewedPostprint (author's final draft

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

First insights on the potential of Sentinel-1 for landslides detection

This paper illustrates the potential of Sentinel-1 for landslide detection, Accepted 23 March 2016 mapping and characterization with the aim of updating inventory maps and monitoring landslide activity. The study area is located in Molise, one of the smallest regions of Italy, where landslide processes are frequent. The results achieved by integrating Differential Synthetic Aperture Radar Interferometry (DInSAR) deformation maps and time series, and Geographical Information System (GIS) multilayer analysis (optical, geological, geomorphological, etc.) are shown. The adopted methodology is described followed by an analysis of future perspectives. Sixty-two landslides have been detected, thus allowing the updating of pre-existing landslide inventory maps. The results of our ongoing research show that Sentinel-1 might represent a significant improvement in terms of exploitation of SAR data for landslide mapping and monitoring due to both the shorter revisit time (up to 6 days in the close future) and the wavelength used, which determine an higher coherence compared to other SAR sensors

Basin scale assessment of landslides geomorphological setting by advanced InSAR analysis

An extensive investigation of more than 90 landslides affecting a small river basin in Central Italy was performed by combining field surveys and remote sensing techniques. We thus defined the geomorphological setting of slope instability processes. Basic information, such as landslides mapping and landslides type definition, have been acquired thanks to geomorphological field investigations and multi-temporal aerial photos interpretation, while satellite SAR archive data (acquired by ERS and Envisat from 1992 to 2010) have been analyzed by means of A-DInSAR (Advanced Differential Interferometric Synthetic Aperture Radar) techniques to evaluate landslides past displacements patterns. Multi-temporal assessment of landslides state of activity has been performed basing on geomorphological evidence criteria and past ground displacement measurements obtained by A-DInSAR. This step has been performed by means of an activity matrix derived from information achieved thanks to double orbital geometry. Thanks to this approach we also achieved more detailed knowledge about the landslides kinematics in time and space

Spatio-temporal quality indicators for differential interferometric Synthetic Aperture Radar data

Satellite-based interferometric synthetic aperture radar (InSAR) is an invaluable technique in the detection and monitoring of changes on the surface of the earth. Its high spatial coverage, weather friendly and remote nature are among the advantages of the tool. The multi-temporal differential InSAR (DInSAR) methods in particular estimate the spatio-temporal evolution of deformation by incorporating information from multiple SAR images. Moreover, opportunities from the DInSAR techniques are accompanied by challenges that affect the final outputs. Resolving the inherent ambiguities of interferometric phases, especially in areas with a high spatio-temporal deformation gradient, represents the main challenge. This brings the necessity of quality indices as important DInSAR data processing tools in achieving ultimate processing outcomes. Often such indices are not provided with the deformation products. In this work, we propose four scores associated with (i) measurement points, (ii) dates of time series, (iii) interferograms and (iv) images involved in the processing. These scores are derived from a redundant set of interferograms and are calculated based on the consistency of the unwrapped interferometric phases in the frame of a least-squares adjustment. The scores reflect the occurrence of phase unwrapping errors and represent valuable input for the analysis and exploitation of the DInSAR results. The proposed tools were tested on 432,311 points, 1795 interferograms and 263 Sentinel-1 single look complex images by employing the small baseline technique in the PSI processing chain, PSIG of the geomatics division of the Centre Tecnològic de Telecomunicacions de Catalunya (CTTC). The results illustrate the importance of the scores—mainly in the interpretation of the DInSAR outputs.This research was partially funded by the Agència de Gestió d’Ajuts Universitaris i de Recerca (AGAUR), Generalitat de Catalunya—through a grant for the recruitment of early-stage research staff (Ref: 2021FI_B2 00186).Peer ReviewedPostprint (published version

Deformation due to the 17 August 1999 Izmit, Turkey, earthquake measured from SPOT images

The geometry of the ruptured areas and the coseismic slip distribution data are key to highlighting the behavior of seismic faults. This information is generally retrieved from field investigations and geodetic measurements or synthetic aperture radar (SAR) interferometry. Here we show that SPOT images can also be used to accurately map the fault zone and to determine the slip distribution by subpixel correlation of images acquired before and after an earthquake. The measured slip includes the contribution of possible distributed shear that might not be clearly expressed in surface ruptures and smoothes out possible along-strike variability due to near-surface fault complexities. We apply the technique to the M_s = 7.4, 1999, Izmit earthquake. Our results reveal a <100-m-wide and very linear fault zone that can be traced for 70 km from Gölcük to Akyazi, along which supershear rupture has been inferred. The obtained slip distribution compares well with the field measurements and is consistent with ground deformation measured at some distance from the fault zone using SAR images. Very little deformation was accommodated off the main fault plane. Maximum slip is observed near Sapanca lake at a small fault jog that has probably influenced rupture propagation

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