Next Article in Journal
Estimating the Leaf Water Status and Grain Yield of Wheat under Different Irrigation Regimes Using Optimized Two- and Three-Band Hyperspectral Indices and Multivariate Regression Models
Previous Article in Journal
Comparison between Deep Learning and Tree-Based Machine Learning Approaches for Landslide Susceptibility Mapping
Previous Article in Special Issue
Forecasting of Debris Flow Using Machine Learning-Based Adjusted Rainfall Information and RAMMS Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modelling Debris Flow Runout: A Case Study on the Mesilau Watershed, Kundasang, Sabah

1
Razak Faculty of Technology and Informatics, Universiti Teknologi Malaysia (UTM), Kuala Lumpur 54100, Malaysia
2
Disaster Preparedness and Prevention Centre (DPPC), Malaysia-Japan International Institute of Technology (MJIIT), Universiti Teknologi Malaysia (UTM), Kuala Lumpur 54100, Malaysia
3
Department of Civil, Environmental and Applied Systems Engineering, Kansai University, Osaka 564-8680, Japan
4
Faculty of Engineering, Tottori University, Tottori 680-8552, Japan
*
Authors to whom correspondence should be addressed.
Water 2021, 13(19), 2667; https://doi.org/10.3390/w13192667
Submission received: 24 August 2021 / Revised: 16 September 2021 / Accepted: 23 September 2021 / Published: 27 September 2021
(This article belongs to the Special Issue Debris Flows Research: Hazard and Risk Assessments)

Abstract

:
Debris flows are among the fatal geological hazards in Malaysia, with 23 incidents recorded in the last two decades. To date, very few studies have been carried out to understand the debris flow processes, causes, and runouts nationwide. This study simulated the debris flow at the Mesilau watershed of Kundasang Sabah caused by the prolonged rainfall after the 2015 Ranau earthquake. Several interrelated processing platforms, such as ArcGIS, HEC-HMS, and HyperKANAKO, were used to extract the parameters, model the debris flow, and perform a sensitivity analysis to achieve the best-fit debris flow runout. The debris flow travelled at least 18.6 km to the Liwagu Dam. The best-fit runout suggested that the average velocity was 12.5 m/s and the lead time to arrive at the Mesilau village was 4.5 min. This high debris flow velocity was probably due to the high-water content from the watershed baseflow with a discharge rate of 563.8 m3/s. The flow depth and depositional thickness were both lower than 5.0 m. This study could provide crucial inputs for designing an early warning system, improving risk communication, and strengthening the local disaster risk reduction and resilience strategy in a tectonically active area in Malaysia.

1. Introduction

Debris flows are among the disastrous geological hazards that occur in both developed and developing countries. This sediment-related disaster strikes quickly without warning, and results in severe consequences, including human and economic losses [1]. In 2021, a number of deaths due to debris flow alone were reported: debris flow in Uttarakhand, India (February 2021), which killed more than 200 people [2]; debris flow (mudslide) in Atami, Japan (July 2021), which killed at least three people, and 80 people remained missing [3]; and debris flow in Gunung Jerai, Kedah, Malaysia (August 2021), which wiped out at least six people with two victims still missing [4,5]. In Malaysia, approximately 23 debris flow incidents were recorded nationwide between 1995 and 2015 (Table 1). These events were mostly induced by rainfall and very few were triggered by the cascading earthquake. Incidentally, Malaysia’s geographic location near the equator with a tropical climate receives an abundant annual rainfall, i.e., an average of 2400 mm [6]. In Sabah, the vicinity of Mount Kinabalu, the first UNESCO World Heritage Site in Malaysia has recorded two recent debris flows at the Mesilau watershed of Kundasang (the southeast flank of Mount Kinabalu) and the Kedamaian watershed of Kota Belud (the northwest flank of Mount Kinabalu), respectively [7,8].
The most fatal of debris flow in Malaysia occurred in 1996 at Keningau, Sabah, which killed 302 people and caused an economic loss of 458.9 million RM [23,25]. The debris flow was induced by Typhoon Gregg; it remains the deadliest geological disaster in the history of Malaysia. The first earthquake-induced debris flow in Malaysia happened in 2015 at the Mesilau watershed of Kundasang, Sabah, after prolonged post-earthquake rainfall. Although the debris flow recorded no fatality along the channelised river, it caused a significant socioeconomic impact and indirect economic losses [8,29,30].
The concept of debris flow initiation at the Mesilau watershed was adapted from the debris flow after the Gorkha earthquake of Nepal in 2015 [31]. After the earthquake, two forms of hazards, known as the primary and secondary hazards, were initiated [31]. The primary hazard encompasses the direct geological failures from the earthquake (i.e., rock avalanches, rock falls, and landslides), and the secondary hazard refers to the reactivated or remobilised landslides dam from the earthquake, for example, debris flows.
Incidentally, the initiation of debris flow of the Gorkha earthquake could also adequately describe the debris flow after the Ranau earthquake at the Mesilau watershed (Figure 1). The primary hazard was seen at the foot slope of Mount Kinabalu, where the earthquake induced many landslides, forming a temporary landslide dam, while the secondary hazard, known as debris flow, occurred along the channelised Mesilau river (Figure 1, the yellow polygon). Figure 2 shows the evidence of accumulated boulders along the channelised river after the debris flow.
The debris flow in the channelised Mesilau river was mainly induced by prolonged rainfall ten days after the 6.0 Mw Ranau earthquake on 5 June 2015 (Figure 3) with a series of aftershocks. Frequent rains were recorded for seven days with a cumulative rainfall of 66.3 mm and the highest rainfall intensity of 14.2 mm/h on 15 June 2015, which triggered the event. Chronologically, the extensive ground shaking from the main earthquake triggered many shallow and deep-seated landslides within the vicinity of Mount Kinabalu on 5 June 2015 and partially accumulated the debris and sediment as the temporary landslide dam [7]. The subsequent prolonged rainfall further reduced the stability of the landslide dam and reactivated the landslide dam as debris flow. Several debris flows impacted the channelised river and the village area, as reported by [8].
Several studies have been conducted to investigate the earthquake at Kundasang [7], which have focused mainly on the causes and impacts of the earthquake [8], landslide inventory [32,33], landslide susceptibility [34], landslide hazard and risk [35], land cover assessment [36], and business continuity plan [30]. However, little is known about the debris flow travelling along the channelised Mesilau river from the upstream to the downstream areas.
Debris flow modelling has become one of the effective methods to understand the debris flow processes [37]. The modelling of debris flow allows the simulation of past events, prediction of future events, and simulation of an area without historical evidence [38]. Therefore, the simulation results can contribute to the analysis of potential hazard scenarios and suitable planning for debris flow risk reduction strategies [39,40]. Generally, three methods for debris flow modelling are available which include physical modelling, empirical modelling, and dynamic modelling [41]. The differences between these methods are that physical modelling is conducted based on the field observation and supported by the controlled laboratory experiment; whilst the empirical modelling is performed based on the well-documented observation and, usually, it is practical to estimate the travel distance without considering the rheology of debris flow [42]. In addition, dynamic modelling is performed by using the numerical methods applying the momentum and energy conservation law [43]. Recently, a variety of debris flow simulation models exist globally, such as TRENT2D [44], AschFLow [45], HyperKANAKO [46], and RAMMS:Debris Flow [47,48]. Each model is distinguished by its applicability, price, parameters, and algorithms within the model [38].
This study examined the debris flow runout processes via simulation using HyperKANAKO and reconstruction of a past debris flow event. This study encompassed the preparation of the base model, extraction of the rainfall-runoff model, and model calibration for generating the best-fit debris flow runout via several interrelated processing platforms. A series of field investigations and interviews with the local community was also conducted to validate the runout results. The simulated result could provide important insight for improving the management of hazards and risks in this tectonically active area.

2. Materials and Methods

2.1. Study Area

The selected study area is an undulating terrain located at the southeast flank of Mount Kinabalu, Kundasang with altitudes ranging from 578 to 2384 m [7]. Figure 4 shows the Mesilau watershed, ranging from the upstream area in the mountain to the downstream area in the Ranau town (the irregular yellow boundary), while the source area of debris flows is located in the eastern plateau of Mount Kinabalu (the red polygon). The latitudes of the study area range between 06°05′02.0″ and 05°57′34.6″, while the longitudes encompass 116°32′53.3″ and 116°41′01.6″.
There are three settlements within the study area (green polygons), i.e., the Mesilau village, Naradau village, and Ranau town. Unlike the other two settlements, the Mesilau village has been exposed to many geological hazards and associated risks due to its proximity to the foot slope of the mountain, two main rivers, and tourist areas with many attractive infrastructures. These infrastructures include Desa Dairy Farm, Mesilau Golf Club, Strawberry Farm Mesilau, Mesilau Cat’s village, Mesilau Nature Resort, homestays, Maragang hill, and Sososdikon hill [7]. The area is also widely known as one of the primary agriculture production zones in Sabah, due to its climate, topography, and tourism activities in the vicinity of Mount Kinabalu.
Geologically, the study area is mainly dominated by five lithologies, i.e., Serpentinite, the Crocker formation, the Trusmadi formation, Granite, and the Pinousok gravel [49,50,51,52]. The Crocker formation consisted of strongly folded and faulted sandstone, siltstone, red and grey shale, mudstone, and argillite; the Trusmadi formation can be described as strongly folded and faulted grey and dark grey argillite, slate, siltstone and sandstone with volcanic; and the Pinousok gravel is characterized by its poorly unconsolidated gravel up to boulder sizes in a sandy to clayey matrix [49,50,51,52,53]. Typically, these lithologies comprised mostly weathered materials [54]. Figure 5 shows the five lithologies within the watershed. The source area at the foot slope of Mount Kinabalu is comprised of the granitic rock and Serpentinite, while the transportation and deposition area is comprised of Pinousok gravel, Crocker formation, and Trusmadi formation. There are two major prominent and active faults in the study area: the Lobou-Lobou fault, a left-lateral strike faulting N20E and the Mensaban fault trending northwest-southeast [55]. The associated fault during the 2015 Ranau earthquake was due to the normal slip of the Lobou-Lobou fault, with the epicentre located at the highland of Kundasang town and a shallow depth of 10 km beneath Mount Kinabalu [8]. Historically, the Ranau district recorded three earthquakes exceeding 5.0 Mw [7,8], i.e., in 1966 (Mw 5.3), 1991 (Mw 5.2), and 2015 (Mw 6.0), with approximate return periods of every 24 to 25 years [7]. Evidence of weak geological materials and earthquake studies have categorised this area as high risk to geological hazards and the associated risks [54,56,57].

2.2. Datasets

This study employed three main datasets, i.e., the orthophoto, Digital Terrain Model (DTM), and rainfall dataset (Table 2). Additionally, the Light Detection and Ranging (LiDAR) point-cloud dataset was obtained from the Department of Mineral and Geoscience Malaysia (JMG Malaysia); it was acquired a year after the earthquake disaster under the Slope Hazard and Risk Mapping (PBRC) project. The LiDAR acquisitions (i.e., orthophoto and DTM) were conducted through airborne laser scanning (ALS) with the combined systems of Laser Scanner 6800-400, IMU-Ilf Inertial Measurement Unit, Trimble R7 GPS Receiver, DigiCAM H60, DigiCAM Lens Cone HC 3,5/50-11, Clear-Protection-Filter, 77 mm, LiteMapper system mount, LiteMapper Data Recorder 680, and AEROcontrol Sensor Management Unit [58]. The Interferometric Synthetic Aperture Radar (IfSAR) dataset was acquired from the Intermap STAR Technologies®’ proprietary X-Band IfSAR, and the hourly rainfall records were retrieved from the Malaysian Meteorological Department (MMD). The coordinate system of the study area was set as WGS 1984 UTM Zone 50N to be congruent with that of the HyperKANAKO model. In this study, the orthophoto and LiDAR-derived DTM after the debris flow event were only used as a reference to evaluate and validate geomorphological changes in the channel and to extract the channel parameters, respectively. Datasets of IfSAR DTM and hourly rainfall were used to model the past debris flow. The overall flowchart of the study is summarized in Figure 6.

2.3. Base Model Extraction

The extraction of the base model prepared the watershed model for processing the rainfall runoff and modelling the debris flow. The base model was mainly extracted using the ArcGIS software version 10.8 with extensions of ArcHydro and HEC-GeoHMS. These extensions enabled the efficient creation of the base model in sequence following the step-wise procedures of other studies [59,60]. Figure 7 summarises the steps used to generate the base model in this study, starting with the extraction of the base model for the rainfall runoff processes. Steps 1 to 5 were conducted using the ArcHydro extension (the pink-blue rectangle), and Steps 6 to 10 were processed with the HEC-GeoHMS extension (the yellow-blue rectangle). Figure 8 shows the extracted base model; it was then exported as *.hms, and later imported into the Hydrologic Modelling System (HEC-HMS) software for processing the rainfall runoffs.
To model the debris flow, the outline shapefile of the overall watershed area (*.shp) was used to clip the raw IfSAR DTM (*.tiff) in order to obtain the boundary of the watershed area. Figure 9 shows the clipped IfSAR DTM (*.tiff) following the shapefile boundary of the watershed area (yellow). Yellow dots along the channelised river represent the distance of the past debris flow runout for every 100 m. The total debris flow runout was 18,600 m (18.6 km), with 187 dots from the initiation to the deposition area. Then, the watershed area (*.tiff) was imported into the open-source QGIS software to model the debris flow.

2.4. Rainfall Runoff Processes

Rainfall runoffs were processed with HEC-HMS software to obtain the hydrograph that triggered the debris flow event; two parameters were required by the software, i.e., the prepared base model by ArcHydro and HEC-GeoHMS extensions in ArcGIS and the hourly rainfall dataset. Figure 10 shows the workspace of the HEC-HMS software highlighting the imported base model (*.hms). Overall, the imported base model defined the basin properties, i.e., basin, sub-basin, junction, outlet, etc.
In the HEC-HMS basin model, the loss method, transform method, and baseflow method were set as the “initial and constant”, “Snyder Unit Hydrograph”, and “constant monthly”, respectively. The 2015 Ranau earthquake altered the hydrological processes within the watershed area. Water infiltrating the soil and fracturing the bedrock drained rapidly, causing large rapid flow [61]. Therefore, the initial and continuous losses were set from 1 to 5 mm based on the impervious properties of geological and land use of the watershed area.
Meanwhile, the baseflow is derived from bedrock water storage near-surface valley soils and riparian zones. Unlike the observable surface flow, the baseflow was estimated using various methods, such as tracers and baseflow separation with streamflow hydrograph. Due to the insufficient streamflow data for validation, the baseflow value was determined on a trial-and-error basis. A baseflow value of 950 m3/s for Kenyir Lake [62], the largest reservoir in Peninsular Malaysia, was used as a reference. Overall, the value for the baseflow in the Mesilau watershed was calibrated with at least ten trial-and-errors, ranging from 100 to 1000 m3/s. The Snyder Unit Hydrograph basin lag was calculated using Equation (1) below [62]:
t p = C C t ( L L c ) 0.3
where C t is the basin coefficient (0.8), L is the length of the mainstream from the farthest outlet to the divide (4.65 km), L c is the length along the mainstream from the outlet to the nearest point in the watershed centroid (2.33 km), and C is the conversion constant (0.75).
The time-series data were set by importing the hourly rainfall dataset retrieved from MMD. The precipitations were assumed constant in all areas in the watershed regardless of the distance from the rainfall station.
Lastly, the control specification was set by adjusting the start/end of the date and time, and the time interval output. Following the HyperKANAKO model format, the time interval output in HEC-HMS software was set for every 60 s (1 min). The control specification allowed the software to simulate the duration of the hydrograph. Later, the ten hydrographs generated from the simulation were imported into the HyperKANAKO model to obtain the best-fit debris flow runout.

2.5. Debris Flow Modelling

The HyperKANAKO model is a debris flow simulation system that uses the debris-flow simulator KANAKO 2D developed by [46], with the flexibility of preparing datasets, user-friendly models, and an enhanced interface. Therefore, this study used this HyperKANAKO model to examine the debris flow processes and runout.
Equations for modelling included the momentum in the x-axis (Equation (2)) and y-axis (Equation (3)), the continuation of debris-flow volumes (Equation (4)), the continuation of debris flow particles (Equation (5)), and changes in the bed surface elevation (Equation (6)) [46,63]. Additionally, this model incorporated the mutual influences of the one-dimensional (1D) simulation area (i.e., gullies), and two-dimensional (2D) simulation area (i.e., alluvial fan) [64].
Equation (2) provides the momentum in the direction of the x-axis:
u t + u u x + v u y = g sin θ w x τ x ρ h
where u is the flow velocity in the direction of the x-axis, θ w x is the flow-surface gradient in the x-axis, τ x is the riverbed shearing stresses in the x-axis, g is the gravity acceleration, ρ is the mass density of fluid phase, and h is the flow depth.
Equation (3) describes the momentum in the direction of the y-axis:
v t + u v x + v v y = g sin θ w x τ y ρ h
where v is the flow velocity in the direction of the y-axis, θ w x is the flow-surface gradients in the y-axis, and τ y is the riverbed shearing stresses in the y-axis.
The continuation of debris-flow volume was computed using Equation (4) as follows:
h t + u h x + v h y = S T
where S T is the sediment erosion or the deposition velocity. The continuation of debris flow particles was calculated by Equation (5):
C h t + u C u h x + v C v h y = S T C *
where C is the sediment concentration by volume in the debris flow and C * is the sediment concentration by volume in the moveable bed layer.
The determination of changes in the bed elevation was based on Equation (6) as follows:
z t + S T = 0
where z is the bed elevation.
Meanwhile, the model setup, began with the streamline within the initiation and deposition area. The initiation line of this study was located in the landslide dam, whereas the deposition line was within the Liwagu Dam of the Ranau town. A green line and a yellow rectangle were displayed on the workspace, representing the expected runout and the deposition area, respectively. Subsequently, the deposition area was adjusted by changing the mesh sizes of the x- and y-axes until the debris flow was projected to be deposited in the Liwagu Dam. Figure 11 shows the model setup within the study area. The type of dam and its heights were set within the initiation and deposition area, respectively. The heights of the post-landslide dam and the Liwagu Dam were measured using the laser-range finder during a field investigation on 31 August 2019. The landslide dam within the initiation area was set to the height of 30.0 m (Figure 12a), while the closed dam within the downstream area was adjusted to a height of 5.0 m (Figure 12b).
Additionally, the modelling parameters were calibrated through the open-source LibreOffice 5 software with three primary tabs, i.e., the “DEM, Dams, and Observation”, “riverbed”, and “hydrograph”. Table 3 shows the calibrated parameters for the “DEM, Dams, and Observation”; these parameters were calibrated by reviewing the available historical records of the study area with modifications based on global events, such as Japan [63,65] and Indonesia [66]. Additionally, several parameters, such as mass density of bed materials, mass density of fluids, coefficient of erosion and deposition rate, and acceleration of gravity of the study area were set as the default. For the “riverbed”, the average channel width was adjusted to 5.0 m upon crosschecking the orthophoto and DEM before and after the debris flow event by changing the value in the riverbed. Table cross-sections were prepared for every 100 m (the yellow dots) to evaluate the channel width. Finally, the hydrograph generated from the HEC-HMS was imported into the relevant tab by manually inserting these values for every timestamp.
Additionally, a sensitivity analysis was conducted to attain the best-fit debris flow runout for the past event with ten simulated scenarios by importing ten hydrographs acquired through the HEC-HMS. Indicators for evaluating the runouts included the flow depth, depositional thickness, and runout distance. The flow depth was evaluated along the Mesilau Nature Resort bridge (upstream area), whereas the depositional thickness and runout distance were observed within the Liwagu Dam (downstream area). The best-fit debris flow runout was presented in a spatial map.

3. Results

3.1. Best-Fit Debris Flow Runout

Table 4 summarizes the simulation results for the ten different scenarios. Overall, the runout distance increased along with the hydrographs. For example, Scenarios 1 to 4 produced short runouts, whereas Scenarios 6 to 10 yielded overestimated runouts. The best-fit runout happened between Scenarios 5 and 6, and both scenarios were nearly similar to the past event. Therefore, the recalibration was conducted for the parameters between Scenarios 5 and 6 to produce the best-fit debris flow runout labelled as “Scenario BF”. The result shows that the best-fit debris flow runout was obtained as a hydrograph of 563.8 m3/s and a baseflow of 550 m3/s.

3.2. Flow Depth and Depositional Thickness of Best-Fit Runout

Figure 13 and Figure 14 show the flow depth and depositional thickness of the best-fit debris flow runout, respectively. As highlighted, a multicoloured line was visualised along the channel, and each colour represented different ranges of flow depths and depositional thickness. These maps also show five different zones, highlighting significant areas, such as Zone A (initiation area), Zone B (Mesilau village), Zone C (transportation area), Zone D (Naradau village), and Zone E (Liwagu Dam). Table 5 and Table 6 summarise the flow depth and depositional thickness, respectively, for each zone. For a detailed evaluation, the flow depth was analysed in Zone B during the transportation process, whereas the depositional thickness was focused on Zone E after the simulation had ended. These two areas were selected based on the field data inputs and evidence collected by interviewing the eyewitnesses of the debris flow event and obtaining input from local communities.
The flow depth within Zone B was most frequently observed and marked by the green line with flow depths of 1.0 to 5.0 m; the least visualised zone was represented by the yellow line with a flow depth ranging from 5.1 to 10.0 m (Figure 13). Additionally, the flow depth did not overflow from the channel bank, indicating that the simulation result was congruent with the past event. Interviews with the local community further supported flow depth by justifying the highest depth within the channel banks during the event, especially along the village area (Figure 13).
The depositional thickness of Zone E, mostly highlighted with the blue to green line, ranged from 0.1 to 5.0 m (Figure 14). The widespread of blue line within the Liwagu Dam was indicative of the overflowing debris from the channel bank (Figure 14).

3.3. Estimated Velocity and Lead Time

The best-fit runout was also used to estimate the debris-flow velocity and the lead time to arrive at the nearest villages, i.e., Mesilau village, Naradau village, and Ranau town. The estimated channel distance from the landslide dam (the initiation area) to the nearest villages were 3.9 km, 11.1 km, and 18.6 km, respectively. The estimated average debris flow velocity derived from the simulation result was 12.5 m/s. The lead time to arrive at each village was 4.5 min, 14 min, and 24 min, respectively. The analysis showed that the local community of Mesilau village would have the shortest time of 4.5 min to evacuate as compared with the Naradau village and Ranau town. This finding would be crucial for designing an early warning system (EWS) to disseminate risk information and the evacuation route.

4. Discussion

This study investigated the debris flow runout triggered by the prolonged rainfall after the 2015 Ranau earthquake in the Mesilau watershed, Kundasang, Sabah. To date, this is the first study conducted within East Malaysia to simulate the past debris flow in the vicinity of Mount Kinabalu. This study provides better insights into debris flow characteristics, including the estimated velocity, lead time, flow depth, and depositional thickness, which would allow the development of debris flow prevention and mitigation measures. Thus far, no EWS has been installed in the area, which increases the likelihood of future debris flows with mild seismic activity and rainfall.
Overall, the debris flow in this area was different from other reported events, such as Genting Sempah, Gunung Pulai, and Pos Dipang [9,14,19] as the source materials (i.e., sediment and debris) of this area were induced by the earthquake and later remobilised as a debris flow following a seven-day prolonged rainfall (66.3 mm) and least amount of rainfall intensity (14.2 mm/h). Consequently, a long debris flow runout of 18.6 km was recorded flowing from the foot slope of Mount Kinabalu until the Liwagu Dam of Ranau town. This debris flow permanently changed the river morphology, disrupting critical infrastructures (e.g., bridges, roads, and dams), degrading the water quality, and affecting socioeconomic activities.
The HyperKANAKO model used in this study was sensitive towards the hydrographs, i.e., the hydrographs increased along with flow depths, depositional thickness, and runout distance. The best-fit debris flow runout, recalibrated between Scenarios 5 and 6, showed that the past debris flow had a discharge of 563.8 m3/s, mostly derived from the baseflow of 550 m3/s, indicating that the least amount of rainfall was enough to trigger the debris flow after the earthquake in the Mesilau watershed. The baseflow value of 550 m3/s for the Mesilau watershed was reasonable as it was less than the baseflow for water bodies [62]. Additionally, the estimated debris flow velocity and the lead time to arrive at the nearest Mesilau village were 12.5 m/s and 4.5 min, respectively. Due to its high debris-flow velocity, this study recommends a community evacuation of 4.5 min before the debris flow arrives in the village. Since the evaluated lead time of 4.5 min is too short for the community to complete the evacuation, there is a need to develop an impact-based multi-hazard EWS for disseminating the warning message as early as possible. Additionally, an early detection system on natural or landslide dams in the highland area is suggested. Once the emergence and growth of landslide dams are detected, a warning should be issued for the community to evacuate from high-risk areas before the commencement of debris flow. The EWS is proposed for the watershed area following the Malaysia’s initiative to achieve the global targets of the Sendai Framework for Disaster Risk Reduction 2015–2030.
The simulation results of this study show that the flow depth within Zone B (Mesilau village) was below 10.0 m, and this result was congruent with past evidence. Then, the flow depth gradually decreased as it travelled to Zones C, D, and E. Meanwhile, the depositional thickness ranged from 0.1 to 5.0 m within Zone E (Liwagu Dam), probably due to the accumulated debris within the closed dam.
These results (i.e., reconstruction of past events, flow depth, depositional thickness, debris flow velocity, and lead time) might serve as baseline data to create greater awareness among local stakeholders and vulnerable communities towards understanding the past debris flow event for preventing future risks. This study is timely and significant as the earthquake return period is approximately 25 years in the Ranau district, which might trigger future debris flows in the area [57]. Thus, further studies are needed to include more parameters for modelling debris flow and to design suitable structural and nonstructural mitigation measures in the localised area. An integrated framework is critically needed for responding to sediment-related disasters, as the area is highly vulnerable to future debris flow given the active tectonic activities, hanging rocks in the mountains due to geodynamic activity, and high exposure to typhoons and their hydrometeorological impact.

5. Conclusions

This study modelled the past debris flow in the Mesilau watershed by applying several interrelated processing platforms, including ArcGIS (ArcHydro and HEC-GeoHMS), HEC-HMS, QGIS, and HyperKANAKO. On the basis of the observation data and the best-fit simulation result, it can be concluded that:
(1)
The debris flow happened on 15 June 2015, i.e., ten days after the 2015 Ranau earthquake (Mw 6.0) with a seven-day cumulative rainfall of 66.3 mm. The maximum rainfall intensity was 14.2 mm/h; it breached the landslide dam and initiated the debris flow. The early identification showed that the least amount of rainfall was sufficient to trigger the debris flow after the earthquake in the Mesilau watershed.
(2)
According to the best-fit simulation, the debris flow velocity was estimated to be 12.5 m/s and the lead time to arrive at the nearest Mesilau village was 4.5 min, representing the required evacuation time by the community to minimise the debris flow impacts and prevent human losses.
(3)
Additionally, the baseflow during the past event was 550 m3/s, yielding a discharge of 563.8 m3/s. According to the reference value of the Kenyir lake [62], the baseflow for the Mesilau watershed should not exceed the baseflow of water bodies, i.e., 950 m3/s.
(4)
The approximate runout distance flowing from the landslide dam to the Liwagu Dam of the Ranau town was 18.6 km, due to a large amount of sediment supply (accumulated debris) generated during the 2015 Ranau earthquake.
(5)
As the debris began flowing, the flow depth gradually decreased from 15.0 m at the initiation area to 5.0 m at the deposition area. In Zone B (the Mesilau village), the flow depth ranged from 1.0 to 10.0 m with no overflowing debris identified. In Zone E (the Liwagu Dam), the depositional thickness was less than 5.0 m; the closed dam blocked the accumulated debris along the channel.
There was some degree of uncertainty while simulating the event since assumptions and empirical laws were used based on the inputs that the HyperKANAKO model requires. Although debris flow runout models have been used with regularity in the past to re-establish past events via the calibration of input parameters, there are still some limitations in the physical description. Therefore, future research is needed for both physical and modelling studies, especially the complex processes and mechanism of debris flow in this earthquake-prone area in order to quantify and to reduce uncertainty.
Predictive assessments may facilitate improved management of future debris flow, and therefore reduce the hazard. This study supports the commitment of the Malaysia Government to achieve the global targets of Sendai Framework for Disaster Risk Reduction 2015–2030 particularly the Global Target G-5 on the accessible, understandable, usable, and relevant disaster risk information, which can later be used to design and develop the multi-hazard EWS. Nonstructural mitigation measures such as early warning systems might allow the local government to estimate the percentage of the population exposed to or at risk from disasters and protect them through preventive evacuation following early warning (Global Targets G1, G3, G4, and G6), and to develop local disaster risk reduction and resilience strategies (Global Target E2) [67].

Author Contributions

Conceptualization, M.I.R. and F.C.R.; writing original draft preparation, M.I.R.; methodology, M.I.R. and F.C.R.; formal analysis, M.I.R.; validation, M.I.R., F.C.R., K.A.R., S.A., S.A.K., A.N.A., A.M., T.T. and Y.O.; writing review and editing, F.C.R., K.A.R., S.A., S.A.K., A.N.A., A.M., T.T. and Y.O.; visualization, M.I.R.; supervision, F.C.R., K.A.R., S.A. and A.N.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Japan International Cooperation Agency (JICA).

Acknowledgments

The authors would like to thank JICA for initiating the debris flow project and the financial support to make the study possible. The authors would also like to thank the international project members from Japan who contributed to the project. Our highest appreciations go to Kana Nakatani from Kyoto University, Japan for the HyperKANAKO dongle and continuous support and advice throughout the simulation. Our great thanks also go to Muhammad Fitri Shahrim and Nur Afiqah Mohd Kamal for their assistance, and to the local stakeholder and community, who provided information on the past debris flow event and validated the best-fit debris flow runout.

Conflicts of Interest

The authors declare no conflict of interest.

References and Note

  1. Santi, P.M.; Hewitt, K.; VanDine, D.F.; Cruz, E.B. Debris-flow impact, vulnerability, and response. Nat. Hazards 2011, 56, 371–402. [Google Scholar] [CrossRef]
  2. Shugar, D.H.; Jacquemart, M.; Shean, D.; Bhushan, S.; Upadhyay, K.; Sattar, A.; Schwanghart, W.; McBride, S.; de Vries, M.V.W.; Mergili, M.; et al. A massive rock and ice avalanche caused the 2021 disaster at Chamoli, Indian Himalaya. Science 2021, 373, 300–306. [Google Scholar] [CrossRef] [PubMed]
  3. Petley, D. The Atami Mudslides in Japan. Available online: https://blogs.agu.org/landslideblog/2021/07/05/atami-mudslides/ (accessed on 2 August 2021).
  4. Radhi, R. 47 Lokasi Tanah Runtuh Di Gunung Jerai Dikenal Pasti. Available online: https://www.getaran.my/artikel/semasa/10261/47-lokasi-tanah-runtuh-di-gunung-jerai-dikenal-pasti (accessed on 10 September 2021).
  5. Jabatan Mineral dan Geosains. The Early Assessment of Debris Flow in Gunung Jerai, Kedah. Kedah. Virtual Interview by Muhammad Iylia Rosli. 13 September 2021.
  6. Mukhlisin, M.; Matlan, S.J.; Ahlan, M.J.; Taha, M.R. Analysis of Rainfall Effect to Slope Stability in Ulu Klang, Malaysia. J. Teknol. 2015, 72, 15–21. [Google Scholar] [CrossRef] [Green Version]
  7. Rosli, M.I.; Kamal, N.A.M.; Razak, K.A. Assessing Earthquake-induced Debris Flow Risk in the first UNESCO World Heritage in Malaysia. Remote Sens. Appl. Soc. Environ. 2021, 23, 100550. [Google Scholar] [CrossRef]
  8. Tongkul, F. The 2015 Ranau Earthquake: Cause and Impact. Sabah Soc. J. 2016, 32, 1–28. [Google Scholar]
  9. Chow, W.; Jadid, M.S.M.; Yaacob, S. Geological and geomorphological investigations of debris flow at Genting Sempah, Selangor. In Proceedings of the GSM, Forum on Geohazards: Landslide & Subsidence, Kuala Lumpur, Malaysia, 22 October 1996; pp. 5.1–5.15. [Google Scholar]
  10. Ghazali, M.A.; Rafek, A.G.; Desa, K.M.; Jamaluddin, S. Effectiveness of Geoelectrical Resistivity Surveys for the Detection of a Debris Flow Causative Water Conducting Zone at KM 9, Gap-Fraser’s Hill Road (FT 148), Fraser’s Hill, Pahang, Malaysia. J. Geol. Res. 2013, 2013, 1–11. [Google Scholar] [CrossRef] [Green Version]
  11. Hashim, K.; Among, H.L. Geological investigation on Ruan Changkul landslide. Bull. Geol. Soc. Malays. 2003, 46, 125–132. [Google Scholar] [CrossRef]
  12. Huat, L.T.; Ali, F.; Osman, A.R.; Rahman, N.A. Web based real time monitoring system along North-South Expressway, Malaysia. Electron. J. Geotech. Eng. 2012, 17, 623–632. [Google Scholar]
  13. Jaafar, K. Heavy rainfall intensity triggering the landslide at Sungai Ruil Cameron Highland Malaysia-A Case study. In Proceedings of the EGU General Assembly Conference Abstracts, Vienna, Austria, 22–27 April 2012; p. 13728. [Google Scholar]
  14. Jamaludin, S.; Abdullah, C.H.; Kasim, N. Rainfall Intensity and Duration for Debris Flow Triggering in Peninsular Malaysia. In Landslide Science for a Safer Geoenvironment; Springer Science and Business Media LLC: Berlin, Germany, 2014; pp. 167–172. [Google Scholar]
  15. Jamaludin, S.; Hussein, A.N. Landslide hazard and risk assessment: The Malaysian experience. Notes 2006, 1–10. Available online: http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.484.88&rep=rep1&type=pdf (accessed on 21 September 2021).
  16. Joe, E.J.; Tongkul, F.; Roslee, R. Engineering Properties of Debris Flow Material at Bundu Tuhan, Ranau, Sabah, Malaysia. Pak. J. Geol. 2018, 2, 22–26. [Google Scholar] [CrossRef]
  17. Kasim, N.; Abu Taib, K.; Mukhlisin, M. Comparison of Debris Flow Simulation Model with Field Event in Kuala Kabu Baru, Malaysia. In Proceedings of the 9th International Conference of Geotechnical & Transportation Engineering (GEOTROPIKA) and 1st International Conference on Construction and Building Engineering (ICONBUILD)—GEOCON2013, At: Persada Johor International Convention Centre, Johor Bahru, Malaysia, 12–13 November 2013; pp. 1–11. [Google Scholar]
  18. Kasim, N.; Taib, K.A.; Mukhlisin, M.; Kasa, A. Triggering mechanism and characteristic of debris flow in Peninsular Malaysia. Am. J. Eng. Res. 2016, 112–119. Available online: http://www.ajer.org/papers/v5(04)/M050401120119.pdf (accessed on 21 September 2021).
  19. Komoo, I. Pembinaan Ambil Kira Aliran Debris Mampu Kurang Risiko. Available online: https://origin.bharian.com.my/taxonomy/term/61/2015/11/97519/pembinaan-ambil-kira-aliran-debris-mampu-kurang-risiko (accessed on 2 July 2019).
  20. Kuraoka, S.; Abdullah, C.H.; Jamaluddin, S.; Kasim, N.; Nazri, M.; Sze, L.T. Study on the Mechanisms of Debris Flows that Damaged Flexible Barriers in the Channel of Fraser Hill, Malaysia. In Proceedings of the Civil Engineering Conference in the Asian Region (CECAR7), Hawaii, HI, USA, 30 August–2 September 2016. [Google Scholar]
  21. Lay, U.S.; Pradhan, B. Identification of Debris Flow Initiation Zones Using Topographic Model and Airborne Laser Scanning Data. In Proceedings of the EECE 2020, St.Petersburg, Russia, 19–20 November 2020. [Google Scholar]
  22. Loh, A.I.; Leen, C.L. Heavy rains cause mud flood in Camerons, one dead. The Star Online 2014. Available online: https://www.thestar.com.my/news/nation/2014/11/06/mud-flood-camerons/ (accessed on 3 July 2019).
  23. Public Work Department. Malaysia: Sectoral Report-National Slope Master Plan 2009–2023; Public Work Department (PWD): Kuala Lumpur, Malaysia, 2009.
  24. Rahman, H.A. An Overview of Environmental Disaster in Malaysia and Prepa-redness Strategies. Iran. J. Public Health 2014, 43, 17–24. [Google Scholar]
  25. Singh, H.; Huat, B.B.; Jamaludin, S. Slope assessment systems: A review and evaluation of current techniques used for cut slopes in the mountainous terrain of West Malaysia. Electron. J. Geotech. Eng. 2008, 13, 1–24. [Google Scholar]
  26. Sum, C.W.; Mohamad, Z. Debris Slide at Kampung Sg. Chinchin, Gombak, Selangor; Geological Society of Malaysia: Kuala Lumpur, Malaysia, 2003; Volume 13, pp. 1–24. [Google Scholar]
  27. Tan, B.K.; Ting, W. Some Case Studies on Debris Flow in Peninsular Malaysia. In Geotechnical Engineering for Disaster Mitigation and Rehabilitation; Springer Science and Business Media LLC: Berlin, Germany, 2009; pp. 231–235. [Google Scholar]
  28. Japan International Cooperation Agency. Region, Natural Disaster Risk Assessment and Area Business Continuity Plan Formulation for Industrial Agglomerated Areas in the ASEAN; Japan International Cooperation Agency: Tokyo, Japan, 2015.
  29. Kamarudin, K.H. Socio-Economic Impacts of Natural Disasters In The Rural Region Of Kundasang, Sabah: A Community Livelihood Analysis. In Proceedings of the RRPG 7th International Conference and Field Study in Malaysia 2016, Johor Bahru, Malaysia, 15–17 August 2016; pp. 492–501. [Google Scholar]
  30. Lehan, N.; Razak, K.; Kamarudin, K. Business continuity and resiliency planning in disaster prone area of Sabah, Malaysia. Disaster Adv. 2020, 13, 25–32. [Google Scholar]
  31. Rosser, N.; Densmore, A.; Oven, K. Landslides following the 2015 Gorkha Earthquake: Monsoon 2016. Available online: http://ewf.nerc.ac.uk/2016/06/15/landslides-following-2015-gorkha-earthquake-monsoon-2016/ (accessed on 2 July 2018).
  32. Bibi, T.; Razak, K.A.; Rahman, A.A.; Latif, A. Spatio Temporal Detection and Virtual Mapping of Landslide Using High-Resolution Airborne Laser Altimetry (Lidar) in Densely Vegetated Areas of Tropics. ISPRS-Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2017, XLII-4/W5, 21–30. [Google Scholar] [CrossRef] [Green Version]
  33. Yusoff, H.H.M.; Razak, K.A.; Yuen, F.; Harun, A.; Talib, J.; Mohamad, Z.; Ramli, Z.; Razab, R.A. Mapping of post-event earthquake induced landslides in Sg. Mesilou using LiDAR. In Proceedings of the IOP Conference Series: Earth and Environmental Science; IOP Publishing: Bristol, UK, 2016; Volume 37, p. 12068. [Google Scholar]
  34. Asmadi, M.A. Landslide Susceptibility Mapping using Remotely Sensed Data in Kundasang, Sabah. Master Thesis, Universiti Teknologi Malaysia, Johor Bahru, Malaysia, 2018. [Google Scholar]
  35. Jabatan Mineral dan Geosains. Slope Hazard and Risk Mapping (Disaster Risk Reduction Initiative); Jabatan Mineral dan Geosains (JMG): Putrajaya, Malaysia, 2016.
  36. Kamal, N.A.M.; Razak, K.A.; Rambat, S. Land Use/Land Cover Assessment in a Seismically ACTIVE Region in Kundasang, Sabah. ISPRS-Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2019, XLII-4/W16, 433–440. [Google Scholar] [CrossRef] [Green Version]
  37. Turnbull, B.; Bowman, E.T.; McElwaine, J.N. Debris flows: Experiments and modelling. Comptes Rendus Phys. 2015, 16, 86–96. [Google Scholar] [CrossRef]
  38. Quan, L. Dynamic Numerical Run-Out Modeling for Quantitative Landslide Risk Assessment. Ph.D. Thesis, University of Twente, Enschede, The Netherlands, 2012. [Google Scholar]
  39. Fan, R.; Zhang, L.; Wang, H.; Fan, X. Evolution of debris flow activities in Gaojiagou Ravine during 2008–2016 after the Wenchuan earthquake. Eng. Geol. 2018, 235, 1–10. [Google Scholar] [CrossRef]
  40. Horton, A.J.; Hales, T.C.; Ouyang, C.; Fan, X. Identifying post-earthquake debris flow hazard using Massflow. Eng. Geol. 2019, 258, 105134. [Google Scholar] [CrossRef]
  41. Chen, H.; Lee, C.F. Numerical simulation of debris flows. Can. Geotech. J. 2000, 37, 146–160. [Google Scholar] [CrossRef]
  42. Luna, B.Q.; Blahut, J.; Camera, C.A.S.; Van Westen, C.; Apuani, T.; Jetten, V.; Sterlacchini, S. Physically based dynamic run-out modelling for quantitative debris flow risk assessment: A case study in Tresenda, northern Italy. Environ. Earth Sci. 2014, 72, 645–661. [Google Scholar] [CrossRef]
  43. Hussin, H.Y. Probabilistic Run-Out Modeling of a Debris Flow in Barcelonnette, France; University of Twente: Enschede, The Netherlands, 2011. [Google Scholar]
  44. Zugliani, D.; Rosatti, G. TRENT2D❄: An accurate numerical approach to the simulation of two-dimensional dense snow avalanches in global coordinate systems. Cold Reg. Sci. Technol. 2021, 190, 103343. [Google Scholar] [CrossRef]
  45. Luna, B.Q.; Blahut, J.; Van Asch, T.; Van Westen, C.; Kappes, M. ASCHFLOW-A dynamic landslide run-out model for medium scale hazard analysis. Geoenvironmental Disasters 2016, 3, 2379. [Google Scholar] [CrossRef] [Green Version]
  46. Nakatani, K.; Iwanami, E.; Horiuchi, S.; Satofuka, Y.; Mizuyama, T. Development of “Hyper KANAKO”, a debris flow simulation system based on Laser Profiler data. In Proceedings of the 12th Congress INTERPRAEVENT, Grenoble, France, 23–26 April 2012; pp. 269–280. [Google Scholar]
  47. Christen, M.; Kowalski, J.; Bartelt, P. RAMMS: Numerical simulation of dense snow avalanches in three-dimensional terrain. Cold Reg. Sci. Technol. 2010, 63, 1–14. [Google Scholar] [CrossRef] [Green Version]
  48. Christen, M.; Bartelt, P.; Kowalski, J. Back calculation of the In den Arelen avalanche with RAMMS: Interpretation of model results. Ann. Glaciol. 2010, 51, 161–168. [Google Scholar] [CrossRef] [Green Version]
  49. Cottam, M.A.; Hall, R.; Sperber, C.; Kohn, B.P.; Forster, M.A.; Batt, G.E. Neogene rock uplift and erosion in northern Borneo: Evidence from the Kinabalu granite, Mount Kinabalu. J. Geol. Soc. 2013, 170, 805–816. [Google Scholar] [CrossRef]
  50. Hall, R.; Cottam, M.; Suggate, S.; Tongkul, F.; Sperber, C.; Batt, G. The Geology of Mount Kinabalu; Sabah Parks: Kota Kinabalu, Malaysia, 2008; Volume 13, pp. 17–20. [Google Scholar]
  51. Jabatan Mineral dan Geosains. Geological Terrain Mapping of the Ranau Area, Sabah; No. Laporan: JMG. Sabah (GBN) 01/2010; Ninth Malaysian Plan: Projek Geobencana Negeri Sabah; Jabatan Mineral dan Geosains: Putrajaya, Malaysia, 2010.
  52. Jacobson, G. Gunung Kinabalu area, Sabah, Malaysia. Geol. Surv. Malays. Rep. 1970, 8, 1–111. [Google Scholar]
  53. Roslee, R.; Tongkul, F. Engineering Geological Mapping on Slope Design in the Mountainous Area of Sabah Western, Malaysia. Pak. J. Geol. 2018, 2, 1–10. [Google Scholar] [CrossRef]
  54. Tjia, H.D. Kundasang (Sabah) at the intersection of regional fault zones of Quaternary age. Bull. Geol. Soc. Malays. 2007, 53, 59–66. [Google Scholar] [CrossRef] [Green Version]
  55. Sharir, K.; Roslee, R.; Lee, K.E.; Simon, N. Landslide Factors and Susceptibility Mapping on Natural and Artificial Slopes in Kundasang, Sabah. Sains Malays. 2017, 46, 1531–1540. [Google Scholar] [CrossRef]
  56. Tongkul, F. Active tectonics in Sabah–seismicity and active faults. Bull. Geol. Soc. Malays. 2017, 64, 27–36. [Google Scholar] [CrossRef] [Green Version]
  57. Jabatan, G. Geological Terrain Mapping of the Kundasang Area, Sabah; Jabatan Mineral dan Geosains (JMG): Putrajaya, Malaysia, 2003.
  58. Bumitouch PLMC. LiDAR LiteMapper 6800-400. Available online: http://www.bumitouchplmc.com/lidar-litemapper-6800-400/ (accessed on 2 July 2019).
  59. Merwade, V. Terrain Processing and HMS-Model Development Using GeoHMS. Lafayette Indiana 2012. Available online: https://web.ics.purdue.edu/~vmerwade/education/terrain_processing.pdf.https://web.ics.purdue.edu/~vmerwade/education/terrain_processing.pdf (accessed on 21 September 2021).
  60. Merwade, V. Watershed and Stream Network Delineation Using ArcHydro Tools; 2012. Available online: https://web.ics.purdue.edu/~vmerwade/education/hechms_sma.pdf (accessed on 21 September 2021).
  61. Salve, R.; Ghezzehei, T.A.; Jones, R. Infiltration into fractured bedrock. Water Resour. Res. 2008, 44, 44. [Google Scholar] [CrossRef] [Green Version]
  62. Ros, F.; Sidek, L.; Ibrahim, N.; Razad, A.A. Probable maximum flood (PMF) for the Kenyir Catchment, Malaysia. In Proceedings of the International Conference on Construction and Building Technology 2008, Kuala Lumpur, Malaysia, 16–20 June 2008; pp. 325–334. [Google Scholar]
  63. Nakatani, K.; Hayami, S.; Satofuka, Y.; Mizuyama, T. Case study of debris flow disaster scenario caused by torrential rain on Kiyomizu-dera, Kyoto, Japan-using Hyper KANAKO system. J. Mt. Sci. 2016, 13, 193–202. [Google Scholar] [CrossRef]
  64. Wada, T.; Satofuka, Y.; Mizuyama, T. Integration of 1-and 2-dimensional models for debris flow simulation. J. Jpn. Soc. Eros. Control. Eng. 2008, 61, 36–40. [Google Scholar] [CrossRef]
  65. Yanagisaki, G.; Aono, M.; Takenaka, H.; Tamamura, M.; Nakatani, K.; Iwanami, E.; Horiuchi, S.; Satofuka, Y.; Mizuyama, T. Debris Flow Simulation by Applying the Hyper KANAKO System for Water and Sediment Runoff from Overtopping Erosion of a Landslide Dam. Int. J. Eros. Control. Eng. 2016, 9, 43–57. [Google Scholar] [CrossRef] [Green Version]
  66. Syarifuddin, M.; Oishi, S.; Legono, D. Lahar Flow Simulation in Merapi Volcanic Area by HyperKANAKO Model. J. Jpn. Soc. Civ. Eng. Ser. B1 Hydraul. Eng. 2016, 72, I_865–I_870. [Google Scholar] [CrossRef] [Green Version]
  67. United Nations for Disaster Risk Reduction. Sendai Framework for Disaster Risk Reduction. Available online: https://www.preventionweb.net/sendai-framework/sendai-framework-monitor/indicators (accessed on 2 July 2019).
Figure 1. The debris flow (secondary hazard) after the Ranau earthquake (primary hazard) of 2015 at the Mesilau watershed of Kundasang, Sabah.
Figure 1. The debris flow (secondary hazard) after the Ranau earthquake (primary hazard) of 2015 at the Mesilau watershed of Kundasang, Sabah.
Water 13 02667 g001
Figure 2. The evidence of boulders accumulated in the channelised Mesilau river after the debris flow event.
Figure 2. The evidence of boulders accumulated in the channelised Mesilau river after the debris flow event.
Water 13 02667 g002
Figure 3. The recorded total rainfall (daily) and cumulative rainfall between earthquake and debris flow events, retrieved from Malaysia Meteorological Department (MMD).
Figure 3. The recorded total rainfall (daily) and cumulative rainfall between earthquake and debris flow events, retrieved from Malaysia Meteorological Department (MMD).
Water 13 02667 g003
Figure 4. The overall debris flow of the study area at the Mesilau watershed, Kundasang, Sabah.
Figure 4. The overall debris flow of the study area at the Mesilau watershed, Kundasang, Sabah.
Water 13 02667 g004
Figure 5. Illustration of the five lithologies within the watershed area synthesised from [8,49,50,51,52,53,55,57].
Figure 5. Illustration of the five lithologies within the watershed area synthesised from [8,49,50,51,52,53,55,57].
Water 13 02667 g005
Figure 6. The overall flowchart of the study.
Figure 6. The overall flowchart of the study.
Water 13 02667 g006
Figure 7. The summarised steps for generating the base model in sequences. Steps 1 to 5 were performed using ArcHydro (pink-blue rectangle), while Steps 6 to 10 were conducted with HEC-GeoHMS (yellow-blue rectangle).
Figure 7. The summarised steps for generating the base model in sequences. Steps 1 to 5 were performed using ArcHydro (pink-blue rectangle), while Steps 6 to 10 were conducted with HEC-GeoHMS (yellow-blue rectangle).
Water 13 02667 g007
Figure 8. The base model extracted using the ArcHydro and Hec-GeoHMS extensions.
Figure 8. The base model extracted using the ArcHydro and Hec-GeoHMS extensions.
Water 13 02667 g008
Figure 9. The extracted DTM of the watershed area covering from the upstream to the downstream area, with yellow dots along the channel representing the distance of debris flow runout for every 100 m.
Figure 9. The extracted DTM of the watershed area covering from the upstream to the downstream area, with yellow dots along the channel representing the distance of debris flow runout for every 100 m.
Water 13 02667 g009
Figure 10. The workspace of HEC-HMS software showing the imported base model using ArcHydro and HEC-GeoHMS extensions in ArcGIS.
Figure 10. The workspace of HEC-HMS software showing the imported base model using ArcHydro and HEC-GeoHMS extensions in ArcGIS.
Water 13 02667 g010
Figure 11. The setup of the HyperKANAKO model (the green line represents the runout area, and the yellow box indicates the deposition area).
Figure 11. The setup of the HyperKANAKO model (the green line represents the runout area, and the yellow box indicates the deposition area).
Water 13 02667 g011
Figure 12. The setup of dams: (a) The landslide dam with a height of 30.0 m; (b) the closed dam with a height of 5.0 m.
Figure 12. The setup of dams: (a) The landslide dam with a height of 30.0 m; (b) the closed dam with a height of 5.0 m.
Water 13 02667 g012
Figure 13. The flow depth simulation results within the Mesilau watershed.
Figure 13. The flow depth simulation results within the Mesilau watershed.
Water 13 02667 g013
Figure 14. Simulation results of the depositional thickness of the Mesilau watershed.
Figure 14. Simulation results of the depositional thickness of the Mesilau watershed.
Water 13 02667 g014
Table 1. Records of debris flows in Malaysia from 1995 to 2015.
Table 1. Records of debris flows in Malaysia from 1995 to 2015.
No.DateLocationTriggeringFatalities
130 June 1995Km 38.6 Kuala Lumpur–Karak Highway, Genting Sempah, SelangorRainfall20
229 August 1996Kg. Orang Asli, Pos Dipang, Kampar, PerakRainfall44
326 December 1996Keningau, SabahTyphoon302
4January 2000Cameron Highland, PahangRainfall6
522 September 2001Chinchin Village, Gombak, SelangorRainfall1
628 December 2001Channelized Pulai river, Gunung Pulai, JohorRainfall5
728 January 2002Ruan Changkul, Simunjan, SarawakRainfall16
8November 2002Hulu Kelang, outskirt of Kuala Lumpur CityRainfall8
910 November 2003Section 23.3 to 24.10 Kuala Kubu Baru, SelangorRainfall0
102 November 2004Km 52.4 Kuala Lumpur–Karak highway, LentangRainfall0
1110 November 2004Km 302 North South Expressway, G. Tempurung, PerakRainfall0
12November 2004Taman Harmonis, Gombak, SelangorRainfall1
1312 April 2005Km 33 Simpang Pulai, Cameron Highland, PahangRainfall0
1415 November 2007Km 4 to 5 Gap, Fraser’s Hill Road, PahangRainfall0
15January 2008Channel of Fraser Hill’sRainfall0
163 January 2009Section 62.4, Lojing Gua Musang Road, KelantanRainfall0
177 August 2011Kampung Orang Asli Sungai Ruil, Cameron Highlands, PahangRainfall7
1823 October 2013Bertam valley, Cameron HighlandRainfall1
195 November 2014Km 28, Jalan Tamparuli, Ranau, SabahRainfall0
2011 Jane 2015Channel of Fraser Hill’sRainfall0
2118 May 2015Km 38.80, Jalan Penampang Tambunan Dongongan, SabahRainfall0
2215 June 2015Channelized Mesilau river, Kundasang, SabahEarthquake and rainfall0
23August 2015Channelized Kedamaian, and Panataran river, Kota Belud, SabahEarthquake and rainfall0
Sources: [7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28].
Table 2. The types of datasets used for this study.
Table 2. The types of datasets used for this study.
Type of DatasetSourceResolution/RecordDate
OrthophotoLiDAR0.07 m2016
DTM (before debris flow)IfSAR5.0 m2008
DTM (after debris flow)LiDAR0.25 m2016
RainfallRain-gaugeHourly, mmJune 2015
Table 3. The calibrated parameters in the “DEM, Dams, and Observation” tab.
Table 3. The calibrated parameters in the “DEM, Dams, and Observation” tab.
Parameters/VariablesValueUnit
Simulation time40min
Time step1s
Diameter of material1m
Mass density of bed material2650kg/m3
Mass density of fluid1200kg/m3
Concentration of moveable bed0.2-
Internal friction angle35degree
Acceleration of gravity9.8m/s2
Coefficient of erosion rate0.0007-
Coefficient of deposition rate0.05-
Manning’s roughness coefficient0.03s/m1/3
Number of calculation points250-
Interval of 1D calculation points5m
Number of 2D calculation points480 × 160-
Interval of 2D calculation points5 × 5m × m
Table 4. The 11 summarized simulation scenarios by varying the baseflow and hydrographs.
Table 4. The 11 summarized simulation scenarios by varying the baseflow and hydrographs.
ScenarioBaseflow
(m3/s)
Hydrograph
(m3/s)
Simulation Runout Distance
Scenario 1100113.8Short
Scenario 2200213.8Short
Scenario 3300313.8Short
Scenario 4400413.8Short
Scenario 5500513.8Almost reach
* Scenario BF550563.8Best-fit
Scenario 6600613.8Over-estimate
Scenario 7700713.8Over-estimate
Scenario 8800813.8Over-estimate
Scenario 9900913.8Over-estimate
Scenario 1010001013.8Over-estimate
* represents the best-fit simulation result.
Table 5. The flow depth of every zone.
Table 5. The flow depth of every zone.
ZonesColour RangeFlow Depth (m)
Zone A (Initiation area)green to orange1.0–15.0
Zone B (Mesilau village)green to yellow1.0–10.0
Zone C (Transportation area)blue to yellow0.1–10.0
Zone D (Naradau village)blue to green0.1–5.0
Zone E (Liwagu Dam)blue to green0.1–5.0
Table 6. The depositional thicknesses for every zone.
Table 6. The depositional thicknesses for every zone.
ZonesColour RangeDepositional Thickness (m)
Zone A (Initiation area)green to orange1.0–15.0
Zone B (Mesilau village)blue to yellow0.1–10.0
Zone C (Transportation area)blue to green0.1–5.0
Zone D (Naradau village)blue to green0.1–5.0
Zone E (Liwagu dam)blue to green0.1–5.0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rosli, M.I.; Che Ros, F.; Razak, K.A.; Ambran, S.; Kamaruddin, S.A.; Nor Anuar, A.; Marto, A.; Tobita, T.; Ono, Y. Modelling Debris Flow Runout: A Case Study on the Mesilau Watershed, Kundasang, Sabah. Water 2021, 13, 2667. https://doi.org/10.3390/w13192667

AMA Style

Rosli MI, Che Ros F, Razak KA, Ambran S, Kamaruddin SA, Nor Anuar A, Marto A, Tobita T, Ono Y. Modelling Debris Flow Runout: A Case Study on the Mesilau Watershed, Kundasang, Sabah. Water. 2021; 13(19):2667. https://doi.org/10.3390/w13192667

Chicago/Turabian Style

Rosli, Muhammad Iylia, Faizah Che Ros, Khamarrul Azahari Razak, Sumiaty Ambran, Samira Albati Kamaruddin, Aznah Nor Anuar, Aminaton Marto, Tetsuo Tobita, and Yusuke Ono. 2021. "Modelling Debris Flow Runout: A Case Study on the Mesilau Watershed, Kundasang, Sabah" Water 13, no. 19: 2667. https://doi.org/10.3390/w13192667

APA Style

Rosli, M. I., Che Ros, F., Razak, K. A., Ambran, S., Kamaruddin, S. A., Nor Anuar, A., Marto, A., Tobita, T., & Ono, Y. (2021). Modelling Debris Flow Runout: A Case Study on the Mesilau Watershed, Kundasang, Sabah. Water, 13(19), 2667. https://doi.org/10.3390/w13192667

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