Permanent Scatterer SAR Interferometry (PSInSAR) aims to identify coherent radar targets exhibiting high phase stability over the entire observation time period. These targets often correspond to ...point-wise, man-made objects widely available over a city, but less present in non-urban areas. To overcome the limits of PSInSAR, analysis of interferometric data-stacks should aim at extracting geophysical parameters not only from point-wise deterministic objects (i.e., PS), but also from distributed scatterers (DS). Rather than developing hybrid processing chains where two or more algorithms are applied to the same data-stack, and results are then combined, in this paper we introduce a new approach, SqueeSAR, to jointly process PS and DS, taking into account their different statistical behavior. As it will be shown, PS and DS can be jointly processed without the need for significant changes to the traditional PSInSAR processing chain and without the need to unwrap hundreds of interferograms, provided that the coherence matrix associated with each DS is properly "squeezed" to provide a vector of optimum (wrapped) phase values. Results on real SAR data, acquired over an Alpine area, challenging for any InSAR analysis, confirm the effectiveness of this new approach.
We present the continuous monitoring of ground deformation at regional scale using ESA (European Space Agency) Sentinel-1constellation of satellites. We discuss this operational monitoring service ...through the case study of the Tuscany Region (Central Italy), selected due to its peculiar geological setting prone to ground instability phenomena. We set up a systematic processing chain of Sentinel-1 acquisitions to create continuously updated ground deformation data to mark the transition from static satellite analysis, based on the analysis of archive images, to dynamic monitoring of ground displacement. Displacement time series, systematically updated with the most recent available Sentinel-1 acquisition, are analysed to identify anomalous points (i.e., points where a change in the dynamic of motion is occurring). The presence of a cluster of persistent anomalies affecting elements at risk determines a significant level of risk, with the necessity of further analysis. Here, we show that the Sentinel-1 constellation can be used for continuous and systematic tracking of ground deformation phenomena at the regional scale. Our results demonstrate how satellite data, acquired with short revisiting times and promptly processed, can contribute to the detection of changes in ground deformation patterns and can act as a key information layer for risk mitigation.
Permanent scatterer InSAR (PSInSAR™) provides a new high‐resolution methodology for detecting and precisely measuring long‐term and seasonal aquifer‐system response to pumping and recharge. In ...contrast to conventional InSAR, the permanent scatterer methodology utilizes coherent radar phase data from thousands of individual radar reflectors on the ground to develop displacement time series and to produce velocity field maps that depict aquifer‐system response with a high degree of spatial detail. In this study, we present the first results of a prototype study in Las Vegas Valley, Nevada, that demonstrate how this methodology can be utilized in heavily pumped groundwater basins to analyze aquifer‐system response to long‐term and seasonal pumping. We have developed a series of velocity field maps of the valley for the 1992–1996, 1996–2000, and 2003–2005 time periods that show that despite rising water levels associated with an artificial recharge program, long‐term, residual, inelastic aquifer‐system compaction (subsidence) is continuing in several parts of the valley. In other areas, however, long‐term subsidence has been arrested and locally reversed. The seasonal, elastic responses to alternating pumping and recharge cycles were segregated from the long‐term trends and analyzed for spatial and temporal patterns. The results show oscillations in which the maximum seasonal responses are associated with the late stages of the annual artificial recharge cycles, and that similar seasonal subsidence signals are related to summer pumping cycles. The differentiation of the seasonal response through the use of time series data further allows the estimation of elastic and inelastic skeletal storage coefficients, providing a basis for future work that could characterize the storage properties of an aquifer system with a high degree of spatial resolution.
This paper presents the results of a blind experiment that is performed using two pairs of dihedral reflectors. The aim of the experiment was to demonstrate that interferometric synthetic aperture ...radar (InSAR) measurements can indeed allow a displacement time series estimation with submillimeter accuracy (both in horizontal and vertical directions), provided that the data are properly processed and the impact of in situ as well as atmospheric effects is minimized. One pair of dihedral reflectors was moved a few millimeters between SAR acquisitions, in the vertical and east-west (EW) directions, and the ground truth was compared with the InSAR data. The experiment was designed to allow a multiplatform and multigeometry analysis, i.e., each reflector was carefully pointed in order to be visible in both Envisat and Radarsat acquisitions. Moreover, two pairs of reflectors were used to allow the combination of data gathered along ascending and descending orbits. The standard deviation of the error is 0.75 mm in the vertical direction and 0.58 mm in the horizontal (EW) direction. GPS data were also collected during this experiment in order to cross-check the SAR results
High-resolution interferometric synthetic aperture radar (InSAR) permanent scatterer data allow us to resolve the rates and variations in the rates of slow-moving landslides. Satellite-to-ground ...distances (range changes) on landslides increase at rates of 5 to 7 millimeters per year, indicating average downslope sliding velocities from 27 to 38 millimeters per year. Time-series analysis shows that displacement occurs mainly during the high-precipitation season; during the 1997-1998 El Niño event, rates of range change increased to as much as 11 millimeters per year. The observed nonlinear relationship of creep and precipitation rates suggests that increased pore fluid pressures within the shallow subsurface may initiate and accelerate these features. Changes in the slope of a hill resulting from increases in the pore pressure and lithostatic stress gradients may then lead to landslides.
This work describes a new procedure aimed to semi-automatically identify clusters of active persistent scatterers and preliminarily associate them with different potential types of deformational ...processes over wide areas. This procedure consists of three main modules: (i) ADAfinder, aimed at the detection of Active Deformation Areas (ADA) using Persistent Scatterer Interferometry (PSI) data; (ii) LOS2HV, focused on the decomposition of Line Of Sight (LOS) displacements from ascending and descending PSI datasets into vertical and east-west components; iii) ADAclassifier, that semi-automatically categorizes each ADA into potential deformational processes using the outputs derived from (i) and (ii), as well as ancillary external information. The proposed procedure enables infrastructures management authorities to identify, classify, monitor and categorize the most critical deformations measured by PSI techniques in order to provide the capacity for implementing prevention and mitigation actions over wide areas against geological threats. Zeri, Campiglia Marittima–Suvereto and Abbadia San Salvatore (Tuscany, central Italy) are used as case studies for illustrating the developed methodology. Three PSI datasets derived from the Sentinel-1 constellation have been used, jointly with the geological map of Italy (scale 1:50,000), the updated Italian landslide and land subsidence maps (scale 1:25,000), a 25 m grid Digital Elevation Model, and a cadastral vector map (scale 1:5000). The application to these cases of the proposed workflow demonstrates its capability to quickly process wide areas in very short times and a high compatibility with Geographical Information System (GIS) environments for data visualization and representation. The derived products are of key interest for infrastructures and land management as well as decision-making at a regional scale.
Deformation in the northern San Francisco Bay area is dominated by a series of sub‐parallel strike‐slip faults. Existing GPS observations provide some constraint on the slip rates of these faults, ...however these have only limited resolution for resolving shallow fault behavior, such as brittle creep. We use a 30 image Permanent Scatterer InSAR (PS‐InSAR) dataset spanning the time interval 1992–2001 to dramatically increase the density of surface deformation observations. We find a discontinuity in observed surface velocities across the Rodgers Creek fault, around Santa Rosa and further north, consistent with shallow creep at rates of up to 6 mm/yr. The creeping segments are located in areas of local transtension, suggesting that lowered normal stresses may play a role in the distribution of creep. The existence of creep could significantly reduce expected moment release in future earthquakes on the Rodgers Creek fault, and thus has implications for seismic hazard assessment.
On 13 January 2012, the Italian vessel, Costa Concordia, wrecked offshore Giglio Island, along the coast of Tuscany (Italy). The ship partially sunk, lying on the starboard side on a 22 degree steep ...rocky seabed, making the stability conditions of the ship critically in danger of sliding, shifting and settling. The tilted position of the ship created also pernicious conditions for the divers involved in the search and rescue operations. It became immediately clear that a continuous monitoring of the position and movements of the ship was of paramount importance to guarantee the security of the people working around and within the wreck. Starting from January 19, the Italian constellation of synthetic aperture radar (SAR) satellites, COSMO-SkyMed (CSK), was tasked to acquire high resolution images of the wreck. Thanks to CSK's short response and revisiting time and its capability to acquire high resolution images in Spotlight mode, satellite data were integrated within the real time, ground-based monitoring system implemented to provide the civil protection authorities with a regular update on the ship stability. Exploitation of both the phase (satellite radar interferometry, InSAR) and amplitude (speckle tracking) information from CSK images, taken along the acquisition orbit, Enhanced Spotlight (ES)-29, revealed a general movement of the translation of the vessel, consistent with sliding toward the east of the hull on the seabed. A total displacement, with respect to the coastline, of 1666 mm and 345 mm of the bow and stern, respectively, was recorded, over the time period of 19 January-23 March 2012.
On May 20 and 29, 2012, two earthquakes of magnitudes 5.9 and 5.8 (Mw), respectively, and their aftershock sequences hit the central Po Plain (Italy), about 40 km north of Bologna. More than 2,000 ...sizable aftershocks were recorded by the Isti-tuto Nazionale di Geofisica e Vulcanologia (INGV; National Institute of Geophysics and Volcanology) National Seismic Network (http://iside.rm.ingv.it/). The sequence was generated by pure compressional faulting over blind thrusts of the western Ferrara Arc, and it involved a 50-km-long stretch of this buried outer front of the northern Apennines. The focal mechanisms of the larger shocks agree with available structural data and with present-day tectonic stress indicators, which show locally a maximum horizontal stress oriented ca. N-S; i.e. oriented perpendicular to the main structural trends. Most of the sequence occurred between 1 km and 12 km in depth, above the local basal detachment of the outer thrust fronts of the northern Apennines. We measured the surface displacement patterns associated with the mainshocks and some of the larger aftershocks (some of which had Mw >5.0) by applying the Interferometric Synthetic Aperture Radar (InSAR) technique to a pair of C-Band Radarsat-1 images. We then used the coseismic motions detected over the epicentral region as input information, to obtain the best-fit model fault for the two largest shocks. …