Generation Of Deformation Time-series Research Articles

Spatiotemporal Evolution of Ground Subsidence and Extensional Basin Bedrock Organization: An Application of Multitemporal Multi-Satellite SAR Interferometry

Since the early 1990s, the European (ESA) and Italian (ASI) space agencies have managed and distributed a huge amount of satellite-recorded SAR data to the research community and private industries. Moreover, the availability of advanced cloud computing services implementing different multi-temporal SAR interferometry techniques allows the generation of deformation time series from massive SAR images. We exploit the information provided by a large PS dataset to determine the temporal trend of ground deformation and the relative deformation rate with millimetric accuracy to analyze the spatial and temporal distribution of land subsidence induced by water pumping from a deep confined aquifer in the Northern Valle Umbra Basin (Central Italy), exploiting 24 years of Permanent Scatterers—interferometric SAR data archives. The SAR images were acquired between 1992 and 2016 by satellites ERS1/2 and ENVISAT, the Sentinel 1 ESA missions and the COSMO-SkyMed ASI mission. We observed ground velocities and deformation geometries between 1992 and 2016, with displacements of more than 70 cm and velocities of up to 55 mm/yr. The results suggest that the shape and position of the surface ground displacement are controlled by the fault activity hidden under the valley deposits.

Automatic Generation of Sentinel-1 Continental Scale DInSAR Deformation Time Series through an Extended P-SBAS Processing Pipeline in a Cloud Computing Environment

We present in this work an advanced processing pipeline for continental scale differential synthetic aperture radar (DInSAR) deformation time series generation, which is based on the parallel small baseline subset (P-SBAS) approach and on the joint exploitation of Sentinel-1 (S-1) interferometric wide swath (IWS) SAR data, continuous global navigation satellite system (GNSS) position time-series, and cloud computing (CC) resources. We first briefly describe the basic rationale of the adopted P-SBAS processing approach, tailored to deal with S-1 IWS SAR data and to be implemented in a CC environment, highlighting the innovative solutions that have been introduced in the processing chain we present. They mainly consist in a series of procedures that properly exploit the available GNSS time series with the aim of identifying and filtering out possible residual atmospheric artifacts that may affect the DInSAR measurements. Moreover, significant efforts have been carried out to improve the P-SBAS processing pipeline automation and robustness, which represent crucial issues for interferometric continental scale analysis. Then, a massive experimental analysis is presented. In this case, we exploit: (i) the whole archive of S-1 IWS SAR images acquired over a large portion of Europe, from descending orbits, (ii) the continuous GNSS position time series provided by the Nevada Geodetic Laboratory at the University of Nevada, Reno, USA (UNR-NGL) available for the investigated area, and (iii) the ONDA platform, one of the Copernicus Data and Information Access Services (DIAS). The achieved results demonstrate the capability of the proposed solution to successfully retrieve the DInSAR time series relevant to such a huge area, opening new scenarios for the analysis and interpretation of these ground deformation measurements.

National Scale Surface Deformation Time Series Generation through Advanced DInSAR Processing of Sentinel-1 Data within a Cloud Computing Environment

We present an automatic pipeline implemented within the Amazon Web Services (AWS) Cloud Computing platform for the interferometric processing of large Sentinel-1 (S1) multi-temporal SAR datasets, aimed at analyzing Earth surface deformation phenomena at wide spatial scale. The developed processing chain is based on the advanced DInSAR approach referred to as Small BAseline Subset (SBAS) technique, which allows producing, with centimeter to millimeter accuracy, surface deformation time series and the corresponding mean velocity maps from a temporal sequence of SAR images. The implemented solution addresses the aspects relevant to i) S1 input data archiving; ii) interferometric processing of S1 data sequences, performed in parallel on the AWS computing nodes through both multi-node and multi-core programming techniques; iii) storage of the generated interferometric products. The experimental results are focused on a national scale DInSAR analysis performed over the whole Italian territory by processing 18 S1 slices acquired from descending orbits between March 2015 and April 2017, corresponding to 2612 S1 acquisitions. Our analysis clearly shows that an effective integration of advanced remote sensing methods and new ICT technologies can successfully contribute to deeply investigate the Earth System processes and to address new challenges within the Big Data EO scenario.

The Parallel SBAS Approach for Sentinel-1 Interferometric Wide Swath Deformation Time-Series Generation: Algorithm Description and Products Quality Assessment

We present an advanced differential synthetic aperture radar (SAR) interferometry (DInSAR) processing chain, based on the Parallel Small BAseline Subset (P-SBAS) technique, for the efficient generation of deformation time series from Sentinel-1 (S-1) interferometric wide (IW) swath SAR data sets. We first discuss an effective solution for the generation of high-quality interferograms, which properly accounts for the peculiarities of the terrain observation with progressive scans (TOPS) acquisition mode used to collect S-1 IW SAR data. These data characteristics are also properly accounted within the developed processing chain, taking full advantage from the burst partitioning. Indeed, such data structure represents a key element in the proposed P-SBAS implementation of the S-1 IW processing chain, whose migration into a cloud computing (CC) environment is also envisaged. An extensive experimental analysis, which allows us to assess the quality of the obtained interferometric products, is presented. To do this, we apply the developed S-1 IW P-SBAS processing chain to the overall archive acquired from descending orbits during the March 2015-April 2017 time span over the whole Italian territory, consisting in 2740 S-1 slices. In particular, the quality of the final results is assessed through a large-scale comparison with the GPS measurements relevant to nearly 500 stations. The mean standard deviation value of the differences between the DInSAR and the GPS time series (projected in the radar line of sight) is less than 0.5 cm, thus confirming the effectiveness of the implemented solution. Finally, a discussion about the performance achieved by migrating the developed processing chain within the Amazon Web Services CC environment is addressed, highlighting that a two-year data set relevant to a standard S-1 IW slice can be reliably processed in about 30 h.The presented results demonstrate the capability of the implemented P-SBAS approach to efficiently and effectively process large S-1 IW data sets relevant to extended portions of the earth surface, paving the way to the systematic generation of advanced DInSAR products to monitor ground displacements at a very wide spatial scale.

The Parallel SBAS Approach for Sentinel-1 Interferometric Wide Swath Deformation Time-Series Generation: Algorithm Description and Products Quality Assessment

The Parallel SBAS Approach for Sentinel-1 Interferometric Wide Swath Deformation Time-Series Generation: Algorithm Description and Products Quality Assessment

Long-term ERS/ENVISAT deformation time-series generation at full spatial resolution via the extended SBAS technique

We extend the small baseline subset (SBAS) differential synthetic aperture radar (SAR) interferometry (DInSAR) approach to allow the generation of deformation time-series by processing, at the full spatial resolution scale, long sequences of European Remote Sensing (ERS-1 and ERS-2) and Environmental Satellite (ENVISAT) SAR data acquired with the same illumination geometry. In particular, we avoid the generation of ERS/ENVISAT cross-interferograms, which are severely affected by noise phenomena due to the carrier frequency separation of the two SAR systems, and we focus on single-platform interferograms only (i.e. ERS/ERS and ENVISAT/ENVISAT interferograms) that are properly combined by applying the singular value decomposition (SVD)-based SBAS approach. Moreover, we exploit the Doppler centroid variations of the post-2000 acquisitions of the ERS-2 sensor and the carrier frequency difference between the ERS-1/2 and the ENVISAT systems, in order to maximize the number of investigated SAR pixels and to improve their geocoding. The presented results, achieved on two data sets relevant to the Napoli Bay area and to the Murge region, both located in southern Italy, confirm the effectiveness of the extended SBAS technique and demonstrate the relevance of deformation analysis carried out at the scale of single buildings or human-made structures with more than 15 years of ERS and ENVISAT acquisitions.

SBAS-Based Satellite Orbit Correction for the Generation of DInSAR Time-Series: Application to RADARSAT-1 Data

We present an algorithm aimed at correcting satellite orbit information for the generation of differential SAR interferometry (DInSAR) deformation time-series. Our approach exploits small baseline differential interferograms, to preserve their spatial coherence, and is directly compatible with the Small BAseline Subset (SBAS) DInSAR technique. In particular, the algorithm investigates the differential phase gradient directly computed from the wrapped interferograms, and is focused on the estimation of the perpendicular baseline and of the parallel baseline azimuth rate components, separately performed along the range and azimuth directions, respectively. Starting from the estimations carried out on the interferograms, we then retrieve the orbit correction associated with each SAR acquisition of our time-series by solving a system of linear equations via the SVD method, extending the SBAS inversion concept also to the orbit estimation problem. Key application of this technique is the generation of deformation time-series from interferometric sequences of RADARSAT-1 SAR acquisitions which are available for several areas in the world, but are characterized by significantly low accuracy of the or bit information. The presented results, obtained by processing a data set consisting of 33 RADARSAT-1 images of Big Island at Hawaii, show that we may retrieve DInSAR time-series with sub centimeter accuracy, thus confirming the effectiveness of the pro posed technique.

New Advances of the Extended Minimum Cost Flow Phase Unwrapping Algorithm for SBAS-DInSAR Analysis at Full Spatial Resolution

We present an efficient space-time phase unwrapping (PhU) algorithm that allows us to process sequences of multitemporal full resolution differential synthetic aperture radar (SAR) interferograms for the generation of deformation time-series. The core of the proposed technique, dealing with sparse data grids, is represented by the extended minimum cost flow (MCF) (EMCF) PhU algorithm that was originally developed for the analysis of sequences of multilook interferograms. In particular, our method relies on the joint analysis of the spatial and temporal relationships among a set of properly selected multitemporal differential interferograms, which are compatible with the Small BAseline subset (SBAS) deformation time-series technique. The key point of the approach is the idea to split the complex MCF network problem, representing the overall PhU operation, into that of simpler subnetworks. More precisely, we start by identifying and solving a primary network that involves a proper selection of coherent pixels of the computed interferograms, representing the backbone structure of the overall network. Subsequently, this result is applied for constraining the solution of the subnetworks connected to the primary one, involving the entire set of analyzed pixels. To achieve this task, we solve a constrained optimization problem based on the computation of a constrained Delaunay triangulation in the azimuth/range domain. The overall procedure is implemented through two successive processing steps that are both carried out by using the EMCF PhU technique, which has been slightly modified to take into account the Doppler centroid differences of the exploited interferometric SAR data pairs. The experimental results, achieved by applying the proposed approach to a data set consisting of European Remote Sensing (ERS) SAR data acquired from June 1992 to August 2007 over the Napoli (Italy) bay area, confirm the effectiveness of the proposed PhU approach.

Space‐borne radar interferometry techniques for the generation of deformation time series: An advanced tool for Earth's surface displacement analysis

This work is focused on advanced differential SAR interferometry (DInSAR) techniques for the generation of deformation time series from sequences of SAR images. We first present the basic rationale of these techniques providing some details of the most well known algorithms. Subsequently, through the analysis of selected case studies focused on the available C‐band SAR data archives, we show the relevance of the retrieved spatially dense deformation time series for the comprehension of several geophysical phenomena. We finally introduce, again via a real case study, the advances brought in by the new generation X‐band space‐borne SAR sensors, highlighting new investigation possibilities for fast varying deformation phenomena.

On the Extension of the Minimum Cost Flow Algorithm for Phase Unwrapping of Multitemporal Differential SAR Interferograms

In this paper, an extension of the minimum cost flow (MCF) algorithm dealing with a sparse data grid, which allows the unwrapping of multitemporal differential synthetic aperture radar (SAR) interferograms for the generation of deformation time series, is presented. The proposed approach exploits both the spatial characteristics and the temporal relationships among multiple interferograms relevant to a properly chosen sequence. In particular, the presented solution involves two main steps: first of all, for each arc connecting neighboring pixels on the interferometric azimuth/range grid, the unwrapped phase gradients are estimated via the MCF technique applied in the temporal/perpendicular baseline plane. Following this step, these estimates are used as a starting point for the spatial-unwrapping operation implemented again via the MCF approach but carried out in the azimuth/range plane. The presented results, achieved on simulated and real European Remote Sensing satellite SAR data, confirm the effectiveness of the extended MCF unwrapping algorithm