Probabilistic Slope Stability Analysis of Mount St. Helens Using Scoops3D and a Hybrid Intelligence Paradigm

: In the past, numerous stratovolcanoes worldwide witnessed catastrophic ﬂank collapses. One of the greatest risks associated with stratovolcanoes is a massive rock failure. On 18 May 1980, we witnessed a rock slope failure due to a volcano eruption, and a 2185.60 m high rock slope of Mount St. Helens was collapsed. Thus, from the serviceability perspective, this work presents an effective computational technique to perform probabilistic analyses of Mount St. Helens situated in Washington, USA. Using the ﬁrst-order second-moment method, probability theory and statistics were employed to map the uncertainties in rock parameters. Initially, Scoops3D was used to perform slope stability analysis followed by probabilistic evaluation using a hybrid computational model of artiﬁcial neural network (ANN) and ﬁreﬂy algorithm (FF), i.e., ANN-FF. The performance of the ANN-FF model was examined and compared with that of conventional ANN and other hybrid ANNs built using seven additional meta-heuristic algorithms. In the validation stage, the proposed ANN-FF model was the best-ﬁtted hybrid model with R 2 = 0.9996 and RMSE = 0.0042. Under seismic and non-seismic situations, the reliability index and the probability of failure were estimated. The suggested method allows for an effective assessment of the failure probability of Mount St. Helens under various earthquake circumstances. The developed MATLAB model is also attached as a supplementary material for future studies.


Introduction
Volcanic edifices have slope failures that range from small rocks falling to massive collapses. Enormous flank (with volume > 0.1 km 3 ) slides have significantly changed over 200 stratovolcanoes across the globe [1] and are among the most sudden, damaging, and potentially hazardous volcanic events. As a result of catastrophic side collapses, many volcanoes and stratovolcanoes have changed in big ways worldwide [2]. The most significant volcanic collapses in history, such as those observed at Bandai (Japan), Mount St. Helens (USA), and Bezymianny (Russia), were initiated by the movement along extensive curved failure surfaces [3]. In the United States, Mount St. Helens had a terrible fall in 1980 [4]. Debris avalanches, which can later mobilize into debris disasters, pose major threats to the structure and areas further downslope or downstream when these enormous collapses occur, often involving more than 0.1 km 3 of material. Moreover, a significant number of approximately 700 stratovolcanoes found on Earth present a potential hazard to individuals

Research Significance
Numerous studies have emphasized the use of ANNs in the domains of engineering and science [28]. The ability of ANNs to represent non-linear problems without considering a functional relationship between input and output is a significant advantage over other ML algorithms. Furthermore, the output generation is unaffected by one or more corrupted cells. Despite these advantages, ANN possesses notable limitations, including challenges related to the entrapment in local minima and the occurrence of overfitting. Moreover, the challenge of accurately determining the precise global minimum can lead to unfavorable outcomes [28]. To address these issues, researchers employed various optimization algorithms (OAs), including PSO, ABC, and genetic algorithm (GA) [28]. Due to their robust global search capabilities, OAs can iteratively optimize the learning parameters of ANNs, resulting in enhanced prediction performance. Over the past decade, there has been a significant increase in the utilization of hybrid ANNs, such as ANN-ABC, ANN-GA, ANN-PSO, etc., for addressing various problems [28,37]. Tun et al. [38] used GA to assess 3D slopes with several failure regions. Nevertheless, it should be noted that there is a lack Mathematics 2023, 11, 3809 4 of 23 of research on the reliability of slope analysis for the Mount St. Helens utilizing hybrid ANNs. Moreover, no study has been conducted on the evaluation of hybrid ANNs that are built using different groups of OAs for slope stability analysis in cone-shaped terrains.
Thus, this study aims to address the gap in the literature by using hybrid ANNs built with different groups of OAs to perform a probabilistic analysis of Mount St. Helens. Specifically, eight distinct OAs viz., ALO, BBO, CPA, DE, EO, FF, GA, and PSO (see Section 5.2 for details) were used to optimize weights and biases ANNs, resulting in eight hybrid ANNs, viz., ANN-ALO, ANN-BBO, ANN-CPA, ANN-DE, ANN-EO, ANN-FF, ANN-GA, and ANN-PSO. The computational findings were used to select the most effective hybrid ANN model for performing RA of Mount St. Helens under seismic and non-seismic conditions. The outcomes of the current study were compared to the findings of Tun et al. [38] and evaluated in the subsequent sections.

Study Area
This study investigates RA of Mount St. Helens, located in Skamania County, Washington, USA, approximately 55 km west of Mount Adams in the western part of the Cascade Range, experienced a catastrophic eruption on 18 May 1980 [4]. Notably, the failure surfaces may be readily apparent for basic geometrical aspects of a slope (such as polygonal or polyhedral-shaped slopes), but this may not be the case for more complex/real topographies like Mount St. Helens. Thus, the present study utilizes the details of Mount St. Helens (obtained from Google Earth) to conduct probabilistic slope analyses. The use of Google Maps and Google Earth images has facilitated the depiction of Mount St. Helen's geographical location and topographical features, presented in Figure 1.
Before the tragic collapse of Mount St. Helens, rock strength and density distributions in 3D were unknown. After the collapse, Voight et al. [4] and Glicken [39] used geologic sections provided by Hopson and Melson [40] to recreate the inner geology structure. The mountain's core comprises earlier dacite dome lavas and flank breccias hydrothermally altered. This unit was covered by andesitic, basaltic, and tephra lava flows. Summit and Goat Rocks domes were built with more modern dacite. The dacitic crypto-dome intruded on each of these components in the year 1980.
Mathematics 2023, 11, x FOR PEER REVIEW 4 of 24 been a significant increase in the utilization of hybrid ANNs, such as ANN-ABC, ANN-GA, ANN-PSO, etc., for addressing various problems [28,37]. Tun et al. [38] used GA to assess 3D slopes with several failure regions. Nevertheless, it should be noted that there is a lack of research on the reliability of slope analysis for the Mount St. Helens utilizing hybrid ANNs. Moreover, no study has been conducted on the evaluation of hybrid ANNs that are built using different groups of OAs for slope stability analysis in cone-shaped terrains. Thus, this study aims to address the gap in the literature by using hybrid ANNs built with different groups of OAs to perform a probabilistic analysis of Mount St. Helens. Specifically, eight distinct OAs viz., ALO, BBO, CPA, DE, EO, FF, GA, and PSO (see Section 5.2 for details) were used to optimize weights and biases ANNs, resulting in eight hybrid ANNs, viz., ANN-ALO, ANN-BBO, ANN-CPA, ANN-DE, ANN-EO, ANN-FF, ANN-GA, and ANN-PSO. The computational findings were used to select the most effective hybrid ANN model for performing RA of Mount St. Helens under seismic and non-seismic conditions. The outcomes of the current study were compared to the findings of Tun et al. [38] and evaluated in the subsequent sections.

Study Area
This study investigates RA of Mount St. Helens, located in Skamania County, Washington, USA, approximately 55 km west of Mount Adams in the western part of the Cascade Range, experienced a catastrophic eruption on 18 May 1980 [4]. Notably, the failure surfaces may be readily apparent for basic geometrical aspects of a slope (such as polygonal or polyhedral-shaped slopes), but this may not be the case for more complex/real topographies like Mount St. Helens. Thus, the present study utilizes the details of Mount St. Helens (obtained from Google Earth) to conduct probabilistic slope analyses. The use of Google Maps and Google Earth images has facilitated the depiction of Mount St. Helen's geographical location and topographical features, presented in Figure 1.
Before the tragic collapse of Mount St. Helens, rock strength and density distributions in 3D were unknown. After the collapse, Voight et al. [4] and Glicken [39] used geologic sections provided by Hopson and Melson [40] to recreate the inner geology structure. The mountain's core comprises earlier dacite dome lavas and flank breccias hydrothermally altered. This unit was covered by andesitic, basaltic, and tephra lava flows. Summit and Goat Rocks domes were built with more modern dacite. The dacitic cryptodome intruded on each of these components in the year 1980.

Methodology
Theoretical detail of the deterministic analysis to calculate the FOS using Bishop's simplified method is presented in this section, followed by the details of probabilistic analysis to conduct reliability analyses.

Deterministic Analysis
The FOS of slope failure in 3D was calculated using Scoops3D, a computer software that employs Ordinary and Bishop simplified methods. As stated above, Bishop's simplified method was employed in this study, wherein spherical failure surfaces were utilized to analyses 3D slopes. The simplified form of Bishop's approach for calculating FOS in the absence of groundwater is given: where c j,k is the effective cohesion; φ j,k is the effective internal friction angle; R j,k is the distance from the j, k column's trial slip region to its axis of rotation, A j,k is the column's trial surface area, W j,k is the column's weight; e j,k is the horizontal driving force moment arm. A free-body diagram of the j, k column is shown in Figure 2. Note that, for a 3D formulation, it is also assumed that the combined normal and shear forces acting along the sides of the columns be zero in both the x and y axes.

Methodology
Theoretical detail of the deterministic analysis to calculate the FOS using Bishop's simplified method is presented in this section, followed by the details of probabilistic analysis to conduct reliability analyses.

Deterministic Analysis
The FOS of slope failure in 3D was calculated using Scoops3D, a computer software that employs Ordinary and Bishop simplified methods. As stated above, Bishop's simplified method was employed in this study, wherein spherical failure surfaces were utilized to analyses 3D slopes. The simplified form of Bishop's approach for calculating FOS in the absence of groundwater is given: ; , = cos , + tan and = sin , (1) where , is the effective cohesion; , is the effective internal friction angle; , is the distance from the , column's trial slip region to its axis of rotation, , is the column's trial surface area, , is the column's weight; e , is the horizontal driving force moment arm. A free-body diagram of the , column is shown in Figure 2. Note that, for a 3D formulation, it is also assumed that the combined normal and shear forces acting along the sides of the columns be zero in both the x and y axes. In Figure 2, W denotes the column's weight; , and , denote inter-column normal forces in x and y axes, respectively; , and , are horizontal shear forces in y-z plane; , and , denote inter-column shear forces occurring in the x-z plane; , and , denote effective normal force and pore water force at the base, respectively; , denotes mobilized shear force on the base; , denotes x-y plane slide angle, and and denote the base inclination of the middle column in the x-z and y-z planes, respectively.

Digital Elevation Modelling
A digital representation of the topographic relief is commonly called a digital elevation model (DEM) [41]. Notably, regular grids are the most commonly encountered type In Figure 2, W denotes the column's weight; E x j,k and E y j,k denote inter-column normal forces in x and y axes, respectively; H x j,k and H y j,k are horizontal shear forces in y-z plane; X x j,k and X y j,k denote inter-column shear forces occurring in the x-z plane; N j,k and U j,k denote effective normal force and pore water force at the base, respectively; S j,k denotes mobilized shear force on the base; α j,k denotes x-y plane slide angle, and α x and α y denote the base inclination of the middle column in the x-z and y-z planes, respectively.

Digital Elevation Modelling
A digital representation of the topographic relief is commonly called a digital elevation model (DEM) [41]. Notably, regular grids are the most commonly encountered type of DEMs, and they are offered in various forms [42]. DEM is widely used in geomorphology, representing various landscape features and offering other benefits, including efficient data storage and processing power. The accuracy of a DEM is influenced by various factors pertaining to the topography, such as landforms, elevations, texture, ruggedness, and vegetation. The accuracy is also influenced by the methods employed for collecting elevation data, the specific process for generating the DEM, the type of grid employed for the DEM, and the resolution of the DEM [43]. Also, the liDAR-DEM (light detection and ranging-DEM) to the DEM was created using 1 × 1 m horizontal resolution and 0.2 m vertical precision from data collected using airborne laser scanning (ALS) [44][45][46][47][48]. In addition, remote sensing allows for massive mountain ranges, deep valleys, and ocean floor magnetic striping [49].
The United States Geological Survey (USGS) developed a free software package, Scoops3D, to analyze the 3D stability of slopes based on DEM [36]. To begin with, precise calculations of the terrain's topography are required. As a result, because DEM can easily employ spatial and picture data from geographic information systems (GIS), it may be used to develop a 3D model for slope concerns [50]. In the context of probabilistic analyses, it is adequate to consider the variations in material properties along with the slope geometry. Also, with the increasing use of satellite images and aerial photography, DEM representations are more accessible than ever. Reid et al. [36] used 100-m resampled DEMs to accurately compute the stability of potential failures with a volume > 0.1 km 3 . However, the DEM provided by Reid et al. [36] was used in this study, generated from photographs taken on 12 May 1980 (see Figure 3). Additionally, Figure 4 depicts the DEM profile of Mount St. Helens and necessary information was extracted from the work of Reid et al. [36] for this study.
Mathematics 2023, 11, x FOR PEER REVIEW 6 of 24 of DEMs, and they are offered in various forms [42]. DEM is widely used in geomorphology, representing various landscape features and offering other benefits, including efficient data storage and processing power. The accuracy of a DEM is influenced by various factors pertaining to the topography, such as landforms, elevations, texture, ruggedness, and vegetation. The accuracy is also influenced by the methods employed for collecting elevation data, the specific process for generating the DEM, the type of grid employed for the DEM, and the resolution of the DEM [43]. Also, the liDAR-DEM (light detection and ranging-DEM) to the DEM was created using 1 × 1 m horizontal resolution and 0.2 m vertical precision from data collected using airborne laser scanning (ALS) [44][45][46][47][48]. In addition, remote sensing allows for massive mountain ranges, deep valleys, and ocean floor magnetic striping [49]. The United States Geological Survey (USGS) developed a free software package, Scoops3D, to analyze the 3D stability of slopes based on DEM [36]. To begin with, precise calculations of the terrain's topography are required. As a result, because DEM can easily employ spatial and picture data from geographic information systems (GIS), it may be used to develop a 3D model for slope concerns [50]. In the context of probabilistic analyses, it is adequate to consider the variations in material properties along with the slope geometry. Also, with the increasing use of satellite images and aerial photography, DEM representations are more accessible than ever. Reid et al. [36] used 100-m resampled DEMs to accurately compute the stability of potential failures with a volume > 0.1 km 3 . However, the DEM provided by Reid et al. [36] was used in this study, generated from photographs taken on 12 May 1980 (see Figure 3). Additionally, Figure 4 depicts the DEM profile of Mount St. Helens and necessary information was extracted from the work of Reid et al. [36] for this study.   of DEMs, and they are offered in various forms [42]. DEM is widely used in geomorphology, representing various landscape features and offering other benefits, including efficient data storage and processing power. The accuracy of a DEM is influenced by various factors pertaining to the topography, such as landforms, elevations, texture, ruggedness, and vegetation. The accuracy is also influenced by the methods employed for collecting elevation data, the specific process for generating the DEM, the type of grid employed for the DEM, and the resolution of the DEM [43]. Also, the liDAR-DEM (light detection and ranging-DEM) to the DEM was created using 1 × 1 m horizontal resolution and 0.2 m vertical precision from data collected using airborne laser scanning (ALS) [44][45][46][47][48]. In addition, remote sensing allows for massive mountain ranges, deep valleys, and ocean floor magnetic striping [49]. The United States Geological Survey (USGS) developed a free software package, Scoops3D, to analyze the 3D stability of slopes based on DEM [36]. To begin with, precise calculations of the terrain's topography are required. As a result, because DEM can easily employ spatial and picture data from geographic information systems (GIS), it may be used to develop a 3D model for slope concerns [50]. In the context of probabilistic analyses, it is adequate to consider the variations in material properties along with the slope geometry. Also, with the increasing use of satellite images and aerial photography, DEM representations are more accessible than ever. Reid et al. [36] used 100-m resampled DEMs to accurately compute the stability of potential failures with a volume > 0.1 km 3 . However, the DEM provided by Reid et al. [36] was used in this study, generated from photographs taken on 12 May 1980 (see Figure 3). Additionally, Figure 4 depicts the DEM profile of Mount St. Helens and necessary information was extracted from the work of Reid et al. [36] for this study.

Probabilistic Analysis
Inaccuracies in modelling, subsurface stratigraphic uncertainty, and intrinsic spatial variability of soil and rock characteristics are some of the uncertainties that can affect a slope stability analysis. Using probability theory and statistics, the intrinsic spatial variability of soils and rocks and other soil and rock characteristics uncertainties can be considered in geotechnical analysis and/or design. This gives a balanced way to avoid these uncertainties in the analysis of slope stability study [23,24]. The reliability of a slope is explained as the probability that the performance standards will not be met. In other terms, the POF is the probability that the minimum FOS is less than one, i.e., POF = P(FOS < 1).
Assume that g(x) is a performance function used to determine the FOS, and that g(x) has a series of random variables x = [x 1 , x 2 , x 3 , . . . . . . , x k ] denoting uncertain model parameters. Thus, the value of β is given by: where µ FOS and σ FOS denote the mean and standard deviation of the FOS, respectively. The µ FOS signify the value of g(x) at mean values of µ 1 , µ 2 , . . . . . . , µ K random variables x 1 , x 2 , x 3 , . . . . . . , x k , given by: The σ FOS is defined as: where σ i denotes the random variable's standard deviations x i ; ρ i,j denotes the correlation coefficient between the two distinct unknowable variables, x i and x j ; ∂g ∂x j are the performance function's partial derivatives with respect to x i . Hence, the P f of slope can be obtained using POF = 1 − Φ(z), where Φ represents cumulative standard distribution function with µ = 0 and σ = 1.

Overview of Employed Models
This section provides a comprehensive overview of the computational models utilized in the present study. Firstly, a concise introduction to ANN, followed by a brief overview of meta-heuristic algorithms. Following this, the methodology for constructing a hybrid ANN is described and explained

Artificial Neural Network
ANNs are a subset of ML models that are constructed based on the principles of neuronal organization observed in biological neural networks found in animal brains. ANNs acquire information from the given dataset and produce predictions in response. ANN consists of one input layer, one or more hidden layers, and one output layer (see Figure 5). Neurons that are not computational collect information in the input layer, whereas computational neurons execute linear/non-linear computations in the hidden and output layers. Notably, weighted connections connect each and every neuron in the input, hidden, and output layers. The hidden and output layers both have biases that are proportional to the number of respective neurons. Input and output layers have a relationship between their respective numbers of input and output variables and their number of neurons. The number of hidden neurons varies from problem to problem and is selected through trial and error. In ANN, weights and biases are calculated iteratively to reduce the discrepancy between estimated and real values. The fitness function most frequently employed is the root mean square error (RMSE) index. Multiple training techniques, including conjugate gradient, gradient descent, Levenberg-Marquardt functions, etc., are used during the training of ANN and to adjust weights and biases. Application of ANNs in different engineering domains can be seen in the literature [51][52][53][54][55].
Mathematics 2023, 11, x FOR PEER REVIEW 8 of 24 In ANN, weights and biases are calculated iteratively to reduce the discrepancy between estimated and real values. The fitness function most frequently employed is the root mean square error (RMSE) index. Multiple training techniques, including conjugate gradient, gradient descent, Levenberg-Marquardt functions, etc., are used during the training of ANN and to adjust weights and biases. Application of ANNs in different engineering domains can be seen in the literature [51][52][53][54][55].

Overview of OAs
The complexity and difficulty of solving challenges in the actual world have increased during the past few decades. As a result, optimization techniques, especially meta-heuristic algorithms, have become important. These techniques, which approximate optimal solutions, rely on randomness. The optimization strategy determines the optimal decision variables through minimizing or maximizing the objective function. Due to their simplicity and ease of implementations, meta-heuristic algorithms have been deployed as an alternative solution. OAs are calorized into two main groups: (a) popbased, which considers the entire population, and (b) single-solution-based, which considers only one possible solution [56]. Notably, most of the pop-based OAs have been derived from natural phenomena. Figure 6 presents a classification of different kinds of pop-based OAs.

Overview of OAs
The complexity and difficulty of solving challenges in the actual world have increased during the past few decades. As a result, optimization techniques, especially meta-heuristic algorithms, have become important. These techniques, which approximate optimal solutions, rely on randomness. The optimization strategy determines the optimal decision variables through minimizing or maximizing the objective function. Due to their simplicity and ease of implementations, meta-heuristic algorithms have been deployed as an alternative solution. OAs are calorized into two main groups: (a) pop-based, which considers the entire population, and (b) single-solution-based, which considers only one possible solution [56]. Notably, most of the pop-based OAs have been derived from natural phenomena. Figure 6 presents a classification of different kinds of pop-based OAs. In ANN, weights and biases are calculated iteratively to reduce the discrepancy between estimated and real values. The fitness function most frequently employed is the root mean square error (RMSE) index. Multiple training techniques, including conjugate gradient, gradient descent, Levenberg-Marquardt functions, etc., are used during the training of ANN and to adjust weights and biases. Application of ANNs in different engineering domains can be seen in the literature [51][52][53][54][55].

Overview of OAs
The complexity and difficulty of solving challenges in the actual world have increased during the past few decades. As a result, optimization techniques, especially meta-heuristic algorithms, have become important. These techniques, which approximate optimal solutions, rely on randomness. The optimization strategy determines the optimal decision variables through minimizing or maximizing the objective function. Due to their simplicity and ease of implementations, meta-heuristic algorithms have been deployed as an alternative solution. OAs are calorized into two main groups: (a) popbased, which considers the entire population, and (b) single-solution-based, which considers only one possible solution [56]. Notably, most of the pop-based OAs have been derived from natural phenomena. Figure 6 presents a classification of different kinds of pop-based OAs.  The categorization of pop-based OAs can be based on their sources of inspiration, which include algorithms inspired by swarm intelligence (SI), algorithms based on human behavior, algorithms derived from evolutionary algorithms (EAs), and physics-based algorithms. SI algorithms replicate the social behaviors of organisms residing in swarms, flocks, and herds. The class of SI algorithms includes various approaches such as ALO [57], CPA [58], FF [59], GWO [60], PSO [61], and SSA [62], among others. The second group of OA pertains to human-based algorithms derived from the collective behavior and collaborative problem-solving approaches exhibited by individuals working in groups. The algorithms encompassed within this category include FDO [63], HS [64], ICA [65], TLBO [66], etc. EAs simulate natural evolutionary processes like mutation, selection, and recombination. BBO [67], DE [68], ES [69], and GA [70] are some of the examples of this group. The final category of popular algorithms draws inspiration from the rules of physics. BBBC [71], EO [72], GSA [73], and SA [74] are some examples of physics-based OAs. As stated above, eight OAs (viz., ALO, BBO, CPA, DE, EO, FF, GA, and PSO) from different groups were considered in this study, which are widely used in different engineering domains. Notably, detailed working principles of these OAs are not presented in this study because they are well established, and the original works of ALO [57], BBO [67], CPA [58], DE [68], EO [72], FF [59], GA [70], and PSO [61] can be referred to for more details.

Hybridization Procedure of ANN and OAs
Numerous studies have been done in the previous few decades on how to use multiple OAs to enhance the performance of conventional ANNs. The ANN may not perform well in some cases since back-propagation isn't great at locating precise global minima [51,75]. ANNs are susceptible to becoming trapped in local minima, unlike OAs, which can adjust their learning parameters to overcome this obstacle. OAs have the flexibility to perform exploration and exploitation techniques and generate optimal values for ANN weights and biases. The construction process of hybrid ANNs involves the following steps:

Data Description and Modelling
The geology of Mount St. Helens appeared more consistent before it collapsed due to volcanic eruptions than that of surrounding stratovolcanoes. Voight et al. [4] described the physical properties of debris avalanche materials. Intact edifice rock was found to have an average unit weight of 24 kN/m 3, and average values of and c were 40° and 1000 kN/m 2 , respectively. This information was used to perform the slope stability analysis. Following the normal distribution sampling technique, a total of 100 samples were generated using mean values of c = 1000 kN/m 2 , ϕ = 40°; and γ = 24 kN/m 3 . Subsequently, the FOS of the slope was determined using Scoops3D for the generated samples with five distinct values of seismic coefficient (ke), viz., 0 (for non-seismic case), 0.05, 0.10., 0.15, and 0.20, referred to as Set 1 to Set 5, respectively. Descriptive details of the 500 samples (i.e., 100 samples against each ke value) are presented in Table 1. The present study incorporates certain simplifications, including the assumption of homogenous material properties, intact rock mass, the exclusion of groundwater effects, and the inclusion of seismic loading. However, according to

Data Description and Modelling
The geology of Mount St. Helens appeared more consistent before it collapsed due to volcanic eruptions than that of surrounding stratovolcanoes. Voight et al. [4] described the physical properties of debris avalanche materials. Intact edifice rock was found to have an average unit weight of 24 kN/m 3, and average values of φ and c were 40 • and 1000 kN/m 2 , respectively. This information was used to perform the slope stability analysis. Following the normal distribution sampling technique, a total of 100 samples were generated using mean values of c = 1000 kN/m 2 , φ = 40 • ; and γ = 24 kN/m 3 . Subsequently, the FOS of the slope was determined using Scoops3D for the generated samples with five distinct values of seismic coefficient (k e ), viz., 0 (for non-seismic case), 0.05, 0.10., 0.15, and 0.20, referred to as Set 1 to Set 5, respectively. Descriptive details of the 500 samples (i.e., 100 samples against each k e value) are presented in Table 1. The present study incorporates certain simplifications, including the assumption of homogenous material properties, intact rock mass, the exclusion of groundwater effects, and the inclusion of seismic loading. However, according to Table 1 Table 1.
Stage I: After FOS estimations, the primary database of 500 records was finalized. This database was normalized randomly between 0 and 1 and then divided into training and testing subsets. Using 5-fold cross-validation, 80% of the entire dataset was used for training, i.e., 400 samples, whereas the remaining 20%, i.e., 100 samples, were used for testing. Notably, the training subgroup was used to construct hybrid ANNs, while the testing subgroup was used for validation. After model construction, multiple performance matrices, namely mean absolute error (MAE), Nash-Sutcliffe efficiency (NSE), performance index (PI), coefficient of determination (R 2 ), RMSE, and weighted mean absolute percentage error (WMAPE), were determined and assessed. Subsequently, the best-forming hybrid ANN was selected for the probabilistic assessment of Mount St. Helens.
Stage II: Subsequent to the selection of the best-performing model, RA was performed in seismic and non-seismic conditions. This was achieved by producing different sets of input parameters with different coefficient of variation (COV) values. The following stage involves the normalization of the new dataset based on the original input variables (see Table 1). The best-obtained paradigm was then used to generate FOS, followed by a probabilistic assessment of the slope. Figure 8 shows the entire process of FOS estimation and probabilistic analyses of Mount St. Helens. Stage I: After FOS estimations, the primary database of 500 records was finalized. This database was normalized randomly between 0 and 1 and then divided into training and testing subsets. Using 5-fold cross-validation, 80% of the entire dataset was used for training, i.e., 400 samples, whereas the remaining 20%, i.e., 100 samples, were used for testing. Notably, the training subgroup was used to construct hybrid ANNs, while the testing subgroup was used for validation. After model construction, multiple performance matrices, namely mean absolute error (MAE), Nash-Sutcliffe efficiency (NSE), performance index (PI), coefficient of determination (R 2 ), RMSE, and weighted mean absolute percentage error (WMAPE), were determined and assessed. Subsequently, the best-forming hybrid ANN was selected for the probabilistic assessment of Mount St. Helens.
Stage II: Subsequent to the selection of the best-performing model, RA was performed in seismic and non-seismic conditions. This was achieved by producing different sets of input parameters with different coefficient of variation (COV) values. The following stage involves the normalization of the new dataset based on the original input variables (see Table 1). The best-obtained paradigm was then used to generate FOS, followed by a probabilistic assessment of the slope. Figure 8 shows the entire process of FOS estimation and probabilistic analyses of Mount St. Helens.

Results and Discussions
This section provides comprehensive assessment of the results obtained from the slope stability analyses, the performance evaluation of hybrid ANNs, and the probabilistic analysis conducted on the slope. However, before presenting the outcomes of probabilistic assessment, the outcomes of slope stability analysis and parametric configurations of hybrid ANNs are presented and discussed.

Slope Stability Analysis
As stated above, the slope stability of Mount St. Helen was carried out using the Scoops3D computer program. For each 100 samples, the FOS was computed for five sets of k e , i.e., Sets 1 to 5. The outcomes of slope stability assessment are demonstrated in Figure 9. Additionally, Figure 10

Results and Discussions
This section provides comprehensive assessment of the results obtained from the slope stability analyses, the performance evaluation of hybrid ANNs, and the probabilistic analysis conducted on the slope. However, before presenting the outcomes of probabilistic assessment, the outcomes of slope stability analysis and parametric configurations of hybrid ANNs are presented and discussed.

Slope Stability Analysis
As stated above, the slope stability of Mount St. Helen was carried out using the Scoops3D computer program. For each 100 samples, the FOS was computed for five sets of ke, i.e., Sets 1 to 5. The outcomes of slope stability assessment are demonstrated in Figure 9. Additionally, Figure 10

Computational Modelling and Performance Assessment
The derived FOS values were used as the output variable, while the rock parameters viz., c, ϕ, γ, and ke. were used as the input variables for computational modelling. It is important to note that ANN`s hyper-parameters (i.e., NH and NHL) must be carefully selected to construct the best possible model. Using Levenberg-Marquardt backpropagation and tansigmoid activation functions and RMSE as the fitness function, the most appropriate value of NH was obtained as 8, i.e., NH = 8 and kept constant for other hybrid ANNs. However, deterministic parameters of OAs, i.e., NS, , ub, and lb, were also selected using trialand-error approaches. Since there are four inputs, the total number of were determined to be 49 (i.e., 4 × 8 + 8 + 8 + 1). The values of NS, , ub, and lb were selected as 50, 500, +1, and −1, respectively. The traditional ANN was also constructed using the above configurations in which RMSE was used as the fitness function.
On the contrary, evaluating the performance of a hybrid model requires examining the convergence behavior of OAs. It exhibits the ability of OAs to conquer local minimum and arrive at a solution more quickly. Convergence curves for all the hybrid ANNs utilized to calculate the FOS are shown in Figure 11. After 500 iterations, the computational cost of the hybrid ANNs were estimated to be 200.

Computational Modelling and Performance Assessment
The derived FOS values were used as the output variable, while the rock parameters viz., c, φ, γ, and k e . were used as the input variables for computational modelling. It is important to note that ANN's hyper-parameters (i.e., N H and N HL ) must be carefully selected to construct the best possible model. Using Levenberg-Marquardt backpropagation and tan-sigmoid activation functions and RMSE as the fitness function, the most appropriate value of N H was obtained as 8, i.e., N H = 8 and kept constant for other hybrid ANNs. However, deterministic parameters of OAs, i.e., N S , itr max , ub, and lb, were also selected using trial-and-error approaches. Since there are four inputs, the total number of O w+b were determined to be 49 (i.e., 4 × 8 + 8 + 8 + 1). The values of N S , itr max , ub, and lb were selected as 50, 500, +1, and −1, respectively. The traditional ANN was also constructed using the above configurations in which RMSE was used as the fitness function.
On the contrary, evaluating the performance of a hybrid model requires examining the convergence behavior of OAs. It exhibits the ability of OAs to conquer local minimum and arrive at a solution more quickly. Convergence curves for all the hybrid ANNs utilized to calculate the FOS are shown in Figure 11. After 500 iterations, the computational cost of the hybrid ANNs were estimated to be 200. ation between the actual and estimated values. Notably, for a perfect prediction model, the value of these indices should be equal to their ideal values as follows: MAE = 0, NS = 1, PI = 2, R 2 = 1, RMSE = 0, and WMAPE = 0. These indices are used widely for assessing the performance of data-driven models [28,76,77]. The detailed mathematical expressions of these indices are not presented in this study because they are well-established, and the studies of Bardhan and Samui [28] and Koopialipoor et al. [75] can be referred to for more details.  Tables 2 and 3 display the outcomes of the hybrid ANNs during the training and testing phases, respectively. As stated above, the training subset was used for model construction, while the testing subset was for validation of the constructed models. According to Table 2, the developed ANN-FF realized the desired level of precision (R 2 = 0.9996 and RMSE = 0.0041) during the training phase, followed by ANN-EO (R 2 = 0.9981 and RMSE = 0.0094), ANN-CPA (R 2 = 0.9970 and RMSE = 0.0117), and so on (see Table 2 for other models). Except for ANN-DE, all hybrid ANNs demonstrated a strong fit (R 2 > 0.95 or 95%), indicating a strong fit to the FOS database. Also, the employed ANN model exhibits satisfac- After model construction, the performance of the developed models was assessed using multiple performance matrices viz., MAE, NSE, PI, R 2 , RMSE, and WMAPE. Among these indices, MAE, RMSE, and WMAPE are the error indices, and NSE, PI, and R 2 are the trend-measuring indices. The use of numerous indices allowed the performance of the models to be assessed in various ways, including the degree of correlation/fitness, associated error in absolute terms, variance in error, and the amount of variation between the actual and estimated values. Notably, for a perfect prediction model, the value of these indices should be equal to their ideal values as follows: MAE = 0, NS = 1, PI = 2, R 2 = 1, RMSE = 0, and WMAPE = 0. These indices are used widely for assessing the performance of data-driven models [28,76,77]. The detailed mathematical expressions of these indices are not presented in this study because they are well-established, and the studies of Bardhan and Samui [28] and Koopialipoor et al. [75] can be referred to for more details.
Tables 2 and 3 display the outcomes of the hybrid ANNs during the training and testing phases, respectively. As stated above, the training subset was used for model construction, while the testing subset was for validation of the constructed models. According to Table 2, the developed ANN-FF realized the desired level of precision (R 2 = 0.9996 and RMSE = 0.0041) during the training phase, followed by ANN-EO (R 2 = 0.9981 and RMSE = 0.0094), ANN-CPA (R 2 = 0.9970 and RMSE = 0.0117), and so on (see Table 2 for other models). Except for ANN-DE, all hybrid ANNs demonstrated a strong fit (R 2 > 0.95 or 95%), indicating a strong fit to the FOS database. Also, the employed ANN model exhibits satisfactory performance (R 2 = 0.9950 and RMSE = 0.0159). However, among the developed ANNs, the ANN-DE was found to be the least effective model with R 2 = 0.9039 and RMSE = 0.0669.  The remaining dataset, referred to as the testing subset, was utilized to assess the model's capacity for generalization after its construction. The performance indicators utilized to evaluate the model's performance during the training phase were also determined for the testing subset and presented in Table 3. The outcomes reveal that the built ANN-FF achieved the highest precision with R 2 = 0.9996 and RMSE = 0.0042, followed by ANN-EO (R 2 = 0.9979 and RMSE = 0.0103), ANN-ALO (R 2 = 0.9969 and RMSE = 0.0124), and so on (see Table 3 for more details). In the testing phase, the performance of the employed ANN was determined to be R 2 = 0.9961 and RMSE = 0.0133, which is quite satisfactory compared to some of the hybrid ANNs. Nonetheless, the developed ANN-FF model attained the most precise performance in both phases of FOS estimation. For better demonstration, scatterplot, error plot, and line plot between the actual and estimated FOS values are presented in Figures 12 and 13 for the training and testing phases, respectively. Herein, the best-performing model, i.e., ANN-FF, is only illustrated. In the following sub-section, assessments of the POF of Mount St. Helen are presented.
ployed ANN was determined to be R 2 = 0.9961 and RMSE = 0.0133, which is quite satisfactory compared to some of the hybrid ANNs. Nonetheless, the developed ANN-FF model attained the most precise performance in both phases of FOS estimation. For better demonstration, scatterplot, error plot, and line plot between the actual and estimated FOS values are presented in Figures 12 and 13 for the training and testing phases, respectively. Herein, the best-performing model, i.e., ANN-FF, is only illustrated. In the following sub-section, assessments of the POF of Mount St. Helen are presented. ployed ANN was determined to be R 2 = 0.9961 and RMSE = 0.0133, which is quite satisfactory compared to some of the hybrid ANNs. Nonetheless, the developed ANN-FF model attained the most precise performance in both phases of FOS estimation. For better demonstration, scatterplot, error plot, and line plot between the actual and estimated FOS values are presented in Figures 12 and 13 for the training and testing phases, respectively. Herein, the best-performing model, i.e., ANN-FF, is only illustrated. In the following sub-section, assessments of the POF of Mount St. Helen are presented.

Assessment of POF
The results of probabilistic analyses of the volcanic slope of Mount St. Helen are presented in this sub-section. As previously indicated, the ANN-FF model, a superior hybrid ANN, was used to undertake the probabilistic assessments. Assuming a specific mean (µ) for different rock parameters, the standard deviation (σ) was calculated using: σ = µ × COV. The details of COV considered in this study are presented in Table 4. Note that the values of COVs considered in this study is in line with the COV range pre-sented in the existing literature. The steps of probabilistic analysis is organized as follows: (i) determination of µ and COVs for different rock parameters viz., c, φ, and γ; (ii) compute the σ FOS of the predicted FOS values; (iii) use rock properties and their standard deviations to produce random values; (iv) select additional input parameters such as k e ; (v) dataset finalization for five different k e values; (vi) dataset normalization as per minimum and maximum values of actual database (as presented in Table 1); (vii) estimation of FOS using the developed ANN-FF model, and (viii) the estimation of β and POF as per Section 4.3.
Note that, five different combinations of COV (Cases 1 to 5) were considered in this study. The COVs of c and φ were set in the range of 10-30% by 5% and 2-10% by 2%, respectively. However, a constant COV of 5% was considered for parameter γ. Table 4 displays COV values for all five cases. To perform probabilistic analysis, a total of 30 samples were generated against each COV case, resulting in a total of 150 (i.e., 30 × 5) samples. These samples were investigated against five k e values. Therefore, 25 instances (i.e., 5 COV cases and 5 k e sets) were investigated. The results of probabilistic analyses are presented in Table 5 According to the information presented in Table 5 and Figure 14, β of the Mount St. Helen reduces when COV and k e increase, indicating higher POF. Also, the impact of COV on the POF was observed throughout a range of k e values. Thus, it can be stated that the selection of appropriate values for the COVs has a major impact on the slope failure probability. Using GA and MCS, Tun et al. [38] conducted RA of the same slope of the Mount St. Helen. In their studies, the authors assumed material homogeneity and ignored the pore pressure ratio and seismic loading. Using a single COV of 23% for c and 7% for φ, Tun et al. [38] reported β = 6.14 and POF = 0.0%, i.e., very low POF.
Similarly, using the MCS approach, the β and POF values were reported as 5.89 and 0.0%, respectively. Nonetheless, details of generated samples and a sample calculation for β and POF are presented in Appendix A. Herein, the samples generated with COV Case 1 are only presented. However, the values of β and POF are shown for all the five k e combinations. The details of weights and biases of the developed ANN-FF model are given in Appendix B, which can be used to assess β and POF for other cases. Moreover, the developed MATLAB model of ANN-FF and DEM profiles are included as Supplementary Materials for future studies.    Figure 14. Illustration of (a) β and (b) POF at different COV levels. Figure 14. Illustration of (a) β and (b) POF at different COV levels.

Summary and Conclusions
RA was carried out on the real topography of Mount St. Helens, which collapsed in 1980 due to an eruption, and the results are described in this study. Initially, the FOS was determined by performing a series of calculations in Scoops3D with various input parameters, such as rock characteristics that vary with location and seismic coefficient. The next step was to construct eight hybrid ANNs, using FOS as an output and the corresponding rock parameters and k e as inputs. After performance assessment, the best-performing model, i.e., ANN-FF (based on performance in the testing phase with R 2 = 0.9996 and RMSE = 0.0042), was chosen for probabilistic analyses the Mount St. Helen in seismic and non-seismic scenarios.
For this purpose, five distinct COV combinations were investigated for five different k e values, varying between 0 and 0.20. According to the experimental results, the POF varies in the range of 1.1 × 10 −44 to 1.7 × 10 −04 and 2.1 × 10 −38 to 0.19 in non-seismic and seismic conditions, respectively. Even at the high COV and seismic levels, the POF was between 9.7 × 10 −20 and 0.19. These results indicate that the failure probability of the slope is negligible even at higher COV levels. Thus, it is deduced that if Mount St. Helens does not erupt, COVs will have little effect on the POF. Given the uncertainty of rock properties, the suggested ANN-FF-based RA is determined to be an appropriate solution for calculating the POF of the Mount St. Helen.
Furthermore, the proposed technique demonstrates an accurate estimation of the FOS of the slope irrespective of seismic conditions. The proposed ANN-FF has the following advantages: (i) improved generalization, (ii) faster convergence, and (iii) higher prediction accuracy in both phases. However, the suggested ANN-FF model has higher computational cost than other hybrid ANNs built in this study. In addition, the above investigation was carried out using the upper and lower ranges of rock parameters (as detailed in Table 1). Hence, it can be deduced that the above analysis may not yield satisfactory results beyond these values, which can be considered as one of the limitations of the present study. However, further assessment is required for this case. Therefore, the following points should be considered as the future scope of the study: (i) implementation of improved mechanism to reduce the computational cost of ANN-FF paradigm; (ii) indepth assessment of slope failure probability at high COV and larger range of c, Ø, and γ parameters, (iii) adoption of Mononobe-Okabe trick combined with seismic actions for an in-depth assessment, and (iv) a comprehensive assessment of other hybrid models of ANN, RVM, and ELM constructed with different group of OAs. Nevertheless, per the author's knowledge and literature review, this study is the primary implementation of the ANN-FF model to perform probabilistic analyses of the Mount St. Helens in seismic and non-seismic conditions.

Appendix A. Sample Calculation
Details of generated samples and estimation of β and POF are presented in Table A1.