Data-Driven Deformation Prediction of Accumulation Landslides in the Middle Qinling-Bashan Mountains Area

: Accurately predicting landslide deformation based on monitoring data is key to successful early warning of landslide disasters. Landslide displacement–time curves offer an intuitive reflection of the landslide motion process and deformation predictions often reference the Saito curve for correlational analysis with cumulative deformation curves. Many scholars have applied machine learning techniques to individual landslide deformation predictions with considerable success. However, most landslide monitoring data lack a full lifecycle, making it challenging to predict unexperienced evolutionary stages. Cross-learning between similar landslide datasets provides a potential solution to issues of data scarcity and accurate prediction. First, this paper proposes a landslide classification and displacement machine learning method, along with predictive performance evaluation metrics. Further, it details a study of 13 landslides with evident deformation signs in the middle Qinling–Bashan Mountains area, conducting refined landslide classification. Based on a data-driven approach, this study conducts an analysis of the importance of characteristics influencing landslide deformation and establishes predictive models for similar-type landslide deformation, mixed-type landslide deformation, and individual landslide deformation using machine learning algorithms. The models trained on the dataset are used to predict the deformation of the West of Yinpo Yard landslide at different periods, with the predictive performance evaluated using two indices. The results indicate that the models trained on similar-type landslide data and those based on individual landslide data yielded comparable predictive performances, substantially addressing challenges such as insufficient early-stage monitoring data and low prediction accuracy.


Introduction
Due to the movements of the Pacific Plate, the Indian Ocean Plate and the Asian-European Plate, China is characterized by strong tectonic and seismic activities, and complex topography and geological conditions, coupled with a variety of climatic types and intense human engineering activities, which make geologic hazards prone to occurring, prolific, and frequent.The middle Qinling-Bashan Mountains area is an area highly prone to geologic hazards, and is located in the key prevention and control area of landslides, avalanches, and mudslides in the Qinling-Bashan Mountains of China.The geological disasters in the region are characterized as cluster, regional, and chain-born.The distribution ent influencing factors.There are many existing classification standards for landslides and many scholars have studied the classification of landslides over the years.Classification of landslides according to material composition is one of the most commonly used means of classification at home and abroad, classifying landslides into rocky and soil landslides, of which soil landslides include stockpiled soil landslides, loess landslides, clayey soil landslides, and landfill landslides, etc. Varnes (1978) [20] categorized slope movements based on rock and soil movement characteristics into falls, topples, slides, lateral spreads, flows, and complex movements.Liu Guangrun (1992) [21], in his study on major geological and seismic issues in the Three Gorges of the Yangtze River, classified slope structures based on the relationship between structural planes and the slope surface into dip-slope normal, dipslope reverse, and inverse slope types.Liu Guangrun and Xu Kaixiang (1993) divided slope deformations into natural and anthropogenic dynamics based on different causative factors.Wang Lansheng, Zhang Zhuoyuan, et al. (1994) [22] proposed five basic combination models for the deformation mechanism of slopes composed of layered or layer-containing rock bodies: creep-shear, slip-induced tension cracking, bending-tension cracking, plastic flow-tension cracking, and slip-bending.Yan Tongzhen et al. (1994,2000) [23] summarized nine types of sliding mechanisms based on the initial conditions, fundamental causes, and apparent sliding modes of landslides: rheological overturning landslide, stress release translational landslide, vibration collapse or liquefaction landslide, subsidence landslide due to erosion, geochemical suspension-sinking landslide, high potential energy leapfrog landslide, pore water pressure floatation landslide, erosion-loading landslide, and giant high-speed long-runout landslide.Liu Guangrun et al. (2002) [24,25] proposed a comprehensive and systematic landslide classification system for the first time, based on the characteristics of the landslide body, deformation dynamics, and activity features, offering a detailed classification of landslides.
In recent years, the landslide classification system has continually evolved and developed, gradually forming a comprehensive classification system based on "category, size, type, mode, and stage".The type of landslide mass (material structural composition) is classified by "category", which includes rock landslides, soil landslides, and accumulation layer landslides, further subdivided based on material composition and the relationship between structural planes and slope surfaces.Rock landslides include rock layer-conformable, rock layer-nonconformable, and rock cross-layer landslides.Soil landslides include loess, clayey soil, and man-made fill landslides.Accumulation layer landslides include colluvium and talus landslides.The scale of the landslide is categorized by "size", dividing landslides into nine types based on landslide volume and landslide mass thickness, ranging from small to gigantic, and shallow to super-deep.The causative forces are classified by "type", with nine types of causative forces divided into natural and anthropogenic.The deformation movement characteristics and failure modes are sorted by "mode", classifying landslides based on kinematic features, dynamic characteristics, and deformation failure modes.The development stage of the landslide is classified by "stage", distinguishing between new and reactivated landslides in terms of their evolutionary stages.The research group based on the aforementioned geological model carries out a refined classification of accumulation layer landslides in the middle Qinling-Bashan Mountains (Figure 1).

Model Construction Method
For on-site landslide monitoring data, processing such as data partitioning is performed.According to the characteristics of time variables, data can be divided into timeseries and non-time-series data, and preprocessing is conducted based on the unique attributes of the two datasets.The landslide basic survey table records non-time-series data, including landslide scale, thickness, volume, quantity, and basic information on deformation movement characteristics.On the other hand, landslide monitoring data, such as surface displacement, rainfall, and tilt, constitute time-series data.The model needs to be capable of processing time-series data while also possessing good interpretability to assist professional technical personnel in monitoring and early warning of geological disasters, and understanding the influencing factors in landslide deformation [26].Deep learning models have difficulty interpreting the relationship between prediction results and input features.

Model Construction Method
For on-site landslide monitoring data, processing such as data partitioning is performed.According to the characteristics of time variables, data can be divided into timeseries and non-time-series data, and preprocessing is conducted based on the unique attributes of the two datasets.The landslide basic survey table records non-time-series data, including landslide scale, thickness, volume, quantity, and basic information on deformation movement characteristics.On the other hand, landslide monitoring data, such as surface displacement, rainfall, and tilt, constitute time-series data.The model needs to be capable of processing time-series data while also possessing good interpretability to assist professional technical personnel in monitoring and early warning of geological disasters, and understanding the influencing factors in landslide deformation [26].Deep learning models have difficulty interpreting the relationship between prediction results and input features.
Predictive methods based on individual slopes are closely related to historical data and struggle to reveal the characteristics of landslide movement over the complete lifecycle.Similar-type landslide prediction models consider the interpretability of the model, the landslide's primary influencing factors, and the characteristics of landslide movement throughout the lifecycle of similar landslides to construct a predictive model based on the concept of similar-type landslides.This predictive model mainly consists of two steps: the first step uses the XGBoost machine learning model for feature importance analysis, with the analyzed results serving as feature input for subsequent deformation prediction.The second step employs the LSTM deep learning model to output predictions.
The XGBoost decision tree algorithm [27] can calculate the importance of each input feature to the final prediction result.The greater the weight of feature importance, the more significant its impact on the prediction outcome.Thus, this model is selected for feature importance analysis.Additionally, XGBoost is known for its efficiency, flexibility, and robustness: it employs parallel processing to rapidly train large datasets, it can handle both numerical and categorical features, and it is capable of managing missing values.
The XGBoost objective function is shown in Equation (1): Predictive methods based on individual slopes are closely related to historical data and struggle to reveal the characteristics of landslide movement over the complete lifecycle.Similar-type landslide prediction models consider the interpretability of the model, the landslide's primary influencing factors, and the characteristics of landslide movement throughout the lifecycle of similar landslides to construct a predictive model based on the concept of similar-type landslides.This predictive model mainly consists of two steps: the first step uses the XGBoost machine learning model for feature importance analysis, with the analyzed results serving as feature input for subsequent deformation prediction.The second step employs the LSTM deep learning model to output predictions.
The XGBoost decision tree algorithm [27] can calculate the importance of each input feature to the final prediction result.The greater the weight of feature importance, the more significant its impact on the prediction outcome.Thus, this model is selected for feature importance analysis.Additionally, XGBoost is known for its efficiency, flexibility, and robustness: it employs parallel processing to rapidly train large datasets, it can handle both numerical and categorical features, and it is capable of managing missing values.
The XGBoost objective function is shown in Equation (1): where Ω(f k ) as shown in Equation (2): contains the regularization term, ŷi is the output of the model, y i is the true label, f k denotes the kth base classifier, T is the number of leaves in the tree, ω denotes the score in the corresponding leaves, and γ is the penalty term.Expanding Equation (1) using Taylor's formula yields Equation (3): The formula g i is the same as h t and Ω(f t ).The expression equations are Equations ( 4)-( 6), respectively, as follows: ) Substituting Equations ( 4)-( 6) into Equation (3) yields the leaf node weight formula shown in Equation (7): G j and H j denote the values of first-order gradient and second-order gradient at the leaf nodes, respectively; the optimal solution of the objective function is obtained by bringing Equation (7) into Equation (3) as shown in Equation ( 8): Landslide deformation prediction primarily focuses on forecasting time-series data for surface displacement [28].The Long Short-Term Memory (LSTM) network algorithm, compared to traditional Recurrent Neural Networks (RNNs), is more adept at addressing gradient vanishing problems in long sequences.It can handle longer sequence data and better capture long-term dependencies within sequential data.Its time sensitivity enables more effective learning of the features in time-series data [29].The plan is to apply a hybrid XGBoost and LSTM model to landslide deformation prediction to enhance the accuracy of predictions while also achieving model interpretability [30][31][32].
LSTM [33] introduces a memory cell that stores and accesses information and controls the flow of information through a gating mechanism.The key parts of LSTM include input gate, forget gate, and output gate (Figure 2).and  it outputs a number between 0 and 1.An output of 1 means "keep this value completely" and 0 means "throw this value away completely".
The forgetting gate is calculated as in Equation ( 9).

• Forget gate
The first step in LSTM is to decide what information we need to throw away from the cell state.This decision is made by a sigmoid layer called the forget gate.With the inputs h t−1 and x t it outputs a number between 0 and 1.An output of 1 means "keep this value completely" and 0 means "throw this value away completely".
The forgetting gate is calculated as in Equation ( 9). •

Input gate
The second step is to decide what kind of information we need to store in the cell state.There are two parts to this problem.First, the sigmoid layer calls an "input gate" to decide what data needs to be updated.Then, a tan h layer creates a vector ∼ C t as new candidate value that can be added to the state.In the next step, we want to merge these two parts to create an update to the state.
The input gate is calculated as in Equations ( 10) and ( 11): The forgetting gate and input gate are combined as in Equation ( 12): • Output gate Finally, we need to decide what we want to output.This output is based on our cell state but will be a filtered value.
The output gate is calculated as in Equations ( 13) and ( 14):

Selection of Evaluation Indicators
In order to analyze and study the prediction ability of the displacement monitoring data model, the model test uses two assessment indicators, MSE and MAPE, to measure the model prediction effect.
MSE squares the difference between the predicted and actual values for each time period, and then averages the squares of the differences over multiple time periods to obtain the average mean square error.The MSE metric has the advantage of being able to highlight extreme errors, as in Equation (15).
MAPE denotes the mean percentage of direct relative error between predicted and measured values in the range [0, +∞); the smaller it is the more accurate the prediction model is, as in Equation (16).

Research Case
The middle Qinling-Bashan Mountains is located in the Qinling Mountains and Bashan Mountain Range in southern Shaanxi, ranges from 31 • 42 ′ N to 34 • 33 ′ N and from 105 • 29 ′ E to 111 • 15 ′ E (Figure 3), and is mainly located in the cities of Hanzhong, Ankang, and Shangluo in southern Shaanxi (Figure 3).The average annual temperature in the study area is 11-14 • C and the annual precipitation is 750-1000 mm, which is mostly concentrated in July-September each year, accounting for about 50% of the annual precipitation, during which torrential rains and continuous rainy weather are predominant.Landslides in the middle Qinling-Bashan Mountains are predominantly mediumand small-scale accumulation layer landslides.The material of these landslide masses mainly consists of clay, silt, silty clay, and gravel of various sizes.The sliding surfaces are generally at the contact between bedrock and accumulation layer material, or along interfaces of bedrock with different degrees of weathering.The landslides are significantly influenced by heavy summer rainfall and human engineering activities.The main disaster evolution pattern of landslides in this central section is creep-shear, with the primary evolutionary process detailed in Table 1.These types of accumulation layer landslides have a relatively long deformation evolution process, allowing for the collection of long-term time-series displacement monitoring data, making them suitable for machine learning applications.
The research objects for this study in the middle Qinling-Bashan Mountains are accumulation layer landslides selected from automated geological hazard monitoring sites within the region.After analyzing and comparing data, 13 landslides equipped with rain gauges and surface displacement devices, and showing clear deformation trends were The geomorphology of the Qinling-Bashan Mountains area belongs to low, middle, and high mountains and wide valley basin landforms, with altitudes ranging from 170 to 3000 m above sea level, and a series of subsidence basins have been formed by the fracture and subsidence action of the landmass.The middle Qinling-Bashan Mountains area is a multiplate convergence zone; the plates squeeze each other, the ground stress is extremely complex, the rock deformation is strong, faults and folds are extremely developed, and the peak acceleration of the ground vibration is generally 0.05 g.The region is subject to intense tectonic movement, resulting in the stacking and interweaving of tectonic plates, which leads to an extremely complex distribution of rock types.The predominant lithologies in the area are Quaternary unconsolidated deposits, slate, and phyllite.Groundwater in the area primarily forms through the vertical infiltration of atmospheric precipitation.
Its movement and storage vary with depth, with shallow layers typically containing unconsolidated sediment pore water and deeper layers hosting bedrock fracture water.Human engineering activities in the area are mainly focused on housing construction, road building, and mining.The region is characterized by high mountains, deep valleys, and steep slopes, and is rife with geological structures such as faults.Additionally, concentrated rainfall and high cumulative precipitation volumes make this region highly susceptible to geological hazards.Rock types such as slate, schist, and phyllite, which are prone to slipping, as well as Quaternary slope debris accumulations and expansive soils, are widely distributed.Investigation and analysis show that geological hazards in the region exhibit interannual and monthly patterns, with more incidents occurring in wet years and fewer during dry years.The vast majority of landslides, collapses, and unstable slope failures are concentrated in the rainy season, with statistical data indicating that about 80% occur between July and September [34][35][36][37].
By the end of 2022, the region had a total of 6084 registered geologic hazard sites, posing a potential threat to the safety of 165,636 people and about CNY 1.4 billion worth of property.Among these hazards, landslides accounted for the highest proportion at 89.6%, followed by collapses at 4.9%, debris flows at 4.7%, and other types at 0.7%.
Landslides in the middle Qinling-Bashan Mountains are predominantly medium-and small-scale accumulation layer landslides.The material of these landslide masses mainly consists of clay, silt, silty clay, and gravel of various sizes.The sliding surfaces are generally at the contact between bedrock and accumulation layer material, or along interfaces of bedrock with different degrees of weathering.The landslides are significantly influenced by heavy summer rainfall and human engineering activities.The main disaster evolution pattern of landslides in this central section is creep-shear, with the primary evolutionary process detailed in Table 1.These types of accumulation layer landslides have a relatively long deformation evolution process, allowing for the collection of long-term time-series displacement monitoring data, making them suitable for machine learning applications.

State of Affairs Schematic Diagram of the Evolution of the Pattern of Disasters Characterization
Natural state

Natural state
The original slope is in its natural state, with the surface of the slope exposed, and the strongly weathered and moderately weathered surfaces of the slope are potential sliding surfaces.
The original slope is in its natural state, with the surface of the slope exposed, and the strongly weathered and moderately weathered surfaces of the slope are potential sliding surfaces.The research objects for this study in the middle Qinling-Bashan Mountains are accumulation layer landslides selected from automated geological hazard monitoring sites within the region.After analyzing and comparing data, 13 landslides equipped with rain gauges and surface displacement devices, and showing clear deformation trends were chosen to construct the dataset.These include the Balipo landslide, Chengouwan landslide, Hujiawan landslide, Koujiayuanzi landslide, Leijiapo landslide, Liujiazhuang landslide, Lvjiawan landslide, Mapoling Group 6 landslide, MeiJiayuanzi landslide, Wangjialaowuchang Group 4 landslide, Wuxuecun Group 5 landslide, Xiajiapo landslide, and West of Yinpoyuanzi landslide (Figures 4-7).The 13 landslides were classified according to the geological model established in Figure 1 (Table 2) and can be subdivided into seven types: new traction type colluvial landslide, new traction type talus landslide, new translational type colluvial landslide, reactivated translational type talus landslide, reactivated translational type broken accumulation landslide, reactivated traction type colluvial landslide, and reactivated traction type talus landslide.The classification codes are as shown in Table 3.
surfaces of the slope are potential sliding surfaces.

Experimental Process
The research team conducted machine learning experiments based on the 13 landslide accumulation layers selected from in the middle Qinling-Bashan mountains area mentioned earlier, carrying out comparative analysis work on deformation prediction for similartype landslides and mixed-type landslides.In the first group of experiments (Table 4), the training dataset for similar-type landslides consisted entirely of Type 7 landslides, which also represented the most numerous type in this dataset.The training dataset for mixed-type landslides included landslides of types other than 7.In the second group of experiments (Table 5), the training dataset for similar-type landslides was entirely composed of Type 7 (including historical monitoring data of West of Yinpo Yuanzi landslide), while the training dataset for individual slope landslides was based on historical monitoring data of the landslide on the west of the Yinpo yard.The prediction set for all experiments used the data from West of Yinpo Yuanzi landslide.

Key Feature Selection
Due to the abnormal power supply of equipment in the field, the failure of the equipment itself, or the weak communication signal causing the monitoring data to jump or be missing, the first step is to carry out the cleaning of landslide deformation monitoring data, including the filling of missing data and correction of abnormal data [38][39][40][41][42]. Due to the diversity of geohazard monitoring data, it is difficult to completely identify abnormal data with the same algorithm [43,44].In order to identify the anomalous data, the long time series of monitoring data is processed using the method of region segmentation and the specific algorithm for different morphology data is used to identify the anomalies.For data with smooth morphology, the 3σ criterion and isolated forest are used to eliminate the anomalous data.For data segments characterized by a clear trend, a least squares polynomial curve is used to fit the trend line.
Landslide deformation is a complex multidimensional nonlinear movement, which is the result of a variety of factors, and it is necessary to deeply explore the features of landslide deformation and the triggering factors affecting the deformation, such as rainfall on the same day, effective rainfall, slope, volume type, thickness type, structural type, deformation movement characteristics, and development stage.Due to the large amount of feature data, if all of them are used as input items of the model, the problem of dimensional catastrophe will occur, which will weaken the accuracy and generalization ability of the model, and lead to a decrease in the accuracy of the model.Our research group selects the random forest XGBoost machine learning algorithm in the embedding method to screen all the features (Table 6) and screens out the features in the feature data that have higher correlation with the stage of landslide displacement change, so as to get the subset of features for the input of the prediction and early warning model, which reduces the dimensionality of the input features of the landslide prediction and early warning model, and improves the model precision and generalizability.

Analysis of Results and Discussion
In order to keep the experiments comparable, the structure, hyperparameters, and amount of training data are the same for the three sets of models, which predict landslide displacements after 12 h, 24 h, and 48 h, respectively.The number of past data referenced for prediction is 18 (6 h unit), learning rate is 0.0002, and number of epochs is 300.The size of hidden size of the LSTM model is 512 and the number of training data in a single batch is 24.The optimal hyperparameters of the model are selected by using the grid search method in the experiments.The evaluation indexes of the two sets of experimental results are shown in the table, and some of the prediction effects are shown in Tables 7 and 8 and Figures 10-12   The first experiment's similar-type landslide prediction model showed optimal performance in predicting the next 12 h on the test set, with an MSE of 5.036 and a MAPE of 0.01.Under the same training conditions and predicting the same duration, the joint training of similar-type landslides yielded better results than the joint training of mixed-type landslides.
In the second experiment, the similar-type landslide prediction model exhibited the best performance in predicting the next 12 h on the test set, with an MSE of 5.099 and a MAPE of 0.01.However, for predicting the next 24 and 48 h, the individual slope landslide prediction model outperformed the similar-type landslide prediction model.

Analysis of Results and Discussion
In order to keep the experiments comparable, the structure, hyperparameters, and amount of training data are the same for the three sets of models, which predict landslide displacements after 12 h, 24 h, and 48 h, respectively.The number of past data referenced for prediction is 18 (6 h unit), learning rate is 0.0002, and number of epochs is 300.The size of hidden size of the LSTM model is 512 and the number of training data in a single batch is 24.The optimal hyperparameters of the model are selected by using the grid search method in the experiments.The evaluation indexes of the two sets of experimental results are shown in the table, and some of the prediction effects are shown in Tables 7 and 8 and Figures 10-12       The first experiment's similar-type landslide prediction model showed optimal performance in predicting the next 12 h on the test set, with an MSE of 5.036 and a MAPE of 0.01.Under the same training conditions and predicting the same duration, the joint training of similar-type landslides yielded better results than the joint training of mixed-type landslides.
In the second experiment, the similar-type landslide prediction model exhibited the best performance in predicting the next 12 h on the test set, with an MSE of 5.099 and a

Discussion
This paper presents a study on the displacement prediction of similar-type landslides in the middle Qinling-Bashan Mountains, focusing on accumulation layer landslides.The enhancement of prediction accuracy necessitates the use of prolonged and continuous monitoring data for model refinement.The methods employed for landslide classification, model construction, and prediction provide a reference for displacement prediction of landslides in other regions and of different types.The use of similar-type landslide prediction models to supplement data enables the acquisition of sliding characteristics at different evolutionary stages of landslides, offering a scientific basis for early landslide prediction and warning.
The study proposes a machine-learning-based method for constructing similar-type landslide deformation prediction models, based on the premise that similar landslides exhibit analogous movement characteristics.This approach includes key steps such as landslide classification methods, feature importance analysis and selection, construction of training sample sets, and sample learning training with optimization and model tuning.A monitoring and early warning geological model based on landslide classification is proposed, encompassing the structure, scale, inducing factors, deformation movement characteristics, failure modes, and development cycles of the landslide body.This refined classification provides a basis of expert knowledge and experience for machine learning.Furthermore, the study introduces a method combining XGBoost and LSTM models for landslide deformation prediction, enhancing the accuracy of the predictions while achieving model interpretability.The XGBoost model is used to analyze the main features of landslide deformation, which are then incorporated as input features in the training of the samples.The contribution of rainfall is approximately 55%, volume around 13%, thickness 12%, slope gradient about 11%, and material structure around 9%, with other features being negligible.
The refined classification of landslides requires manual intervention and is closely related to the experience of the researchers, leading to variability in results among different researchers, which can influence the final landslide displacement predictions.The predictive results for different landslides are closely related to the developmental stage of the landslides and are significantly correlated with the developmental stage and quality of monitoring data of similar landslides in the study area.Regions with high landslide susceptibility and widespread displacement monitoring tend to have better predictive outcomes, suggesting the potential for adoption and wider application of these methods.The predictive method proposed in this paper requires a substantial amount of long-term continuous observational data as samples for machine learning.The requirements for deformation prediction are relatively high and the applicability of these methods is somewhat limited.The predictive results are highly correlated with the duration and quality of the monitoring data, and the predictive models require continuous refinement to be truly effective.

Conclusions
This paper presents a method and process for constructing a prediction model for deformation in similar-type landslides.Taking the landslide accumulation layers in the middle Qinling-Bashan Mountains area as an example, a training sample set for these landslides was constructed.A data-driven deformation prediction model for the middle section of the Qinling-Bashan Mountains area's accumulation layers was developed, utilizing historical monitoring data and key feature parameters of similar-type landslides for predicting future 12 h, 24 h, and 48 h deformations.The model's predictive performance was also evaluated.The main conclusions of this study are as follows: 1.The displacement prediction model for similar-type landslides performs better in displacement prediction than the mixed-type landslide prediction model.

Water 2024 , 23 Figure 1 .
Figure 1.Geological model for classifying landslides in the middle Qinling-Bashan Mountains area with stacked layers.

Figure 1 .
Figure 1.Geological model for classifying landslides in the middle Qinling-Bashan Mountains area with stacked layers.

Water 2024 , 23 Figure 2 .
Figure 2. LSTM model.•ForgetgateThe first step in LSTM is to decide what information we need to throw away from the cell state.This decision is made by a sigmoid layer called the forget gate.With the inputs ℎ and  it outputs a number between 0 and 1.An output of 1 means "keep this value completely" and 0 means "throw this value away completely".The forgetting gate is calculated as in Equation (9).

Table 1 .
Creep-slip-pull-type evolutionary model of landslides in the middle portion of the accumulation layer in the Qinling-Bashan Mountains area.
at the foot of slopes and gradual decrease in stability Landslide formation Landslide formation Sliding surface penetrates, slide falls, landslide forms.Sliding surface penetrates, slide falls, landslide forms.Water 2024, 16, 464 9 of 19

Figure 10 .
Figure 10.Subgroups predicting effects over the next 12 h.(a) First set of similar slopes projected for the next 12 h.(b) First set of mixed slopes projected for the next 12 h.(c) Second set of similar slopes projections for the next 12 h.(d) Second set of single-slope projected for the next 12 h.

Figure 10 .
Figure 10.Subgroups predicting effects over the next 12 h.(a) First set of similar slopes projected for the next 12 h.(b) First set of mixed slopes projected for the next 12 h.(c) Second set of similar slopes projections for the next 12 h.(d) Second set of single-slope projected for the next 12 h.

Figure 10 .
Figure 10.Subgroups predicting effects over the next 12 h.(a) First set of similar slopes projected for the next 12 h.(b) First set of mixed slopes projected for the next 12 h.(c) Second set of similar slopes projections for the next 12 h.(d) Second set of single-slope projected for the next 12 h.

Figure 11 .
Figure 11.Subgroups predicting effects over the next 24 h.(a) First set of similar slopes projected for the next 24 h.(b) First set of mixed slopes projected for the next 24 h.(c) Second set of similar slopes projections for the next 24 h.(d) Second set of single-slope projected for the next 24 h.

Figure 11 .Figure 12 .
Figure 11.Subgroups predicting effects over the next 24 h.(a) First set of similar slopes projected for the next 24 h.(b) First set of mixed slopes projected for the next 24 h.(c) Second set of similar slopes projections for the next 24 h.(d) Second set of single-slope projected for the next 24 h.A comparison of the two sets of experiments reveals that (1) under the same amount of training data, in most cases, the shorter the prediction duration of the same model, the higher the accuracy; (2) for a 12 h prediction, the similar-type landslide prediction model performs best, followed by the individual slope landslide model, and then the mixed-type landslide model; (3) for 24 h and 48 h predictions, the individual slope landslide prediction model performs best, followed by the similar-type landslide model and then the mixed-type landslide model; (4) both the similar-type landslide prediction model and the individual slope landslide prediction model perform relatively well.Especially in terms of the MAPE metric, the prediction accuracy of the similar-type landslide model is consistent with that of the individual slope landslide model.Water 2024, 16, x FOR PEER REVIEW 19 of 23

Figure 12 .
Figure 12.Subgroups predicting effects over the next 48 h.(a) First set of similar slopes projected for the next 48 h.(b) First set of mixed slopes projected for the next 48 h.(c) Second set of similar slopes projections for the next 48 h.(d) Second set of single-slope projected for the next 48 h.

Figure 12 .
Figure 12.Subgroups predicting effects over the next 48 h.(a) First set of similar slopes projected for the next 48 h.(b) First set of mixed slopes projected for the next 48 h.(c) Second set of similar slopes projections for the next 48 h.(d) Second set of single-slope projected for the next 48 h.Analyzing the different results for predicting 12 h versus 24 and 48 h in the second group, the reason may be that the individual slope landslide has already gone through three stages: initial deformation, uniform deformation, and accelerated deformation.With a sufficient amount of training sample data, the individual slope landslide prediction model has learned the characteristics of different evolutionary stages, hence the better prediction performance, not fully demonstrating the advantage of joint training of similartype landslides.

2 .
The displacement prediction models for similar-type landslides and individual slope landslides tend to converge in performance.For predicting displacement in the next 12 h, the similar-type landslide model outperforms the individual slope model.The effectiveness of longer-duration predictions is closely related to the developmental stage of the predicted landslide.3. The similar-type landslide displacement prediction model can learn the sliding characteristics of landslides at different evolutionary stages, providing a scientific basis for early prediction and warning of landslides.It effectively addresses the issues of insufficient early monitoring data and low prediction accuracy in landslide monitoring.Author Contributions: Conceptualization, J.M., Q.Y. and M.Z.; Methodology, J.M., Q.Y. and M.Z.; Software, C.O.; Validation, Q.Y. and M.Z.; Formal analysis, Q.Y. and M.Z.; Investigation, J.M., Q.Y., M.Z. and D.M.; Resources, J.M., Q.Y. and M.Z.; Data curation, J.M., Q.Y., M.Z. and W.Z.; Writingoriginal draft, J.M., Q.Y. and M.Z.; Writing-review & editing, J.M., Q.Y., M.Z. and Y.C.; Visualization, Q.Y.; Supervision, J.M., Q.Y. and M.Z.; Project administration, J.M., Q.Y. and M.Z.; Funding acquisition, M.Z.All authors have read and agreed to the published version of the manuscript.

Table 2 .
Refined classification of landslides in the middle Qinling-Bashan Mountains region of the accumulation layer.

Table 3 .
Refined classification codes for landslides in the middle Qinling-Bashan mountains area.

Table 4 .
First group of experiments.

Table 5 .
Second group of experiments.

Table 7 .
Evaluation indicators of the results of the first group of experiments.

Table 8 .
Evaluation indicators of the results of the second group of experiments.

Table 7 .
Evaluation indicators of the results of the first group of experiments.

Table 8 .
Evaluation indicators of the results of the second group of experiments.