<p>Comprehensive landslide risk zoning in open-pit mines is fundamental to precise safety monitoring, early warning, and disaster prevention and control. Existing approaches often rely on single static indicators and exhibit limited capacity for dynamic risk regulation. To address these limitations, this study proposes a scientifically grounded slope zoning method that couples multiple factors across "engineering-geological-environmental (seasonal variations, distribution of infrastructure)". First, a hierarchical structure model comprising four key influencing degree of crack development, slope angle, slope height, and rock and soil properties—is constructed using the Analytic Hierarchy Process (AHP). By computing the weight coefficients of these evaluation indices and deriving a slope hazard index, we classify slope hazard zones. Additionally, the two-dimensional rigid-body Limit Equilibrium Method (LEM) is applied to calculate slope factor of safety for multiple cross-sections, facilitating the delineation of slope stability zones. To integrate these results, we employ a combined cross-matrix analysis and simple weighted averaging method, further refining the comprehensive risk zoning by dynamically incorporating factors such as slope along-strike surface shape, seasonal (climatic) variations, fault characteristics, and distribution of infrastructure. This approach is validated through application to a case study in an open-pit mine, where one higher-risk, five medium-risk, and six low-risk zones are identified. Notably, the higher-risk zone exhibits strong spatial agreement with actual landslide occurrences. The results demonstrate that the proposed method significantly enhances zoning accuracy, providing a theoretical basis for the optimized allocation of monitoring resources, the development of differentiated early-warning models, and the full-lifecycle management of slopes in open-pit mining operations.</p>