The article presents the rainfall threshold development framework and its employment for sediment-related disaster from Malaysia’s perspective. The open issues and challenges on deriving the rain threshold are discussed in details. The outcomes of this review could serve as references for the application in sediment-related disaster management worldwide, especially for developing countries.
1. Introduction
Landslides, debris and mudflows (DMF), and rockfalls are types of sediment-related disasters whereby rainfall is the main triggering factor [
1] The debris flow slide downslope under the influence of gravity consists of a mixture of coarse material, air and water to create a slurry and its ability to wash away anything along its path in high velocity depends on the slope gradient. Mudflow, also known as earth flow, contains at least 50% sand, silt and clay particles and the materials flow downslope rapidly. It is often referred to as a mud flood when the mixture reaches the equilibrium state, as well as when the landslide dam collapses. These sediment-related disasters occur in a sporadic and unpredictable manner. However, based on the extent, regularity and consequences of historical events, forecasting these phenomena has become a global subject of considerable interest. As Guzzetti et al. [
2] defined, the rainfall threshold is the condition or amount of rainfall that is likely to initiate sediment-related disasters. Therefore, it is agreed that the rainfall threshold is the most significant tool to predict impending occurrences, and this has been enhanced to facilitate an adaptable and reliable prediction [
3].
Referring to the work of Piciullo et al. [
4], the rainfall threshold can be implemented as a warning model for landslide early warning systems (LEWS). LEWS are one of the passive countermeasures in the framework of landslide risk management at both local and regional scales [
5]. Local LEWS (Lo-LEWS) are more specific in the monitoring of a single landslide at the slope scale, while territorial LEWS (Te-LEWS) monitor multiple sites at a regional scale.
The issues with and challenges of assessing the rain threshold for the alert, the forecast and the warning rain greatly concern the development of the landslide model. The effectiveness of the warning model also depends on the false alert rate and time for evacuation. Therefore, this kind of effective evaluation would only be possible once the warning model is well established with the specified or identified threshold. However, the hydrological threshold does not provide the spatial information and impact of the upcoming event.
In Malaysia, sediment-related disasters are associated with hilly and mountainous landscapes, coupled with intense or prolonged rainfall [
6,
7,
8,
9,
10,
11,
12,
13].
Figure 1 graphically summarizes the major disaster events recorded in Malaysian history. The first-ever massive landslip took place in 1961, which injured 35 people and claimed another 16 lives in Ringlet, Cameron Highlands. The 1996 debris flow incident in Keningau, Sabah, was perhaps the deadliest disaster with 302 recorded deaths and 4925 houses destroyed. The threat of landslides was further witnessed as recently as 2020, when three people were killed. Over the past decades, DMF and landslide disasters have occurred on a frequent basis with a total loss of more than 600 lives (
Figure 1 and
Table 1), and several research efforts have been carried out, focusing on the high-risk areas.
Figure 1. Location of sediment-related disaster events (involving death) in Malaysia.
Table 1. Details of the historical summary of sediment-related disaster events in Malaysia.
2. Rainfall Threshold
A warning model for LEWS can be based on physical or empirical models. The empirical model is mainly based on the correlations between rainfall events that have caused landslides [
14], whereas the physical method involves physically computing the force balance between shear-force and shear-force resistance at each slope. Papa et al. [
15] highlighted that the established physical model requires instrumentation to compute certain inputs such matric suction within the soil layers, hydraulic conductivity, soil moisture and pore water pressure. It also includes the infinite slope stability analysis, of which the main parameter is the groundwater table that, when fully saturated, leads to failure. Therefore, the physical-based model is highly suitable for a site-specific slope, but it is very costly if one has to apply it to the regional threshold. The conceptual model is also considered as a process-based model that represent the physical analogy, representing the quantitative understanding between rainfall and landslides. It is more widespread than the instrumentation installation requirement. The example of the conceptual model is the leaky barrel model [
16] and the soil water index (SWI) [
17,
18]. On the other hand, the work of Yamazaki et al. [
19] has successfully identified the critical rainfall condition to simulate the disaster temporally and spatially. The advantage of this method is that it takes into account that the event may occur during the observation with no rain forecast. The rainfall condition curve is segregated according to their slope classification or spatial condition when the current rainfall condition breaches the curve, hence indicating when the disaster is likely to occur.
While the empirical model is more widespread, the application of empirical rainfall thresholds has been in place since 1987, with many studies conducted in many countries, especially in Italy, and central and southern Europe [
20]. There are several types and characteristics of rainfall thresholds, namely the intensity–duration (I-D), accumulated rainfall–duration (E-D), intensity–accumulated rainfall (I-E), accumulated rainfall–accumulated rainfall (E-E) and antecedent or working rainfall or rainfall index. The most common and pioneering is the I-D threshold established by Caine [
21] Most of the threshold correlating the rainfall characteristic to landslide events is presented in terms of the power law relationship [
22]. The main goal in the empirical model development is to identify the rain that has triggered the disaster despite other additional factors, if present. Therefore, the definition of rainfall intensity plays the most crucial part in this manner [
23]. Rainfall intensity refers to the rainfall depth over a certain duration, depending on the observation period. The timescale is normally in an hourly format [
3]. For a short period, it may represent the current intensity of the rainfall during the occurrence of the disaster. Likewise, for a long period, it may represent the average value of the rainfall over hours, i.e., 2, 3 or 24 h. The definition should be consistent during the construction of the numerical rainfall characteristic of the historical events, and the forecasting rainfall (forecasting events) and the observed rainfall (monitoring events). Due to the diversity of the definition, one should state clearly the definition of rainfall intensity in their work for others to compare and adopt the method accordingly. The disaster is likely to occur when the current rainfall characteristic based on the observed rainfall data (or real-time data) are beyond the rainfall characteristic that resulted in the disaster in the past, referred to as the rainfall threshold. The framework is adopted in the early warning system as a warning model that the impending disaster can be forecasted, in order to alert the respective stakeholder (authority or community).
Figure 2 illustrates the conceptual description of the I-D threshold under ideal conditions, where there is a clear separation between the sediment-related disaster-triggering and non-triggering rainfall conditions [
24].
Figure 2. Conceptual description of intensity–duration of rainfall threshold under ideal conditions [
24].
Both the physical and empirical model are significant and can be used as warning models employed in any early warning system in order to reduce the disaster risk. However, the chosen model should be tailored according to the local environment, climate and geological condition of the area under surveillance. Since sediment-related disaster is caused by many factors—either spatial or climatic factors and types of behavior—it is pertinent to compile all of the important information related to the possible causes [
25]. In general, the process or step of the threshold development is illustrated in the figure below.
Based on
Figure 3, the general steps of rainfall threshold development are as follows: firstly, the collection of the dataset features. This process includes the landslide data inventory and the determination of the rainfall data source. The catalogue of the data should be accurate or almost accurate, including items such as the date, time, type of landslide, the location and the data source. Other information such as the slope gradient, topography and lithology and details of the landslide are essential for analysis and comparison. The selection of rainfall data is very crucial as it is correlated with the risk of landslide. Therefore, details on the criteria of the rainfall data to be used in the analysis are very important. The method of criteria selection of rainfall data such as a single reference, manual/expert judgement, the reference rain gauge, the nearest rainfall station and the Theisen polygon method should be specified clearly in the study [
3]. The rainfall parameter that governs the relationship of the landslide is determined, whereby the governing parameter is considered to be the causal rain. Next is the identification of the threshold. The outcome of this step is the established relationship between the causal rain and the disaster with the threshold line. This process requires the modeler to draw the separation line and once it is breached, the disaster is likely to occur. Therefore, the drawn line is very subjective; if it is too conservative then the warning could be missed; on the contrary, if it is too low then it could end up resulting in a false alarm. The practical method, as suggested by Nikolopoulos et al. [
24], is based upon visual observation of the modeler to separate the safe zone and unsafe zone at the heuristic lower boundary of the empirical data. Another method is based on the
frequentist method [
26]. Lastly, is the validation of the model. This step is categorized into two parts. Firstly, is the model validation and secondly, is the model performance. In order to reach the second category, the model should be implemented or applied in the early warning system. As per the first category, there is no standard method in place worldwide as yet, but Segoni et al. [
3] suggested that the validation can be made using the same dataset or independent datasets. The warning model is not static but has a dynamic manner.
This entry is adapted from the peer-reviewed paper 10.3390/w14203212