1. Introduction
Differential SAR (Synthetic Aperture Radar) Interferometry techniques can offer great support to natural hazard detection and monitoring since they provide ground displacement data characterized by rapid and easily updatable measurements, millimetric and centimetric precision, high spatial resolution, good temporal sampling, and wide area coverage, consequently, reducing efforts and costs. In particular, multi-temporal interferometric techniques,
i.e., Persistent Scatterer Interferometry (PSI), have recently demonstrated to be a valuable and successful tool to detect slope instability and map geomorphological processes, such as slow-moving landslides [
1–
12].
At present, the availability of large stacks of SAR data from both historical and new-generation satellite sensors permits to spatially and temporally detect landslide displacements at regional and local scale with improved resolution and precision [
13–
20], especially on wide areas that are too large or inaccessible to be investigated through conventional monitoring techniques.
Nevertheless, PSI techniques hold some limitations, mainly due to the acquisition geometry of the satellite systems. PSI-based displacement data are mono-dimensional along the satellite LOS (Line Of Sight); consequently, they are dependent on the relationships between sensor acquisition geometry, the local topography, and the real direction of movement. This is critical when dealing with landslides, as only a fraction of slope movement component can be captured in the LOS. Additionally, the same ground movement can be estimated with opposite sign and different module from ascending and descending orbits, making the slope dynamics interpretation not immediately intelligible [
6,
8,
21–
23].
In this paper, a methodology to fully exploit PS data acquired in L-band by ALOS (Advanced Land Observing Satellite) satellite equipped with PALSAR (Phased Array type L-band Synthetic Aperture Radar) sensor is presented for detecting and representing landslide displacements in the Tramuntana range in the island of Majorca (Spain), taking into account both satellite parameters and geomorphological features. The main aim is assessing the ground velocity, the geometrical visibility, and radar detection suitability of the phenomena, and generating landslide activity maps, which can represent alert maps of the study area.
The Tramuntana range is the main mountainous alignment of Majorca, where a wide range of slope movements occurs, among which the most frequent ones are rockfalls, rotational, and complex landslides [
24]. The Tramuntana range is also characterized by a high exposure of elements at risk, due to the presence of several tourist resorts and apartments, and the Ma-10 road that represents the main road of the island and suffered several damages across time, due to both rockfalls and landslides.
After elaborating landslide and ground motion activity maps, an innovative confidence degree evaluation is finally carried out, in order to attest the reliability of these provided products, resting upon the availability of PS and further external data, i.e., damages inventory.
2. Study Area: The Tramuntana Range, Majorca (Spain)
The study area is the Tramuntana Range in Majorca (Spain), the main mountainous alignment of the island, extended up 1,480 km
2 in the northwestern part of Majorca [
25] (
Figure 1a). This Alpine range spreads parallel to the north coast, being approximately 90 km maximum long and meanly 15 km wide. The mountain line exceeds 600 m a.s.l. (above sea level) everywhere, and the central sector (Lluc area) is the highest, with several peaks higher than 1,000 m, including Puig Major (1,445 m a.s.l.), which is the highest peak on Majorca island [
26]. The north coastal side of Tramuntana mountain range is the most abrupt and rugged, with prevalent very high cliffs overlooking the sea and steep slopes with 30–40% gradient percentage, while the southern side is smoother, characterized by a gentler gradient (
Figure 1c). Flat regions cover 4% of the total area. Slopes S-SE oriented are slightly frequent than the other orientations, due to the structural outline of the range (
Figure 1d).
The present geomorphologic setting is determined by the geological structure of the chain and was formed during the Miocenic structuring linked to the Alpine fold, characterized by a series of NW overlapping thrusts [
27]. From the geological point of view, the main framework of the Tramuntana range is constituted by Mesozoic carbonate rocks (limestones and dolostones) (
Figure 1b). The lowest units of the stratigraphic sequence of the Tramuntana range are composed of Triassic siliceous sandstone (Buntsandstein, Lower Triassic), dolostones and marls (Rhaetian, Upper Triassic), and in some places clay and gypsum (Keuper, Upper Triassic). Jurassic limestones and Lias dolostones overlap the Triassic series. The upper units are conglomerate and marls (Paleogene-Miocene) and recent colluvial sediments composed of Quaternary coarse debris, which cover most of the slopes with several meters of thickness (
Figure 1b).
From a structural point of view, the Tramuntana range is an NE-SW directed imbricated system of overlapping thrusts [
28]. The regional detachment level of the nappe thrusts is located in the Keuper sediments (clays with gypsum and altered volcanic rocks) [
27,
29].
2.1. Landslide Movements on Tramuntana Range
The landslides of the Tramuntana range have been studied and described by Mateos [
30] and Mateos
et al. [
25,
28]. The variety of lithologies cropping out in this mountain chain determines a wide range of slope movements [
24]. Therefore, the geological structure of the Tramuntana range conditions the typology and the spatial distribution of the recent mass movements, as well as the presence of faults and discontinuities that often control the planes and the local failure surfaces of ground motions. Shallow landslides are frequent phenomena: most of the movements are translational slides and earth flows or complex movements, mostly involving the soft sediments from the Upper Triassic (clays with gypsum), and the Paleogene-Neogene loam deposits that occasionally outcrop along the mountain range. Movements, such as debris slides and debris flows, are also widespread as most of the slopes in the Tramuntana range are covered with these coarse colluviums. Nevertheless, the most frequent slope movements in the Tramuntana range are rockfalls due to the predominance of Jurassic rocky massifs made of limestone and dolostones [
30].
The steep topography of the chain, related to its geological and structural complexity, determines an intense slope dynamics [
25], quickened by the Mediterranean climate. The island of Majorca has a typically Mediterranean climate, with mild winters and warm-dry summers. The maximum precipitation takes place during the autumn months due to the arrival of the first cold air masses, which contrast with the high temperature retained by the sea. This phenomenon is known as “gota fría” (cold rain), in which heavy storms are accompanied by intense rainfall with episodes of up to 250 mm in 24 h [
25]. The orography of the island clearly controls the distribution of precipitation. The central sector of the Tramuntana range registers an average annual precipitation of 1,200 mm, which gradually drops towards the extreme SW of the range, where average annual precipitation is no more than 300 mm. Mateos
et al. [
25] reveals that recent rockfalls have also occurred after several freeze-thaw cycles, being a determining and unusual factor in this warm region.
The economy of the Tramuntana range is based exclusively around tourism, which accounts for 95% of its income. In particular, the dwellings and the road network are of great importance in this area as it welcomes 8.5 million visitors every year [
25].
The occurrence of many shallow landslides cause huge damages on the island, e.g., the Estellencs landslide (1971, 1979, 2010), the Fornalutx landslide (1924, 1974, 2002), and the Biniarroi landslide (1721, 1813, 1943) that are the most important mass movements that occurred in the Balearic islands, due to both their dimensions and caused damages [
26].
A total of 918 events have been recorded in the Tramuntana range, being 65% rockfalls and 35% landslides. From this amount, 19% of the events have produced damage to buildings, structures and infrastructure. Recently, during the years 2008 to 2010, Majorca experienced one of the coldest and wettest winters in living memory, leading to the occurrence of 34 mass movements that produced damages to dwellings and road networks. Extensive damages occurred, especially on the Ma-10 Road, which is the main road of the northern part of Tramuntana range with an intense vehicular traffic, causing, therefore, significant economic losses, estimated at about 11M Euro [
25].
The rockfall and landslide damage inventory dated up to 2011 is pointed out in
Figure 1e: it was elaborated at a 1:1,000 scale from different sources, and accounts 172 records (56 of them corresponding to landslide occurrences and 116 to rockfalls).
The available pre-existing inventory of ground movements in Tramuntana range includes rockfalls and landslides, covering the period 2008–2010 (
Figure 1e). It was produced at a 1:5,000 scale by means of orthophoto and field-work [
25,
30].
The PSI technique for mapping and monitoring mass movements can only be applied to slow-moving landslide phenomena that are classified, according to Cruden and Varnes [
31], as extremely slow (velocity < 16 mm/yr) and very slow landslides (16 mm/yr < velocity < 1.6 m/yr), due to the acquisition parameters (
i.e., signal wavelength, revisit cycle) of the satellite systems. Thus, the comparison of the PS spatial distribution with the mass movement inventory was primarily focused on the inventoried slow-moving phenomena,
i.e., landslides, and on the related damage inventory due to landslides. The pre-existing landslide inventory contains 206 landslides covering an area of 11 km
2.
4. Methodology
The landslide activity map integrates the PSI data, the post-processing derived products and the improved landslide inventory.
Input data (grey panels in
Figure 2) include morphological and satellite PS data, the pre-existing landslide inventory (LSI), as well as ground truth data,
i.e., damage database and field survey of the study area. Within the PSI post processing stage (
Figure 2a), the visibility of the targeted study area is assessed and the LOS PS displacements are projected downslope. The visibility assessment, calculated through the computation of two indexes (R- and LU-index), is useful to identify slopes, PS targets, or landslides with an adequate geometry for satellite measurements (
i.e., slope direction parallel to the satellite LOS direction), to be distinguished from those where geometrical distortions will be introduced in the measurement. These geometrical distortions, characteristic of mountainous environments, come out from the combination of the local topography with respect to the satellite acquisition parameters [
8,
34] and can be reduced through the projection of satellite Line Of Sight velocities (V
LOS) along the local steepest slope (V
SLOPE). In a second stage (
Figure 2b), PSI data and post-processing derived products are integrated with the landslide inventory to generate a “landslide activity map”. Therefore, for every landslide phenomenon, both the landslide visibility assessments as well as the displacement measurements are provided.
Separately, we also analyze PS benchmarks that were discarded through the downslope projection and that are not referred to any landslide of LSI, but even moving, since they could deal with unmapped landslides or other kinds of terrain motion processes. These clusters of moving PS are considered for the generation of a “ground motion activity map”, being the target of further necessary field investigations that could explain the mechanisms generating measured displacements. In the last step (
Figure 2c), the confidence degree of the activity maps is evaluated (high, medium, or low). This evaluation aims at assessing how reliable the PSI measured displacement is to represent a landslide movement. It is important to stress that the reliability on the PSI measurement itself is not evaluated, but whether this measurement is related to landslide dynamic or not.
5. PSI Post-Processing Phase
5.1. PS Visibility Analysis
The geometrical visibility of the area of interest is evaluated with respect to the morphology and land-use, and to the acquisition geometry of the available satellite system [
8,
22,
34,
35]. The R- and LU-indexes [
8] are calculated for the whole area and assigned to each PS. Then, these indexes are also assigned to every landslide by computing the average or the moda value of all the PS included within its boundary, in order to assess how good the ground surface covered by each phenomenon is to backscatter the satellite radar signal with respect to its relief and land use.
The R-Index (RI) is based on the approach proposed by Notti
et al. [
8] and it is calculated through
Equation (1) that represents the ratio between the slant range and the ground range, taking into account the acquisition geometry of the radar (
i.e., incidence angle = θ; angle from North of the satellite track = α) and the geometry of the ground surface (slope and aspect layers, S and A, respectively, derived from a Digital Eleveation Model –DEM- with 20 m resolution):
The RI ranges from the maximum value +1 to −1 and four classes are defined (
Figure 3). In particular, it assumes negative or around 0 values (class 1) when no PS can be detected or when geometrical distortions (
i.e., layover and shadow effects for negative values and foreshortening effect for RI ≈ 0) strongly affect the area. Values from 0.66 to +1 are estimated when the suitability of the area to be investigated through the PS dataset is the best: this means that the slope direction of the area is approximately parallel to the satellite LOS direction (class 4). RI values from 0.33 to 0.66 include either PS located on slopes with an acceptable geometry with respect to the satellite LOS or PS located in flat areas (column portion with darker colors in
Figure 3c). RI values between 0.00 and 0.33 indicate that the slope geometry is not favorable for the satellite acquisition geometry and a lower PS number is retrieved. As it expected, in the Tramuntana range, a greater number of PS is found for an increasing RI, as long as the PS located in flat areas (darker color columns in
Figure 3c) are removed.
In order to assess the suitability of every landslide to be monitored by ALOS PALSAR satellite, the average RI of all the PS included within the landslide geometry is calculated. From this analysis, 78 landslides exhibit the best geometry to be monitored (i.e., slope direction approximately parallel to the satellite LOS direction). Assuming that the landslide propagation follows the steepest slope, the measured LOS displacements will represent, almost completely, the real motion of these landslides. On the other hand, there are 133 landslides with RI < 0.66, where the geometry of the slope will prevent LOS displacement from fully representing the real landslide motion.
The Land Use Index (LU-index) calculation permits to analyze the effect of the land coverage on PS detectability [
8]. Five main categories of land use, based on the theoretical radar signal penetration properties, are defined: (i) areas with buildings and structures; (ii) rocks, talus of wide dimensions, rocky areas with sparse vegetation; (iii) cultivated fields; (iv) scrublands with bushes and woods; and (v) water reservoirs. Land-use categories on Tramuntana range, extracted from the Corine Land Cover of the Balearic Islands, are homogenized and reduced to the five classes (
Figure 3b). In
Figure 3d, the PS number in each LU class is also shown. It can be seen that rocky outcrops with sparse vegetation show the highest number of PS, followed by scrubland with bushes and woods. This is due to the dominant presence of massive rock outcrops in vegetated areas all through the Tramuntana range. The calculation of the LU moda value for all the PS detected within each landslide reveals that 180 landslides (85% of the total) are monitored in vegetated environments with a dominant presence of massif rocks (classes 2 and 4) (
Figure 3d).
5.2. LOS Displacement Projection along the Steepest Slope
PS displacement represents the one-dimensional projection in the satellite LOS direction of the real movement that actually occurs in all three dimensions. In order to overcome the differences given by ground geometry and satellite LOS, all the PS average yearly V
LOS (mm/yr) are projected into the same direction of the steepest slope, by means of
Equation (2), obtaining V
SLOPE (mm/yr) [
22,
23]. Thus, this downslope projection is designed to compare landslide velocities with different slope orientations, enriching V
LOS information.
The C coefficient is the fraction of the 3D displacement that can be measured by PS targets and β the angle between the steepest slope and the LOS direction (
Figure 4b). Several limitations of this method must be taken into account: when β is almost 90°, C is close to 0 and V
SLOPE tends to infinity. Following the work done by Herrera
et al. [
21], an absolute maximum value of β = 72° corresponding to cos β = 0.3 is fixed and, as a result, V
SLOPE cannot be higher than 3.33 times the V
LOS. This threshold corresponds to the condition number of 15 proposed by Cascini
et al. [
22] as the number for the inversion matrix solving the algebraic system used for the projection process. In order to reduce any V
SLOPE exaggeration, it is assumed C = −0.3 when −0.3 < C< 0 and C = 0.3 when 0 < C < 0.3. Whereas PS V
SLOPE values turn positive (V
SLOPE > 0) they are discarded [
21]. This is because a positive V
SLOPE would represent uphill movement. Although positive displacements may be present at the toe of landslides where vertical movements can occur, the horizontal displacement vector should remain oriented downhill [
21].
Additionally, the velocity for only PS localized over ≥5° slopes is projected, since in this study area landslides are unreported when the slope is smaller. Finally, the downslope projection assumes that the local steepest slope is the most probable direction of movement and so the landslide motion is parallel to the slope. Nevertheless the surface runoff of some mass movements could be not parallel to the slope, e.g., at the crown of a rotational landslide where also vertical motions can take place.
In the Tramuntana range the V
SLOPE calculation reduces the V
LOS PS population to 45% (
Table 2).
Even though the reduction of PS data is significant, this transformation improves their integration with the landslides inventory. Exploiting the V
LOS measurements, only 71 landslides (with almost 1 PS) show a negative displacement (
i.e., downhill displacement or away from the satellite), whereas 95 of them exhibit a positive V
LOS (
i.e., uplift or towards the satellite). As it was explained in the previous section, this is due to the distortion introduced by the satellite acquisition geometry with respect to the slope relief. Once the V
LOS is projected along the steepest slope, a total of 148 landslides reveal a negative displacement, among which 77 landslides (with almost 1 PS) show an average V
SLOPE < −5 mm/yr (
Table 3).
5.3. Velocity Thresholds Determination
Velocity thresholds need to be established for distinguishing moving from not moving PS. Some authors, such as Herrera
et al. [
21] and Righini
et al. [
9], consider 2 mm/yr as the V
LOS stability threshold for PSI C-band data. Other works [
11] use a threshold of 5 mm/yr for V
SLOPE obtained by a calibration procedure. However, there are no previously published thresholds for PSI L-band data results, due to the reduced number of applications.
Herein, authors use velocity threshold criteria based on the statistical analysis of PS velocity distributions. As V
LOS and V
SLOPE PS populations show different velocity distribution patterns, being the first approximately normal (Gaussian) and the second negatively skewed due to the PS reduction, the two PS distributions are classified separately (
Figure 4).
The classification of V
LOS into velocity classes is based on the Inter-Quantile Range (IQR) around zero (
Figure 4a). The central position (Q2) is zero, while Q1 and Q3 are the upper and lower positions, found by adding/detracting the standard deviation (σ) value to Q2 (0 ± σ). Thus, the IQR, which represents the “dispersion” of data around zero, is considered the stable threshold range. The further lower and upper velocity classes are derived using the σ value.
Thus, the velocity stable range for L-band radar displacement measurements (VLOS) is basically chosen taking into account the statistical velocity distribution and the standard deviation of the ALOS dataset in the study area.
Most of the published PSI analyses establish the stability V
LOS threshold at ±2 mm/yr [
9,
21]. Taking into account that V
LOS sensitivity increases with the radar frequency, a greater threshold would be expected for L-band (1.2 GHz). Evidently, ALOS L-band PS data have a higher standard deviation than C-band in similar conditions, but this outcome also depends on the number of images, on the temporal period and quality of selected points. In this case, L-band results are approaching C-band map quality with a reduced number of images (
i.e., 14 images) over a long time (
i.e., some years) period.
In
Figure 5a total of 7,265 PS points (1% of the whole dataset) were selected within some different rock outcrops assumed as good stable scattering areas over the Tramuntana Range. These PS targets are characterized by a good coherence (values from 0.60 to 1.00) and show a V
LOS velocity distribution ranging from −2.9 to +3 mm/yr that is consequently assumed as a reasonable motionless velocity range.
Taking into account that the V
SLOPE PS population is negatively skewed, the classification of V
SLOPE into velocity classes is based on determining the V
SLOPE stable range where a certain percentage (
i.e., 68%) of the PS population is found from zero [
21]. In the study area, the velocity ranges of the V
LOS and V
SLOPE PS distributions are ±2.3 mm/yr and −4.5 mm/yr, respectively (
Table 2).
In order to classify landslides as active or no-active, considering landslide average V
LOS and V
SLOPE velocities, the thresholds of ±3 mm/yr for V
LOS and −5 mm/yr for V
SLOPE are chosen, being slightly greater than those identified for the PS population (
i.e., rounded up to the next integer) (
Table 2).
Being this facet a critical and uncertain issue for radar landslide mapping and monitoring, another threshold is set at −10 mm/yr. This threshold, proposed by Mansour
et al. [
36] from a literature review, establish that moderate damages are produced to buildings and infrastructures when the displacement rate is between −10 and −100 mm. This threshold is adopted as the activity maps are elaborated, as final purpose, for supporting environmental planning and management strategies in order to highlight areas characterized by higher hydrogeological hazard, hence, they need to be addressed to potential damages detection and prevention.
7. Confidence Degree Evaluation
The confidence degree of landslide and ground motion activity maps can be evaluated throughout their comparison with external independent sources, such as damages inventories, in situ measurements, field checks, etc. This evaluation aims at assessing if measured displacement represents landslides dynamics. It is important to highlight that the reliability on the PSI measurement itself is not evaluated, but whether this measurement is related to landslide activity or not. Three confidence degree levels are established—low, medium, and high—depending upon ground truth data availability.
In the Tramuntana range, the confidence degree determination of the activity maps was performed through the use of Ma-10 road damage database, which is the most relevant element at risk with high exposure in the area of interest. Comparing LSI and damages to the road, it can be seen that 18 landslides intersect the Road M-10: five of them caused damages in 2007–2010, and 13 produced damages before 2007 (
Table 4).
ALOS PS data acquired in 2007–2010 monitored 11 landslides intersecting the Road, among which two phenomena produced damages in this period and nine produced damages before 2007, classified as “True” and “False” landslides, respectively (
Table 4). The two landslide polygons that intersect the road and produced damages in 2007–2010, detected by ALOS data, show an average V
SLOPE of −8 and −11 mm/yr. According to the 10 mm/yr landslide damage threshold adopted by Mansour
et al. [
36], the former would be a “True Negative” (
i.e., inactive landslide producing damages), whereas the latter would be a “True Positive” (
i.e., active landslide producing damages). Among the nine landslides detected by ALOS that did not produce damages to the road, three of them are “False Positives” (
i.e., active landslides not producing damages), and six of them are “False Negatives” (
i.e., inactive landslides not producing damages).
According to this analysis, our success rate is 63% over the 11 landslides that were monitored with ALOS data in the period 2007–2010. Consequently, a high confidence degree is attributed to seven landslides, specifically, one “True Positive” and six “False Negative” phenomena, as measured displacements (active/inactive) can be correlated with reported damages (occurred/not occurred).
A medium confidence degree is assigned to the rest of detected landslides (93 phenomena) with more than 4 PS or more than 50% of their spatial extension covered by PS pixel area. In these landslides, it cannot be confirmed that measured displacement (active/inactive) is related to landslide activity. However it is the most probable explanation since displacement measurements were gathered within the landslide boundaries.
A low confidence degree is attributed to those landslides where neither sufficient number of PS data (or no PS at all), nor ground information, provide evidence of terrain displacement. This confidence degree is also attributed to all the active clusters derived from moving PS out-points, since the mechanism of the triggering cause is unknown. Thus, a low confidence degree is allocated to 51 detected landslides with less than 4 PS or less than 50% of their surface covered by PS pixel area, and to all the 615 active clusters with an average V
SLOPE < −10 mm/yr, 18 of which intersect the Ma-10 Road (
Figure 11).
8. Discussion
Overall, the PSI analysis of ALOS PALSAR images in the Tramuntana Range reveals a good performance for mapping and monitoring landslides in a rural environment, due to the high penetration capacity of the L-band sensor and a sufficient sampling of SAR images for a three-year temporal span (2007–2010). This time interval coincides with the most critical rainy and cold period of the study area, when several landslide phenomena were triggered.
The reduced geometrical decorrelation combined with the long satellite revisiting time in L-band (
i.e., 46 days for ALOS system) leads to successful applications of L-band data especially for conventional differential interferometry, e.g., [
39–
42]. However, the use of L-band PSI allows a feasible investigation on natural areas such as Tramuntana range covered by rocky outcrops and scrubland, due to its lower volumetric and temporal decorrelation in vegetated zones with respect to C-band data, which are usually used in literature [
6,
9,
37] for landslide detection and mapping. Furthermore, L-band SAR data yields an intrinsic capability to detect faster ground movements with respect to C-band, just due to the longer orbital repeat cycle and wavelength.
The PSI post-processing stage permits to evaluate how good is the area to backscatter the radar signal, permitting to assess the suitability of the detected PS in terms of the satellite acquisition geometry, the relief, and the land-use. In the Tramuntana range, 37% of the landslides exhibit the best geometry to be monitored (i.e., average RI > 0.66; slope direction approximately parallel to the satellite LOS direction), whereas for 63% a certain geometrical distortion will be introduced in the PS measurement (average RI < 0.66). Concerning the land use, rocky outcrops with sparse vegetation and scrubland with bushes or woods show the highest number of PS, because this land use class represents the most widespread one in the Tramuntana range and because of the presence of many good stable reflectors retrieved in L-band. In particular, 180 landslides (85% of the total) are successfully monitored in a vegetated environment with a dominant presence of massif rocks (23 PS/landslide on average), confirming the high penetration capacity on the ground of ALOS PALSAR sensor.
The projection of the VLOS along the steepest slope permits to homogenize landslide velocities with different slope orientations. This procedure significantly reduces the PS population (45% in the Tramuntana range), but enriches the VLOS interpretation. The number of landslides monitored with a sufficient number of PS (> than 4 PS or more than 50% of the landslide area covered by PS) using the VLOS dataset is 96, being reduced to 80% (77 landslides) when VSLOPE dataset is considered.
However, it is worth to highlight that the downslope projection of VLOS values is only valid when the landslide movement is parallel to the slope. This requirement typically holds for planar slides with slow-moving flow, but not for rotational landslides, which generally possess a vertical movement at the crown and horizontal movement at the toe. In the Tramuntana range, the projection procedure turns out to be suitable since the most frequent phenomena, besides rockfalls, are translational slides, debris slides, and earth/debris flows.
In order to classify landslides velocity from L-band radar displacement measurements, the thresholds of ±3 mm/yr for landslide average V
LOS and −5mm/yr for landslide average V
SLOPE are chosen. These values are slightly greater than those identified for the PS population (
i.e., rounded up to the next integer) (
Table 2), and they account for 77% and 68%, respectively, of the PS population. Most of the published landslide PSI analyses have been conducted on C-band (5.3 GHz) and establish the stability V
LOS threshold at ±2 mm/yr [
9,
21]. Taking into account that the V
LOS sensitivity increases with the radar frequency, a proportionally greater threshold would be expected for L-band (1.2 GHz),
i.e., higher by a factor of more than four. ALOS L-band data have a higher standard deviation than C-band in similar conditions, but this feature also depends on the number of images, the temporal period and the quality of selected points. In this case, in the Tramuntana range L-band results are close to C-band map quality with a reduced number of images over a long acquisition time period. Moreover, it is worth to underline that ALOS PS data in the study area show a high coherence and hence a low standard deviation is achieved even within L-band dataset. Therefore, the chosen ±3 mm/yr threshold that is only 50% greater than the commonly adopted C-band threshold is assumed as a reasonable range for landslide movements based on L-band V
LOS distribution.
By the way, being a critical issue for landslide monitoring, and the adopted threshold potentially small, a −10 mm/yr threshold is set to identify active landslides. This threshold, proposed by Mansour
et al. [
36] from a literature review, establishes that moderate damages are produced to buildings and infrastructures when the displacement rate is between −10 and −100 mm/yr. The adoption of this threshold for landslide activity maps elaboration is adequate, as the final purpose of these maps is to highlight areas where potential damages could be produced, as a support to environmental planning and management strategies.
The landslide activity map integrates the PSI data, the post-processing derived products (RI, LU, C, and VSLOPE) with the landslide inventory (LSI), permitting to improve the interpretation of measured displacements with respect to landslide dynamics. In the Tramuntana range, 35 landslides were identified as inactive and 42 of them as active, from which seven of them overpass the −10 mm/yr threshold.
The “Ground motion activity map” consists of a set of active PS clusters that are not included within the mapped landslide boundaries, and that can be related to either unknown landslides or to other types of geomorphological processes. In the Tramuntana range 615 active clusters show a velocity above the 10 mm/yr threshold (18 of them intersecting the Ma-10) representing potentially hazardous areas.
The confidence degree of landslide and ground motion activity maps is evaluated throughout their comparison with external independent sources, aiming to assess if measured displacement represents the landslide dynamics or not. In the Tramuntana range, the confidence degree determination of the activity maps was performed through the use of the damages database on Ma-10 road, which is the most relevant element at risk with high exposure in the area of interest. A high confidence degree is attributed to seven landslides as measured displacements (active/inactive) can be correlated with reported damages (occurred/not occurred). Particular attention is paid to Estellencs landslide, where the acceleration trend observed between August 2009 and March 2010 could be related to the acceleration of the moving mass before the total failure of the road occurred on 8 March 2010. A medium confidence degree is assigned to 93 landslides as PSI measured displacement (active/inactive) coincide with the landslide boundaries, but there are no external sources to confirm landslide activity. A low confidence degree is attributed to 51 landslides where an insufficient number of PS (or no PS at all) provides evidence of terrain displacement. This confidence degree is also attributed to the 615 detected active clusters since their mechanism of the triggering cause is unknown.
Overall, the landslide database on the Tramuntana range (
Figure 1e) has been enhanced by detecting and analyzing landslide phenomena and further unstable areas (
Figure 11) characterized by higher motion rates. The reliability of the delivered maps has been assessed through the confidence degree evaluation, making this “added value” interesting and useful for the policy decisions and strategies concerning risk mitigation and urban-environmental design.
9. Conclusions
In this paper landslide and ground motion activity maps are elaborated through the use of PS (Persistent Scatterers) ALOS data in L-band, over the Tramuntana Range in the island of Majorca (Spain), extensively affected by mass movements across time, especially during the last years.
The landslide activity map provides, for every monitored landslide, an assessment of the PS visibility according to the relief, land use and satellite acquisition parameters. The average VLOS and VSLOPE velocities are also attributed to every landslide, providing an estimate of their state of activity and their potential to cause damages. Additionally, a ground motion activity map is also generated, based on active PS clusters not included within any mapped landslide phenomenon. Finally, a confidence degree evaluation is carried out to attest the reliability of measured displacements to represent landslide dynamics. In the Tramuntana range 42 landslides were identified as active (VSLOPE < −5mm/yr) and seven of them with a potential to produce moderate damage (VSLOPE < −10mm/yr). A high confidence degree was attributed to seven landslides, a medium one to 93, and a low confidence degree to 51 landslide phenomena and to 615 active clusters, 18 of them intersecting the Ma-10 and thus representing further potentially hazardous areas.
Overall, the outcomes of this work reveal the usefulness of landslide activity maps for environmental planning activities, being exportable to other radar data and different geomorphological settings.