Toward Sustainable Geohazard Assessment: Dynamic Response and Failure Characteristics of Layered Rock Slopes Under Earthquakes via DEM Simulations

Understanding the dynamic response and failure mechanisms of rock slopes during earthquakes is crucial in sustainable geohazard prevention and mitigation engineering. The initiation of landslides involves complex interactions between seismic wave propagation, dynamic rock mass behavior, and crack network evolution, and these interactions are heavily influenced by the slope geometry, lithology, and structural parameters of the slope. However, systematic studies remain limited due to experimental challenges and the inherent variability of landslide scenarios. This study employs Discrete Element Method (DEM) modeling to comprehensively investigate how geological structure parameters control the dynamic amplification and deformation characteristic of typical bedding/anti-dip layered slopes consist of parallel distributed rock masses and joint faces, with calibrated mechanical properties. A soft-bond model (SBM) is utilized to accurately simulate the quasi-brittle rock behavior. Numerical results reveal distinct dynamic responses between bedding and anti-dip slopes, where local amplification zones (LAZs) act as seismic energy concentrators, while potential sliding zones (PSZs) exhibit hindering effects. Parametric analyses of strata dip angles and thicknesses identify a critical dip range where slope stability drastically decreases, highlighting high-risk configurations for earthquake-induced landslides. By linking the slope failure mechanism to seismic risk reduction strategies, this work provides practical guidelines for sustainable slope design and landslide mitigation in tectonically active regions.