Brief Literature Review and Classification System of Reliability Methods for Evaluating the Stability of Earth Slopes

: The issue of slope stability is one of the most important and yet most difficult geotechnical problems. Assessing slope stability is particularly difficult because of the many uncertainties in ‐ volved in the process. To take these uncertainties into account, probabilistic methods are used, and the reliability approach is adopted. There are many methods for reliability assessment of earth slope stability. However, there is no system that would organize all of these methods in an unambiguous way. In fact, these methods can be classified in different ways: by assignment to a deterministic classification of methods, by description of uncertainties of soil parameters, by level of reliability according to the theory of reliability, etc. The huge number of articles summarizing the research in this field, but in various “disordered” directions, certainly do not facilitate the understanding or ultimately the practical application of the reliability approach by the engineer. The paper proposes a universal classification system of reliability methods for evaluating the stability of earth slopes. This proposal is preceded by a brief literature review of both historical background and contempo ‐ rary research on reliability analysis of earth slope stability.


Introduction
The assessment of slope stability in terms of reliability analysis is still an active research topic being investigated by many scientific centres around the world. In only the last dozen years or so, in the leading geological and geotechnical journals as well as proceedings summarising international conferences, several thousand publications have appeared on this subject. Research is being conducted in different directions covering a wide spectrum of problems related to the deterministic and probabilistic methods used, comparative analyses, modelling of uncertain soil parameters (random variable, random field), influence of randomness of various factors affecting slope stability, types of problem, acting loading, calculation procedures with the possibility of simplifying these calculations, etc. In the present form, such a huge number of articles, of which only some contribute to the development of the discipline, is a disorderly "multidirectional" collection. It seems obvious that it should be systematised, because even for most researchers specialising in this subject, the quantity of disordered information can lead to confusion and possibly even incorrect conclusions. Of course, such a collection still would not give the engineer a universal tool for work. On the contrary, it would consistently discourage them from reliably analysing slope stability. In fact, the unwillingness of engineers to apply reliability to slope analysis and other geotechnical problems results from the lack of basic knowledge of statistics and probability theory or from common misconceptions with respect to the requirements regarding probabilistic methods, especially from the insufficiency of convincing literature illustrating the implementation and benefits of such analyses. To facilitate putting the reliability approach to slope stability into practice, the applied methodologies and procedures should be simple, well-known to engineers and able to solve real slope problems. In order to achieve this, it is first necessary to set in order and systematise the extensive literature on the subject, comprehend its particular aspects and then classify it. Also, some concepts and terms, not only those appearing in recent literature, should be clarified. Thus, a brief overview of earth slope stability reliability approaches starting from the historical background and moving to contemporary research is presented in this paper. In the latter, apart from the methods of analysis, particular attention is paid to the issues developed in the last dozen or so years such as system reliability and description of random soil properties, especially spatial variability. Some comments on applications of reliability approaches to slope stability analysis are discussed as well. The state of research in the subject presented in an orderly manner, as well as the explanation of a number of terms and concepts occurring here, will bring the engineer closer to understanding this difficult issue and will contribute to the wider application of the reliability approach in geotechnical practice.
This paper proposes a classification system of reliability methods for earth slope stability assessment that integrates deterministic slope stability methods, modelling of uncertain soil parameters and reliability level (commonly used in the structural safety analysis). Also, in the case of the most sophisticated approaches, further divisions related to improvement of computation efficiency are suggested.

Historical Background
In general, there are two types of analyses for evaluating the stability of slopes against failure: deterministic and probabilistic (reliability). In the first case, three basic approaches have been developed: the limit equilibrium method (LEM), the displacementbased finite element method (FEM) used by the strength reduction method (SRM) (also called shear strength reduction (SSR)) and limit analysis (LA) by lower and upper bound solutions theorems. The limit equilibrium method utilises an assumed slip surface and determines its static equilibrium, usually by discretising the assumed failing soil mass into slices. The forces are then summed for each slice, creating a statically determinate problem following some assumptions. By introducing the factor of safety for the entire sliding mass, global equilibrium is maintained for a system on the verge of failure. The LEM is the oldest technique for evaluating slope stability and the most commonly used in practice. However, it is restricted by its arbitrary choice of failure mechanism and by interslice forces. Several variants of this method have been proposed [1][2][3][4][5][6]. In the FEM by SRM used in evaluating slope stability, soil strength parameters continuously decrease until the first indications of failure appear. The safety factor is defined as the ratio of the real shear strength of the soil to the reduced shear strength. This method seems to be superior to the LEM because there is no need for the primary guess at determining the critical failure surface. In addition, this method does not require any assumptions about interslice forces. FEM analysis is a more rigorous and universal technique, but often less attractive due to its dependence on mesh density and the available computational capacity. The primary advantage of this method is that the critical slip surface is found automatically from the shear strain, which increases as the shear strength decreases. Unfortunately, other ''slip'' surfaces (i.e., local minima) are omitted. However, because of the high speed of modern computer systems, analysis by FEM is used today more often than before. The finite element method by strength reduction method was first proposed and applied to slope stability by Zienkiewicz et al. [7] and then used to assess slope stability in, among others, [8][9][10]. Limit analysis models soil as a material that is perfectly plastic and obeys an associated flow rule. This method employs a dichotomy of theorems to provide a solution: either upper bound or lower bound plasticity. The upper bound theorem of limit analysis is predominantly used in solving slope stability problems. Unfortunately, the application of LA is still limited, since most of the research findings are chart-based and prepared for particular cases, and there is no stability chart available to cover a wide range of different slope material properties, geometries, etc. The concept of limit analysis was proposed by Drucker and Prager [11] and was utilised in slope stability in [12][13][14] and others. A review of the three basic deterministic approaches of slope stability, including their shortcomings and possible errors, is discussed in [15].
In general, the probabilistic methods used in geotechnical engineering can be divided into two groups, depending on the description of the uncertainties of the soil parameters: the single random variable (SRV) approach and the random field (RF) approach. The former method is commonly used in practice due to its concise concepts and simplicity in analysis. The latter reliably analyses the spatial variability of the soil properties. Regarding slope stability, these methods make use both of the LEM and FEM and rarely LA. Interest in the analysis of slope stability using the reliability approach began over 50 years ago. The number of papers on this subject has increased considerably since 1975, when a second ICASP International Congress was held in Aachen. After that, the correction factor methods, the first-order second-moment (FOSM), second-order second-moment (SOSM) and Monte Carlo (MC) methods, dominated. However, the majority of probabilistic or reliability methods made use of traditional slope stability analysis techniques, i.e., LEM [16][17][18][19][20][21][22][23][24][25]. In the case of a highly nonlinear function of factor of safety, computations of the derivatives are impossible or inconvenient, thus rendering FOSM results inaccurate. In addition, different results can be obtained depending on how the limit state function is formulated. To avoid the main drawback of the FOSM and SOSM methods (in which the results unduly depend on the mean value), as observed in [26,27], the first-order reliability method (FORM) has begun to be widely used [28][29][30][31]. It is well known that FORM works only for slopes with a small probability of failure or a high reliability index. Otherwise, this method underestimates the probability of failure of slopes. A summary of research on probabilistic analysis of slope stability was given in the monograph [32]. A review of the literature on this topic can also be found in [24] or [33].

Contemporary Research
Recently, numerical methods have found particularly significant value in reliability analysis, mainly due to the rapid development of computerisation. Their use, however, was associated with further research problems. Some of them are briefly presented below and others are only mentioned.

Methodology
A significant change in approach to probabilistic slope analysis occurred with the use of the finite element method (FEM) for geotechnical problems. Different probabilistic methods related to FEM have been proposed, such as the perturbation method, the Neumann expansion method, the partial differential method, the spectral stochastic finite element method (SSFEM), etc. Along with the development of computers and software, the MC simulation became dominant because of its relative ease of application. This method is a conceptually simple tool for reliability analysis of slope stability regardless of the form of the performance function or the number of scenario failure events. It employs statistical averaging over random samples generated from the probability density function of the parameters to evaluate the probability of failure. It is the easiest to apply; however, its simulations are usually time consuming and computation demanding. The MC method is also robust to various deterministic analysis methods for slope stability analysis, such as LEM or FEM/FDM [30,53]. Ali et al. [54] combined finite element limit analysis with random fields to deliver a valuable tool for probabilistic analysis of stability problems. Using a two-layer soil slope as an example, they showed that their proposed approach performed much better than equivalent uniform meshes (used in deterministic analysis) in reducing the gap between the upper and lower bounds of probabilistic solutions. The random finite element/difference method (RFEM/RFDM) was developed to account for the effect of spatial variability [25,55]. This method combines random field theory with the finite element/difference method in a Monte Carlo framework. However, it does not offer insight into the relative contributions of various uncertainties to the failure probability. This is of particular interest in engineering practice. The critical slip surface varies spatially and needs to be located for each random sample generated during the MCS. The RFEM/RFDM also suffers from a lack of efficiency at small probability levels, especially in spatially variable soils. A comparative analysis of different variants of probabilistic and reliability methods based on FEM and applied to slope reliability assessment problems was given in [55][56][57][58].
Unfortunately, in the case of complicated slopes, RFEM requires extensive computational effort. In order to improve the computation efficiency of slope reliability analysis, research has been conducted in two basic directions: simplifying the performance function of slope stability and reducing the sample size of Monte Carlo simulations. In the first case, the response surface method [59], Kriging methodology [60], surrogate models [61], artificial neural network [62], support vector machine [63] or genetic algorithms [64] are usually utilised. The latter includes advanced sampling techniques such as subset simulation (SS), importance sampling (IS), etc. For example, the importance sampling probability density function is much closer to the failure region; thus, fewer realisations are needed for slopes with a low probability of failure. On this basis, Wang et al. [53] incorporated the SS method and Ching et al. [65] the IS method into the LEM. Li et al. [66] proposed the deaggregation approach, whereas Li et al. [67] combined MCS with both the LEM (in a first step) and the FEM. Hung et al. [68] combined method SS with the RFEM and Li et al. [69] proposed an advanced MCS method called generalised subset simulation (GSS). In GSS, the system failure event is decomposed into a number of scenario failure events through fault tree analysis, and the system failure probability and scenario failure probabilities are calculated by a single run of GSS. Thus, repeatedly performing the simulation for different failure events becomes unnecessary. This is of special importance when the number of slope failure events is large. An efficient machine learning (ML)aided stochastic reliability analysis technique for spatially variable slopes was introduced in [70]. An attempt to classify the finite element method for slope stability analysis by comparing random fields was made in [71,72].
Recently, the material point method (MPM), which is a variant of the finite element method, has been applied in the reliability analysis of earth slopes. This method has been shown to be a robust spatial discretisation method for simulating multiphase interactions. It has distinct advantages in solving extremely large deformation problems. Wang et al. [73] proposed the random material point method (RMPM), which combines random field theory and the MPM. It differs from the random finite-element method (RFEM) by assigning random field (cell) values to material points that are free to move relative to the computational grid rather than to Gauss points in a conventional finite-element mesh. The authors demonstrated the method for an idealised slope in a strain softening soil. The occurrence and evolution of various slope failure modes during large deformation in spatially variable soils using Monte Carlo simulation combined with the MPM was investigated in [74].
Most of the research regarding the reliability analysis of slope stability (in the case of simplifying the performance function) has been devoted to the RSM. The fundamental concept of this method is the approximation of the relationship between the factor of safety and soil model parameters (such as cohesion and friction angle) by a simple and explicit polynomial function (approximation model) or by artificial neural networks. In problems of slope stability, nonlinear problems dominate. Therefore, the quadratic response surface method is most popular for more complicated system reliability analyses. The RSM was introduced into the reliability analysis of slope stability, with implicit performance functions, by Wong [75]. Several papers have been devoted to this problem [55,59,[76][77][78][79][80][81][82][83]. The RSM may be applicable to some relatively simple cases, but such an approximation may not be possible for complicated multilayered slopes. The stratified response surface method was proposed in [84] in order to tackle the complexity problem. In this method, the number of stratified response surfaces is related to the number of soil layers that have uncertain soil parameters. A review of RSM can be found in [85]. A comparative analysis on the computational efficiency and accuracy of four commonly-used RSM approaches was performed in [67] for several cases of slope reliability problems with or without taking the spatial variability of soil into account.
In the past few decades, the point estimate method (PEM) proposed by Rosenblueth [86] has found increasing use in many engineering fields. There are some limitations to this proposal for multiple variables, but a number of researchers have modified the procedure to optimise computational time and accuracy in these cases [87,88]. In fact, PEM has become one of the most popular methods in practical geotechnical reliability analysis. Many authors have applied this method for probabilistic analysis of slope design, indicating its efficiency and accuracy [89][90][91][92][93].
In recent years, various modifications of the "classical" reliability methods as well as their combinations have been developed. For example, it is worth mentioning a novel simplified method, the so-called quantile-based first-order second-moment method (QFOSM). This method has a level of simplicity close to the FOSM and yet has accuracy close to the FORM. It does not require performing large number of performance function evaluations (as in MCS, IS) nor operating in the standard normal space (as in FORM/SORM). The search for the critical quantile position is a one-dimensional root finding problem, which can be readily implemented [94]. A second-order orthogonal experimental design, SOED-based RSM, was proposed in [78]. The SOED is constructed by changing the length of star points, and the main characteristic of the SOED is that the design matrix is diagonal. The so-called combined response surface Method (CRSM), a simple, straightforward hybrid computational procedure incorporating the response surface method (RSM) linked with standard Monte Carlo (MC) simulations and the point estimate method (PEM), was proposed by Winkelmann et al. [95]. In order to address the multicollinearity existing in slope reliability analysis, a modified version of the RSM utilising the least absolute shrinkage and selection operator (Lasso), ridge regression, elastic net regression and stepwise regression was given in [96]. A novel multidimensional cloud model coupled with connection numbers theory for evaluation of slope stability was proposed in [97]. The cloud model is a tool for depicting fuzziness and randomness of evaluation indicators in a unified way. A modified one-dimensional conditional Markov chain model describing stratigraphic boundary uncertainty and MCS was introduced by [98].

System Reliability
There may be many possible slip surfaces in a slope reliability analysis. However, two basic types of failure probability should be distinguished: the failure probability of the slope along an individual slip surface and the overall failure probability of the slope. The critical deterministic surface is one of many failure surfaces for which the safety factor is minimal and the reliability index is calculated. However, this index is not necessarily the minimum value. The minimum reliability index occurs on some critical probabilistic surface that does not, in general, coincide with the critical deterministic one. Locating the critical probabilistic surface may require additional computations. The critical slip surface with the minimum safety factor varies spatially when spatial variability is considered. Identification of that surface among a large number of potential slip surfaces is an elementary step in slope stability analysis. In real soil slopes, various slope failure modes (i.e., slip surfaces) are caused by stratification (i.e., layered soils). In general, there may exist many potential failure modes or slip surfaces. Thus, it is justified to analyse the slope reliability problem rather in the framework of system reliability. In such a case, each potential slip surface is a component, and the critical slip surface is the weakest component. System failure takes place when a landslide occurs along the critical slip surface. Slope reliability analysis is then defined as a system reliability analysis. It is obvious that the overall failure probability is greater than the failure probability of any individual potential slip surface because of system effects. The difference depends on the correlation between the failure probabilities of the different slip surfaces, associated with various slip surfaces, for which no general formulation is available. Although a slope may have many potential slip surfaces, its system failure probability is managed by only a few critical slip surfaces, which are called representative slip surfaces (RSS). The system reliability can be conveniently calculated based on these slip surfaces if the RSS are identified. An important aspect of reliability slope stability analysis is identifying the RSS. The failure probability of a system is generally more difficult to calculate than that of a single mode. First, instead of calculating the system failure probability directly, several researchers assessed the bounds of system failure probability [99,100]. However, these bounds could be wide when the single-mode failure probabilities are all large. The system reliability of slope stability has found much greater use along with the increasing development of simulation methods. These methods have allowed system failure probability to be predicted more accurately. For example, Ching et al. [65] applied the IS method, Huang et al. [101] combined the MCS method with the FEM, Wang et al. [53] used the SS and Zhang et al. [36] utilised the Hassan and Wolff method, which identifies the most critical slip surface of soil slopes based on sensitivity analysis, to assess the system reliability of a slope. It is also worth mentioning the papers [84] and [31], where a local search method and a barrier method for identification of RSS were developed, respectively. A system reliability analysis approach applied for layered soil slopes based on multivariate adaptive regression splines and MCS was given in [102]. A relatively new method to identify the RSS of arbitrary shapes based on the shear strength reduction method, for system reliability of soil slopes, was proposed in [103]. In that paper, a Kriging-based response surface to approximate the deterministic slope stability model was used. An efficient RSM-based MCS for risk assessment considering multiple failure surfaces was proposed in [104], whereas its extension by the effects of the stratigraphic boundary uncertainty was given in [105]. Duan et al. [106], viewing the task of representative slip surfaces as a multi-design-point identification problem, proposed a barrier-based optimisation method based on the shear strength reduction method to identify the RSS of arbitrary shapes.

Uncertainties of Soil Parameters
The uncertainties of the soil properties are a major contributor to the uncertainty in the stability of slopes. They comprise both a random measurement error (a result of a statistical uncertainty in the expected values and the effects of bias) and real spatial variability (variations in the trend of the parameters). Research on the influence of systematic error of soil parameters on slope stability has been carried out almost from the beginning of the application of probabilistic methods [24,26,93,[107][108][109][110][111]. The copula theory to evaluate slope reliability in the presence of incomplete probability information was given in [112][113][114]. Application of a coupled Markov chain model that considered geological uncertainty for slope reliability was introduced in [115]. That model was also used in [116], where the inherent variability of the soil parameters was included. Among the available extensive literature on slope analysis, a significant part concerns hypothetical cases, usually simple slopes with a single soil layer or with two soil layers. The spatial variability of soil properties is one of the main sources of geotechnical uncertainty and significantly affects the reliability of slope stability. Thus, it is important that the spatially varying soil properties should be modelled appropriately. However, in such a case, the computing time consumed can be significantly longer. Ignoring this variability leads to an overestimation of the variance of the slope safety factor that can result either under-or overestimating the probability of slope failure [25,53,117,118]. Research on the influence of spatial variability of soil properties on the probability of slope failure has already begun; it was addressed by Vanmarcke [20], and a significant contribution was made in [119]. The effect of the spatial variability of soil properties on slope stability has been studied using the LEM in [24,53,120,121] and using the RFEM in [55,76,101,122,123]. The influence of spatially correlated soil properties on the probability of failure was analysed in [124]. The cross-correlation among Markov random fields leading to a generic approach for modelling multivariate cross-correlated geotechnical random fields based on vine copulas was proposed in [114]. Griffith et al. [55] showed that the probability of failure of slopes computed by RFEM may differ significantly from that obtained using a random field and the LEM. The RSM was used for this purpose in [59]. Application of the RFEM requires an understanding the relationships between correlation length, mesh coarseness, and the necessary number of simulations. Fundamental aspects in the employment of random fields in the numerical analysis of geotechnical structures were investigated in [125]. Application of Bayesian analysis to account for spatially variable soil in the reliability slope approach was proposed in [126]. Simplified approaches to incorporating the spatial variability of soil properties in the reliability analysis of slope stability should also be noted. Li et al. [127] used equivalent homogeneous random soil parameters in the form of single random variables applied to the RFDM, whereas Liu et al. [128] utilised a simplified framework based on the multiple response surface method (MRSM) and Monte Carlo simulation. It is also worth mentioning a simplified approach [129], proposing a method that combines the advantages of the sequential search and bisection search of the strength reduction method (SRM).

Applications
Use of reliability slope stability analysis in practical applications is still limited. Most engineers are not familiar with probabilistic concepts. The methods used here are unknown or seem to be too difficult for them. Even if they are not, the data needed for a detailed statistical evaluation of the soil properties is usually unavailable. There are too few papers presenting applications of the reliability approach in a comprehensible and simple way.
Practical spreadsheet techniques in the Microsoft Excel software for a probabilistic slope stability analysis that obtains the same reliability index as FORM, based on Janbu's generalised procedure of slices, were proposed by Low and Tang [28]. Those authors also modified spreadsheets by including cases with correlated non-normals and explicit and implicit performance functions [130]. The same technique was applied to the generalised Morgenstern-Price method of slices [29]. Furthermore, a spreadsheet approach based on Monte Carlo simulation and using Bishop's method, including spatial variability of the input variables, was proposed in [24]. Wang et al. [53] implemented an advanced MCS method called subset simulation. Low and Phoon [131] illustrated geotechnical reliabilitybased analyses, among others for a soil slope, using the Excel spreadsheet platform. They focused their attention on practical procedures available in commercial FORM/SORM packages and suggested such reliability-based analyses to fill a complementary role to the Eurocode 7 design approach. An iterative algorithm for FORM analysis involving correlated non-normal variables and their spatial variability was proposed in [123]. The authors demonstrated its usefulness for practicing engineers for geotechnical reliability analysis where deterministic software is used. A logical framework for the load and resistance factor design (LRFD) of slopes based on reliability analysis was proposed in [110]. MCS has been adopted into commercial software packages such as SVSlope, FLAC, PLAXIS and GEO5 [64]. In addition, Cao et al. [111] presented a practical approach of reliability slope stability analysis that implements an advanced MCS method called "subset simulation" in a spreadsheet environment.
There are also simplified probabilistic slope stability design charts. For example, the chart for purely cohesive soil [25] and for cohesive-frictional soil [109].

Classification System
In the case of rock slopes, the slope stability probability classifications (SSPC) proposed by Hack et al. [132] is commonly used. It is based on a three-step approach and on the probabilistic assessment of independently different failure mechanisms in a slope. It is also known as the modified Hack classification system [133] and was applied in, among others, [134,135]. Unfortunately, there is still no such system for the reliability stability methods of earth slopes. Thus, a classification proposal is presented below.
In a structural safety analysis, there are different levels of reliability, depending on the importance of the structure, grouped under four basic levels [136]:  level 1-partial factor approach-employs only one "characteristic" value of each uncertainty parameter;  level 2-estimates two values of each uncertainty parameter, usually the mean value, standard deviation, and the correlation between these parameters;  level 3-best estimate of the probability of failure-knowledge of the join distribution of all uncertain parameters is required;  level 4-reliability methods appropriate for structures of major economic importance, taking into account the structures' economic value, including the consequences of their failure.
The level 2 reliability methods are included in limit state design codes [137]. Methods of level 2 include a range of approximate or iterative procedures such as the perturbation method, FOSM, SOSM, FORM, SORM, PEM, etc. Methods of level 3, in the strict sense, require determining the mathematically exact probability of structural failure as a result of integrating the joint probability density function of random variables. In the case of slope stability, they can only be used for simplified, idealised cases. However, in a broader sense, level 3 methods require estimates of all probabilistic measures. The Monte Carlo simulation method has become the dominant procedure here as a result of the rapid development of computer techniques that have taken place in recent years. In order to improve the efficiency of this method while maintaining the accuracy of calculations, various reduction techniques have been developed (e.g., stratified sampling, Latin Hypercube simulation, importance sampling and Russian roulette and splitting). The response surface method and the methods of artificial neural networks have also grown in popularity in the analysis of the reliability of slopes. Both methods allow all probabilistic measures to be estimated and can also be qualified to the level 3 reliability method.
In the deterministic approach to slope stability, three basic approaches-the limit equilibrium method (LEM), the displacement-based finite/different element method (FEM/DEM), and limit analysis (LA)-have been developed. Probabilistic methods used in geotechnical engineering are commonly divided into two groups, depending on the description of the uncertainties of the soil parameters: the single random variable (SRV) approach and the random field (RF) approach.
The proposal of a classification of reliability methods of slope stability is presented the Table 1. It combines deterministic methods of slope stability with random modelling of the soil medium and includes levels of reliability. The abbreviations in Table 1 refer to the random soil model, deterministic slope stability method and reliability level. For example: SRVLEM1-soil modelled as a single random variable (SRV), limit equilibrium method (LEM) and Level 1 reliability method; RFFEM3-soil modelled as a random field, finite/different element method and Level 3 reliability method; etc.
Currently, the FEM/DEM is the predominant deterministic method of stability assessment, soil is usually modelled as a random field and level 3 reliability methods are usually applied. Most of the research focus here is on computation efficiency. Thus, the RFFEM3 methods can be divided into two groups depending on how the performance function of slope stability is simplified and on the reduction techniques ( Figure 1). Slopes can be subjected both to static (gravity) and dynamic (earthquake, waves) loads as well as environmental loads (rainfall, temperature changes). Also, slope failure modelling can be carried out in a two-dimensional or three-dimensional analysis. However, these factors should not affect the proposed classification system of reliability methods. Instead, new subdivisions could be introduced on their basis.

Conclusions
In reviewing the reliability methods for evaluating earth slope stability, several important conclusions emerged. First of all, this issue has probably been the most published (as a geotechnical problem in terms of reliability) both because of its importance and level of difficulty. Many methods of probabilistic analysis have been developed, and thanks to the rapid progress of computer techniques, there has also been rapid development and acceleration of reliability analysis methods. It was also important to improve the method of ground investigation and to introduce the spatial variation of soil parameters into the calculation model. Unfortunately, the papers on reliability analysis of slopes actually cover basically all possible aspects, while the journals or proceedings in which they are published are not thematically assigned. Thus an attempt was made to organise the huge number of publications by assigning them to several research directions. In this paper, a review of the literature on reliability assessment of earth slope stability is presented in a fairly concise manner.
Generally, the engineer is quite conservative when it comes to applying new methods, especially when they are not simple-and such are probabilistic methods and reliability analysis. The terms and concepts explained in this paper as well as the general characteristics of the methods used should help the engineer to understand the benefits of the reliability approach to geotechnical problems and contribute to its application in practice.
Many probabilistic methods of slope stability assessment are known, and they can be divided according to different criteria. However, the classification of these methods for earth slopes has not yet been developed. The disorder that occurs as a result may make it difficult or discouraging for an engineer to read publications. The classification system proposed in this article organises these methods and will contribute to a better perception of this difficult issue.
This literature review shows that although reliability methods are extremely advanced and the computational possibilities almost unlimited, in practice, the simplest methods are usually used. This is mainly due to the fact that engineers are not familiar with probabilistic concepts; thus, it is difficult to incorporate them into practice.
The author is aware of the fact that the literature review included in the paper is incomplete. Certainly, a number of publications unavailable to him were omitted here. First of all, in accordance with the title of the work, the study area was narrowed thematically only to earth slopes. Environmental impacts, e.g., those caused by earthquakes, excessive precipitation or temperature changes, were also ignored. However, the reliability methods of slope stability analysis used in these cases are included in the proposed computing classification system. It is also possible that among the great number of recent proposals, there are methods that do not fit into the presented classification system, and it could be necessary to modify it accordingly.