Discrete element modeling of the mass movement and loose material supplying the gully process of a debris avalanche in the Bayi Gully, Southwest China

• Three-dimensional mass movement of debris avalanche is simulated by DEM.
• FORTRAN and AutoCAD are used for the secondary development of PFC3D.
• Evolution of landslide travel distance, velocity and volume supply to gully is presented.
• Concept of loose materials supply gully process is presented and discussed.
• Mass movement influenced by landslide volume, topography and parameters are analyzed.

The dynamic process of a debris avalanche in mountainous areas is influenced by the landslide volume, topographical conditions, mass-material composition, mechanical properties and other factors. A good understanding of the mass movement and loose material supplying the gully process is very important for understanding the dynamic properties of debris avalanches. Three-dimensional particle flow code (PFC3D) was used to simulate a debris avalanche in Quaternary deposits at the Bayi Gully, Southwest China. FORTRAN and AutoCAD were used for the secondary development to display the mass movement process and to quantitatively describe the mass movement and loose material supplying the gully process. The simulated results show that after the landslide is initiated, the gravitational potential energy is converted into kinetic energy with a variation velocity for the sliding masses. Two stages exist for the average-movement velocity: the acceleration stage and the slowdown stage, which are influenced by the topographical conditions. For the loose materials supplying the gully process, the cumulative volume of the sliding masses into the gully gradually increases over the time. When the landslide volume is not large enough, the increasing landslide volume does not obviously influence the movement process of the sliding masses. The travel distance and movement velocity increase with the decreasing numerical parameters, and the mass-movement process is finished more quickly using low-value parameters. The deposition area of the sliding masses decreases with the increasing numerical parameters and the corresponding deposition thickness increases. The mass movement of the debris avalanche is not only influenced by the mechanical parameters but is also controlled by the topographical conditions.