Application of a Probabilistic Method Based on Neutrosophic Number in Rock Slope Stability Assessment

Featured Application: A vector similarity measure under neutrosophic number environment was developed for rock slope stability assessments with multiple levels of attributes. This method does not require sophisticated modeling of slopes and the data demand is relatively low. It could serve as a quick preliminary method to identify the key inﬂuential parameters and the risk of natural hazards, such as rock slope failures. Abstract: The stability of natural rock slopes is influenced by a wide spectrum of factors, such as mechanical properties of bedrocks and spatial distribution of discontinuities. Their specific values are typically incomplete, due mainly to the lack of effective and comprehensive methods to accurately characterize these factors, especially those inside of the slopes. The neutrosophic number is a useful tool to solve problems in indeterminate environment. This study introduces the neutrosophic theory into slope stability assessment. A vector similarity measure developed under neutrosophic environment was employed to establish a stability assessment method considering multilevel attributes of slopes. Using this method, the level of stability for studied slopes, i.e., stable, mostly stable, less stable, and instable, was determined by computing the relation indices. The method was applied to a group of rock slopes located in Zhejiang province, China, and the calculated results were compared with the reality of in situ survey. The field application showed that the developed method has a good efficiency and precision in assessing the stability of rock slopes. The obtained weight vector can reveal the key influential parameters that inherently control the stability of rock slopes.


Introduction
Rock slope failures are a typical type of geological hazard that feature large scale and serious damage, frequently occurring in mountainous regions. Rock slopes can be considered as a complex system that consists in random, discrete, and nonlinear information, such as the characteristics of geological structures, geomorphic types, and environmental impacts [1][2][3]. For the better protection of civil engineering, environmental conservation, and efficient operation, the assessment of rock slope stability has always been a challenging topic in the world.
In the past several decades, different approaches have been developed to assess the stability of rock slopes, which can be primarily categorized into deterministic and probabilistic methods. Limit equilibrium and numerical analysis are the most common forms of deterministic method. The limit equilibrium method utilizes the static equilibrium principle to analyze the stress state under various failure modes of slopes [4], including Bishop's, Janbu's, Morgenstern-Price's and Fellenius' techniques, etc. [5][6][7]. Meanwhile, different numerical methods, such as finite element method, discrete element method, discontinuous deformation analysis, and numerical manifold method, have been employed in slope stability assessments and landslide evolution, which have achieved great success [8][9][10][11][12]. In engineering practices, the amount of available data for analysis is commonly limited. It has been realized that the deterministic approaches frequently fail to give reliable assessment when complex slope systems are encountered [13,14], owing to the inadequacy of information regarding site characteristics as well as the inherent variability and measurement errors in geological and geotechnical parameters [15].
With the need to consider the risk and consequence of slope instability and parameter uncertainty, probabilistic approaches have been developed and applied to slope stability assessment. Different mathematics methods have been employed for rock slope engineering, such as the reliability theory, the fuzzy mathematics, the grey theory, and the artificial neural network (ANN). Early studies developed assessment models directly based on a single probabilistic method [16][17][18][19][20]. It has been realized later that the combination of different probabilistic methods can yield better assessments, such as the hybrid model of the projection pursuit (PP), particle swarm optimization (PSO), and the logistic curve function (LCF) [21]. The probabilistic methods have also been introduced into different failure criterions and verified against physical models to overcome the limitation of conventional methods without considering randomness and uncertainty of the slope status and mechanical parameters [22][23][24][25][26][27][28]. With the effort for probabilistic sensitivity analysis, more sophisticated models led to improved accuracy of assessment [29]; however, huge challenges remain in the high computational effort, complex modeling procedures, and the large amount of data required to yield results with acceptable accuracy.
In recent years, the neutrosophic theory, as a new theory of uncertainty, has been accepted and rapidly applied to different fields, such as electrical engineering, medical diagnosis, and roughness estimation [30][31][32]. Neutrosophic number (NN) is a branch of neutrosophy that evolved from classic philosophy. Smarandache [33][34][35] originally proposed the concept of neutrosophic number, which consists of both indeterminate and determinate information for a problem with uncertainty. Ye [36,37] developed three vector similarity measures, i.e., Jaccard, Dice, and cosine similarity measures, and incorporated them into the neutrosophic number to solve problems related to decision-making. Ye [38] later combined possibility degree ranking method and ordered weighted aggregation operators of interval neutrosophic numbers to handle decision-making problems. More recently, several extensions of neutrosophic theory have been put forward to deal with different types of problems. Roy and Das [39] solved multicriteria production planning problem by linear programming approaches. Abdel-Baseet et al. [40] proposed a new approach of big data analysis by neutrosophic mining algorithm, which increased the number of generated association rules. The neutrosophic number provides a simple and reliable way for assessment of uncertainty in the real world, such as probability analysis of various failure phenomena.
Extensive and detailed geological survey on slopes is essential to establish an effective regional plan for risk assessment and management of slope disasters, which, however, is not always practical in many mountainous regions, where only the geological data of a few selected slopes typically with greater risks are available. Given that most probabilistic methods require a sufficient amount of data to construct a baseline or a rule for stability assessment, the insufficiency of geological data has been a significant obstacle faced by researchers in the applications of probabilistic methods. In the present study, a new probabilistic methodology of rock slope stability assessment is developed based on the theory of neutrosophic number, aiming to provide a fast and reliable approach to assess the rock slope stability with limited geological data.

General Concept of Neutrosophic Number
Some scientific problems with uncertainty cannot be directly estimated in a quantitative manner in the real world. The uncertain parts need to be translated into mathematical forms for utilization. Hence, Smarandache [33] originally proposed a concept of NN to express real information. In reality, different kinds of data are used to denote different attributes, which are distinct and can be expressed by crisp numbers, such as 0 and 1. Under neutrosophic environment, the crisp number can be denoted as N = Q + EI, and this expression contains two parts, i.e., the determinate part Q and the indeterminate part EI, where the values Q and E represent the real number. The indeterminacy I is a dynamic interval number owing to the domain, which is variable according to the actual situation. The indeterminacy I has the following properties: I 2 = I, 0 × I= 0, and I/I = undefined.
Taking a neutrosophic number N = 6 + 4I, as an example, if I ∈ [0, 0.5], it is equivalent to N ∈ [6,8] and therefore N ≥ 6. This number consists of the determinate part 6 and the indeterminate part 4I. The domain I ∈ I L , I U can be selected according to the practical situation. Here, L and U stand for the lower and upper bounds for the domain, respectively. When an event is considered to be true, the value of indeterminacy becomes 0. For example, if the value of a true event equals to 8, the neutrosophic number will be N ∈ [8,8], and it is defined as N = 8 + 0 × I.

Development of a Multiattribute Similarity Measure
The Dice similarity measure is a kind of similarity evaluation method that is used to calculate the similarity degree between two objectives [41]. In regard to statistics, this similarity degree indicates the amount of association between two given objectives compared to the amount of association between them expected by chance [42]. Compared with other similarity measures, it can evaluate the similarity when one vector is equal to zero and its calculation process is relatively convenient. Here, the Dice similarity measure was incorporated into the neutrosophic number to develop a method for slope stability assessment, since slope is a complex system that consists of a series of variables with different ranges.
Considering two vectors A = (a 1 , a 2 , · · · , a n ) and B = (b 1 , b 2 , · · · , b n ) that feature length n and positive coordinates. The Dice similarity measure can be defined as Appl. Sci. 2019, 9, 2309 4 of 18 Presume that there are two neutrosophic numbers N Aj = a Aj + b Aj I and N Bj = a Bj + b Bj I. For convenience, those numbers can be written as two sets A = {N A1 , N A2 , · · · , N An } and B = {N B1 , N B2 , · · · , N Bn }. Then, the neutrosophic Dice similarity measure is defined as D(A, B) = 1 n n j=1 2[(a A j +inf(b Aj I))×(a Bj +inf(b Bj I))+(a Aj +sup(b Aj I))×(a Bj +sup(b Bj I))] (a Aj +inf(b Aj I)) 2 +(a Bj +inf(b Bj I)) 2 + (a Aj +sup(b Aj I)) 2 +(a Bj +sup(b B j I)) 2 (8) Obviously, a single value can be denoted as a special interval between the lower and upper bounds. For example, a single value 6 can be transformed into an interval value [6,6] that is computable in Equation (8). This measure D(A, B) has the following properties.
The measure with a single attribute does not satisfy the requirement of slope stability assessment because a slope possesses attributes of various types and levels. Multiple attribute method has to be put forward to solve this problem.
The neutrosophic numbers N Aj and N Bj were further divided into N Ajk and N Bjk , which restructured the evaluation system to two levels. By doing so, both attribute and subattributes of slopes can be considered in the developed method. Let N An = N Aj1 , N Aj2 , · · · , N Ajr and N Bn = N Bj1 , N Bj2 , · · · , N Bjr be two sets of subattributes that can be expressed as N Ajk = a Ajk + b Ajk I and N Bjk = a Bjk + b Bjk I. Then, the weight of elements N Aj and N Bj are given as w j ∈ [0, 1] with n j=1 w j = 1. The weight of subelements is evaluated as r k=1 w k = 1. and w k ∈ [0, 1]. Hence, the neutrosophic Dice similarity measure of multiple elements is denoted as D W (A, B)= n j=1 w j r k=1 w k 2[(a Ajk +in f (b A jk I))×(a B jk +inf(b Bjk I))+(a A jk +sup(b Ajk I))×(a B jk +sup(b Bjk I))] (a Ajk +in f (b Ajk I)) 2 +(a Bjk +in f (b Bjk I)) 2 + (a A jk +sup(b Ajk I)) 2 +(a Bjk +sup(b B jk I)) 2 (12) D W (A,B) also satisfies the properties denoted in Equations (9)-(11).

Classification of Slope Properties
Before applying the developed method, it is required to classify the data of slopes to fit the algorithm. First, slope stability was classified into four grades: Stable (I), mostly stable (II), less stable (III), and instable (IV), which is commonly employed in previous studies and is recommended by related technical codes [43]. Then, according to the geological survey results, as presented in the next section, a set of predominant influencing factors was selected for stability classification analysis. The multiattribute classification method was adopted, which classified the property of slopes into three attributes and each attribute consisted of three subattributes. The slope classification is described as follows. In Table 1, the subattributes have corresponding units, such as degree and meter, because these values were directly extracted from the in situ survey data. These real values are not suitable for a mathematical computation. They were thereafter transformed to normalized values by dividing each value by the maximum in each index. The normalized values are tabulated in Table 2. Note that each subattribute has slightly different intervals because the original data have different ranges, units, and characters. Variable intervals were therefore defined to accommodate their characters [44]. In the next step, the range for each subattribute was transformed into neutrosophic number to indicate its incomplete and indeterminate information, as tabulated in Table 3. The indeterminacy is set to I ∈ [0, 0.24]. Here, the range of I could be assigned by any value between 0 and 1 without changing the assessment results.
When applying the neutrosophic Dice similarity measure to the collected data, the weight for each attribute and subattribute should be assigned as shown in Equation (12). The appropriate value of weight depends on the characteristics of all collected data and may be determined in a trial-and-error manner. To investigate the influence of weight on the assessment performance, all collected slopes were randomly divided into two groups (A and B), ten for each one. Group A was used to determine the weight vector and group B was used to validate the performance of the developed method. We selected three samples, five samples, six samples, eight samples, and ten samples from group A to calculate a set of weight vectors, respectively. These weight vectors were employed to assess the stability of slopes in group B. Because slopes in group B did not participate in the determination of weight, such a procedure could effectively appraise whether reasonable weight vectors exist in the developed method. The minimum sample number required to calculate a reasonable set of weight vectors can be determined.

Classification of Slope Mass Rating
To improve the understanding on the inherent correlation between the geological features of the slope with the stability, the following slope mass rating (SMR) is used to quantitatively estimate the stability of the slope [45][46][47].
where CSMR is the classification of slope mass rating, RMR (rock mass rating) is an empirical classification system for rock masses developed by Bieniawski [48]; F 1 -F 4 are the coefficients of correction for the dip direction and dip angle of joint sets for different rock mass type λ, the attitude correlation between slope and the controlling joint sets, and the excavation method, respectively; ξ is the altitude correction factor according to the height of slope H. After CSMR estimation, the defined classes can be linked to the instability level of the slope as tabulated by Romana [45], where a total of 100 points were divided into equally spaced four levels (I-IV), indicating different degrees of stability.

Geological Setting
The study area is located in Zhejiang province, the southeastern coastal region of China with a continental area of 101,800 km 2 . As shown in Figure 1, the hilly and mountainous terrain with elevation above 300 m accounts for 70.6% of its total area. The altitude decreases gradually from the southwest to the northeast and drops in a step form. The southwest part is a mountainous region with the existence of several peaks around 1000 m high. The central part is dominated by hills and scattered basins. The northeast portion is dominated by an alluvial plain. Appl. Sci. 2019, 9, x 7 of 18 This region suffers from complex geological and extreme climatic settings. Figure 2 shows the exposed bedrock in the mountainous area that consists mainly of intensely weathered granite and tuff [49]. With a geochemical process, the bedrock received strong weathering in this region. The degree of fragmentation alters from intact rock to gravel-like based on the degree of weathering. Generally, these weathering profiles have more complex engineering properties than other sedimentary soils [50]. Located in the coastal region, rainstorms are frequently encountered from May to October. The annual precipitation is around 1380 mm, and the maximum monthly precipitation of over 200 mm is typically recorded in June [51]. The weak geological properties and heavy rainfall brought by typhoon lead to geological disasters of great frequency and scale, which cost significant damages to lives and properties [52]. According to the geological investigation report in the research area from Department of Land and Resources of Zhejiang Province, 2832 geological disasters occurred in the recent five years, resulting in more than 50 civilian casualties as well as a direct economic loss of around US$40 million [53]. Prevention of geohazards is therefore of great importance in Zhejiang province. Similar disaster events have been reported worldwide, showing a strong demand for reliable field survey and instability assessment methods for disaster prevention [54,55].
To study the applicability of the developed method, a representative area in the central north of Zhejiang province was selected, which was located on a fault and was frequently struck by geohazards. A dataset of 20 rock slopes from 12 locations was collected through geologic survey as shown in Figure 3, where the data meet the requirements for the developed method. The study area is divided into two zones by the NE-striking Jiangshan-Shaoxing fault. The northwest area belongs to the Yangtze paraplatform, where the strata are well developed from middle Proterozoic group to Quaternary of the Cenozoic. The magmatic activity is weak, and sedimentary deposits are well exposed. In contrast, the southeast portion belongs to the Southern China fold system, which undergoes strong magmatism, leading to obvious losses of formation. It is covered with igneous rocks and the outcrops are mainly formed by the Mesozoic and Cenozoic continental sedimentary This region suffers from complex geological and extreme climatic settings. Figure 2 shows the exposed bedrock in the mountainous area that consists mainly of intensely weathered granite and tuff [49]. With a geochemical process, the bedrock received strong weathering in this region. The degree of fragmentation alters from intact rock to gravel-like based on the degree of weathering. Generally, these weathering profiles have more complex engineering properties than other sedimentary soils [50]. Located in the coastal region, rainstorms are frequently encountered from May to October. The annual precipitation is around 1380 mm, and the maximum monthly precipitation of over 200 mm is typically recorded in June [51]. The weak geological properties and heavy rainfall brought by typhoon lead to geological disasters of great frequency and scale, which cost significant damages to lives and properties [52]. According to the geological investigation report in the research area from Department of Land and Resources of Zhejiang Province, 2832 geological disasters occurred in the recent five years, resulting in more than 50 civilian casualties as well as a direct economic loss of around US$40 million [53]. Prevention of geohazards is therefore of great importance in Zhejiang province. Similar disaster events have been reported worldwide, showing a strong demand for reliable field survey and instability assessment methods for disaster prevention [54,55].
To study the applicability of the developed method, a representative area in the central north of Zhejiang province was selected, which was located on a fault and was frequently struck by geohazards. A dataset of 20 rock slopes from 12 locations was collected through geologic survey as shown in Figure 3, where the data meet the requirements for the developed method. The study area is divided into two zones by the NE-striking Jiangshan-Shaoxing fault. The northwest area belongs to the Yangtze paraplatform, where the strata are well developed from middle Proterozoic group to Quaternary of the Cenozoic. The magmatic activity is weak, and sedimentary deposits are well exposed. In contrast, the southeast portion belongs to the Southern China fold system, which undergoes strong magmatism, leading to obvious losses of formation. It is covered with igneous rocks and the outcrops are mainly formed by the Mesozoic and Cenozoic continental sedimentary rocks. There are several distinct geological outcrops in this region, which belong to various geological periods, such as Proterozoic, Paleozoic, Mesozoic, and Cenozoic. rocks. There are several distinct geological outcrops in this region, which belong to various geological periods, such as Proterozoic, Paleozoic, Mesozoic, and Cenozoic.  rocks. There are several distinct geological outcrops in this region, which belong to various geological periods, such as Proterozoic, Paleozoic, Mesozoic, and Cenozoic.

Characteristics of Studied Rock Slopes
A detailed engineering geological survey was conducted on the studied slopes to obtain the main status and characteristics. According to the suggested method by ISRM (2007-2014), point loading test (PL) and Schmidt hammer rebound test (SH) were employed to estimate the characteristics of weathering profiles [56]. Fractured and loose structures dominate the shallow rock masses in the range of 10-20 m combining means of geological logging and drilling data. The dimension of most slopes is small to medium (less than 150 m) and the amount of debris of failed slopes is therefore relatively small. It is indicated from the statistics of geological disasters in the study area that slope failures are bound up with rainfall [49,57,58]. Seismic activities are extremely low, while over half slopes suffered from human activities such as quarrying and road constructions in different degrees. Five of the total 20 slopes have failed in the past, and one example of stable slopes and one example of collapsed slopes are described below to highlight their geological characteristics. Given the close distance, other slopes in the studied area exhibit similar characteristics and the details are quantified by the nine subattributes (Table 1) shown below.

Jidong Town Quarry Slope
Jidong Town rock slope is an abandoned quarry located in the southwest piedmont of Mountain Dabu. It is an excavated slope located in the leading edge of the mountain. The maximum height is 65.8 m, and the slope dip orientation is 265 • N. The upper part and the lower part exhibit different geological characteristics as shown in Figure 4. The upper part is covered by debris deposit, including sandy soil mixed with gravels. The measured point loading strength is less than 2 MPa and the rebound number is less than 15, which together with the visual investigation indicate that these rocks are completely weathered. This part has an average height of 18.6 m and an average slope angle of 45 • . The lower terrain is steeper (about 85 • ) than that of the upper part, and the rocks are slightly weathered according to the test results. At the toe of the slope, there is an artificial excavation platform with a width of 54 m. Construction facilities and quarrying induced deposits can be observed on Figure 4.

Characteristics of Studied Rock Slopes
A detailed engineering geological survey was conducted on the studied slopes to obtain the main status and characteristics. According to the suggested method by ISRM (2007-2014), point loading test (PL) and Schmidt hammer rebound test (SH) were employed to estimate the characteristics of weathering profiles [56]. Fractured and loose structures dominate the shallow rock masses in the range of 10-20 m combining means of geological logging and drilling data. The dimension of most slopes is small to medium (less than 150 m) and the amount of debris of failed slopes is therefore relatively small. It is indicated from the statistics of geological disasters in the study area that slope failures are bound up with rainfall [49,57,58]. Seismic activities are extremely low, while over half slopes suffered from human activities such as quarrying and road constructions in different degrees. Five of the total 20 slopes have failed in the past, and one example of stable slopes and one example of collapsed slopes are described below to highlight their geological characteristics. Given the close distance, other slopes in the studied area exhibit similar characteristics and the details are quantified by the nine subattributes (Table 1) shown below.

Jidong Town Quarry Slope
Jidong Town rock slope is an abandoned quarry located in the southwest piedmont of Mountain Dabu. It is an excavated slope located in the leading edge of the mountain. The maximum height is 65.8 m, and the slope dip orientation is 265°N. The upper part and the lower part exhibit different geological characteristics as shown in Figure 4. The upper part is covered by debris deposit, including sandy soil mixed with gravels. The measured point loading strength is less than 2 MPa and the rebound number is less than 15, which together with the visual investigation indicate that these rocks are completely weathered. This part has an average height of 18.6 m and an average slope angle of 45°. The lower terrain is steeper (about 85°) than that of the upper part, and the rocks are slightly weathered according to the test results. At the toe of the slope, there is an artificial excavation platform with a width of 54 m. Construction facilities and quarrying induced deposits can be observed on Figure 4.  Geological logging by fine line method was adopted in the field survey with a total length of 220 m. The results indicated that the slope mainly comprises three joint sets. The orientations of these joints, given as dip direction and dip angle in degrees, are N205∠30-50 • (J 1 ), N340∠65-75 • (J 2 ), and N245∠30-40 • (J 3 ). The spacing of joint sets J 1 and J 2 ranges between 0.5 m and 1 m, while the spacing of joint set J 3 ranges from 1 m to 2 m. The orientation of inclination of the major joint sets is close to that of the slope surface, while their inclination angles are smaller than the slope angle. The tested PL index is over 400 MPa and the SH rebound is over 50 for the bedrock, which indicate that the bedrock has a good quality. Although three sets of joints exist in the slope, they do not form wedges or other structures that are unfavorable to stability. Given these geological features, the Jidong town quarry slope is presently in a stable state, which is also confirmed by the judgment of experts based on in situ survey.

Daoxu Town Quarry Slope
Daoxu Town rock slope is a quarry situated towards the north piedmont of Mountain Chen. The major lithological unit is the lower Cretaceous Gaowu formation. It comprises low-grade ignimbrite in dark grey color and has a loosen structure with lots of cranny. The slope orientation is about N52 • ∠70 • and the average height is 65.8 m. The estimated slope area is about 9000 m 2 . As shown in Figure 5, a sudden gradient change appears in the middle of the slope, above which a steep bedrock cliff exists and under which a gentle accumulative zone exists. According to the results of PL and SH tests, the rocks on the top and surfaces are completely decomposed, and the strength is weak. m. The results indicated that the slope mainly comprises three joint sets. The orientations of these joints, given as dip direction and dip angle in degrees, are N205∠30-50° (J1), N340∠65-75° (J2), and N245∠30-40° (J3). The spacing of joint sets J1 and J2 ranges between 0.5 m and 1 m, while the spacing of joint set J3 ranges from 1 m to 2 m. The orientation of inclination of the major joint sets is close to that of the slope surface, while their inclination angles are smaller than the slope angle. The tested PL index is over 400 MPa and the SH rebound is over 50 for the bedrock, which indicate that the bedrock has a good quality. Although three sets of joints exist in the slope, they do not form wedges or other structures that are unfavorable to stability. Given these geological features, the Jidong town quarry slope is presently in a stable state, which is also confirmed by the judgment of experts based on in situ survey.

Daoxu Town Quarry Slope
Daoxu Town rock slope is a quarry situated towards the north piedmont of Mountain Chen. The major lithological unit is the lower Cretaceous Gaowu formation. It comprises low-grade ignimbrite in dark grey color and has a loosen structure with lots of cranny. The slope orientation is about N52°∠ 70° and the average height is 65.8 m. The estimated slope area is about 9000 m 2 . As shown in Figure  5, a sudden gradient change appears in the middle of the slope, above which a steep bedrock cliff exists and under which a gentle accumulative zone exists. According to the results of PL and SH tests, the rocks on the top and surfaces are completely decomposed, and the strength is weak.
Four major joint sets were identified and the orientations of them are N90°∠18° (J1), N75°∠75° (J2), N270°∠80° (J3), and N130°∠65° (J4). Joint sets J1 and J2 run parallel to the bedding plane with the spacing varying from 0.5 to 1 m. Because of the large inclination of both joint sets J3 and J4, and their spatial correlation with the slope face, wedge-like key blocks hanging on the slope were formed through the combination of these joints ( Figure 5). In addition, a great amount of fractures exists in the slope due to the blasting for quarrying in the past. The hanging rocks of Daoxu Town rock slope collapsed on June 2017 when a monthly precipitation of 208.2 mm and a maximum daily precipitation of 373.4 mm were recorded, and formed the accumulative zone in the toe.  Four major joint sets were identified and the orientations of them are N90 • ∠18 • (J 1 ), N75 • ∠75 • (J 2 ), N270 • ∠80 • (J 3 ), and N130 • ∠65 • (J 4 ). Joint sets J 1 and J 2 run parallel to the bedding plane with the spacing varying from 0.5 to 1 m. Because of the large inclination of both joint sets J 3 and J 4 , and their spatial correlation with the slope face, wedge-like key blocks hanging on the slope were formed through the combination of these joints ( Figure 5). In addition, a great amount of fractures exists in the slope due to the blasting for quarrying in the past. The hanging rocks of Daoxu Town rock slope collapsed on June 2017 when a monthly precipitation of 208.2 mm and a maximum daily precipitation of 373.4 mm were recorded, and formed the accumulative zone in the toe.

Results
Taking Daoxu Town Quarry slope as an example, Tables 4 and 5 tabulate the values for RMR and CSMR estimations, respectively. A value of 30.8 was obtained by Equations (13) and (14) of the slope, and therefore the slope is classified to grade III, which is in the range from 25 to 50 in CSMR. This result revealed the unstable state of the slope, however, given the actual collapsed state, a grade IV should be estimated to fit the reality. It is considered that this discrepancy was yielded due to the incomplete and imprecise information being used in the estimation, which again suggested the necessity of the introduction of probabilistic methods to slope stability assessments. The stability grades for slopes of group A and group B are tabulated in Tables 6 and 7, respectively. These grades were initially determined based on in situ survey data using ordinary estimation methods such as CSMR, followed by cross-examination by experienced experts. After modification and validation, these grades mostly represented the actual status of the studied slopes, which are labelled as T 1 -T 10 . Here, Jidong Town quarry slope corresponds to T 4 , and Daoxu Town quarry slope corresponds to T 9 in Table 6. There are five slopes that are ranked as IV, which mean that collapse or failure had happened on them.  III  IV  III  III  I  II  IV  IV  I  II Three rock slope samples (T 1 , T 2 , and T 3 ) from group A were firstly selected to determine the weight w for the attributes, and w 1 , w 2 , and w 3 for the subattributes. As shown in Table 1, there are three indices for attribute or subattribute, therefore, the weight vector has a form of W = (x 1 , x 2 , x 3 ), where x 1 + x 2 + x 3 = 1. A simple algorithm based on MATLAB was developed, which tested different combinations of weight values to identify the suitable vector under that the calculated stability grades of T 1 , T 2 , and T 3 agree with their actual values (III, IV and II) as tabulated in Table 6.
The calculated weight vectors are as follows. This set of weight vectors was employed to estimate the stability of slopes of group B via Equation (12). Note that in calculation, all the single values tabulated in Table 6 should be expressed as neutrosophic number. For example, a single value 0.3 is transformed into [0.3, 0.3], and its neutrosophic number is The calculated possibilities D W of each grade for the slopes in group B are tabulated in Table 8. The maximum value of D W for each slope represents its calculated stability grade. Comparison between the calculated and the actual grades shows that the stability grades of slopes T 2 , T 7 , and T 8 are underestimated using the current set of weight vector. The accuracy is 70% when the weight vectors are determined based on the data of three slopes.  Tables 9 and 10. In this case, the accuracy increased to 80%, but the stability grades of T 2 and T 7 were still mismatched.

Discussions
Generally, a greater weight value indicates that an attribute or a subattribute is more intimately correlated with the stability of a slope. The determined weight vector for the three attributes (0.6, 0.2, 0.2) shows that the rock mass characteristics play the most important role in the stability of slopes, which agrees with the common understanding that the physical, mechanical, and geological properties of rock masses govern the stability of rock slopes [4,8,14,17,25,46,47]. Meanwhile, among the three subattributes of the rock mass characteristics, rock structure has greater importance over other two indexes, which is in line with the actual situation that most studied instable slopes are formed by bedding structures with multiple joint sets [59]. Moreover, the plant cover is particularly important in the topographic features, which agrees with the fact that most failed slopes in the studied area are barely covered by plants mainly due to human activities such as quarrying, and road and house constructions as demonstrated by the Daoxu Town quarry slope. Although the precipitation and seismic intensity are almost identical for all the studied slopes given the limited area selected, the precipitation seems to have greater importance than the seismic intensity. In situ survey showed that the majority of failure incidents of studied slopes happened after heavy rainfalls, such as the Daoxu Town quarry slope, demonstrating the important role precipitation plays in the stability of slopes. The maximum daily precipitation may be a better parameter for characterizing the influence of rainfall as long as the data are available. It is noted that the seismic intensity in Zhejiang Province is remarkably low, and no strong earthquakes have been recorded in recent 50 years.
These results show that the rock slopes formed by bedding structures of fractured rocks disturbed by human activities and undergone heavy rainfall are most unstable, which are in concert with our common knowledge on slope stability [60][61][62]. The developed method can help quantitatively reveal the relative importance of different influencing factors, and help pinpoint the key factors that trigger instability of slopes. This could contribute to the policy-making regarding maintenance works of slopes.
The above results and analysis indicate that the developed method holds a straightforward computational manner. Compared with other probabilistic methods such as ANN that typically require a complex learning process for updating stability knowledge [63,64], the present method only requires a simple trial-and-error approach to determine the weight vector. Meanwhile, ANN requires a large amount of data that possess the main characteristics of a problem for learning and training, which need to be repeated when new data are included in the assessment. In contrast, the requirement of data amount is significantly small by the present method, making it suitable for application to regional slope stability assessment with limited geological survey data. Once an effective weight vector is determined, the method can be applied to the slopes of a whole region possessing similar geological background. Given the simplicity in computation and low demand on input data, this method, supported by expert opinion, can serve as a fast and reliable tool for regional stability assessment of rock slopes.

Conclusions
The present study developed a vector similarity measure under neutrosophic number environment for rock slope stability assessments with multiple levels of attributes. The advantage of the method lies on that it significantly simplifies the computational process of fuzzy problems with indeterminate information while maintaining the capability of similarity analysis on multiattribute information based on limited input data.
When applying this method, the weight vector has to be firstly determined based on the data of at least eight slope samples. This amount of data supply is much lower than other representative probabilistic methods. The determined weight vector reflects the relative importance of different parameters on the stability of slopes in a studied region. It could help understand the key influential parameters for the stability of a slope and take corresponding measures to improve the stability. By incorporating this vector into the multiattribute classification of the influential parameters of slope stability, the developed method can accurately estimate the stability degree of studied slopes. This method does not require sophisticated modeling of slopes and the data demand is relatively low. It could serve as a quick preliminary method to identify the key influential parameters and the risk of rock slope failures. Thorough in situ survey and real-time monitoring could then be applied to the identified slopes at risk to save cost and time.
In the present study, the weight vector was determined via a trial-and-error manner using different amounts of data. In the future, this procedure needs to be improved by developing an algorithm for weight determination. This method will be applied to multiple regions with insufficient or incomplete geological survey data to verify its applicability. The attributes and subattributes need to be continuously improved to establish a more efficient classification system that fits the neutrosophic theory and could well represent the geological characteristics of studied regions.