Next Article in Journal
Uptake and Recovery of Gold from Simulated Hydrometallurgical Liquors by Adsorption on Pine Bark Tannin Resin
Next Article in Special Issue
Rainfall-Induced Shallow Landslide Detachment, Transit and Runout Susceptibility Mapping by Integrating Machine Learning Techniques and GIS-Based Approaches
Previous Article in Journal
Coupling Sediment Transport Dynamics with Sediment and Discharge Sources in a Glacial Andean Basin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Using Field-Based Monitoring to Enhance the Performance of Rainfall Thresholds for Landslide Warning

by
Minu Treesa Abraham
1,
Neelima Satyam
1,
Maria Alexandra Bulzinetti
2,
Biswajeet Pradhan
3,4,5,
Binh Thai Pham
6 and
Samuele Segoni
2,*
1
Discipline of Civil Engineering, Indian Institute of Technology Indore, Madhya Pradesh 452020, India
2
Department of Earth Sciences, University of Florence, 50121 Florence, Italy
3
Centre for Advanced Modelling and Geospatial Information Systems (CAMGIS), Faculty of Engineering and Information Technology, University of Technology Sydney, Sydney P.O. Box 123, Australia
4
Center of Excellence for Climate Change Research, King Abdulaziz University, P.O. Box 80234, Jeddah 21589, Saudi Arabia
5
Earth Observation Center, Institute of Climate Change, Universiti Kebangsaan Malaysia, Bangi 43600, Malaysia
6
Department of Geotechnical Engineering, University of Transport Technology, Hanoi 100000, Vietnam
*
Author to whom correspondence should be addressed.
Water 2020, 12(12), 3453; https://doi.org/10.3390/w12123453
Submission received: 3 November 2020 / Revised: 5 December 2020 / Accepted: 7 December 2020 / Published: 9 December 2020
(This article belongs to the Special Issue Rainfall-Induced Shallow Landslides Modeling and Warning)

Abstract

:
Landslides are natural disasters which can create major setbacks to the socioeconomic of a region. Destructive landslides may happen in a quick time, resulting in severe loss of lives and properties. Landslide Early Warning Systems (LEWS) can reduce the risk associated with landslides by providing enough time for the authorities and the public to take necessary decisions and actions. LEWS are usually based on statistical rainfall thresholds, but this approach is often associated to high false alarms rates. This manuscript discusses the development of an integrated approach, considering both rainfall thresholds and field monitoring data. The method was implemented in Kalimpong, a town in the Darjeeling Himalayas, India. In this work, a decisional algorithm is proposed using rainfall and real-time field monitoring data as inputs. The tilting angles measured using MicroElectroMechanical Systems (MEMS) tilt sensors were used to reduce the false alarms issued by the empirical rainfall thresholds. When critical conditions are exceeded for both components of the systems (rainfall thresholds and tiltmeters), authorities can issue an alert to the public regarding a possible slope failure. This approach was found effective in improving the performance of the conventional rainfall thresholds. We improved the efficiency of the model from 84% (model based solely on rainfall thresholds) to 92% (model with the integration of field monitoring data). This conceptual improvement in the rainfall thresholds enhances the performance of the system significantly and makes it a potential tool that can be used in LEWS for the study area.

1. Introduction

Landslides are geomorphological processes, which can sometimes lead to severe loss of lives and properties [1,2]. Due to the increase in population, most hilly terrains are being urbanised, and the risk associated with landslides is increasing worldwide in similar environments, leading to the need of establishing early warning systems. The definition of rainfall thresholds is the most widely followed approach for the setting up of landslide early warning systems [3,4]. A rainfall threshold can be defined as the rainfall conditions that, when reached or exceeded, are likely to trigger landslides [5]. While at the local scale rainfall thresholds can be defined with physically-based approaches [5,6,7,8], at regional scale empirical approaches are used and the threshold definition is usually based on establishing a statistical correlation between the occurrence of landslides in the study area and some parameters describing the triggering rainfall [9,10,11,12]. Generally speaking, three parameters are considered to define a rainfall event: (i) The cumulative event rainfall (E), which is the total precipitation recorded from the beginning of continuous rainfall till the time of occurrence of failure; (ii) the rainfall duration (D) and (iii) the rainfall intensity (I), which is calculated as the ratio of cumulated event rainfall to duration. ID thresholds are defined on a 2D plane with I on y axis and D on x axis, and similarly ED thresholds with E and D on y and x axes, respectively. Thresholds are calibrated according to historical data and since the pioneer work of Caine [10], they are mathematically represented by power law equations, with the parameter on y axis (I or E) expressed in terms of the one on x axis (D). If the forecasted rainfall is above this threshold line, a slope failure can be expected in the region. The thresholds define a lower bound to the critical rainfall parameters, and no landslides are expected to occur below this minimum threshold line [9,13]. The conventional approach that defines rainfall thresholds is the ID thresholds [9,10,14,15,16], but more recently the functional relationship between I and D has led to the development of ED thresholds [17,18,19]. Moreover, automatic algorithms were recently developed for the objective identification of responsible rainfall events [18,20,21,22]. All these methods assume that the hydro-geological conditions of the study area remain unchanged, with respect to time, and the rainfall conditions that have triggered landslides in the past can trigger landslides in the future as well.
However, the technique of regional scale rainfall thresholds has some limitations, associated mainly with the intrinsic uncertainties of the empirical approach and with the reliability of historical rainfall and landslide data, which usually lead to define conservative thresholds, resulting in a very high number of false alarms [23,24,25,26,27,28]. To overcome this limitation, several authors tried to strengthen the classical rainfall threshold approach by proposing hydro-meteorological thresholds [29], integrating soil moisture measurements [30,31], using weighted indexes [32], considering a third additional parameter [33] or accounting for additional physical processes [34,35]. However, the approach of improving rainfall thresholds effectiveness by an integration with instrumental monitoring remains relatively unexplored in regional scale studies [36]. This is because getting information about slope deformation with good spatial and temporal accuracy on wide areas is usually very challenging, as with remote sensing-based monitoring methods: Ground-based and satellite-based Synthetic Aperture Radar (SAR), Light Detection and Ranging (LiDAR), Permanent Scatterer Synthetic Aperture Radar Interferometry (PSInSAR) and many others [37,38,39,40,41,42,43,44,45].
The advances in Internet of Things (IoT) and wireless networks, data transmission and real-time monitoring have become cost-effective approaches [46]. Hence, instrumental monitoring using low-cost sensors [47] is now becoming a popular approach for monitoring unstable slopes. Measuring the displacements of moving soil mass using extensometers is a physical-based monitoring technique which has been in practice for a long time [43]. Another approach is the use of deep hole inclinometers, which is effective in detecting the deep slip surfaces [48]. Such sensors can be connected via wireless networks for real-time data transmission, opening a wide range of possibilities in monitoring of unstable areas [49,50,51].
In this study, in a test site in India, a rainfall threshold approach was coupled with instrumental monitoring data provided by MicroElectroMechanical Systems (MEMS) tilt sensors, which have already proved to be a reliable and low-cost method [47,52]. MEMS devices have become an essential component in a wide range of applications, ranging from medical, engineering, mechanical, biological, hydraulic sciences, proving to be useful allies also in landslide prediction [53,54]. The sensors used in this study were originally developed in Japan and since then have been found effective in China and Japan [55,56]. The climatic and geological settings of the study area in India differ from the locations where MEMS tilt sensors have been already tested. In a global database of fatal landslides, it has been mentioned that Asia contributes 75% of the world’s fatal landslides, with a significant number of events in the Himalayan Arc [57]. The study also states that landslide events are increasing in in India, Nepal, Bangladesh, Northern Pakistan, and Bhutan during the summer monsoon. When rainfall induced landslides are considered, India and Nepal contribute 16% and 10% of the events in the global dataset [57]. Several studies are conducted for deriving the empirical rainfall thresholds for the occurrence of landslides in this region [16,58,59,60]. The Himalayan Arc [57,61] and Western Ghats [62] are critical landslide prone areas in India, contributing a major share of the total landslides in the country. Hence, it is important to verify its applicability in the study area, by comparing the readings from sensors with field observations. The primary triggering factor of landslides is rainfall, and hence several rainfall thresholds are being used to forecast the occurrence of landslides in this region [22,63,64,65,66,67,68]. This study evaluates the use of field-based monitoring data, for enhancing the forecasting performance of empirical rainfall thresholds defined for the study area in Darjeeling Himalayas.
An algorithm is proposed and validated in this study, to integrate the field monitoring data with the rainfall thresholds. The algorithm can be used by the authorities to issue alerts to the public for possible occurrence of landslides. This study, making use of both historical data and real time observations for landslide forecasting, is the first of its kind for Darjeeling Himalayas. The proposed approach signifies the need for a network of real-time field monitoring sensors on critical slopes. For a country like India, the investment in such low-cost sensor networks can help in reducing the risk associated with landslides.

2. Study Area

Kalimpong town is situated in the state of West Bengal, India (Figure 1). The town is close to Sikkim, West Bengal border in the Darjeeling Himalayas. Kalimpong town lies in between two rivers, Tista and Relli. The western slopes towards Tista River are very steep, and the slopes towards Relli River are gentle. The ground inclination varies from 20° near Tista River and increases to around 40o at higher elevations towards the town. The morphometric classification ranges from gentle slopes to escarpment. The topsoil consists of soil of size ranging from silt to sand. Colluvium, older debris, and young debris deposits are found in the region, which gets eroded during heavy rains.
Kalimpong is a part of Darjeeling Himalayas, with intra-thrusted Fold-Thrust Belt (FTB) rock formation aged from Precambrian to Quaternary [69]. The bedrock lithology consists mainly of golden to silver coloured quartz mica schist of Daling series [70]. The quartz mica schist at some portions grades to dull white mica quartz schist and bright silver white mica schist. The bedrock is largely covered with thick layers of weathered material. The rocks are highly jointed and weathered, resulting in intense land sliding and erosion rates during the monsoon season and the formation of loose soil deposits, especially in the upper portion of the young folded mountain sectors. The general landslide susceptibility of the area is also increased by land use practices (urbanisation and slope cuts for roads and buildings) and lateral erosion of streams. The most common landslide types in the region are debris flows and slides, often influenced by slope cuts and riverbank erosion. During monsoon season, the soil gets saturated due to rainwater. The joints in the rocks also allow the rainwater to percolate. The saturation of soil further reduces the shear strength and leads to slope failure. The region, with high drainage density, has the most critical locations along the mountain rivulets (jhoras in local dialect). The jhoras flowing in high velocity erodes soil from its banks and eventually result in slope failure. One of the most hazardous areas of the test site is the Pyarieni jhora, which is located in Chibo village of Kalimpong (Figure 1c): It develops through deposits of an ancient debris slide consisting of a sandy silty matrix with phyllite fragments, phyllitic quartzite and mica schist [69], with a regolith thickness varying from 4 to 5 m along the banks. The upper portions of the location have scanty exposures of schist, phyllite and sheared gneiss, while the exposures of silvery grey mica schist are found along the jhora banks. The rocks are traversed by two to three sets of joints and are foliated. Seepage of rainwater through the loose topsoil and the jointed rocks makes the slopes unstable during monsoon season.
The rainfall data used for this study is of daily resolution, measured at Tirpai, Kalimpong [71]. The average monthly rainfall data for the monitoring period, from July 2017 to August 2020, is reported in Table 1 below and it can be observed that the major share of the total annual rainfall starts in the month of July, during the monsoon season.
The daily rainfall data during the study period has been plotted in Figure 2. It can be observed that rainfall with an intensity greater than 100 mm per day has happened multiple times during the monsoon. The average annual rainfall of the region is around 2000 mm and during the investigated period the least rainfall was received in 2018, and the maximum amount was received in 2020. From January to August 2020, a total rainfall of 2246.3 mm is recorded in the region.

3. Methodology

The study evaluates the performance of empirical rainfall thresholds, field monitoring data and the combined use of both, in forecasting the slope failures near a mountain rivulet, Pyarieni jhora in Chibo village of Kalimpong. The unstable slopes near the rivulet are being monitored using three tilt sensors since July 2017. The empirical thresholds for the region have been derived using the rainfall and landslide data from 2010 to 2016 and has been validated in this study using the data from July 2017 to August 2020. The field monitoring data and the combined approach are calibrated and verified using the data from July 2017 to August 2020.
Both empirical thresholds and field monitoring are combined using an algorithm, which first considers the rainfall thresholds and verify the in situ conditions using the data from tilt meters. All the three methods are explained in the following sections in detail.

3.1. Empirical Thresholds

Defining the triggering rainfall of landslide events and identifying the associated rainfall parameters are the most critical part in deriving a rainfall threshold [21]. The parameters associated with the responsible rainfall event are then plotted on a two-dimensional plane in logarithmic scale. For example, while defining ID thresholds, the values of I and D of the rainfall events that have resulted in landslides are plotted in y and x axes, respectively. The logarithmic scale is used to plot data of multiple orders. The power law in the original scale will get converted to a straight line in logarithmic scale and the best fit line can be derived using the method of least squares. The definition of empirical thresholds considers only those events which have resulted in landslides. The effect of rainfall events that have not triggered landslides are usually ignored in the analysis.
The thresholds in this study are defined using a frequentist approach, according to the distribution of data around the best fit line. The data plotted on the logarithmic scale are of multiple orders, and the best fit line can be calculated using the method of least squares. On a logarithmic scale, the ID thresholds are generally straight lines with negative slope, and ED thresholds are straight lines with positive slopes. The distance of each data point in the vertical axis from the best fit line is the shift of each data point. The best fit line represents a threshold line with 50% exceedance probability. It is assumed that the data follows normal distribution, and the best fit line represents the mean value. The threshold lines can be defined using the shift of data points from the best fit line. Hence the shift is fitted using a standard Gaussian distribution and thresholds can be defined based on the distribution of data.
The standard form of ID thresholds is given by:
I = α D β
where α and β are empirical parameters in which α   is the intercept and β is the scaling parameter. In the case of ED thresholds, the power law is in of the similar form, given by:
E = a D b
where a and b are empirical parameters, similar to the case of ID thresholds. For the same data, β and b are related by the following equation:
b   =   β + 1
The rainfall thresholds were derived for the region using rainfall and landslide from 2010 to 2016 [22,58] in previous studies conducted by the same research group. The landslide data has been collected from various government sources and media reports to derive a relationship between the occurrence of rainfall and landslides in the region. For the definition of ID thresholds, responsible rainfall events were identified by expert judgement based on the time and location of landslide. Then, the intensity and duration of the rainfall has been calculated from the daily rainfall data, which were used for the definition of minimum threshold, in the form of a power law. The ID threshold (Figure 3a) with an exceedance probability of 5% is derived for Kalimpong [58] as:
I = 3.52   D 0.41 .
To overcome the limitations of the defined ID threshold due to uncertainty in data and the spatial distribution of location of landslides and rain gauge, an automatic algorithm based approach [18] was used to define the ED thresholds (Figure 3b) with 5% exceedance probability for Kalimpong [22].
E = 4.2   D 0.56 .
The algorithm defines the threshold using frequentist approach, but the rainfall event for the initiation of landslide is identified by the reconstruction of all rainfall events based on the rainfall series data. Next, the Maximum Probability Rainfall Conditions (MPRC) were identified using the rainfall condition with highest weight. The weight assigned is directly proportional to the cumulative event rainfall and intensity of the rainfall conditions and inversely proportional to the square of distance between rainfall and landslide. The thresholds were calibrated using the rainfall and landslide data from 2010 to 2016 and have been validated in this study using the data from July 2017 to August 2020.

3.2. Field Monitoring System

For monitoring the unstable slopes near Pyarieni jhora in Chibo, three sensor units were installed in July 2017. Each field monitoring unit consists of an MEMS tilt sensor, a volumetric water content sensor and a control unit as shown in Figure 4.
The MEMS tilt sensor measures the tilt angle in two directions: One parallel to the sloping surface, and one in the perpendicular direction. The abscissa of the sensors was kept parallel to the direction of slope and ordinate in the perpendicular direction. The angles were recorded with an accuracy of 0.017°. The sensitivity of tilt sensors was 4 V/g, and the output was measured as digital voltage. This voltage value was converted to the tilt angle using the conversion equation provided by the manufacturer. The volumetric water content sensors measured the water content with a precision of ±3%, with a response time of 10 ms. The control unit consists of a wireless communication module which transfers the data from each sensor unit to a data logger. The data received from all sensors, were collected by the data logger and transferred in real-time through internet. The power to the sensors was provided by four C size alkaline batteries. Such batteries can work continuously in field conditions, for a long time. The sensor sleeps for 10 min after sending a signal, which increases the battery life in the field. The configuration of the sensor is shown in Figure 5. The sensor for this research has been procured from Chuo Kaihatsu Corporation, Japan [72]. The manufacturing cost of sensors are less when compared to the conventional approaches, due to the integration of MEMS with integrated circuits, with very small dimensions [73].
The sensors can be installed easily at shallow depths. The tilt sensors are placed at a depth of 1 m or slightly more, and volumetric moisture content sensors at around 3 cm depth from the ground surface. Thus, tilt sensors overcome the disadvantages of inclinometers and extensometers, with the easy installation and quick response time [47], making them an easy-to-use tool for providing early warning. The low cost of the sensor unit helps in installing many sensors in a single slope. The major disadvantage associated with the sensors is the possibility of false alarms due to human on animal interventions. Any impact on the peg or pole attached to the sensor can result in tilting of the sensor. Hence, utmost care has to be taken while interpreting the data from sensors. The data from sensors should always be correlated with the rainfall and moisture data before arriving at conclusions.
The locations of the sensors were selected after detailed field investigations, consultation with Geological Survey of India, and discussions with local people [69]. The sensors S1, S2, and S3 were placed near Pyarieni jhora. The location of sensors, along with the drainage map, is shown in Figure 6.
The sensors are used to monitor the slopes since July 2017 and the readings are matched with field observations, to assess the reliability of the data [46,69]. The recordings from the tilt meters during monitoring period is plotted in Figure 7 below.
It can be observed from Figure 7 that the maximum variation in tilt has been observed in sensor 2, which is the most critical location. Also, there are many abrupt changes, observed in the monitoring data, which are mostly vertical, due to sudden impacts. From field observations, it has been noted that such changes are not always associated with slope failures. The slope instabilities are occurring during monsoon seasons, due to heavy rainfall, during which the tilt meters record variations in tilting angles in accordance with the directions of failing mass.
During the 2017 monsoon, the tilting angles of sensors 3 and 2 showed a notable increase in the tilt angles. The displacement periods were 28 and 29 July 2017, and 13 to 17 August 2017. During monsoon 2018, maximum variations in tilting angles were recorded in sensor 2. From 3 August to 15 September 2018, significant changes were observed in tilt angles, with a flat period from 14 to 24 August 2018. During 2019 monsoon, both sensors 1 and 2 showed minor variations. There were two displacement periods, one from 12 to 20 July 2019 and the second was from 9 to 18 August 2019. The displacement periods during 2020 monsoon were from 17 to 22 June and 26 July to 27 August. The second displacement period during 2020 was the most critical (Figure 8) where cracks were observed in the village roads near the jhora (Figure 8). The cracks are close to the location of sensor 2, which has recorded the maximum variations in tilt angles. The recorded data were in accordance with the field observations. The data from all four monsoon seasons were used to calibrate and validate the thresholds for tilt meters.
The peak hourly tilt rate is used for issuing the alert and after evaluating the field monitoring data, it was set to 0.03°/h. All three sensors were considered for analysis, and if a reading crossing the threshold is recorded in any of the three sensors, alert should be issued. The method of measuring tilting rate has been modified from the previous studies, for ease of calculation. The common approach of calculating tilting rate is to calculate the slope of tilt angle vs. time graph, for a small constant interval of tilt angle. But in this study, we have kept the time interval constant instead of the change tilt angle. The change in tilt angle for every hour has been recorded as tilting rate. It has been verified if the rate is above 0.03°/h continuously for at least an hour. This has been done to minimize the errors due to noises in tilt angle readings. When the threshold is crossed, minor displacements can be expected near the jhora, which may result in very slow sinking of the uphill area. The process is very slow and happens in a duration of a few days. If the threshold value crossed is 0.1°/h, the alert is critical and there are chances for sudden slope failures, leading to landslides.

3.3. Combined Use of Field Monitoring Data and Empirical Thresholds

The abovementioned rainfall thresholds and field monitoring data were integrated in a procedure that is depicted in Figure 9 and represents a prototype of warning systems combining both techniques.
The warnings are issued only when a precursor event has happened in the study region. The precursor event is a rainfall event which is above the threshold line. From the historical data, such an event has the potential to trigger a slope failure in the region. It is not necessary that the slope failure should occur along with the event or immediately after the event. When the soil is moist, there are chances that the failure occurs within a short span, after the occurrence of the precursor event. The lag time defines the possibility of occurrence of landslides, within 72 h after the occurrence of a rainfall event, as observed during the test period. The value of 72 h has been chosen for Kalimpong, using a trial-and-error approach to minimize the missed alerts. The algorithm first checks if the rainfall threshold is crossed or not, and then checks if the peak hourly tilting rate has exceeded the tilting angle thresholds on the previous day. An alert is issued only if both the conditions are satisfied. The term Ts in Figure 9 represents the time for which alert is issued. If the rainfall forecast for a day is less than the empirical threshold, the algorithm checks if the time span is less than the lag time (72 h in this study) or not. If the lag time is not crossed, it checks for the tilt rate threshold again and issue alert if the tilt rate threshold is crossed.
The variable T5 in Figure 9 indicates the empirical rainfall threshold with 5% exceedance probability. The lag time is the time of issuing alert, after the ending of the precursor event, rainfall. The algorithm was applied to both ID and ED thresholds, with a lag time of 72 h, and the resulting performance were evaluated.

4. Results

The rainfall [71], monitoring data and landslide data for Kalimpong town from 1 July 2017 to 31 August 2020 were used for evaluation of the performance of the different methods considered. A confusion matrix was used, and some commonly used skill scores were calculated to quantify the prediction performance. The outcome or prediction by the thresholds were classified into four categories as true positives, false positives, false negatives, and true negatives for quantitative performance evaluation. The days in which the threshold predicts the possibility of a slope failure are a positive outcome. If landslide has happened on the day, it is true positive and if not; it is defined as a false positive. Negative outcomes are days on which the rainfall is below the threshold condition. If no landslides or movements are reported on the day, the outcome is a true negative and otherwise it is a false negative. By combining these four attributes, the performance of a model can be quantitatively evaluated through statistical indexes.
Efficiency is an indicator which can be used to understand the ratio of true outcomes to the total outcomes [74]. The ratio of true outcomes to the false outcomes is termed as odds ratio. For landslide forecasting models, the number of true negatives is expected to be of a higher order than the other three outcomes, and hence the terms like efficiency and odds ratio are dominated by the value of true negatives. Hence, it is important to consider this aspect and make use of terms in which the effect of true negatives is minimum. The terms sensitivity and specificity and the attributes derived from these variables are helpful in this aspect. The term sensitivity denotes the ratio of landslide events, which was correctly predicted to the total number of landslide events, and specificity denotes the ratio of non-landslide days, which was correctly predicted to the total number of non-landslide days. The value of both these variables should be 1 for a perfect model. Likelihood ratio is another term to evaluate the overall performance, the ratio of sensitivity to 1-specificity. While analysing the performance of a model, all these attributes should be evaluated thoroughly, to identify which is the best model according to the needs of the user.

4.1. Performance of Empirical Thresholds

The skill scores of the empirical thresholds are calculated using the confusion matrix and are mentioned in Table 2. The ID thresholds were 84% efficient and ED thresholds were 86% efficient in predicting true outcomes. From Table 2, it can be understood that the number of false positives predicted by the empirical thresholds are much higher than that of true positives, thus the performance of the model is highly conservative and, if used in a LEWS, it would result in a very high number of false alarms. The values of odds ratio are 5.36 and 5.98 for ID and ED thresholds, respectively. Even though the false positives are larger than true positives, odds ratio has a value greater than 1 due to the higher number of true negatives.
In our case of study, both rainfall thresholds fail to predict almost 50% of the alerts and the number of false alerts is almost twice the number of true alerts. The sensitivity of the thresholds is 0.51, which is too less for an operational LEWS. These drawbacks affect the reliability of the model and false predictions will negatively affect the response of public towards the alerts issued. Hence it is important to overcome these two disadvantages, to develop an effective LEWS.

4.2. Performance of Tilt Meters

From the data received from sensors, it was observed that the sensors near Pyarieni jhora (sensors 1, 2, and 3) have recorded variations during monsoon.
The prediction performance, by using the field monitoring data alone, is given in Table 3. The data shows a total of 273 false positives on days without landslides. The same can be observed as abrupt changes in tilt angle during the months without rainfall. The variation in tilt angles were associated with external interference and noises, and not related to slope failure. Hence it is important to correlate the data from the sensor with the in-situ moisture content and daily rainfall data for better interpretation and issuing alarms, as done by the algorithm proposed in this work and explained in Section 3.3. As can be seen from Table 3, the field monitoring data alone cannot help in developing a good early warning system, as many false alarms would be issued. Even though the use of field monitoring is more sensitive than the rainfall thresholds, the higher number of false alarms make it less specific, resulting in a lesser likelihood ration of 3.41.

4.3. Performance of Proposed Algorithm

The peak daily tilting rates were employed to enhance the forecasting performance of conventional empirical rainfall thresholds. As mentioned in Table 2 and Table 3, the major disadvantage associated with the rainfall thresholds and tilt meters when used alone, is the high number of false positives. Another important aspect is that such thresholds do not consider a lag time, which is the time between ending of the precursor event (rainfall) and the occurrence of landslides. Hence this study proposes an algorithm which uses both rainfall thresholds and field monitoring data to issue an alert to the public.
The performance statistics obtained using the proposed algorithm are mentioned in Table 4.
From Table 4, it can be observed that the performance of both ID and ED thresholds have been improved considerably, using the algorithm-based approach employed in this research. It can be understood from Table 4 that the major variation after the use of field monitoring data has happened in the decrease of false alarms. The number of false positives has been reduced from 273 to 63, using the algorithm combining ED thresholds and tilt meters.

5. Discussion

The empirical rainfall thresholds, ID and ED were observed to be in good agreement with each other. Both the thresholds had similar prediction performance with many false alerts and missed alerts. The tilt meters, when used alone, are able to predict most of the events, but with the downside of a very high number of false alerts. This drawback would limit the potential use of such methods in an operational LEWS: When false alarms are issued multiple times, the response to future warnings might get affected. Since the social response to the early warning is an important aspect in the implementation of an LEWS, it is crucial to lower the number of false alarms to the best possible extent. This is done by the combined approach proposed in this study. Figure 10 shows the variations in alerts issued by different methods considered. Another important component of the proposed warning system is the use of lag time, which helps in minimising the number of false alarms on days without rainfall. The lag time explores the effect of antecedent rainfall conditions on the initiation of landslides [75]. When compared with empirical thresholds, the combined approach shows an improved performance in terms of all the attributes. The number of true alerts is increased and the number of both false alerts and missed alerts are decreased (Figure 10), making the model more sensitive and specific. It should also be noted that there is a slight variation in the performance of the algorithm, depending upon the initial rainfall threshold used, ID or ED. The performance of ED thresholds is found to be better than that of ID thresholds when used individually and when used in algorithm. The difference comes in the decreased number of false alerts issued when ED thresholds are used (Figure 10). The reason for this improvement is the method of selection of rainfall parameters considered for the definition of thresholds. The use of automatic algorithm has considered both spatial and temporal variations in the rainfall and landslide data and has reduced the number of false alarms issued.
The very high number of false alerts provided by tilt sensors can be explained with the local conditions of the test site. Seepage happens at different locations within the old slide debris and through the joints of rocks. Erosion of debris material during rain may sometimes result in tilting the sensor, and such movements can be less critical and not always associated with landslide processes. Second, when sensors are placed in locations where chances of external factors cannot be neglected, they may record tilt angles which are not related to the slope failure, but due to micro-scale deformations not connected with slope instability or can be caused by external factors like disturbance from animals or human activities. This strengthens the initial idea of the proposed approach, for which data from the sensors must be interpreted in connection with the rainfall data, to take decisions regarding early warning.
Figure 10 depicts that the number of false alarms issued by the field monitoring has been reduced considerably by the combined approach. But this comes with the cost of missing some of the true alerts, which have happened long after the occurrence of the precursor event (rainfall). Such events have started during the occurrence of rainfall, but have continued at very slow rates a few more days after the rainfall. The algorithm fails to predict such events if the lag time is more than 72 h. Increasing the lag time further will result in a greater number of false alerts. Thus, the use of tilt meters alone is a more sensitive approach with less missed alerts, but the higher number of false alerts makes it unsuitable for direct use in LEWS. The combined approach has a better overall performance with higher efficiency and likelihood ratio. The method proposed in this study is conceived to issue alerts with a lead time of 24 h. That lead time would allow the local community to implement some counteractions, which should be carefully evaluated in subsequent stages of the research, based also on the expected risk scenarios (yet to be defined) and the reliability of the warning system. The prototype LEWS should be tested first for a longer duration and detailed action plans should be decided, after considering different risk scenarios.
The development of LEWS demands strong scientific research in data collection, interpretation, and integration of different modules. The proposed algorithm overcomes the disadvantages associated with both rainfall thresholds and field monitoring data. The method is simple and easy to export, with historical rainfall and landslide data and a network of slope monitoring sensors. Cost of installation of a large number of sensors is a major concern for developing countries like India, hence a cost-effective approach is used in this study. The use of MEMS sensors along with the rainfall thresholds is thus a simple and economical approach for the prediction of landslide events.
Indeed, the methodology is far to get perfect predictions and needs further improvements: A purely statistical rainfall threshold approach is too simplistic to fully cope with the complex geological setting of this area. The combined use of instrumental monitoring helps getting some more insights on the behaviour of the slopes, but the error rates returned in the validation demonstrated that the complex geology of the area remains not adequately caught by the model. This conclusion opens other research directions to be explored in the near future such as: Extending the monitoring network to gather more data at different locations; or combining the rainfall threshold method with a zonation of landslide susceptibility [15] based on morphometry and geotechnical characteristics (e.g., weathering degree, orientation of joints, residual strength of the soil) of the area. Such improvements could help to better relate the triggering factors (i.e., rainfall) with the variability of the geological setting of the area.
Another limitation of the methodology lies in the intrinsic nature of empirical rainfall threshold method, which does not directly consider the failure mechanism. Since the major rock type of the study area is schist, the chances for purely translational failure cannot be neglected. Theoretically speaking, that leaves open the possibility of tilt angles undetected by the sensors. Since such circumstances have not happened during the study period, the response of tilt sensors to purely translational failure, remains not yet examined.
The proposed system should be considered prototypal as a high rate of errors is still present; nonetheless the testes improvements could be evaluated ameliorative with respect to the rainfall thresholds and tilt sensors used alone, thus they represent a step in the right direction for a more balanced hazard management in the study area. The methodology proposed in this study, which is the first of its kind for the Darjeeling Himalayas, can be easily exported to the other parts of the world, to test the possibility to improve the forecasting performance of existing empirical thresholds in different hydrological and geological settings.
To understand the applicability of the proposed algorithm for different meteo-hydrological settings, experimental investigations can also be conducted on laboratory scale. This could verify the possibility of exporting the model to different regions, by simulating different rainfall conditions on different types of soil and topography [76,77]. The prototype developed in this study is a site-specific case from the Pyarieni Jhora of Kalimpong, where failure occurs at a very slow rates and observed as cracks and displacements in the uphill areas. In a laboratory scale study, it is possible to simulate different types of landslides and define different thresholds for each case. Such studies can help in the calibration of model and critical tilt rate for different conditions of precursors and slope failures and can aid in field scale studies.

6. Conclusions

The rainfall thresholds defined for Kalimpong region in Darjeeling Himalayas has been evaluated quantitatively, using the recent rainfall and landslide data from 2017 to 2020 monsoons. The analysis proves that the large number of false alarms makes the conventional empirical rainfall thresholds unsuitable for use in an LEWS. Even though the thresholds predict the most landslides events correctly, the false alarms affect the reliability of the system and makes difficult to issue a trustworthy early warning to the public. To overcome this limitation, real-time field monitoring data from Kalimpong was utilized to enhance the forecasting performance of empirical thresholds.
The field monitoring system consists of MEMS units integrating a tilt sensor, a soil moisture meter and a real-time wireless transmitter. The sensors which recorded peak variations were used to modify the empirical rainfall thresholds, using an algorithm that combines rainfall thresholds, lag time, and recorded tilting rates. According to the algorithm, the rainfall thresholds were taken as a first line of action and the criticality was verified using the field monitoring data, to issue alerts. The results showed that the combined approach has an efficiency of 92% and a likelihood ratio of 11.33, which makes them a more suitable part of the proposed LEWS for Kalimpong respect to both the original version based only on rainfall thresholds and a monitoring system based solely on MEMS. The proposed method is also found to have better performance than an early warning system using field monitoring data only. The algorithm overcomes the limitations of both rainfall thresholds and field monitoring data when used separately and can be used by the authorities to warn the public regarding a possible slope failure. The study discusses the first prototype of LEWS for this highly vulnerable landslide zone of Kalimpong and is found to be promising to be implemented as an operational LEWS.

Author Contributions

Conceptualization, M.T.A. and N.S.; formal analysis, M.T.A.; supervision, N.S., B.P., B.T.P., and S.S.; writing—literature review and original draft, M.T.A. and M.A.B.; writing—review and editing, N.S., B.P., and S.S. All authors have read and agreed to the published version of the manuscript.

Funding

The work has been carried out with the financial support of Department of Science & Technology (DST), New Delhi, for funding the research project Grant No. (NRDMS/02/31/015(G)) and Università degli Studi di Firenze under the Grant 58517_INTERNAZIONALIZZAZIONE.

Acknowledgments

The authors express our sincere gratitude to both the funding agencies. We are also thankful to Praful Rao and Save the Hills organization for their constant support throughout the study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dowling, C.A.; Santi, P.M. Debris flows and their toll on human life: A global analysis of debris-flow fatalities from 1950 to 2011. Nat. Hazards 2014, 71, 203–227. [Google Scholar] [CrossRef]
  2. Petley, D. Global patterns of loss of life from landslides. Geology 2012, 40, 927–930. [Google Scholar] [CrossRef]
  3. Devoli, G.; Tiranti, D.; Cremonini, R.; Sund, M.; Boje, S. Comparison of landslide forecasting services in Piedmont (Italy) and Norway, illustrated by events in late spring 2013. Nat. Hazards 2018, 19, 1351–1372. [Google Scholar] [CrossRef] [Green Version]
  4. Baum, R.L.; Godt, J.W. Early warning of rainfall-induced shallow landslides and debris flows in the USA. Landslides 2010, 7, 259–272. [Google Scholar] [CrossRef]
  5. Segoni, S.; Piciullo, L.; Gariano, S.L. A review of the recent literature on rainfall thresholds for landslide occurrence. Landslides 2018, 15, 1483–1501. [Google Scholar] [CrossRef]
  6. Fusco, F.; De Vita, P.; Mirus, B.B.; Baum, R.L.; Allocca, V.; Tufano, R.; Clemente, E.D.I.; Calcaterra, D. Physically based estimation of rainfall thresholds triggering shallow landslides in volcanic slopes of Southern Italy. Water 2019, 11, 1915. [Google Scholar] [CrossRef] [Green Version]
  7. De Vita, P.; Fusco, F.; Napolitano, E.; Tufano, R. Physically-based models for estimating rainfall triggering debris flows in Campania (Southern Italy). In Advancing Culture of Living with Landslides; Springer International Publishing: Cham, Switzerland, 2017; pp. 289–297. [Google Scholar]
  8. Napolitano, E.; Fusco, F.; Baum, R.L.; Godt, J.W.; De Vita, P. Effect of antecedent-hydrological conditions on rainfall triggering of debris flows in ash-fall pyroclastic mantled slopes of Campania (southern Italy). Landslides 2016, 13, 967–983. [Google Scholar] [CrossRef]
  9. Guzzetti, F.; Peruccacci, S.; Rossi, M.; Stark, C.P. Rainfall thresholds for the initiation of landslides in central and southern Europe. Meteorol. Atmos. Phys. 2007, 98, 239–267. [Google Scholar] [CrossRef]
  10. Caine, N. The rainfall intensity-duration control of shallow landslides and debris flows: An update. Geogr. Ann. Ser. A Phys. Geogr. 1980, 62, 23–27. [Google Scholar]
  11. Aleotti, P. A warning system for rainfall-induced shallow failures. Eng. Geol. 2004, 73, 247–265. [Google Scholar] [CrossRef]
  12. Bordoni, M.; Corradini, B.; Lucchelli, L.; Valentino, R.; Bittelli, M.; Vivaldi, V.; Meisina, C. Empirical and physically based thresholds for the occurrence of shallow landslides in a prone area of Northern Italian Apennines. Water 2019, 11, 2653. [Google Scholar] [CrossRef] [Green Version]
  13. Berti, M.; Martina, M.L.V.; Franceschini, S.; Pignone, S.; Simoni, A.; Pizziolo, M. Probabilistic rainfall thresholds for landslide occurrence using a Bayesian approach. J. Geophys. Res. Earth Surf. 2012, 117, 1–20. [Google Scholar] [CrossRef] [Green Version]
  14. Guzzetti, F.; Peruccacci, S.; Rossi, M.; Stark, C.P. The rainfall intensity-duration control of shallow landslides and debris flows: An update. Landslides 2008, 5, 3–17. [Google Scholar] [CrossRef]
  15. Brunetti, M.T.; Peruccacci, S.; Rossi, M.; Luciani, S.; Valigi, D.; Guzzetti, F. Rainfall thresholds for the possible occurrence of landslides in Italy. Nat. Hazards Earth Syst. Sci. 2010, 10, 447–458. [Google Scholar] [CrossRef]
  16. Dahal, R.K.; Hasegawa, S. Representative rainfall thresholds for landslides in the Nepal Himalaya. Geomorphology 2008, 100, 429–443. [Google Scholar] [CrossRef]
  17. Peruccacci, S.; Brunetti, M.T.; Luciani, S.; Vennari, C.; Guzzetti, F. Lithological and seasonal control on rainfall thresholds for the possible initiation of landslides in central Italy. Geomorphology 2012, 139–140, 79–90. [Google Scholar] [CrossRef]
  18. Melillo, M.; Brunetti, M.T.; Peruccacci, S.; Gariano, S.L.; Guzzetti, F. An Algorithm for the objective reconstruction of rainfall events responsible for landslides. In Landslide Dynamics: ISDR-ICL Landslide Interactive Teaching Tools; Fundamentals, Mapping and Monitoring; Springer International Publishing: Cham, Switzerland, 2014; Volume 1, pp. 433–447. [Google Scholar]
  19. Melillo, M.; Brunetti, M.T.; Peruccacci, S.; Gariano, S.L.; Guzzetti, F. Rainfall thresholds for the possible landslide occurrence in Sicily (Southern Italy) based on the automatic reconstruction of rainfall events. Landslides 2016, 13, 165–172. [Google Scholar] [CrossRef]
  20. Zhao, B.; Dai, Q.; Han, D.; Dai, H.; Mao, J.; Zhuo, L. Probabilistic thresholds for landslides warning by integrating soil moisture conditions with rainfall thresholds. J. Hydrol. 2019, 574, 276–287. [Google Scholar] [CrossRef]
  21. Segoni, S.; Rossi, G.; Rosi, A.; Catani, F. Landslides triggered by rainfall: A semi-automated procedure to define consistent intensity-duration thresholds. Comput. Geosci. 2014, 63, 123–131. [Google Scholar] [CrossRef]
  22. Teja, T.S.; Dikshit, A.; Satyam, N. Determination of rainfall thresholds for landslide prediction using an algorithm-based approach: Case study in the Darjeeling Himalayas, India. Geoscience 2019, 9, 302. [Google Scholar] [CrossRef] [Green Version]
  23. Abraham, M.T.; Satyam, N.; Rosi, A.; Pradhan, B.; Segoni, S. The selection of rain gauges and rainfall parameters in estimating intensity-duration thresholds for landslide occurrence: Case study from Wayanad (India). Water 2020, 12, 1000. [Google Scholar] [CrossRef] [Green Version]
  24. Abraham, M.T.; Pothuraju, D.; Satyam, N. Rainfall thresholds for prediction of landslides in Idukki, India: An empirical approach. Water 2019, 11, 2113. [Google Scholar] [CrossRef] [Green Version]
  25. Gariano, S.L.; Melillo, M.; Peruccacci, S.; Brunetti, M.T. How much does the rainfall temporal resolution affect rainfall thresholds for landslide triggering? Nat. Hazards 2020, 100, 655–670. [Google Scholar] [CrossRef] [Green Version]
  26. Segoni, S.; Rosi, A.; Fanti, R.; Gallucci, A.; Monni, A.; Casagli, N. A regional-scale landslide warning system based on 20 years of operational experience. Water 2018, 10, 1297. [Google Scholar] [CrossRef] [Green Version]
  27. Rosi, A.; Canavesi, V.; Segoni, S.; Dias Nery, T.; Catani, F.; Casagli, N. Landslides in the mountain region of rio de Janeiro: A proposal for the semi-automated definition of multiple rainfall thresholds. Geosciences 2019, 9, 203. [Google Scholar] [CrossRef] [Green Version]
  28. Nikolopoulos, E.I.; Crema, S.; Marchi, L.; Marra, F.; Guzzetti, F.; Borga, M. Impact of uncertainty in rainfall estimation on the identification of rainfall thresholds for debris flow occurrence. Geomorphology 2014, 221, 286–297. [Google Scholar] [CrossRef]
  29. Bogaard, T.; Greco, R. Invited perspectives: Hydrological perspectives on precipitation intensity-duration thresholds for landslide initiation: Proposing hydro-meteorological thresholds. Nat. Hazards Earth Syst. Sci. 2018, 18, 31–39. [Google Scholar] [CrossRef] [Green Version]
  30. Wicki, A.; Lehmann, P.; Hauck, C.; Seneviratne, S.I.; Waldner, P.; Stähli, M. Assessing the potential of soil moisture measurements for regional landslide early warning. Landslides 2020, 17, 1881–1896. [Google Scholar] [CrossRef] [Green Version]
  31. Mirus, B.B.; Becker, R.E.; Baum, R.L.; Smith, J.B. Integrating real-time subsurface hydrologic monitoring with empirical rainfall thresholds to improve landslide early warning. Landslides 2018, 15, 1909–1919. [Google Scholar] [CrossRef]
  32. Glade, T.; Crozier, M.; Smith, P. Applying probability determination to refine landslide-triggering rainfall thresholds using an empirical “Antecedent Daily Rainfall Model.”. Pure Appl. Geophys. 2000, 157, 1059–1079. [Google Scholar] [CrossRef]
  33. Rosi, A.; Segoni, S.; Canavesi, V.; Monni, A.; Gallucci, A.; Casagli, N. Definition of 3D rainfall thresholds to increase operative landslide early warning system performances. Landslides 2020, 1–13. [Google Scholar] [CrossRef]
  34. Martelloni, G.; Segoni, S.; Lagomarsino, D.; Fanti, R.; Catani, F. Snow accumulation/melting model (SAMM) for integrated use in regional scale landslide early warning systems. Hydrol. Earth Syst. Sci. 2013, 17, 1229–1240. [Google Scholar] [CrossRef] [Green Version]
  35. Wei, X.; Fan, W.; Cao, Y.; Chai, X.; Bordoni, M.; Meisina, C.; Li, J. Integrated experiments on field monitoring and hydro-mechanical modeling for determination of a triggering threshold of rainfall-induced shallow landslides. A case study in Ren River catchment, China. Bull. Eng. Geol. Environ. 2020, 79, 513–532. [Google Scholar] [CrossRef]
  36. Pecoraro, G.; Calvello, M.; Piciullo, L. Monitoring strategies for local landslide early warning systems. Landslides 2018, 16, 213–231. [Google Scholar] [CrossRef]
  37. Casagli, N.; Catani, F.; Del Ventisette, C.; Luzi, G. Monitoring, prediction, and early warning using ground-based radar interferometry. Landslides 2010, 7, 291–301. [Google Scholar] [CrossRef]
  38. Yin, Y.; Zheng, W.; Liu, Y.; Zhang, J.; Li, X. Integration of GPS with InSAR to monitoring of the Jiaju landslide in Sichuan, China. Landslides 2010, 7, 359–365. [Google Scholar] [CrossRef]
  39. Derron, M.H.; Jaboyedoff, M. LIDAR and DEM techniques for landslides monitoring and characterization. Nat. Hazards Earth Syst. Sci. 2010, 10, 1877–1879. [Google Scholar] [CrossRef]
  40. Joyce, K.E.; Samsonov, S.V.; Levick, S.R.; Engelbrecht, J.; Belliss, S. Mapping and monitoring geological hazards using optical, LiDAR, and synthetic aperture RADAR image data. Nat. Hazards 2014, 73, 137–163. [Google Scholar] [CrossRef]
  41. Jaboyedoff, M.; Oppikofer, T.; Abella, A.; Derron, M.; Loye, A.; Metzger, R.; Pedrazzini, A. Use of LIDAR in landslide investigations: A review. Nat. Hazards 2010, 61, 5–28. [Google Scholar] [CrossRef] [Green Version]
  42. Raspini, F.; Bianchini, S.; Ciampalini, A.; Del Soldato, M.; Montalti, R.; Solari, L.; Tofani, V.; Casagli, N. Persistent Scatterers continuous streaming for landslide monitoring and mapping: The case of the Tuscany region (Italy). Landslides 2019, 16, 2033–2044. [Google Scholar] [CrossRef] [Green Version]
  43. Angeli, M.G.; Pasuto, A.; Silvano, S. A critical review of landslide monitoring experiences. Eng. Geol. 2000, 55, 133–147. [Google Scholar] [CrossRef]
  44. Bhattacharya, A.; Mukherjee, K. Review on InSAR based displacement monitoring of Indian Himalayas: Issues, challenges and possible advanced alternatives. Geocarto Int. 2017, 32, 298–321. [Google Scholar] [CrossRef]
  45. Meisina, C.; Zucca, F.; Notti, D.; Colombo, A.; Cucchi, A.; Savio, G.; Giannico, C.; Bianchi, M. Geological interpretation of PSInSAR data at regional scale. Sensors 2008, 8, 7469–7492. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Abraham, M.T.; Satyam, N.; Pradhan, B.; Alamri, A.M. IoT-based geotechnical monitoring of unstable slopes for landslide early warning in the Darjeeling Himalayas. Sensors 2020, 20, 2611. [Google Scholar] [CrossRef]
  47. Towhata, I.; Uchimura, T.; Seko, I.; Wang, L. Monitoring of unstable slopes by MEMS tilting sensors and its application to early warning. IOP Conf. Ser. Earth Environ. Sci. 2015, 26, 012049. [Google Scholar] [CrossRef] [Green Version]
  48. Zhang, Y.; Tang, H.; Li, C.; Lu, G.; Cai, Y.; Zhang, J.; Tan, F. Design and testing of a flexible inclinometer probe for model tests of landslide deep displacement measurement. Sensors 2018, 18, 224. [Google Scholar] [CrossRef] [Green Version]
  49. Chen, Y.; Irfan, M.; Uchimura, T.; Zhang, K. Feasibility of using elastic wave velocity monitoring for early warning of rainfall-induced slope failure. Sensors 2018, 18, 997. [Google Scholar] [CrossRef] [Green Version]
  50. Tao, S.; Uchimura, T.; Fukuhara, M.; Tang, J.; Chen, Y.; Huang, D. Evaluation of soil moisture and shear deformation based on compression wave velocities in a shallow slope surface layer. Sensors 2019, 19, 3406. [Google Scholar] [CrossRef] [Green Version]
  51. Le Cozannet, G.; Hosford, S.; Douglas, J.; Serrano, J.J.; Coraboeuf, D.; Comte, J. Connecting hazard analysts and risk managers to sensor information. Sensors 2008, 8, 3932–3937. [Google Scholar] [CrossRef]
  52. Uchimura, T.; Towhata, I.; Wang, L. Miniature ground inclinometer for slope monitoring. In Proceedings of the 14th Asian Regional Conference on Soil Mechanics and Geotechnical Engineering, Hong Kong, China, 23–27 May 2011. [Google Scholar]
  53. Huang, Y.; Sai, A.; Vasan, S.; Doraiswami, R.; Osterman, M.; Pecht, M. MEMS reliability review. IEEE Trans. Device Mater. Reliabil. 2012, 12, 482–493. [Google Scholar] [CrossRef] [Green Version]
  54. Jung, H.; Kim, C.J.; Kong, S.H. An optimized MEMS-based electrolytic tilt sensor. Sens. Actuators A Phys. 2007, 139, 23–30. [Google Scholar] [CrossRef]
  55. Yang, Z.; Cai, H.; Shao, W.; Huang, D.; Uchimura, T.; Lei, X.; Tian, H.; Qiao, J. Clarifying the hydrological mechanisms and thresholds for rainfall-induced landslide: In situ monitoring of big data to unsaturated slope stability analysis. Bull. Eng. Geol. Environ. 2019, 78, 2139–2150. [Google Scholar] [CrossRef]
  56. Yang, Z.; Shao, W.; Qiao, J.; Huang, D.; Tian, H.; Lei, X.; Uchimura, T. A multi-source early warning system of MEMS based wireless monitoring for rainfall-induced landslides. Appl. Sci. 2017, 7, 1234. [Google Scholar] [CrossRef] [Green Version]
  57. Froude, M.J.; Petley, D.N. Global fatal landslide occurrence from 2004 to 2016. Nat. Hazards Earth Syst. Sci. 2018, 18, 2161–2181. [Google Scholar] [CrossRef] [Green Version]
  58. Dikshit, A.; Satyam, N. Rainfall thresholds for the prediction of landslides using empirical methods in Kalimpong, Darjeeling, India. In Proceedings of the JTC1 Workshop on Advances in Landslide Understanding, Barcelona, Spain, 24–26 May 2017. [Google Scholar]
  59. Dahal, R.K. Rainfall-induced Landslides in Nepal. Int. J. Jpn. Eros. Control Eng. 2012, 5, 1–8. [Google Scholar] [CrossRef] [Green Version]
  60. Gabet, E.J.; Burbank, D.W.; Putkonen, J.K.; Pratt-sitaula, B.A.; Ojha, T. Rainfall thresholds for landsliding in the Himalayas of Nepal. Geomorphology 2004, 63, 131–143. [Google Scholar] [CrossRef]
  61. Dikshit, A.; Sarkar, R.; Pradhan, B.; Segoni, S.; Alamri, A.M. Rainfall induced landslide studies in Indian Himalayan region: A critical review. Appl. Sci. 2020, 10, 2466. [Google Scholar] [CrossRef] [Green Version]
  62. Abraham, M.T.; Satyam, N.; Reddy, S.K.P.; Pradhan, B. Runout modeling and calibration of friction parameters of Kurichermala debris flow, India. Landslides 2020, 1–8. [Google Scholar] [CrossRef]
  63. Dikshit, A.; Satyam, D.N. Estimation of rainfall thresholds for landslide occurrences in Kalimpong, India. Innov. Infrastruct. Solut. 2018, 3, 24. [Google Scholar] [CrossRef]
  64. Dikshit, A.; Satyam, N. Probabilistic rainfall thresholds in Chibo, India: Estimation and validation using monitoring system. J. Mt. Sci. 2019, 16, 870–883. [Google Scholar] [CrossRef]
  65. Abraham, M.T.; Satyam, N.; Pradhan, B.; Alamri, A.M. Forecasting of landslides using rainfall severity and soil wetness: A probabilistic approach for Darjeeling Himalayas. Water 2020, 12, 804. [Google Scholar] [CrossRef] [Green Version]
  66. Abraham, M.T.; Satyam, N.; Kushal, S.; Rosi, A.; Pradhan, B.; Segoni, S. Rainfall threshold estimation and landslide forecasting for Kalimpong, India using SIGMA model. Water 2020, 12, 1195. [Google Scholar] [CrossRef] [Green Version]
  67. Dikshit, A.; Sarkar, R.; Satyam, N. Probabilistic approach toward Darjeeling Himalayas landslides—A case study. Cogent Eng. 2018, 5, 1–11. [Google Scholar] [CrossRef]
  68. Dortch, J.M.; Owen, L.A.; Haneberg, W.C.; Caffee, M.W.; Dietsch, C.; Kamp, U. Nature and timing of large landslides in the Himalaya and Transhimalaya of northern India. Quat. Sci. Rev. 2009, 28, 1037–1054. [Google Scholar] [CrossRef]
  69. Dikshit, A.; Satyam, D.N.; Towhata, I. Early warning system using tilt sensors in Chibo, Kalimpong, Darjeeling Himalayas, India. Nat. Hazards 2018, 94, 727–741. [Google Scholar] [CrossRef]
  70. Mukherjee, A.; Mitra, A. Geotechnical study of mass movements along the kalimpong approach road in the Eastern Himalayas. Indian J. Geol. 2001, 73, 271–279. [Google Scholar]
  71. Save The Hills Blog. Available online: http://savethehills.blogspot.com/ (accessed on 15 September 2020).
  72. Chuo Kaihatsu Corporation—Company Profile. Available online: http://www.ckcnet.co.jp/ (accessed on 18 October 2020).
  73. Lin, W.; Shunsaku, N.; Taro, U.; Ikuo, T.; Ling, S.; Shangning, T. An early warning system of unstable slopes by multi-point MEMS tilting sensors and water contents. In Advancing Culture of Living with Landslides; Springer International Publishing: Cham, Switzerland, 2017; pp. 147–154. [Google Scholar]
  74. Lagomarsino, D.; Segoni, S.; Rosi, A.; Rossi, G.; Battistini, A.; Catani, F.; Casagli, N. Quantitative comparison between two different methodologies to define rainfall thresholds for landslide forecasting. Nat. Hazards Earth Syst. Sci. 2015, 15, 2413–2423. [Google Scholar] [CrossRef] [Green Version]
  75. Melillo, M.; Brunetti, M.T.; Peruccacci, S.; Gariano, S.L.; Roccati, A.; Guzzetti, F. A tool for the automatic calculation of rainfall thresholds for landslide occurrence. Env. Model. Softw. 2018, 105, 230–243. [Google Scholar] [CrossRef]
  76. Ivanov, V.; Arosio, D.; Tresoldi, G.; Hojat, A.; Zanzi, L.; Papini, M.; Longoni, L. Investigation on the role of water for the stability of shallow landslides—Insights from experimental tests. Water 2020, 12, 1203. [Google Scholar] [CrossRef] [Green Version]
  77. Cogan, J.; Gratchev, I. A study on the effect of rainfall and slope characteristics on landslide initiation by means of flume tests. Landslides 2019, 16, 2369–2379. [Google Scholar] [CrossRef]
Figure 1. Location details of study area (a) India, (b) West Bengal and (c) Geological map of Kalimpong town (Modified after [22]).
Figure 1. Location details of study area (a) India, (b) West Bengal and (c) Geological map of Kalimpong town (Modified after [22]).
Water 12 03453 g001
Figure 2. Daily rainfall data for Kalimpong (July 2017–August 2020).
Figure 2. Daily rainfall data for Kalimpong (July 2017–August 2020).
Water 12 03453 g002
Figure 3. Empirical thresholds derived for Kalimpong along with the rainfall events considered for derivation of thresholds (a) ID threshold [58] and (b) ED threshold [22].
Figure 3. Empirical thresholds derived for Kalimpong along with the rainfall events considered for derivation of thresholds (a) ID threshold [58] and (b) ED threshold [22].
Water 12 03453 g003
Figure 4. Schematic arrangement of sensor unit [46].
Figure 4. Schematic arrangement of sensor unit [46].
Water 12 03453 g004
Figure 5. Configuration of a sensor unit.
Figure 5. Configuration of a sensor unit.
Water 12 03453 g005
Figure 6. Location of sensors in Chibo, Kalimpong.
Figure 6. Location of sensors in Chibo, Kalimpong.
Water 12 03453 g006
Figure 7. The recordings from sensors along with daily rainfall data and displacement periods (a) tilting angle in the direction parallel to slope, and (b) tilting angle in the direction perpendicular to slope.
Figure 7. The recordings from sensors along with daily rainfall data and displacement periods (a) tilting angle in the direction parallel to slope, and (b) tilting angle in the direction perpendicular to slope.
Water 12 03453 g007
Figure 8. Cracks developed in road near Pyarienii jhora during 2020 monsoon (a) on 27 July 2020 and (b) on 6 August 2020.
Figure 8. Cracks developed in road near Pyarienii jhora during 2020 monsoon (a) on 27 July 2020 and (b) on 6 August 2020.
Water 12 03453 g008
Figure 9. Proposed algorithm for issuing slope failure alert in Chibo area, Kalimpong.
Figure 9. Proposed algorithm for issuing slope failure alert in Chibo area, Kalimpong.
Water 12 03453 g009
Figure 10. Comparison between different methods considered.
Figure 10. Comparison between different methods considered.
Water 12 03453 g010
Table 1. Average monthly rainfall during July 2017–August 2020.
Table 1. Average monthly rainfall during July 2017–August 2020.
MonthJanuaryFebruaryMarchAprilMayJuneJulyAugustSeptemberOctoberNovemberDecember
Rainfall (mm)14.6030.129.1074.87105.87292.20675.85512.68347.5010.330.002.00
Table 2. Prediction performance of empirical thresholds.
Table 2. Prediction performance of empirical thresholds.
Statistical AttributesID Threshold [58]ED Threshold [22]
a = True positives5050
b = False positives134118
c = False negatives4848
d = True negatives926942
Efficiency = (a + d)/(a + b + c + d)0.840.86
Odds ratio = (a + d)/(b + c)5.365.98
Specificity = d/(b + d)0.870.89
Sensitivity = a/(a + c)0.510.51
Likelihood ratio = Sensitivity/(1 − Specificity)4.044.58
Table 3. Prediction performance of field monitoring data.
Table 3. Prediction performance of field monitoring data.
Statistical AttributesField Monitoring
True positives86
False positives273
False negatives12
True negatives787
Efficiency 0.75
Odds ratio 3.06
Specificity 0.74
Sensitivity 0.88
Likelihood ratio 3.41
Table 4. Prediction performance of algorithm-based approach, using empirical rainfall thresholds.
Table 4. Prediction performance of algorithm-based approach, using empirical rainfall thresholds.
Statistical AttributesID + TiltVariation from IDED + TiltVariation from ED
True positives66+1666+16
False positives71−6363−55
False negatives32−1632−16
True negatives989+63997+55
Efficiency 0.91+0.070.92+0.06
Odds ratio 10.24+4.8811.19+5.21
Specificity 0.93+0.060.94+0.05
Sensitivity 0.67+0.160.67+0.16
Likelihood ratio 10.05+6.0211.33+6.75
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abraham, M.T.; Satyam, N.; Bulzinetti, M.A.; Pradhan, B.; Pham, B.T.; Segoni, S. Using Field-Based Monitoring to Enhance the Performance of Rainfall Thresholds for Landslide Warning. Water 2020, 12, 3453. https://doi.org/10.3390/w12123453

AMA Style

Abraham MT, Satyam N, Bulzinetti MA, Pradhan B, Pham BT, Segoni S. Using Field-Based Monitoring to Enhance the Performance of Rainfall Thresholds for Landslide Warning. Water. 2020; 12(12):3453. https://doi.org/10.3390/w12123453

Chicago/Turabian Style

Abraham, Minu Treesa, Neelima Satyam, Maria Alexandra Bulzinetti, Biswajeet Pradhan, Binh Thai Pham, and Samuele Segoni. 2020. "Using Field-Based Monitoring to Enhance the Performance of Rainfall Thresholds for Landslide Warning" Water 12, no. 12: 3453. https://doi.org/10.3390/w12123453

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

Article Metrics

Back to TopTop