Critical rainfall conditions for the initiation of torrential flows: results from the Rebaixader catchment (Central Pyrenees)

Torrential flows like debris flows or debris floods are fast movements formed by a mix of water and different amounts of unsorted solid material. They generally occur in steep torrents and pose high risk in mountainous areas. Rainfall is their most common triggering factor and the analysis of the critical rainfall conditions is a fundamental research task. Due to their wide use in warning systems, rainfall thresholds for the triggering of torrential flows are an important outcome of such analysis and are empirically derived using data from past events. In 2009, a monitoring system was installed in the Rebaixader catchment, Central Pyrenees (Spain). Since then, rainfall data of 25 torrential flows (“TRIG rainfalls”) were recorded, with a 5-min sampling frequency. Other 142 rainfalls that did not trigger torrential flows (“NonTRIG rainfalls”) were also collected and analyzed. The goal of this work was threefold: (i) characterize rainfall episodes in the Rebaixader catchment and compare rainfall data that triggered torrential flows and others that did not; (ii) define and test Intensity–Duration (ID) thresholds using rainfall data measured inside the catchment by with different techniques; (iii) analyze how the criterion used for defining the rainfall duration and the spatial variability of rainfall influences the value obtained for the thresholds. The statistical analysis of the rainfall characteristics showed that the parameters that discriminate better the TRIG and NonTRIG rainfalls are the rainfall intensities, the mean rainfall and the total rainfall amount. The antecedent rainfall was not significantly different between TRIG and NonTRIG rainfalls, as it can be expected when the source material is very pervious (a sandy glacial soil in the study site). Thresholds were derived from data collected at one rain gauge located inside the catchment. Two different methods were applied to calculate the duration and intensity of rainfall: (i) using total duration, Dtot, and mean intensity, Imean, of the rainfall event, and (ii) using floating durations, D, and intensities, Ifl, based on the maximum values over floating periods of different duration. The resulting thresholds are considerably different (Imean = 6.20 Dtot-0.36 and Ifl_90% = 5.49 D-0.75, respectively) showing a strong dependence on the applied methodology. On the other hand, the definition of the thresholds is affected by several types of uncertainties. Data from both rain gauges and weather radar were used to analyze the uncertainty associated with the spatial variability of the triggering rainfalls. The analysis indicates that the precipitation recorded by the nearby rain gauges can introduce major uncertainties, especially for convective summer storms. Thus, incorporating radar rainfall can significantly improve the accuracy of the measured triggering rainfall. Finally, thresholds were also derived according to three different criteria for the definition of the duration of the triggering rainfall: (i) the duration until the peak intensity, (ii) the duration until the end of the rainfall; and, (iii) the duration until the trigger of the torrential flow. An important contribution of this work is the assessment of the threshold relationships obtained using the third definition of duration. Moreover, important differences are observed in the obtained thresholds, showing that ID relationships are significantly dependent on the applied methodology.Peer ReviewedPostprint (author's final draft

Digital image correlation (DIC) analysis of the 3 December 2013 Montescaglioso landslide (Basilicata, Southern Italy). Results from a multi-dataset investigation

Image correlation remote sensing monitoring techniques are becoming key tools for providing effective qualitative and quantitative information suitable for natural hazard assessments, specifically for landslide investigation and monitoring. In recent years, these techniques have been successfully integrated and shown to be complementary and competitive with more standard remote sensing techniques, such as satellite or terrestrial Synthetic Aperture Radar interferometry. The objective of this article is to apply the proposed in-depth calibration and validation analysis, referred to as the Digital Image Correlation technique, to measure landslide displacement. The availability of a multi-dataset for the 3 December 2013 Montescaglioso landslide, characterized by different types of imagery, such as LANDSAT 8 OLI (Operational Land Imager) and TIRS (Thermal Infrared Sensor), high-resolution airborne optical orthophotos, Digital Terrain Models and COSMO-SkyMed Synthetic Aperture Radar, allows for the retrieval of the actual landslide displacement field at values ranging from a few meters (2–3 m in the north-eastern sector of the landslide) to 20–21 m (local peaks on the central body of the landslide). Furthermore, comprehensive sensitivity analyses and statistics-based processing approaches are used to identify the role of the background noise that affects the whole dataset. This noise has a directly proportional relationship to the different geometric and temporal resolutions of the processed imagery. Moreover, the accuracy of the environmental-instrumental background noise evaluation allowed the actual displacement measurements to be correctly calibrated and validated, thereby leading to a better definition of the threshold values of the maximum Digital Image Correlation sub-pixel accuracy and reliability (ranging from 1/10 to 8/10 pixel) for each processed dataset

Post-failure evolution analysis of a rainfall-triggered landslide by multi-temporal interferometry SAR approaches integrated with geotechnical analysis

Persistent Scatterers Interferometry (PSI) represents one of the most powerful techniques for Earth's surface deformation processes' monitoring, especially for long-term evolution phenomena. In this work, a dataset of 34 TerraSAR-X StripMap images (October 2013–October 2014) has been processed by two PSI techniques - Coherent Pixel Technique-Temporal Sublook Coherence (CPT-TSC) and Small Baseline Subset (SBAS) - in order to study the evolution of a slow-moving landslide which occurred on February 23, 2012 in the Papanice hamlet (Crotone municipality, southern Italy) and induced by a significant rainfall event (185 mm in three days). The mass movement caused structural damage (buildings' collapse), and destruction of utility lines (gas, water and electricity) and roads. The results showed analogous displacement rates (30–40 mm/yr along the Line of Sight – LOS-of the satellite) with respect to the pre-failure phase (2008–2010) analyzed in previous works. Both approaches allowed detect the landslide-affected area, however the higher density of targets identified by means of CPT-TSC enabled to analyze in detail the slope behavior in order to design possible mitigation interventions. For this aim, a slope stability analysis has been carried out, considering the comparison between groundwater oscillations and time-series of displacement. Hence, the crucial role of the interaction between rainfall and groundwater level has been inferred for the landslide triggering. In conclusion, we showed that the integration of geotechnical and remote sensing approaches can be seen as the best practice to support stakeholders to design remedial works.Peer ReviewedPostprint (author's final draft

Basin scale assessment of landslides geomorphological setting by advanced InSAR analysis

An extensive investigation of more than 90 landslides affecting a small river basin in Central Italy was performed by combining field surveys and remote sensing techniques. We thus defined the geomorphological setting of slope instability processes. Basic information, such as landslides mapping and landslides type definition, have been acquired thanks to geomorphological field investigations and multi-temporal aerial photos interpretation, while satellite SAR archive data (acquired by ERS and Envisat from 1992 to 2010) have been analyzed by means of A-DInSAR (Advanced Differential Interferometric Synthetic Aperture Radar) techniques to evaluate landslides past displacements patterns. Multi-temporal assessment of landslides state of activity has been performed basing on geomorphological evidence criteria and past ground displacement measurements obtained by A-DInSAR. This step has been performed by means of an activity matrix derived from information achieved thanks to double orbital geometry. Thanks to this approach we also achieved more detailed knowledge about the landslides kinematics in time and space

Automated Satellite-Based Landslide Identification Product for Nepal

Landslide event inventories are a vital resource for landslide susceptibility and forecasting applications. However, landslide inventories can vary in accuracy, availability, and timeliness as a result of varying detection methods, reporting, and data availability. This study presents an approach to use publicly available satellite data and open source software to automate a landslide detection process called the Sudden Landslide Identification Product (SLIP). SLIP utilizes optical data from the Landsat 8 OLI sensor, elevation data from the Shuttle Radar Topography Mission (SRTM), and precipitation data from the Global Precipitation Measurement (GPM) mission to create a reproducible and spatially customizable landslide identification product. The SLIP software applies change detection algorithms to identify areas of new bare-earth exposures that may be landslide events. The study also presents a precipitation monitoring tool that runs alongside SLIP called the Detecting Real-time Increased Precipitation (DRIP) model that helps identify the timing of potential landslide events detected by SLIP. Using SLIP and DRIP together, landslide detection is improved by reducing problems related to accuracy, availability, and timeliness that are prevalent in the state-of-the-art of landslide detection. A case study and validation exercise was performed in Nepal for images acquired between 2014 and 2015. Preliminary validation results suggest 56% model accuracy, with errors of commission often resulting from newly cleared agricultural areas. These results suggest that SLIP is an important first attempt in an automated framework that can be used for medium resolution regional landslide detection, although it requires refinement before being fully realized as an operational tool

Bayesian Analysis of the Impact of Rainfall Data Product on Simulated Slope Failure for North Carolina Locations

In the past decades, many different approaches have been developed in the literature to quantify the load-carrying capacity and geotechnical stability (or the factor of safety, Fs) of variably saturated hillslopes. Much of this work has focused on a deterministic characterization of hillslope stability. Yet, simulated Fs values are subject to considerable uncertainty due to our inability to characterize accurately the soil mantles properties (hydraulic, geotechnical, and geomorphologic) and spatiotemporal variability of the moisture content of the hillslope interior. This is particularly true at larger spatial scales. Thus, uncertainty-incorporating analyses of physically based models of rain-induced landslides are rare in the literature. Such landslide modeling is typically conducted at the hillslope scale using gauge-based rainfall forcing data with rather poor spatiotemporal coverage. For regional landslide modeling, the specific advantages and/or disadvantages of gauge-only, radar-merged and satellite-based rainfall products are not clearly established. Here, we compare and evaluate the performance of the Transient Rainfall Infiltration and Grid-based Regional Slope-stability analysis (TRIGRS) model for three different rainfall products using 112 observed landslides in the period between 2004 and 2011 from the North Carolina Geological Survey database. Our study includes the Tropical Rainfall Measuring Mission (TRMM) Multi-satellite Precipitation Analysis Version 7 (TMPA V7), the North American Land Data Assimilation System Phase 2 (NLDAS-2) analysis, and the reference truth Stage IV precipitation. TRIGRS model performance was rather inferior with the use of literature values of the geotechnical parameters and soil hydraulic properties from ROSETTA using soil textural and bulk density data from SSURGO (Soil Survey Geographic database). The performance of TRIGRS improved considerably after Bayesian estimation of the parameters with the DiffeRential Evolution Adaptive Metropolis (DREAM) algorithm using Stage IV precipitation data. Hereto, we use a likelihood function that combines binary slope failure information from landslide event and null periods using multivariate frequency distribution-based metrics such as the false discovery and false omission rates. Our results demonstrate that the Stage IV-inferred TRIGRS parameter distributions generalize well to TMPA and NLDAS-2 precipitation data, particularly at sites with considerably larger TMPA and NLDAS-2 rainfall amounts during landslide events than null periods. TRIGRS model performance is then rather similar for all three rainfall products. At higher elevations, however, the TMPA and NLDAS-2 precipitation volumes are insufficient and their performance with the Stage IV-derived parameter distributions indicates their inability to accurately characterize hillslope stability

ANALYZING THE LIFE-CYCLE OF UNSTABLE SLOPES USING APPLIED REMOTE SENSING WITHIN AN ASSET MANAGEMENT FRAMEWORK

An asset management framework provides a methodology for monitoring and maintaining assets, which include anthropogenic infrastructure (e.g., dams, embankments, and retaining structures) and natural geological features (e.g., soil and rock slopes). It is imperative that these assets operate efficiently, effectively, safely, and at a high standard since many assets are located along transportation corridors (highways, railways, and waterways) and can cause severe damage if compromised. Assets built on or around regions prone to natural hazards are at an increased risk of deterioration and failure. The objective of this study is to utilize remote sensing techniques such as InSAR, LiDAR, and optical photogrammetry to identify assets, assess past and current conditions, and perform long-term monitoring in transportation corridors and urbanized areas prone to natural hazards. Provided are examples of remote sensing techniques successfully applied to various asset management procedures: the characterization of rock slopes (Chapter 2), identification of potentially hazardous slopes along a railroad corridor (Chapter 3), monitoring subsidence rates of buildings in San Pedro, California (Chapter 4), and mapping displacement rates on dams in India (Chapter 5) and California (Chapter 6). A demonstration of how InSAR can be used to map slow landslides (those with a displacement rate \u3c 16 mm/year and may be undetectable without sensitive instrumentation) and update the California Landslide Inventory on the Palos Verdes Peninsula is provided in Chapter 7. Long-term landslide monitoring using optical photogrammetry, GPS, and InSAR measurements is also used to map landslide activity at three orders of magnitude (meter to millimeter scales) in Chapter 8. Remote sensing has proven to be an effective tool at measuring ground deformation, which is an implicit indicator of how geotechnical asset condition changes (e.g., deteriorates) over time. Incorporating these techniques into a geotechnical asset management framework will provide greater spatial and temporal data for preventative approaches towards natural hazards

Multi-technique approach to rockfall monitoring in the Montserrat massif (Catalonia, NE Spain)

Montserrat Mountain is located near Barcelona in Catalonia, in the northeast of Spain, and its massif is formed by conglomerate interleaved by siltstone/sandstone with steep slopes very prone to rockfalls. The increasing number of visitors in the monastery area, reaching 2.4 million per year, has highlighted the risk derived from rockfalls for this building area and also for the terrestrial accesses, both roads and the rack railway. A risk mitigation plan has been launched, and its first phase during 2014-2016 has been focused largely on testing several monitoring techniques for their later implementation. The results of the pilot tests, performed as a development from previous sparse experiences and data, are presented together with the first insights obtained. These tests combine four monitoring techniques under different conditions of continuity in space and time domains, which are: displacement monitoring with Ground-based Synthetic Aperture Radar and characterization at slope scale, with an extremely non-uniform atmospheric phase screen due to the stepped topography and atmosphere stratification; Terrestrial Laser Scanner surveys quantifying the frequency of small or even previously unnoticed rockfalls, and monitoring rock block centimetre scale displacements; the monitoring of rock joints implemented through a wireless sensor network with an ad hoc design of ZigBee loggers developed by ICGC; and, finally, monitoring singular rock needles with Total Station.Peer ReviewedPostprint (author's final draft