Comment on egusphere-2022-1398

Debris flows are one of the perilous landslide-related hazards due to their fast flow velocity, large impact force, and long runout, in association with poor predictability. Debris-flow barriers that can minimize the energy of debris flows have been widely constructed to mitigate potential damages. However, the interactions between debris flows and barriers remain poorly understood, which hampers the optimal barrier installation against debris flows. Therefore, this study examined the effect of barrier locations, in particular source-to-barrier distance, on velocity and volume of debris flows via the numerical approach based on smoothed particle hydrodynamics (SPH). A debris-flow event was simulated on a 3D terrain, in which a closed-type barrier was numerically created at predetermined locations along a debris-flow channel, varying the source-to-barrier distance from the initiation point. In all cases, the closed-type barrier significantly reduced the velocity and volume of the debris flows, compared to the cases without a barrier. When the initial volume of source debris was small, or when the flow path was short, the barriers effectively blocked the debris flow regardless of the source-to-barrier distance. However, with a long flow path, installation of the barrier closer to the initiation location appeared more effective by preventing the debris volume from growing by entrainment. Our results contribute to a better understanding of how source-to-barrier distance influences debris-flow behavior, and show that the methodology presented herein can be further used to determine optimum and efficient designs for debris-flow barriers.

To further the understanding of how rivers process an instantaneous increase in sediment supply, a study was conducted on a gravel bed channel affected by post-fire debris flows. Sleeping Child Creek (SCC), a tributary to the Bitterroot River, is located in the Sapphire Mountains of west central Montana. In the summer of 2000, 1500 acres of Sleeping Child Creek's watershed were burned with 65% of the area classified as a high severity burn. The following summer, intense rainfall triggered debris flows originating from tributaries of SCC with headwaters toward the north. Along a 7-km study reach, cross sections, longitudinal profiles, and pebble counts were taken near six debris flow fans in the summer of 2012. In addition, a survey of large woody debris was conducted throughout the study reach. Comparing the cross sections with previous measurements taken in 2005 demonstrated that the river has aggraded throughout most of the site. Moreover, the bed material has become finer and the amount of large woody debris has doubled in seven years. The aggradation and the finer bed material can be attributed to the increase in large woody debris in the channel. In 2005, the large woody debris in the channel was from the debris flows. Since then, the large woody debris has come from the severely burned slopes adjacent to the channel. The response documented in Sleeping Child Creek demonstrates that, in burned mountainous landscapes, the large woody debris can be instrumental in modulating the storage and release of the post-fire sediment pulse.

Frequency–magnitude distribution of debris flows compiled from global data, and comparison with post-fire debris flows in the western U.S.

Due to their potentially long runout, debris flows are a major hazard and an important geomorphic process in mountainous environments. Understanding runout is therefore essential to minimize risk in the near‐term and interpret the pace and pattern of debris flow erosion and deposition over geomorphic timescales. Many debris flows occur in forested landscapes where they mobilize large volumes of large woody debris (LWD) in addition to sediment, but few studies have quantitatively documented the effects of LWD on runout. Here, we analyze recent and historic debris flows in southeast Alaska, a mountainous, forested system with minimal human alteration. Sixteen debris flows near Sitka triggered on August 18, 2015 or more recently had volumes of 80 to 25 000 m3 and limited mobility compared to a global compilation of similarly‐sized debris flows. Their deposits inundated 31% of the planimetric area, and their runout lengths were 48% of that predicted by the global dataset. Depositional slopes were 6°–26°, and mobility index, defined as the ratio of horizontal runout to vertical elevation change, ranged from 1.2 to 3, further indicating low mobility. In the broader southeast Alaskan region consisting of Chichagof and Baranof Islands, remote sensing‐based analysis of 1061 historic debris flows showed that mobility index decreased from 2.3–2.5 to 1.4–1.8 as average forest age increased from 0 to 416 years. We therefore interpret that the presence of LWD within a debris flow and standing trees, stumps, and logs in the deposition zone inhibit runout, primarily through granular phenomena such as jamming due to force chains. Calibration of debris flow runout models should therefore incorporate the ecologic as well as geologic setting, and feedbacks between debris flows and vegetation likely control the transport of sediment and organic material through steep, forested catchments over geomorphic time. © 2020 John Wiley & Sons, Ltd.

Volumes of erosional and depositional landforms were estimated by reconnaissance methods in a steep, forested drainage basin in the Coast Mountains, British Columbia, to examine what insight might be gained into the sediment budget in the absence of direct measurements of sediment transfer processes. Sediment transfers in postglacial and contemporary time were inferred on hillslopes, from hillslopes to stream channels, and within stream channels. Data were collected by aerial photographic analysis and field observation. Postglacial sediment sources are major gullies created by debris slides and flows, failures in glaciolacustrine terraces, and stream channel degradation. Depositional landforms consist of talus slopes and colluvial and alluvial fans. Contemporary sediment sources include debris slides and flows in established gullies and minor processes on hillslopes. Debris slide and debris flow volumes were calculated, and other processes were estimated from regional values. Erosion rate averaged over postglacial time is 276 t·km–2·a–1(0.15 mm·a–1surface lowering), with gullies and stream channel degradation contributing 170 and 82 t·km–2·a–1, respectively. A terminal alluvial fan provides an independent check of the results. In contemporary time, erosion rates are calculated to be 350 t·km–2·a–1, with debris flows and slides contributing nearly all of this sediment. The contemporary rate is probably perturbed by recent land use history. Mass-movement processes appear to be the dominant mechanism of sediment transfer and, contributions from Pleistocene valley deposits have declined significantly during Holocene time.

Predicting debris flow runout is of major importance for hazard mitigation. Apart from topography and volume, runout distance and area depends on debris flow composition and rheology, but how is poorly understood. We experimentally investigated effects of composition on debris flow runout, depositional mechanisms, and deposit geometry. The small‐scale experimental debris flows were largely similar to natural debris flows in terms of flow behavior, deposit morphology, grain size sorting, channel width‐depth ratio, and runout. Deposit geometry (lobe thickness and width) in our experimental debris flows is largely determined by composition, while the effects of initial conditions of topography (i.e., outflow plain slope and channel slope and width) and volume are negligible. We find a clear optimum in the relations of runout with coarse‐material fraction and clay fraction. Increasing coarse‐material concentration leads to larger runout. However, excess coarse material results in a large accumulation of coarse debris at the flow front and enhances diffusivity, increasing frontal friction and decreasing runout. Increasing clay content initially enhances runout, but too much clay leads to very viscous flows, reducing runout. Runout increases with channel slope and width, outflow plain slope, debris flow volume, and water fraction. These results imply that debris flow runout depends at least as much on composition as on topography. This study improves understanding of the effects of debris flow composition on runout and may aid future debris flow hazard assessments.

The authors predicted mudflow and debris flow velocities with 350 field and laboratory measurements and obtained a slight trend for V=u to increase with h=d50, and the ratio of V=u is approximately 10 and rarely exceeds 30. The logarithmic relationship of turbulent resistance agrees reasonably well with the measurements. Good results also are obtained with the Manning-Strickler approach. The dispersive stress equation also compares well with the measurements, but only when h=d50 < 50. For individual series data such as Hashimoto and Puerco, there is a trend for V=u to increase with h=d50. But for other individual series data, such as Davies, Rickenmann, Paris, St. Helens, Nevado del Ruiz, Wanglin, Wang, and Wenhai, there is a trend for V=u to decrease with h=d50 (Fig. 1 of the original paper). The discussers have some field data that also show the same phenomena, which are presented in Fig. 1. The database includes a total of 119 flow velocity measurements, in which each point includes flow depth, density (volume concentration), median grain diameter d50, grain size of less than 10% d10, and slope. The data of Jiangjia were obtained from the field measurements in Jiangjia Ravine, Yunnan, China, in 1999; the data of Hunshui were obtained from the field measurements in Hunshui Ravine, Yunnan, China, in 1976–1978; and the data of Liuwan were obtained from the field measurements in Liuwan Ravine, Gansu, China, in 1963–1964. Debris flow is a gravity flow, and the gravity force plays an important role in the movement of debris flow. So debris flow moves on a large slope and deposits on a small slope. Fei and Su (2004) showed that the main driving force of debris flow is provided by particles, not by water. When the volume concentration C 1⁄4 0:27 (or density ρ 1⁄4 1:46 g=cm3), the particle driving force is equal to the water driving force in debris flow. They defined that the minimum density of debris flow is 1:46 g=cm3, and 1:46–1:80 g=cm3 is the range of less viscous debris flow. Therefore, the minimum density of viscous debris flow is 1:80 g=cm3 (C 1⁄4 0:47) (Fei and Su 2004; Yu 2008a). Two types of shear stresses describe these two kinds of debris flows: (1) the viscous stress for viscous debris flow, and (2) the turbulent stress for less viscous debris flow (Fei and Su 2004). As the driving force of debris flow is provided primarily by particles, the coarse particle plays an important role in the movement of debris flow, especially for viscous debris flow with a large volume concentration. So the larger particle diameter, the larger is the velocity. There is a trend for V to increase with d50, and there is a trend for V=u to increase with d50=h. Two empirical equations of mean velocity of viscous debris flow described this relationship (Wu et al. 1993; Yu 2001): V u 1⁄4 27:57 d50 h 0:245 ð1Þ

Experience shows that debris flows containing large woody debris (driftwood) can be more damaging than debris flows without driftwood. In this study, the deposition process of debris flows carrying driftwood was investigated using numerical simulations and flume experiments. Debris-flow trapping due to driftwood jamming in a slit-check dam was also investigated. A numerical model was developed with an interacting combination of Eulerian expression of the debris flow and Lagrangian expression of the driftwood, in which the fluctuating coordinates and rotation of the driftwood were treated stochastically. The calculated shapes and thicknesses of a debris-flow fan and the positions and orientations of the deposited driftwood on a debris-flow fan were consistent with experimental flume results. The jamming of driftwood in a slit-check dam was evaluated based on geometry and probability. The simulated results of outflow discharge and the proportion of driftwood passed through the slit-check dam also agreed with the experimental results.

Gravitational mass movements, such as debris and mud flows, are among the most destructive natural hazards, leading to substantial fatalities and extensive economic damage worldwide. Improving our understanding and modeling of these processes is essential for developing effective risk management and early warning strategies. When debris and mud flows pass through a curved channel, centrifugal forces may cause a difference in flow height between the inner and outer banks of the channel. This height difference, known as superelevation, can be described using analytical models that establish a relationship between the superelevation height and the flow velocity.Analytical models often employ a forced vortex approach, incorporating parameters such as the cross-sectional slope of the flow surface, flow width, and bend radius. These models, however, rely on assumptions such as a linear flow surface between mud deposits on the banks, a rectangular cross-section, and neglect both complex rheological behaviors and solid-fluid interactions. As a result, an empirically determined correction factor is required within the formula. The absence of a clear mechanical rationale for this correction factor presents challenges, as it is currently derived only through field investigations and laboratory experiments.This study presents an enhancement to the existing forced vortex approach by incorporating insights from numerical modeling. A coupled SPH-DEM numerical model is employed, where DEM particles represent coarse solid particles, and SPH accounts for the fluid phase, comprising fines and water. The SPH-DEM coupling is based on the no-slip interaction model, with simulations performed using a GPU-based solver to ensure enhanced computational efficiency. To validate the approach, a parametric test is conducted, initially back-calculating laboratory-scale experiments. The study further involves varying the water content in debris and mud flows to examine its impact on flow behavior and superelevation. Larger water contents lead to an increased superelevation angle. Results from the parametric test reveal a clear correlation between water content and the flow surface shape in curved channels. Specifically, mud flows are characterized by convex upward surface shapes, whereas more granular debris flows typically exhibit concave downward shapes.The distribution of material within the cross-section of the flow is governed by the equilibrium between boundary forces and centrifugal forces acting on the flow, which directly influences the superelevation. Numerical investigations are conducted to determine a correction factor and assess the extent to which inertial effects contribute to this correction factor for different material mixtures. Furthermore, we demonstrate that the effect of the flow surface shape is significant and is currently only accounted for by the empirical correction factor. This study offers new physical insights for the back-calculation of debris flow velocities in curved sections marked by mud deposits. Large-scale SPH-DEM simulations of a real debris flow event at Illgraben (Switzerland) are performed, showing good agreement with field data and its potential for further real-scale modeling.

The role of large woody debris in modulating the dispersal of a post-fire sediment pulse

Our aim is to determine the run-out distance of the debris flow that is crucial in the assessment, prevention and control of the debris flow hazard. Based on the variation characteristic of debris flow velocity in the alluvial fan, this paper proposes the calculation method of the velocity attenuation coefficient of the debris flow. By defining the velocity attenuation coefficient and deducing its calculating formula, this paper puts forward a new method to determine the run-out distance of the debris flow based on the velocity attenuation coefficient, and Gangou debris flow in Luding County, Sichuan Province is selected as a case for calculation and verification. Having 10 m as its measuring spacing, this paper measured 19 sections at the alluvial fan of the Gangou debris flow (among them, 11 sets of data are valid). And based on the measurement, this paper analyzes the characteristic of the velocity attenuation and calculates its velocity attenuation coefficient after the 2005 debris flow. The study indicates that when the velocity of Gangou debris flow at the alluvial fan is greater than 12 % of the initial velocity (at the mouth of gully), the attenuation is quite remarkable. But when the velocity at the alluvial fan is less than 12 % of the initial velocity, the attenuation is quite slow. Besides, when Gangou debris flow rushes out of the gully mouth (the initial velocity is 10 m/s) and when it attenuates to the 32 time, its velocity is less than 0.1 m/s, the debris flow is considered to stop flowing, and the run-out distance of Gangou debris flow is calculated to be 320 m. But the present alluvial fan of Gangou debris flow is measured to be 285 m in length, and the calculated run-out distance is 320 m, which is 35 m longer than its present length. This means when the debris flow runs out in 2005, it blocked up the main river (Dadu River) in some extent. And this finding is generally in accordance with that from the field survey. The findings can be of theoretical and practical significance in the debris flow hazard assessment, as well as its prevention and mitigation.

In order to effectively prevent and control debris flow disasters, it is one of the most effective engineering measures to construct a blocking project in the valley section of the circulation section to actively reduce the flow and velocity of debris flow. The check dams is taken as a control structure to intercept the solid materials of flush flood and debris flow in the channel, which play an important role in debris flow prevention and control project. However, the huge impact force of debris flow often causes the check dam to be destroyed. Based on the theory of fluid–solid coupling, by means of structural dynamics simulation software, fluid dynamics simulation software and coupling plane data transfer software, numerical simulation of the interaction between debris flow and the check dam has been carried out. Research results indicate that (1) the overflow section is subjected to the impact of the debris flow greater than the abutment; (2) debris flow impact only has a greater impact on a certain depth of the check dam, and the impact beyond this range is less. The numerical simulation results can provide data support for the design and optimization of the check dam.

Debris flows frequently occurred in Wenchuan earthquake region from 2008 to 2010, resulting in great damage to localities and being a prolonged threat to reconstruction. Forty three events' data including debris-flow volume, sediment volume and watershed area are analyzed and compared with other debris-flow events in Eastern Italian Alps, burned areas in USA and in Taiwan. The analysis reveals that there is a strong empirical relationship between debris-flow volume and loose materials volume in the earthquake region. In addition, the relationship between debris-flow volume and watershed area in the earthquake region has a wider variation range than that in other three regions while the debris volume also appears to be larger than that in the other three regions, which implies the volume of debris flows with strong influence of earthquakes is larger than that with no such influence and it is hard to predict the post-quake volume only by the watershed area. The comparison of the maximal debris-flow erosion modulus in the Wenchuan region and in Taiwan indicates that debris flows will be very active in a short time after strong earthquake.

Experimental investigation of the effects of shrub filter strips on debris flow trapping and interception

Abstract. This article presents MassMov2D, a two-dimensional model of mud and debris flow dynamics over complex topography, based on a numerical integration of the depth-averaged motion equations using a shallow water approximation. The core part of the model was implemented using the GIS scripting language PCRaster. This environment provides visualization of the results through map animations and time series, and a user-friendly interface. The constitutive equations and the numerical solution adopted in MassMov2D are presented in this article. The model was applied to two field case studies of mud flows on torrential alluvial fans, one in the Austrian Tyrol (Wartschenbach torrent) and the other in the French Alps (Faucon torrent). Existing data on the debris flow volume, input discharge and deposits were used to back-analyze those events and estimate the values of the leading parameters. The results were compared with modeling codes used by other authors for the same case studies. The results obtained with MassMov2D matched well with the observed debris flow deposits, and are in agreement with those obtained using alternative codes.