Next Article in Journal
Sustainable Micro, Small, and Medium Enterprises: A Developed Model of Entrepreneurial Marketing Behaviors for Thailand’s MSMEs
Previous Article in Journal
Can Green Finance Drive the Development of the Green Building Industry?—Based on the Evolutionary Game Theory
Previous Article in Special Issue
Investigating the Relationship between the Manning Coefficients (n) of a Perforated Subsurface Stormwater Drainage Pipe and the Hydraulic Parameters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessment of Debris Flow Impact Based on Experimental Analysis along a Deposition Area

by
Muhammad Khairi A.Wahab
1,*,
Mohd Remy Rozainy Mohd Arif Zainol
2,*,
Jazaul Ikhsan
3,
Mohd Hafiz Zawawi
4,
Mohamad Aizat Abas
5,
Norazian Mohamed Noor
6,
Norizham Abdul Razak
7 and
Moh Sholichin
8
1
School of Civil Engineering, Engineering Campus, Universiti Sains Malaysia, Nibong Tebal 14300, Penang, Malaysia
2
River Engineering and Urban Drainage Research Centre (REDAC), Engineering Campus, Universiti Sains Malaysia, Nibong Tebal 14300, Penang, Malaysia
3
Department of Civil Engineering, Universitas Muhammadiyah Yogyakarta, Yogyakarta 55183, Indonesia
4
College of Engineering, Universiti Tenaga Nasional, Bandar Baru Bangi 43650, Selangor, Malaysia
5
School of Mechanical Engineering, Engineering Campus, Universiti Sains Malaysia, Nibong Tebal 14300, Penang, Malaysia
6
Centre of Excellence Geopolymer and Green Technology (CEGeoGTech), Faculty of Civil Engineering & Technology, Universiti Malaysia Perlis, Arau 02600, Perlis, Malaysia
7
School of Aerospace Engineering, Engineering Campus, Universiti Sains Malaysia, Nibong Tebal 14300, Penang, Malaysia
8
Department of Chemical Engineering, Faculty of Engineering, Universitas Brawijaya, Malang 65145, Indonesia
*
Authors to whom correspondence should be addressed.
Sustainability 2023, 15(17), 13132; https://doi.org/10.3390/su151713132
Submission received: 9 March 2023 / Revised: 5 June 2023 / Accepted: 17 June 2023 / Published: 31 August 2023

Abstract

:
Debris flow is a devastating phenomenon that happens in hilly and mountainous regions and has a serious impact on affected areas. It causes casualties and serious damage to the environment and society. Therefore, a susceptible assessment is necessary to prevent, mitigate, and raise awareness of the impact of debris flows. This paper focuses on evaluating the deposition area along the deposition board. The methodology involved an experiment on a physical model by demonstrating the debris flow based on the steepness of the flume slope at 15°, 20°, and 25° angles. The limestone particles with a total volume of 2.5 × 106 mm3 acted as debris and were released with water from the tank to the deposition board with an area of 10 × 105 mm2. The volume, area, and length of particle distribution carried from the flume to the deposition board were then determined. Based on the experimental results, the deposition board is covered with particles of about 696.19 × 103 mm3, 748.29 × 103 mm3, and 505.19 × 103 mm3 volume for each 15°, 20°, and 25° angle, respectively. In actual situations, debris flow is capable of causing significant risk to the affected area. This study can be deemed useful for a risk assessment approach, to help develop guidelines, and to mitigate the regions where debris flows are most probable to occur.

1. Introduction

In mountainous areas where uphill valleys or streams are fed by heavy annual rainfall, debris flow has been recognized as a natural phenomenon in recent years. Due to its size (volume of transported debris) and long runout, the phenomenon is among the most dangerous natural processes that affect the mountainous environment [1]. Debris flow contains multiple mixtures of solids, soil, wood, gravel, rocks, snow, ice, and water that slide down slopes under the driving force of gravity [2]. These flows are frequently reported to result in major property damage, loss of life, and morphological changes along riverbeds and mountain slopes [1,3,4].
In many areas of the country, major flooding and landslides were caused by significant amounts of rain that were recorded as a result of prolonged rainfall events [5,6]. Debris flows are often caused by many different mechanisms, one of which is heavy rain, which turns landslides into fast-moving liquid sediments [6,7,8].
Debris flow formation consists of three stages: the initiation area, the propagation zone, and the deposition area [1,9].
The phase of initiation is where the initial mass is released. The most typical type is debris flow, which is mobilized from a single large landslide or from a large series of smaller landslides. They occur when a landslide or debris slide turns into a debris flow. Mobilization is the process of creating debris flows from a static mass of liquid soil, silt, or rock. A major contributing component to the development of a debris flow is the impact of surface-water runoff [9,10,11,12,13]. Debris flows can also be caused by steep channels with an abundance of sediment, shallow landslides, and frequent in high-mountain regions [14,15,16]. The triggered process depends on several characteristics, including the geomorphology, geotechnical characteristics of slope inclination, hydrological factors, and geological setting [8,17]. It is of utmost importance to investigate and understand the debris flow initiation mechanisms with variously sophisticated circumstances for developing debris flow mitigation [18,19].
Once debris flows initiate, they travel down the channel. The zone of propagation is the path where the debris flow after initiation has been triggered, usually by a few factors. One of the main characteristics of many debris flows is the entrainment of channel path and torrent flank material as well as sediment deposition during runout. Through abrupt changes in the flow volume and its rheological behavior, such entrainment mechanisms have the ability to drastically change the mobility of the flow [1]. For researchers to predict the behavior of a debris flow in this zone, parameters such as runout distance are essential [20]. Efforts to analyze, predict, and evaluate the runout flow are necessary for identifying the hazardous region that can be potentially impacted and determining the intensity of the disaster for hazard and risk assessment. Prevention strategies can also be conducted, as can mitigation designs, through dynamic analysis [21,22].
Deposition zones, also known as alluvial fans, typically have many fan levels that alternate between aggradation and incision. There are situations where it could develop hundreds of kilometers away from its source, and deposition of the debris is the final stage of debris flows [23]. Channelized debris flows can rapidly transport large volumes of sediment to the depositional zone, endangering nearby residents’ infrastructure, such as buildings and roads, as well as the environment [24,25,26].
Mitigation and prevention are necessary as they involve lives and the destruction of property. Therefore, it is very important to understand the initiation, propagation, and deposition of the debris flow, and it can be a dominant tool to reduce and minimize the effects of disasters [27].
Debris flow incidents are only investigated and reported when they have already resulted in damage to people or property. Currently, research on debris flows is limited to post-disaster analysis at the catastrophe site. The catastrophes ultimately lead to huge financial losses and fatalities [15]. There have been comparatively few observations of debris flow events in the field. The mechanisms and characteristics of debris flows have long been poorly understood [1]. Given a specific shape and collection of materials, observations of laboratory experiments utilizing model-scale flumes can help us understand the mechanics of debris flows better and, as a result, help us quantify the spatial component of risk. Furthermore, parametric investigations can be carried out in a flume experiment, which is simpler and highly controlled.

2. Materials and Methods

2.1. Physical Model Test

The model for this study was developed with three main components: a water tank, a flume, and a deposition board. In this design, the water tank is filled manually using a pump system, and flow to the tank is controlled by a butterfly valve connected to the pipe. The actual model is 4.55 m long, 1.1 m wide, and 1.5 m high. The water tank is 1 m long, 1 m wide, and 1 m high. The total length of the flume is 2.5 m long, 0.1 m wide, and 0.01 m high. The deposition board dimension is 1 m long and 1 m wide. The model is constructed using a polyvinyl chloride (PVC) sheet with a 10 mm thickness. Figure 1 shows the geometry of the physical model, and Figure 2 shows the physical model in the laboratory.

2.2. Physical Model Configuration

The determination of the design parameters of the model needs to be considered, such as the flume angle and water level in the water tank.
Twelve cases were included in the study and tested at different water levels, flume angles, and water gate openings. Table 1 summarizes the case studies and operations used in this section.
In this research, the water and sediment are separated in the experiment. In order to study these processes in isolation, we systematically varied the water content (height of water in the tank) of the bed sediment. Through experiments, we were able to observe and measure the effect of different water contents on the process of entrainment. According to Calhoun and Clague [28], water-rich mass flows are one of the causes of debris flow. This process involves the movement of a continuous mass of debris, resembling a flood, which is directed by an inflow. The primary mode of transportation for the majority of sediments within this flow is particle-by-particle movement, forming what is known as a bed load. The water in the tank was used to reproduce the conditions of a debris flow. In their research, [29] utilized this method to study the segregation of particle routing in debris flow mechanisms.

2.2.1. Sediment Particle and Preparation

The limestone used in this study was obtained as a by-product of the marble industry in Ipoh, Malaysia. Prior to undertaking the experiments, a general characterization of the limestone was done to determine its particle size. In order to get the necessary particle size for the media, about 1 kg of limestone was sieved using a sieve shaker. The particles were sieved to pass through a 1.18 mm aperture and be retained in a 2 mm size range. The limestone used in the experiment had a density of 2573 kg/m3. Figure 3 shows the total sample collection as well as the sieve shaker that was used to sieve the samples. The selection of debris flow material offers valuable data insights, especially when focusing on limestone. This is particularly significant in the Malaysian region, where limestone is known to have extensive formation and prevalence [30].
Sediment is filling behind the sediment block in the flume, which is 50 mm high, 500 mm long, and 100 mm wide. The distance between the bottom sediment block and the bottom of the flume is 1250 mm. The sediment block is higher than the sample to maintain its position before the test is carried out. The details of the sediment are shown in Figure 4.

2.2.2. Water Tank

The water in the tank will be filled based on the test that will be conducted according to the height of the water, either 120 mm or 150 mm. The water in this tank resembles the pore water pressure in the soil. The increase in water content leads to a decrease in the shear strength of the soil in the process of rainfall infiltration. Eventually, the balance of the debris source fails, triggering a debris flow [19]. The pump will be turned on to fill the tank with water and controlled using a valve to reach the desired water level. The marker for the water level in the water tank is shown in Figure 5.

2.2.3. Water Gate Opening

For this study, there are two openings (50 mm and 100 mm) controlled by a water gate using pneumatic actuators. The pneumatic actuators converted the energy of compressed air into a mechanical motion that regulates the water gate. This opening is to see the influence of the amount of water on the sediment on the dynamics of debris flow, where water has an extensive influence on it [31]. Figure 6 shows the position of the water gate system and water opening on the physical model.
The pneumatic actuator is connected to the control box with a push button to control the opening and closing of the water gate. The control box comprises a solenoid valve, circuit board, power supply, and push button that have their respective functions to ensure that it works as desired. Figure 7 indicates the electronic device on the control box.

2.2.4. Inclination of Flume

The degree of the slope on the flume is adjusted depending on the case. Prior to conducting the test procedure, the slope of the flume is set. In this study, three levels of flume angle (°) were studied to represent the initiation level of debris flow: high (25°), intermediate (20°), and low (15°) levels. This is based on the fact that the average slope angle of the channel bed at the prospective area for debris flows typically ranges from 10° to 25° [8,20,32,33,34]. The inclination of the flume is shown in Figure 8.

2.2.5. Deposition Board

A deposition board is located at the end of the flume to receive particle deposits from the flow current. An alluvial fan will be formed, and the contour will be measured as a grid. The boards are kept in clean condition, in good alignment, and positioned properly to avoid data distortion.
Equation (1), developed by [20], was used to determine the position and volume of the sediment in the flume and the size of the deposition board. Three runout distances calculated for three different slopes met the equation with a runout distance of less than 1000 mm.
L = 1.9V0.16 H0.83
where L is the propagation distance, H is the height where the debris flow occurred, and V is the volume of the debris flow. Based on the height and volume of the debris flow occurrence, the behavior of the debris flow was studied. As a result, the outcomes were based on a topographical viewpoint.

2.3. Data Collection

2.3.1. Shape and Thickness

In order to have an accurate condition of the sediment deposit, the grid method is used. In this method, a rigid bamboo stick with a diameter of 4 mm and a length of 15 mm with red and white colors is used as the interval measurement. Each interval is 20 mm to facilitate the process of retrieving data for thickness. The stick is placed as a grid on a board, 50 × 50 mm for the middle part and 100 × 100 mm for the right, left, and end parts of the board. The middle part of the grid is smaller because more deposits occur in this part than in other parts. Figure 9 shows the grid points on the deposition board and the labels on the deposition board for data retrieval purposes. In research conducted by researchers [35], it was determined that bamboo sticks have no impact on the runout and pattern of debris flows. The study aimed to investigate the efficacy of tree stems in intercepting debris flows in forested fan areas. The researchers observed that the high-forest management type, which was represented by bamboo sticks in this study, did not appear to have a significant impact on reducing debris-flow motion for solid concentrations greater than or equal to 0.60.
A measuring tape was used to measure the shape deposit of each grid on the x and y axes. A digital caliper depth rod was used to measure the thickness of the deposit on the deposition board and lower jaws were used to measure the deposit shape for parts that required high-accuracy readings to get the actual shape.

2.3.2. Contour Sketch and Visualize

Once the data were collected, they were processed to generate a contour map using the Kriging method. One of the advantages of kriging is that it allows for quantifying the estimation uncertainty, which is important for decision-making processes. Kriging is often regarded as an effective approach because it can produce accurate maps for a wide range of datasets. Additionally, it has the capability to compensate for clustered data by assigning less weight to the cluster in the overall prediction [36].
The resulting contour maps were used to further understand deposition shapes and patterns. Cross-sectional profiles, volumes, areas, and distances can also be computed and exported. This can simplify the analysis process. In this study, Surfer was utilized as a software tool for post-processing and analysis of the collected data. Surfer is a comprehensive and powerful software package developed by Golden Software, which specializes in providing advanced scientific data visualization tools. Specifically, Surfer is designed to create and visualize high-quality 2D and 3D maps, models, and analyses of spatial data, making it a suitable choice for a wide range of scientific applications [37,38]. Figure 10 shows some examples of post-processing results.

2.4. Test Procedure

This section explains the process of the experiment. The stages of the experimental process chart are shown in Figure 11. This chart shows how the experiment is carried out.
The experiment is conducted as follows. Firstly, the water gate was closed using the pneumatic actuator with an on/off switch button. The flume is tilted to the desired inclination according to the case to be carried out. The water tank is filled with water at the desired level, according to the case. The deposition board is ensured to be in its proper position and in good alignment by using a spirit level. The sediment is then prepared in the flume. Once everything is ready, the switch is pressed to open the water gate for 5 s. Then, the switch is pressed again to close back the water gate.
The next step is to take data for the sediment contour on the deposition board using an electronic caliper and measuring tape. Contour data are recorded in x, y, and z coordinate formats. The recorded data will then visualize the contour in 3D.

3. Results and Discussions

3.1. Deposition Observation

A total of 12 cases were analyzed. The contouring was performed with a measuring tape and digital calipers. As a result of the data analysis, a 3D contour map was generated to represent the spatial distribution of the sampled data. Then, a planar (2D) view of the cross-section of the 3D contour was created to provide a more detailed view of the distribution pattern along a specific axis or plane. Each color in the contour line represents an elevation value, which is measured in millimeters (mm). According to [21], the deposition pattern is unaffected by any obstructions or discrete channels on the fan. Typically, the topography causes the material to rise over the channel and begin to spread over the fan. Table 2 displays the deposition pattern as classified by the observations of [21]. Table 2 categorizes the results of the deposition pattern from the experimental observation. This process is necessary to comprehend the particle characteristics, particle distribution, and physical data of the deposited material.
The researchers [21] made an important observation regarding the debris-flow deposition pattern known as type A. They identified a distinct correlation between the properties of the circle sector (specifically the radius, denoted as Lf, and the aperture angle, denoted as Ψ) and the volume of the event (V). This correlation allowed them to propose a theoretical approximation for the radius, Lf in Equation (2). Moreover, they conducted data analysis and successfully fitted the data with a power-law function, resulting in a coefficient of determination (R2) value of 0.78. As depicted in Figure 12, it was determined that type A represents an ideal shape for achieving uniform runout patterns.
L fpred = 65.3 Ψ V 2 / 3 0.49 65.3 Ψ 1 / 2 V 1 / 3
Due to the observational outcomes being nearly identical for all four cases, only two of the four cases on each slope of the flume (15°, 20°, and 25°) representing the water level (120 mm, 150 mm) with half gate opening are displayed in Figure 13.
As shown in Figure 13, the deposition trajectory reaches the end of the deposition board for all 6 cases. From case A1 to case C11, the deposition pattern reveals that the fan size is decreasing. According to the details, A1 recorded the highest estimated total volume of deposition, which is 1.677518 × 103 mm3, with a surface area of 537.07 × 103 mm2. The total volume of particles on the deposition board decreases as the inclination angle increases, and this trend can be seen in all conditions.
In every case, the pattern type is categorized as type A because it represents the ideal shape for uniform runout patterns because there are no obstacles or distinct channels on the fan affecting deposition.
In the conducted research, Figure 14 visually presents the outcome following the completion of each case. Mainly, no sediment remains either behind the sediment block or in front of it once the water gate is closed, as clearly depicted in Figure 14a. Furthermore, Figure 14b effectively illustrates the condition of the deposition board subsequent to the complete deposition of sediment. The deposit’s shape on the board is clearly observable without any obstructions caused by the bamboo stick. Additionally, the total volume of sediment that traverses over the board is visibly captured and measured through the retention on the strainer, with the corresponding values recorded in Table 2.
Table 3 shows the summary results for all experimental cases. The table consists of the total volume of deposited sediment, total volume of sediment retained in the strainer, planar area, surface area, aperture angle, highest point, and percentage of particles on the deposition board. The effect of the slope gradient on the flume and water level inside the tank can be seen in Table 3.
From Table 3, the maximum percentage of particles covering the sedimentation board is 53.68% for A1 and the minimum percentage is 29.00% for B8, respectively, at angles of 15° and 25°. The decrease in the percentage of particles on the deposition board from case A1 to case C12 can be attributed to the change in slope gradient (15° to 25°) at the flume and water level (120 to 150 mm) in the tank. This decrease in the percentage of particles on the deposition board could be related to an increase in pore water pressure, causing entrainment to occur and separating the particles, making them liquefy and accelerate in motion during runout. This causes the sediment to flow out of the deposition board, as supported by the total volume results. This situation shows a higher mobility and results in flatter deposits, with only one case recording the highest point on the contour at 11 mm and the rest below 7 mm. Additionally, larger flume slopes produce larger runout distances and reduce the area due to greater gravitational potential energy.
According to our observations, the deposition fan in case B5 exhibits the widest aperture angle, measuring 124.88°, while the smallest aperture angle of 65.86° is found in case C12. In case C5, there is a flume inclined at 20° with a water level of 120 mm and the water gate partially open. Similarly, in Case C12, the flume has a 25° inclination, a water level of 150 mm, and the gate is fully open. By examining the case with a 15° inclination, we can observe a reduction in angle from Case A1 to A4. This pattern of angle reduction is also evident in cases with 20° and 25° of inclination.
As a result of predicting the debris flow length using Equation (2), it was found that there were only two cases where the runout length exceeded 1000 mm compared to all the other cases. During the observation, it was found that some of the sediment that has been deposited will be carried to the end of the board by the excess water flowing from the flume after the gate is open for 5 s. This indicates that the duration of the open gate affects the movement of sediment in the flume.

3.2. Statistical Analysis

This analysis explores the relationships between flume inclination (θd) and fluid volume (VW) as independent variables and total volume of deposited sediment (VT), total volume of sediment retained in strainer (Vout), planar area (PA), surface area (SA), aperture angle (Ψ), highest point (HP) and percentage particle on deposition board (PB) as dependent variables using Pearson’s correlation coefficient (r) and regression analysis.
Several variables exhibited a significant correlation with the inclination of the flume. For instance, a weak positive correlation was observed between flume inclination and the aperture angle (r = 0.302, p < 0.05). Furthermore, the inclination demonstrated a weak negative correlation with the total volume of sediment retained in the strainer (r = −0.005, p < 0.05), as shown in Table 4. Notably, a strong negative correlation was found between the flume inclination and the total volume of deposited sediment (r = −0.637, p < 0.05). It was evident that the distribution of sediment deposition in terms of total volume was influenced by the inclination of the flume and indicated that when θd is increased, there is a clear correlation with a decrease in the volume of deposits present on the board. In a study conducted by [39], the experimental findings provided support for the idea that as the channel slope, which represents the flume inclination, increased, the total travel distance of debris flow also increased. Aside from flume inclination, fluid volume also showed a significant association with VT (r = −0.260, p < 0.05), Vout (r = 0.783, p < 0.01), PA (r = −0.210, p < 0.05), SA (r = −0.210, p < 0.05), Ψ (r = −0.285, p < 0.05), Hp (r = 0.107, p < 0.05), and PB (r = −0.210, p < 0.05). A significant positive correlation between Vout (fluid volume) and the runout distance of debris flow, suggesting that higher fluid volumes contribute to increased travel distances. According to a study carried out by [40], it has been demonstrated that the phenomenon of entrainment exerts a noteworthy impact on flow depths and runout. The study reveals a substantial augmentation in both maximum flow heights and the distance covered by the flow when entrainment is taken into consideration. In fact, the calculated maximum height of the flow, when entrainment is absent, can be multiplied by two or even three times when compared to the estimated maximum heights that incorporate the entrainment process. Similarly, the length of the runout distance traveled can be expanded by nearly one-third. Moreover, this study reveals that a higher fluid volume also leads to a greater retention of sediment on the strainer.
In [41], the researchers mentioned that the polynomial regression model has the potential to accurately predict responses in certain data segments. Therefore, for the analysis of θd and Vw, this research explored alternative regression models, including linear, quadratic, and nonlinear approaches. The subsequent equations were formulated to account for each type of response.
y = a + b   ×   t ,
y = a + b   ×   t + c   × t 2 ,
y = exp b 1   ×   x / b 2 + b 3   ×   x
Several datasets exhibit poor conditioning when applied to a quadratic model. Fits for the linear and quadratic are reasonably comparable for the θd with VT (R2 = 40.6, 40.7). The comparative analysis of the nonlinear model reveals its superiority in estimating the values of θd and VW for the aperture angle data, as indicated by its significantly lower standard error of the regression (S) when compared to alternative models. However, it should be noted that the nonlinear model does not provide a coefficient of determination (R2) to assess the goodness of fit. Table 5 shows the regression goodness of fit for several variables.

4. Conclusions

Data collected on all 12 cases were used as input to be analyzed and visualized. The results from the analysis of all 12 cases show the maximum percentage of particles covering the deposition board is 1.677518 × 103 mm3 (53.68%) for A1 and the minimum percentage is 1.048050 × 103 mm3 (29.00%) for B8. The decrease in the percentage of particles on the deposition board from case A1 to case C12 can be attributed to the change in slope gradient (15° to 25°) at the flume and water levels (120 to 150 mm) in the tank. Evidently, an increase in water content and pore water pressure leads to entrainment, which separates the particles, causes them to liquefy, and accelerates the motion during runout. As a result, the sediment was driven off the deposition board.
The research concludes that even though the affected area comprises only 29% or the minimum with a small aperture angle, the strainer retains a larger amount of sediment. This suggests that the debris flow has the potential to travel longer distances. These findings have significant implications for the development of guidelines, mitigation plans, and risk assessment techniques in regions with extensive limestone formations and where debris flows are expected to occur.

Author Contributions

Conceptualization, M.K.A.W. and M.R.R.M.A.Z.; methodology, M.K.A.W. and N.A.R.; software, M.K.A.W. and M.A.A.; validation, N.A.R., J.I. and M.H.Z.; formal analysis, M.K.A.W., N.M.N. and M.S.; investigation, M.K.A.W., N.A.R., J.I., M.A.A. and M.H.Z.; resources, M.K.A.W.; data curation, M.K.A.W.; writing—original draft preparation, M.K.A.W.; writing—review and editing, M.K.A.W.; visualization, M.K.A.W.; supervision, M.R.R.M.A.Z.; project administration, M.R.R.M.A.Z., N.A.R., J.I., M.H.Z. and M.A.A.; funding acquisition, M.R.R.M.A.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Ministry of Higher Education Malaysia for Fundamental Research Grant Scheme with Project Code: FRGS/1/2016/TK01/USM/02/4.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Acknowledgement to “Ministry of Higher Education Malaysia under the Highest Institution Center of Excellence (HICOE) research grant (311.PREDAC.4403901) and Ministry of Higher Education Malaysia for Fundamental Research Grant Scheme with Project Code: FRGS/1/2016/TK01/USM/02/4”.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Takahashi, T. A Review of Japanese Debris Flow Research. Int. J. Eros. Control Eng. 2013, 2, 1–14. [Google Scholar] [CrossRef]
  2. Zhou, G.D.; Law, R.P.H.; Ng, C.W.W. The mechanisms of debris flow: A preliminary study. In Proceedings of the 17th International Conference on Soil Mechanics and Geotechnical Engineering: The Academia and Practice of Geotechnical Engineering, Alexandria, Egypt, 5–9 October 2009; Volume 2, pp. 1570–1573. [Google Scholar]
  3. Mohd Arif Zainol, M.R.R.; AWahab, M.K. Hydraulic Physical Model of Debris Flow for Malaysia Case Study. In Proceedings of the IOP Conference Series: Materials Science and Engineering, Iasi, Romania, 17–18 May 2018; Institute of Physics Publishing: Bristol, UK, 2018; Volume 374. [Google Scholar]
  4. Tiranti, D.; Crema, S.; Cavalli, M.; Deangeli, C. An Integrated Study to Evaluate Debris Flow Hazard in Alpine Environment. Front. Earth Sci. 2018, 6, 60. [Google Scholar] [CrossRef]
  5. Muhammad, N.S.; Julien, P.Y.; Salas, J.D. Probability Structure and Return Period of Multiday Monsoon Rainfall. J. Hydrol. Eng. 2016, 21, 0401548. [Google Scholar] [CrossRef]
  6. Rahman, H.A.; Mapjabil, J. Landslides Disaster in Malaysia: An Overview. Health Environ. J. 2017, 8, 58–71. [Google Scholar]
  7. Joe, E., Jr.; 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]
  8. Kasim, N.; Taib, K.A.; Mukhlisin, M.; Kasa, A. Triggering Mechanism and Characteristic of Debris Flow in Peninsular Malaysia. Am. J. Eng. Res. 2016, 5, 112–119. [Google Scholar]
  9. Liu, K.F.; Huang, M.C. Numerical simulation of debris flows. In Proceedings of the International Conference on Offshore Mechanics and Arctic Engineering—OMAE, Honolulu, HI, USA, 31 May–5 June 2009; Volume 6, pp. 375–382. [Google Scholar]
  10. Chang, T.C.; Wang, Z.Y.; Chien, Y.H. Hazard assessment model for debris flow prediction. Environ. Earth Sci. 2010, 60, 1619–1630. [Google Scholar] [CrossRef]
  11. Li, M.L.; Jiang, Y.J.; Yang, T.; Huang, Q.B.; Qiao, J.P.; Yang, Z.J. Early warning model for slope debris flow initiation. J. Mt. Sci. 2018, 15, 1342–1353. [Google Scholar] [CrossRef]
  12. Vargas-Cuervo, G.; Rotigliano, E.; Conoscenti, C. Prediction of debris-avalanches and -flows triggered by a tropical storm by using a stochastic approach: An application to the events occurred in Mocoa (Colombia) on 1 April 2017. Geomorphology 2019, 339, 31–43. [Google Scholar] [CrossRef]
  13. Zhou, J.W.; Cui, P.; Yang, X.G.; Su, Z.M.; Guo, X.J. Debris flows introduced in landslide deposits under rainfall conditions: The case of Wenjiagou gully. J. Mt. Sci. 2013, 10, 249–260. [Google Scholar] [CrossRef]
  14. Chen, S.C.; Wu, C.Y. Debris flow disaster prevention and mitigation of non-structural strategies in Taiwan. J. Mt. Sci. 2014, 11, 308–322. [Google Scholar] [CrossRef]
  15. Lay, U.S.; Pradhan, B. Identification of Debris Flow Initiation Zones Using Topographic Model and Airborne Laser Scanning Data. In Lecture Notes in Civil Engineering, Proceedings of the GCEC 2017: The 1st Global Civil Engineering Conference 1, Kuala Lumpur, Malaysia, 25–28 July; Pradhan, B., Ed.; Springer: Singapore, 2019; Volume 9. [Google Scholar]
  16. Pradhan, B.; Bakar, S.B.A. Debris Flow Source Identification in Tropical Dense Forest Using Airborne Laser Scanning Data and Flow-R Model. In Laser Scanning Applications in Landslide Assessment; Pradhan, B., Ed.; Springer: New York, NY, USA, 2017; pp. 85–112. [Google Scholar]
  17. Borga, M.; Stoffel, M.; Marchi, L.; Marra, F.; Jakob, M. Hydrogeomorphic response to extreme rainfall in headwater systems: Flash floods and debris flows. J. Hydrol. 2014, 518, 194–205. [Google Scholar] [CrossRef]
  18. Liu, X.; Wang, F.; Nawnit, K.; Lv, X.; Wang, S. Experimental study on debris flow initiation. Bull. Eng. Geol. Env. 2020, 79, 1565–1580. [Google Scholar] [CrossRef]
  19. Zhuang, J.Q.; Cui, P.; Peng, J.B.; Hu, K.H.; Iqbal, J. Initiation process of debris flows on different slopes due to surface flow and trigger-specific strategies for mitigating post-earthquake in old Beichuan County, China. Environ. Earth Sci. 2013, 68, 1391–1403. [Google Scholar] [CrossRef]
  20. Jeong, S.; Lee, J.; Song, C.G.; Lee, S.O. Propagation characteristics of debris flows considering entrainment effect. Int. J. Eng. Technol. 2018, 7, 122–128. [Google Scholar] [CrossRef]
  21. Rickenmann, D.; Scheidl, C. Debris-Flow Runout and Deposition on the Fan. In Dating Torrential Processes on Fans and Cones; Springer: Dordrecht, The Netherlands, 2012; Volume 47, pp. 75–93. [Google Scholar]
  22. Wang, W.; Chen, G.; Han, Z.; Zhou, S.; Zhang, H.; Jing, P. 3D numerical simulation of debris-flow motion using SPH method incorporating non-Newtonian fluid behavior. Nat. Hazards 2016, 81, 1981–1998. [Google Scholar] [CrossRef]
  23. Talling, P.J.; Wynn, R.B.; Masson, D.G.; Frenz, M.; Cronin, B.T.; Schiebel, R.; Amy, L.A. Onset of submarine debris flow deposition far from original giant landslide. Nature 2007, 450, 541–544. [Google Scholar] [CrossRef]
  24. Bouchut, F.; Vilotte, J.P.; Lajeunesse, E.; Aubertin, A.; Pirulli, M. On the use of Saint Venant equations to simulate the spreading of a granular mass. J. Geophys. Res. Solid Earth 2005, 110, 1–17. [Google Scholar]
  25. Jakob, M.; Stein, D.; Ulmi, M. Vulnerability of buildings to debris flow impact. Nat. Hazards 2012, 60, 241–261. [Google Scholar] [CrossRef]
  26. Stancanelli, L.M.; Musumeci, R.E. Geometrical characterization of sediment deposits at the confluence of mountain streams. Water 2018, 10, 401. [Google Scholar] [CrossRef]
  27. Liu, J.; Nakatani, K.; Mizuyama, T. Effect assessment of debris flow mitigation works based on numerical simulation by using Kanako 2D. Landslides 2013, 10, 161–173. [Google Scholar] [CrossRef]
  28. Calhoun, N.C.; Clague, J.J. Distinguishing between debris flows and hyperconcentrated flows: An example from the eastern Swiss Alps. Earth Surf. Process. Landf. 2018, 43, 1280–1294. [Google Scholar] [CrossRef]
  29. Yamashiki, Y.; Rozainy, M.A.Z.M.R.; Matsumoto, T.; Takahashi, T.; Takara, K. Particle Routing Segregation of Debris Flow Mechanisms Near the Erodible Bed. APCBEE Procedia 2013, 5, 527–534. [Google Scholar] [CrossRef]
  30. Tan, B.K. Environmental Geology of Limestone in Malaysia. IGCP 2002, 448, 54–67. [Google Scholar]
  31. Turnbull, B.; Bowman, E.T.; McElwaine, J.N. Debris flows: Experiments and modelling. Comptes Rendus Phys. 2015, 16, 86–96. [Google Scholar] [CrossRef]
  32. Chang, T.C.; Chien, Y.H. The application of genetic algorithm in debris flows prediction. Environ. Geol. 2007, 53, 339–347. [Google Scholar] [CrossRef]
  33. Lin, D.G.; Hsu, S.Y.; Chang, K.T. Numerical simulations of flow motion and deposition characteristics of granular debris flows. Nat. Hazards 2009, 50, 623–650. [Google Scholar] [CrossRef]
  34. Reid, M.E.; Iverson, R.M.; Logan, M.; Lahusen, R.G.; Godt, J.W.; Griswold, J.P. Entrainment of bed sediment by debris flows: Results from large-scale experiments. Ital. J. Eng. Geol. Environ. 2011, 367–374. [Google Scholar]
  35. Bettella, F.; Michelini, T.; D’Agostino, V.; Bischetti, G.B. The ability of tree stems to intercept debris flows in forested fan areas: A laboratory modelling study. J. Agric. Eng. 2018, 49, 42–51. [Google Scholar] [CrossRef]
  36. Maris, F.; Potouridis, S.; Angelidis, P. Kriging Interpolation Method for Estimation of Continuous Spatial Distribution of Precipitation in Cyprus. Br. J. Appl. Sci. Technol. 2013, 3, 1286–1300. [Google Scholar] [CrossRef]
  37. Marschallinger, R. Interface programs to enable full 3-D geological modeling with a combination of AutoCAD and SURFER. Comput. Geosci. 1991, 17, 1383–1394. [Google Scholar] [CrossRef]
  38. Yang, C.-S.; Kao, S.-P.; Lee, F.-B.; Hung, P.-S. Twelve different interpolation methods: A case study of Surfer 8.0. Proc. XXth ISPRS Congr. 2004, 35, 778–785. [Google Scholar]
  39. Chang, H.; Ryou, K.; Lee, H. Debris flow characteristics in flume experiments considering berm installation. Applied Sciences 2021, 11, 2336. [Google Scholar] [CrossRef]
  40. Luna, B.Q.; Remaître, A.; van Asch, T.W.J.; Malet, J.P.; van Westen, C.J. Analysis of debris flow behavior with a one dimensional run-out model incorporating entrainment. Eng. Geol. 2012, 128, 63–75. [Google Scholar] [CrossRef]
  41. Song, E.; Oh, M.S.; Billard, L.; Moss, A.; Pesti, G.M. A statistical analysis of nonlinear regression models for different treatments for layers. Poult. Sci. 2002, 101, 102004. [Google Scholar] [CrossRef]
Figure 1. The geometry of the physical model: (a) three-dimensional view of the model; (b) geometry model. (Unit: mm).
Figure 1. The geometry of the physical model: (a) three-dimensional view of the model; (b) geometry model. (Unit: mm).
Sustainability 15 13132 g001
Figure 2. Physical model in the laboratory.
Figure 2. Physical model in the laboratory.
Sustainability 15 13132 g002
Figure 3. (a) Sieve shaker; (b) total sample collection.
Figure 3. (a) Sieve shaker; (b) total sample collection.
Sustainability 15 13132 g003
Figure 4. The dimension details of sediment.
Figure 4. The dimension details of sediment.
Sustainability 15 13132 g004
Figure 5. Water level marker in water tank.
Figure 5. Water level marker in water tank.
Sustainability 15 13132 g005
Figure 6. Water gate opening system on model.
Figure 6. Water gate opening system on model.
Sustainability 15 13132 g006
Figure 7. Electronic device on the control box.
Figure 7. Electronic device on the control box.
Sustainability 15 13132 g007
Figure 8. The angle of the flume is changed according to the case to be carried out.
Figure 8. The angle of the flume is changed according to the case to be carried out.
Sustainability 15 13132 g008
Figure 9. Deposition board: (a) grid points; (b) labels.
Figure 9. Deposition board: (a) grid points; (b) labels.
Sustainability 15 13132 g009
Figure 10. Example of post-processing results: (a) generated contour map; (b) cross-sectional profile from the contour maps.
Figure 10. Example of post-processing results: (a) generated contour map; (b) cross-sectional profile from the contour maps.
Sustainability 15 13132 g010
Figure 11. Overview of the test workflow.
Figure 11. Overview of the test workflow.
Sustainability 15 13132 g011
Figure 12. Pattern A in the runout exhibits a Ψ symbol representing the aperture angle of a circular sector. Lf corresponds to the radius from the initial point to the farthest point of the observed deposition. The start point of the deposition and commonly known as the fan apex (modified after [21]).
Figure 12. Pattern A in the runout exhibits a Ψ symbol representing the aperture angle of a circular sector. Lf corresponds to the radius from the initial point to the farthest point of the observed deposition. The start point of the deposition and commonly known as the fan apex (modified after [21]).
Sustainability 15 13132 g012
Figure 13. Post-processing results: (af) the visualized contour map and a cross-sectional profile for the deposition pattern.
Figure 13. Post-processing results: (af) the visualized contour map and a cross-sectional profile for the deposition pattern.
Sustainability 15 13132 g013aSustainability 15 13132 g013b
Figure 14. The final condition after the water gate is closed: (a) no sediment left at the flume behind the block; (b) retained sediment in strainer.
Figure 14. The final condition after the water gate is closed: (a) no sediment left at the flume behind the block; (b) retained sediment in strainer.
Sustainability 15 13132 g014
Table 1. Investigated cases.
Table 1. Investigated cases.
Case StudyFlume Degree (°)Water Level (mm)Gate Opening
A115120Half
A215120Full
A315150Half
A415150Full
B520120Half
B620120Full
B720150Half
B820150Full
C925120Half
C1025120Full
C1125150Half
C1225150Full
Table 2. Debris flow deposition pattern categorization [21].
Table 2. Debris flow deposition pattern categorization [21].
TypeRemarks
ARepresents an ideal shape for uniform runout patterns with no obstructions or distinct channels on the fan affect deposition. Typically, the topography causes the material to rise over the channel and begin to spread over the fan.
B1Obstructions such as structures or mitigation measures in the lower part of the deposition region affect the deposition pattern.
B2
C1Linear structures crossing the deposition zone, such as roads, railway tracks, or receiving streams that caused a pronounced lateral spreading in the lower part of the deposition region.
C2
D
Table 3. Summary of the results.
Table 3. Summary of the results.
Case StudyTotal
Volume of Deposited Sediment
(mm3)
Total
Volume of Sediment Retained in Strainer
(mm3)
Planar Area
(mm2)
Surface Area
(mm2)
Aperture Angle (°)Highest Point (mm)Percentage Particle on Deposition Board (%)Type of PatternLfpred [21] (mm)
A11.677518 × 1032498.322 × 103536.79 × 103537.07 × 10390.197.553.68%A884.87
A21.508610 × 1032498.491 × 103436.93 × 103437.13 × 10386.636.543.69%A902.50
A31.656639 × 1032498.343 × 103488.42 × 103488.65 × 10381.836.548.84%A928.07
A41.544908 × 1032498.455 × 103402.34 × 103402.61 × 10367.877.540.23%A1017.15
B51.497587 × 1032498.502 × 103461.38 × 103461.57 × 103124.885.546.14%A754.44
B61.730672 × 1032498.269 × 103447.06 × 103447.57 × 10398.87744.71%A845.91
B71.54 × 1032498.460 × 103415.86 × 103416.17 × 10385.46741.59%A908.54
B81.048050 × 1032498.951 × 103290.03 × 103290.33 × 10384.006.529.00%A916.24
C91.312914 × 1032498.687 × 103367.62 × 103367.87 × 10386.526.536.75%A903.07
C101.251489 × 1032498.748 × 103321.46 × 103321.93 × 10373.171132.15%A980.36
C111.377437 × 1032498.622 × 103485.07 × 103485.23 × 103102.454.448.51%A831.30
C121.179189 × 1032498.820 × 103298.70 × 103299.02 × 10365.86729.87%A1032.24
Table 4. Statistical analysis for correlation coefficient, θd, Vw, and other variables based on Pearson analysis.
Table 4. Statistical analysis for correlation coefficient, θd, Vw, and other variables based on Pearson analysis.
θd (°)Vw (mm3)VT (mm3)Vout (mm3)PA (mm2)SA (mm2)Ψ (°)HP (mm)PB (%)
θd (°)1
Vw (mm3)0.0001
VT (mm3)−0.637 *−0.2601
Vout (mm3)−0.0050.783 **−0.3151
PA (mm2)−0.527−0.2100.844 **−0.4101
SA (mm2)−0.527−0.2100.844 **−0.4111.000 **1
Ψ (°)0.302−0.285−0.078−0.3050.1300.1291
HP (mm)−0.2480.107−0.3080.107−0.418−0.418−0.5651
PB (%)−0.527−0.2100.844 **−0.4101.000 **1.000 **0.129−0.4181
* Correlation is significant at the 0.05 level (2-tailed). ** Correlation is significant at the 0.01 level (2-tailed).
Table 5. Regression goodness of fit for θd, VW with other variables.
Table 5. Regression goodness of fit for θd, VW with other variables.
θd (°) Vw (mm3)
VariableRegression TypeR2 (%)SEquationR2 (%)SEquation
VT (mm3)Linear40.6171.288VT = 2077 − 31.67θd 6.75214.616VT = 1918 − 3.514Vw
Quadratic40.7180.357VT = 1840 − 69θd − 0.620θd2N/AN/AN/A
NonlinearN/A180.835VT = exp(−0.0218546θd)/(0.000449157 + 1.72044 × 10−12θd)N/A1157.87VT = exp(−0.0172781Vw)/(−1.1024 + 0.00918734Vw)
Vout (mm3)Linear0.0222.246Vout = 2,498,562 − 0.27θd138.1861.34Vout = 2,497,126 + 10.59Vw
Quadratic0.01234.261Vout = 2,498,510 + 5.2θd − 0.136θd2N/AN/AN/A
NonlinearN/A234.268Vout = exp(−1.0799 × 10−7θd)/(4.0023 × 10−7 − 3.31555 × 10−16θd)N/A2,040,192Vout = exp(−0.0296687Vw)/(−1.10221 + 0.00918512Vw)
PA (mm2)Linear27.7870,606.7PA = 608,453 − 9791θd4.481,234.5PA = 555,753 − 1060Vw
Quadratic28.4974,057.8PA = 816,737 – 31,525θd + 543θd2N/AN/AN/A
NonlinearN/A74,057.8PA = exp(0.0223781θd)/(3.74916 × 10−7 + 1.7508 × 10−7θd)N/A335,061PA = exp(−0.0277149Vw)/(−1.10245 + 0.00918709Vw)
SA (mm2)Linear27.7770,574.4SA = 608,634 − 9785θd4.4181,187.1SA = 556,277 − 1062Vw
Quadratic28.4874,026.5SA = 816,075 − 31,431θd + 541θd2N/AN/AN/A
NonlinearN/A74,026.5SA = exp(0.0222817θd)/(3.83081 × 10−7 + 1.74142 × 10−7θd)N/A335,257SA = exp(−0.0276883Vw)/(−1.10245 + 0.00918709Vw)
Ψ (°)Linear9.1416.1373Ψ = 64.42 + 1.145θd8.1116.2287Ψ = 126.9 − 0.2934Vw
Quadratic29.9214.9388Ψ = 293.6 − 22.77θd + 0.5978θd2N/AN/AN/A
NonlinearN/A14.9388Ψ = exp(0.134118θd)/(−0.221576 + 0.0205293θd)N/A17.1065Ψ = exp(0.0274097Vw)/(−1.48245 + 0.0147908Vw)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

A.Wahab, M.K.; Mohd Arif Zainol, M.R.R.; Ikhsan, J.; Zawawi, M.H.; Abas, M.A.; Mohamed Noor, N.; Abdul Razak, N.; Sholichin, M. Assessment of Debris Flow Impact Based on Experimental Analysis along a Deposition Area. Sustainability 2023, 15, 13132. https://doi.org/10.3390/su151713132

AMA Style

A.Wahab MK, Mohd Arif Zainol MRR, Ikhsan J, Zawawi MH, Abas MA, Mohamed Noor N, Abdul Razak N, Sholichin M. Assessment of Debris Flow Impact Based on Experimental Analysis along a Deposition Area. Sustainability. 2023; 15(17):13132. https://doi.org/10.3390/su151713132

Chicago/Turabian Style

A.Wahab, Muhammad Khairi, Mohd Remy Rozainy Mohd Arif Zainol, Jazaul Ikhsan, Mohd Hafiz Zawawi, Mohamad Aizat Abas, Norazian Mohamed Noor, Norizham Abdul Razak, and Moh Sholichin. 2023. "Assessment of Debris Flow Impact Based on Experimental Analysis along a Deposition Area" Sustainability 15, no. 17: 13132. https://doi.org/10.3390/su151713132

APA Style

A.Wahab, M. K., Mohd Arif Zainol, M. R. R., Ikhsan, J., Zawawi, M. H., Abas, M. A., Mohamed Noor, N., Abdul Razak, N., & Sholichin, M. (2023). Assessment of Debris Flow Impact Based on Experimental Analysis along a Deposition Area. Sustainability, 15(17), 13132. https://doi.org/10.3390/su151713132

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