## 1 Introduction
They employed the model proposed by [PERSON] [19] and [PERSON] and [PERSON] [20] to compare the results with those found by other authors for the same location [19; 21; 22; 23; 24; 25]. Two types of slope stability analyses were performed, one considering the existence of a tension crack, and the other not considering it. In both cases, cohesion, friction angle, tension crack depth, and water column height were modified in order to determine the influence of height on the design of rock slopes. The authors concluded that it is effective to include the variability and uncertainty of rock properties and slope geometry in the design. [PERSON] et al. [26] and [PERSON] and [PERSON] [27] used the control variable method (CVM) and the orthogonal design method (ODM) to study the influence of freeze-thaw cycles on slope stability in sandstones and dry cohesionless slopes in the Chinese provinces of Xinjiang and Yingpan, respectively. [PERSON] et al. [26] found that safety factor decreased as density, cohesion, and friction angle decreased. The authors concluded that the friction angle presented the highest sensitivity index in their case study. Meanwhile, [PERSON] and [PERSON] [27] used a simplified analytical model with a plane failure surface and concluded that the parameters with the greatest sensitivity index were both friction angle and cohesion. [PERSON] et al. [28] and [PERSON] et al. [29] implemented the Grey relational method (GRM) [30] combined with the Universal Distinct Element method (UDEM) to determine the factors influencing toppling and bending of layers of rocks inclined perpendicular to the angle of the slope. [PERSON] et al. [28] carried out a theoretical study based on normal parameters, while [PERSON] et al. [29] applied a similar methodology to the study of a slope located on the Yalong River in southeast China. This method requires fewer tests than the orthogonal test method (OTM) and allows the influence of sub-factors on the main factor to be determined. These sub-factors are the geometry, the physical and mechanical parameters of the materials, and the mechanical parameters of the discontinuities or joints. [PERSON] et al. [12] applied the combined slope hydrology and stability model (CHASM) to analyze the slope stability of sandy and clayery materials under the influence of the water table position, using cases found in the literature [31; 32]. They performed modeling using the limit equilibrium method, specifically [PERSON]'s method [33]. The soil strength parameters were analyzed using the variance-based sensitivity analysis (VBSA) method of [PERSON]'s [34] and the replicated Latin hypercube sampling (rLHB) [35]. In the analyzed cases, the authors reported that the parameter with the greatest influence in sandy clay is cohesion, while in sand, it is the friction angle. [PERSON] et al. [36] used the least angle sensibility regression (LASR) [37] and [PERSON]'s method combined with Monte Carlo simulation to determine the range of input factor values. The authors used the simplified analytical expression for infinite slopes, and stated that LASR presents better performance and that it is not possible to distinguish between sensitivity to slope geometry and soil parameters, with slope inclination being the factor with the highest sensitivity index. [PERSON] et al. [15] implemented the sparse polynomial chaos expansion (SPCE) model to quantify uncertainty, combined with the upper bound limit analysis theory and a pseudo-static analysis using a horn-type failure surface. The probabilistic analysis was performed using Monte Carlo simulation, defining the performance function and failure probability. Latin hypercube sampling was used to obtain the cohesion, friction angle, horizontal seismic coefficient, specific weight, slope angle, and width-height ratio parameters required for the 3D model.