Study of critical failure surface influencing factors for loose rock slope

In current research, a novel geometrical technique is introduced to estimate the length (L) of sliding arc by considering slip surface entry point distance (de). Morgenstern-Price limiting equilibrium and strength reduction methods are selected within the framework of SLOPE/W and FLAC3D softwares respectively, to study the sliding mechanism and stability-effecting factor for a loose rock slope comprise of two different material. Also, the variability influence of strength parameters (cohesion force and friction angle) and the depth of pore water pressure on safety factor, failure depth (D) and distribution range of critical failure surface (CFS) are investigated through software. The results show that the shear strength parameters and the depth of pore water pressure plays a significant role in altering overall stability and the distribution range of CFS for a loose rock slope. Shallow failure with composite failure surface is much easier to occur for loose rock slope subjected to pore water pressure and material heterogeneity. These findings can be effectively used for loose rock slope reinforcement design.


Introduction
Engineers mostly investigate rock slope stability to precisely classify the most distinct 2D/3D possible critical slip surface [1]. Engineers also use several analytical methods to evaluate the stability factor of natural or artificial slopes [2]. Limit equilibrium techniques are leading among these methods, in which the empirically based established sitestability charts may be applicable for slope stability calculation in certain conditions [3]. In general, nevertheless, laboratory tests, field assessment and more precise numerical solutions are normally suggested [4]. The outcomes of above studies are usually reported as the probability of failure and FS for the geotechnical parameters [4].
Although the limiting equilibrium (LE) method extensively used for stability analysis, is lacking in definite respects: predominantly, the incapability of the technique to consider the multifaceted constitutive relations of the slope materials [5,6]. Furthermore, during LE calculations assumptions are made by which water influences are isolated, because the influences of preliminary stresses are not straightforwardly incorporated, and this is impossible to compute local safety factors for slope [7,8]. So, the position of the CFS is found capably and comparatively easily if the resultant explanation of the slope is adequate with the range of the available geotechnical data.
Limit-equilibrium approaches are also widely selected for slope stability investigation, mainly due to their ease of use and simplicity. Among these techniques, the technique of slices selected with a circular sliding surface (CSS) and frequently employed when the circles are often approximate and convenient to analyse the studied failure surface [2]. Inappropriately, there are several cases, in which simple slip-circle technique cannot be applied for stability analysis and any choice of the nature of the failure surface cannot be commonly useable for further study. This is particularly due to the heterogeneous profiles of slope wherein a critical sliding surface (CSS) is anticipated to make parallel with the line of minimal resistance for example anisotropy, bedding, etc. In the site, examinations of actual sliding surfaces show that circular failure surfaces can indicate the variety of characters with severe changes at the boundary between dissimilar geological materials [9,10]. Hence, circular failure surface can referred to the distinct case of noncircular sliding exploration where the previous is considered a systematic method. Therefore, the conventional limit equilibrium study for the slope stability factor assessment includes two phases: (a) development of a law of resemblance in between the safety factor and a probable sliding surface. (b) Searching for the minimum value of safety factor (FS) and over all realistic sliding surface, which shows minimum value of FS, is refer to critical sliding surface (CSS).
Significant effort already has been made to the first aspect; the techniques of Bishop [11], Baker [12], extended Spencer [13], and a global procedure of Morgenstern and Price [14] are well known today. In the second aspect, various techniques have been used to analyze the location of the critical sliding surfaces [15,16]. These methods are not particularly laborious [17]. But, numerical methods are useful for slope stability assessment [18], and key block analytical method can be used for rock slope failure analysis [19]. Blocky structural sliding of rock slope is analyzed by using fuzzy key-block method [20]. Simple genetic algorithm [21] and ant colony optimization [22] are applicable for critical non-circular failure surface of rock slope stability analysis. A suitable method of stability calculation can be determined by a combination of the above logical methods. The position and size of the critical failure surface cannot be determined well with limit-equilibrium techniques. But, nural network and optimization algorithm could be easily selected to define its position [23,24].
In current research, a simple geometrical scheme is introduced first time to calculate the failure arc (L) of a loose rock slope. SLOPE/W computer code is used to analyze the failure mechanism of CFS of a real case of instability. In addition, the effect of cohesion, internal friction and pore water pressure on the distribution of CFS is studied by software. At last, FLAC 3D software is also selected to check the correctness of results evaluated from SLOPE/W. In these findings, the CFS allocation will provide concrete guidelines for choosing the more typical sliding surface in risk assessments and system reliability analyses of heterogeneous loose rock slope beside road.

Methodology
There are numerous approaches selected for slope stability investigation. Among these approaches, the limiting equilibrium (LE) technique of slices [25][26][27] have attracted significant attention, due to its accuracy and simplicity. In this process, the ratio between the resisting forces to driving forces on a possible slip surface is welldefined as the probability of slid. Only if the intended value of FS more than one then slope is thought-out safe. Limiting equilibrium calculations consider the slope failure surface as a rigid body [28]. Sliding body or surface delivers necessary information and has been extensively adopted previously to evaluate the stability of all types of slopes. For example the stability of heterogeneous slopes was calculated with the help CFS to take safety measures [27], for dam slopes [29], rock slopes [30], highway slopes [31] and rock cliff face [32].
The limiting equilibrium procedures are frequently used as the systematic approaches to determine the stability of weathered rockslides. In geotechnical engineering field, investigating the stability of slopes are also one of the oldest kinds of study. In the 20th century many researcher presented the knowledge of dividing the probable sliding mass into several number of slices for stability determination purpose [27,33]. In the 1960's due to the dawn of computers technology, which made it conceivable to extra voluntarily switch for iterative measures which are essential scheme to control the statistically more laborious equations [33]. On the other hand, the idea of numerically discretizing a large slippery mass into small parts for investigation objectives is relatively innovative at the time. Moreover, in LE calculations assumptions are made first to generalize the effects of water pressure and comprehensively incorporate the influence of primary stresses. In addition, this is impossible to evaluate potential FS. However, if the resultant generalization of the problem is appropriate or acceptable with the scope of the existing geotechnical information than the position of the CSS could be found quite efficiently easily. SLOPE/W computer program is a product of Geo-Studio software system. One of the more commanding topographies of this joined tool is that this provides access to sorts of investigates of a complex and broader field of difficulties, together with the selection of stresses and calculated water pressures in the stability study. It is possible in limit equilibrium (LE) computer programs such as SLOPE/W to deal with extremely irregular water pressure, complex stratigraphy, variety of nonlinear and linear strength models, material heterogeneity, concentrated loads, structural reinforcement and practically any type of sliding surface shape [34]. Then search method is prolonged to compare the CSS for a loose slope.

Determination of shear strength parameters
To estimate the unknown material constants, such as cohesion (c) and internal friction angle (ϕ), a new C-tanφ relationship chart is first created from multiple sliding surfaces within the framework of SLOPE/W tool as shown in Fig. 1. Two curves showing the C/ϕ relationships can be plotted by using the method of slice [35]. Set the hypothetical safety factor of the said slope intended for repeated calculation to be one. The following lists of strata whose soil constants are; sandy soil: unit weight γ t = 18.0 (kN/m 3 ), friction angle ϕ (°) and cohesion force c (kN/m 2 ) were determined from reverse calculation by using density of water γ w is 9.80 (kN/m 3 ). Weathered rock: unit weight γ t = 20.0 (kN/m 3 ), friction angle ϕ = 34° and cohesion c = 15 (kN/m 2 ). The unknown constants (c and ϕ) for sandy soil were then determined from C-tanϕ relational expression diagram ( Fig. 1), by assuming that 120 m 2 area of slope have c = 1 KN/m 2 . With reference to the Fig. 1, the value of c and tanϕ is noted c = 15.3 KN/m 2 and tanϕ = 0.41 (ϕ = 25°). Figure 2 shows a representative failure surface (AB), failure arc (L), slices discretization scheme and the forces acting on a slice of the slope. Normal forces (E) and shear forces (X) are acting on the base and one sides of the slice, respectively (Fig. 2). Limit equilibrium (LE) calculation is consist on two safety factors calculations and permits for a boundary of inter-slice normal-shear force circumstances. One calculation gives the FS regarding horizontal force (F f ) while the other calculation provides the FS with regard to moment (F m ). An equation to handle inter-slice shear forces can be written as:

General limit equilibrium formulation
where X =shear force, E =normal force, λ =dimensionless function and f(x) =a function. In current analysis, a distinctive half-sine function is chosen in SLOPE/W software as shown in Fig. 3. In this figure, the upper and lower curves are actual specified and used functions. Lambda (λ) is the ratio between the two curves. One of the basic problems in LE calculation is knowing how to describe interslice and FS equation relating to moment equilibrium is: The safety factor (FS) equation regarding horizontal force equilibrium is: Base normal is defined as: In above equations, ϕ´ = effective internal angle of friction, c´ = effective cohesion, N = base normal force, u = water pressure, W = slice weight, α = slice base inclination, D = point load and β, x, R, f, ω, d are geometric parameters. Equation 4 is acquired by the summation of all vertical forces. F = Fm, when base normal force N is substituted in the Eq. 2 and F = F f when N is substituted in Eq. 3. The normal force (N) is dependent on the inter-slice shear forces (X R and X L ) on both sides of the slice.
Limit equalibrium formulation calculates F f and F m a choice of dimensionless function (λ) values. A plot is shown in Fig. 4, which was computed to illustrate how F f and F m vary with lambda (λ). At the end, for the whole slope there is only one FS. F f and F m are similar when both force and moment equilibrium were satisfied and FS was same for each slice.

Determination of failure arc
A new geometrical procedure is developed to calculate L. Figure 5 shows a possible geometry of a failure surface with entry point E and exit point M. sliding surface entry and exit points are assumed in this particular research. From Fig. 2, x n+1 and x 1 are the x-coordinates of entrance and exit points on the slope face. The y-coordinates of a point on the sliding surface have entrance and exit points, which can be used to find out the length of failure arc. To solve this problem, it is assumed that the following Eq. 5 satisfied the failure circle formula. In Eq. 5, x 0 , y 0 and r are unknown variables. To find these variables three equations are needed. Since, entry and exit points can satisfy the Eq. 5. Consequently, two of the equations created by inserting their coordinates in the Eq. 5 as: and the total length of sliding cure is given as: From Eq. 6 L is length of failure arc which satisfy the Eqs. 8 and 9 as: The three equations, which are needed, can be evaluated with the help of Eq. 9 as follow: The Eq. 10 can be solved by putting the value of geometric parameters (a, de, α, β, and L) and a = h/cosβ, and h is the height of slope and de is entry point distance. In Eq. 10, r is radius and 0 , and 0 are the coordinates of slip circle. However, this geometric methodology is suitable for circular slip surface, where there are only three different control variables, so this method cannot be used to composite and non-circular slip surface. In addition, this method is not used material shear strength properties, thus there is no need to conduct laboratory test. Therefore, proposed geometrical method is simple and time saving as compare to other methods of slices.

Illustrative example
The promotion of China Pakistan Economic Corridor (CPEC) strategy has led to a large number of construction and exploitation in Pakistan. The stability of such largescale engineering activities on hills or mountains along road side directly related to the progress of the project, economy and the safety of the workers [36]. Hence, the slope stability analyses are conducted by engineers to calculate the stability of road slopes, dam slopes, retaining wall and excavated slopes with critical sliding surface [35]. A loose rock cut slope is located along National Highway (N70) in Pakistan considered as a case study. The slope material is composed of boulder. Terrace is distributed in vast area of the project. It consists of 20-1000 mm diameter of rounded rock pieces (Fig. 6). This formation is valley floor sediment and supposed to consist of two sandy layers. Water pressure measured with piezometer. The coordinates of water line are presented in Table 1.

Model parameters
The width and height of 70 m 35 m respectively, measured with total station (Fig. 6). Slope angle of 45° was recorded.
The shear strength parameters of slope material assigned to corresponding region (Fig. 7a). Half-sin function and Mohr-Coulomb failure criterion is selected in current study. Trial slip surface entry and exit point are located as shown in Fig. 7a. Stability factor and the location of the CFS are then simulated by SLOPE/W. An appropriate manner of realistically showing the result of the CFS is a safety map (Fig. 7b). All trial failure surfaces are shown with corresponding factor of safety. The location and distribution of all possible CFS   Fig. 7b. Figure 7b shows that critical failure surface initiate from crest and cut the slope near the toe.

Sliding mechanism for typical failure surfaces
In this sub-section the failure mechanism and the location of failure surfaces is analyzed at the various value of cohesion, friction angle and unit weight. . Sliding surface location and its sliding style for three sets of input parameters (c, ϕ, γ) are plotted in Fig. 8. The representative identifications according to Fig. 8 have many characteristics of CFSs location by varying the strength parameters. In Fig. 8, CFSs are located between the crest and the toe of loose rock slope and the sliding of the slope is a local failure at average value of strength parameters. The representative insights in Fig. 8 have different appearances of CFS failure mechanism and its location, which is categorized into three parts as follows: (1) In first case, the CFS was shallow and cut the toe of slope as shown in Fig. 8a. The weak zone is comparatively small and is situated around the surface of the slope. The location of the CFS, in this case, is well dependable with the position of the weak zone and overall failure style is a local failure. Because the factor of safety (FS) for the CFS passing through the weak zone is smaller than that for other big CFSs representing an overall slope failure. When the failure of a slope surface happens in such a manner that the sliding surface cuts the slope surface or passes from slope toe, it is named slope failure as shown in Fig. 8a, b. In this case, the slip surface exit and entry points are respectively situated near the top and on the toe of slope surface. (2) When the strength parameters were kept relatively large, the CFS was deep (Fig. 8b). In second case, the weak zones are relatively large and situated at the upper most portion of slope. The associated failure mechanism for a typical CFS is a narrow failure as shown in Fig. 8b. The observation for this particular case is that, the value of safety factor was lesser for big sliding surfaces and passing through the weak area of the slope that showing an overall circular failure. (3) When the large value of strength parameters selected the CFS is located at the little distance from the toe and the slip mode of a loose rock slope is global slip (Fig. 8c). Because safety factor (FS) for a CFS is bigger as compared to other all-valid slip surfaces in the weak areas of slope. In this particular example, the loose rock slope sliding is a deep-seated failure and the position of CFS not dependable with the actual situation of the weaker areas. Because the weight of slope material is the part of the gravitational force and this is related to the size of the sliding body. The resistance shear force of slope material contains on two components, cohesion part and internal friction part. Apparently, as the volume of the sliding body decreased, the diminishing amount of related CFS was minor than that of sliding body. Especially, if the input strength properties of slope kept constant and the volume of slippery body reduced, the weight of slope reduced more rapidly than the slope resistance. For a comparatively large CFS as displayed in Fig. 8c, the cohesion force is small enough to produce a safety factor for the weak area lesser than the safety factor for a global CFS.
It can be evidently seen from the above analysis, that the sliding characteristics for a heterogeneous loose rock slope (formed by colluvium material) with a moderately small value of strength parameters is a global slid in most cases. In certain cases, although potential slid may occur, but the related CFS have entry and exit points situated at the crest and near the toe of slope, respectively.

Discussion
In this section, the influence of cohesion force (c), internal friction (ϕ) and the depth of pore water pressure (PWP) on safety factor (FS), depth of failure surface (D), length of failure arc (L) and its distribution is discussed by using SLOPE/W software. The general failure surface with an exit or an entry points on the slope face are assigned to horizontal coordinates of lowermost and uppermost points of CFS. The distribution range of all valid slip surfaces show similar phenomena for this particular case.

Influence of cohesion
In this sub-section, the influence of cohesion (c) is analyzed for same slope model configuration. For this purpose, unit weight and angle of internal friction (ϕ) are chosen as 20 KN/m 2 and 30°, respectively. The value of cohesion (c) ranging from 10 KN/m 2 to 20 KN/m 2 is assigned in SLOPE/W software for analysis. Three different slope model are tested by varying the values of c as 10, 15 and 20 (KN/ m 2 ), respectively. The summary of computed safety factors are graphically portrayed and all-valid failure surfaces fall in the limit of all possible trial slips. The distribution range of failure surfaces are consistent with safety map (Fig. 9). Cohesion force has a substantial effect on the safety factor of loose rock cut slope principally for larger values.
As clearly displayed in Fig. 9a-c, when the cohesion force is increased, than there are enormous gapes among the allocations of all valid failure surfaces. In addition, the FS and maximum D increased linearly (Fig. 9). The reason is that the cohesion is a strength parameter increasing this value increase in FS. In other words, very large or local slip surface are more likely to occur, as the cohesion is kept enormous.
As shown in Fig. 9a, b, there are no significant differences in-between FS and the allocations of failure surfaces when c = 10 and c = 20 (KN/m 2 ) for upper and lower layer, respectively. Hence, the cohesion is set to be 20 and 25 (KN/m 2 ) for upper and lower layer, respectively. In addition, all valid sliding surfaces are a band of curves with different sizes, shape and locations. Safety factor and the distribution range of all valid failure surfaces suddenly increased with an increasing cohesion parameter as shown in Fig. 9c. The normal value of c with any probable failure surface is near to the average value of cohesion and a possible CFS is likely to run across both layers of slope material as shown in Fig. 9b, c. Hence, when the safety factors are comparatively small than the affected area and the distribution range of failure surfaces were larger for large safety factors as compare to small safety factors. This mechanism could be also examine by comparing the all-valid slip surfaces in Figs. 9 and 10.
The spatial changeability effect of cohesion on maximum failure length (L) is also examine in this study. The failure arc was supposed to be initiated from the crest (entry point) and cut the slope from or near the toe (Fig. 9). Figure 9 also shows the effect of c on L. Changing the value of c also change in failure arc length and position (Fig. 9). When cohesion was 15, 20 and 25, (KN/m 2 ) to be set the maximum length of failure arc is 48.8 m, 50.2 m and 50.4 m recorded, respectively. Failure arc length evaluated by the equation system presented in Sect. 2.3 was 50.0 m, which is consistent the result determined from SLOPE/W software. It is clearly seen in Fig. 9 that by increasing the value of c L increased significantly.

Influence of friction angle
This sub-section mainly discussed the changeability influence of angle of internal friction (ϕ) on the FS, maximum D, L and distribution range of all valid failure surfaces. As compare to cohesion parameters the friction parameter usually has a smaller effect on FS, depth and the slip surface distribution (e.g., [37]). Same procedure as for cohesion is adopted to investigate the effect of internal friction. By varying the value of friction angle, for both layers, three cases were investigated through SLOPE/W. The values of friction angle are set to be 20°, 25°, 30° and 25°, 30° and 35° for upper layer and lower layer, respectively. The reason is that the value of ϕ is varying from 20° to 35° for studied case. Unit weight and cohesion are set to be 20 KN/m 3 and 15 KN/m 2 , respectively. Figure 10a-c plot the value of D and the distributions range of all valid trial failure surfaces for slope considering different values of friction angle. As the value of ϕ increased, the FS and distribution range of failure surfaces also increased significantly, but decrease in the depth of slip surface (Fig. 10). Local failures happened in first two cases as ϕ = 20° and 25°, and 25° and 30° selected for both layers, respectively (Fig. 10a, b). When the value of ϕ is 30° and 35° chosen, no local failure occurs as shown in Fig. 10c. The reason is that ϕ is a strength parameter increasing this increase in material resistance force, which produced greater FS (Fig. 10c). In other words, when the value of ϕ is large, all trial failure surfaces are extremely likely to enter the slope from the crest (entry point) and cut the slope near the toe. By contrast, when the value of ϕ is comparatively small, local failures have an exit and entry points situated on the crest and around the toe of the loose rock slope.
There were no large variations among the distributions of sliding arcs as c and ϕ are modeled in software. The reason is that there are no clear weak zones were founded near the slope surface. In addition, the more representative CFS is deep in both cases (Fig. 10a, b). While, as friction increased the distribution range of possible failure surfaces also increased (Fig. 10c), and these insights agree well with the results published by Ahmed et al. [37]. Nevertheless, in case of friction angle CFS is not cross the lower layer of slope material as shown in Fig. 10. The reason is that, for angle of internal friction, the whole slope area is almost homogeneous. Figure 10 also shows the variability effect of internal friction (ϕ) on the L. The length of failure

Influence of pore water pressure
The pore water pressure is one of the key factor affecting slope stability. When the whole slip surface cuts the water zone, and or some slices of the sliding mass lays within the water zone then this leads to sudden rise in the water pressure on the failure surface, which decreases the value of FS and also change the depth (D) of the possible sliding surface of the slope. The shape of CFS of the slope might be circular or spherical and might be shallow or deep-seated depends on the pore water distribution and material profile in the slope [38].
In this sub-section, the effect of the depth of the water zone on the failure mode is also investigated choosing the minimum material strength properties for same slope model configuration. Figure 11 illustrates the locations of the potential CFSs when the depth of the water zone was 15 m, 20 m, and 25 m, respectively. The overall trend is that the CFS for composite type is deep than that of circular type (Fig. 11).
Two cases, based on material maximum and minimum shear strength properties (SSP), are considered to evaluate the different effects of pore water pressure on the  Fig. 12. FS of the loose rock slope, apart of the simulation technique, decreased by increasing the height of the water zone considering maximum and minimum material shear strength properties (Fig. 12). At both material strength properties, the safety factor of the cut slope showed a decreasing trend (linear) by increasing the height of the water zone. Which means that the increase in water zone depth plays a big role in altering the slope stability. The reason for this is that the slip surface is situated within the water zone and thus the variability effect of unit weight and other strength parameters are not much significant for the stability of the slope. At the material minimum strength parameters (SSP), the safety factor of the slope drops slowly (Fig. 12). As the depth of the water zone reached at 30 m, the FS maintain the value of 2.42 (Fig. 12). When the depth of the water zone was 1 m, the safety factor was almost consistent with first case (Fig. 12). The reason is that increasing the height of pore water increase in the weight of material, result in decrease in FS. In addition, c and ϕ are resisting forces by reducing these decreases in FS. It can be concluded that by ignoring the pore water leads to the overestimate of the FS of the slope.
The effect of water zone height on the depth of the slope failure modes (critical or composite) investigated as shown in Fig. 13. The depth of both CFS and composite failure surface suddenly decreased first and then increased by increasing the depth of the water zone (Fig. 13). At maximum and minimum material properties, there is no big difference between the depths of CFS and composite failure surface. The reason is obvious and not repeated herein. Figure 14 shows the relationship between safety factors (SF) of the slope and the depth of water zone under different failure modes e.g., circular or composite. In the Fig. 14, the safety factors corresponding to circular slip surface are consistently lower than those corresponding to composite slip surface, which is consistent with the observation in Fig. 12. This indicates that under the situation when water table exists, the slope is more exposed to landslide with composite slip surface.  The CFS is mostly to be affected by the mishmash of c and ϕ. It means that cohesion and friction parameters are resisting forces, decreasing these value leads to decrease in the value of FS, D, and de and slip surface distribution. Deep-seated or local failures may take place if the spatial (longitudinal) changeability of cohesion and friction is considered. This is worth observing that potential slope slips can take place in authenticity (e.g. [37],). This detail shows the significance of choosing the spatial changeability influence of strength properties in slope stability analysis [39]. These results indicate that under the circumstances when pore water pressure exists, the loose slope is more exposed to slip with critical failure surface.
Therefore, it can be concluded that the effect of strength parameters for a loose rock material are helpful to choose appropriate value of these input parameters (i.e., cohesion and friction) for slope reinforcement design. Meanwhile, the factors of safety (FS) at c = 20 KN/ m 2 and ϕ = 30° for CFS is obviously more than one. These values can be enhanced, based on the actual situation, using related rock reinforcement measures, such as anchor length.

Comparison of failure surface of slope by SLOPE/W and FLAC
In this sub-section, the SLOPE/W and FLAC 3D software was selected to compare the prediction of sliding surface of loose rock slope. The boundary and geometry conditions in FLAC 3D software were kept same as those in SLOPE/W software. The same material properties which are presented in Sect. 3.3 and other material constants were Poisson's ratio μ = 0.33 and Young's modulus E = 60 MPa. FLAC 3D uses strength reduction technique to search possible value of safety factor (FS). In this method, c and ϕ are reduced with altered values of SRF until slope slip take place as: where c f is reduced cohesion, φ f reduced internal friction angle and f is safety factor. With increasing reduction factor the displacement also increases [40]. The safety factor FS = 0.86, 1.02 and 1.2 obtained via SLOPE/W (shown in Fig. 10) is consistent with FS = 0.88, 1.06 and 1.16 attained via FLAC 3D (Fig. 15). Regarding the prediction of the possible sliding mechanism of slope, (11) C f = c∕f (12) f = arctan tan f the distribution of the shear strain attained via FLAC 3D is presented in Fig. 15. This can be clearly observed from the chart that the shear strain localization band take place firstly followed by the development of the clear sliding surface of the loose rock slope. As can be observed form the Fig. 15, the numerical findings are same as far as the size and shape of the sliding surface of the studied slope is considered, however the positions are considerably different. On the other hand, SLOPE/W has certain benefits over FLAC 3D with respect to material heterogeneity problems. It can be concluded that material heterogeneity problems can be captured practically well by the SLOPE/W software to avoid an extremely skewed mesh in FLAC 3D . For that reason, this recommends that SLOPE/W is more appropriate computer code to FLAC 3D when this comes to forecasting the sliding surface of heterogeneous material subjected to water zone.

Conclusions
In this study, a new search procedure based on the geometrical parameters is presented to estimate the length of failure arc. The failure mechanism of critical failure surface (CFS) and variability effect of strength parameters, and the depth of water zone on the distribution of CFS is examine by SLOPE/W software. The correctness of results is checked by FLAC 3D software. Results show that the rock slope formed by loose material can merely produce a global failure by ignoring the strength properties but with an increasing these properties the distribution range of CFSs increases. The pore water pressure (PWP) also plays a dominant role in altering the slope stability. Comparing the SFs and the locations of the failure surfaces, the landslides with composite slip surfaces are much easier to occur for loose rock slope subjected to water pressure.
The method is presented to estimate the length of sliding arc (L) provides the new researcher with a valuable superfluous procedure of analyzing loose rock slope stability. It is also noted that the current study underestimate the length of failure arc for loose rock slope due to the assumptions of sliding surface entry and exit points. The current geometrical method cannot estimate the value of FS, so this method should be modified in future research.