Phase unwrapping and integration of finite differences are key problems in several technical fields, among which is synthetic aperture radar (SAR) interferometry. In this paper, we propose a general ...formulation for robust and efficient integration of finite differences and for phase unwrapping, which includes standard techniques (e.g., minimum cost flow and least squares phase unwrapping) as subcases. The proposed approach allows obtaining more reliable and accurate solutions by exploiting redundant differential estimates (not only between nearest neighboring points) and multidimensional information (e.g., multitemporal). In addition, a model of the signal (e.g., multibaseline or multifrequency) or external data (e.g., GPS or leveling measurements) can be integrated. The method requires the solution of linear or quadratic programming problems, for which computationally efficient algorithms exist. The validation tests performed on real and simulated SAR data confirm the validity of the method, which was integrated in our production chain and successfully used also in massive productions.
Persistent scatterer interferometry is a widely used technique to detect and monitor slow terrain movements, with millimetric accuracy, from satellite synthetic aperture radar (SAR) data. We have ...recently proposed a method, named persistent scatterer pair (PSP), aimed at overcoming some limitations of standard techniques. The PSP method is characterized by the fact of exploiting only the relative properties of neighboring pairs of points for both detection and analysis of persistent scatterers (PSs), intended in the general sense of scatterers that exhibit interferometric coherence for the time period and baseline span of the acquisitions, including both point-like and distributed scatterers. Thanks to the pair-of-point approach, the PSP technique is intrinsically not affected by artifacts slowly variable in space, like those depending on atmosphere or orbits. Moreover, by exploiting a very redundant set of pair-of-point connections, the PSP approach guarantees extremely dense and accurate displacement and elevation measurements, both in correspondence of structures and when the backscattering is weak or distributed as in the case of natural terrains. In all cases, the measurements keep the full resolution of the input SAR images. In this work, the qualifying characteristics of the PSP technique are described, and several application examples and validation tests based on COSMO-SkyMed data are reported, which demonstrate the validity of the proposed approach.
Kinematic characterization of a landslide at large, small, and detailed scale is today still rare and challenging, especially for long periods, due to the difficulty in implementing demanding ground ...surveys with adequate spatiotemporal coverage. In this work, the suitability of space-borne synthetic aperture radar sub-pixel offset tracking for the long-term monitoring of the Slumgullion landslide in Colorado (US) is investigated. This landslide is classified as a debris slide and has so far been monitored through ground surveys and, more recently, airborne remote sensing, while satellite images are scarcely exploited. The peculiarity of this landslide is that it is subject to displacements of several meters per year. Therefore, it cannot be monitored with traditional synthetic aperture radar differential interferometry, as this technique has limitations related to the loss of interferometric coherence and to the maximum observable displacement gradient/rate. In order to overcome these limitations, space-borne synthetic aperture radar sub-pixel offset tracking is applied to pairs of images acquired with a time span of one year between August 2011 and August 2013. The obtained results are compared with those available in the literature, both at landslide scale, retrieved through field surveys, and at point scale, using airborne synthetic aperture radar imaging and GPS. The comparison showed full congruence with the past literature. A consistency check covering the full observation period is also implemented to confirm the reliability of the technique, which results in a cheap and effective methodology for the long-term monitoring of large landslide-induced movements.
Interferometric processing of series of data acquired over time by synthetic aperture radar (SAR) satellites makes it possible to measure millimetric deformations (typically due to landslides, ...subsidence and earthquake or volcanic phenomena) and to monitor the stability of terrain and infrastructures. Despite the unique capability to observe very large areas, this technology has been typically applied to the analysis of relatively small sites or specific geophysical phenomena. In this work, we present the first application of this technology to a national scale project, which required the processing, through advanced persistent scatterer interferometry (PSI) techniques, of about 20,000 SAR images acquired from 1992 to 2014 over the whole Italian territory by the ERS, Envisat, and COSMO-SkyMed satellites. The obtained results provide a huge database of surface deformation measurements covering the last 20years and the whole Italian territory, which represents an extremely useful and pioneering service for geo-hazards mapping and prevention.
•Millimetric ground deformation measurements over the whole Italian territory from 1992 to 2014•Monitoring of landslides, subsidence, and earthquake or volcanic phenomena•First application of synthetic aperture radar interferometry to national scale•Processed about 20,000 images acquired by ERS, Envisat and COSMO-SkyMed satellites
We measured ground displacements before and after the 2009 L'Aquila earthquake using multi-temporal InSAR techniques to identify seismic precursor signals. We estimated the ground deformation and its ...temporal evolution by exploiting a large dataset of SAR imagery that spans seventy-two months before and sixteen months after the mainshock. These satellite data show that up to 15 mm of subsidence occurred beginning three years before the mainshock. This deformation occurred within two Quaternary basins that are located close to the epicentral area and are filled with sediments hosting multi-layer aquifers. After the earthquake, the same basins experienced up to 12 mm of uplift over approximately nine months. Before the earthquake, the rocks at depth dilated, and fractures opened. Consequently, fluids migrated into the dilated volume, thereby lowering the groundwater table in the carbonate hydrostructures and in the hydrologically connected multi-layer aquifers within the basins. This process caused the elastic consolidation of the fine-grained sediments within the basins, resulting in the detected subsidence. After the earthquake, the fractures closed, and the deep fluids were squeezed out. The pre-seismic ground displacements were then recovered because the groundwater table rose and natural recharge of the shallow multi-layer aquifers occurred, which caused the observed uplift.
The present work focuses on the postseismic deformation observed in the region of L'Aquila (central Italy) following the Mw 6.3 earthquake that occurred on April 6, 2009. A new, 16-month-long dataset ...of COSMO-SkyMed SAR images was analysed using the Persistent Scatterer Pairs interferometric technique. The analysis revealed the existence of postseismic ground subsidence in the mountainous rocky area of Mt Ocre ridge, contiguous to the sedimentary plain that experienced coseismic subsidence. The postseismic subsidence was characterized by displacements of 10 to 35 mm along the SAR line of sight. In the Mt Ocre ridge, widespread morphological elements associated with gravitational spreading have been previously mapped. We tested the hypothesis that the postseismic subsidence of the Mt Ocre ridge compensates the loss of equilibrium induced by the nearby coseismic subsidence. Therefore, we simulated the coseismic and postseismic displacement fields via the finite element method. We included the gravitational load and fault slip and accounted for the geometrical and rheological characteristics of the area. We found that the elastoplastic behaviour of the material under gravitational loading best explains the observed postseismic displacement. These findings emphasize the role of gravity in the postseismic processes at the fault scale.
The present work focuses on the postseismic deformation observed in the region of L’Aquila (central Italy) following the Mw 6.3 earthquake that occurred on April 6, 2009. A new, 16-month-long dataset ...of COSMO-SkyMed SAR images was analysed using the Persistent Scatterer Pairs interferometric technique. The analysis revealed the existence of postseismic ground subsidence in the mountainous rocky area of Mt Ocre ridge, contiguous to the sedimentary plain that experienced coseismic subsidence. The postseismic subsidence was characterized by displacements of 10 to 35 mm along the SAR line of sight. In the Mt Ocre ridge, widespread morphological elements associated with gravitational spreading have been previously mapped. We tested the hypothesis that the postseismic subsidence of the Mt Ocre ridge compensates the loss of equilibrium induced by the nearby coseismic subsidence. Therefore, we simulated the coseismic and postseismic displacement fields via the finite element method. We included the gravitational load and fault slip and accounted for the geometrical and rheological characteristics of the area. We found that the elastoplastic behaviour of the material under gravitational loading best explains the observed postseismic displacement. These findings emphasize the role of gravity in the postseismic processes at the fault scale.
The availability of several sensors with complementarity characteristics, in terms of spatial resolution and tasking flexibility, offers new opportunities for SAR interferometry applications. In this ...work, we quantitatively discuss the complementarity of C-band low-resolution sensors (Sentinel-1 and Envisat) and the X-band high-resolution COSMO-SkyMed SAR constellation, both with theoretical analyses and experiments, focusing the comparison on interferometric measurements on real targets, in terms of density, deformation precision, and 3D positioning precision. Moreover, we provide a characterization of the measurement precision as a function of the interferometric baseline values. The obtained results confirm that whereas Sentinel-1 can systematically cover very large areas, COSMO-SkyMed can guarantee a more detailed analysis, which can be fundamental for building and infrastructure monitoring.
Synthetic aperture radar (SAR) interferometry is a powerful technology for measuring slow terrain movements. The extraction of this information is a complex task, because the phase of the signal is ...measured only modulo 2pi and is affected by noise and systematic terms. The persistent scatterer (PS) approach brought important advances in the solution of this problem. In this work, we present a new method, named persistent scatterer pairs (PSP) method, for the identification and the analysis of PS in series of full resolution SAR images. The problems coming from orbital and atmosphere phase artifacts are effectively overcome by exploiting their spatial correlation, without using model based interpolations or fits, which can be advantageous when the atmospheric artifacts or the displacement to be retrieved are not very well described by the models used in the standard PS approach. Moreover, the proposed method does not need a preprocessing to calibrate the data and is insensitive to the density of PS candidates, it is able to identify PS in natural terrains and PS characterized by non linear movements, is computationally efficient and highly parallelizable. The results obtained on real ERS SAR data confirm the validity of the proposed approach.