The research focuses on the Baspa Valley in Kinnaur district of Himachal Pradesh, one of the most landslide-prone corridors in the Indian Himalaya. The valley is traversed by the Sangla Detachment Fault, a major tectonic structure that fractures and weakens the rock through which it passes, making road cuts along the valley exceptionally susceptible to recurring slope failures. The region also falls within Seismic Zone IV of India’s seismic zonation, meaning engineers must consider strong earthquake shaking, represented in the study by a horizontal seismic coefficient of 0.12 g, in addition to the ordinary pull of gravity.

To build a trustworthy dataset, the researchers combined field investigations, laboratory testing and published literature to characterise the rock masses along the valley. Nine key geotechnical parameters went into their models: slope geometry, the Geological Strength Index, or GSI, the uniaxial compressive strength of the intact rock, Young’s modulus, Poisson’s ratio, the Hoek-Brown material constant, the disturbance factor and the seismic loading condition. These parameters describe everything from how fractured the rock mass is to how much blasting or excavation damage it has suffered, and they form the vocabulary in which any slope-stability problem can be expressed.

Using the finite element method, a numerical technique that divides a slope into thousands of small elements and solves the equations of stress and strain across them, the team ran 494 separate slope stability simulations under both static and seismic conditions. Each simulation produced a factor of safety, the classic engineering ratio that compares the resisting forces in a slope to the driving forces trying to pull it down, with values below 1.0 signalling imminent danger. A parallel set of seismic factors of safety was generated for the same slopes under pseudo-static earthquake loading, giving the researchers two rich datasets on which to train their predictive models.

With the simulation data in hand, the team benchmarked a suite of conventional machine-learning algorithms. Gradient Boosting emerged as the strongest of the traditional approaches, achieving a coefficient of determination on the test set of 0.8366 for the static factor of safety and 0.8349 for the seismic version. Random Forest followed with scores of 0.7973 and 0.7299 respectively, while Support Vector Regression and Decision Trees performed moderately, with the Decision Tree in particular showing signs of overfitting, and the K-nearest neighbours method lagging behind the rest of the field.

Feature-importance analyses drawn from the tree-based models delivered one of the study’s most consequential findings: the Geological Strength Index accounted for more than half of the predictive power across the models. The GSI is a field-based visual assessment of how blocky, fractured and disturbed a rock mass is, and its dominance means that careful geological mapping remains the single most valuable input an engineer can provide. Slope angle proved to be the critical secondary variable for static stability, while slope height gained prominence once seismic loading was introduced, a shift that mirrors how taller cuts amplify inertial forces during an earthquake.

The real leap in accuracy, however, came from artificial neural networks. A baseline ANN reached a test R-squared of 0.9019 for the static factor of safety and 0.7864 for the seismic case, already outperforming every conventional method. The researchers then turned to metaheuristics, nature-inspired optimisation algorithms that search for the best possible set of neural network weights and hyperparameters, in the same way that a bee colony optimises its foraging or an ant colony finds efficient paths. Three optimisers were tested: the Artificial Bee Colony, Ant Colony Optimization and the Genetic Algorithm.

Each metaheuristic pushed performance beyond the baseline network. The Artificial Bee Colony-tuned ANN achieved the best static prediction, with a test R-squared of 0.9376, while also reaching 0.8983 under seismic loading. Ant Colony Optimization delivered 0.9316 and 0.8758 for the two cases, and the Genetic Algorithm-tuned network posted 0.9301 for static conditions and an outstanding 0.9178 for seismic conditions, making it the top performer where earthquake shaking matters most. These figures indicate that the optimised networks reproduce the output of full finite element simulations with errors of only a few percent, while evaluating a new slope in a fraction of a second rather than the minutes to hours a numerical model requires.

That speed-to-accuracy trade-off is the heart of GeoRockSlope. The application, written in Python and deployed on the Streamlit platform, is freely available on GitHub and lets practitioners enter the nine geotechnical parameters for any road cut and receive both a static and a seismic factor of safety almost instantly. Because the underlying models were trained on physics-based finite element results rather than on sparse case histories, the tool effectively distils the rigour of numerical modelling into a field-ready package. It eliminates the need for primitive empirical classification charts on one hand and for repeated, computationally intensive FEM analyses on the other, and the study notes that the machine-learning code and the web application are both openly accessible for verification and reuse.

The stakes of such a tool are far from academic. Landslides are among the deadliest natural hazards in mountainous regions worldwide, and roughly 12 percent of India’s land area is considered prone to them, with the Himalaya bearing a disproportionate share of the losses. In a tectonically fractured valley like the Baspa, where new road cuts are carved regularly and every monsoon brings fresh failures, a prediction engine that can screen hundreds of slope geometries in an afternoon could reshape how highway engineers prioritise stabilisation work, where to place rockfall barriers, and which sections to monitor during seismic events. For communities perched along these fragile corridors, the difference between a slow expert assessment and an instant, well-calibrated forecast may ultimately be measured in lives saved.

Subject of Research: Machine learning prediction of rock-slope stability along seismically active Himalayan road cuts

Article Title: GeoRockSlope: an open access application for prediction of FEM-ANN guided road cut rock-slope instability

Article References: Pandey, V. H. R., Kainthola, A., & Kushwaha, G. (2026). GeoRockSlope: an open access application for prediction of FEM-ANN guided road cut rock-slope instability. Natural Hazards, 122(20), Article 654. https://doi.org/10.1007/s11069-026-08407-z

Image Credits: AI Generated

DOI: 10.1007/s11069-026-08407-z

Keywords: GeoRockSlope, slope stability, machine learning, artificial neural networks, finite element method, factor of safety, Himalaya, seismic hazard, artificial bee colony, Genetic Algorithm, Gradient Boosting, road cuts