Probabilistic analysis of tunnel face seismic stability in layered rock masses using Polynomial Chaos Kriging metamodel

Face stability is an essential issue in tunnel design and construction.Layered rock masses are typical and ubiquitous;uncertainties in rock properties always exist.In view of this,a comprehensive method,which combines the Upper bound Limit analysis of Tunnel face stability,the Polynomial Chaos Kriging,the Monte-Carlo Simulation and Analysis of Covariance method(ULT-PCK-MA),is proposed to investigate the seismic stability of tunnel faces.A two-dimensional analytical model of ULT is developed to evaluate the virtual support force based on the upper bound limit analysis.An efficient probabilistic analysis method PCK-MA based on the adaptive Polynomial Chaos Kriging metamodel is then implemented to investigate the parameter uncertainty effects.Ten input parameters,including geological strength indices,uniaxial compressive strengths and constants for three rock formations,and the horizontal seismic coefficients,are treated as random variables.The effects of these parameter uncertainties on the failure probability and sensitivity indices are discussed.In addition,the effects of weak layer position,the middle layer thickness and quality,the tunnel diameter,the parameters correlation,and the seismic loadings are investigated,respectively.The results show that the layer distributions significantly influence the tunnel face probabilistic stability,particularly when the weak rock is present in the bottom layer.The efficiency of the proposed ULT-PCK-MA is validated,which is expected to facilitate the engineering design and construction.

Influence of joint spacing and rock characteristics on the toppling stability of cut rock slope through a simplified limit equilibrium method

Toppling failure of rock mass/soil slope is an important geological and environmental problem.Clarifying its failure mechanism under different conditions has great significance in engineering.The toppling failure of a cutting slope occurred in a hydropower station in Kyushu,Japan illustrates that the joint characteristic played a significant role in the occurrence of rock slope tipping failure.Thus,in order to consider the mechanical properties of jointed rock mass and the influence of geometric conditions,a simplified analytical approach based on the limit equilibrium method for modeling the flexural toppling of cut rock slopes is proposed to consider the influence of the mechanical properties and geometry condition of jointed rock mass.The theoretical solution is compared with the numerical solution taking Kyushu Hydropower Station in Japan as one case,and it is found that the theoretical solution obtained by the simplified analysis method is consistent with the numerical analytical solution,thus verifying the accuracy of the simplified method.Meanwhile,the Goodman-Bray approach conventionally used in engineering practice is improved according to the analytical results.The results show that the allowable slope angle may be obtained by the improved Goodman-Bray approach considering the joint spacing,the joint frictional angle and the tensile strength of rock mass together.

Evaluation of slope stability through rock mass classification and kinematic analysis of some major slopes along NH-1A from Ramban to Banihal, North Western Himalayas

The network of Himalayan roadways and highways connects some remote regions of valleys or hill slopes,which is vital for India’s socio-economic growth.Due to natural and artificial factors,frequency of slope instabilities along the networks has been increasing over last few decades.Assessment of stability of natural and artificial slopes due to construction of these connecting road networks is significant in safely executing these roads throughout the year.Several rock mass classification methods are generally used to assess the strength and deformability of rock mass.This study assesses slope stability along the NH-1A of Ramban district of North Western Himalayas.Various structurally and non-structurally controlled rock mass classification systems have been applied to assess the stability conditions of 14 slopes.For evaluating the stability of these slopes,kinematic analysis was performed along with geological strength index(GSI),rock mass rating(RMR),continuous slope mass rating(CoSMR),slope mass rating(SMR),and Q-slope in the present study.The SMR gives three slopes as completely unstable while CoSMR suggests four slopes as completely unstable.The stability of all slopes was also analyzed using a design chart under dynamic and static conditions by slope stability rating(SSR)for the factor of safety(FoS)of 1.2 and 1 respectively.Q-slope with probability of failure(PoF)1%gives two slopes as stable slopes.Stable slope angle has been determined based on the Q-slope safe angle equation and SSR design chart based on the FoS.The value ranges given by different empirical classifications were RMR(37-74),GSI(27.3-58.5),SMR(11-59),and CoSMR(3.39-74.56).Good relationship was found among RMR&SSR and RMR&GSI with correlation coefficient(R 2)value of 0.815 and 0.6866,respectively.Lastly,a comparative stability of all these slopes based on the above classification has been performed to identify the most critical slope along this road.

Borehole stability in naturally fractured rocks with drilling mud intrusion and associated fracture strength weakening:A coupled DFN-DEM approach

Borehole instability in naturally fractured rocks poses significant challenges to drilling.Drilling mud invades the surrounding formations through natural fractures under the difference between the wellbore pressure(P w)and pore pressure(P p)during drilling,which may cause wellbore instability.However,the weakening of fracture strength due to mud intrusion is not considered in most existing borehole stability analyses,which may yield significant errors and misleading predictions.In addition,only limited factors were analyzed,and the fracture distribution was oversimplified.In this paper,the impacts of mud intrusion and associated fracture strength weakening on borehole stability in fractured rocks under both isotropic and anisotropic stress states are investigated using a coupled DEM(distinct element method)and DFN(discrete fracture network)method.It provides estimates of the effect of fracture strength weakening,wellbore pressure,in situ stresses,and sealing efficiency on borehole stability.The results show that mud intrusion and weakening of fracture strength can damage the borehole.This is demonstrated by the large displacement around the borehole,shear displacement on natural fractures,and the generation of fracture at shear limit.Mud intrusion reduces the shear strength of the fracture surface and leads to shear failure,which explains that the increase in mud weight may worsen borehole stability during overbalanced drilling in fractured formations.A higher in situ stress anisotropy exerts a significant influence on the mechanism of shear failure distribution around the wellbore.Moreover,the effect of sealing natural fractures on maintaining borehole stability is verified in this study,and the increase in sealing efficiency reduces the radial invasion distance of drilling mud.This study provides a directly quantitative prediction method of borehole instability in naturally fractured formations,which can consider the discrete fracture network,mud intrusion,and associated weakening of fracture strength.The information provided by the numerical approach(e.g.displacement around the borehole,shear displacement on fracture,and fracture at shear limit)is helpful for managing wellbore stability and designing wellbore-strengthening operations.

Heterogeneous information phase space reconstruction and stability prediction of filling body–surrounding rock combination

Traditional research believes that the filling body can effectively control stress concentration while ignoring the problems of unknown stability and the complex and changeable stress distribution of the filling body–surrounding rock combination under high-stress conditions.Current monitoring data processing methods cannot fully consider the complexity of monitoring objects,the diversity of monitoring methods,and the dynamics of monitoring data.To solve this problem,this paper proposes a phase space reconstruction and stability prediction method to process heterogeneous information of backfill–surrounding rock combinations.The three-dimensional monitoring system of a large-area filling body–surrounding rock combination in Longshou Mine was constructed by using drilling stress,multipoint displacement meter,and inclinometer.Varied information,such as the stress and displacement of the filling body–surrounding rock combination,was continuously obtained.Combined with the average mutual information method and the false nearest neighbor point method,the phase space of the heterogeneous information of the filling body–surrounding rock combination was then constructed.In this paper,the distance between the phase point and its nearest point was used as the index evaluation distance to evaluate the stability of the filling body–surrounding rock combination.The evaluated distances(ED)revealed a high sensitivity to the stability of the filling body–surrounding rock combination.The new method was then applied to calculate the time series of historically ED for 12 measuring points located at Longshou Mine.The moments of mutation in these time series were at least 3 months ahead of the roadway return dates.In the ED prediction experiments,the autoregressive integrated moving average model showed a higher prediction accuracy than the deep learning models(long short-term memory and Transformer).Furthermore,the root-mean-square error distribution of the prediction results peaked at 0.26,thus outperforming the no-prediction method in 70%of the cases.

Stability prediction of hard rock pillar using support vector machine optimized by three metaheuristic algorithms

Hard rock pillar is one of the important structures in engineering design and excavation in underground mines.Accurate and convenient prediction of pillar stability is of great significance for underground space safety.This paper aims to develop hybrid support vector machine(SVM)models improved by three metaheuristic algorithms known as grey wolf optimizer(GWO),whale optimization algorithm(WOA)and sparrow search algorithm(SSA)for predicting the hard rock pillar stability.An integrated dataset containing 306 hard rock pillars was established to generate hybrid SVM models.Five parameters including pillar height,pillar width,ratio of pillar width to height,uniaxial compressive strength and pillar stress were set as input parameters.Two global indices,three local indices and the receiver operating characteristic(ROC)curve with the area under the ROC curve(AUC)were utilized to evaluate all hybrid models’performance.The results confirmed that the SSA-SVM model is the best prediction model with the highest values of all global indices and local indices.Nevertheless,the performance of the SSASVM model for predicting the unstable pillar(AUC:0.899)is not as good as those for stable(AUC:0.975)and failed pillars(AUC:0.990).To verify the effectiveness of the proposed models,5 field cases were investigated in a metal mine and other 5 cases were collected from several published works.The validation results indicated that the SSA-SVM model obtained a considerable accuracy,which means that the combination of SVM and metaheuristic algorithms is a feasible approach to predict the pillar stability.

A stability evaluation method for deep-seated toppling in the upper Lancang river,Southwestern China

Deep-seated toppling in the upper reaches of the Lancang River,southwest China involves deformations exceeding 100 m in depth.The slope deformation is initiated by river downcutting and evolves distinctive characteristics with a depth of river incision.In this study,we propose a system for evaluating the stability of deep-seated toppled slopes in different evolutionary stages.This system contains identification criteria for each evolutionary stage and provides the corresponding stability evaluation methods.Based on the mechanical and kinematic analysis of slope blocks,the specific stage of slope movement can be identified in the field through outcrop mapping,in situ tests,surface displacement monitoring,and adit and borehole explorations.The stability evaluation methods are established based on the limiting equilibrium theory and the strain compatibility between the undisturbed zone and the toppled zone.Finally,several sample slopes in different evolution stages have been investigated to verify the applicability and accuracy of the proposed stability evaluation system.The results indicate that intense tectonic activity and rapid river incision lead to a maximum principal stress ratio exceeding 10 near the slope surface,thus triggering widespread toppling deformations along the river valley.When considering the losses of joint cohesion during the further rotation process,the safety factor of the slope drops by 7%e28%.The self-stabilization of toppling deformation can be recognized by the layer symmetry configuration after the free rotation of the deflected layers.Intensely toppled rock blocks mainly suffer sliding failures beyond the layer symmetry condition.The factor of safety of the K73 rockslide decreased from 1.17 to 0.87 by considering the development of the potential sliding surface and the toesaturated zone.

Automatic identification of rock discontinuity and stability analysis of tunnel rock blocks using terrestrial laser scanning

Local geometric information and discontinuity features are key aspects of the analysis of the evolution and failure mechanisms of unstable rock blocks in rock tunnels.This study demonstrates the integration of terrestrial laser scanning(TLS)with distinct element method for rock mass characterization and stability analysis in tunnels.TLS records detailed geometric information of the surrounding rock mass by scanning and collecting the positions of millions of rock surface points without contact.By conducting a fuzzy K-means method,a discontinuity automatic identification algorithm was developed,and a method for obtaining the geometric parameters of discontinuities was proposed.This method permits the user to visually identify each discontinuity and acquire its spatial distribution features(e.g.occurrences,spac-ings,trace lengths)in great detail.Compared with hand mapping in conventional geotechnical surveys,the geometric information of discontinuities obtained by this approach is more accurate and the iden-tification is more efficient.Then,a discrete fracture network with the same statistical characteristics as the actual discontinuities was generated with the distinct element method,and a representative nu-merical model of the jointed surrounding rock mass was established.By means of numerical simulation,potential unstable rock blocks were assessed,and failure mechanisms were analyzed.This method was applied to detection and assessment of unstable rock blocks in the spillway and sand flushing tunnel of the Hongshiyan hydropower project after a collapse.The results show that the noncontact detection of blocks was more labor-saving with lower safety risks compared with manual surveys,and the stability assessment was more reliable since the numerical model built by this method was more consistent with the distribution characteristics of actual joints.This study can provide a reference for geological survey and unstable rock block hazard mitigation in tunnels subjected to complex geology and active rockfalls.

Stability mechanism and control of the pumpable supports in longwall recovery room

The load-bearing performance(LBP)of pumpable supports(PPS)is crucial for the stability of longwall pre-driven recovery room(PRR)surrounding rock.However,the unbalanced bearing coefficient(UBC)of the PPS(undertaking unequal load along the mining direction)has not been investigated.A mechanical model of the PRR was established,considering the main roof cantilever beam structure,to derive an assessment formula for the load,the failure criteria,and the UBC of the PPS.Subsequently,the generation mechanisms,and influencing factors of the UBC were revealed.Global sensitivity analysis shows that the main roof hanging length(l_(2))and the spacing between the PPS(r)significantly impact the UBC.A novel design of the PPS and the coupling control technology were proposed and applied to reduce the UBC of the PPS in the adjacent longwall PRR.Monitor results showed no failure of the PPS at the test site,with the UBC(ζ)reduced to 1.1 consistent with the design value(1.15)basically,fully utilizing the collaborative LBP of the PPS.Finally,the maximum roof-to-floor convergence of the PRR was 234 mm,effectively controlling the stability of the surrounding rock of the PRR and ensuring the mining equipment recovery.

Classification and rating of disintegrated dolomite strata for slope stability analysis

Although disintegrated dolomite,widely distributed across the globe,has conventionally been a focus of research in underground engineering,the issue of slope stability issues in disintegrated dolomite strata is gaining increasing prominence.This is primarily due to their unique properties,including low strength and loose structure.Current methods for evaluating slope stability,such as basic quality(BQ)and slope stability probability classification(SSPC),do not adequately account for the poor integrity and structural fragmentation characteristic of disintegrated dolomite.To address this challenge,an analysis of the applicability of the limit equilibrium method(LEM),BQ,and SSPC methods was conducted on eight disintegrated dolomite slopes located in Baoshan,Southwest China.However,conflicting results were obtained.Therefore,this paper introduces a novel method,SMRDDS,to provide rapid and accurate assessment of disintegrated dolomite slope stability.This method incorporates parameters such as disintegrated grade,joint state,groundwater conditions,and excavation methods.The findings reveal that six slopes exhibit stability,while two are considered partially unstable.Notably,the proposed method demonstrates a closer match with the actual conditions and is more time-efficient compared with the BQ and SSPC methods.However,due to the limited research on disintegrated dolomite slopes,the results of the SMRDDS method tend to be conservative as a safety precaution.In conclusion,the SMRDDS method can quickly evaluate the current situation of disintegrated dolomite slopes in the field.This contributes significantly to disaster risk reduction for disintegrated dolomite slopes.

Effect of surface retaining elements on rock stability:laboratory investigation with sand powder 3D printing

This study aims to investigate the benefcial efects of surface retaining elements (SREs) on the mechanical behaviors of bolted rock and roadway stability. 3D printing (3DP) technology is utilized to create rock analogue prismatic specimens for conducting this investigation. Uniaxial compression tests with acoustic emission (AE) and digital image correlation techniques have been conducted on 3DP specimens bolted with diferent SREs. The results demonstrate that the strength and modulus of elasticity of the bolted specimens show a positive correlation with the area of the SRE;the AE characteristics of the bolted specimens are higher than those of the unbolted specimen, but they decrease with an increase in SRE area, thus further improving the integrity of the bolted specimens. The reinforcement efect of SREs on the surrounding rock of roadways is further analyzed using numerical modelling and feld test. The results provide a better understanding of the role of SREs in rock bolting and the optimization of rock bolting design. Furthermore, they verify the feasibility of 3DP for rock analogues in rock mechanics tests.

Theoretical study on stability evolution of soft and hard interbedded bedding reservoir slopes

Soft and hard interbedded bedding rock slopes,which is prone to failure,are widely distributed in the Three Gorges Reservoir,China.Limit equilibrium method(LEM)is commonly used to analyze the stability of bedding rock slopes that have a single failure plane.However,this method cannot accurately estimate the stability of soft and hard interbedded bedding reservoir slopes because the strength parameters of a soft and hard interbedded rock mass vary spatially along the bedding plane and deteriorate with time due to periodic fluctuations of reservoir level.A modified LEM is proposed to evaluate the stability evolution of soft and hard interbedded bedding reservoir slopes considering the spatial variation and temporal deterioration of shear strength parameters of rock masses and bedding planes.In the modified LEM,the S-curve model is used to define the spatial variation of shear strength parameters,and general deterioration equations of shear strength parameters with the increasing number of wettingdrying cycles(WDC)are proposed to describe the temporal deterioration.Also,this method is applied to evaluate the stability evolution of a soft and hard interbedded bedding reservoir slope,located at the Three Gorges Reservoir.The results show that neglecting the spatial variation and temporal deterioration of shear strength parameters may overestimate slope stability.Finally,the modified LEM provides useful guidance to reasonably evaluate the long-term stability of soft and hard interbedded bedding reservoir slopes in reservoir area.

Probabilistic assessment of rock slopes stability using the response surface approach——A case study

This paper examines the stability condition of a jointed rock slope in the south western region of Saudi Arabia using deterministic and probabilistic approaches,under both dry and wet conditions.The study area is characterized by complex geology in rugged terrains.The stability analysis is carried out using the code FLAC3D to generate a 3-dimensional,ubiquitous joint model,to determine the influence of the dominant,unfavourable discontinuity orientation with respect to the slope face.The deterministic analysis is first implemented using the mean values of the selected random variables,namely the dip,dip direction and friction angle of the dominant discontinuity set,and the stability condition is assessed with a factor of safety based on the classical frictional joint constitutive model.A Box-Behnken design(BBD)approach is then adopted to create the surface response function as a second order polynomial for the factor of safety.To do so,fifteen FLAC3D models are generated in accordance with the BBD.Based on this,10,000 simulations of different slope realizations are carried out using Monte-Carlo simulation technique,and the probability of unsatisfactory of performance of the rock slope is assessed.It is shown that the probabilistic approach provides more insight and confidence in the stability condition of the rock slope,both under dry and steady state heavy rainfall conditions.A discussion is presented on the significance of accepting lower safety factors when heavy rainfall conditions are encountered.

Microseismic monitoring and numerical simulation on the stability of high-steep rock slopes in hydropower engineering

For high-steep slopes in hydropower engineering, damage can be induced or accumulated due to a seriesof human or natural activities, including excavation, dam construction, earthquake, rainstorm, rapid riseor drop of water level in the service lifetime of slopes. According to the concept that the progressivedamage （microseismicity） of rock slope is the essence of the precursor of slope instability, a microseismicmonitoring system for high-steep rock slopes is established. Positioning accuracy of the monitoringsystem is tested by fixed-position blasting method. Based on waveform and cluster analyses of microseismicevents recorded during test, the tempo-spatial distribution of microseismic events is analyzed.The deformation zone in the deep rock masses induced by the microseismic events is preliminarilydelimited. Based on the physical information measured by in situ microseismic monitoring, an evaluationmethod for the dynamic stability of rock slopes is proposed and preliminarily implemented bycombining microseismic monitoring and numerical modeling. Based on the rock mass damage modelobtained by back analysis of microseismic information, the rock mass elements within the microseismicdamage zone are automatically searched by finite element program. Then the stiffness and strengthreductions are performed on these damaged elements accordingly. Attempts are made to establish thecorrelation between microseismic event, strength deterioration and slope dynamic instability, so as toquantitatively evaluate the dynamic stability of slope. The case studies about two practical slopes indicatethat the proposed method can reflect the factor of safety of rock slope more objectively. Numericalanalysis can help to understand the characteristics and modes of the monitored microseismic events inrock slopes. Microseismic monitoring data and simulation results can be used to mutually modify thesensitive rock parameters and calibrate the model. Combination of microseismic monitoring and numericalsimulation provides a more objective basis for the numerical model and parameters and a solidmechanical foundation for the microseismic monitoring.

Application of rock mass classification systems to rock slope stability assessment:A case study

The stability of rock slopes is considered crucial to public safety in highways passing through rock cuts, as well as to personnel and equipment safety in open pit mines. Slope instability and failures occur due to many factors such as adverse slope geometries, geological discontinuities, weak or weathered slope materials as well as severe weather conditions. External loads like heavy precipitation and seismicity could play a significant role in slope failure. In this paper, several rock mass classification systems developed for rock slope stability assessment are evaluated against known rock slope conditions in a region of Saudi Arabia, where slopes located in rugged terrains with complex geometry serve as highway road cuts. Selected empirical methods have been applied to 22 rock cuts that are selected based on their failure mechanisms and slope materials. The stability conditions are identified, and the results of each rock slope classification system are compared. The paper also highlights the limitations of the empirical classification methods used in the study and proposes future research directions.

Deformation Mechanism and Stability of a Rocky Slope

A high slope is located on the side of the spillway at a hydropower station in Southwest China, which has some weak inter-layers inclining outwards. Parts of the slope show heavy weathering and unloading. There appeared deformation and tensile crack either on the surface or on the afteredge of the slope during excavation, and under a platform （elev. 488 m）, two levels of slopes collapsed on the downriver side. Based on the investigation in situ and the analysis of the geological structure, the conceptual model of deformation and failure mechanism was erected for this slope. Furthermore, the deformation characteristics were studied with FLAC^3D numerical simulation. Comprehensive analysis shows that the whole deformation of the slope is unloading rebound in certain depth scope and the whole body does not slide along any weak interlayer. In addition, two parts with prominent local deformation in the shallow layer of the slope show the models of ＂creep sfiding-tensile cracking＂ and ＂slidlng-tensile cracking＂, respectively. Based on the above analysis, the corresponding project of support and reinforcement is proposed to make the slope more stable.

Model test and numerical simulation on the dynamic stability of the bedding rock slope under frequent microseisms

Shake table testing was performed to investigate the dynamic stability of a mid-dip bedding rock slope under frequent earthquakes. Then, numerical modelling was established to further study the slope dynamic stability under purely microseisms and the influence of five factors, including seismic amplitude, slope height, slope angle, strata inclination and strata thickness, were considered. The experimental results show that the natural frequency of the slope decreases and damping ratio increases as the earthquake loading times increase. The dynamic strength reduction method is adopted for the stability evaluation of the bedding rock slope in numerical simulation, and the slope stability decreases with the increase of seismic amplitude, increase of slope height, reduction of strata thickness and increase of slope angle. The failure mode of a mid-dip bedding rock slope in the shaking table test is integral slipping along the bedding surface with dipping tensile cracks at the slope rear edge going through the bedding surfaces. In the numerical simulation, the long-term stability of a mid-dip bedding slope is worst under frequent microseisms and the slope is at risk of integral sliding instability, whereas the slope rock mass is more broken than shown in the shaking table test. The research results are of practical significance to better understand the formation mechanism of reservoir landslides and prevent future landslide disasters.

Stability classification model of mine-lane surrounding rock based on distance discriminant analysis method

Based on the principle of Mahalanobis distance discriminant analysis (DDA) theory, a stability classification model for mine-lane surrounding rock was established, including six indexes of discriminant factors that reflect the engineering quality of surrounding rock: lane depth below surface, span of lane, ratio of directly top layer thickness to coal thickness, uniaxial comprehensive strength of surrounding rock, development degree coefficient of surrounding rock joint and range of broken surrounding rock zone. A DDA model was obtained through training 15 practical measuring samples. The re-substitution method was introduced to verify the stability of DDA model and the ratio of mis-discrimination is zero. The DDA model was used to discriminate 3 new samples and the results are identical with actual rock kind. Compared with the artificial neural network method and support vector mechanic method, the results show that this model has high prediction accuracy and can be used in practical engineering.

Seismic stability of jointed rock slopes under obliquely incident earthquake waves

Seismic stability of slopes has been traditionally analyzed with vertically propagated earthquake waves.However,for rock slopes,the earthquake waves might approach the outcrop still with a evidently oblique direction.To investigate the impact of obliquely incident earthquake excitations,the input method for SV and P waves with arbitrary incident angles is conducted,respectively,by adopting the equivalent nodal force method together with a viscous-spring boundary.Then,the input method is introduced within the framework of ABAQUS software and verified by a numerical example.Both SV and P waves input are considered herein for a 2 D jointed rock slope.For the jointed rock mass,the jointed material model in ABAQUS software is employed to simulate its behavior as a continuum.Results of the study show that the earthquake incident angles have significance on the seismic stability of jointed rock slopes.The larger the incident angle,the greater the risk of slope instability.Furthermore,the stability of the jointed rock slopes also is affected by wave types of earthquakes heavily.P waves induce weaker responses and SV waves are shown to be more critical.

Assessment of rock slope stability by probabilistic-based Slope Stability Probability Classification method along highway cut slopes in Adilcevaz-Bitlis(Turkey)

Rock slope stability is of great concern along highway routes as stability problems on cut slopes may cause fatal events as well as loss of property.In rock slope engineering,stability evaluations are commonly performed by means of analytical or numerical analyses,principally considering the factor of safety concept.As a matter of fact,the probabilistic assessment of slope stability is progressively getting popularity due to difficulties in assigning the most appropriate values to design parameters in analytical or numerical methods.Additionally,the effect of heterogeneities in rock masses and discontinuities on the analysis results is minimized through the probabilistic concept.In this study,slope stability of high and steep sedimentary rock cut slopes along a state highway in AdilcevazBitlis(Turkey) was evaluated on the basis of probabilistic approach using the Slope Stability Probability Classification(SSPC) system.The probabilistic assessment indicates major slope stability problems because of discontinuity controlled and discontinuity orientation independent mass movements.Almost all studied cut slopes suffer from orientation-independent stability problems with very low stability probabilities.Additionally,the probability of planar and toppling failures issignificantly high with respect to the SSPC system.The stability problems along the investigated rock slopes were also verified by field reconnaissance.Remedial measures such as slope re-design and reinforcement at the studied locations should be taken to prevent hazardous events along the highway.On the other hand,the probabilistic approach may be a useful tool during rock slope engineering to overcome numerous uncertainties when probabilistic and analytic results are compared.