Spaceborne, UAV and ground-based remote sensing techniques for landslide mapping, monitoring and early warning

Background

The current availability of advanced remote sensing technologies in the field of landslide analysis allows for rapid and easily updatable data acquisitions, improving the traditional capabilities of detection, mapping and monitoring, as well as optimizing fieldwork and investigating hazardous or inaccessible areas, while granting at the same time the safety of the operators. Among Earth Observation (EO) techniques in the last decades optical Very High Resolution (VHR) and Synthetic Aperture Radar (SAR) imagery represent very effective tools for these implementations, since very high spatial resolution can be obtained by means of optical systems, and by the new generations of sensors designed for interferometric applications. Although these spaceborne platforms have revisiting times of few days they still cannot match the spatial detail or time resolution achievable by means of Unmanned Aerial Vehicles (UAV) Digital Photogrammetry (DP), and ground-based devices, such as Ground-Based Interferometric SAR (GB-InSAR), Terrestrial Laser Scanning (TLS) and InfraRed Thermography (IRT), which in the recent years have undergone a significant increase of usage, thanks to their technological development and data quality improvement, fast measurement and processing times, portability and cost-effectiveness. In this paper the potential of the abovementioned techniques and the effectiveness of their synergic use is explored in the field of landslide analysis by analyzing various case studies, characterized by different slope instability processes, spatial scales and risk management phases.

Results

Spaceborne optical Very High Resolution (VHR) and SAR data were applied at a basin scale for analysing shallow rapid-moving and slow-moving landslides in the emergency management and post- disaster phases, demonstrating their effectiveness for post-disaster damage assessment, landslide detection and rapid mapping, the definition of states of activity and updating of landslide inventory maps. The potential of UAV-DP for very high resolution periodical checks of instability phenomena was explored at a slope-scale in a selected test site; two shallow landslides were detected and characterized, in terms of areal extension, volume and temporal evolution. The combined use of GB-InSAR, TLS and IRT ground based methods, was applied for the surveying, monitoring and characterization of rock slides, unstable cliffs and translational slides. These applications were evaluated in the framework of successful rapid risk scenario evaluation, long term monitoring and emergency management activities. All of the results were validated by means of field surveying activities.

Conclusion

The attempt of this work is to give a contribution to the current state of the art of advanced spaceborne and ground based techniques applied to landslide studies, with the aim of improving and extending their investigative capacity in the framework of a growing demand for effective Civil Protection procedures in pre- and post-disaster initiatives. Advantages and limitations of the proposed methods, as well as further fields of applications are evaluated for landslide-prone areas.

Landslides play an important role in the evolution and shaping of aerial/subaerial landscapes (Brunetti et al., 2015), representing a major cause of loss of life, injuries, property damage, socio-economic disruption and environmental degradation (WP/WLI, 1993; Canuti et al., 2004; Petley et al., 2005; Petley, 2012), especially if they are associated with other natural hazards (like earthquakes, volcanic eruptions, meteorological events and wildfires). Because of such habitual combinations, reliable numbers for the social impact only due to landslides are difficult to obtain on a global scale and the economic losses are certainly underestimated (or not quoted at all). This general condition often contributes to reducing the concern individuals and authorities have about landslide risk (Kjekstad, and Highland 2009). Although in most of the disaster-prone areas the consideration of the social-cultural and socio-economic conditions in relation to their physical safety is still very confused, the application of appropriate technologies for landslide detection, monitoring and early warning systems are increasingly considered crucial by local authorities in reducing the risk of landslide disasters. EO from space has found many uses in the natural sciences, but it is only in the last decades that technological advances have also extended to landslides analysis (Mantovani et al., 1996; Ferretti et al., 2001; Canuti et al., 2004; Metternicht et al., 2005; van Westen et al., 2008; Casagli et al., 2010; Martha et al., 2010; Guzzetti et al., 2012; Lu et al., 2012; Tofani et al., 2013a). Nowadays rapid advances are making EO techniques more effective for landslide detection, mapping, monitoring and hazard assessment. Applications are originating from nearly all types of sensors available today (Tofani et al. 2013b). Rapid developments in this field are fostered by the very high spatial resolution obtained by optical systems (currently in the order of tens of centimeters) and by the launching of SAR sensors, purposely built for interferometric applications with revisiting times of few days, such as TerraSAR X and COSMO-SkyMed (Tofani et al., 2013a). Landslide detection and mapping benefit from both optical (Hervas et al. 2003, Cheng et al., 2004, Marcelino et al., 2009, Martha and Kerle 2012, Lu et al., 2011) and radar imagery (Singhroy, 1995; Fruneau et al., 1996; Massonnet and Feigl, 1998; Kimura and Yamaguchi, 2000; Hilley et al., 2004; Hanssen, 2005; Colesanti and Wasowski, 2006; Meisina et al., 2008; Herrera et al., 2009; 2011; Bardi et al., 2014; Crosetto et al., 2016) to study slow moving landslides. The ability to make numerous point measurements of displacement over the landslide body allows one the detection and mapping of the actively deforming slopes (e.g. Righini et al., 2012), the characterization and monitoring of landslide mechanism (Tofani et al., 2013b) and, through the analysis of time series of deformation, the identification of velocity changes in the landslide evolution (Berti et al., 2013), as well as the modeling of large slope instability (Berardino et al., 2003). Advanced terrestrial remote sensing technologies, such as GB-InSAR, TLS, IRT and digital photogrammetry (DP) are nowadays applied in the field of slope instability detection, mapping and monitoring, for short/long term landslide management (real time, near real time and deferred time) (Lillesand et al., 2014). They are characterized by operational efficiency and accuracy of data not reached by traditional methods: high-resolution acquisition, multifunction versatility, device portability, low cost sensors, easy and fast data processing. Such equipment allows for systematic and easily updatable acquisitions of data that may also enhance the implementation of effective early warning systems at slope scale. In this paper the potential of the abovementioned remote sensing techniques (both spaceborne and ground-based), and their applications for landslide detection and mapping are evaluated.

The presented techniques are described by means of their main technical features and applicability in different observed scenarios, typology of landslide (Cruden and Varnes, 1996; Hungr et al., 2014) and geomorphological setting. Some case studies are also shown and discussed in order to exhibit good practices in landslide characterization and prediction by means of different techniques and sensors in synergic action. The main advantages and disadvantages of the presented techniques are described in the text and in a tabular form.

Spaceborne platforms

Optical VHR imagery

The most important active optical satellites are reported in Fig. 1. Optical data are usually used for landslide detection and mapping through visual inspection or analytical methods (Metternicht et al., 2005; Fiorucci et al., 2011; Parker et al. 2014; Guzzetti et al., 2012; Mondini et al., 2014). For example, several optical derivative products (panchromatic, pan sharpen, false colour composits, rationing) can help in visual interpretation for landslide mapping (Casagli et al., 2005; Marcelino et al., 2009; Ma et al. 2016). In image fusion procedures, multispectral channels, characterized by a coarser spatial resolution than the panchromatic, are downscaled through analytical models based on the panchromatic-derived spatial information (Eyers et al., 1998; Chini et al., 2011; Martha and Kerle, 2012; Kurtz et al., 2014). The False Colour Composites (FCCs) of the VHR images are often used to discriminate lithologies or terrain having different characteristics (weathering, water content, vegetation cover) (Ciampalini et al., 2012; Lamri et al., 2016).

Active optical and SAR satellites for landslide mapping and monitoring. The numbers on the right of the figure report the revisting time of each satellite. RCM: Radarsat constellation mission, CSK: COSMO-Skymed, CSK–SG: COSMO-Skymed Second Generation

Full size image

The Normalized Difference Vegetation Index (NDVI) derived from optical images, is another index widely used to map landslides by means of evaluating the vegetation cover rate (Lin et al., 2004). Higher values of NDVI can be related to a wide vegetation cover, whereas lower values can represents areas affected by landslides. Furthermore, multispectral images can be enhanced to detect landslides by means of analytical methods based on the spectral characteristics of the land surface and automatic approaches focus on the classification of image pixels (Martha et al., 2010; Mondini et al., 2011). Few studies have described the use of hyperspectral data for recognition and classification of landslides based on Earth surface characteristics since most of the hyperspectral satellite sensors are still under development (Scaioni et al., 2014).

SAR data

The family of SAR satellite sensors (Fig. 1) orbits the Earth at an altitude ranging from 500 to 800 km, following sun-synchronous, near-polar orbits, slightly inclined with respect of Earth meridians. The most commonly used bands in SAR applications are C-band (5–6 GHz, ~5,6 cm wavelength), X-band (8–12 GHz, ~3,1 cm wavelength) and L-band (1–2 GHz ~23 cm wavelength) with a temporal resolution depending on the satellite revisiting time (Fig. 1). A SAR image is composed of pixel characterized by amplitude and phase values. Phase values of a single SAR image is partly depends on the sensor-target distance and is the key element to detect ground displacement. SAR Interferometry is the technique focused on the measure changes of signal phase over time through the analysis of at least two SAR images (Fruneau et al., 1996; Singhroy et al., 1998). A suitable approach to exploit phase variation between two consecutive radar images acquired over the same target is the Differential Interferometric SAR (D-InSAR) (Bamler and Hartl, 1998; Rosen et al., 2000). Geometrical and temporal decorrelation and atmospheric effects caused by the variation of the phase reflectivity value of some radar targets reduce the reliability of the D-InSAR technique (Berardino et al., 2002). In order to overcome these limitations InSAR-based information can be enhanced through multi-temporal interferometric techniques (MIT), based on analysis of long stacks of co-registered SAR imagery (Ferretti et al. 2001; Crosetto et al, 2016). In the past years, several MIT approaches have been developed such as: the Permanent Scatterers Interferometry, named PSInSAR™ (Ferretti et al., 2011; Colesanti et al., 2003), the SqueeSAR™ (Ferretti et al., 2011), the Stanford Method for Persistent Scatterers StaMPS (Hooper et al., 2004; Hooper et al., 2007), the Interferometric Point Target Analysis IPTA (Werner et al., 2003; Strozzi et al., 2006), the Coherence Pixel Technique CPT (Mora et al., 2006), the Small Baseline Subset SBAS (Lanari et al., 2004; Berardino et al., 2003), the Stable Point Network SPN (Casu et al., 2006; Crosetto et al., 2008), the Persistent Scatterer Pairs PSP (Herrera et al., 2011) and the Quasi PS technique QPS (Costantini et al., 2008). Signal analysis of a network of coherent radar targets (Permanent Scatterers, PS) allows estimating occurred displacement, acquisitions by acquisition. Line of Sight (LOS) deformation rate can be estimated with an accuracy theoretically better than 0.1 mm/yr. Each measurement is referred temporally and spatially to a unique reference image and to a stable reference point. MIT analysis is designed to generate time-series of ground deformations for individual PS, assuming different types of deformation models (e.g., linear, nonlinear or hybrid). .In the field of landslide investigations the potential of SAR data has been exploited at different scales: from national (Adam et al., 2011) to regional (Meisina et al., 2008; 2013; Ciampalini et al. 2016a, b) basin (Lu et al., 2012) slope (Frodella et al., 2016) and building scale (Ciampalini et al., 2014; Bianchini et al., 2015; Nolesini et al., 2016), as well as in different phases of landslide response (Canuti et al., 2007) and Civil Protection practice (Farina et al., 2008). Other application fields include subsidence phenomena (Raspini et al., 2012; 2014; Rosi et al. 2014; 2016), earthquakes (Bűrgmann et al., 2005; Sousa et al., 2010) and volcanic activities (Hooper et al., 2004; Vilardo et al., 2010; Parker et al., 2014).

UAV and Ground-based methods

UAV-DP

DP is a well-established technique for acquiring dense 3D geometric information in slopes from stereoscopic overlaps of photo sequences captured by a calibrated digital camera (Chandler, 1999; Lane et al., 2000; Sturzenegger and Stead, 2009; Zhang et al., 2004). During past few years, with the rapid development of DP techniques and the availability of ease-using, focusable and relatively cheap digital cameras, this technique gained wide applications in many fields, such as 3D building reconstruction, heritage protection and landslides studies (Grussenmeyer et al., 2008; Scaioni et al., 2015; Fan et al., 2016). In this latter field, depending on the camera lens-setting, DP can be divided into two fields of activity (Gopi, 2007): far range, usually more exploited for landslide characterization and general mapping (Wolter et al., 2014), and close range, having a wide use in high precision metrological and deformation monitoring applications (Cardenal et al., 2008; Scaioni et al., 2015). More recently the combination of rapid development of low cost and small UAVs and the improvements of conventional sensors in terms of cost and size, led to new, promising scenarios in environmental remote sensing, surface modelling and monitoring (Colomina and Molina, 2014; James and Robson, 2012; Remondino et al., 2011; Eisenbeiss and Sauerbier, 2011).

GB-InSAR

GB-InSAR system consists of a computer-controlled microwave transceiver, characterized by a transmitting and receiving antennas, which by moving along a mechanical linear rail is capable to synthesize a linear aperture along the azimuth direction (Tarchi et al., 1997; Rudolf et al., 1999; Pieraccini et al., 2002). The obtained SAR image contains amplitude and phase information of the observed scenario backscattered echo in the acquiring time interval (from few to less than 1 min with the most modern systems) (Luzi et al., 2004; 2010; Monserrat et al., 2014). In a GB-InSAR interferogram the displacement obtained from the phase difference calculation can be represented in 2D maps, in which the chromatic scale covers a total value corresponding to half of the wavelength used. However, since the phase is periodic, it cyclically assumes the same values crating image-interpreting problems. This issue, known as phase ambiguity, and can be solved through interpretation based on field geological knowledge or by adopting apposite phase unwrapping algorithms (Ghiglia & Romero, 1994), which count the number of cycles performed by the wave obtaining cumulated displacement maps. Given the relative short distances at which GB-InSAR apparatuses usually operate (typically less than 3 km), they work in Ku band (1.67–2.5 cm). The main research applications of GB-InSAR soon became focused on slope monitoring (Tarchi et al., 2003; Pieraccini et al., 2002; 2003), for civil protection purposes (Del Ventisette et al., 2011; Intrieri et al., 2012; Bardi et al., 2014; 2016; Lombardi et al., 2016) and, more recently, for mining safety (Farina et al., 2011; Severin et al., 2014). Other fields include volcanoes monitoring (Di Traglia et al., 2013; 2014a; 2014b; Intrieri et al., 2013; Nolesini et al., 2013; Calvari et al., 2016), cultural heritage sites (Tapete et al., 2013; Pratesi et al., 2015; Nolesini et al., 2016; Frodella et al., 2016), glaciers and snowpack sinkholes (Intrieri et al., 2015).

TLS

A TLS device produces and emits a beam characterized by a directional, coherent and in-phase electromagnetic radiation (Jaboyedoff et al., 2012). The laser scanner by measuring with high accuracy (millimeter or centimeter) the back-scattered laser signal, is capable of obtaining the exact position of a mesh of points (point cloud), characterized by (x, y, z) cartesian coordinates (Slob et al., 2002; Frohlich and Mettenleiter, 2004; Turner et al., 2006; Slob et al., 2007). The device high acquisition rate (up to hundreds of thousands points per second) makes the detailed 3D shape of the object available in a short operating time. By defining the coordinates of specific laser reflectors within the surveyed area through a Differential Global Positioning System in Real Time Kinematic mode (DGPS-RTK; Morelli et al., 2012; Tapete et al., 2015; Pazzi et al., 2016), it is possible to link the obtained high-resolution 3D surface digital model to a global reference system. In landslide studies TLS has been increasingly used for the geometrical and geostructural characterization and unstable rock cliffs monitoring applications (Abellán et al. 2006; 2011; Jaboyedoff et al., 2007; Ferrero et al., 2009; Oppikofer et al. 2009; Gigli et al. 2014a, b, c). Thanks to the high resolution of the laser scanning survey it is also possible to extract even the smallest features, such as the structural crack pattern, the crack opening direction (Gigli et al., 2009; 2012), and the orientation of critical discontinuities within the rock mass (Gigli and Casagli, 2011; Gigli et al., 2014b). Furthermore, this technique is capable of measuring ground 3D temporal displacements by comparing sequential datasets of the same scenario (Rosser et al., 2005; Abellán et al., 2011). The intensity data can also provide some information about the type of material and the soil moisture content of the targets, which can add information regarding the landslide main geomorphologic features (Voegtle et al., 2008; Franceschi et al., 2009).

IRT

IRT is the branch of remote sensing dealing with measuring the radiant temperature of Earth’s surface features from a distance (Spampinato et al. 2011). The product of an infrared thermographic survey is a pixel matrix (thermogram), collected through the thermal camera array detector (Maldague, 2001), which following the correction of the sensitive parameters (object emissivity, path length, air temperature and humidity) represents a radiant temperature map of the investigated object. The presence within the observed surface of fractures, subsurface voids, moisture and seepage zones, will influence the material thermal characteristics (density, thermal capacity and conductivity) modifying its heat transfer (Teza et al., 2012). Therefore, the presence of an inhomogeneity within the observed scenario will be displayed in the corresponding radiant temperature map as an irregular thermal pattern with respect to the surroundings (a “thermal anomaly”) (Frodella et al., 2014b). In recent years IRT has undergone a significant increase of applications in the field of geosciences (Spampinato et al. 2011), nevertheless in the study of slope instability processes it is still experimentally used, except for a few interesting experimental studies (Wu et al., 2005; Baroň et al. 2012; Frodella et al., 2014b). In particular, IRT (often coupled with laser scanning) is applied with the following purposes: i) obtain information about the rock mass fracturing (Squarzoni et al. 2008); ii) detect shallow surface weakness in rock walls (Teza et al. 2012); iii) perform rockfall/slide susceptibility assessment (Gigli et al. 2014a, c; Teza et al. 2015); iv) map ephemeral drainage patterns (Frodella et al., 2014a; 2015); v) integrate traditional geo-structural and geomechanical surveys (Mineo et al., 2015; Mineo and Pappalardo 2016; Pappalardo et al., 2016).

In this section, the potential of the presented techniques and their synergic use is explored for the detection, mapping, and monitoring of landslides, through various case studies characterized by different types and scales of instability phenomena, hazardous scenarios and operational modes (Fig. 2). The employed remote sensing systems are listed in Table 1.

a Landslide case studies location; b Giampilieri (debris flows); c Nebrodi area (complex landslides); d Ricasoli (shallow landslides); e Western Elba Island (unstable rock masses); f San Leo (collapsed rock cliff); g Santa Trada (translational slide)

Full size image

Spaceborne platforms

The study areas for spaceborne applications comprise two different mountain chain sectors (the Peloritani and the Nebrodi mountains; Sicily Island, Southern Italy; Fig. 2), in which the geological features are characterized by the typical features of recently uplifted areas, developed on a crystalline basement with steep slopes and shallow clayey soil cover. In the late afternoon of October 1^st 2009, an intense storm affected the area between the Peloritani Mountains ridge and the Ionian coastline (Ciampalini et al., 2015a; Del Ventisette et al., 2012), where the main villages are located. During the same night, the persisting rainfall triggered more than 600 landslides, such as shallow soil slides and debris flows, on an area of about 50 km². The assessed number of fatalities caused by landslides and inundation was 37 (including 31 deaths and 6 missing persons), with 122 injured people and 2019 evacuated people (Ardizzone et al., 2012; Del Ventisette et al., 2012; Raspini et al., 2013); the worst damages were reported in the village of Giampilieri (Fig. 2). Furthermore, between 2009 and 2010, following heavy and persisting heavy rainfall, several municipalities in the Nebrodi Mountains were strongly affected by several complex, rotational and deep-seated landslides which damaged buildings and infrastructures (Ciampalini et al., 2014; Bardi et al., 2014; Bianchini et al., 2015; Ciampalini et al. 2015a, b; 2016a,b).

Object-Oriented Analysis (OOA) for mapping shallow rapid-moving landslides: the Giampilieri case study

The purpose of this case study is to introduce a new approach for a rapid mapping of newly-triggered landslides using an objected-oriented change detection technique. The methodology aims at a semi-automatic and rapid analysis with a minimum of operator involvement and manual analysis steps. Compared to conventional approaches for landslide mapping, this approach benefits from (i) an image segmentation with problem-specified scale optimization, and (ii) a multi-temporal analysis at object level with several systemized spectral and textural metrics. This procedure has been applied to the two of the most damaged areas of Giampilieri, including a training area (ca. 1.8 km²) for algorithm development, and a larger independent testing area (ca. 8.1 km²). The latter allows the robustness and transferability of the algorithm (without any change of ruleset and threshold) and the corresponding accuracy to be assessed by comparison with a manually mapped landslide inventory prepared from fieldworks and subsequent modifications from image interpretation. Two Quickbird images acquired on September 6^th 2006 and October 8^th 2009, with 0.3% and zero cloud cover respectively, were used in the study (Table 1). The application with the optical data is based upon the OOA (Lu et al., 2011). OOA is mainly dealing with the measuring unit of ‘object’, which can be defined as ‘individually resolvable entities located within a digital image which are perceptually generated from high-resolution pixel groups’ (Hay et al., 2003). Detailed information on the methodology and the algorithm developed can be found in Lu et al. (2011). The algorithm developed based on the training area was directly applied in the testing area.

The final outputs for the testing area are shown in Fig. 3: in order to evaluate the accuracy of this approach, OOA-derived landslides were compared with a manually-mapped landslide inventory. The accuracy assessment was carried out for the number and the spatial extent of mapped landslides. For the spatial extent of landslides a user’s accuracy of 75.9% and a producer’s accuracy of 69.9% were achieved. In terms of the number of landslides, user’s and producer’s accuracies of 81.8 and 69.5%, respectively, were reached.

The used Quickbird imagery in the optical VHR analysis of Giampilieri area (Peloritani mountains): a pre-event QuickBird imagery; b post-event QuickBird imagery (false color 4-3-2); c The result of OOA landslide mapping in the independent testing area (yellow areas = mapped shallow landslides)

Full size image

Detection and mapping of slow-moving landslides with SAR data: the Sicily case studies

Spaceborne SAR analysis of ground deformation in the Peloritani (nearby the Giampilieri village) and Nebrodi area (Fig. 2) was performed using the SqueeSAR approach (Ferretti et al., 2011). The SqueeSAR algorithm has been applied to C-band SAR dataset acquired by ERS (08/09/1992 - 24/11/2000) and Envisat (22/01/2003-20/05/2009) missions along ascending orbits. Following the approaches proposed by Farina et al. (2008) and Bianchini et al. (2012), deformation measurements extracted by means of SqueeSAR technique have been coupled and integrated with thematic maps (topographic and geological maps), optical data (ortophoto, optical satellite VHR images and multi-temporal aerial photos) and available landslides inventory maps, in order to identify the areas characterized by high hydro-geological hazard (hotspot mapping), related to the occurrence of extremely and very slow moving landslide (according to the classification of Cruden & Varnes, 1996). Twenty-six sites have been identified, for which landslides have been detected and mapped (Fig. 4). On the basis of available multi–interferometric data these sites were assessed as the most critical in terms of hydro-geological hazard, both for the type of instability detected and/or the extent of the mapped phenomena and/or the measured deformation velocities and/or the presence of elements at risk. In the inset of Fig. 4 the hotspot analysis for the area of the village of Croce (identified with number 9) is reported. The area is located on the right bank of the Guidomandri creek and is characterized by the presence of several elements at risk, including, beside the village itself, isolated buildings and minor settlements. Pre-existing landslide inventory maps do not report slope instability in the study area.

Location of the twenty-six sites in the Nebrodi area which are characterized by high hydro-geological risk (hotspot mapping) according to the D-InSAR analysis. In the inset the landslide mapping of the village of Croce is reported

Full size image

The SqueeSAR results also show a large sector of the slope characterized by a displacement with velocity ranging from 1.6 to 4.8 mm/yr Envisat dataset (2003–2009). Photo-interpretation of stereoscopic colour images (1:3500 scale) and analysis of information provided by SqueeSAR results allow to detect and map a large complex system of active continuous slides affecting the area. Such deformation rates do not pose threat to population, but can cause, persisting for many years, damages to buildings and manmade infrastructures. Cracks and damages have been surveyed during field validation, which also helped to confirm the presence and the extension of the active movements through the identification of tension cracks, scarps and counterscarps.

After the 2009 and 2010 events that affected the Nebrodi Mountains, the SqueeSAR technique was applied to characterize the triggered hillslope phenomena both at the basin and at the local scale (Fig. 5). At the basin scale, SqueeSAR PSI data was applied to update the available Landslide Inventory Map (LIM) including information on typology and state of activity of each identified landslide. The updating procedure has been performed using: (i) radar interpretation of four different available SAR datasets; (ii) photo-interpretation of 1:33000 scale aerial photographs flown in 1954, 1955 and 2005; and (iii) field surveys. InSAR displacement measurements were acquired in different periods (2006–2009, RADARSAT-1 scenes and 2011–2012, COSMO-SkyMed images). Considering the limitation of the adopted technique, the updating of the pre-existing LIM was limited to the extremely slow and very slow moving landslides (faster phenomena have been excluded due to their rapid kinematics).

Ground deformation velocity maps in the surroundings of Militello Rosmarino town (Nebrodi) obtained using Radarsat-1 (a) and COSMO-SkyMed (b) PSI data. Pre-existing LIM (c) and updated LIM (d)

Full size image

The new LIM (Fig. 5) includes 566 events: 15 (2.7%) rockfalls and topples, 136 (24.0%) complex landslides, 188 (33.2%) flows and 227 (40.1%) slides, covering an area of 74.1 km². The comparison between the pre-existing and the new LIMs has been performed using three classes: (i) confirmed, (ii) enlarged with respect to the pre-existing LIM and (iii) new (landslide not included in the pre-existing LIM). This approach led to the enlargement of 120 events (21.2%) of the pre-existing LIM, to the confirmation of 155 events (27.4%), and to the recognition of 291 (51.4%) new phenomena.

UAV and Ground based methods

UAV-DP for landslide characterization and mapping: The Ricasoli case study

A periodical check was performed in the Ricasoli village (Upper Arno river Valley,Tuscany, Italy; Fig. 2), in order to evaluate the potential of UAV-DP to characterize and to monitor landslides. In particular, a multitemporal photogrammetric survey, carried out for the northern slope of Ricasoli, are compared to define at very high resolution, morphologic features of the slope and their evolution in time. The survey was performed using a multicopter drone (Saturn) with an innovative perimetric chassis, fully designed, built and patented by the Department of Earth Science of the University of Florence (Fig. 6). The images were processed using Agisoft Photoscan Professional (Agisoft LLC, 2016) software and the resulting data were implemented in a GIS environment using the ESRI ArcGIS package. Three 3D point clouds, acquired at a few months one to each other and filtered in order to remove all the points processed on trees and high vegetation, were used to build high-resolution DTMs (0.05 m/pix) (Fig. 6). The DTMs were compared to detect any morphological change between the three acquisitions, to characterize the landslides and, in addition, to precisely point out features as indicators of landslide-prone areas on the slope (Fig. 6). As a result, two landslides were detected and characterized, in terms of areal extension, volume and temporal evolution. The overall extent and volume of the mass movements detected in Ricasoli are summarized in Table 2. The slope is currently being monitored by performing repeated aerial surveys, and the assessed landslide evolution is being used as an input for current mitigation works planning.

a, b and c Orthophotos acquired during the three surveys in Ricasoli village by means of a special camera mounted on the Saturn multicopter drone (d) and (e and f) the differences in height calculated using the high resolution DTMs. The main scarps and landslide features are detected and analyzed based on the results of the DTM comparison

Full size image

TLS and IRT for risk scenario assessment: the Elba Island case study

The investigated area is located on the western Elba Island coastline (Central Italy), along a 250 m stretch of a local panoramic roadway (provincial roadway n° 25) (Fig. 2). The area is characterized by very steep rock slopes overlooking the roadway, which due to their complex geostructural setting and degree of fracturing (Gigli et al., 2014a), in 2009 underwent the detachment of rock mass portions and rock debris. In order to define the risk scenarios for the roadway transportation security conditions, the slope instability occurrences were investigated through a methodology based on the integration of accurate geological and geomechanical field surveys and terrestrial remote sensing techniques, such as TLS and IRT (Frodella and Morelli, 2013; Gigli et al., 2014a). IRT surveys in particular were carried out in correspondence of rock mass most critical sectors, in order to detect thermal anomalies connected to open fractures, water seepage and moisture zones, validate the unstable block volume calculation, and rapidly assess the hydraulic conditions along the more critical rock mass discontinuities. The obtained TLS 3D surface model contributed to characterize the morphological variability of the investigated area: a rough morphology, characterized by creek erosion gullies isolating jutting rock mass portions (Fig. 7c). Figure 6c shows the stereographic projection of the collected field survey structural data: five main discontinuity sets were identified, JN3 set in particular, including high persistent decimetric-spaced discontinuity planes (=exfoliation joints, EJ) dipping parallel with respect to the slope, represent slipping planes isolating large rock mass portions.

a High-definition 3D surface of the western Elba coastline (dots mark the different TLS scan positions, the square delimitates semiautomatic geomechanical surveyed area; b 3D representation of all the joint sets extracted; stereographic projection of discontinuity poles and modal planes of the main sets collected in the investigated area by means of traditional field surveys (c), and the semi-automatic analysis (d) (modified after Gigli et al., 2014a)

Full size image

Furthermore a semi-automatic geo-structural survey was performed by means of a Matlab tool (DiAna = Discontinuity Analysis; Gigli and Casagli 2011), on a limited sector of the rock mass not covered by nets, rock bolts, and fences. Figure 7d reports the poles of the semi-automatically extracted discontinuities (labeled from D1 to D7, and represented in 3D in Fig. 7b). Given the geological setting of the investigated area, and the most probable detected failure mechanism occurring (planar failure along JN3 discontinuity set), an iterative procedure was applied with the aim of identifying the maximum credible scenario. A Matlab routine was built for this purpose by moving on the 3D surface a plane with the same orientation of JN3 set. By selecting a volume threshold value of 1000 m³, three protruding rock masses were detected and labeled from north to south as M1, M2, and M3 (Figs. 7a and 8a, b, c).

Optical images of the unstable rock masses of the western Elba coastline (a = M1; c = M2; e = M3); related 3D digital model with the detected basal and lateral slipping planes (b, d, f); mosaicked thermograms from IRT (g = M1, h = M2, and i = M3) acquired around 1 p.m., November 2011 (dotted lines mark the basal slipping planes; white squares on the thermogram allow a comparison with the correspondent sectors in the optical images, acquired by the built-in digital camera) (modified after Gigli et al., 2014a)

Full size image

M3 rock mass, in addition to the basal slipping plane, is also delimited southeastward from the stable portion of the rock slope by a second sub-vertical plane (belonging to JN2 set in Fig. 5c and D3 in Fig. 5d). The obtained surface temperature maps highlighted warm thermal anomalies connected to air circulation were detected in correspondence of the open portions of the JN3 discontinuities delimiting the detected M1, M2, and M3 masses (Fig. 8g, h, i). The abovementioned discontinuities detected on the thermograms follow closely the EJ basal planes; this interpretation was strengthened by the comparison of the thermograms with the optical images that confirmed no evidence of water flow along the detected discontinuities. For these reasons, dry conditions were diagnosed for all M1, M2, and M3 basal slipping planes, and the absence of water pressure was considered in the carried out stability analysis (Gigli et al., 2014a). The resulting assessed rock mass volumes (expressed in cubic meter) are 3706 (M1), 4359 (M2), and 1293 (M3) respectively (Fig. 8c, d, e).

Long-term monitoring of collapse-affected rock wall by means of GB-InSAR, TLS and IRT: The San Leo case study

The town of San Leo is located in the southwestern sector of the Emilia Romagna Region (northern Italy; Fig. 2), on top of a limestone isolated rock massif overlying clayey slopes, which is historically affected by instability phenomena. On February 27^th, 2014 an entire portion of the rock plate north-eastern sector collapsed, causing a huge rockfall. Following the event a GB-InSAR monitoring activity coupled with TLS surveys was carried out, in order to manage the post-event emergency phase and evaluate the residual risk (Frodella et al., 2016). Furthermore, IRT surveys were performed in order to integrate the TLS and GB-InSAR data for the rock wall characterization. The obtained 3D terrain model revealed a rock wall surface, characterized by criticalities such as overhanging sectors, ledges and niches (Fig. 9a). 3D temporal variations of the terrain model were detected by comparing sequential datasets acquired in the carried out different laser scanning surveys; the resulting 3D rock wall temporal variations, from March 7^th 2014 to December 18^th 2014 (Fig. 9b, c). The deformational field analysis provided evidence of an ongoing rock block toppling (with an estimated volume of 450 m³), which displacement evolution reached values ranging from 12 to about 50 cm (Fig. 9b, c). The scan comparison also provided the detection of minor rockfall phenomena (areas colored in blue in Fig. 9b, c) which volumes are listed in Table 3.

a High-resolution 3D surface of the surveyed rock wall in Sal Leo rock cliff (red dashed square corresponds to the area affected by block detachments). b Comparison between March 7^th 2014-April 9th 2014 TLS scans: white oval showing the first monitored rock block detachment; orange-yellow areas enhance rock wall sector characterized by displacements. c Comparison between March 7^th 2014-December 18th 2014 scans, enhancing the occurred detached rock block sectors (in blue); yellow ovals enhance the minor block detachments. d Correspondent sectors in optical image (black ovals enhance major detachments; dashed line delimits the rock wall displaced sector); (after Frodella et al., 2016)

Full size image

The TLS 3D model was merged with the GB-InSAR data obtaining a 3D GB-InSAR cumulative displacement map, which allowed to read detected LOS displacements directly on the observed scenario 3D representation, and therefore both to better localize the most critical areas, and compare the different techniques displacement data (Fig. 10).

GBInSAR and TLS data integration. a San Leo rock wall March 7th 2014 – 2015 3D GB-InSAR cumulative displacement map (after Frodella et al., 2016). b Thermogram acquired during April 9^th 2014 (after Barla et al., 2016); c correspondent photo (seepage sectors 1-4 are characterized by lower temperatures, due to local rock wall surface cooling caused by water evaporation)

Full size image

The GB-InSAR data acquired during the first monitoring year allowed to assess a general stability of the rock cliff and the observed town structures, and to detect critical areas, corresponding to: i) a detensioned rock block located at the foot of the monitored rock wall central sector (confirming the TLS analysis of an ongoing rock block toppling); and ii) the rockfall deposits (metric and decametric size boulders and blocks in a coarse sandy-clayey matrix, corresponding to the maximum recorded cumulative displacement in the investigated area. Surface temperature maps collected on April 9^th 2014 (following a period characterized by local intense rainfall) allowed to detect widespread seepage sectors in correspondence of a rock mass key discontinuity, corresponding to a high persistent normal fault dissecting the whole rock massif (oval 1 in Fig. 10b, c). In this geological, morphological and structural context discontinuities affected by seepage represent potential criticalities with respect to instability phenomena, as confirmed by minor seepage sectors (ovals 4 in Fig. 9b, c), which are located in correspondence of rock wall sectors affected by a widespread fracture network and block detachments (ovals 1–2 in Fig. 10c, d).

Short-term GB-InSAR monitoring for emergency management: Santa Trada case study

The Santa Trada landslide (Calabria Region, Southern Italy; Fig. 2) occurred on January 30^th 2009, after a period characterized by heavy rainfall (Del Ventisette et al., 2011), putting at high risk a viaduct sector along the A3 national motorway, and could have also dammed the stream below. It is a 100 m high, 90 m wide translational slide (estimated thickness is between 3 and 5 m), developed in sand and conglomerates originating from metamorphic weathered rocks (Fig. 11). For safety reasons this tract of the motorway was closed to traffic and on 31^st January a ground-based had been installed. Already on February 2^nd 2010, after a short-monitoring campaign the motorway was partially reopened, thanks to the structure stability assessment performed by means of the first GB-InSAR monitoring data (Figs. 11 and 12). During the mid-term monitoring campaign (lasted until April 24^th 2010) two approaches for calculating the interferograms and displacement maps were adopted:

a Photo of the Santa Trada landslide, 1^st February 2009; b cumulated displacement map calculated using GB-InSAR data with incremental approach spanning from 2nd February to 29th April 2009. The letters indicate the corresponding points between the two images (after Del Ventisette et al., 2011)

Full size image

Displacement and velocity time series of Santa Trada landslide obtained from the GB-InSAR system (after Del Ventisette et al., 2011)

Full size image

• Differential: the time span between the first and last image composing the interferogram is kept constant (e.g. 1 day) and the interferograms represent sequential moments (e.g. day-by-day displacement). This approach is used in particular for kinematic analyses as it permits to identify acceleration phases because it enables to evaluate two comparable time periods.

• Incremental: the first image is taken as a reference and the interferograms are all calculated between the reference and the last image; in this way the time span is increasing with time. This approach is useful to evaluate the total cumulative displacement and to measure the displacement even in the slower portions of the landslide. Furthermore, it is suitable for spatial analyses as the total extension of the unstable area can be assessed. On the other hand, long time intervals can cause phase ambiguity and loss of coherence.

For the Santa Trada landslide, thanks to the displacement maps provided by the GB-InSAR system, it was possible to delimitate the area affected by the movement (Fig. 11) and to identify some temporal phases characterized by different activity levels and to assess the risk scenarios temporal evolution (Fig. 12).

The reliability and effectiveness of the described remote sensing techniques, as well as their synergic use, have been enhanced, providing a wide range of surveying and monitoring activities for different landslide types (Fig. 13).

Schematised workflow of the applied techniques showing the different degree of connection between the advanced products (coming from the available instrumentation) which are the basis for addressing the various landslides issues

Full size image

A brief overview of applications (individual or combined technique) has been shown through some selected case studies in section 3. In this section the abovementioned case studies are discussed in order to show, for the employed techniques their main advantages and limitations and the possibility of a synergic use for different applications (Table 4).

Spaceborne platforms

With the increasing large constellation of VHR satellites, imagery can be acquired timely after major landslide events and with daily temporal resolution at nearly global coverage. The main advantage of VHR imagery is the great density of spatial information, whereas, with more competing satellite operators entering the market, prices constantly decrease. A main advantage of optical datasets is their synergetic values for several other applications such as:

• post-disaster damage assessment;

• updating of land cover and landslide inventory maps;

• corresponding archives are often available over a given area.

A greater diversity of platforms increases the chance to acquire cloud free imagery of a given area with a specified time frame; nevertheless atmospheric conditions remain an important factor that, depending on the climate zone and the season, may delay the acquisition of suitable images. For these reasons, satellite tasked for images with sub-meter resolution can still be associated with considerable costs. Higher spatial, spectral (e.g. WorldView-2) and temporal resolutions strongly increase the computational load for the storage and analysis of the datasets, especially for mapping over wide areas. This can considerably slow down the analysis and may need for further investments in hardware and software. This is closely related to the desirable exploitation of spatial context which is typically computational intensive. Pixel-based change detection (typically image differencing) is relatively easy to apply and can be accurate when most of surface changes are caused by landslides. In many cases, it might also be possible to account for sensor and illumination differences by cross calibration and image transformation. However, only limited accuracy can be expected from such approaches in situation where other similar surface changes such as deforestation or barren fields are present in the same scene. The selection of an appropriate threshold to distinguish between changed and unchanged areas remains as a general difficulty for the application of pixel-based methods. Further problems are usually encountered when pixel-base change detection is applied on VHR imagery because of the higher spectral variance and stronger impacts of small co-registration errors. Due to a better exploitation of the spatial context within remote sensing images, OOA approaches generally yield better results than could be achieved with per pixel analyses. It has been demonstrated that OOA rule sets are not only capable to accurately delineate areas affected by landslides but can also be used to distinguish among different landslide types (Martha et al. 2010). In the study area of Giampilieri (Lu et al., 2011) the OOA approach has proved an effective tool to map rapid landslides, and support the local authorities and civil protection department for the emergency management. For both the obtained number and spatial extent of detected landslides, the results show a lower producer’s accuracy than user’s accuracy: specifically, ca. 31% of all manually mapped landslides were omitted in the OOA based detection. This indicates an overestimation of false positives during their classification, accompanied with an underestimation of true positives obtained from the membership function of the selected samples. Further improvements should include a more accurate definition of these thresholds for classifying false positives and a more careful selection of representative samples.

As discussed by Colesanti & Wasowski (2006), due to the inherent limitations of current space observation systems and relevant data processing techniques, the practical applicability of Persistent Scatterer Interferometry (PSI) approaches is usually limited to two landslide classes of the Cruden & Varnes (1996) classification: extremely slow and very slow movements (vel < 16 mm/year and 16 mm/year ≤ vel < 1.6 m/year, respectively). These phenomena are suitable for analysis based on PSI techniques, as long as they evolve with very low displacement rates (few tens of centimeters per year) and their velocities do not exceed the intrinsic limits of the techniques (related to the radar wavelength, revisiting time of the platform and the spatial density of measurement points). Moreover, landslide-induced displacements, detectable through PSI techniques, are restricted to “coherent” landslides with very slow dynamics, i.e. with little internal deformation such as deep-seated deformations (García-Davalillo et al, 2014), creep (Cascini et al., 2010), and, in some cases, slides (Raspini et al., 2015b) roto-translational slides (Tofani et al. 2013a, b), rockslide (Lauknes et al., 2010), complex landslides (Bardi et al., 2014), slow earth flows (Herrera et al., 2011) and badlands (Herrera et al., 2009). InSAR data can provide useful information about pre-event movements, often characterized by low displacement rates (few cm/year) persisting over long time periods (Bardi et al., 2016; Frodella et al., 2016). This deformation regime is quite different to failure events, which occur suddenly and may produce ground displacements of several meters (Raspini et al., 2015a).

With reference to landslide detection/mapping the most advantageous aspects of the multipass D-InSAR approach are (Colesanti and Wasowski, 2006):

• The cost-effectiveness for wide-area (hundreds and thousands of km²) applications, typical of spaceborne remotely sensed data.

• The high density of benchmarks (up to several hundreds per km²).

• The use of “natural” benchmarks not requiring deployment and maintenance.

• The possibility of geo-locating the benchmarks with a precision in the order of 1–5 meters.

• The availability of the extremely valuable ESA (European Space Agency) ERS archive spanning about 20 years, which enables to carry out retrospective studies.

Recent studies proved the feasibility of combining (stitching) SAR data acquired by different sensors (e.g., ENVISAT with ERS, or RADARSAT-1 with RADARSAT-2), despite slight differences in critical image acquisition parameters. Regular revisiting time in the order of 20–40 days (up to 6 days with the new Sentinel-1 ESA mission).

On the whole, the case studies described in the scientific literature highlight that with reference to the detection/mapping of slow-moving landslide phenomena the main benefits regard:

• the definition of the boundaries of already detected mass movements;

• the definition of the states of activity;

• the detection of previously unmapped unstable areas.

However, several limiting factors need to be properly taken into account (Colesanti and Wasowski, 2006):

• Displacement data represent the one-dimensional projection along the LOS projection of a deformation that can actually occur in all three dimensions.

• The ambiguity of phase measurements implies the impossibility to track correctly (i.e., unambiguously) the relative LOS displacement between two scatterers exceeding λ/4 (=1.4 cm for ERS) within one revisiting time interval (35 days for ERS), i.e. approximately 14.5 cm/yr. In practice, it is extremely difficult to detect LOS displacement rates exceeding 8 – 10 cm/yr in the presence of low density of stable scatterers, such as in the case of landslides where topography and vegetation introduce a limitation in the number of detected scatterers. This limits the use of multi-interferometric approaches only to landslides ranging from extremely to very slow phenomena according to the velocity classification of Cruden and Varnes (1996).

• Limited versatility in terms of (a) positioning of the measurement points and (b) revisiting time. Both factors (a) and (b) cannot be optimized as degrees of freedom while planning an analysis.

• Finally, it is still difficult to forecast the coherent pixel density in rural areas without carrying out at least several processing steps on a significant number (15–20) of SAR images.

UAV and Ground-based methods

In terrestrial applications it is still necessary to perform time-consuming and hazardous activities, such as placing ground control points on the analysed scenario (Stavroulaki et al., 2016), if not supported by other technologies (Forlani et al., 2014). In this framework UAV photogrammetry has the following advantages: real-time, flexibility, high-resolution, low costs, as it allows the collection of information in dangerous environments without risk (Chang-chun et al. 2011). Furthermore, the recent development of new algorithms for digital photogrammetry, based on Structure from Motion (SfM) (Westoby et al. 2012) and Multi-View Stereo (MVS) (James and Robson, 2012) techniques, allows obtaining high-resolution 3D models, even by using compact and consumer-grade digital cameras (Lucieer et al., 2013; Rossi et al., 2016). In the case of landslide monitoring and characterization, acquiring aerial imagery using drones permits to overcome some limits of ground-based photogrammetric surveying, such as shadowing effects, which can drastically reduce the accuracy of the resulting digital models.

The synergic use of TLS and IRT has been applied at Elba island case study (section 3.2.1). The TLS survey yielded a detailed 3-D remote structural, geometrical, and geomechanical characterization of the investigated rock masses. In particular, a semiautomatic geomechanical survey made possible the automatic calculation of six of the ten parameters suggested by ISRM (1985) for the quantitative description of discontinuities (orientation, spacing, persistence, roughness, number of sets, and block size). A total of 1359 planes were recognized and clustered according to seven different discontinuity sets, adding two more discontinuity sets to the five detected by means of the traditional field survey, therefore improving the rock mass structural-geomechanical characterization (Fig. 7). The obtained TLS 3-D products also provided reference morphological maps useful for both further detailed field inspections and the design of possible future restoration works. The proposed approach proved to be an effective tool in the field of emergency management, when it is often urgently necessary and gather all the required information (characterization and mapping) as fast as possible in dangerous environments.

The investigated area showed favorable logistic conditions; in fact, the roadway at the foot of the investigated rock slope was fundamental in carrying out up-close the field inspections, the TLS, and the TIR surveys. Had this condition not existed, the point cloud resolution would not probably have been high enough for such detailed analyses. In San Leo an integrated use of GB-InSAR, TLS and IRT has been applied for landslide monitoring (section 3.2.2). The GB-InSAR one year monitoring campaign allowed analysing the short-term behavior of the 2014 rockfall event deposits (Frodella et al., 2016). A general stability of the town monitored structures and analysed rock wall was assessed, and 4 critical sectors were detected in the monitored area by means of 2D displacement maps, while 7 acceleration events were detected from the GB-InSAR control points time series. The accuracy in locating a GBInSAR control point is controlled by the system azimuth and range resolutions, which are in terms related to the distance between the sensor and the backscattering objects. Regarding the San Leo case study, the investigated rock wall has a subvertical geometry which in some portions shows overhanging sectors; in this framework, the GB-InSAR system installation frontal with respect to the surveyed scenario, lead to a different range resolution of scenario sectors located at different heights along the surveyed rock wall (Fig. 10). The obtained 3D displacement map in fact shows an ambiguity in locating deformation sectors along the rock wall height; therefore, in this specific case study sectors characterized by displacements are displayed as vertical zones instead of pixel clusters located at the cliff bottom (Fig. 9). The 3D GB-InSAR displacement map was also used for a comparison between TLS and GB-InSAR displacement data. The difference in the recorded displacements detected by the two monitoring systems in correspondence of the detected rock block toppling is related to the different displacement components of the recorded movements, due to the two different monitoring systems LOSs. The analyzed rock cliff is characterized by complex geomorphological and geometric features, different ongoing landslide processes with various state of activity. Each single employed monitoring technique can be considered not adequate for their intrinsic limitations. The combined use of the abovementioned techniques provided an effective monitoring system for landslide characterization and state of activity monitoring, thanks to the different instrument characteristics (LOS, range of detectable velocity, repetition time), which allowed to overcome the limitations of each single employed technique. Compared to GB-InSAR, TLS does not suffer from problematics such as loss of coherence, decorrelation, and displacement detection capability only along the sensor LOS; on the other hand, GB-InSAR single measure can reach sub-millimeter accuracy, while using a TLS it is not possible to easily detect displacements smaller than 10 mm. Therefore, TLS was considered more suitable for the detection, characterization and volume assessment of the minor rockfall events affecting the newly formed cliff (too fast phenomena to be detected by means of the GB-InSAR system). In both the abovementioned case studies, the capability to remotely collect the scenario surface temperature at a detailed spatial resolution proved that IRT can be usefully applied in the field landslide analyses. Nevertheless, IRT alone is insufficient for a complete landslide characterization; in order to obtain a more accurate interpretation of the results IRT could be more profitably used as an ancillary low cost technique through the integration with other ground based remote sensing techniques, such as TLS and GB-InSAR. Future developments should include the application of fixed IRT installations for gathering continuous, high-resolution, real-time data to be compared with those of the integrated GB-InSAR and TLS monitoring systems. Furthermore, the portability of modern thermal cameras and the rapid evolution of IRT technology opens up future scenarios of automated inspections, which could be quickened by using remotely controlled UAV platforms. The Santa Trada landslide offers a good example of GB-InSAR application in emergency conditions in order to assess the risk impending on a critical infrastructure. In Santa Trada area, a GB-InSAR device was promptly installed in order to understand the temporal evolution of a landslide that seriously threatened the functionality and the safety of a strategic road infrastructure. This technique worked with all weather conditions and with a continuous surveillance for all the time of emergency, allowing the rapid assessment of the overall dynamics of the instable slope and related risks scenarios. This application was among the first to demonstrate the full effectiveness of this system in managing landslides emergencies since it greatly facilitated the intervention operations by designated authorities that aimed to restore a normal service in the shortest possible time after a precautionary closure of the motorway. It also granted the possibility to acquire data during precipitations, which represent the most critical moment in stability terms; indeed, usually traditional monitoring instruments are not able to work in such conditions and to provide such useful information real-time, since they normally require a longer time for the installation and a direct access to the unstable and unsafe sites. The GB-InSAR technique proved to be a suitable and versatile tool to assess the actual hazard of the landslide in order to enable the re-opening of the motorway.

In the recent years, remote sensing techniques for landslide analysis have been interested by several technical and scientific improvements. These techniques can play an important role in landslide risk management, as they allow the representation of large surfaces with dense spatial sampling, offering clear advantages with respect to traditional topographical systems (such as GPSs and robotized total stations), which on the contrary provide data that are accurate but necessarily limited to a small number of control points (Teza et al., 2008). Currently landslide analysis operators can select the most proper methodology with respect to their specific needs, which can be related to the different civil protection phases, technical issues, financial budget, environmental factors and specific features of the studied instability phenomena. The selection of the proper methodology to be adopted can be related to technical issues, economic budget, environmental factors and specific feature of the landslides to be monitored. Because of a growing demand for effective Civil Protection procedures in pre- and post-disaster initiatives in landslide-prone areas, the purpose of researchers in the near future is to improve the investigative capacity of the such instruments and consequently to extend their fields of application. In fact, landslides managing in order to reduce vulnerability is currently considered more feasible (in terms of faster and extensive results) than governing all the natural conditions leading to instability, such as the spatial distribution of geology and geomorphology and the climatic influence.

Authors’ contributions

NC conceived and structured the whole manuscript. He also supervised the writing, especially during the organization of contributions coming from different techniques. WF, SM, EI were responsible for the preparation of the ground-based sections. VT, AC, FR, PL were responsible for the preparation of the spaceborne sections. GR and LT were responsible for the preparation of the UAV sections. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Authors and Affiliations

1. Department of Earth Sciences, University of Florence, Via G. La Pira 4, 50121, Florence, Italy

   Nicola Casagli, William Frodella, Stefano Morelli, Veronica Tofani, Andrea Ciampalini, Emanuele Intrieri, Federico Raspini, Guglielmo Rossi & Luca Tanteri

2. College of Surveying and Geo-Informatics, Tongji University, Shanghai, 200092, China

   Ping Lu

Corresponding author

Correspondence to Stefano Morelli.

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions