Next Article in Journal
Changes in Glaciers and Glacial Lakes in the Bosula Mountain Range, Southeast Tibet, over the past Two Decades
Previous Article in Journal
Coherence of Eddy Kinetic Energy Variation during Eddy Life Span to Low-Frequency Ageostrophic Energy
Previous Article in Special Issue
Susceptibility Analysis of Land Subsidence along the Transmission Line in the Salt Lake Area Based on Remote Sensing Interpretation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

CAPS: A New Method for the Identification of Different Surface Displacements in Landslide and Subsidence Environments through Correlation Analysis on Persistent Scatterers Time-Series from PSI

1
Dipartimento di Scienze della Terra dell’Ambiente e della Vita (DISTAV), Università degli Studi di Genova, 16132 Genova, Italy
2
Research Institute for Geo-Hydrological Protection, National Research Council, Strada delle Cacce 73, 10135 Torino, Italy
3
Dipartimento di Fisica (DIFI), Università degli Studi di Genova, 16146 Genova, Italy
4
Sistemi Informativi Territoriali, Comune di Genova, Via di Francia 3, 16149 Genova, Italy
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(15), 3791; https://doi.org/10.3390/rs14153791
Submission received: 23 June 2022 / Revised: 2 August 2022 / Accepted: 4 August 2022 / Published: 6 August 2022
(This article belongs to the Special Issue Applications of Remote Sensing in Geological Engineering)

Abstract

:
Persistent Scatterer Interferometry (PSI) is one of the most powerful tools for identifying and monitoring areas exposed to surface deformations such as landslides or subsidence. In this work, we propose a new method that we named CAPS (Correlation Analysis on Persistent Scatterers), to extend the capability of PSI in recognizing and characterising areas influenced by complex ground deformations and differential motions. CAPS must be applied to both ascending and descending orbits separately and comprises three major steps: (i) calculating the cross-correlation matrix on detrended PS time-series; (ii) extracting PS pairs with similarity greater than a given threshold; (iii) grouping PS in families by sorting and classification. Thus, in both orbits, PS Families identify groups of PS with similar movements. This allows distinguishing sectors characterised by different displacements over time even in areas with similar LOS (Line of Sight) velocities. As test sites, we considered four different known geological scenarios: two representing landslide environments (Santo Stefano d’Aveto and Arzeno, both in Liguria, NW Italy) and two subsidence environments (Rome and Venice, urban and surrounding areas). This method proved to be versatile, applicable to different geological situations and at different scales of observation, for recognizing both regional and local differential deformations.

Graphical Abstract

1. Introduction

Ground movements, such as land subsidence and slope instabilities, represent a major threat to population, buildings and infrastructures, especially in a period of increasing frequency of drastic meteo-hydrologic phenomena [1,2,3]. Conducting a careful prevention campaign to identify and monitor risk areas is not always easy and is often very expensive with common in situ methods. In the last decades, the Interferometric Analysis of Synthetic Aperture Radar images (InSAR) has been proven to be a very important complementary tool to detect, quantify, and monitor landslides and urban subsidence, as well as deformation of volcanic edifices or seismogenic displacements [4,5,6,7]. The launch of the Sentinel-1A (2014) and Sentinel-1B (2016) satellites within the frame of the European Copernicus Programme [8] improved the capability of InSAR analysis, also thanks to the small repeat cycle of the satellites, which is 12 days each or 6 days considering both together, due to 180° shifting of their orbits.
Three InSAR techniques have been developed so far:
-
Differential Interferometry (DInSAR): ground deformation is analysed by comparing two successive SAR acquisitions;
-
Small Baseline Subset (SBAS): provides time-series of movements and velocity maps by long-term, multi-temporal analysis of pairs of SAR acquisitions with small spatial and temporal baselines;
-
Persistent Scatterer Interferometry (PSI): provides time-series of movements over natural or anthropogenic-based targets (e.g., outcrops, building, pylons), called Persistent Scatterers (PS), and/or over homogenous areas (e.g., non-cultivated lands, desert areas), called Distributed Scatterers (DS), by analysing at least 30–35 SAR acquisitions.
In this study, we developed a new method to improve the capability of PSI. PSI results are generally displayed in terms of time-series of displacement observed along the Line of Sight (LOS) of the satellite in both Ascending and Descending orbits. The angular coefficient of the line interpolating all the points of the time-series represents the mean LOS velocity (mm/y) of the PS. However, similar LOS velocities can derive from trends that could even be very different from each other. An example is reported in Figure 1 where two groups of PS with different time-series and similar LOS velocities (~0 mm/y) are shown. They likely refer to two different geologic situations that cannot be distinguished by exclusively considering their velocities. With our method, which we named CAPS (Correlation Analysis on Persistent Scatterers), we aim at a better characterisation of complex geological scenarios through an automatic identification of groups of PS (i.e., PS Families) with similar deformation trends without considering their velocity values. CAPS is based on the Waveform Similarity Analysis (WSA), a method commonly used in seismology to find similar seismic events through the calculation of a cross-correlation matrix [9,10,11,12,13,14]. After extrapolating PS for both orbits, in the R software environment we detrended and removed the offset from the time-series. Then, we calculated the cross-correlation matrix from which we selected PS pairs with similarity greater than a chosen threshold. Finally, we recognised the PS Families by a sorting and classification technique.
We applied CAPS to both orbits to four different areas: Santo Stefano d’Aveto and Arzeno landslides, and Rome and Venice, where urban and surrounding areas are subject to differential subsidence motions. Furthermore, for Rome and Venice we cut a series of subscenes, thus also applying CAPS at different scales.
Results indicate that the PS Families found with CAPS are composed of neighbour PS characterised by time-series that are similar in shape. This allowed recognition of sectors with homogeneous deformation in both landslides and subsidence environments, from larger areas to localised sectors, down to single districts or groups of edifices. The application of CAPS proved to be a powerful and versatile tool, useful to better understand the kinematics of larger as well as more localised areas where several different movements can occur together.
The method presented in this paper is an improvement of the procedure used by [15] to investigate the landslide affecting Santo Stefano d’Aveto. In [15], PS Families were identified through a visual interpretation of the cross-correlation matrix. Specifically, the cross-correlation matrix was drawn using a suitable colour scale in order to visually identify PS Families. CAPS makes the analysis of the cross-correlation results more rigorous, faster and easier to apply to both orbits and to any case study. The families are defined through a multi-step algorithm (see Section 3.2) that automatically extracts groups of PS from the cross-correlation matrix, that have time series with a similarity greater than a given threshold. The results obtained by applying CAPS to the Santo Stefano d’Aveto landslide confirm what was already found by [15]; considering the descending orbit only, three sectors in the landslide characterised by different movements are clearly recognizable, in accordance with previous interpretations of in situ data [16,17].

2. Test Sites

In this work, we selected four different geological scenarios already analysed in previous studies (citations are in the next sub-paragraphs): two represent active landslide processes (Santo Stefano d’Aveto and Arzeno villages, both in Liguria, NW Italy) and two are active subsidence environments (Rome and Venice, urban and surrounding areas). Our scope is to present and test the CAPS method; thus, we will not discuss in depth the geological implication of our results and we just report a brief setting for each test site.

2.1. Santo Stefano d’Aveto

The village of Santo Stefano d’Aveto is located in the Ligurian Apennine, near the border between Liguria and Emilia-Romagna (NW Italy). The town centre is found in a lateral branch of the Aveto Valley, at about 1000 m a.s.l.; nearby, up to 1300 m a.s.l., there are some smaller hamlets, the most important being Roncolongo and Rocca d’Aveto. The valley of Santo Stefano is ENE–WSW oriented, and it is surrounded by some of the highest peaks of the Ligurian Apennine, reaching 1700–1800 m a.s.l.
The geology of the area is complex, with marly limestone flysch and sedimentary melanges comprising sandstones, heterogeneous breccias, and ophiolitic olistoliths (basalts, serpentinites) up to km in size [18,19]. The bedrock is covered by colluvial and landslide deposits.
The main geomorphological feature is a complex earth slide/earth flow landslide which extends along the valley floor [20] (Figure 2). According to the Italian Inventory of Landslide Phenomena (IFFI) [21], the landslide is divided into an active upper sector and a dormant lower sector. Recent studies involving different monitoring methodologies, both in situ and by remote sensing, show that both sectors are active; the upper sector is moving at about 40–50 mm/y, while the estimated speed of the lower sector is 20–30 mm/y [15,21,22,23].

2.2. Arzeno

Arzeno is a small village of the municipality of Ne (Liguria, NW Italy), located in the upper Graveglia valley, in the central-eastern Ligurian Apennine. The main town stands halfway up the eastern slope of the Graveglia valley, at about 600 m a.s.l. The slope goes from about 400 m a.s.l. on the valley floor to 1127 m a.s.l. by the top of Mt. Chiappozzo.
Bedrock formations outcropping in the Arzeno area belong to the Bracco-Val Graveglia Unit. They are composed of a Jurassic ophiolitic sequence (serpentinites, gabbros, basalts and associated ophiolitic breccias) overlayed by a Jurassic-Cretaceous deep-sea sedimentary sequence comprising cherts, limestones, and shales [25,26].
From the geomorphological point of view, evidence suggests the presence of a deep-seated gravitational slope deformation (DSGSD) that involves the entire slope from the Mt. Chiappozzo ridge to the valley floor [27]. This area is also affected by several complex landslides, among which the most relevant and well-known are the Arzeno landslide and the adjacent Prato landslide. The Arzeno landslide extends from about 850 m a.s.l. to the valley floor, where it causes a deviation of the Reppia stream. IFFI [21] identifies several sectors with different activity (Figure 3): active sectors have been identified by the Arzeno settlement and in the lower part of the landslide, near the valley floor; other sectors of the landslide are classified as dormant or stabilised. The Prato landslide extends from about 800 m a.s.l. to 450 m a.s.l., almost reaching the valley floor. The IFFI classifies it as a stabilised landslide.
Recent technical and scientific reports, by using several remote sensing and in situ monitoring methodologies, have estimated average speed of about 20 mm/y for the active sectors of both Arzeno and Prato landslides, with peaks of about 30–40 mm/y in correspondence of the Arzeno settlement [28,29].

2.3. Venice

Venice is the regional capital of Veneto, in north-eastern Italy. One of the most famous cities in the world, it is built on several small, sandy islands in the central part of the Venice Lagoon, near the coast of the Adriatic Sea. The city and the surrounding areas rest on Quaternary alluvial and tidal deposits (Figure 4). The thickness of the Quaternary deposits exceeds 700 m [30].
As it is located at sea level, the city of Venice is susceptible to sea level rise, which increases the risk of flooding due to tides and adverse meteorological conditions. Relative sea level rise has been calculated to be 0.23 m in the XX century; 0.11 m can be attributed to true sea level rise caused by climate change, while the rest is due to subsidence, caused by natural and anthropogenic causes such as the groundwater pumping in past decades and construction of new buildings [31].
As a consequence, subsidence in the area of Venice and the Venice Lagoon has been extensively studied [32,33]. Vertical displacement in the city averages is of the order of 0.8–1 mm/y; the ancient, pre-1500 part of the city is stabler compared to the more recent outskirts. Parts of the city, interested by reconstructions, restorations of buildings and other urban maintenance processes, show higher displacement rates that can be up to 10 mm/y for short time intervals [33]. Different sectors of the Venice Lagoon subside at different rates: the northern lagoon has a 2–3 mm/y rate of vertical displacement, while the southern lagoon has a 3–4 mm/y rate [32].

2.4. Rome

Rome, the Italian capital, is a city with 2.8 million inhabitants, located in the Lazio region in central Italy. The city stands on the lower part of the valley of the Tiber river, which flows across the urban area from NNE to SSW. In the northern part of the city, the Aniene river, coming from the E, flows into the Tiber. On the two sides of the main river, the landscape is characterised by low hills, which represent the lower foothills of the Colli Albani on the SE, the Sabatini Mountains on the NW and the Central Apennines on the NE. Southwest of the city lies the Tiber alluvial plain, which extends to the coast of the Tyrrhenian Sea.
The geology of the area is characterised by volcanic deposits belonging to the Sabatino volcanic district and the Albano volcanic district, Pleistocene marine deposits, and Meso-Cenozoic limestones outcropping on the Apennine foothills on the NE (Figure 5). The Tiber valley and the narrower Aniene valley are filled with alluvial deposits [35].
Ground subsidence in the city of Rome is well known and many studies have dealt with its quantification using InSAR techniques [4,36]. In most cases, subsidence is caused by the weight of construction on unconsolidated sediments. Thus, it is more evident on the alluvial material of the Tiber valley [37,38]. There are some cases, such as the Acque Albule plain, where groundwater extraction has played a major role [39].
Vertical motion due to subsidence in the urban area has been estimated between 2 to 5 mm/y, with peaks of 7–9 mm/y; the higher values are found along the Tiber valley. Even higher values, up to 20 mm/y have been observed for the area of the Fiumicino Airport, which is located in the Tiber alluvial plain near the Tyrrhenian coast. The nearby Lido di Ostia coastal zone shows vertical motion up to 12 mm/y [4].

3. Materials and Methods

3.1. PS Processing: Software and Data

We used Copernicus Sentinel-1A acquisitions from 2015 to 2022 for each test site, downloading 180-to-200 images for each orbit via the Alaska Satellite Facility (ASF) web repository (Table 1). They are Single Look Complex (SLC) TOPSAR data acquired in Interferometric Wide (IW) mode with VV Polarisation. The use of only Sentinel-1A acquisitions proved to be adequate for the purpose of this work; adding the Sentinel-1B dataset would have employed an excessive amount of computing time and storage space without giving appreciable improvements on the final results.
We conducted the PS processing using the latest available versions of the following open-source software:
(1)
SNAP (SeNtinel Application Platform) [41] that allows (i) selecting the Master image, (ii) correcting and splitting both Master and Slave images, (iii) co-registering each Slave with the Master and creating the relative interferograms, and finally (iv) exporting the resulting products in StaMPS format;
(2)
snap2stamps [42] that allows performing the previous ii to iv steps automatically for the entire stack of images with a series of python scripts;
(3)
StaMPS (Stanford Method for Persistent Scatterers) [43] that processes with multiple steps the products resulting from SNAP/snap2stamps to extrapolate PS velocity maps.
A detailed description of the processing chain and the StaMPS parameterization can be found in [4,42,43]. It is important to highlight that the StaMPS parameters greatly influence the final results, also depending on the analysed geological scenario; thus, a proper StaMPS setting is of utmost importance. The parameterizations used in this work for the landslide and subsidence test sites are shown in Table 2. The choice has been made according to previous studies [15,44,45] on similar geological scenarios and by autonomously trying different settings. Below, the parameters we changed are described according to [43]:
-
max_topo_error: value in metre of the maximum uncorrelated DEM error accepted. Changed only for the subsidence cases;
-
unwrap_time_win: length in days of the time window used to filter/smooth the phase in time by estimating the noise for each pair of neighbouring pixels. Changed for both landslide and subsidence cases;
-
unwrap_grid_size: spacing of the resampling grid. Changed for both landslide and subsidence cases;
-
unwrap_gold_n_win: size of the window for the Goldstein filter [46]. Changed for both landslide and subsidence cases;
-
scla_deramp: it estimates the phase ramp for each interferogram (if set to “y”). Changed for both landslide and subsidence cases;
-
scn_time_win: window size of the low-pass temporal filter. Changed for both landslide and subsidence cases.

3.2. Waveform Similarity Analysis (WSA) and Correlation Analysis on Persistent Scatterers (CAPS)

In seismology, in order to better understand the seismic features of an area, it is useful to identify seismic events that present some degree of correlation. In particular, a group of earthquakes occurring very close in space and having a similar source-time function, propagation path, station site and recording instrument are defined as a “Family”. One of the most-used methods to identify earthquake families is the so-called Waveform Similarity Analysis, or WSA [12,47], which is based on seismogram cross-correlation.
Such a method can obviously be applied to other research fields whenever time-series need to be analysed. In this paper, we propose the application of the WSA for the identification of groups of PS showing similar time-series of ground movement and, therefore, to define areas characterised by similar deformation patterns at the surface.
In particular, the similarity between two PS time-series, a1(t) and a2(t), is estimated through the normalised cross-correlation function (NCCF) defined as (1):
N C C F   τ = C 12 τ C 11 0 C 22 0
where τ is the time shift, and C 12 τ is defined by (2):
C 12 τ = + a 1 t a 2 t + τ d t
The maximum value assumed by NCCF (i.e., the cross-correlation coefficient) is indicative of the similarity between the signals. Therefore, the more two PS, a1(t) and a2(t), show similar time-series, the more the function NCCF will approach 1; 1 is the value assumed by the cross-correlation coefficient when signals are clones, −1 when signals are uncorrelated and 0 when signals are completely different.
After calculation of the cross-correlation coefficient for each pair of PS in a set, it is possible to define a cross-correlation matrix that can be visual inspected to define groups of PS having a similarity greater than a predefined threshold (Figure 6).
In order to identify PS Families within the cross-correlation matrix in a more objective way, we have developed a two-step algorithm. First, all pairs of PS showing a similarity greater than a certain threshold are identified and extracted from the matrix, then they are grouped in families through a sorting and classifying technique.
In the first step, the most important parameter that must be defined is the cross-correlation threshold, or index of similarity. The choice of the threshold value is crucial and is essentially based on subjective expert judgement. In many seismological applications [11,13,14], several decreasing values of the threshold are iteratively applied until the best compromise between the number of families and the population of each family is achieved. In this work we tuned such threshold following the previous rule and also considering the characteristics of the investigated PS dataset and the target of the analysis. Specifically, when dealing with small areas and a limited number of PS (e.g., <2000) and the goal is to identify significant heterogeneities in small-scale ground displacements, the threshold should be chosen considering low values (e.g., <0.90). Figure 7a shows an example of time-series of three PS with similarity greater than 0.83. If the analysis is performed on extended areas with a large number of PS (e.g., >10,000) and the goal is to identify slightly different deformations at larger scale, the threshold should be increased to values higher than the previous case (e.g., >0.90). Figure 7b shows an example of PS time-series with similarity greater than 0.95. In other words, the choice of the cross-correlation threshold can be used to guide the analysis towards different levels of resolution. In the next paragraph, we show some examples to explain the influence of the threshold value on the results.
In the second step, the algorithm used for defining the families is based on the Equivalence Class approach [48], already tested and applied to earthquake data sets by many authors [11,49,50].
Given the list of PS pairs with a level of similarity greater than the threshold, it is necessary to group elements that are in the same equivalence class of “sameness”. The adopted algorithm works like a “bridging technique”; if two pairs of PS (A, B) and (B, C) share a common PS (B) then all three PS are attributed to the same family even if the similarity between A and C is lower than the selected threshold value. Event B therefore represents the “bridge” between these pairs.
In summary, the proposed CAPS method is based on the following operative steps:
  • Removal of the offset (offset removal) and of the linear trend (detrending) from the PS time-series using the “detrend” function of the package “pracma” of R;
  • Calculation of the cross-correlation matrix through the “cor” function of the package “stats” of R;
  • Selection of the PS pairs with a level of similarity greater than a certain cross-correlation threshold;
  • Definition of PS Families through the “eclass” subroutine in Fortran 77 [48].
The outcomes resulting from the WSA are strongly dependent on the number of data points considered, which in our application corresponds to the number of SAR images processed. In general, for such analyses the more data points/images there are, the more reliable the results. However, it is not a simple task to objectively define the minimum number of data points in the time-series that guarantees the reliability of the results. To this end, the minimum number of images required for reliable PSI analysis, which conventionally consists of at least 30–35 SAR images, can be roughly assumed as the lower limit. For each case study described in the paper, we consider approximately 180–200 images in each orbit (Table 1) and, therefore, the correlation analysis is performed on time-series with 180–200 points.

4. Results

4.1. Santo Stefano d’Aveto

The Santo Stefano d’Aveto landslide is characterised by LOS velocities ranging between 5 and 25 mm/y in ascending mode and between −5 and −45 mm/y in descending mode. In both orbits, absolute velocity values decrease from north-east to south-west with a maximum in correspondence of the village of Rocca d’Aveto (Figure 8a,e).
We applied CAPS to ascending and descending PS by using various threshold values, from 0.83 to 0.95. As specified in the previous paragraph, there are no objective parameters to be used to choose the best threshold a priori. The choice must be made considering the best compromise between number of lost PS, families distribution and most important, geologically realistic output. For the Santo Stefano d’Aveto landslide, results show that in both orbits the more the threshold increases, the more the landslide is subdivided in multiple Families (Figure 8). In particular, by using a threshold of 0.83, a high number of PS is preserved, and the landslide results subdivided into two families for both orbits. However, in the central-southwestern area, PS within and outside the perimeter of the landslide are grouped together, meaning that this threshold is too low to highlight differences. On the other hand, with 0.95, the number of PS lost is too high and the subdivision of the landslide appears to be excessive. Thus, the best values seem to be 0.91 for the ascending mode and 0.87 for the descending mode. Both maintain a sufficient number of PS (higher in descending) and subdivide similarly the landslide body by identifying four and three families, respectively. The number and percentage of remaining PS together with the number of identified families by varying the threshold values are reported in Appendix A for each test site.
Figure 9 shows the time-series of a reference PS for each family identified within the landslide with the chosen thresholds. To better highlight the differences in trend the velocity information (i.e., the inclination of the time-series) is also shown. It can be noted that the trend of the time-series of each family effectively differs from the others and, in particular, these differences are more evident considering the descending orbit. This can be related to the western orientation of the landslide that makes the descending orbit the best to analyse the area. Analysing the distribution of the families within the landslide perimeter for the ascending and descending modes, the correspondence between families SA1 and SD1, SA2 and SD2 and SA3 + SA4 and SD3 appears clear (Figure 9). This distribution suggests a subdivision of the landslide into several sectors characterized by different kinematics and dynamics likely related to a complex internal structure of the landslide body and to the geological and geomorphological characteristics of the area. This confirms what was previously proposed in [15] based on the descending orbit only.

4.2. Arzeno

The landslide affecting the village of Arzeno, considering the perimeter from IFFI [21], is characterised by LOS velocity ranging between 2.5 and 15 mm/y in ascending mode and between −10 and −25 mm/y in descending mode (Figure 10). The highest values of displacement along the LOS for both orbits are observed in the topographically higher part of the Arzeno village whilst the western-lower part is affected by slower movements. The perimeter made by IFFI excludes from the landslide the village of Prato (see Figure 3) to the north. However, our results show that also this area is affected by important values of displacements observed along the LOS. These results may suggest that the village of Prato should also be included within the area of the landslide.
We analysed these data with CAPS by using threshold values from 0.83 to 0.92 in ascending mode and from 0.75 to 0.87 in descending mode. By increasing the threshold, the area of the landslide is not divided into an increasing number of families such as in the case of Santo Stefano d’Aveto. Again considering the village of Prato, the maximum number of families here identified is two in ascending (0.83 and 0.87) and three in descending (0.82); however, due to the high loss of PS (i.e., PS not grouped into families), high threshold values appear not meaningful. In this case, results that describe a similar scenario in both orbits, with a good number of PS preserved, are obtained with ascending 0.87 and descending 0.75. However, by comparing the time-series of the three PS families obtained with descending 0.82, it is possible to note important differences between family ArD2 and family ArD3 (Figure 11). Thus, the subdivision into three families seems realistic and, in addition, reflects the clear change in LOS velocity highlighted in Figure 10a,e.
For this case study, the most appropriate threshold values are ascending 0.87 and descending 0.82. This last, in particular, despite preserving a low number of PS, especially in the topographically lower part of Arzeno, highlights an interesting subdivision of the landslide. Similarly to the Santo Stefano d’Aveto test site, this subdivision can be referred to relevant heterogeneities in structure and geotechnical properties of the landslide body. The site of Arzeno is currently the subject of numerous studies and in situ analyses, not yet published. Our results from CAPS could therefore provide useful constraints to gather better knowledge on the kinematics of the landslide.

4.3. Venice

The third case study concerns the Venice Lagoon and surrounding areas, where we conducted two analyses. The first one (Figure 12) is a large scale overview that considers Venice and neighbouring islands, Lido and the continental area of Mestre. By applying CAPS, we aimed at distinguishing macro-families that identify sector characterised by similar movements or by high heterogeneity. However, considering large areas implies the analysis of a huge amount of PS. In this case, for instance, we extrapolated more than 180,000 PS for each orbit, and this would imply calculating two cross-correlation matrices made of >32 billion of elements (180,000 × 180,000). This requires extreme computational resources that are not always accessible. In order to reduce the necessary resources and calculation time, we randomly selected 1 PS every 10 reducing the total number of PS to ~18,000 and the number of elements of the matrices to ~320 million. This is an acceptable approximation that does not undermine the reliability of the results.
Clearly, more detailed analysis on smaller areas requires maintenance of the highest possible number of PS. This is highlighted by the results of our second analysis focused only on the main group of islands of Venice and aimed at distinguishing deformations related to single buildings or groups of edifices. In this case, we maintained all the PS extrapolated that are ~20,000 for each orbit.
Results of the application of CAPS at large overview reported in Figure 12 show that for both orbits the more the threshold is increased, the more the continental area to the west is subject to multiple subdivisions and to a growing loss of PS. To the east, Venice and Lido (see Figure 4) are included in a single family with all the thresholds considered with the exception of descending 0.75 that groups together both the western and eastern areas. Thus, in this case, ascending 0.90 and descending 0.86 can be considered the best threshold values to identify macro-areas characterised by similar movements. In particular, with these thresholds, it is possible to identify the north-westernmost part of the study area, i.e., the farthest from the lagoon, as a single macro-area whereas the anthropized area of Mestre (see Figure 4) as characterised by multiple families. These results could help to consider and plan further analysis on the area of Mestre in order to clarify the nature of the identified families such as, for instance, relationships with different behaviours of the substrate related to different types of infrastructures and foundations.
The results of the application of CAPS to the main group of islands of Venice (Figure 13) highlights in both orbits the presence of one large family that includes almost all of the PS. A series of smaller families are dispersed within the main one and referred to single buildings or groups of edifices as shown in Figure 13b,e. For this case, we present only the most appropriate threshold values that are ascending 0.93 and descending 0.88. Figure 13b,e show a further close-up on the eastern sector of Venice where five PS families in ascending mode and four PS families in descending mode have been identified. In Figure 13c,f the time-series of families VeA1, VeA2, VeA3 and VeD1, VeD2, VeD3, VeD4 are shown. The differences in trend between the main family (VeA1/VeD1) and the other families are evident. In particular, in ascending mode the three time-series differ starting from October 2017. This evidence could be compared to the local geophysical or mechanical properties in order to identify the mechanisms controlling the movement in correspondence of the edifices highlighted with CAPS (e.g., foundation instabilities, structural distress). The family distribution can thus help to plan in-depth analyses of specific targets.

4.4. Rome

The analysis at large overview of the metropolitan area of Rome is shown in Figure 14. Even in this case, we conducted two analyses: in the first one we focused on the entire Metropolitan area of Rome considering 1 PS every 65 (in both orbits: total initial PS >1,300,000; remaining PS ~20,000); the second one is focused on Tivoli and surrounding areas considering 1 PS every 5 in ascending mode (Total initial PS >101,000; remaining PS ~20,000) and 1 PS every 6 in descending mode (Total initial PS >115,000; remaining PS ~19,000).
PS maps (Figure 14a,c) show the presence of important subsidence deformations mainly to the west, in the area of Ponte Galeria, and to the east, in the area of Tivoli (see Figure 5). Here, velocities observed along the LOS reach values of −7 to −9 mm/y.
Similarly to the other test sites, PS families identified with CAPS (Figure 14b,d) do not reflect the velocity distribution. Areas characterised by higher velocities are included in macro-families that refer to different behaviours of the surface likely related to different environments. For instance, the green family to the east, very similar in both orbits, could be related to a coastal environment influenced by the delta of the Tiber River. The urban area of the city of Rome is poorly represented probably because of the high presence of heterogeneities. This suggests that, as in the case of Mestre, highly anthropized areas should be analysed with a more local-scale overview.
The analysis of the area of Tivoli through the application of CAPS (Figure 15) shows that the PS in the central part are included in a single family (blue).
This PS family does not reflect the distribution of the velocities suggesting that subsidence could interest an area larger compared to that could be derived from the only velocity distribution. In particular, by overlaying the central family to the lithologic map (Figure 15b,d) it is possible to observe that most of the PS belonging to the central blue family (~60%) are in correspondence of the Travertines (grey) and the Eluvial-colluvial deposits (light green). This agrees with the analyses conducted by [51,52] on the characteristics of the Acque Albule aquifer in Tivoli with different in situ methods. In addition, the blue family also includes PS from the eastern valley where the Aniene River flows from east to west arising new question marks about the dynamics of the area and more specifically, about the geometry and the extension of the aquifer.
As stated before, the aim of this paper is not to explain the causes related to each single family. However, by referring to [4], we can argue that in the area of Rome the subdivision into families is controlled by the variation in geological and hydrogeological characteristics and, more specifically, in the groundwater table geometry and geotechnical properties of the subsoil.

5. Discussion

The method proposed in this paper, CAPS (Correlation Analysis on Persistent Scatterers), proved to be a powerful tool for processing and analysing PSI results. Through a correlation analysis on PS time-series, CAPS allows the identification of groups (or families) of PS showing similar movement over time. The association of a PS to a family is independent of its velocity but only depends on the observed movement in ascending or descending orbit. Therefore, CAPS is also useful to highlight differences in surface deformation in areas that are characterized by similar velocities.
In the examples proposed in the present study for testing the method, CAPS not only allows us to confirm the results of previous studies, but also to provide new findings. For the Santo Stefano d’Aveto and Arzeno landslides, the geographic distribution of the PS families highlights different sectors affected by different kinematics, thus confirming what pointed out by several authors [12,17,24]. For Venice and Rome, the PS family distribution shows macro-areas characterized by similar subsidence movement, but also small areas characterized by small deformation differences (as those highlighted in Venice main island and Tivoli area).
The advantages of the proposed method are twofold. The first benefit is versatility. CAPS can be applied to different geological scenarios. In this paper we have applied it to landslides, which are characterized by significant variation of superficial movement and by complex geological–geomorphological setting, and to subsidence environments, where the deformation rate is lower and less heterogeneous than the previous cases. Moreover, CAPS can work at both large and small scale of observation depending on the purpose of the analysis and/or on the required spatial resolution. It can be used to identify similar environments in terms of ground surface deformation in very large areas (as in Roma and Venice Lagoon test sites), in small areas (as in the main group of islands of Venice, Santo Stefano d’Aveto and Arzeno test sites), or to explain anomalous movements of buildings and infrastructures.
The second benefit is adaptability. As already stated, the identification of PS Families is mainly driven by the value of the cross-correlation threshold used to quantify the similarity between data. On the one hand, such value can be appropriately chosen as a function of the scope of the work and the characteristics of the investigated area. On the other hand, it can be also selected in order to satisfy a priori information. For instance, the value of the threshold can be defined in order to reproduce or confirm findings of previous studies or to point out the deformation in correspondence of specific targets subjected to known movements.
In Appendix A, a comparison table to highlight similarities and differences between each test site is reported.

6. Conclusions

In this research, we proposed a new method, called CAPS (Correlation Analysis on Persistent Scatterers), to analyse time-series of displacement derived from PSI. CAPS is a tool for post-processing ascending and descending PS that groups them into families as a function of the level of similarity of the times series of displacements regardless of the LOS velocity values. This makes it possible to distinguish areas characterized by similar LOS velocity but different ground displacement over time, providing useful constraints for relating the deformations observed at the surface with the geological, geophysical or structural properties of the study area. To this aim, we combined PSI results from both ascending and descending orbits separately with WSA, a procedure generally used in seismology for investigating seismicity patterns. The distribution of the PS families provides information about the kinematics of a targeted area, allowing the identification of sectors with similar displacement. Of note, CAPS can be applied to ascending and descending PS data separately and the comparison between the resulting families can also provide useful information about intrinsic errors and uncertainties in PSI results (due, for instance, to the relative orientation of the satellite respect to the direction of the movement).
The most important parameter that has to be chosen when using CAPS is the value of the cross-correlation threshold, which specifies the level of similarity between two PS time-series. The choice of that value depends on:
-
the focus of the study (e.g., highlighting structural anomalies of single edifices, heterogeneities within a landslide body or different environments in extended areas);
-
the characteristics of the investigated area or phenomena (e.g., orientation compared to the LOS; geological, geotechnical or geophysical heterogeneities that drive kinematics and dynamics of the area);
-
the scale of the investigated area or phenomena (e.g., local or regional)
-
the characteristics of the PS dataset (e.g., number and distribution of PS).
In the application presented in the foregoing, we have provided useful suggestions to choose the value of the threshold as a function of the characteristics of the investigated area or phenomena. When dealing with landslides characterized by high spatially variable deformation rates (>15 mm/y), which are generally related to highly heterogeneous geomorphology as well as complex kinematics and dynamic features, the best threshold value should be sought for the interval 0.85–0.90. Since landslides often develop in little anthropized areas, where few numbers of PS can be extrapolated, the use of relatively low thresholds allows maintaining a sufficient number of PS to be grouped in families. When dealing with areas characterized by low spatial heterogeneities and low displacements such as subsidence environments (absolute LOS velocity, generally, lower than 10 mm/y) the most appropriate threshold should be chosen in the interval 0.90–0.95 or even higher. These are often spatially extended areas where the ground movements are generally related to geological and hydrogeological properties, variations in the groundwater table, and geotechnical characteristics of the subsoil. In these cases, a very large number of PS can be extrapolated and, even using very high threshold values (e.g., 0.98), a sufficient number of PS is maintained at the end of the analysis, thus making possible the identification of subtle differences. Finally, when focusing attention on single buildings, infrastructures or small areas where the deformation is related to local geotechnical proprieties of soil foundations and/or structural distress, the thresholds should be set around 0.90. Obviously, a more helpful guide to setting the most appropriate value of the cross-correlation threshold may be compiled in the future by applying the CAPS tool to many other case studies (also considering the application of CAPS to active volcanoes or glaciers). This will also be helpful to better relate the identified PS families to the geological and geophysical properties that control the ground deformation.
Another benefit of this method is that, when grouping PS into a limited number of families, the analysis of the movement observed at each PS become much simpler. For each family, the analyst can select a reference PS whose time-series can be analysed and compared to the others in order to highlight changes over time in the deformation rate of each homogeneous sector.
To summarise, CAPS has been proven to be a powerful tool for gathering a better interpretation of the deformation maps derived from PS observations. The identified PS families can be related to the differences in terms of ground displacement in areas subject to movement, such as landslides or subsidence, and can be used for identifying the geological and geophysical properties that lead to that deformation. We have shown that the number of families depends on the value of the cross-correlation threshold, whose setting is crucial. Even if the suggestions provided in this study can be helpful for future applications, a trial-and-error procedure is still the best way to determine its value. When increasing the threshold, the total number of PS grouped into families decreases as well as the number of PS in each family, while the number of families tends to increase. Hence, the “best” value should lead to a compromise between the number of families and their population in terms of number of PS. In addition, when geological and geophysical information of the area under study is available, the spatial distribution of the PS families obtained by varying the value of the cross-correlation threshold should be compared to such information. This comparison can thus provide important constraints to identify homogeneous sectors characterised by similar geological and geophysical properties (which determine their displacement) such as depth of the water table, morphology, subsoil geotechnical properties, and geo-lithological composition.
Obviously, CAPS does not allow us to overcome the intrinsic uncertainties of the PSI method (and of InSAR in general) related to the viewing geometry. This remains a crucial parameter that can have a huge impact on the PSI results. The more the real ground movements are perpendicular to the LOS, the more they are underestimated. The viewing geometry also affects the final value of the cross-correlation threshold chosen for the ascending or the descending orbit. As an example, in our landslide test sites (that are better viewed in descending mode due to their orientation), a higher threshold value for the ascending orbit was adopted in order to achieve a distribution of PS families similar to that obtained with the descending orbit.

Author Contributions

Conceptualization, E.B. and G.F.; methodology, E.B., G.F., M.T. and A.F.; software and data curation, G.F., E.B. and D.S.; validation, L.C., P.C., S.B., S.T. and F.F.; writing—original draft preparation, E.B., G.F. and A.F.; writing—review and editing by all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Copernicus Sentinel-1A data 2015–2022 retrieved from ASF DAAC. Available online: https://asf.alaska.edu/ (accessed on 20 June 2022), processed by ESA.

Acknowledgments

The authors would like to acknowledge the open-source software and packages SNAP, StaMPS, snap2stamps, QGIS 3.22 and StaMPS_Visualizer (https://doi.org/10.5281/zenodo.4407188) used for data processing and for visualization and interpretation of results and figures. The Step Forum, available online: https://forum.step.esa.int (accessed on 20 June 2022), a place to share problems and ideas, is also acknowledged. Anonymous reviewers are acknowledged as well, as their suggestions certainly improved this work.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Number of PS, percentage of remaining PS and number of families by varying the threshold are reported for each test site.
Table A1. Number of PS, percentage of remaining PS and number of families by varying the threshold are reported for each test site.
Santo Stefano d’Aveto
PS tot0.830.870.910.95
Ascending324310 PS (96%)
5 Fam.
278 PS (86%)
6 Fam.
208 PS (64%)
6 Fam.
114 PS (35%)
9 Fam.
Descending384357 PS (93%)
5 Fam.
343 PS (89%)
9 Fam.
288 PS (75%)
12 Fam.
138 PS (36%)
8 Fam.
Arzeno
PS tot0.750.820.870.92
Ascending168141 PS (84%)
2 Fam.
117 PS (70%)
4 Fam.
90 PS (54%)
6 Fam.
37 PS (22%)
3 Fam.
Descending12898 PS (77%)
5 Fam.
80 PS (63%)
6 Fam.
38 PS (30%)
3 Fam.
25 PS (20%)
2 Fam.
Venice
PS tot0.800.860.900.95
Ascending1863716183 PS (87%)
11 Fam.
13116 PS (70%)
21 Fam.
990 PS (53%)
45 Fam.
3973 PS (21%)
67 Fam.
Descending182079458 PS (52%)
17 Fam.
5659 PS (31%)
34 Fam.
3116 PS (17%)
48 Fam.
661 PS (4%)
28 Fam.
Venice zoom
PS tot0.880.900.930.95
Ascending20798-19735 PS (95%)
22 Fam.
18056 PS (87%)
22 Fam.
15189 PS (73%)
43 Fam.
Descending1958212004 PS (61%)
43 Fam.
9644 PS (49%)
51 Fam.
5487 PS (28%)
110 Fam.
2896 PS (15%)
133 Fam.
Rome
PS tot0.910.930.950.96
Ascending2003818335 PS (92%)
6 Fam.
16773 PS (84%)
27 Fam.
13803 PS (69%)
111 Fam.
11459 PS (57%)
181 Fam.
Descending1990914486 PS (73%)
7 Fam.
12282 PS (62%)
37 Fam.
9184 PS (46%)
68 Fam.
7287 PS (37%)
98 Fam.
Rome Tivoli
PS tot0.940.960.980.99
Ascending2035119775 PS (97%)
3 Fam.
18577 PS (91%)
5 Fam.
13177 PS (65%)
47 Fam.
6690 PS (33%)
137 Fam.
Descending1924516981 PS (88%)
3 Fam.
14091 PS (73%)
5 Fam.
7591 PS (39%)
37 Fam.
2959 PS (15%)
78 Fam.

References

  1. Guzzetti, F.; Peruccacci, S.; Rossi, M.; Stark, C.P. Rainfall thresholds for the initiation of landslides in central and southern Europe. Meteorol. Atmos. Phys. 2007, 98, 239–267. [Google Scholar] [CrossRef]
  2. Kalubowila, P.; Lokupitiya, E.; Halwatura, D.; Jayathissa, G. Threshold rainfall ranges for landslide occurrence in Matara district of Sri Lanka and findings on community emergency preparedness. Int. J. Disaster Risk Reduct. 2021, 52, 101944. [Google Scholar] [CrossRef]
  3. Casagli, N.; Tofani, V.; Sassa, K.; Bobrowsky, P.T.; Takara, K. (Eds.) Understanding and Reducing Landslide Disaster Risk: Volume 3 Monitoring and Early Warning; Springer Nature: Berlin, Germany, 2020. [Google Scholar] [CrossRef]
  4. Blasco, J.M.D.; Foumelis, M.; Stewart, C.; Hooper, A. Measuring urban subsidence in the Rome metropolitan area (Italy) with sentinel-1 SNAP-StaMPS Persistent Scatterer Interferometry. Remote Sens. 2019, 11, 129. [Google Scholar] [CrossRef] [Green Version]
  5. Hooper, A.; Zebker, H.; Segall, P.; Kampes, B. A new method for measuring deformation on volcanoes and other natural terrains using InSAR persistent scatterers. Geophys. Res. Lett. 2004, 31, 1–5. [Google Scholar] [CrossRef]
  6. Bürgmann, R.; Rosen, P.A.; Fielding, E.J. Synthetic aperture radar interferometry to measure Earth’s surface topography and its deformation. Annu. Rev. Earth Planet. Sci. 2000, 28, 169–209. [Google Scholar] [CrossRef]
  7. Orellana, F.; Hormazábal, J.; Montalva, G.; Moreno, M. Measuring Coastal Subsidence after Recent Earthquakes in Chile Central Using SAR Interferometry and GNSS Data. Remote Sens. 2022, 14, 1611. [Google Scholar] [CrossRef]
  8. Copernicus: Europe’s Eye on Earth. Available online: https://www.copernicus.eu (accessed on 20 June 2022).
  9. Shearer, P.M. Evidence from a cluster of small earthquakes for a fault at 18 km depth beneath Oak Ridge, Southern California. Bull. Seismol. Soc. Am. 1998, 88, 1327–1336. [Google Scholar] [CrossRef]
  10. Astiz, L.; Shearer, P.M.; Agnew, D.C. Precise relocations and stress change calculations for the Upland earthquake sequence in Southern California. J. Geophys. Res. 2000, 105, 2937–2953. [Google Scholar] [CrossRef]
  11. Ferretti, G.; Massa, M.; Solarino, S. An improved method for the identification of seismic families: Application to the Garfagnana–Lunigiana area (Italy). Bull. Seismol. Soc. Am. 2005, 95, 1903–1915. [Google Scholar] [CrossRef] [Green Version]
  12. Barani, S.; Ferretti, G.; Massa, M.; Spallarossa, D. The waveform similarity approach to identify dependent events in instrumental seismic catalogues. Geophys. J. Int. 2007, 168, 100–108. [Google Scholar] [CrossRef] [Green Version]
  13. Barani, S.; Ferretti, G.; Scafidi, D.; Spallarossa, D. Analysis of seismicity and micro-seismicity associated with the October–November 2010 Sampeye swarm, Southwestern Alps. Tectonophysics 2014, 611, 130–140. [Google Scholar] [CrossRef]
  14. Cattaneo, M.; Augliera, P.; Spallarossa, D.; Lanza, V. A waveform similarity approach to investigate seismicity patterns. Nat. Hazards 1999, 19, 123–138. [Google Scholar] [CrossRef]
  15. Balbi, E.; Terrone, M.; Faccini, F.; Scafidi, D.; Barani, S.; Tosi, S.; Crispini, L.; Cianfarra, P.; Poggi, F.; Ferretti, G. Persistent Scatterer Interferometry and Statistical Analysis of Time-Series for Landslide Monitoring: Application to Santo Stefano d’Aveto (Liguria, NW Italy). Remote Sens. 2021, 13, 3348. [Google Scholar] [CrossRef]
  16. Tofani, V.; Raspini, F.; Catani, F.; Casagli, N. Persistent Scatterer Interferometry (PSI) technique for landslide characterization and monitoring. Remote Sens. 2013, 5, 1045–1065. [Google Scholar] [CrossRef] [Green Version]
  17. Intrieri, E.; Frodella, W.; Raspini, F.; Bardi, F.; Tofani, V. Using satellite interferometry to infer landslide sliding surface depth and geometry. Remote Sens. 2020, 12, 1462. [Google Scholar] [CrossRef]
  18. Marino, M.; Terranova, R. Ophiolithifer Complexes of Mt.Aiona and Mt.Penna and their relation with sedimentary series (Ligurian-Parmisan Appenine). Boll. Soc. Geol. It. 1980, 99, 183–203. [Google Scholar]
  19. Elter, P.; Lasagna, S.; Marroni, M.; Pandolfi, L.; Vescovi, P.; Zanzucchi, G. Note Illustrative della Carta Geologica d’Italia alla Scala 1:50.000, Foglio 215 Bedonia; Regione Emilia Romagna; S.EL.CA. srl: Firenze, Italy, 2005; 116p. [Google Scholar]
  20. Faccini, F.; Robbiano, A.; Roccati, A. Geomorphological characterization, site investigation and management of the S Stefano d’Aveto Landslide (Ligurian Apennine Italy). In Proceedings of the 33rd International Geological Society Congress, GSM-01 General Contributions to Geomorphology, Oslo, Norway, 6–14 August 2008. [Google Scholar]
  21. Progetto IFFI (Inventario dei Fenomeni Franosi in Italia) [Inventory of Landslide Phenomena in Italy]. Available online: https://www.isprambiente.gov.it/it/progetti/cartella-progetti-in-corso/suolo-e-territorio-1/iffi-inventario-dei-fenomeni-franosi-in-italia (accessed on 17 June 2022).
  22. ARPAL, Remover Project 2008–2017. Available online: http://www.vincolimap.it/img/Remover/Commenti_Siti/GE013_commento_tot.pdf (accessed on 17 June 2022).
  23. Bottero, B.; Poggi, F. The integration of different monitoring techniques for the analysis of slow moving landlside hazard: S. Stefano d’Aveto case study (Genoa-Liguria). Geoing. E Attività Estrattiva 2014, 1, 5–18. [Google Scholar]
  24. Regione Liguria. Litologia. Available online: https://srvcarto.regione.liguria.it/geoviewer2/pages/apps/geoportale/index.html?id=1907 (accessed on 17 June 2022).
  25. Bortolotti, V.; Principi, G. The Bargonasco-Upper Val Graveglia ophiolitic succession, Northern Apennines, Italy. Ofioliti 2003, 28, 137–140. [Google Scholar]
  26. Bortolotti, V.; Mannori, G.; Piccini, I.; Aiello, W.; Chiari, M.; Principi, G.; Reale, V. Note Illustrative della Carta Geologica d’Italia alla Scala 1:50.000, Foglio 232 Sestri Levante; S.EL.CA. srl: Firenze, Italy, 2015; 117p. [Google Scholar]
  27. Brandolini, P.; Canepa, G.; Faccini, F.; Robbiano, A.; Terranova, R. Geomorphological and geo-environmental features of the Graveglia Valley (Ligurian Apennines, Italy). Geogr. Fis. Dinam. Quat. 2007, 30, 99–116. [Google Scholar]
  28. Allasia, P.; Baldo, M.; Faccini, F.; Godone, D.; Notti, D.; Poggi, F. The Role of Measure of Deep-Seated Displacements in the Monitoring Networks on Large-Scale Landslide. In Workshop on World Landslide Forum; Springer: Cham, Switzerland, 2020; pp. 49–57. [Google Scholar]
  29. Paliaga, G.; Donadio, C.; Bernardi, M.; Faccini, F. High-Resolution Lightning Detection and Possible Relationship with Rainfall Events over the Central Mediterranean Area. Remote Sens. 2019, 11, 1601. [Google Scholar] [CrossRef] [Green Version]
  30. Tosi, L.; Rizzetto, F.; Bonardi, M.; Donnici, S.; Serandrei-Barbero, R.; Toffoletto, F. Note Illustrative della Carta Geologica d’Italia alla Scala 1:50.000, Foglio 128 Venezia; Regione Veneto: Venezia, Italy, 2007; 164p. [Google Scholar]
  31. Gatto, P.; Carbognin, L. The lagoon of Venice: Natural environmental trend and man-induced modification. Hydrol. Sci. Bull. 1981, 26, 379–391. [Google Scholar] [CrossRef] [Green Version]
  32. Bock, Y.; Wdowinski, S.; Ferretti, A.; Novali, F.; Fumagalli, A. Recent subsidence of the Venice Lagoon from continuous GPS and interferometric synthetic aperture radar. Geochem. Geophys. Geosyst. 2012, 13, 1–13. [Google Scholar] [CrossRef]
  33. Tosi, L.; Teatini, P.; Strozzi, T. Natural versus anthropogenic subsidence of Venice. Sci. Rep. 2013, 3, 2710. [Google Scholar] [CrossRef] [Green Version]
  34. Regione Veneto. Database delle Diverse Litologie che Compongono il Territorio della Regione Veneto Scala 1:250.000. Available online: https://idt2.regione.veneto.it/idt/webgis/viewer?previewLayerId=474 (accessed on 15 June 2022).
  35. Funiciello, R.; Giordano, G. Note Illustrative della Carta Geologica d’Italia alla Scala 1:50.000, Foglio 374 Roma; S.EL.CA. srl: Firenze, Italy, 2008; 164p. [Google Scholar]
  36. Cigna, F.; Lasaponara, R.; Masini, N.; Milillo, P.; Tapete, D. Persistent scatterer interferometry processing of COSMO-skymed stripmap HIMAGE time series to depict deformation of the historic centre of Rome, Italy. Remote Sens. 2014, 6, 12593–12618. [Google Scholar] [CrossRef] [Green Version]
  37. Stramondo, S.; Bozzano, F.; Marra, F.; Wegmuller, U.; Cinti, F.R.; Moro, M.; Saroli, M. Subsidence induced by urbanisation in the city of Rome detected by advanced InSAR technique and geotechnical investigations. Remote Sens. Environ. 2008, 112, 3160–3172. [Google Scholar] [CrossRef]
  38. Bozzano, F.; Esposito, C.; Mazzanti, P.; Patti, M.; Scancella, S. Imaging Multi-Age Construction Settlement Behaviour by Advanced SAR Interferometry. Remote Sens. 2018, 10, 1137. [Google Scholar] [CrossRef] [Green Version]
  39. Bozzano, F.; Esposito, C.; Franchi, S.; Mazzanti, P.; Perissin, D.; Rocca, A.; Romano, E. Analysis of a Subsidence Process by Integrating Geological and Hydrogeological Modelling with Satellite InSAR Data. In Engineering Geology for Society and Territory; Springer: New York, NY, USA, 2015; Volume 5, pp. 155–159. [Google Scholar]
  40. Città Metropolitana di Roma Capitale. Carta Geologica e Litologica. Available online: https://g3w-suite.cittametropolitanaroma.it/it/map/carta-geologica-litologica/ (accessed on 15 June 2022).
  41. Veci, L.; Lu, J.; Prats-Iraola, P.; Scheiber, R.; Collard, F.; Fomferra, N.; Engdahl, M. The Sentinel-1 Toolbox. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Quebec City, QC, Canada, 13–18 July 2014; pp. 1–3. [Google Scholar]
  42. Foumelis, M.; Blasco, J.M.D.; Desnos, Y.-L.; Engdahl, M.; Fernandez, D.; Veci, L.; Lu, J.; Wong, C. ESA SNAP-StaMPS integrated processing for sentinel-1 persistent scatterer interferometry. In Proceedings of the 2018 IEEE International Geoscience and Remote Sensing Symposium, Valencia, Spain, 22–27 July 2018; pp. 1364–1367. [Google Scholar] [CrossRef]
  43. Hooper, A.; Bekaert, D.; Spaans, K.; Arıkan, M. Recent advances in SAR interferometry time-series analysis for measuring crustal deformation. Tectonophysics 2012, 514–517, 1–13. [Google Scholar] [CrossRef]
  44. Höser, T. Analysing the Capabilities and Limitations of InSAR Using Sentinel-1 Data for Landslide Detection and Monitoring. Master’s Thesis, Department of Geography, University of Bonn, Bonn, Germany, 2018. [Google Scholar]
  45. Shimosato, K.; Ukita, N. Multi-Modal Data Fusion for Land-Subsidence Image Improvement in PSInSAR Analysis. IEEE Access 2021, 9, 141970–141980. [Google Scholar] [CrossRef]
  46. Goldstein, R.M.; Werner, C.L. Radar interferogram filtering for geophysical applications. Geophys. Res. Lett. 1998, 25, 4035–4038. [Google Scholar] [CrossRef] [Green Version]
  47. Maurer, H.; Deichmann, N. Microearthquake cluster detection based on waveform similarities, with an application to the western Swiss Alps. Geophys. J. Int. 1995, 123, 588–600. [Google Scholar] [CrossRef] [Green Version]
  48. Press, W.H.; Teukolsky, S.A.; Wetterling, W.T.; Flannery, P.B. Numerical Recipes in in Fortran 77. In The Art of Scientific Computing, Second Edition Volume 1 of Fortran Numerical Recipes; Cambridge University Press: Cambridge, UK, 1992; ISBN 0-521-43064-X. [Google Scholar]
  49. Aster, R.C.; Scott, J. Comprehensive characterization of waveform similarity in microearthquake data set. Bull. Seismol. Soc. Am. 1993, 83, 1307–1314. [Google Scholar] [CrossRef]
  50. Massa, M.; Ferretti, G.; Spallarossa, D.; Eva, C. Improving automatic location procedure by waveform similarity analysis: An application in the South Western Alps (Italy). Phys. Earth Planet. Inter. 2006, 154, 18–29. [Google Scholar] [CrossRef]
  51. Billi, A.; De Filippis, L.; Poncia, P.P.; Sella, P.; Faccenna, C. Hidden sinkholes and karst cavities in the travertine plateau of a highly-populated geothermal seismic territory (Tivoli, central Italy). Geomorphology 2016, 255, 63–80. [Google Scholar] [CrossRef]
  52. De Filippis, L.; Anzalone, E.; Billi, A.; Faccenna, C.; Poncia, P.P.; Sella, P. The origin and growth of a recently-active fissure ridge travertine over a seismic fault, Tivoli, Italy. Geomorphology 2013, 195, 13–26. [Google Scholar] [CrossRef]
Figure 1. Example of six different PS with similar LOS velocities of ~0 mm/y but time-series in (a) strongly differ from those in (b).
Figure 1. Example of six different PS with similar LOS velocities of ~0 mm/y but time-series in (a) strongly differ from those in (b).
Remotesensing 14 03791 g001
Figure 2. Geological sketch map of the area of the Santo Stefano d’Aveto landslide (red dot in the upper image), modified from [24]: (1a) Active landslide; (1b) Inactive landslide; (1c) Undifferentiated landslide deposits; (2) Other quaternary deposits: slope deposits, alluvial deposits, lacustrine deposits; (3) Shales; (4) Marly limestone flysch; (5) Sandstones; (6) Polygenic breccias; (7) Basalts; (8) Serpentinites; (9) Perimeter of the Santo Stefano landslide; (10) Buildings. Landslide activity according to IFFI [21].
Figure 2. Geological sketch map of the area of the Santo Stefano d’Aveto landslide (red dot in the upper image), modified from [24]: (1a) Active landslide; (1b) Inactive landslide; (1c) Undifferentiated landslide deposits; (2) Other quaternary deposits: slope deposits, alluvial deposits, lacustrine deposits; (3) Shales; (4) Marly limestone flysch; (5) Sandstones; (6) Polygenic breccias; (7) Basalts; (8) Serpentinites; (9) Perimeter of the Santo Stefano landslide; (10) Buildings. Landslide activity according to IFFI [21].
Remotesensing 14 03791 g002
Figure 3. Geological sketch map of the area of the Arzeno landslide (red dot in the upper image), modified from [24]: (1a) Active landslide; (1b) Inactive landslide; (2) Other quaternary deposits: slope, alluvial, lacustrine; (3) Shales with limestones; (4) Limestones; (5) Cherts; (6) Basalts; (7) Ophicalcites and ophiolitic breccias; (8) Serpentinites; (9) Perimeter of the Arzeno landslide; (10) Buildings. Landslide activity according to IFFI [21].
Figure 3. Geological sketch map of the area of the Arzeno landslide (red dot in the upper image), modified from [24]: (1a) Active landslide; (1b) Inactive landslide; (2) Other quaternary deposits: slope, alluvial, lacustrine; (3) Shales with limestones; (4) Limestones; (5) Cherts; (6) Basalts; (7) Ophicalcites and ophiolitic breccias; (8) Serpentinites; (9) Perimeter of the Arzeno landslide; (10) Buildings. Landslide activity according to IFFI [21].
Remotesensing 14 03791 g003
Figure 4. Geological sketch map of the study area of Venice Lagoon (red dot in the upper image): (1) Alluvial deposits; (2) Tidal deposits; (3) Littoral deposits; (4) Venice Lagoon and Adriatic Sea; (5) Buildings. Modified from [34].
Figure 4. Geological sketch map of the study area of Venice Lagoon (red dot in the upper image): (1) Alluvial deposits; (2) Tidal deposits; (3) Littoral deposits; (4) Venice Lagoon and Adriatic Sea; (5) Buildings. Modified from [34].
Remotesensing 14 03791 g004
Figure 5. Geological sketch map of the Urban area of Rome and surroundings (red dot in the upper image), modified from [40]. (1) Slope deposits; (2) Eluvial-colluvial deposits; (3) Alluvial and lacustrine deposits; (4) Coastal deposits; (5) Travertines; (6) Volcanic deposits; (7) Pleistocene marine deposits; (8) Shales, marls, siltites; (9) Limestones; (10) Buildings.
Figure 5. Geological sketch map of the Urban area of Rome and surroundings (red dot in the upper image), modified from [40]. (1) Slope deposits; (2) Eluvial-colluvial deposits; (3) Alluvial and lacustrine deposits; (4) Coastal deposits; (5) Travertines; (6) Volcanic deposits; (7) Pleistocene marine deposits; (8) Shales, marls, siltites; (9) Limestones; (10) Buildings.
Remotesensing 14 03791 g005
Figure 6. Cross-correlation matrix for ascending (left) and descending (right) orbits of the Santo Stefano d’Aveto case study. The X and Y axes show the number of PS, the colour scale indicates the similarity between each pair (−1 maximum uncorrelation; 0 different signals; 1 maximum correlation). In a first visual interpretation, if a threshold of 0.8 is considered, each orange-to-red square would represent a single PS Family. A variable number of PS Families can be recognised in descending and ascending orbits due to the different PS extrapolated during the processing and, most importantly, to the orientation of the landslide that, in this case, is better visualised by the descending orbit.
Figure 6. Cross-correlation matrix for ascending (left) and descending (right) orbits of the Santo Stefano d’Aveto case study. The X and Y axes show the number of PS, the colour scale indicates the similarity between each pair (−1 maximum uncorrelation; 0 different signals; 1 maximum correlation). In a first visual interpretation, if a threshold of 0.8 is considered, each orange-to-red square would represent a single PS Family. A variable number of PS Families can be recognised in descending and ascending orbits due to the different PS extrapolated during the processing and, most importantly, to the orientation of the landslide that, in this case, is better visualised by the descending orbit.
Remotesensing 14 03791 g006
Figure 7. Example of time-series grouped in the same family with CAPS using two different thresholds: (a) 0.83 and (b) 0.95. Time-series in (a) are less similar to each other than those in (b) due to the lower threshold.
Figure 7. Example of time-series grouped in the same family with CAPS using two different thresholds: (a) 0.83 and (b) 0.95. Time-series in (a) are less similar to each other than those in (b) due to the lower threshold.
Remotesensing 14 03791 g007
Figure 8. The Santo Stefano d’Aveto landslide, perimeter in black line according to IFFI [21], with: PS distribution in (a) ascending mode and (e) descending mode (colour code indicates different LOS velocity ranges), and their subdivision into families with CAPS method by using different threshold values for (bd) ascending and (fh) descending (different colours indicate different PS families).
Figure 8. The Santo Stefano d’Aveto landslide, perimeter in black line according to IFFI [21], with: PS distribution in (a) ascending mode and (e) descending mode (colour code indicates different LOS velocity ranges), and their subdivision into families with CAPS method by using different threshold values for (bd) ascending and (fh) descending (different colours indicate different PS families).
Remotesensing 14 03791 g008
Figure 9. In (a,c) the chosen thresholds for ascending and descending orbit, respectively, are reported for the Santo Stefano d’Aveto test site. In (b,d) the time-series of a reference PS for each family identified with CAPS in both orbits are shown.
Figure 9. In (a,c) the chosen thresholds for ascending and descending orbit, respectively, are reported for the Santo Stefano d’Aveto test site. In (b,d) the time-series of a reference PS for each family identified with CAPS in both orbits are shown.
Remotesensing 14 03791 g009
Figure 10. The Arzeno landslide, perimeter in black line according to IFFI [21], with: PS distribution in (a) ascending mode and (e) descending mode (colour code indicates different LOS velocity ranges), and their subdivision into families with CAPS method by using different threshold values for (bd) ascending and (fh) descending (different colours indicate different PS Families).
Figure 10. The Arzeno landslide, perimeter in black line according to IFFI [21], with: PS distribution in (a) ascending mode and (e) descending mode (colour code indicates different LOS velocity ranges), and their subdivision into families with CAPS method by using different threshold values for (bd) ascending and (fh) descending (different colours indicate different PS Families).
Remotesensing 14 03791 g010
Figure 11. In (a,c) the chosen thresholds for ascending and descending orbits, respectively, are reported for the Arzeno test site. In (b,d) the time-series of a reference PS for each family identified with CAPS in both orbits are shown.
Figure 11. In (a,c) the chosen thresholds for ascending and descending orbits, respectively, are reported for the Arzeno test site. In (b,d) the time-series of a reference PS for each family identified with CAPS in both orbits are shown.
Remotesensing 14 03791 g011
Figure 12. The Venice Lagoon and surrounding areas with: PS distribution in (a) ascending mode and (e) descending mode (colour code indicates different LOS velocity ranges), and their subdivision in families with CAPS method by using different threshold values for (bd) ascending and (fh) descending (different colours indicate different PS Families).
Figure 12. The Venice Lagoon and surrounding areas with: PS distribution in (a) ascending mode and (e) descending mode (colour code indicates different LOS velocity ranges), and their subdivision in families with CAPS method by using different threshold values for (bd) ascending and (fh) descending (different colours indicate different PS Families).
Remotesensing 14 03791 g012
Figure 13. In (a,d) the subdivision of families of Venice at smaller scale. In (b,e) a further zoom showing the main family and a series of smaller ones corresponding to single or small groups of edifices. In (c,f) the time-series of the identified families within the close-up are compared.
Figure 13. In (a,d) the subdivision of families of Venice at smaller scale. In (b,e) a further zoom showing the main family and a series of smaller ones corresponding to single or small groups of edifices. In (c,f) the time-series of the identified families within the close-up are compared.
Remotesensing 14 03791 g013
Figure 14. The metropolitan area of Rome with: PS distribution in (a) ascending mode and (c) descending mode (colour code indicates different LOS velocity ranges), and their subdivision into families with CAPS method by using as threshold 0.95 in (b) ascending and (d) descending (different colours indicate different PS families).
Figure 14. The metropolitan area of Rome with: PS distribution in (a) ascending mode and (c) descending mode (colour code indicates different LOS velocity ranges), and their subdivision into families with CAPS method by using as threshold 0.95 in (b) ascending and (d) descending (different colours indicate different PS families).
Remotesensing 14 03791 g014
Figure 15. The area of Tivoli to the east of Rome. In (a,c) the PS distribution in ascending and descending mode, respectively. In (b,d) the distribution of the identified PS Families overlaid to the lithologic chart [40].
Figure 15. The area of Tivoli to the east of Rome. In (a,c) the PS distribution in ascending and descending mode, respectively. In (b,d) the distribution of the identified PS Families overlaid to the lithologic chart [40].
Remotesensing 14 03791 g015
Table 1. Information about the Sentinel-1A acquisitions used for each test site.
Table 1. Information about the Sentinel-1A acquisitions used for each test site.
TrackFrameFirst ImageLast ImageImagesMasterSub-Swath
Santo Stefano d’Aveto
Ascending15139, 141-2-3-4-523 March 20152 February 202220517 August 2019IW2
Descending168441-2-3-4-5-6-722 March 20151 February 202220325 March 2019IW3
Arzeno
Ascending15139, 141-2-3-4-518 November 20152 February 20221847 April 2019IW2
Descending66441-2-3-4-5-610 November 20156 February 202218917 January 2019IW1
Venice
Ascending117142-3, 145-6-730 March 20159 February 202220428 October 2018IW3
Descending95438-9, 440, 442-329 March 20158 February 202219721 September 2018IW2
Rome
Ascending117130-1, 133, 135-630 March 20159 February 202220428 October 2018IW2
Descending22450-1-2-3-4-5-6-724 March 201522 January 202219221 April 2018IW3
Table 2. List of StaMPS parameters used in this research. The Default and the Used values together with the involved StaMPS steps are shown. Note that max_topo_error has been changed only for the two subsidence test sites.
Table 2. List of StaMPS parameters used in this research. The Default and the Used values together with the involved StaMPS steps are shown. Note that max_topo_error has been changed only for the two subsidence test sites.
ParametersStaMPS StepDefaultUsedEnvironment
max_topo_error220 m10 mSubsidence
unwrap_time_win6730 days24 daysBoth
unwrap_grid_size6200 m10 mBoth
unwrap_gold_n_win6328Both
scla_deramp7nyBoth
scn_time_win8365 days50 daysBoth
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Balbi, E.; Ferretti, G.; Ferrando, A.; Faccini, F.; Crispini, L.; Cianfarra, P.; Scafidi, D.; Barani, S.; Tosi, S.; Terrone, M. CAPS: A New Method for the Identification of Different Surface Displacements in Landslide and Subsidence Environments through Correlation Analysis on Persistent Scatterers Time-Series from PSI. Remote Sens. 2022, 14, 3791. https://doi.org/10.3390/rs14153791

AMA Style

Balbi E, Ferretti G, Ferrando A, Faccini F, Crispini L, Cianfarra P, Scafidi D, Barani S, Tosi S, Terrone M. CAPS: A New Method for the Identification of Different Surface Displacements in Landslide and Subsidence Environments through Correlation Analysis on Persistent Scatterers Time-Series from PSI. Remote Sensing. 2022; 14(15):3791. https://doi.org/10.3390/rs14153791

Chicago/Turabian Style

Balbi, Evandro, Gabriele Ferretti, Andrea Ferrando, Francesco Faccini, Laura Crispini, Paola Cianfarra, Davide Scafidi, Simone Barani, Silvano Tosi, and Martino Terrone. 2022. "CAPS: A New Method for the Identification of Different Surface Displacements in Landslide and Subsidence Environments through Correlation Analysis on Persistent Scatterers Time-Series from PSI" Remote Sensing 14, no. 15: 3791. https://doi.org/10.3390/rs14153791

APA Style

Balbi, E., Ferretti, G., Ferrando, A., Faccini, F., Crispini, L., Cianfarra, P., Scafidi, D., Barani, S., Tosi, S., & Terrone, M. (2022). CAPS: A New Method for the Identification of Different Surface Displacements in Landslide and Subsidence Environments through Correlation Analysis on Persistent Scatterers Time-Series from PSI. Remote Sensing, 14(15), 3791. https://doi.org/10.3390/rs14153791

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop