Next Article in Journal
SST Anywhere—A Portable Solution for Wide Field Low Earth Orbit Surveillance
Next Article in Special Issue
Gully Morphological Characteristics and Topographic Threshold Determined by UAV in a Small Watershed on the Loess Plateau
Previous Article in Journal
Remote Sensing Estimation and Spatiotemporal Pattern Analysis of Terrestrial Net Ecosystem Productivity in China
Previous Article in Special Issue
Quantitative Evaluation of Gully Erosion Using Multitemporal UAV Data in the Southern Black Soil Region of Northeast China: A Case Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Simple Deposition Model for Debris Flow Simulation Considering the Erosion–Entrainment–Deposition Process

1
Department of Agricultural and Rural Engineering, Chungnam National University, Daejeon 34134, Korea
2
Geologic Environment Division, Korea Institute of Geoscience and Mineral Resources, Daejeon 34132, Korea
3
Department of Regional Infrastructure Engineering, Kangwon National University, Chuncheon 24341, Korea
*
Author to whom correspondence should be addressed.
Submission received: 8 February 2022 / Revised: 1 April 2022 / Accepted: 13 April 2022 / Published: 14 April 2022

Abstract

:
This study aimed to determine the depositional effect and improve the identification of debris flow risk zones. To accomplish this goal, we developed a two-dimensional debris flow model (Deb2D) based on a hyperbolic conservation form of the mass and the momentum balance equation with consideration of the erosion–entrainment effect as well as the depositional effect. In this model, we implemented a widely-used rheological equation—the Voellmy equation—and a quadtree adaptive grid-based shallow-water equation. This model was applied to two study sites to assess the depositional effect. The impact area, volume of soil loss, maximum velocity, inundated depth, and erosion depth resulting from the debris-flow modeling were compared with the field data. The simulation results with/without deposition were evaluated using the receiver operating characteristic method. The implementation results of the erosion–entrainment model with deposition showed superior accuracy when estimating the damage range and flow time.

1. Introduction

Debris flows, defined as gravity-driven sediment mixtures with surface and subsurface flow in mountainous areas, are one of the most serious natural disasters. Additionally, it can be difficult to predict the extent and scale of their resultant damage. In particular, debris flows are capable of rapidly transporting large volumes of sediment and large boulders over long distances, making them destructive and dangerous [1,2] according to local properties such as the geology, topography, and saturated soil by rainfall. Therefore, several studies have attempted to analyze and predict the flow and accumulation process of debris flows using various methods including experimentation [3,4], monitoring [5,6,7,8,9], and numerical modeling [10,11,12,13,14,15]. In particular, many notable numerical studies of debris flow have been performed in recent decades [16,17] and they have provided a better and more detailed understanding of debris flow mechanisms, and their results reflect the suggested theory.
Numerical models are useful for analyzing the flow of difficult-to-implement mixtures such as debris flows. Various phenomena occur during the debris flow process, of which erosion, entrainment, and deposition are typical. Here, the erosion process refers to the removal of the topsoil, and the entrainment refers to the absorption of the topsoil removed through the erosion process. Therefore, erosion–entrainment is an essential factor when predicting and estimating the magnitude of damage because this effect increases the initial volume of the debris flow by 10–50 times [18,19]. In addition, deposition refers to a phenomenon in which flowing sediment movement stops and accumulation begins after reaching a specific condition; deposition is closely linked to the travel distance [2,20]. Therefore, deposition is a significant factor in predicting and calculating the extent of damage caused by debris flows. Therefore, erosion, entrainment, and deposition must all be considered in numerical modeling to accurately analyze and predict the extent and scale of damage caused by debris flows.
Several one-dimensional (1D) and two-dimensional (2D) models have been proposed, and 2D simulations have provided clear insights into the dynamic behavior of debris flows [21,22,23,24,25]. In addition, several erosion–entrainment models have been developed based on different theories and research on the erosion and entrainment phenomena that occur during the debris flow process has been actively conducted in recent decades [10,19,26,27,28,29].
In several studies [30,31,32,33], the erosion–deposition process was implemented using the sediment concentration of debris flows and physical parameters. Pudasaini and Fischer [34] developed the erosion and deposition algorithm through a full physical mechanism based on a two-phase system. However, setting these physical parameters (e.g., grain size, particle diameter, internal friction angle) requires many field surveys. In addition, these characteristics are different for each basal, and it is limited to being defined as a representative value for one basin. Therefore, the applicability of models using physical parameters is still challenging. Medina et al. [10] developed a “stop-and-go” mechanism to estimate damaged areas. However, this mechanism is limited in its ability to implement the deposition process because it simply determines at what point a debris flow will stop, and in addition, its applicability is limited because a detailed explanation of the involved mechanisms was not included. Since the deposition process is relevant in predicting the extent of damage caused by debris flow, several studies on the deposition process have been conducted. However, since most algorithms use physical parameters based on physical mechanics, their usefulness is low for various basal and rainfall situations. Therefore, to predict/respond instantaneously to debris flow that occurs in landslide-prone areas, additional research on the deposition process algorithm based on the conceptual system is needed.
This study developed a simple deposition model for debris flow simulation with a combined erosion–entrainment model. This model refers to the algorithms developed in Medina et al. [10] and Frank et al. [29] that assume that the erosion–entrainment process occurs when the basal shear stress exceeds a specific value, and these algorithms well-simulate the erosion–entrainment process at real-scale debris flow events. Therefore, the developed model simulates the erosion, entrainment, and deposition process based on the algorithms of Medina et al. [10] and Frank et al. [29]. To verify the necessity and importance of the deposition effect, the past two debris flow events (the 2011 Mt. Umyeon debris flow and the 2020 Gokseong-gun debris flow) were selected and analyzed.
Since the field survey data are abundant, the 2011 Mt. Umyeon landslides, which are suitable for the analysis of debris flow, were selected as the study events and the parameters of the deposition algorithm were analyzed. In addition, a recent debris flow event, which occurred in 2020 in Gokseong-gun and killed six people, was simulated to verify the versatility and accuracy of the algorithm developed in this paper. To evaluate the analytical results, the simple method and the receiver operating characteristic (ROC) methods were used to quantitatively analyze the area damaged by the debris flow.

2. Methods

2.1. Governing Equation and Numerical Model

This study used the Deb2D model developed by An et al. [25] to test the deposition model’s performance. The Deb2D model’s governing equation is a hyperbolic conservation form of the mass and the momentum balance equation is expressed as follows:
q t + f x + g y = s ,
where t denotes time; x and y are Cartesian coordinates; and q, f, g, and s are vectors representing conserved variables, fluxes in the x and y directions, and source terms, respectively. The vectors can be written as:
q = ( h h u h v ) ,   f = ( h u h u 2 + g h 2 / 2 h u v ) ,   g = ( h v h u v h v 2 + g h 2 / 2 ) ,   s = ( E D S g x S f x S g y S f y ) ,
where h is the depth of the debris-flow mixture; u and v are the depth-averaged velocity components in the x and y directions, respectively; g is the acceleration due to gravity; ED is the entrainment (ED > 0) and deposition rate (ED < 0); Sgx and Sgy represent the gravitational acceleration in the x and y directions, respectively; and Sfx and Sfy represent the driving friction in the x and y directions, respectively.
The rheological equation is one of the most significant factors and it dominates the behavior of debris flows. Several rheological equations have been used to simulate debris flows and Deb2D can consider the Voellmy, Bingham, and Coulomb-viscous rheological equations. This study used the Voellmy equation because it demonstrated good performance in several previous studies (e.g., [10,35]) while having fewer parameters than other rheological equations. The Voellmy rheology equation is expressed as follows:
S f x = ( μ g h + g | u | 2 ξ ) ,   S f y = ( μ g h + g | v | 2 ξ ) ,
where μ, which represents the Coulomb friction coefficient, dominates the deceleration behavior when the flow is slow and ξ, which represents the turbulent friction coefficient, prevails when the flow is rapid [29,36].
Numerical models based on shallow-water-type governing equations frequently suffer from the imbalance that occurs between the slope of the water depth and the slope on irregular terrain. This imbalance may cause unphysical perturbations and fluctuations in the simulation near shocks or wet–dry transitions. Deb2D implements the hydrostatic reconstruction technique proposed by Audusse et al. [37], which has successfully been applied to the quadtree adaptive grid-based shallow-water equation. For a detailed numerical discretization and scheme, refer to An et al. [25,38]. In addition, Figure 1 depicts a flowchart of the Deb2D model.

2.2. Deposition Model with Erosion–Entrainment Phenomena

Many previous studies have identified the importance of the erosion–entrainment phenomena that occur as debris flows progress [10,25,29,39,40]. The treatment of the erosion–entrainment process in numerical models can be classified into two methods. The first approach assumes that the erosion–entrainment occurs at a constant rate [10,29] and the second approach considers that the erosion-entrainment process’s velocity is determined by the state of the debris flow [10,27,28]. This study used the first approach, which demonstrates the erosion–entrainment process based on physics using basal shear stress. On the other hand, the second approach simulated the erosion–entrainment process using full conceptual/experimental equations that have been limited to demonstrate the erosion–entrainment action reasonably. According to Medina et al.’s [10] static and Frank et al.’s [29] algorithms, the erosion and entrainment process occurred at a constant velocity when the basal shear stress was higher than a critical condition during the flow. This study assumed that the deposition occurred at a constant velocity when the shear stress was lower than a specific condition. The erosion–entrainment–deposition mechanism ED is expressed as follows:
E D ( x , y , t ) = { d z d t e r o s . if   τ > τ e r o s . d z d t d e p o s . else   if   τ < τ d e p o s . ,
h max ( x , y , 0 ) = { d z d τ ( τ τ e r o s . ) = d z d τ ( ρ g h s τ e r o s . ) if   τ > τ e r o s . 0 else ,
where dz/dteros. is a constant erosion–entrainment rate; dz/dtdepos. is a constant deposition rate; τ is the basal shear stress; τeros. and τdepos. are the critical shear stress of erosion and deposition, respectively; hmax(x, y, 0) is the initial height of the erosion–entrainment layer at position (x, y). dz/dτ is the average potential erosion depth; ρ is the mass density; and s is the channel slope. Figure 2 briefly illustrates this algorithm. In addition, this algorithm’s performance, which well simulated the erosion–entrainment–deposition process, has been verified based on laboratory experiments conducted by Lim et al. [41].
The Voellmy rheological equation, which simulates the sediment flow process, uses the parameters calibrated in Lee et al. [40], who analyzed the parameters of the Voellmy equation on the 2011 debris flow in Mt. Umyeon. To determine the importance of deposition, only the deposition algorithm’s parameter was changed to compare the simulation results.

2.3. Analysis Method

Godt et al. [42] and Cepeda et al. [43] used the ROC method to quantitatively analyze the debris flow’s process and damage. Therefore, we used the ROC method to analyze qualitative data including the range of the impact areas from the simulations. In particular, the ROC method’s accuracy factor was used. The accuracy analysis of the ROC is as follows:
Acc . R = D . A + nD . A C . A
where Acc., which is the accuracy, is a value in the range 0–1 and 1 denotes 100% accuracy and 0 denotes 0% accuracy. D.A is the area that successfully implemented the actual damaged area; nD.A denotes the area that was not damaged and also indicates that the simulation results showed that there was no damage; and C.A is the calculation area (Figure 3). For details regarding the ROC method, please refer to Godt et al. [42] and Cepeda et al. [43].
The accuracy of the quantitative results, which include the sediment height, flow velocity, and relocated sediment during the debris flow process, were calculated as follows [40]:
Acc . L = { 0 if   Sim . val . 2 × Obs . val . 1 | Obs . val . Sim . val . | Obs . val . else ,
where Obs.val. and Sim.val. are the observed and simulated values, which include the inundated depth, flow velocity, and volume of eroded sediment. The overall accuracy TotalAcc. was calculated as follows:
TotalAcc . = i = 1 N ( Acc . i × w i ) i = 1 N w i .
where Acc.i is the accuracy calculated using Equations (6) and (7); N is the number of all comparisons; and wi is the weight of each comparison. Here, the weight of each comparison was determined by referring to the reliability of the field data and the coefficient proposed by Cepeda et al. [43]. Therefore, using Equation (8), we attempted to comprehensively evaluate the flow process and impact areas that appeared after the flow.

3. Study Event

Two debris flow events that occurred in the Republic of Korea were selected to validate the deposition model proposed in this study. The first event was the debris flow that occurred in Mt. Umyeon in 2011, as the field data have been investigated extensively by many previous researchers [25,40,44,45,46,47] due to its severe damage and location at the center of the capital city (Figure 4). In particular, field survey data are essential for parameter correction and back-analysis and are also useful when validating a numerical model. The second chosen event was the debris flow that occurred in Gokseong-gun in August 2020. The extent of the damage after the debris flow has been investigated directly and analyzed, so this event was used to test the model performance.

3.1. 2011 Mt. Umyeon Debris Flows

In the Umyeon Mountain area located in Seoul City, 16 shallow landslides and two runout debris flows occurred due to intense rainfall (rainfall amount: 500 mm/day; maximum rainfall intensity: 80 mm/h) on 26 and 27 July 2011 [48]. The bedrock in this area is primarily composed of Precambrian banded biotite gneiss and granitic gneiss and the study area had an average 34° slope [48]. The soil depth within the hillslopes in the study area was in the range of 1–4 m (average 2 m) [48]. When the debris flow occurred due to torrential rainfall, there was significant damage to people and property in the residential area near Mt. Umyeon and the sediment deposited on the nearby road paralyzed traffic.
Debris flows translated by shallow landslides occurred at four points near the top of the mountain within the Raemian catchment (Figure 4a) and at three points in the Sindonga catchment (Figure 4b). The observed channel lengths in the Raemian and Sindonga basins were reportedly 606 m and 664 m, respectively. The volume of soil loss was evaluated by comparing the digital elevation models using light detection and ranging (LiDAR DEMs) before and after the debris flow events. The soil loss volumes were 25,940 m3 and 21,070 m3 in the Raemian and Sindonga basins, respectively. The maximum velocity of the debris flows, which were analyzed using CCTV and a car black box, was estimated to be approximately 28 m/s and 18 m/s at the Raemian and Sindonga apartment blocks, respectively. It was confirmed that the Raemian and Sindonga apartments were directly damaged to the 3rd floor (approximately 10 m) and the 2nd floor (8 m), respectively, due to the debris flows.
In this study, a LiDAR DEM (1 m × 1 m) was used to input data for the topography. The scale of the collapse spots was determined based on the difference in the LiDAR DEMs before and after the landslide. Therefore, the initially measured landslide volumes were 350 m3 for both catchments.
Using Equations (6)–(8), the results were compared and analyzed using field data. The criteria were (1) the impact area; (2) the volume of soil loss; (3) the inundated depth observed near the apartment blocks; and (4) the debris flows’ maximum velocity.

3.2. 2020 Gokseong-Gun Debris Flows

In the Seongdeok reservoir, Gokseong-gun, which is located in South Jeolla Province, debris flows recently occurred due to intense rainfall (rainfall amount: 277 mm/3 days; maximum rainfall intensity: 52 mm/h) 5–7 August 2020 (Figure 4c). The length of the Seongdeok basin is 678 m, which is similar to the two Mt. Umyeon basins. The debris flow caused by the torrential rainfall was rapid, whereas the flow velocity of the debris flow was slower than that of the 2011 Mt. Umyeon debris flow. Nevertheless, the volume of the slope collapse that occurred near the road was quite large and it caused considerable damage to rural areas and rice paddies in the lower part of the watershed, and five fatalities occurred [49].
According to the on-site survey, a large-scale collapse occurred on the slope facing the road and the collapsed area had an average slope of 35°. According to Choi et al. [49] and the field survey, the erosion depth was <2.5 m and the volume of soil loss was evaluated at ~30,000 m3. There is a river in the lower region, therefore, some debris flows into the river. However, since this study simulated the debris flow based on a one-phase/-layer, the implementation for such inflow into the river was limited. The principal flat bottom valleys in Gokseong-gun are composed of deeply weathered soils [50]. The characteristic weathering profile of these granite and binary rocks are gravel and laminated gravel tens of meters thick with irregularly oriented crack areas under the loose weathering material [49]. A DEM (5 m × 5 m) was used as the input data for the topography. The data on the collapsed area were constructed by referring to the digital surface model (DSM) built by drones after the debris flow event and the volume was set to approximately 14,500 m2. Since this basin lacks field data, precise comparative analysis is limited so we attempted to confirm the implementation rate of the impact area using Equation (6).

4. Results

4.1. Simulation Results: Mt. Umyeon Area

Table 1 summarizes the parameters used in the catchment of the Raemian and Sindonga apartments. The Voellmy rheological equation and the erosion–entrainment–deposition model’s parameters were calibrated for each basin. In addition, to reveal the importance of deposition, the parameter used in the deposition algorithm was the only one we changed (dz/dtdepos. = 0 or 0.01).
Figure 5 presents the simulated maximum flow height of the debris flow according to the algorithm with/without the deposition process in both catchments. As this figure shows, there is a clear difference in the extent of the damage depending on the deposition. When the deposition process was considered, the impact area observed through the field surveys was adequately demonstrated in both catchments (Figure 5a,c). However, when the deposition process was not considered, the simulated damage caused by the debris flow occurred over a wider area because the impact area was overestimated; it appeared that it went outside of the calculation range set in the numerical model (Figure 5b,d).
To precisely analyze this phenomenon, Figure 6 and Figure 7 show the height of the debris flow over time in the Raemian apartment basin. Figure 6 shows the results when the deposition was considered; when the flow time reached 2 min, most of the debris flow had stopped and was deposited. The damage range did not spread any farther. However, according to the simulation results in Figure 7, which did not consider the deposition process, the flow of sediment did not stop over time and the impact area was unphysically expanded. In addition, we attempted to determine the effect of the parameter change in the deposition algorithm on the estimation of the damaged area (Figure 8). A larger dz/dtdepos. meant that the deposition process was larger and the simulated impact area was very small. Conversely, when the parameter was closer to 0, the impact area was wider.
Figure 9 shows the simulation results of the erosion depth for the two catchments. Figure 9a,d depict the observations for each catchment, which were constructed using LiDAR DEMs before and after the occurrence of the debris flow. According to the results, in both basins, there was no critical difference with or without the deposition process. However, when the deposition was accounted for, as shown in Figure 9b,e, the phenomenon of widespread erosion that was demonstrated in Figure 9c,f was reduced, which is contrast from the results from the field data. In addition, it showed that there was a considerable difference in the volume of soil loss according to the with/without deposition process. Therefore, this result revealed that the erosion depth, which was equivalent to the topography change, was simulated as relatively similar to the results shown in Figure 9b,e in which the deposition process was considered.
Table 2 summarizes the simulation results for the debris flow in the two catchments. In the Raemian apartment basin, the volume of soil loss analyzed with/without the deposition process was 25,260 m3 and 47,390 m3, respectively. When deposition was considered, the observed value of 25,940 m3 was well-simulated. The amount of soil loss analyzed with/without deposition in the Sindonga apartment basin was 21,900 m3 and 48,170 m3, respectively. In addition, when the deposition process was accounted for, the observation of 21,070 m3 was well-simulated. As such, the analysis results for the volume of soil loss revealed large differences depending on the deposition process. However, the analysis results for the maximum damage height of the debris flow observed near the Raemian and Sindonga apartments and the maximum flow velocity of the debris flow observed near the roads were similar, regardless of the deposition process. In the case of the Raemian apartment basin, the inundated depths analyzed with/without deposition were 10 m and 11 m, respectively, and they matched the observed values of 9–11 m. In the Sindonga apartment basin, the inundated depth was implemented to be 6 m, regardless of the deposition process and this value was consistent with the observation of 6–8 m. In addition, the maximum flow velocity was analyzed to be 16 m/s, regardless of the deposition process, which was consistent with the observed value of 18 m/s. However, in the Raemian apartment basin, the maximum flow velocity was analyzed to be 19 m/s regardless of the deposition process, which was slow compared to the observation of 28 m/s.
Table 3 summarizes the implementation rate results of the flow process and impacted area analyzed using Equations (6)–(8). As a result, when the deposition was considered at the Raemian apartment basin in Table 3, TotalAcc. was 0.914, and when the deposition was not considered, it was analyzed to be 0.779. At the Sindonga apartment basin, when the deposition was considered, TotalAcc. was 0.917; when the deposition was not considered, it was analyzed as 0.753. Overall, the simulation accuracy was 0.13–0.17 higher when considering the deposition process. In particular, considering the deposition, impact area and the amount of soil loss resulted in a more accurate prediction. Moreover, the simulation without deposition demonstrated the worst performance in terms of “volume of soil loss”.

4.2. Simulation Results: Seongdeok Area

Table 1 summarizes the parameters used to analyze the Seongdeok basin. The Voellmy rheological equation and the erosion–entrainment–deposition model’s parameters were also calibrated for this basin based on field survey data. The Voellmy rheology parameters that control the fluidity of debris flow were calibrated considering that the rainfall intensity and total rainfall were less than that for the previous two basins. Similar to the Mt. Umyeon landslides, they were implemented with/without the deposition process (dz/dtdepos. = 0 or 0.01). In the case of the Seongdeok basin, the simulated results and the observation data that could be compared and analyzed included only the damaged range. Therefore, using Equation (6), we assessed the implementation rate (accuracy) of the damage range that was analyzed using the algorithm developed in this study.
Figure 10 shows the damaged area observed by drones and on-site surveys as well as the simulation results for the maximum flow depth of the debris flow after 5 min, and the river was briefly displayed. As shown in Figure 10a, the impact area was well-depicted when the deposition was considered. However, when the deposition process was not considered, the damaged area was much wider than what was observed, as shown in Figure 10b. In the numerical model, it was difficult to implement the effect of the river and rice fields. Therefore, regardless of the deposition process, the damage that occurred across the river was not observed in the field. This phenomenon was demonstrated precisely in Figure 11, which analyzed the debris flow height according to time. As shown in Figure 11, the river blocked some sediments. However, the system by which the sediments flowed into the river was not simulated, and it was difficult to implement the sediment flow into the rice fields that could have reduced the damage of the debris flow. However, it was sufficient to simulate and analyze the impact area of the debris flow that occurred in the Seongdeok basin.
According to the field survey and Choi et al. [49], the amount of soil loss was analyzed to be about ~30,000 m3, which is the volume calculated by considering both erosion and deposition. For the simulation result in the Seongdeok basin, when the deposition effect was considered, the analyzed soil loss was 29,600 m3, which is consistent with the volume calculated in the field survey and previous studies. However, when deposition was not considered, the amount of soil loss was analyzed to be 41,980 m3, which is an overestimate.
Table 4 summarizes the simulation results of the debris flow. When simulating with deposition, the impact area was analyzed with an accuracy of 0.952. The simulation without deposition demonstrated the impact area with high accuracy of 0.909. However, the volume of soil loss in the basal was overestimated considerably. In addition, the maximum flow velocity observed in the residential area (rural area) was analyzed to be 10 m/s, and in the simulation, the debris flow affected the village up to 6 m. These results revealed a lower value than that of the Mt. Umyeon debris flow event, but it was estimated to be sufficient to cause damage to the village.

5. Discussion

Figure 5 shows the importance of the deposition process. When deposition was considered, the extent of the damage was simulated similarly to the field survey, as shown in Figure 5a,c. However, when the deposition process was not considered, as shown in Figure 5b,d, the impact area was overestimated compared to the result from the field survey. To analyze this phenomenon more precisely, the flow height of the debris flow over time was analyzed for the Raemian apartment basin, as shown in Figure 6 and Figure 7. According to these simulation results, the deposition process stopped the flow of sediments under certain conditions, as shown in Figure 6, and caused them to be deposited; additionally, the analysis that accounted for deposition helped to estimate the flow time of the debris flow. However, when simulating the debris flow event without the deposition process, the damaged area was extended unnaturally over time, as shown in Figure 7. Therefore, to accurately estimate the area impacted by the debris flow, Figure 5 and Figure 6 demonstrate that it is essential to consider the deposition process. However, as shown in Figure 8, it revealed that parameter calibration was necessary for the deposition algorithm, as in the Voellmy rheological equation and the erosion–entrainment model. This phenomenon, which stopped the flow of sediments, was implemented in Medina et al. [10] using the “stop-and-go” system. Since the “stop-and-go” system controls the flow and stoppage of debris flow, this system made it possible to precisely analyze the damaged area in Median et al.’s [10] study. However, a detailed explanation of the “stop-and-go” mechanism was omitted and it played a strictly different role from that in the deposition process, so its applicability was poorer than the deposition algorithm introduced in this paper.
Figure 9 shows the simulation results of the erosion depths in both catchments. Visually, there was no critical difference according to the presence/absence of the deposition process. However, when analyzed with the “volume of soil loss” in Table 2, it was revealed that the deposition process was necessary to accurately calculate the volume of soil loss and topography changes. Therefore, debris flow analysis using the numerical model and considering the deposition process would help not only to analyze the impact area, but also to accurately estimate the flow time and topography change.
As shown in Table 2 and Table 3, the debris flow simulations for the Raemian and Sindonga apartment basins were successfully carried out. In addition, regardless of whether deposition was considered, the inundated depth and maximum velocity near both apartments were similarly reproduced. However, it was difficult to implement the flow velocity of the debris flow, and in particular, the sediment flowed at 19 m/s in the Raemian apartment basin, which was very low compared to the observed value of 28 m/s. Its rapid velocity was due to the fact that the debris flow event at Mt. Umyeon in 2011 contained considerably more water than the others [40,44,46]. However, the flow velocity in the Sindonga apartment basin was well-implemented. Therefore, it revealed that the Voellmy rheological equation struggles to demonstrate a flow velocity of ≥20 m/s [40].
In addition, the debris flow simulation for the Seongdeok basin was successfully carried out. The soil loss and impact area were well-implemented through the deposition process. Figure 10 also demonstrates the importance of the deposition process in the analysis of the debris flow. In particular, the erosion–entrainment–deposition model had a considerable influence on the modeling of the deposition that occurred in the upper middle region of the basin as well as deposition in the lower region of the basin, which is a residential area. However, the damage across the river, which did not appear in the actual event, was revealed by the simulation (Figure 10 and Figure 11). This result was assessed to be an error caused by the fact that the DEM data are low-resolution to implement the curvature of terrain and was also caused by a limitation in the numerical model in which the sediments flowing into the river could not be implemented. However, using the simulation results in Table 4, it was estimated that the approximate damage scale and damage range were well-implemented.
In this paper, the deposition process was implemented to simulate the debris flow with the erosion and entrainment processes. In addition, deposition plays an important role in the analysis of debris flows as well as erosion and entrainment. However, the current body of research regarding deposition is relatively incomplete. Therefore, studies that explore deposition need to make better predictions of debris flows. In addition, the analysis method in this paper is not practical for unmeasured catchments because data regarding the collapse point is essential for the analysis of debris flows. Therefore, in future research, we plan to comprehensively implement slope stability and debris flow analysis to demonstrate and predict debris flow events that occur in unmeasured areas and areas with a high risk of landslides and debris flows.

6. Conclusions

In this study, we attempted to analyze the importance of depositional processes on the flow and accumulation process of debris flows. Among the debris flow events that occurred in Mt. Umyeon, Seoul in 2011, the Raemian and Sindonga apartment basins, and the Gokseong-gun, South Jeolla Prefecture in 2020, the Seongdeok basin was selected as the study area and the debris flows were well-implemented. To analyze the simulation results such as the flow velocity, amount of soil loss, and inundated depth, the ROC method was used. Implementing the debris flows while considering the erosion–entrainment–deposition process was successful and the importance of the deposition process was revealed through the simulation results.
In particular, deposition played a significant role in estimating the range of the damaged area caused by debris flows. In addition, the debris flows over time revealed that the deposition process helped assess the flow time. Debris flow analyses that consider erosion, entrainment, and deposition should help analyze the topography change more precisely. The overall analysis results demonstrated an approximately 0.13–0.17 (13–17%) higher accuracy when considering deposition. According to this result, analyses of debris flow events that account for the erosion–entrainment–deposition process could calculate the impact area, flow time, and topography changes. In other words, the deposition process helps simulate the debris flow so that the results are closer to the natural phenomenon along with the erosion–entrainment process.
In the Seongdeok basin, it was not possible to accurately implement the debris flow event due to the lack of on-site investigation. Therefore, it was analyzed and verified based on the damaged area captured by drones after the debris flow. In the simulation results, we attempted to precisely analyze the effect of reducing the impact area due to the rice fields and river in the field survey, but this effort was limited to implementing the rice fields and river due to the limitations of the numerical model. However, the debris flow that occurred in the Seongdeok basin was successfully simulated. In addition, in the analysis that considered the deposition process, the impact area was simulated with approximately 0.05 (5%) higher accuracy. However, as observed in the two Mt. Umyeon basins, the deposition did not affect the flow height or flow velocity.

Author Contributions

Conceptualization, S.L. and H.A.; Methodology, S.L. and H.A.; Software, H.A. and M.K.; Validation, S.L., H.L. and Y.K.; Formal analysis, S.L. and H.A.; Investigation, S.L., H.A. and M.K.; Data curation, M.K.; Writing—original draft preparation, S.L.; Writing—review and editing, H.A. and M.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Korea Ministry of Environment as “The SS projects; 2019002830001” and the Basic Research Project of the Korea Institute of Geoscience and Mineral Resources (Project code: 22-3412-1).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Coe, J.A.; Kean, J.W.; Godt, J.W.; Baum, R.L.; Jones, E.S.; Gochis, D.J.; Anderson, G.S. New insights into debris-flow hazards from an extraordinary event in the Colorado Front Range. GSA Today 2014, 24, 4–10. [Google Scholar] [CrossRef]
  2. Takahashi, T. Debris Flow: Mechanics, Prediction and Countermeasures; Taylor & Francis: Abingdon, UK, 2007. [Google Scholar]
  3. Iverson, R.M.; Reid, M.E.; Logan, M.; LaHusen, R.G.; Godt, J.W.; Griswold, J.P. Positive feedback and momentum growth during debris-flow entrainment of wet bed sediment. Nat. Geosci. 2011, 4, 116–121. [Google Scholar] [CrossRef]
  4. Zhou, G.G.; Li, S.; Song, D.; Choi, C.E.; Chen, X. Depositional mechanisms and morphology of debris flow: Physical modelling. Landslides 2019, 16, 315–332. [Google Scholar] [CrossRef]
  5. Oguchi, T.; Oguchi, C.T. Late Quaternary rapid talus dissection and debris flow deposition on an alluvial fan in Syria. Catena 2004, 55, 125–140. [Google Scholar] [CrossRef]
  6. Blasone, G.; Cavalli, M.; Marchi, L.; Cazorzi, F. Monitoring sediment source areas in a debris-flow catchment using terrestrial laser scanning. Catena 2014, 123, 23–36. [Google Scholar] [CrossRef]
  7. Hürlimann, M.; Abancó, C.; Moya, J.; Vilajosana, I. Results and experiences gathered at the Rebaixader debris-flow monitoring site, Central Pyrenees, Spain. Landslides 2014, 11, 939–953. [Google Scholar] [CrossRef]
  8. Gao, Y.; Yin, Y.; Li, B.; Feng, Z.; Wang, W.; Zhang, N.; Xing, A. Characteristics and numerical runout modeling of the heavy rainfall-induced catastrophic landslide–debris flow at Sanxicun, Dujiangyan, China, following the Wenchuan Ms 8.0 earthquake. Landslides 2017, 14, 1361–1374. [Google Scholar] [CrossRef]
  9. Lyu, L.; Xu, M.; Wang, Z.; Cui, Y.; Blanckaert, K. A field investigation on debris flows in the incised Tongde sedimentary basin on the northeastern edge of the Tibetan Plateau. Catena 2022, 208, 105727. [Google Scholar] [CrossRef]
  10. Medina, V.; Hürlimann, M.; Bateman, A. Application of FLATModel, a 2D finite volume code, to debris flows in the northeastern part of the Iberian Peninsula. Landslides 2008, 5, 127–142. [Google Scholar] [CrossRef]
  11. Gregoretti, C.; Degetto, M.; Boreggio, M. GIS-based cell model for simulating debris flow runout on a fan. J. Hydrol. 2016, 534, 326–340. [Google Scholar] [CrossRef]
  12. Han, Z.; Su, B.; Li, Y.; Wang, W.; Wang, W.; Huang, J.; Chen, G. Numerical simulation of debris-flow behavior based on the SPH method incorporating the Herschel-Bulkley-Papanastasiou rheology model. Eng. Geol. 2019, 255, 26–36. [Google Scholar] [CrossRef]
  13. Ouyang, C.; Wang, Z.; An, H.; Liu, X.; Wang, D. An example of a hazard and risk assessment for debris flows—A case study of Niwan Gully, Wudu, China. Eng. Geol. 2019, 263, 105351. [Google Scholar] [CrossRef]
  14. Chang, M.; Liu, Y.; Zhou, C.; Che, H. Hazard assessment of a catastrophic mine waste debris flow of Hou Gully, Shimian, China. Eng. Geol. 2020, 275, 105733. [Google Scholar] [CrossRef]
  15. Liu, W.; He, S. Comprehensive modelling of runoff-generated debris flow from formation to propagation in a catchment. Landslides 2020, 17, 1529–1544. [Google Scholar] [CrossRef]
  16. Hungr, O.; Corominas, J.; Eberhardt, E. Estimating landslide motion mechanism, travel distance and velocity. In Landslide Risk Management; CRC Press: Boca Raton, FL, USA, 2005; pp. 109–138. [Google Scholar]
  17. Mangeney, A.; Bouchut, F.; Thomas, N.; Vilotte, J.P.; Bristeau, M.O. Numerical modeling of self-channeling granular flows and of their levee-channel deposits. J. Geophys. Res. Earth Surf. 2007, 112. [Google Scholar] [CrossRef] [Green Version]
  18. Remaître, A.; Malet, J.P.; Maquaire, O. Morphology and sedimentology of a complex debris flow in a clay-shale basin. Earth Surf. Process. Landf. J. Br. Geomorphol. Res. Gr. 2005, 30, 339–348. [Google Scholar] [CrossRef]
  19. Pudasaini, S.P.; Krautblatter, M. The Mechanics of Landslide Mobility with Erosion. Nat. Commun. 2021, 12, 6793. [Google Scholar] [CrossRef]
  20. Schuerch, P. Debris-Flow Erosion and Deposition Dynamics; Durham University: Durham, UK, 2011. [Google Scholar]
  21. Laigle, D.; Coussot, P. Numerical modeling of mudflows. J. Hydraul. Eng. 1997, 123, 617–623. [Google Scholar] [CrossRef]
  22. Denlinger, R.P.; Iverson, R.M. Flow of variably fluidized granular masses across three-dimensional terrain: 2. Numerical predictions and experimental tests. J. Geophys. Res. Solid Earth 2001, 106, 553–566. [Google Scholar] [CrossRef]
  23. Chen, H.; Lee, C.F. A dynamic model for rainfall-induced landslides on natural slopes. Geomorphology 2003, 51, 269–288. [Google Scholar] [CrossRef]
  24. Rickenmann, D.; Laigle, D.; McArdell, B.; Hübl, J. Comparison of 2D debris-flow simulation models with field events. Comput. Geosci. 2006, 10, 241–264. [Google Scholar] [CrossRef] [Green Version]
  25. An, H.; Kim, M.; Lee, G.; Kim, Y.; Lim, H. Estimation of the area of sediment deposition by debris flow using a physical-based modeling approach. Quat. Int. 2019, 503, 59–69. [Google Scholar] [CrossRef]
  26. Takahashi, T.; Nakagawa, H. Prediction of stony debris flow induced by severe rainfall. J. Jpn. Soc. Eros. Control Eng. 1991, 44, 12–19. [Google Scholar] [CrossRef]
  27. McDougall, S.; Hungr, O. Dynamic modelling of entrainment in rapid landslides. Can. Geotech. J. 2005, 42, 1437–1448. [Google Scholar] [CrossRef]
  28. Sovilla, B.; Burlando, P.; Bartelt, P. Field experiments and numerical modeling of mass entrainment in snow avalanches. J. Geophys. Res. Earth Surf. 2006, 111. [Google Scholar] [CrossRef]
  29. Frank, F.; McArdell, B.W.; Huggel, C.; Vieli, A. The importance of entrainment and bulking on debris flow runout modeling: Examples from the Swiss Alps. Nat. Hazards Earth Syst. Sci. 2015, 15, 2569–2583. [Google Scholar] [CrossRef] [Green Version]
  30. Takahashi, T. Dynamics of Debris Flows in the Inertial Regime. In Debris-Flow Hazards Mitigation: Mechanics, Prediction, and Assessment; ASCE: Reston, VA, USA, 1997; pp. 239–248. [Google Scholar]
  31. Egashira, S. Constitutive Equations of Debris Flow and Their Applicability. In Proceedings of the 1st International Conference Water Resources Engineering Division/ASCE at San Francisco, San Antonio, TX, USA, 14–18 August 1995; ASCE: Reston, VA, USA, 1997; pp. 340–349. [Google Scholar]
  32. Satofuka, Y.; Mizuyama, T. Numerical simulation on debris flow control by a grid dam. J. Jpn. Soc. Eros. Control Eng. 2005, 57, 21–27. [Google Scholar]
  33. Shrestha, B.B.; Nakagawa, H.; Kawaike, K.; Baba, Y. Numerical simulation on debris-flow deposition and erosion processes upstream of a check dam with experimental verification. Disaster Prev. Res. Inst. Annu. 2008, 51, 613–624. [Google Scholar]
  34. Pudasaini, S.P.; Fischer, J.-T. A mechanical erosion model for two-phase mass flows. Int. J. Multiph. Flows 2020, 132, 103416. [Google Scholar] [CrossRef]
  35. Naef, D.; Rickenmann, D.; Rutschmann, P.; McArdell, B. Comparison of flow resistance relations for debris flows using a one-dimensional finite element simulation model. Nat. Hazards Earth Syst. Sci. 2006, 6, 155–165. [Google Scholar] [CrossRef] [Green Version]
  36. Bartelt, P.; Buehler, Y.; Christen, M.; Deubelbeiss, Y.; Graf, C.; McArdell, B. RAMMS–Rapid Mass Movement Simulation, A Modeling System for Debris Flows in Research and Practice, User Manual v1.5, Debris Flow, Manuscript Update: 31 January 2013; WSL Institute for Snow and Avalanche Research SLF: Davos, Switzerland, 2013; Available online: http://ramms.slf.ch/ramms/downloads/RAMMS_DBF_Manual.pdf (accessed on 27 February 2015).
  37. Audusse, E.; Bouchut, F.; Bristeau, M.O.; Klein, R.; Perthame, B.T. A Fast and Stable Well-Balanced Scheme with Hydrostatic Reconstruction for Shallow Water Flows. SIAM J. Sci. Comput. 2004, 25, 2050–2065. [Google Scholar] [CrossRef]
  38. An, H.; Ichikawa, Y.; Tachikawa, Y.; Shiiba, M. Comparison between iteration schemes for three-dimensional coordinate-transformed saturated–unsaturated flow model. J. Hydrol. 2012, 470–471, 212–226. [Google Scholar] [CrossRef]
  39. Shen, W.; Wang, D.; Qu, H.; Li, T. The effect of check dams on the dynamic and bed entrainment processes of debris flows. Landslides 2019, 16, 2201–2217. [Google Scholar] [CrossRef]
  40. Lee, S.; An, H.; Kim, M.; Lim, H. Analysis of debris flow simulation parameters with entrainment effect: A case study in the Mt. Umyeon. J. Korea Water Resour. Assoc. 2020, 53, 637–646. [Google Scholar] [CrossRef]
  41. Lim, H. Erosion-Deposition Characteristics of Debris Flow Using Laboratory Experiment and Numerical Analysis; Kangwon National University: Chuncheon, Korea, 2021. [Google Scholar]
  42. Godt, J.; Baum, R.; Savage, W.; Salciarini, D.; Schulz, W.; Harp, E. Transient deterministic shallow landslide modeling: Requirements for susceptibility and hazard assessments in a GIS framework. Eng. Geol. 2008, 102, 214–226. [Google Scholar] [CrossRef]
  43. Cepeda, J.; Chávez, J.A.; Cruz Martínez, C. Procedure for the selection of runout model parameters from landslide back-analyses: Application to the Metropolitan Area of San Salvador, El Salvador. Landslides 2010, 7, 105–116. [Google Scholar] [CrossRef]
  44. Kim, S.; Paik, J.; Kim, K.S. Run-out modeling of debris flows in Mt. Umyeon using FLO-2D. J. Korean Soc. Civ. Eng. 2013, 33, 965–974. [Google Scholar] [CrossRef] [Green Version]
  45. Park, D.W.; Nikhil, N.; Lee, S. Landslide and debris flow susceptibility zonation using TRIGRS for the 2011 Seoul landslide event. Nat. Hazards Earth Syst. Sci. 2013, 13, 2833–2849. [Google Scholar] [CrossRef] [Green Version]
  46. Lee, D.H.; Lee, S.R.; Park, J.Y. Numerical Simulation of Debris Flow Behavior at Mt. Umyeon using the DAN3D Model. J. Korean Soc. Hazard. Mitig. 2019, 19, 195–202. [Google Scholar] [CrossRef]
  47. Hong, M.; Jeong, S.; Kim, J. A combined method for modeling the triggering and propagation of debris flows. Landslides 2020, 17, 805–824. [Google Scholar] [CrossRef]
  48. Korean Geotechnical Society. Research Contract Report: Addition and Complement Causes Survey of Mt. Woomyeon Landslide; Research Report, 51-6110000-000649-01; Korean Geotechnical Society: Seoul, Korea, 2014; Volume 2011. [Google Scholar]
  49. Choi, S.K.; Ramirez, R.A.; Kwon, T.H. Preliminary report of a catastrophic landslide that occurred in Gokseong County, South Jeolla Province, South Korea, on 7 August 2020. Landslides 2021, 18, 1465–1469. [Google Scholar] [CrossRef]
  50. Cho, M.; Choi, Y.; Ha, K.; Kee, W.; Lachassagne, P.; Wyns, R. Paleoweathering covers in Korean hard rocks: A methodology for mapping their spatial distribution and the thickness of their constituting horizons. Applications to identify brittle deformation and to hard rock hydrogeology. KIGAM Bull. 2002, 6, 12–25. [Google Scholar]
Figure 1. Flowchart of the Deb2D model (modified from An et al. [25]).
Figure 1. Flowchart of the Deb2D model (modified from An et al. [25]).
Remotesensing 14 01904 g001
Figure 2. Processes of (a) erosion–entrainment and (b) deposition.
Figure 2. Processes of (a) erosion–entrainment and (b) deposition.
Remotesensing 14 01904 g002
Figure 3. Definition of variables for the estimation of discrete classifiers for spatially distributed values by the ROC method (modified from Godt et al. [42] and Cepeda et al. [43]).
Figure 3. Definition of variables for the estimation of discrete classifiers for spatially distributed values by the ROC method (modified from Godt et al. [42] and Cepeda et al. [43]).
Remotesensing 14 01904 g003
Figure 4. The study event location: the debris flows at Mt. Umyeon in 2011, (a) Raemian and (b) Sindonga apartment basins, and (c) Gokseong-gun, Seongdeok basin in 2020.
Figure 4. The study event location: the debris flows at Mt. Umyeon in 2011, (a) Raemian and (b) Sindonga apartment basins, and (c) Gokseong-gun, Seongdeok basin in 2020.
Remotesensing 14 01904 g004
Figure 5. The maximum depth of the Mt. Umyeon debris flows after 5 min on the Raemian apartment basin simulation (a) with deposition and (b) without deposition; and the Sindonga apartment basin simulation (c) with deposition and (d) without deposition.
Figure 5. The maximum depth of the Mt. Umyeon debris flows after 5 min on the Raemian apartment basin simulation (a) with deposition and (b) without deposition; and the Sindonga apartment basin simulation (c) with deposition and (d) without deposition.
Remotesensing 14 01904 g005
Figure 6. Simulated flow depth at the Raemian apartment basin over time with deposition where dz/dtdepos. = 0.01.
Figure 6. Simulated flow depth at the Raemian apartment basin over time with deposition where dz/dtdepos. = 0.01.
Remotesensing 14 01904 g006
Figure 7. Simulated flow depth at the Raemian apartment basin over time without deposition where dz/dtdepos. = 0.
Figure 7. Simulated flow depth at the Raemian apartment basin over time without deposition where dz/dtdepos. = 0.
Remotesensing 14 01904 g007
Figure 8. Spreadability of debris flow according to the change in parameter dz/dtdepos. at the Raemian apartment basin.
Figure 8. Spreadability of debris flow according to the change in parameter dz/dtdepos. at the Raemian apartment basin.
Remotesensing 14 01904 g008
Figure 9. Erosion depth of the Mt. Umyeon debris flows after 5 min in the Raemian apartment basin (a) observed, (b) with deposition, where dz/dtdepos. = 0.01 and (c) without deposition, where dz/dtdepos. = 0; and the Sindonga apartment basin (d) observed, (e) with deposition, where dz/dtdepos. = 0.01 and (f) without deposition, where dz/dtdepos. = 0.
Figure 9. Erosion depth of the Mt. Umyeon debris flows after 5 min in the Raemian apartment basin (a) observed, (b) with deposition, where dz/dtdepos. = 0.01 and (c) without deposition, where dz/dtdepos. = 0; and the Sindonga apartment basin (d) observed, (e) with deposition, where dz/dtdepos. = 0.01 and (f) without deposition, where dz/dtdepos. = 0.
Remotesensing 14 01904 g009
Figure 10. Maximum depth of the debris flow at Seongdeok basin after 5 min (a) with deposition and (b) without deposition.
Figure 10. Maximum depth of the debris flow at Seongdeok basin after 5 min (a) with deposition and (b) without deposition.
Remotesensing 14 01904 g010
Figure 11. Simulated flow depth at the Seongdeok basin over time with deposition.
Figure 11. Simulated flow depth at the Seongdeok basin over time with deposition.
Remotesensing 14 01904 g011
Table 1. Calibrated parameters for the Deb2D debris-flow model (Voellmy rheological equation with erosion–entrainment algorithm considering deposition) at the Mt. Umyeon and Gokseung-gun debris flow events.
Table 1. Calibrated parameters for the Deb2D debris-flow model (Voellmy rheological equation with erosion–entrainment algorithm considering deposition) at the Mt. Umyeon and Gokseung-gun debris flow events.
Study EventVoellmy RheologyErosion–EntrainmentDeposition
μξ (m/s2)dz/dteros. (m/s)τeros. (kPa)dz/ (m/kPa)ρ (kg/m3)dz/dtdepos. (m/s)
Raemian apartment basin0.0430000.081.00.319000.01
Sindonga apartment basin0.0430000.071.50.319000.01
Seongdeok basin0.0520000.071.00.219000.01
Table 2. 2011 debris flow simulation results for the Raemian and Sindonga apartment basins with/without deposition process.
Table 2. 2011 debris flow simulation results for the Raemian and Sindonga apartment basins with/without deposition process.
Study EventWith/Without DepositionVolume of Soil Loss (m3)Inundated Depth (m)Maximum Velocity (m/s)
Raemian apartment basinwith-25,2601019
without-47,3901119
observation25,9409–1128
Sindonga apartment basinwith-21,900616
without-48,170616
observation21,0706–818
Table 3. Analysis results reported by the simple and ROC methods when analyzing the debris flows in the 2011 Raemian and Sindonga apartment basin simulations according to with/without deposition action.
Table 3. Analysis results reported by the simple and ROC methods when analyzing the debris flows in the 2011 Raemian and Sindonga apartment basin simulations according to with/without deposition action.
Study EventWith/Without
Deposition
Acc.Total Acc.
Impact AreaVolume of Soil LossInundated DepthMaximum Velocity
wi3121
Raemian apartment basinwith-0.9140.9741.0000.6790.914
without-0.8680.1731.0000.6790.779
Sindonga apartment basinwith-0.8570.9611.0000.8890.917
without-0.7940.0001.0000.8890.753
Table 4. 2020 Seongdeok basin debris flow simulation results with deposition action.
Table 4. 2020 Seongdeok basin debris flow simulation results with deposition action.
Study EventWith/Without DepositionVolume of Soil Loss (m3)Inundated Depth(m)Maximum Velocity (m/s)Impact Area (Acc.)
Seongdeok basinwith-29,6006100.952
without-41,9806100.909
observation~30,000UnmeasuredUnmeasuredField survey
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lee, S.; An, H.; Kim, M.; Lim, H.; Kim, Y. A Simple Deposition Model for Debris Flow Simulation Considering the Erosion–Entrainment–Deposition Process. Remote Sens. 2022, 14, 1904. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14081904

AMA Style

Lee S, An H, Kim M, Lim H, Kim Y. A Simple Deposition Model for Debris Flow Simulation Considering the Erosion–Entrainment–Deposition Process. Remote Sensing. 2022; 14(8):1904. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14081904

Chicago/Turabian Style

Lee, Seungjun, Hyunuk An, Minseok Kim, Hyuntaek Lim, and Yongseong Kim. 2022. "A Simple Deposition Model for Debris Flow Simulation Considering the Erosion–Entrainment–Deposition Process" Remote Sensing 14, no. 8: 1904. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14081904

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