Abstract

A sudden downward movement of the geomaterial, either composed of soil, rock, or a mixture of both, along the mountain slopes due to various natural or anthropogenic factors is known as a landslide. The Himalayan Mountain slopes are either made up of residual soil or rocks. Residual soil is formed from weathering of the bedrock and mainly occurs in gentle-to-moderate slope inclinations. In contrast, steep slopes are mostly devoid of soil cover and are primarily rocky. A stability prediction system that can analyse the slope under both the condition of the soil or rock surface is missing. In this study, artificial neural network technology has been utilised to predict the stability of jointed rock and residual soil slope of the Himalayan region. The database for the artificial neural network was obtained from numerical simulation of several residual soils and rock slope models. Nonlinear equations have been formulated by coding the artificial neural network algorithm. An android application has also been developed to predict the stability of residual soil and rock slope instantly. It was observed that the developed android app provides promising results in predicting the factor of safety and stability state of the slopes.

1. Introduction

Landslide is defined as the downslope movement of a large mass of rocks, rubbles, soil, or a mixture of all, owing to various natural and anthropogenic factors. This downslope movement of soil and rock through the influence of gravity leads to a kind of “mass wasting” [1]. Landslide consists of five various types of slope movement: falls, topples, slides, spreads, and flows [2]. These types are segmented based on their geological make (bedrock, debris, or Earth). Some of the common landslide types are debris flows and rockfalls. Nearly every single landslide has numerous reasons for its occurrence [3, 4]. Slope movement occurs once forces operating down the slope under the influence of gravity surpass the strength of the ground materials that hold the slope [5]. Some of the major factors increasing the landslide occurrence are the impact of downslope forces and the reduced strength of the bedrock or the soil constituting the slope. Landslides are set off on the slopes, which are already on the rim of movement by rainfall, snowmelt, changes in water level, stream erosion, changes in groundwater, earthquakes, volcanic activity, disruption by human activities, or any permutation of these factors [6].

Landslides mostly result in property damage, injuries, and many fatalities [7]. It also negatively alters a variety of natural resources. Water bodies, fisheries, sewage disposal systems, forests, dams, and roadways the whole lot have an effect for a long time after a landslide occurrence [8]. These adverse economic effects of landslides incorporate the expenses for the restoration of structures, loss of property, interruption of transportation routes, and medical expenses in the instance of injuries, and unforeseen costs in the form of lost goods, freshwater accessibility, level, and condition are also impacted [9]. People might decrease their risks of landslides by studying probable landslide hazards and undertaking measures to mitigate those hazards. Due to the isolated nature of the landslide event, they usually occur in an instance without any public forewarning.

Landslide risk mitigation is generally categorised into direct and indirect methods [10]. The direct method includes geometric methods, where the geometry of the hillside is altered (in general the slope); hydrogeological methods, where an attempt is made to decrease the groundwater level or to reduce the water content of the slope material; chemical and mechanical methods, where efforts are made to build up the shear strength of the unsteady mass or else to establish dynamic external forces or passive to counter the destabilising forces. These methods adjust to a certain extent with the kind of material that makes up the mountain slope [1113]. The indirect method of landslide risk mitigation includes landslide hazard zonation maps, landslide hazard charts, and early warning techniques, which is generally used to avoid vulnerable areas for developmental works [1417].

The Himalayan Mountain slopes are either made up of residual soil or rocks. Residual soil is formed from weathering of the bedrock and mostly occurs in gentle-to-moderate slope inclinations (<50°–60° slope angle). In contrast, steep slopes (>60° slope angle) are mostly devoid of soil cover and are primarily rocky [1]. Many times, there is a sudden change in the slope geomaterial (soil or rock) within a few hundred meters of distance. The stability analyses of both cases of slope material are different from each other and cannot be used uniformly. A stability prediction system that can analyse the slope under both the condition of the soil or rock surface is missing.

A detailed study of the residual soil and rock slopes of the Himalayan region was carried out by Ray et al. [18] and Ray et al. [19], respectively. In the present study, the artificial neural network analysis was used for stability prediction of residual soil and rock slope using the numerical simulation data of Ray et al. [18] and Ray et al. [19]. Nonlinear equations have been formulated by coding the artificial neural network algorithm. An android application has also been developed to predict the stability of residual soil and rock slope instantly. The mobile application was rendered user friendly with a basic GUI and is available for public use to quickly identify the stability state of the residual soil and rock slope of the Himalayan region.

2. Area of Study

The study area is located in the Himalayan region in the Indian states of Uttarakhand and Himachal Pradesh, which lies at the convergence zone of two lithospheric plates. One is the Indian plate in the south, and the other is the Eurasian plate in the north. The prevalence of a large number of discontinuities makes this region geologically very fragile and susceptible to frequent landslides. A slight imbalance in elements of shear stress and strength factors might trigger landslides. The Himalayan Mountain slopes can be classified into two distinct categories based on stratigraphy, namely, rock slopes and soil slopes. The rock slopes have highly fractured rock mass at the surface and are devoid of any soil cover. The failure in rock slopes is governed mainly by the discontinuities present within the rock mass. The soil slopes, on the other hand, are formed from the weathering of bedrock and are also referred to as residual soil slope. Slope failure in residual soil slope is a complex phenomenon involving various factors such as slope topography, depth of soil cover, the grain size distribution of the soil, inherent heterogeneity in the geo-mechanical properties of the soil, and the presence of water. Ray et al. [1] have carried out extensive literature studies of the Himalayan region and concluded that up to moderate slope (<50°–60°), the thickness of residual soil varies from 2 to 10 m, with some places exceeding 10 m before encountering the weathered bedrock. Previous studies concerning slope stability were limited either to residual soil slopes or rock slopes. The current study focuses on the slope stability problem associated with jointed rock slopes and residual soil slopes in the Himalayan Range.

A database of 400 slope models, which were previously analysed by Ray et al. [18] using a numerical simulation technique, has been used for the stability analysis of the residual soil slope. Thirteen major influencing parameters have been considered during numerical simulation, which include Young’s modulus of residual soil (Es), unit weight of soil (ɣs), shear strength parameter of residual soil (cohesion (Cs) and angle of internal friction (Φs)), Young’s modulus of the weathered rock mass (Er), shear strength parameter of the weathered rock mass (cohesion (Cr) and angle of internal friction (Φr)), unit weight of weathered rock (ɣr), strength parameter of the soil-rock joint interface (cohesion (Cj) and angle of internal friction (Φj)), average slope angle (α), slope height (H), and residual soil depth (D). The obtained results indicate that all the parameters are affecting the FOS of the residual soil slopes. However, the impact of all those thirteen parameters was not uniform as can be seen from correlation analysis (Table 1). For developing the ANN models, all those parameters whose correlation coefficient is more than 0.2 are considered.

Similarly, a database of 1200 slope models, which were previously analysed by Ray et al. [19], using a numerical simulation technique, has been used for the stability analysis of the rock slope. Thirteen major influencing parameters have been considered during numerical simulation, which include Young’s modulus of the weathered rock mass (Er), the tensile strength of the weathered rock mass (Tr), Poisson's ratio of weathered rock mass, unit weight of weathered rock mass (ɣr), shear strength parameter of the weathered rock mass (cohesion (Cr) and angle of internal friction (Φr)), shear strength parameter of the bedrock (cohesion (Cbr) and angle of internal friction (Φbr)), slope angle (α), slope height (H), width of weathered layer (D), main joint orientation (Dip1), and cross joint orientation (Dip2). The obtained results indicate that all the input parameters affect the FOS of the rock slopes. However, the impact of all those thirteen parameters was not uniform as can be seen from correlation analysis (Table 2). For developing the ANN models, all those parameters whose correlation coefficient is more than 0.2 are considered.

3. Methodology

Four hundred numerical simulation models were analysed by Ray et al. [18] using various permutation combinations of the height of the slope, the angle of the slope, and the thickness of the layer of residual soil. The height of the slope varied from 50 m to 500 m, with an interval of 50 m. The overall slope angle varied from 15° to 60°, with an interval of 15°. The thickness of the layer of the residual soil varied from 0.5 m to 15 m (0.5 m, 1 m, 2 m, 3 m, 4 m, 5 m, 7 m, 9 m, 12 m, and 15 m). Thus, the height of the slope resulted in ten variations, the thickness of residual soil resulted in ten variations, and the angle of slope caused four variations, leading to 400 models (i.e., 10 ∗ 10 ∗ 4 = 400). About 1200 numerical simulation models were analysed by Ray et al. [19] using various permutation combinations of slope height, slope angle, weathered layer thickness, and variations of each cross joint set (models were prepared with two joint sets). There were five variations in slope height, three variations in slope angle, five variations in weathered layer thickness, and four variations of each cross joint set. Using the numerical simulation data obtained from analysing the residual soil slopes and the rock slopes, ANN models were prepared to predict the FOS of the investigated slopes.

3.1. Artificial Neural Network Model

The natural nervous system, such as the brain processing information, is the main inspiration behind the information processing paradigm of an ANN (artificial neural network). It comprises many highly interconnected processing elements (neurons) functioning in concurrence to resolve a particular problem [20]. The artificial neural network is a framework instead of an algorithm for various distinct machine learning algorithms, through which they might work together and process complex computations [21]. These types of systems tend to “learn” in order to implement tasks by ruminating the examples, mostly lacking the need of being programmed with any task-specific instructions. An ANN primarily moves through a training phase, where it figures out how it might perceive the patterns in data, whether using visual, audio, or literal aid [22]. The disparity between the two outputs will be corrected by utilising backpropagation [23]. This further implies that the network will work in reverse, going from the output nodes to the input nodes in order to alter the weight of the connections amongst the nodes up until the variation between the actual and desired outcome yields the minimal conceivable error. Throughout the training and testing phases, the ANN is trained on what things to search for and what its output ought to look like, utilising binary yes/no query types [24].

The activation function is imperious for an ANN to learn and build the logic of some genuinely complex problem. The crucial notion here is to convert an input signal of a node in an ANN into an output signal of the node. The output signals might be utilised as an input to the next layer in the neural network stack. The activation function decides if a neuron ought to be activated by assessing the weighted sum and then further adding up the bias [25]. This is done to get nonlinearity into the output of a node. If there is no activation function, then the output signal will be a more straightforward linear function, which is limited in its complexity and has much less learning power [26]. Deficient in a proper activation function, the model will not learn complex data such as images, videos, audio, and speech. Nonlinear functions have a higher degree, that is, more than one, and they follow a polynomial curve. Some essential types of activation functions are threshold activation function, binary step function, sigmoid activation function, hyperbolic tangent function, and rectified linear units [27].

3.2. Network Architecture

The neural network is generally formulated into various layers. In Figure 1, the input layer is represented in yellow. It directly takes the input from the external source. Each of the nodes in the input layer is then interconnected with the nodes of the next layer. Blue and green colours denote the hidden layers. Hidden layers are interconnected networks of nodes that exist between the output and the input nodes. Hidden layers remain hidden from the outside world and are responsible for the computation done by the neural network. Edges are formed by the connection of nodes between two different layers, and each edge is allotted some weight based on the relative significance of the connection between the nodes. All the node values are based on the inputs provided by the user. These are then multiplied with the assigned weight of their edges, and then it is all summed up to get the value of the particular, mostly done for the hidden layers [21, 28]. Based on this summed-up value, nodes of the hidden layer will be “activated” or not, which is decided by predefined activation functions. Thus, the hidden layers take the input from the input layer and compute based on the weights and the activation function. Then, that result is propagated to the next layer of nodes until it reaches the final output node, where one can view the outcome [28]. Various performance indices such as the coefficient of determination (R2), residual error, root mean square error (RMSE), the variance accounted for (VAF), and the learning rate are used for assessing the efficacy of the developed ANN models [21, 23, 27, 29, 30].

3.3. Mobile Application

Developing a mobile application means developing a software application for smartphones, tablets, and other handheld devices. This application can be preinstalled in mobiles or supplied on a server or client-side as web applications to give a “mobile application-like” familiarity in a web browser portal. This is possible with an application-like processing experience. Mobile application developers consider a long arrangement of screen sizes, hardware specifications, and configurations before proceeding with the development process. The mobile application developers should also consider the extensive diversity of screen sizes, hardware specs, and settings due to fierce competitiveness in mobile applications and changes across all platforms. Being a part of the mobile app development process, mobile user interface (UI) designing will be a vital step in the making of mobile apps. The mobile app UI design's objective is primarily to make the experience of the interface more comprehensible. Much more functionalities will be followed by inculcating mobile enterprise application platforms or integrated development environments (IDEs).

Utilising an artificial neural network (ANN) technique, a microstrip patch antenna will be constructed. In order to design a proper mobile application, both the investigation and creation problems for constructing microstrip patch antennas are taken into. An analysis problem denotes the estimate of resonant frequency, while a synthesis problem denotes the estimate of dimensions. Both these are reciprocal of each other. The outcomes are put into operation utilising a mobile application. Backpropagation of the artificial neural networks is mainly utilised for training the neural network in order to minimise the error and computation time. So, the geometric dimensions of the micro antenna patch will be obtained with higher accuracy in lesser computation time in comparison to the simulation software. Sometimes, the size of the antenna might be a limitation, especially in mobile and wireless applications.

The work required dealing with low-level graphics, high accuracy prediction using ANN, and the vast portion of the app, React Native and Flutter, did not fit the use case, and it is just too much overhead to build the app with C++. Thus, the viable option is to create separate apps for each platform or use Kotlin Multiplatform to share most code. In order to create separate apps, an Android app will be created using Kotlin, and then the code will be shared with iOS, too (Figure 2). It has low risk and high returns. Therefore, Kotlin Multiplatform (KMP) is the obvious choice. The KMP-based app can scale well and would be beneficial in the long run, even if some problems with tooling and language might occur. Kotlin Multiplatform is created by Jetbrains and uses Kotlin as the programming language. It is more a set of technologies than a framework. It provides common business logic, high performance, compiles to bytecode for Android and native for iOS, good interop with Java, and very little to no change for Android developers.

Kotlin Multiplatform Mobile or KMM is a subset of KMP focusing on making the tech stable or used on mobile platforms, namely, Android and iOS and provides improved and better tooling for this use case. Kotlin Multiplatform provides an advantage over other platforms in making scalable and maintainable apps using a single codebase. Although the code-sharing is slightly less than that with other cross-platform solutions, UI cannot be shared; this gives the advantage of creating a user experience that aligns with the platform’s guidelines and best practices. KMP can target JVM for a faster iteration cycle as the whole code can be run and tested on the PC rather than a mobile device. It also provides the future possibility of using the same code for other platforms such as Web and Desktop with minimal changes.

4. Prediction of Stability Class and Factor of Safety

The sample dataset comprises information to build the predictive model. The inputs are independent variables, the targets are the dependent variables, and the unused variables will not be used as inputs or targets. Furthermore, some instances might be training instances that are utilised to create the model, selection instances are utilised for picking the optimum order, testing instances are utilised to validate the model's working, and unused instances are not contemplated at all.

4.1. Prediction of Stability Class or Risk Factor for Residual Soil Slope

The total number of instances is 400. The number of training instances is 240 (60%), the number of selection instances is 80 (20%), the number of testing instances is 80 (20%), and the number of unused instances is 0 (0%). Table 3 lists the number of variables and their type. Table 4 lists the minimum, maximum, mean, and standard deviation for all of the variables in the dataset. In this model, the FOS is the output. As a consequence, no mean and standard deviation are shown for them. The minimum column shows the class of the nominal variable with the least number of appearances. In contrast, the maximum column shows the class of the nominal variable with a significant number of impressions.

Slope angle, slope height, soil cohesion, soil friction angle, Young’s modulus of soil, and soil depth are used as inputs in the artificial neural network. FOS class four, FOS class three, FOS class two, and FOS class one are used as target outputs. The output class is defined based on the FOS range/class for a particular slope profile. This range/class is listed in Table 5.

The most suitable neural network architecture was deduced by training, testing, and validating different combinations of hidden layers and associated neurons. A single hidden layer having a range of 2–15 neurons has been opted to train the neural network model. The number of neurons in the rest two layers, that is, the input layer and output layer, are constrained to the number of input variables and output variables, respectively [31]. The normalising of the input data was performed and the training process is initiated by randomly varying the number of neurons and their associated weights. The network with 6-3-4 architecture (six neurons in the input layer, three neurons in a single hidden layer, and four neurons in the output layer) is selected based on maximum training performance and minimum training error for the developed ANN model after simulation of various network combinations (Table 6). The selected ANN model has a learning rate and momentum of 0.91 and 0.014, respectively, and its architecture is shown in Figure 3.

4.2. Prediction of Factor of Safety for Jointed Rock Slope

The data set includes valuable information for constructing a predictive model. There is a data matrix where columns denote variables and rows denote instances. The total number of instances is 1200. The number of training instances is 720 (60%), the number of selection instances is 240 (20%), the number of testing instances is 240 (20%), and the number of unused instances is 0 (0%). If the data are very irregularly distributed, the model will probably become inferior quality. Table 7 lists the number of variables and their type. Table 8 lists the minimum, maximum, mean, and standard deviation of all the variables in the data set.

The most suitable neural network architecture was deduced by training, testing, and validating different combinations of hidden layers and associated neurons. A single hidden layer having a range of 2–15 neurons has been opted to train the neural network model. The number of neurons in the rest two layers, that is, the input layer and output layer, are constrained to the number of input variables and output variables, respectively [31]. The normalising of the input data was performed and the training process is initiated by randomly varying the number of neurons and their associated weights. The network with 5-10-1 architecture (five neurons in the input layer, ten neurons in a single hidden layer, and one neuron in the output layer) is selected based on maximum training performance and minimum training error for the developed ANN model after simulation of various network combinations (Table 9). The selected ANN model has a learning rate and momentum of 0.96 and 0.011, respectively, and its architecture is shown in Figure 4.

4.3. Efficiency of the Developed ANN Models

Several performance indices, as discussed in detail by Ray et al. [21] and Karir et al. [28], were used to measure the efficiency of the developed ANN models. The coefficient of determination (R2) between the target and predicted values indicates a good prediction performance of the developed models. The RMSE and variance account for (VAF) were computed for studying the performance and the prediction capacity of the developed predictive models. Various performance indices obtained from the developed ANN1 model are presented in Table 10 and Figure 5.

5. Android Application

The developed app has two modules: one is residual soil classification, and the other is rock slope failure prediction using an artificial neural network. The neural network simulation has been done using neuro software. The simulation has been exported into code, and an equation has been used in developing the android app. Simulation of ANN also takes time and memory depending upon model complexity. Therefore, the equation generated from ANN simulation to predict classification and factor of safety was used. Smartphones have turned popular for the research for mobile applications because of their enhanced performance. This mobile application will be developed to facilitate the professionals to get a quick understanding of residual soil and rock slope behaviour to ascertain the potentially dangerous slopes. The smart miner app predicts the risk of slope failure for residual soil and rock, the factor of safety of rock slope, and the classification of residual soil (https://play.google.com/store/apps/details?id=com.maheshpaliwal.smartmineapp).

5.1. Debugging the Android Application

The android application offers dedicated tools to debug the app in Android Studio (IDE). There are two ways an Android app will be debugged: first is the android simulator on the IDE itself and then the second is by testing the app on an actual Android smartphone. In this case, debugging the apps will be done in both ways and will them make stable and as reliable as possible.

5.2. Overview and Design of GUI

The app is built user friendly with a very basic GUI (Figure 6). The sample results are compared with various situations, and a very nice accord is achieved amongst them. The app helps the professionals to get a quick idea about residual soil and rock slope behaviour in order to identify the potential hazardous slopes.

5.3. Analysis of Residual Soil Classification for Safety Class

The inputs are slope angle, slope height, shear strength of residual soil, Young’s modulus, and soil depth. The first step is that all parameter is scaled-down, and then y_1, y_2, and y_3 are calculated using the equation. The output is the probability of failure for all classes, that is, very low, low, medium, and high, is determined using the below-mentioned equation. The risk was predicted using the if-else rule and displayed in GUI using the Android application:

The generated parametric equations are as follows:

The final output equations are as follows:

5.4. Analysis of Stability of Rock Slope

The inputs are slope angle, slope height, weathered layer thickness, joint one, and joint two. The first step is that all parameter is scaled-down, and then y_1 to y_10 are calculated using the equation generated. The final dependent variable is calculated using the tanh function and the set of generated equations. Equation 1 has been used to calculate the factor of safety of rock slope. The factor of safety is displayed in GUI using the Android application:

The generated parametric equations are as follows:

The final generated parametric equations are as follows:-

The final output equations are as follows:

5.5. Prediction of Stability Class and Factor of Safety Using Android Application

After verifying and debugging, the developed app is tested on three real android mobiles; Samsung (model: SM-N750), Asus (model: ASUS T_00J), and Xiaomi (model: Mi Note 8 Pro); it is concluded here that a good agreement is achieved between computed results and their reference counterparts for all seven cases associated with the example. The values for a particular case are listed in Tables 11 and 12. A particular case is presented in these tables. The app results are in normalised form by default, which can be converted before use. The app has two modules: one is residual soil classification and the other is rock slope failure prediction using an artificial neural network. Table 11 lists the prediction classification of residual soil using the ANN algorithm and android app. Table 12 lists the prediction of FOS of rock slope using the ANN algorithm and android app.

6. Conclusions

There are various problems associated with calculating the factor of safety of the hills in the Himalayan region. Various tests are to be performed, which are very time consuming and expensive. As by nature, the slopes are not uniform. There can be a sudden change in the geomaterial (soil to rock and vice versa) within a few hundred meters. Thus, carrying out conventional geotechnical investigations over a stretch having several rocks and soil slope is very cumbersome. An effort has been made to quickly identify the stability of both the rock slope and residual soil slope. Numerical simulation results from previous studies on residual soil and rock slope were used to develop the ANN models. Various performance indices were used to highlight the high efficiency of the developed ANN models. The developed ANN models were further used to develop a mobile application for slope stability assessment.

The proposed software is user friendly and economical (virtually no cost) and has been constructed on an artificial neural network application intended in place of microstrip antenna applications. The performance viability of the developed mobile application has been verified using the virtual android platform in the IDE itself and then utilising various real-world android mobiles possessing unique chipset configurations ranging from higher performance to lower performance. All the mobile devices performed quite positively, providing the same desired results for all use cases. Thus, the proposed solution delivers a roadmap for the development of an android app utilising artificial neural networks due to which the proposed methodology is reasonably faster and has a lower manufacturing cost.

Data Availability

Data available on request

Conflicts of Interest

The authors declare that they have no conflicts of interest.