Next Article in Journal
Earth Observation Contribution to Cultural Heritage Disaster Risk Management: Case Study of Eastern Mediterranean Open Air Archaeological Monuments and Sites
Next Article in Special Issue
A Case Study of Novel Landslide Activity Recognition Using ALOS-1 InSAR within the Ragged Mountain Western Hillslope in Gunnison County, Colorado, USA
Previous Article in Journal
Exploring the Combined Effect of Urbanization and Climate Variability on Urban Vegetation: A Multi-Perspective Study Based on More than 3000 Cities in China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unraveling Spatial and Temporal Heterogeneities of Very Slow Rock-Slope Deformations with Targeted DInSAR Analyses

1
Department of Earth and Environmental Sciences, University of Milano-Bicocca, 20126 Milano, Italy
2
Department of Earth Sciences, ETH Zurich, 8092 Zurich, Switzerland
3
Centro di Monitoraggio Geologico, ARPA Lombardia, 23100 Sondrio, Italy
*
Author to whom correspondence should be addressed.
Submission received: 9 March 2020 / Revised: 9 April 2020 / Accepted: 17 April 2020 / Published: 22 April 2020
(This article belongs to the Special Issue Leveraging on SAR Imagery for Landslide Detection and Monitoring)

Abstract

:
Spaceborne radar interferometry is a powerful tool to characterize landslides at local and regional scales. However, its application to very slow rock slope deformations in alpine environments (displacement rates < 5 cm/year) remains challenging, mainly due to low signal to noise ratio, atmospheric disturbances, snow cover effects, and complexities resulting from heterogeneous displacement in space and time. Here we combine SqueeSARTM data, targeted multi-temporal baseline DInSAR, GPS data, and detailed field morpho-structural mapping, to unravel the kinematics, internal segmentation, and style of activity of the Mt. Mater deep-seated gravitational slope deformation (DSGSD) in Valle Spluga (Italy). We retrieve slope kinematics by performing 2D decomposition (2D InSAR) of SqueeSARTM products derived from Sentinel-1 data acquired in ascending and descending orbits. To achieve a spatially-distributed characterization of DSGSD displacement patterns and activity, we process Sentinel-1 A/B images (2016-2019) with increasing temporal baselines (ranging from 24-days to 1-year) and generate several multi-temporal interferograms. Unwrapped displacement maps are validated using ground-based GPS data. Interferograms derived with different temporal baselines reveal a strong kinematic and morpho-structural heterogeneity and outline nested rockslides and active sectors, that arise from the background displacement signal of the main DSGSD. Seasonal interferograms, supported by GPS displacement measurements, reveal non-linear displacement trends suggesting a complex response of different slope sectors to rainfall and snowmelt. Our analyses clearly outline a composite slope instability with different nested sectors possibly undergoing different evolutionary trends towards failure. The results herein outline the potential of a targeted use of DInSAR for the detailed investigation of very slow rock slope deformations in different geological and geomorphological settings.

Graphical Abstract

1. Introduction

Giant, very slow rock slope deformations known as deep-seated gravitational slope deformations (DSGSD) are widespread in high mountain environments worldwide [1,2,3,4,5]. The evolution of these phenomena is strongly constrained by inherited geological structure and long-term climatic forcing and occurs through progressive damage accumulation. This can lead to progressive strain localization and nucleation of secondary landslides, nested within the main DSGSD mass, that become faster and more sensitive to hydrological processes as damage accumulates [6,7,8].
DSGSDs creep slowly over thousands of years, affecting entire valley slopes and displacing rock volumes up to billions of cubic meters [4,9,10]. Therefore, even if displacement rates are low (usually smaller than 5 cm/year [11,12,13]) and apparently steady, DSGSDs result in large cumulative deformations that are mirrored by different morpho-structural field evidence. These include extensional features as double-crested ridges, trenches, scarps, counterscarps, and half-grabens, usually dominating in the upper slope sectors, and compressional features like bulging, thrusting and folding, and nested large landslides in the middle-lower sectors. The kinematic interpretation of individual features and their associations provide critical information on the overall slope deformation mechanism and the patterns of both distributed and localized strain, reflecting deep-seated deformation mechanisms [14].
Slow rock slope deformations cause severe degradation of rock mass strength, damage to infrastructures [15,16,17] and possible evolution of slope sectors towards faster rates [6]. A reliable characterization of DSGSD mechanisms and activity is thus crucial to identify possible failure scenarios and for risk reduction. Nevertheless, DSGSDs are usually characterized by: (a) strongly heterogeneous displacement patterns, associated with complex mechanisms (strain partitioning, damage localization and/or secondary landslides nested at different depths); (b) variable and often unknown trends of activity and sensitivity to external forcing [4,6,8]. In order to investigate these points, displacements must be characterized in a spatially-distributed fashion and with a sufficiently high rate of temporal sampling.
Among the available in situ and remote landslide mapping and monitoring techniques, spaceborne differential radar interferometry (DInSAR) allows characterizing ground deformation rates from a few millimeters to centimeters per year [18], maximizing the spatial and temporal coverage at a relatively low cost [12,19]. Persistent-scatterer interferometry (PSI) techniques allow extracting coherent information on line-of-sight (LOS) displacement rates at point targets, even when spatial baselines and revisit times of satellite platforms were not favorable [20]. Thus, PSI became a standard tool for landslide investigation, also applied to a regional appraisal of DSGSDs [12,13]. However, the application of point-like data extracted by PSI processing chains to the detailed analysis of slow-moving DSGSDs in high mountain environments remains challenging. This is mainly due to low signal to noise ratio, severe atmospheric disturbances and snow cover effects that limit the occurrence of scatterers with sufficiently high coherence. This results in a sparse point cloud that is not dense enough to fully analyze the spatial displacement pattern.
Such specific problems can be overcome only by targeting DInSAR on the specific features and typical spatial-temporal scales of very slow rock slope deformations, integrating: (a) DInSAR processing over multiple temporal baselines; (b) detailed field morpho-structural constraints; (c) ground-based measurements.
In this perspective, we investigate the Mt. Mater DSGSD in Valle Spluga (Italian Central Alps), which affects a 1300 m-high glacial valley flank over an area of 3 km2. The slope impends over Madesimo village and ski resort and, according to published data [12,13], it is one of the fastest DSGSD in the Lombardia region. In order to unravel the kinematics, segmentation, and style of activity of Mt. Mater DSGSD, we integrate available SqueeSARTM products, site-specific DInSAR processing of Sentinel-1 radar images, GPS data, and detailed morpho-structural information from aerial photointerpretation and field surveys. Our results outline a strong spatial segmentation and suggest a complex seasonal pattern of slope deformations at Mt. Mater, providing important guidance to further modeling activities aimed at risk reduction.

2. Case Study: Mt. Mater DSGSD

2.1. Geology and Geomorphology

The slope is located on the eastern flank of the upper Valle Spluga (Lombardia, Central Italian Alps) and impends over Madesimo village and tourist resort (Figure 1).
The high relief slope (1550 to 3000 m a.s.l.) ranges in inclination between 33° (<2500 m a.s.l.) and 25° (>2500 m a.s.l.). The slope is made of metamorphic rocks of the Suretta and Tambò Penninic nappes [21], with dominant mica schist and paragneiss of the Stella-Timun complex (Suretta nappe) and limited outcrop of Mesozoic metasedimentary rocks at the slope toe, marking the contact with the underlying Tambò nappe. Rocks underwent a polyphase alpine tectono-metamorphic evolution with four main deformation stages since the Paleocene [22,23]. The first two stages involved the regional nappe structure and resulted in tight isoclinal folds associated with a pervasive foliation moderately dipping to the East. Locally, the third stage refolded previous structures with a subvertical N-verging axial plane [23,24]. Later ductile-brittle deformation resulted in the development of NNW-SSE trending normal faults cutting nappe boundaries [22] and in N-S and E-W trending master fractures. Tectonic strain and fabric are strongly heterogeneous due to the widespread occurrence of ductile shear zones, affecting the variability of rock strength and deformability in the brittle field.
The geomorphology of the upper Valle Spluga was strongly imprinted by Quaternary glaciations, especially by the Last Glacial Maximum (LGM [25]), that carved steep slopes in basement rocks with local reliefs up to 1500 m. After post-LGM deglaciation, landscape development was controlled by Lateglacial glacier re-advances [26], testified by well-preserved glacial cirques, steps and moraine ridges below 2350 m a.s.l. (Figure 1). During and following Lateglacial, extensive periglacial processes strongly affected rock degradation, slope stability, and sediment redistribution, as testified by widespread scree deposits and intact rock glaciers above 2400 m a.s.l. [27].

2.2. Gravitational Morpho-Structures

The post-glacial destabilization of the entire slope by means of deep-seated gravitational slope deformation (DSGSD) is testified by evident morpho-structural features [28,29], that affect the slope from the toe to the crest over an area exceeding 3 km2. These features have been recognized, mapped and interpreted through a geological and geomorphological investigation of the area, integrating 1:5000 field surveys, aerial photo-interpretation (stereo-photo coverage; TEM 1981-83; ortho-photos: 2000, 2003, 2007, 2012, 2015) and high-resolution digital elevation models (5 m resolution).
Above 2500 m a.s.l., these features include dominant scarps, suggesting limited or negligible rotational movements within the DSGSD mass. At 2900 m a.s.l., the slope is cut by a sharp triangular head scarp with a vertical downthrow of about 40 m (Figure 1a). Laterally, the scarp trace can be followed continuously to the N and S, suggesting that the DSGSD is bounded by a relatively well-developed basal shear zone. Between 2900 and 2500 m a.s.l., the slope is cut by several steep persistent N-S trending scarps, that extend to the south crossing the ridge with the Val Groppera. Here, NE-SE trenches and NW-SE counterscarps with a maximum length of few tens of meters define a small gravitational graben (Figure 1b).
Moving downslope, shallower arcuate scarps mark the transition to two nested large landslides, affecting the slope between 2400 m a.s.l. to the toe (Figure 1). In these sectors, rock masses are progressively more damaged moving downslope and become almost entirely crushed at 1900–2000 m a.s.l. Nested landslides are bounded by curve-shaped, highly fractured main scarps with vertical downthrow up to 20–25 m (Figure 1c,d). Several secondary scarps occur between 2200 and 2450 m a.s.l. within the landslide masses, suggest their internal segmentation in sectors with different degrees of damage (Figure 1a). Both the nested landslides have highly displaced toes, pushed out towards the valley floor through a series of slope-parallel lobes associated with debris sliding and a localized collapse of large landslide toes (Figure 1e). The main scarp of the northern nested landslide partially dismantled a moraine ridge at about 2300 m a.s.l., while DSGSD morpho-structures interact with intact rock-glaciers. This suggests that large-scale rock-slope instability was active in different stages after the Lateglacial.

3. Materials and Methods

We characterized the DSGSD using products of different spaceborne differential SAR interferometry (DInSAR) techniques. To analyze different aspects of the phenomenon, we used both point-like information from SqueeSARTM data (2015–2017) and spatially distributed information from specially processed multi-temporal baseline interferograms. We integrated this information with GPS ground measurements covering a period of 4 years from 2014 to 2019 (Figure 2).

3.1. SqueeSARTM Data

In order to test the potential of commercially available PSI data processed over large areas for the detailed characterization of a very slow rock slope deformation, we used SqueeSARTM data (TRE Altamira; Figure 3a,b). These were derived from Sentinel 1A/B SLC radar images (C-band, wavelength: 5.56 cm), acquired in TOPS Interferometric Wide swath (IW) mode [30] between March 2015 and July 2017 (revisit: 12 days, 6 days since March 2017; Supplementary Tables S1 and S2) in both ascending (track = 15; orbit azimuth δ = 10.23°; mean LOS angle θ = 41.99°) and descending geometries (track = 168; δ = 8.99°; θ = 41.78°). SqueeSARTM data report LOS displacement rates (VLOS). When the true displacement rate vector deviates from the LOS, InSAR sensitivity decreases making the interpretation of slow movements challenging. The Mt. Mater slope, facing to the W (mean aspect: 262°) with a mean slope gradient of 28°, is favorably oriented to the satellite LOS of descending orbits, able to catch slope movements with the highest sensitivity and without significant geometrical distortions.
To maximize the exploitation of SqueeSARTM data (Figure 3a,b), we integrated the information provided by ascending and descending acquisitions by decomposing the 2D velocity vector in the E-W trending vertical plane (2DInSAR; [14,31,32]). We developed a MatlabTM script to discretize the study area into a regular square grid (cell size: 25 m) and, for each cell, we averaged the LOS velocity values of PS/DS belonging to the same acquisition geometry. The average LOS velocity value computed in each cell for the two geometries was then assigned to pseudo-PS points, corresponding to the cell centroids (Figure 3c). For each pseudo-PS with both ascending and descending mean LOS velocity values, we extracted the E-W horizontal (Ve, Equation (1)) and vertical (Vz, Equation (2)) displacement rate components, as well as the magnitude (VT, Equation (3)) and inclination (dipT, Equation (4)) of the 2D displacement rate vector T.
V e = V d cos θ a V a cos θ d sin ( θ a + θ d )
V z = V d sin θ a V a sin θ d sin ( θ a + θ d )
V T = V z 2 + V e 2
d i p T = cos 1 V e V T
where: Va and Vd are the ascending and descending LOS velocities (mm/year), and θa and θd are the mean LOS angles for Sentinel 1 platform in the two acquisitions geometries.

3.2. DInSAR Processing

To further investigate the spatial and temporal patterns and heterogeneities of DSGSD displacements, we used the software SNAP (ESA Sentinel Application Platform 7.0.0 [33]) to generate 61 interferograms with increasing temporal baselines, ranging from 24 days to 1 year in the period between June 2016 and October 2019 (Figure 2 and Supplementary Table S3). Shorter baseline interferograms (6, 12 days) have been discarded because they are only able to capture fast displacements of scree deposits or periglacial features. To reduce speckle noise and better outline phase signatures in the interferograms, we applied multi-looking factors of 1 (azimuth) and 4 (range) and phase filtering techniques (Goldstein phase filtering [34]). For those interferograms with good signal to noise ratio and high coherence, we obtained phase-unwrapped displacement maps using the minimum cost flow algorithm (MCF, [35]) implemented in the SNAPHU software plugin [36].
Short temporal baseline (24-day) interferograms are frequently affected by atmospheric disturbances in the early June and September-October periods, due to the association of sharp topography and daily temperature variations. Thus, they were carefully selected to avoid submitting misleading information to the unwrapping procedure.
Seasonal interferograms with temporal baselines of several months (Figure 2) were processed considering snow cover occurrence from nivometric data (Campodolcino Meteo Station) and Landsat 8 (OLI) C1 images. Suitable pairs of snow-free images were selected to generate interferograms representative of displacements accumulated from June to October (hereafter summer interferograms) and from October to June (hereafter snow cover period interferograms).
Finally, annual interferograms were targeted to outline persistent (pluri-annual) displacement signals. We considered pairs of images with temporal baseline spanning 1-year (±6 days), rolling over the period June 2016–October 2019.

3.3. GPS Data

We considered 8 GPS benchmarks (operated by ARPA Lombardia), installed in the upper sector of the slope above 2500 m.a.s.l. (Figure 1). Periodical differential GNSS measurements (about 3 surveys per year), carried out using double-frequency GNSS Leica GS09 receivers, started in October 2014 for station 3 and in June 2015 for the other ones. Considering instrumentation specifications and the distance to the master control station (~1.5 Km), measurements are affected by a static (post-processed) horizontal standard error of 3–5 mm + 0.5 ppm and a vertical one of 6–10 mm + 1 ppm. No further corrections have been applied for environmental sources of error affecting measurements in field conditions.

4. Data and Results

We characterized the global kinematics (i.e., sliding mechanism) of the Mt. Mater DSGSD by exploiting the available SqueeSARTM dataset based on Sentinel-1 radar images (Figure 3a,b), from which we derived the products of 2D displacement rate vector decomposition (2DInSAR). The spatial distribution of the dipT variable along the slope allows identifying sectors undergoing lowering, slope-parallel translation or bulging/uplift, thus providing a consistent figure of global slope kinematics [31].
For Mt. Mater, our results indicate an average dipT value of about 30° over the slope, while higher values (around 80°) can be detected only close to the steep DSGSD upper scarps and to the headscarps of nested large landslides. These values and their spatial distribution provide robust insights into the global kinematics of the DSGSD. This is characterized by dominant translational sliding, with rotational sliding components associated with the head sectors of the nested large landslides mapped in the field (Figure 1 and Figure 3d).
VT values (Figure 3d) are generally higher in the upper DSGSD sectors, reaching 80 mm/year above the curved scarp at around 2600 m a.s.l. (Figure 1a) and testifying the ongoing active deformation of the entire slope. At and above the head sectors of the two nested large landslides, VT values range between 30 and 45 mm/year (head of northern landslide; Figure 1a,b) and 25–30 mm/year (above the head of southern landslide; Figure 1a,d). Displacement rates fade to 1–5 mm/year moving towards the slope toe, where landslide materials become more crushed and crumpled and/or the true displacement vector may significantly deviate from the radar line-of-sight.
However, the sparse nature of these data hampers a sound interpretation of the spatial pattern of measured displacements with respect to the scale (i.e., shallow vs deep-seated) and heterogeneity of slope failure mechanisms characteristic of different slope sectors. Thus, to better investigate displacement patterns along the slope, we analyzed the multi-temporal interferograms generated through targeted DInSAR processing across the period June 2016-October 2019 (Figure 4).
In all the processed interferograms, a low coherence area (<0.4) due to patchy vegetation and thermal stratification of the valley air, hampers the possibility to detect deformations close to the slope toe. This effect becomes more severe when increasing the temporal baseline from 24 days to 1-year, in which, the low coherence sector extends up to the middle slope sector because of increasing temporal decorrelation.
Interferograms with a baseline shorter than 24 days (i.e., 6, 12 days) do not show significant displacements, except for localized superficial debris covers or periglacial forms. From 24 days on, a triangular-shaped moving sector emerges from an almost homogeneous ground displacement signal (Figure 4a). It extends between the N-S curvilinear scarp trending at 2500 m a.s.l. and the DSGSD headscarp at 2900 m a.s.l. (Figure 1) and is made of damaged rock with extensive debris cover and periglacial features. This is the fastest slope sector, as it always decorrelates in longer baseline interferograms and is characterized by average displacement rates of 10 mm/month (Figure 4b).
These values are consistent with data provided by GPS benchmarks, all located within the upper slope sector between 2500 and 2650 m a.s.l. (Figure 1 and Figure 5a). As few GPS measurements were carried out during each year, for this comparison we referred to displacement rates averaged over the entire time span covered by GPS (2015 to 2019), that were projected along the descending Sentinel-1 satellite LOS.
GPS data show evident seasonal trends, characterized by increased displacement rates in the early summer-autumn period (June to October) as an effect of combined snowmelt and rainfall (Figure 5b). In fact, due to the high elevation (2800 m a.s.l.), snowmelt in this sector usually starts around the end of May, later than usually observed at the valley floor. In the snowmelt period, major portions of the slopes are still covered by snow, hampering DInSAR processing [12,37,38] and accurate quantification of snowmelt contributions to displacements.
Seasonal interferograms (Figure 4c) with longer temporal baselines clearly highlight continuous displacement fields, corresponding with: (a) localized deformation along the major DSGSD scarps at 2600 m and 2700 m a.s.l. (Figure 1 and Figure 4d), occurring at rates of several millimeters per seasonal cycle. This signal is very different from that of shallow debris movement detected in 24-day interferograms (Figure 4b) and reflects deep-seated slow deformations; (b) displacements of the head and internal sectors of the northern nested landslide below 2400 m a.s.l. (Figure 1 and Figure 4d), at rates up to 5 mm/month depending on the period (Figure 6); (c) movements up to 3–4 mm/month (Figure 6) of the southern nested landside head at 2100 m a.s.l. and the strongly segmented above sector, up to a neat curved scarp at 2500 m a.s.l. (Figure 1 and Figure 4d). DInSAR-derived displacement maps report for these sectors velocities ranging from a mean value of 1.5 to 4 mm/month, depending on the considered seasonal period (Figure 6). This suggests a differential seasonal response also for the middle slope sector, with variable trends conditioned by the snowmelt timing and rainfall input in the across winter and springtime.
To better constrain this observation, we selected three slope sectors (Figure 4f) corresponding to the nested landslides heads and above scarps (sectors 1 and 2) and the area in between (sector 3). We compared displacement rate distributions between the three sectors during consecutive seasonal periods, from 2016 to 2019. Boxplots in Figure 6 show that displacement rates change seasonally in response to the magnitude and timing of different hydrological inputs, with variable trends through the considered years. In 2016–2017, when snowmelt started early in April (Figure 5b), displacement rates recorded across the snow cover periods (Oct-Jun; about 4 mm/month) were higher than in the summer period (Jun-Oct; about 1.5 mm/month).
An opposite behavior is observed for 2017–2018. In this case, snowmelt occurred later (Figure 5b) and its contribution is recorded by the summer interferogram, with larger displacement rates between June and October (Figure 6). In the 2018–2019 period, despite late snowmelt (May; Figure 5a), abundant rainfall during autumn and spring triggered an anticipated acceleration, with mean LOS displacement rates of about 3.5 mm/month. While consistently reflecting seasonal changes, displacement rates in the three considered sectors confirm a sharp internal segmentation of the DSGSD. The northern nested landslide (sector 1) is always faster than the southern one (sector 2) and both are faster than sector 3, which is interpreted as representative of the background DSGSD activity.
Because of the low signal to noise ratio and poor coherence, we could not derive reliable unwrapped displacement maps from 1-year interferograms. However, wrapped phase maps (Figure 4e) still provide valuable information on the persistent, pluri-annual displacement signal of the DSGSD. In fact, two fringe cycles are evident in the nested landslide sectors and a coherent signal turns out between them in the central part, corresponding to the previously identified sector 3. Such signal unlikely results from shallow movements in the debris cover, that would decorrelate with such a long temporal baseline.

5. Discussion

A detailed slope-scale characterization of very slow rock slope deformations like deep-seated gravitational slope deformations [4,29] remains challenging. Successful slope-scale DInSAR applications have been presented in the literature for slow-moving landslides characterized by displacement rates far exceeding 5 cm/year and relatively homogeneous displacement fields [39,40,41,42]. Slower movements typical of DSGSD [13,29] are close to the limits of detection of DInSAR. In fact, DSGSD is characterized by complex mechanisms, reflected by heterogeneous displacement fields and creep at very slow rates. Moreover, these phenomena occur in high mountain environments, where atmospheric and snow-cover disturbances strongly affect the quality of remote sensing products. Although persistent-scatterers interferometry (PSI) techniques proved their ability to measure displacement rates of a few mm/year for coherent targets [12,13] due to their point-like nature, they often fail at capturing a complete picture of the spatial variability associated to slope failure mechanisms.
Our study shows that integrating field morpho-structural observations, remote sensing data from different InSAR techniques and ground-based monitoring allows a successful detail-scale investigation of heterogeneous, very slow-moving rock slope deformations. In particular, we characterized the activity and kinematics of the Mt. Mater DSGSD (Figure 7), wisely using remote sensing techniques to cope with practical difficulties due to: (a) displacement rates close to the limits of remote sensing techniques; (b) strong segmentation of displacement fields, due to secondary large landslides nested at different depths and characterized by different degree of both localized and distributed material damage; (c) lack of subsurface borehole or geophysical investigation data providing constraints to the interpretation of the deep structure of DSGSD.
In the first stage, we tested the potential of available SqueeSARTM PSI data, that are now used as a standard in the scientific community, for the characterization of very slow rock slope deformations. Performing 2DInSAR analysis to combine information from Sentinel 1A/B ascending and descending orbits, we fully exploited PSI data to retrieve a sound interpretation of the global kinematics of the DSGSD (Figure 3). However, these data have a point-like, sparse nature, making them unsuitable to capture the spatial heterogeneity of the phenomenon. Moreover, SqueeSARTM data processed over large areas are usually unable to account for the temporal patterns of activity of individual very slow-moving DSGSD and the site-specific atmospheric and snow cover disturbances, which greatly affect the local signal to noise ratio. Thus, we adopted a targeted DInSAR processing approach integrating multi-temporal baseline interferograms (2016–2019) constrained by detailed field morpho-structural observations and GPS data to unravel the deformation patterns and trends of activity in a spatially distributed fashion (Figure 4).
DInSAR processing, specifically designed on the spatial and temporal scales of slope processes associated with DSGSD, provides data consistent with ground-based GPS measurements. As the ability of GPS to measure very small displacements is limited by several error sources such as satellite coverage, temperature changes, and ground dilation, for each GPS benchmark we considered displacement rates averaged over the entire 4-year measurement period. Despite smoothing non-linearities of the GPS time series, this reduces the errors associated with individual measurements and provides a consistent figure of persistent displacement rates. GPS and DInSAR displacement rates consistently settle within about ±2 mm/24 days, according to the specific 24-day period covered by interferograms. Despite keeping the same spatial trend of DInSAR products obtained with targeted temporal baselines, SqueeSARTM data processed over large areas [37] always underestimate displacement rates (Figure 5a). This might be caused by reduced spatial coherence (quality) of the interferograms included in the SqueeSARTM analysis (mainly associated with long temporal baselines considered and/or the occurrence of relatively large spatial/ temporal phase gradients), which might result in underestimations of local surface velocities [42].
Selecting image pairs covering different temporal windows (i.e., 24-days, seasonal, annual) allows us to maximize the sensor detection capability for complex phenomena, extracting slope displacement rates, that in the studied case range between more than −15 mm/month and less than −3.5 mm/month. These very different values are typical of associated slope instability processes acting on different scales (e.g., shallow vs deep-seated sliding, nested landslides). Failing to discriminate among these different processes can lead to a misleading interpretation of the overall slope instability mechanisms with major impacts on the assessment of risk components (e.g., scenario volumes, intensity), that can be strongly over- or underestimated. Moreover, 24-day and seasonal interferograms effectively highlighted active sectors characterized by different spatial patterns and rates of displacements.
Combining mapping and morpho-structural information with remote sensing constraints on kinematic and spatial segmentation of the slope, we were able to provide an interpretative geological model of Mt. Mater slope (cross-section in Figure 7), even without the support of investigations. Different morpho-structures and their associations are witnesses of different deformation mechanisms and provide a first insight into the kinematic behavior of the landslide. We mapped several orders of persistent scarps that dissect the slope from the edge to the toe and only a few minor counterscarps that cut the southern border of the DSGSD towards the Groppera valley. This morpho-structural assemblage suggests a translational sliding mechanism characterized by synthetic structures developing at different depth levels and bounding discrete sectors. Field surveys and aerial photo-interpretation revealed that the slope deformation is sliced by secondary large landslides, nested at different depths within the main DSGSD. These are characterized by rotational movements in the head sectors, where the highest vertical displacements are recorded, and by a variable amount of internal damage. The presence of nested sectors reflects a deep complexity of slope deformation mechanisms, controlled by multiple shear zones almost parallel to the slope profile (Figure 7).
While 1-year interferograms provide a picture of long-term background DSGSD displacement signals, the combined analysis of seasonal interferograms and GPS data outline a sensitivity of the different slope sectors to hydrological forcing. Displacement rates in the middle-upper slope are dominantly sensitive to snowmelt, modulated by spring rainfall depending on the relative timing and magnitude of the two contributions. This induces differential responses across the snow cover period and the summer season. The large landslides nested in the main DSGSD mass exhibit similar styles of activity, but deform at displacement rates higher than the background signal of the DSGSD.
This may be due to the development of different shear zones at shallower depth from the toe to the top and their different degree of evolution that condition the response of the relative slope sectors to hydrological inputs.
Our approach, that can be applied also to other slow rock deformations in different geological and geomorphological settings, proved to supply key information (i.e., internal segmentation, style of activity, forcing) required to define reference scenarios for risk analysis and mitigation of a widespread, yet challenging class of slope instabilities.

6. Conclusions

Very slow rock slope deformations (e.g., DSGSD) in high mountain environments are challenging phenomena for InSAR applications, due to slow movements (usually less than 5 cm/year), spatial and temporal heterogeneity and strong impacts of atmospheric and snow-cover disturbances. We show how these issues can be handled through an integrated, targeted study approach combining field morpho-structural mapping, DInSAR analyses, and GPS data.
In this context, PSI data are very useful for a first-order characterization of slope activity and to retrieve its kinematics, but often fail to capture spatial segmentation, temporal trends and associated mechanisms for site-specific applications because of their point-like nature. DInSAR processing, targeted on multiple temporal baselines, allows us to successfully unravel the mechanisms, temporal trends of activity and forcing factors of heterogeneous DSGSDs characterized by discrete slope sectors possibly undergoing different evolution. Nevertheless, a sound interpretation of these remote sensing results requires detailed geological and morpho-structural field constraints.
Our observations depict a complex phenomenon characterized by nested sliding surfaces bounding discrete sectors that show different sensitivity to hydrological triggers and possibly undergoing different trends towards failure.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-4292/12/8/1329/s1, Table S1: List of images—SqueeSAR ascending dataset, Table S2: List of images—SqueeSAR descending dataset, Table S3: Interferograms derived by Sentinel-1 SLC image processing.

Author Contributions

Conceptualization, F.A. and C.C.; Data curation, C.C., F.F. and F.A.; Investigation, M.Z. and F.A.; Validation, G.B.C. and L.D.C.; Writing—original draft, C.C. and F.A.; Writing—review & editing, A.M., F.A., C.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Fondazione Cariplo, Grant 2016-0757 Slow2Fast to F. Agliardi.

Acknowledgments

We thank ARPA Lombardia for providing GPS data and Margherita Spreafico and Elena Valbuzzi for the fruitful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bovis, M.J.; Evans, S.G. Extensive deformations of rock slopes in southern Coast Mountains, southwest British Columbia, Canada. Eng. Geol. 1996, 44, 163–182. [Google Scholar] [CrossRef]
  2. Audemard, F.A.; Beck, C.; Carrillo, E. Deep-seated gravitational slope deformations along the active Boconó Fault in the central portion of the Mérida Andes, western Venezuela. Geomorphology 2010, 124, 164–177. [Google Scholar] [CrossRef]
  3. Agliardi, F.; Crosta, G.B.; Frattini, P.; Malusà, M.G. Giant non-catastrophic landslides and the long-term exhumation of the European Alps. Earth Planet. Sci. Lett. 2013, 365, 263–274. [Google Scholar] [CrossRef]
  4. Crosta, G.B.; Frattini, P.; Agliardi, F. Deep seated gravitational slope deformations in the European Alps. Tectonophysics 2013, 605, 13–33. [Google Scholar] [CrossRef]
  5. Lin, Z.; Kaneda, H.; Mukoyama, S.; Asada, N.; Chiba, T. Detection of subtle tectonic-geomorphic features in densely forested mountains by very high-resolution airborne LiDAR survey. Geomorphology 2013, 182, 104–115. [Google Scholar] [CrossRef]
  6. Crosta, G.B.; Agliardi, F.; Rivolta, C.; Alberti, S.; Dei Cas, L. Long-term evolution and early warning strategies for complex rockslides by real-time monitoring. Landslides 2017, 14, 1615–1632. [Google Scholar] [CrossRef]
  7. Preisig, G.; Eberhardt, E.; Smithyman, M.; Preh, A.; Bonzanigo, L. Hydromechanical rock mass fatigue in deep-seated landslides accompanying seasonal variations in pore pressures. Rock Mech. Rock Eng. 2016, 49, 2333–2351. [Google Scholar] [CrossRef]
  8. Riva, F.; Agliardi, F.; Amitrano, D.; Crosta, G.B. Damage-Based Time-Dependent Modeling of Paraglacial to Postglacial Progressive Failure of Large Rock Slopes. J. Geophys. Res. Earth Surf. 2018, 123, 124–141. [Google Scholar] [CrossRef] [Green Version]
  9. Pánek, T.; Klimeš, J. Temporal behavior of deep-seated gravitational slope deformations: A review. Earth Sci. Rev. 2016, 156, 14–38. [Google Scholar] [CrossRef]
  10. Agliardi, F.; Riva, F.; Barbarano, M.; Zanchetta, S.; Scotti, R.; Zanchi, A. Effects of tectonic structures and long-term seismicity on paraglacial giant slope deformations: Piz Dora (Switzerland). Eng. Geol. 2019, 263, 105353. [Google Scholar] [CrossRef]
  11. Rott, H.; Scheuchl, B.; Siegel, A.; Grasemann, B. Monitoring very slow slope movements by means of SAR interferometry: A case study from a mass waste above a reservoir in the Otztal Alps, Austria. Geophys. Res. Lett. 1999, 26, 1629–1632. [Google Scholar] [CrossRef]
  12. Wasowski, J.; Bovenga, F. Investigating landslides and unstable slopes with satellite Multi Temporal Interferometry: Current issues and future perspectives. Eng. Geol. 2014, 174, 103–138. [Google Scholar] [CrossRef]
  13. Frattini, P.; Crosta, G.B.; Rossini, M.; Allievi, J. Activity and kinematic behaviour of deep-seated landslides from PS-InSAR displacement rate measurements. Landslides 2018, 15, 1053–1070. [Google Scholar] [CrossRef]
  14. Eriksen, H.Ø.; Lauknes, T.R.; Larsen, Y.; Corner, G.D.; Bergh, S.G.; Dehls, J.; Kierulf, H.P. Visualizing and interpreting surface displacement patterns on unstable slopes using multi-geometry satellite SAR interferometry (2D InSAR). Remote Sens. Environ. 2017, 191, 297–312. [Google Scholar] [CrossRef] [Green Version]
  15. Frattini, P.; Crosta, G.B.; Allievi, J. Damage to buildings in large slope rock instabilities monitored with the PSinSARTM technique. Remote Sens. 2013, 5, 4753–4773. [Google Scholar] [CrossRef] [Green Version]
  16. Macciotta, R.; Martin, C.D.; Edwards, T.; Cruden, D.M.; Keegan, T. Quantifying weather conditions for rock fall hazard management. Georisk 2015, 9, 171–186. [Google Scholar] [CrossRef]
  17. Béjar-Pizarro, M.; Notti, D.; Mateos, R.M.; Ezquerro, P.; Centolanza, G.; Herrera, G.; Bru, G.; Sanabria, M.; Solari, L.; Duro, J.; et al. Mapping vulnerable urban areas affected by slow-moving landslides using Sentinel-1InSAR data. Remote Sens. 2017, 9, 876. [Google Scholar] [CrossRef] [Green Version]
  18. Gabriel, A.K.; Goldstein, R.M.; Zebker, H.A. Mapping small elevation changes over large areas: Differential radar interferometry. J. Geophys. Res. 1989, 94, 9183–9191. [Google Scholar] [CrossRef]
  19. Colesanti, C.; Ferretti, A.; Prati, C.; Rocca, F. Monitoring landslides and tectonic motions with the Permanent Scatterers Technique. Eng. Geol. 2003, 68, 3–14. [Google Scholar] [CrossRef]
  20. Ferretti, A.; Prati, C.; Rocca, F. Permanent scatterers in SAR interferometry. IEEE Trans. Geosci. Remote Sens. 2001, 39, 8–20. [Google Scholar] [CrossRef]
  21. Schmid, S.M.; Fügenschuh, B.; Kissling, E.; Schuster, R. Tectonic map and overall architecture of the Alpine orogen. Eclogae Geol. Helv. 2004, 97, 93–117. [Google Scholar] [CrossRef]
  22. Baudin, T.; Marquer, D.; Persoz, F. Basement-cover relationships in the Tambo nappe (Central Alps, Switzerland): Geometry, structure and kinematics. J. Struct. Geol. 1993, 15, 543–553. [Google Scholar] [CrossRef]
  23. Marquer, D.; Challandes, N.; Baudin, T. Shear zone patterns and strain distribution at the scale of a Penninic nappe: The Suretta nappe (Eastern Swiss Alps). J. Struct. Geol. 1996, 18, 735–764. [Google Scholar] [CrossRef] [Green Version]
  24. Scheiber, T.; Pfiffner, O.A.; Schreurs, G. Strain accumulation during basal accretion in continental collision—A case study from the Suretta nappe (eastern Swiss Alps). Tectonophysics 2012, 579, 56–73. [Google Scholar] [CrossRef]
  25. Bini, A.; Buoncristiani, J.-F.; Couterrand, S.; Ellwanger, D.; Felber, M.; Florineth, D.; Graf, H.R.; Keller, O.; Kelly, M.; Schlüchter, C.; et al. Die Schweiz während des letzteiszeitlichen Maximums (LGM); Bundesamt für Landestopographie, Swisstopo: Knitz, Switzerland, 2009.
  26. Ivy-Ochs, S.; Kerschner, H.; Reuther, A.; Preusser, F.; Heine, K.; Maisch, M.; Kubik, P.W.; Schlüchter, C. Chronology of the last glacial cycle in the European Alps. J. Quat. Sci. 2008, 23, 559–573. [Google Scholar] [CrossRef]
  27. Scotti, R.; Brardinoni, F.; Alberti, S.; Frattini, P.; Crosta, G.B. A regional inventory of rock glaciers and protalus ramparts in the central Italian Alps. Geomorphology 2013, 186, 136–149. [Google Scholar] [CrossRef]
  28. Agliardi, F.; Crosta, G.; Zanchi, A. Structural constraints on deep-seated slope deformation kinematics. Eng. Geol. 2001, 59, 83–102. [Google Scholar] [CrossRef]
  29. Agliardi, F.; Crosta, G.B.; Frattini, P. Slow rock-slope deformation. Landslides Types Mech. Model. 2012, 207–221. [Google Scholar]
  30. De Zan, F.; Guarnieri, A.M. TOPSAR: Terrain observation by progressive scans. IEEE Trans. Geosci. Remote Sens. 2006, 44, 2352–2360. [Google Scholar] [CrossRef]
  31. Manzo, M.; Ricciardi, G.P.; Casu, F.; Ventura, G.; Zeni, G.; Borgström, S.; Berardino, P.; Del Gaudio, C.; Lanari, R. Surface deformation analysis in the Ischia Island (Italy) based on spaceborne radar interferometry. J. Volcanol. Geotherm. Res. 2006, 151, 399–416. [Google Scholar] [CrossRef]
  32. Dalla Via, G.; Crosetto, M.; Crippa, B. Resolving vertical and east-west horizontal motion from differential interferometric synthetic aperture radar: The L’Aquila earthquake. J. Geophys. Res. Solid Earth 2012, 117. [Google Scholar] [CrossRef] [Green Version]
  33. ESA. SNAP|STEP. 2018. Available online: http://step.esa.int/main/toolboxes/snap/ (accessed on 10 November 2019).
  34. Goldstein, R.M.; Werner, C.L. Radar interferogram filtering for geophysical applications. Geophys. Res. Lett. 1998, 36, 813–821. [Google Scholar] [CrossRef] [Green Version]
  35. Costantini, M. A novel phase unwrapping method based on network programming. IEEE Trans. Geosci. Remote Sens. 1998, 36, 813–821. [Google Scholar] [CrossRef]
  36. Chen, C.W.; Zebker, H.A. Network approaches to two-dimensional phase unwrapping: Intractability and two new algorithms. J. Opt. Soc. Am. A 2000, 17, 401–414. [Google Scholar] [CrossRef]
  37. Ferretti, A.; Fumagalli, A.; Novali, F.; Prati, C.; Rocca, F.; Rucci, A. A new algorithm for processing interferometric data-stacks: SqueeSAR. Proc. IEEE Trans. Geosci. Remote Sens. 2011, 49, 3460–3470. [Google Scholar] [CrossRef]
  38. Bovenga, F.; Belmonte, A.; Refice, A.; Pasquariello, G.; Nutricato, R.; Nitti, D.O.; Chiaradia, M.T. Performance analysis of satellite missions for multi-temporal SAR interferometry. Sensors 2018, 18, 1359. [Google Scholar] [CrossRef] [Green Version]
  39. Handwerger, A.L.; Roering, J.J.; Schmidt, D.A. Controls on the seasonal deformation of slow-moving landslides. Earth Planet. Sci. Lett. 2013, 377, 239–247. [Google Scholar] [CrossRef]
  40. Schlögel, R.; Doubre, C.; Malet, J.P.; Masson, F. Landslide deformation monitoring with ALOS/PALSAR imagery: A D-InSAR geomorphological interpretation method. Geomorphology 2015, 231, 314–330. [Google Scholar] [CrossRef]
  41. Kos, A.; Amann, F.; Strozzi, T.; Delaloye, R.; von Ruette, J.; Springman, S. Contemporary glacier retreat triggers a rapid landslide response, Great Aletsch Glacier, Switzerland. Geophys. Res. Lett. 2016, 43, 12–466. [Google Scholar] [CrossRef]
  42. Manconi, A.; Kourkouli, P.; Caduff, R.; Strozzi, T.; Loew, S. Monitoring surface deformation over a failing rock slope with the ESA sentinels: Insights from Moosfluh instability, Swiss Alps. Remote Sens. 2018, 10, 672. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Main features of the Mt. Mater deep-seated gravitational slope deformation (DSGSD). (a) Simplified map portraying the main morpho-structural features associated with DSGSD and nested large landslides. GPS benchmark locations are outlined. (b) tranches on the S side of the main DSGSD body; (c) main scarp of the northern nested landslide (cross-section in Figure 6); (d) damaged rock mass at the top of the southern nested landslide; (e) blocky accumulation at the southern nested landslide toe.
Figure 1. Main features of the Mt. Mater deep-seated gravitational slope deformation (DSGSD). (a) Simplified map portraying the main morpho-structural features associated with DSGSD and nested large landslides. GPS benchmark locations are outlined. (b) tranches on the S side of the main DSGSD body; (c) main scarp of the northern nested landslide (cross-section in Figure 6); (d) damaged rock mass at the top of the southern nested landslide; (e) blocky accumulation at the southern nested landslide toe.
Remotesensing 12 01329 g001
Figure 2. Time windows and temporal baselines covered by remote sensing and monitoring data.
Figure 2. Time windows and temporal baselines covered by remote sensing and monitoring data.
Remotesensing 12 01329 g002
Figure 3. Kinematics of the Mt. Mater DSGSD derived from SqueeSARTM data analysis. (a,b) permanent (PS) and distributed scatterers (DS), classified by line-of-sight (LOS) displacement rate (VLOS); (c) pseudo-PS derived by combining data from ascending and descending orbit; (d) 3D visualization of the 2D total displacement vector T, derived by the 2D decomposition of LOS velocities at pseudo-PS locations. Vector color and length are classified by T vector velocity (VT).
Figure 3. Kinematics of the Mt. Mater DSGSD derived from SqueeSARTM data analysis. (a,b) permanent (PS) and distributed scatterers (DS), classified by line-of-sight (LOS) displacement rate (VLOS); (c) pseudo-PS derived by combining data from ascending and descending orbit; (d) 3D visualization of the 2D total displacement vector T, derived by the 2D decomposition of LOS velocities at pseudo-PS locations. Vector color and length are classified by T vector velocity (VT).
Remotesensing 12 01329 g003
Figure 4. Examples of interferograms and displacement maps obtained by DInSAR with different temporal baselines (images: Sentinel-1 (Track 66, descending). Left column: wrapped phase interferograms with increasing temporal baseline: (a) 24 days (2–26 July 2017); (c) seasonal (June–October 2019); (e) 1-year (2018–2019). Right column: (b) and (d) unwrapped phase converted to displacements (mm). The unwrapped results for the 1-year interferogram were affected by a high noise level and unreliable values, and for reason are omitted; (f) sectors 1, 2 and 3, representative of the activity of northern and southern nested landslide heads and the background DSGSD, respectively. 24-days and seasonal wrapped and unwrapped interferograms are masked where the phase coherence of the filtered interferogram is < 0.4 and displacement >0. On 1-year interferograms pixels with coherence <0.3 are masked out. The reference point used for phase to displacement conversion is shown.
Figure 4. Examples of interferograms and displacement maps obtained by DInSAR with different temporal baselines (images: Sentinel-1 (Track 66, descending). Left column: wrapped phase interferograms with increasing temporal baseline: (a) 24 days (2–26 July 2017); (c) seasonal (June–October 2019); (e) 1-year (2018–2019). Right column: (b) and (d) unwrapped phase converted to displacements (mm). The unwrapped results for the 1-year interferogram were affected by a high noise level and unreliable values, and for reason are omitted; (f) sectors 1, 2 and 3, representative of the activity of northern and southern nested landslide heads and the background DSGSD, respectively. 24-days and seasonal wrapped and unwrapped interferograms are masked where the phase coherence of the filtered interferogram is < 0.4 and displacement >0. On 1-year interferograms pixels with coherence <0.3 are masked out. The reference point used for phase to displacement conversion is shown.
Remotesensing 12 01329 g004
Figure 5. Spatial and temporal trends of activity from GPS and DInSAR data. (a) spatial trends of displacement rates derived from SqueeSARTM, DInSAR and GPS data, extracted at GPS benchmark locations and scaled to a 24-day period. SqueeSARTM data were processed over a 2-year acquisition period; GPS displacement rates are derived from cumulative displacements averaged over the 4-year observation period and projected along the Sentinel-1 LOS; DInSAR displacement rates were extracted from three unwrapped 24-day interferograms (P1: 14/06/2017–14/07/2017; P2: 02/07/2017–26/07/2017; P3: 12/09/2017–06/10/2017); (b) time series of cumulative GPS 3D displacements compared to daily rainfall and snow cover (data: ARPA Lombardia).
Figure 5. Spatial and temporal trends of activity from GPS and DInSAR data. (a) spatial trends of displacement rates derived from SqueeSARTM, DInSAR and GPS data, extracted at GPS benchmark locations and scaled to a 24-day period. SqueeSARTM data were processed over a 2-year acquisition period; GPS displacement rates are derived from cumulative displacements averaged over the 4-year observation period and projected along the Sentinel-1 LOS; DInSAR displacement rates were extracted from three unwrapped 24-day interferograms (P1: 14/06/2017–14/07/2017; P2: 02/07/2017–26/07/2017; P3: 12/09/2017–06/10/2017); (b) time series of cumulative GPS 3D displacements compared to daily rainfall and snow cover (data: ARPA Lombardia).
Remotesensing 12 01329 g005
Figure 6. Seasonal response of slope sectors. Distributions of LOS displacement rates extracted from seasonal interferograms (June–October and October–June) for the 3 sectors highlighted in Figure 4f. Whisker lengths: upper and lower inner fences based on interquartile range (IQR).
Figure 6. Seasonal response of slope sectors. Distributions of LOS displacement rates extracted from seasonal interferograms (June–October and October–June) for the 3 sectors highlighted in Figure 4f. Whisker lengths: upper and lower inner fences based on interquartile range (IQR).
Remotesensing 12 01329 g006
Figure 7. Interpretative cross-section across the northern sector of the Mt. Mater slope. DSGSD morpho-structures, basal shear zones and nested landslides are reconstructed from field evidence and interpretation of DInSAR displacement patterns. Cross-section trace is shown in the inset.
Figure 7. Interpretative cross-section across the northern sector of the Mt. Mater slope. DSGSD morpho-structures, basal shear zones and nested landslides are reconstructed from field evidence and interpretation of DInSAR displacement patterns. Cross-section trace is shown in the inset.
Remotesensing 12 01329 g007

Share and Cite

MDPI and ACS Style

Crippa, C.; Franzosi, F.; Zonca, M.; Manconi, A.; Crosta, G.B.; Dei Cas, L.; Agliardi, F. Unraveling Spatial and Temporal Heterogeneities of Very Slow Rock-Slope Deformations with Targeted DInSAR Analyses. Remote Sens. 2020, 12, 1329. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12081329

AMA Style

Crippa C, Franzosi F, Zonca M, Manconi A, Crosta GB, Dei Cas L, Agliardi F. Unraveling Spatial and Temporal Heterogeneities of Very Slow Rock-Slope Deformations with Targeted DInSAR Analyses. Remote Sensing. 2020; 12(8):1329. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12081329

Chicago/Turabian Style

Crippa, Chiara, Federico Franzosi, Mattia Zonca, Andrea Manconi, Giovanni B. Crosta, Luca Dei Cas, and Federico Agliardi. 2020. "Unraveling Spatial and Temporal Heterogeneities of Very Slow Rock-Slope Deformations with Targeted DInSAR Analyses" Remote Sensing 12, no. 8: 1329. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12081329

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