Next Article in Journal
2D FEM Numerical Prediction of Local Seismic Effects at San Salvador Municipality (El Salvador) Induced by 2001 Earthquakes
Previous Article in Journal
Rock Magnetism and Magnetic Fabric Study of the Icelandite and Rhyodacite Long Volcanic Sequence at Mauna Kuwale, Wai’anae Volcano, Oahu, Hawaii, USA
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Twisting Theory: A New Artificial Adaptive System for Landslide Prediction

1
Semeion Research Center of Sciences of Communication, 00128 Rome, Italy
2
Department of Mathematical and Statistical Sciences, University of Colorado, Denver, CO 80204, USA
3
Department of Civil, Construction-Architectural and Environmental Engineering, University of L’Aquila, 67100 L’Aquila, Italy
4
Department of Engineering, University of Perugia, 06123 Perugia, Italy
*
Author to whom correspondence should be addressed.
Geosciences 2023, 13(4), 115; https://doi.org/10.3390/geosciences13040115
Submission received: 29 December 2022 / Revised: 22 March 2023 / Accepted: 24 March 2023 / Published: 12 April 2023

Abstract

:
Landslides pose a significant risk to human life. The Twisting Theory (TWT) and Crown Clustering Algorithm (CCA) are innovative adaptive algorithms that can determine the shape of a landslide and predict its future evolution based on the movement of position sensors located in the affected area. In the first part of this study, the TWT and CCA will be thoroughly explained from a mathematical and theoretical perspective. In the second part, these algorithms will be applied to real-life cases, the Assisi landslide (1995–2008) and the Corvara landslide (2000–2008). A correlation of 0.9997 was attained between the model estimates and the expert’s posterior measurements at both examined sites. The results of these applications reveal that the TWT can accurately identify the overall shape of the landslides and predict their progression, while the CCA identifies complex cause-and-effect relationships among the sensors and represents them in a clear, weighted graph. To apply this model to a wider area and secure regions at risk of landslides, it is important to emphasize its operational feasibility as it only requires the installation of GNSS sensors in a predetermined grid in the target area.

1. Introduction

The analysis of landslide susceptibility is complicated by the presence of environmental factors that are either non-correlated or exhibit nonlinear correlations (Huang et al., 2020 [1]). These factors make it difficult to approach the study through physically or mechanically based models. Data-driven models based on advanced learning algorithms, fuzzy logic, and Machine Learning, such as Support Vector Machines (SVM), random forest (RF), boosting methods, and especially Deep Neural Nets (DNNs), have achieved superior performance compared to statistics that, thanks to advanced modeling of complex connections of interlayer interdependencies, obtain satisfactory internal representations of the data (Conforti et al., 2014 [2]; Gomez & Kavzoglu, 2005 [3]; Bui et al., 2020 [4]; Van Dao et al., 2020 [5]; Rudin, 2019 [6]; Gunning et al., 2019 [7]; Adadi & Berrada, 2018 [8]; Wang et al., 2020 [9]; Luti et al., 2020 [10]). However, the advanced data driven approaches utilized to date (especially Deep Neural Network models, given the high number of layers and hyperparameters used) have a significant drawback—not in terms of performance, but in terms of interpretability. The extreme difficulty in understanding the precise relationship between inputs and outputs often hinders experts from comprehending the underlying dynamics, resulting in a lack of trust in the results (Rudin, 2019 [6]; Gunning et al., 2019 [7]; Adadi & Berrada, 2018 [8]). Moreover, local territory and geomorphologies are often not adequately taken into account (Xiao et al. 2020 [11]). This implies that in landslide mitigation application forecasts generated by these systems are frequently minimized or disregarded, particularly if the decision-making process has significant consequences for human lives and involves high insurance and reconstruction costs (Cui et al., 2019 [12]; Froude & Petley, 2018 [13]; Huang & Fan, [14]).
One of the most commonly used approaches in the field of Landslide Susceptibility Assessment (LSA) is to employ Convolutional Neural Networks (CNNs) to create zoning maps in order to identify areas susceptible to slope failures. CNNs are particularly effective in analyzing geospatial data and are capable of automatically learning the features of the data, without the need for explicit knowledge of the physics or geology underlying landslide formation.
In Azarafza et al., 2021 [15] a deep convolutional neural network (CNN-DNN) was developed for mapping landslide susceptibility in the Isfahan province, Iran. The model was trained and validated using historical landslide data, field records, remote sensing images, and geomorphological and environmental factors. Upon evaluation, the model was found to be the most accurate among other machine learning techniques, with a prediction accuracy of 90.9%. The resulting map showed high-susceptibility areas in the western and southwestern regions of the province, which can aid in managing landslide risk and shaping land-use planning in the Isfahan province.
Nikoobakht et al., 2022 [16] used a GIS-based approach to assess landslide susceptibility in the Gorzineh-khil region of northeastern Iran using a convolutional neural network (CNN). The 15-layer CNN was programmed in Python and was trained using data on the main triggering factors for landslides including geomorphologic/topographic and water condition parameters. The results were cross-validated using loss functions and other classifiers like SVM, k-NN, and DT, and the CNN was found to have the best accuracy with 79.0% accuracy.
Nanehkaran et al., 2021 [17] used a fuzzy, logic-based, multicriteria decision-making method with five factors affecting landslides (climate, geomorphology, tectonic and seismic parameters, geological and hydrological conditions, and human activities) to create a hazard map. Through the identification of landslide influencing factors using satellite images, DEM data, and field survey, a Convolutional Neural Network (CNN) was trained; the model was used to classify the hazard potential of the region into five susceptibility classes, showing better performance compared to other machine learning classifiers (SVM, k-NN, and DT).
Therefore, it’s crucial to continuously monitor landslide hazards. Typically, GNSS sensors are deployed at various points in the risk zone to track their movement over time with millimeter or centimeter accuracy. During observations, data on land movement is collected. In this study, utilizing data from two recorded observation efforts, we aim to demonstrate how to achieve two main objectives.
  • Reconstruct the cause-and-effect relationships that evolve over time between the monitored points. In other words, determine “which displacement” is the most probable cause of “which other displacement,” thereby identifying the points that are the center of the landslide and those that are its by-products. To achieve this objective, we will utilize a patented and published algorithm named the Crowd Clustering Algorithm (CCA) [18,19];
  • Determine the shape of the entire landslide area based on the monitored points in a risk zone. To attain this goal, we will utilize an updated version of a previously patented and published algorithm named Twisting Theory (TWT) [18,19]. After reconstructing the landslide shape at time T(n), we can predict its new shape at a subsequent time T(n + 1). To accomplish this forecasting task, we will utilize a novel Deep Neural Network (DNN) [20,21,22,23,24,25,26,27,28,29,30,31,32]. In other words, this innovative DNN has the capability to learn from the images of the complete landslide grids produced by TWT at every observation stage, and predict the location of the sensors and the shape of the entire landslide at time T(n + 1).
To accomplish these two objectives, we only require the information of the GNSS coordinates over the duration of the recorded observations. Now, let us mathematically represent this challenge.
D a t a s e t = x i t , y i t i = 1 P t = 1 T ; w h e r e : x , y =   Coordinates   of   GNSS ; P =   Number   og   GNSS   points ; T =   Number   of   observations .
Reconstructing the form/shape of the landslide is based on the information collected from the monitored points, which is a crucial yet complex task. Merely monitoring the displacement of these points in millimeters and/or centimeters over time does not provide much insight, as the Earth’s surface naturally shifts. For example, Italy moves imperceptibly towards Albania and Croatia each year. Thus, it is essential to go beyond this information to accurately determine the existence of a landslide.
The challenge in this study is to go beyond just monitoring the movement of a few GNSS sensor-monitored points, but to understand the movement of the entire terrain. The scientific challenge is to determine the behavior of all unmonitored points in relation to the few points under observation. The goal of this research is to reconstruct the overall shape of the ground movement to determine the presence, size, and evolution of a landslide.

2. Materials and Methods

This paper will describe the theory and benefits of two adaptive algorithms: Twisting Theory and Crowd Clustering Algorithm. These algorithms belong to the field of machine learning as they are based purely on data, and they evolve and learn patterns present in the data.
“Twisting Theory—TWT” is an innovative algorithm patented by Semeion (Twisting Theory (TWT): a new theory and a new class of algorithms able to model the global deformations of the space, considering the trajectories of only a little sample of points along the time flow. Applicant: Semeion Research Center & CSI. Inventor: M Buscema. USA Patent: 12/969,887. Deposited 16 December 2010) and published in international scientific journals [18,19]. The publication has shown how, starting from the data produced over time from a few points monitored by GNSS sensors, the TWT model was able to reliably reconstruct the behaviour of all “invisible” areas. It has a completely different architecture from classic neural networks: instead of reasoning on record of features, these algorithms directly process the movements of all the position sensors over time. Table 1 shows a summary comparison of advantages and disadvantages between the TWT-CCA method and the Classic LSA ANNs. The TWT-CCA solution expresses the dynamics of the landslide for the entire convex hull box in which the sensors are positioned, through a grid deformation. The widening grid shows a greater speed of the landslide; i.e., it indicates a greater susceptibility, reveals a critical point, and develops greater kinetic energy. Conversely, a narrowing grid shows the places where the landslide stops spontaneously and acquires potential energy. Therefore, the TWT–CCA algorithms are able to simultaneously display multiple susceptibility points for the entire area.
The model has already been applied to several landslides and has undergone independent validation. TWT recently was updated with a special kind of ANNs able to make predictions about the next steps of the landslide analysed.
This paper introduces the fundamental ideas behind two algorithms for landslide analysis: the Crowd Clustering Algorithm (CCA) and Twisting Theory (TWT). The main focus of the paper is to present the concepts of CCA for determining cause-and-effect relationships and TWT for reconstructing the shape of the landslide. The Deep Neural Network (DNN) is not covered in this exposition, as the emphasis is on CCA and TWT.

2.1. Crowd Clustering Algorithm (CCA)

We introduce a novel algorithm called Crowd Clustering (CC) (Patent: An algorithm to model the causation process of a discrete process. Applicant: Semeion Research Centre & CSI. Inventor: P.M. Buscema. USA Patent: 8,521,668/B2. Filed 24 March 2011), which aims to understand the mechanics of the cause-and-effect relationships between entities, in our case sensors, over time in a two-dimensional discrete process. This algorithm provides a new perspective to analyze the interactions between these entities in a dynamic environment. CC is a generalization for data georeferenced in the two-dimensional space of a previous algorithm, capable of managing multiple one-dimensional temporal signals simultaneously [18,19]. The main goal of the Crowd Clustering algorithm is to construct a plausible model that clearly shows the most probable cause-and-effect process, in which each entity at each time step affects other entities in the following time step.
It is assumed that the dataset to be processed by CC has the following:
  • Space Position of each entity (Space—“S”).
  • Amount, quantity, of displacement of each entity (Quantity—“Q”)
  • These data (space position, S, and the amount of displacement, Q) are present for each entity at any time step (Time “T”).
The Causation Model (M) that CC tries to approximate has the general form:
0.1 M = ψ S , T , Q .
Specified in greater detail, a local model (M[t, t + 1]), describing the likelihood of transition of all the entities (i, j) from time (t) to time (t + 1) has the following representation:
0.2   M i , j t , t + 1 = ψ d x i t , t + 1 , d y i t , t + 1 , d x j t , t + 1 , d y j t , t + 1 , σ m i t , t + 1 , σ m j t , t + 1 , Δ ϕ i t , t + 1 , Δ ϕ j t , t + 1 w h e r e : d x i t , t + 1 , d x j   t , t + 1   are   the   x   displacement   of   the   i - th   and   the   j - th   points   from   time   t   to   t + 1 ; d y i t , t + 1 , d y j t , t + 1   are   the   y   displacement   of   the   i - th   and   the   j - th   points   from   time   t   to   t + 1 ; σ m i t , t + 1 , σ m j t , t + 1   are   the   summation   metric   displacement   of   the   i - th   and   the   j - th     points   in   time   t   and   t + 1 ; Δ ϕ i t , t + 1 , Δ ϕ j t , t + 1   are   the   difference   of   phase   displacements   of   the   i - th   and   the   j - th     points   from   time   t   to   t + 1 .
CC needs data coming from a set of points placed in a two-dimensional space (like GNSS sensors), whose position is dynamically sampled in different temporal steps, see Table 2.
The outputs of the CC algorithm applied to any dataset of this type are the following:
  • A table of the most likely and prevalent cause-and-effect relationship among the given sensors (points) during and at the end of the temporal flow.
  • A sparse and a directed graph of the most effective excitatory relationship among the sensors during the whole process.
  • A sparse and a directed graph of the most effective inhibitory relationship among the sensors during the whole process.
The cause-and-effect relationship table plays a crucial role in comprehending the dynamics of a landslide. It highlights the areas that are key to the expansion of the landslide from a local to a larger region. Understanding the distinction between excitatory and inhibitory relationships is vital in determining the direction of the landslide’s development. Excitatory relationships occur when the movement of a sensor, and the area it represents, influences the movement of other areas. In contrast, inhibitory relationships occur when an area impedes the movement of the landslide. Furthermore, the indication of the strength of the excitatory or inhibitory relationship provides insights into the stability of the ongoing landslide phenomenon.
CC algorithm works according to the following steps:
  • The displacement of each sensor at each temporal observation is calculated (Equations (1) and (2):
    d x n = x n x n 1 ;
    d y n = y n y n 1 ;
    where n = t , t 1 .
  • The calculation of the phase and of the momentum between any couple of sequential displacement of all the sensors is represented by Equations (3) and (4):
    ϕ i n = 0 d x i n > 0 d y i n > 0 π d x i n < 0 d y i n < 0 2 π d x i n > 0 d y i n < 0 π d x i n < 0 d y i n > 0 ;
    q i n = a r c t a n d y i n d x i n + ϕ i n ;   q i n 0 , 2 π ;
  • Then the calculation of the strength of each displacement of any sensor is calculated (Equation (5)):
    m i n = 1 m * d x i n 2 + d y i n 2 2 ;   m * = M a x i , n m i n ;   m i n 0 , 1 .
  • The calculation of the resonance of the phase, the momentum, and the strength of any pair of sensors are inferred by Equations (6)–(8):
    q 0 , 2 π ;   Phase ;   m 0 , 1 ;   Module α i , j n , n + 1 = 2 π q i n q j n + 1 q i n q j n + 1 > π q j n q j n + 1 q j n q j n + 1 π ;
    Q i , j n , n + 1 = 1 2 α i , j n , n + 1 π ;
    M i , j n , n + 1 = m i n m j n + 1 ;
  • Finally, the global strength with which the displacement of each sensor acts on the displacement of any other, also considering their reciprocal geographical distance is calculated:
    S i , j = n = 1 N 1 Q i , j n , n + 1 M i , j n , n + 1 e d i , j n σ ;
At this point two other steps are computed:
f.
The calculation specified by Equations (10) and (11) is performed for each sensor in order to determine its maximum excitatory strength over time. This step enables the algorithm to identify the sensor that is most likely to induce the displacement of any other sensor during the temporal flow.
E W i n + , j = ArgMax i S i , j ;
C i , j + = 1 i = W i n + C i , j + = 0 i W i n + ;
g.
The calculation, for each sensor, of the maximum inhibitory strength it exerts on other sensors in time is provided by Equations (12) and (13). This step enables the algorithm to identify which sensor is most likely to exhibit the greatest inhibitory strength among all the sensors during the temporal progression.
I W i n , j = ArgMin i S i , j ;
C i , j = 1 i = W i n C i , j = 0 i W i n ;
At the end of the process there are two square matrices representing the cause-and-effect relationship among sensors during the temporal flow. These two matrices are:
  • A matrix of the excitatory relations, C + (i, j), among the sensors indicating which sensor imposes the direction of its displacement on any other sensor;
  • A matrix of inhibitory relations, C − (i, j), among the sensors indicating which sensor constrains the direction of the displacement on any other sensor.
These matrices can be visualized as direct and weighted graphs, which cluster the sensors and their asymmetrical relationships. These graph representations of the CC algorithm hypotheses, regarding the dynamics of the process, are easily verifiable through experimentation.

2.2. Twisting Theory (TWT)

The second algorithm is TWIST. To understand TWT, we must assume that the displacement of each GNSS in the territory over time will be the entities (x, y) that the algorithm analyzes. The movement of these entities over time represents given data points [18,19], which serve as input to the algorithm. Imagine a set of entities located in a two-dimensional space that are moving due to the influence of some underlying force. This force can be interpreted as the quantity associated with each entity, which imparts momentum to the entity proportional to its magnitude. The motion of each entity results in a trajectory in the environmental space, which is assumed to be linear for simplicity. Additionally, the space is divided into a grid and the forces acting on the entities also affect the grid, causing a deformation of the original grid. As a result, the field of forces generated by the quantities assigned to the entities can be visualized as modifications to the grid. The set of all points on the grid is referred to as geometrical points. What is the best way to describe the actual effect of the forces acting on the entities on the grid?
The Twisting Theory (TWT) algorithm is a model that predicts how each point on the grid will change its coordinates at each step in time, as entities move to their new positions. TWT divides trajectories into N equal substeps, allowing each entity to be identified by its starting coordinates, the local target during each sub-step, and its final destination coordinates (Figure 1):
TWT focuses on the distances between a particular geometrical point on the grid and the positions of the entities as they move along their trajectory. At each substep, TWT measures how the distance between the geometrical point and the entity changes according to the following.
d p j S n = x p G n x j S 2 + y p G n y j S 2 d p j T n = x p G n x j T n 2 + y p G n y j T n 2 . w h e r e x p G n , y p G n   are   the   coordinates   of   a   generic   grid   point   p   of   the   grid   at   substep   n ;   when   n = 0 ,   grid   points   all   sit   upon   the   regular   grid . x j S , y j S                               is   the   origin   Coordinates   of   each   entity   j ,   i . e . ,   assigned   point . x j T n , y j T n   are   the   local   target   of   the   coordinates   of   each   entity   j   at   any   substep   n ;   when   n = 0 ,   the   entity   lies   at   its   Origin ,   whereas   when   n = N   the   entity     has   completed   its   trajectory . d p j S n                             is   the   distance   of   a   generic   grid   point   p   from   the   Origin   of   any   entity     at   substep   n ; d p j T n       is   the   istance   of   a   generic   grid   point   p   from   the   Local   Target   n     of   any   entity   at   substep   n ;
The change in the distance between a given geometrical point and the entity’s position as it moves along its trajectory can be interpreted as the buildup of potential energy that is then able to act on the grid, altering the position of the geometrical points. The more the entity’s position changes relative to the given geometrical point from one step to the next, the more the grid is stretched or compressed in response (Figure 2).
In mathematical terms, the variation of potential energy ΔE across the whole trajectory can be expressed as:
d p , j n S = x p n G x j S 2 + y p n G y j S 2
d p , j n T = x p n G x j n T 2 + y p n G y j n T 2
Δ E p n = j = 1 N e x p d p , j n S + d , j n T α d p , j n S d p , j n T .
However, we also need to consider the side effects, in Equations (14)–(16), of the grid points among themselves. The following Equations (17)–(20) take into account all the side effects of the energy that accumulates at each grid point, absorbed from its neighbour.
R =   neighbour   of   ± Δ E p n   dimension   of   p   grid   point ;
p n + z n = a   grid   point   belongs   to   the   neighbour   R .
Δ K p n = z = 1 R e x p d p n , p + z n G + d p , p + z n 1 G α d p n , p + z n G d p n , p + z n 1 G ;
Δ F p n = Δ E p n + Δ K p n .
Thus, Equations (21)–(30) show how the potential energy is transformed into kinetic energy, that is, how it produces a deformation of the grid as determined by the movement of the source entities:
δ E x p n G = j N e x p d p , j n S + d p , j n T α x j S x j n T ;
δ K x p n G = z = 1 R e x p d p , z n G + d p , z n 1 G α x p n 1 G x p n G ;
δ E y p n G = j N e x p d p , j n S + d p , j n T α y j S y j n T ;
δ K y p n G = z = 1 R e x p d p , z n G + d p , z n 1 G α y p n 1 G y p n G ;
δ F x p n G = δ E x p n + δ K x p n ;
δ F y p n G = δ E y p n + δ K y p n ;
x p n + 1 G = x p n G + Δ F p n   δ F x p n < 0 ;  
x p n + 1 G = x p n G Δ F p n   δ F x p n 0 ;
y p n + 1 G = y p n G + Δ F p n   δ F y p n < 0 ;  
y p n + 1 G = y p n G Δ F p n   δ F y p n 0
.
TWT may be considered as a special kind of ANN whose weights are manifested by the grid deformation during its learning phase, where the dynamics implemented by the equations reaches a point of convergence. A contraction of a grid point indicates an inhibitory force, while an expansion of a grid point indicates an excitatory force. The number of iterations in this type of artificial neural network is determined beforehand by the number of observations of the movement of the entities and the number of substeps in each observation.
It is possible to analogize TWT with another type of neural network: the well-known Self Organizing Map, SOM, designed by Kohonen around the 1970s [22]. The SOM is an Unsupervised Neural Network whose termination criterion is defined a priori before to begin the training phase. This ANN generates a special kind of projection of the input space into a two-dimensional space, based on the multiple similarities of the patterns of the training set. At the end of the training phase, the resulting weights (codebook) define the optimal parameters through which the SOM classifies, in a blind way, the patterns of the testing set, taking into account the global similarity that each one of them has with the already trained codebooks. Under this point of view, TWT works in a very similar way to the SOM:
  • Both the algorithms are unsupervised.
  • The number of the epochs for TWT is also defined a priori, according to the number of observation campaigns sampled to monitor the landslide.
  • During the training phase the TWT codebooks are represented by the progressive twisting of the regular grid defined at the beginning. These codebook at the end of the training represent the weights of TWT, able to provide a picture of how the GNSSs movement has deformed the entire space of the land slide.
  • These codebook are the weights that TWT show at the end of the training phase. They are locale and explainable because the deformation of the entire plane is defined only by the local weights that deform the landslide plane in each specific part of the space. So, TWT does not work as a black box, but is more like the ANNs whose learning law is vectorial quantization (Learning Vector Quantization [23], Adaptive Vector Quantization [24]).
  • TWT is different form a SOM because during its training phase it works only on temporal series (series of observational campaigns); under this profile TWT is a special kind of recurrent ANN, because any updating of its weights considers all the previous states of the algorithm.
  • While most classic ANNs start their training with a random initialization of their weight matrices, TWT initializes its weights with the same values (regular grid—see point “c”), in the same way in which a special type of ANN, the Auto Contractive Map [6], used to do.
  • The only big difference that the TWT algorithm has with the classic ANNs is that TWT has only the training phase, because TWT modifies its weights (codebook) every time that a new observation campaign is added to the old data set. The prediction TWT was add the TWT in order to allow to the TWT to make predictions about the next step of any landslide, after the training phase. The prediction TWT will be explained in the next paragraphers, where TWT will need the support of the classic supervised ANNs to complete its predictive job.
Taking account all these considerations, we are able to assert that TWT belongs to the Artificial Adaptive System family and that, in particular, TWT is a new and unusual kind of ANN.

2.3. Predictive Twisting Theory (P-TWT)

TWT generates a new grid at each time step of its processing task. If the assigned dataset is composed of T observations along the time, TWT generates a series of grids, where each one of them will show the estimated deformation of the entire plane at that time step, according to the displacements (shifts) recorded in the assigned data set. These grids can be used as a new dataset to enable TWT to make predictions about the future shifts of the GNSS at the next steps.
To achieve this new prediction task TWT needs two new methodological operations:
  • A specifying coding of the data generated by each grid;
  • An ANN (deep or shallow) able to learn from the data already generated from the grids, in order to make estimations about the new positions of the GNSSs at the next observation. In this study, classic neural networks (both deep and shallow) as well as non-classical networks such as SVCM (deep) are utilized. Additionally, the MLP with multiple hidden layers was preinitialized with autoencoders and RBMs to prepare the weights of each layer (see [25,26]).
Further, to control the accuracy of these estimations we need to use the first grids generated by TWT as Training set, and the last grids as Testing set, in order to make a blind verification of the accuracy of the predictions of the ANN already trained.
The coding strategy: The organization of the data of each grid generated by TWT in a new data set is quite simple: a pattern will consist in the latitude and longitude of each GNSS in each grid concatenated with the coordinates of its first eight neighbors.
Let us imagine a dataset of 51 GNSS (N = 51) and 13 observations (T = 13). TWT will generate 13 grids (t1, t2, …, t13). In each grid there is the actual position (latitude and longitude) for each GNSS, and the longitude and the latitude of each of the eight first neighbors to each GNSS. Thus, each input pattern will consist of a vector of 18 components: x and y of the considered GNSS, and x and y of each of its first eight neighbors. The target to be learnt of each input pattern will be the coordinates (x and y) of that GNSS in the next TWT grid. Consequently, the topology of the ANN to be trained will present 18 inputs and 2 outputs. Figure 3 shows how the topology of each GNSS of a grid must be coded in a pattern representing that GNSS, and the target that the ANN has to learn to associate at such patterns.
If the number of GNSS of the assigned dataset is N = 51, each grid generated by TWT will be composed of 51 couples of input = 18 and output = 2 patterns.
Consequently, if the assigned dataset is composed of T = 13 observations, TWT will generate 13 grids; twelve of these will be used as training set for ANN learning process, while the last one (grid) will be used as testing set in order to estimate the ANNs prediction accuracy. This means that the training set in this exemplificative case will present 612 patterns (12 × 51), while the testing set will be composed of the 51 patterns of the thirteenth grid.
The ANN learning strategy: Once the training set and the testing set have been prepared, a specific ANN has to be selected for the learning phase. Many different ANNs topologies and Learning Laws were compared to optimize this regression task; classic Multilayer Perceptron, with different learning gradient strategies and topologies (number of hidden layers), has shown many problems of convergence (time and loss function). The best results were represented by a new supervised ANN, not based on gradient descend methodology: the Supervised Contractive Map, SVCm for short [27]. This ANN has already been successfully used in several application areas (earthquake prediction [28], atrial fibrillation detection [29]) where it has demonstrated not only superior performance compared to classical models, but also a great ability to converge in a few epochs with large and highly nonlinear data sets. We have reported the main equations of SVCm, able to suggest the differences and the similarities that this specific ANN has with the more classic and known supervised ANNs.
Legend:
  • l = number or name of the ANN layer;
  • u i l =   values   of   the   all   i t h   nodes   of   the   l t h     layer ;
  • w i , j l =   weight   matrix   connecting   the   layer   l 1   to   the   layer   l ;
  • C l =   number   of   nodes   of   l t h   layer ;
    t i =   value   of   i t h   of   the   dependent   variable ;
  • L C o e f =   ANN   learning   rate .
Signal Transfer from Input layer to Output layer (passing through each hidden layer):
C N e t i l = j C l 1 u j l 1 1 w i j l C l 1 ;
I N e t i l = j C l 1 u j l 1 w i j l ;  
u i l = sin I N e t i l 1 sin ( C N e t i l ) C l 1 .
Weights update:
δ i o u t = t i u i o u t cos I N e t i o u t 1 sin ( C N e t i o u t ) C o u t 1 ;
δ i h i d = k N u m h i d + 1 δ k h i d + 1 w k i h i d + 1 cos I N e t i h i d 1 sin ( C N e t i h i d ) C h i d 1 ;  
Δ w i j l = L C o e f δ i l u j l 1   1 w i j l C l 1 .
The supervised Contractive Map calculates two net inputs for each destination node: a classic weighted input, see Equation (32); and a contractive input, see Equation (31). The first net input tends to decay or to increase when the positive or negative value of the weight (w) becomes close to a specific constant (C).
Equation (33) activates each node according to a sine function of the two net inputs (the contractive input works as a harmonic modulation of the weighted input). The advantages and the disadvantages of sine transfer function to work properly into the topology of Multilayer Perceptron was already analyzed in the scientific literature (Le Cun et al., 1991 [30]; Le Cun et al.,1998 [31]).
Equation (34) shows a typical error calculation using the distance between the desiderate output and the estimated output, times the first derivative of sine transfer function.
Equation (35) works in the same way as Equation (4), but using the chain rule to calculate the local error of each hidden unit.
Equation (36) updates the weight matrices, using the chain rules, with a contractive factor useful to limit an extreme growth of each weight value.
In the next paragraphs (Section 3.3 and 4.2) we show two different applications of Predictive TWT (P-TWT): the Assisi and the Corvara landslide. In both cases, a Supervised Contractive Map was applied to the grid data after the coding processing, already described. In this section we show the topology implemented of SVCm for these two applications. Below are the topology and the starting parameters for the Assisi application (Table 3a), and the same information for the Corvara application (Table 3b) are shown.

3. Application 1: The Assisi Landslide

The Assisi landslide (see map in Figure 4) is described in Bovenga et al., 2013 [32]. It states: “This case refers to a landslide affecting an urban expansion area of Assisi, built starting from the years 1950–1960. The area is located completely outside the historical town centre, at a distance of some hundreds of metres from the ancient city wall. … Since 1995, the Perugia University … has established a GNSS control network over the area, integrated from 1999 on a levelling network. GNSS and levelling monitoring have been performed since then, with an approximately annual data collections up to 2008.” See Figure 5.
The final data set of the Assisi landslide is so far composed of 13 GNSS points collected during 8 observational times. In Figure 6 we have renamed the GNSS points to make clearer the following application.
The vertical displacements (1995–2008) of the points from the reference are shown in blue, while the horizontal displacements by the sensors in the landslide are shown in red.

3.1. CCA Algorithm

The CCA algorithm generates 6 direct and weighted graphs (N-2), considering 8 time steps of observation. The algorithm analyzes the delta of delta displacements of the GNSS points. The strength, displacement, and angle degree that each GNSS exercises on each of the others is made explicit in each graph. Table 4 shows this information after the first two observations, and Figure 7 displays the corresponding cause-and-effect relationship graph at this step.
We have only displayed the Excitatory Cause-and-Effect table, not the Inhibitory Cause-and-Effect table. Let us take the first row as an example. The influence of P2 on P3 is shown in the first row, with a strength of approximately 0.015, as indicated in the “Strength” column. The “Cause Module T_1” column shows the displacement of P2 from its previous delta position, which in this case is the beginning of the landslide. This value is constant for each row of P2, i.e., row 1, 2, and 3. The “Effect Module T_2” column displays the displacement of P3 due to P2’s influence in the previous delta. “Cause Degree T_1” indicates the direction of P2’s displacement in degrees, while “Effect Degree T_2” indicates the direction of P3’s displacement in degrees.
We created a topological graph from Table 4 (Figure 7) by analyzing only the causes and effects of the points. In the table, P2 acts as a force on P3, represented in the graph as a connection from P2 to P3 with an outgoing arrow. P2 also imparts a force on P4, but P4 reciprocates with a force on P2, depicted as a connection between P2 and P4 with a double-headed arrow. This results in the formation of two graphs and clusters. The first graph on the left shows that the interaction between P2 and P4 gave rise to all other cause-and-effect relationships, as indicated by the direction of the arrows. The second graph on the right demonstrates that it was the interaction between P13 and P8 that formed all other cause-and-effect relationships.
CCA cumulates the effects of the previous steps. Consequently, Table 5 and Figure 8 show respectably the effects of step 3 and step 4, taking account the steps 1 and 2.
Accordingly, Table 5 depicts the delta following the one analyzed above. Here, it can be noted that P9 exerts more force than all other points and also experiences a significant amount of force from P13, while P3 experiences minimal force compared to the others. As a result, P9 is marked with the symbol “++”, while P3 is marked with the symbol “--” representing the point of resilience.
In Figure 8, P9 is positioned in the center of the left graph and the arrows indicate that the interaction between P9 and P13 resulted in all other cause-and-effect interactions.
From step 4 to step 5 the landslide exhibits a distinct structure. Figure 9 illustrates the change in the cause-and-effect relationships graph and Table 6 provides additional details about this transformation. It is important to remember that this new configuration of the landslide is a result of the cumulative effects of steps 1, 2, 3, and 4, and the CCA outputs consider the overall dynamics of the landslide.
Step 5 and step 6 represent the final structure of this landslide. Table 7 and Figure 10 display the cumulative information gathered throughout the process and the resulting directed and weighted graph.
At this final stage, P3 is identified as the resilient point marked “--”. P1 exerts very little force on P3, while P5 experiences little force from P12 but exerts force on P1. P3 does not exert any force on any other point. P13 exerts significantly more force than the other points and is also subjected to a large amount of force from P8, therefore, marked with the symbol “++”.
The left graph of Figure 10 shows that the interaction between P8 and P13 generated all other interactions, with P3 being the furthest point from P13. The arrow connecting P12 and P5 is depicted as dashed, as the force exerted by P12 on P5 is close to zero, as previously mentioned.

3.2. The Twisting Theory Algorithm (TWT)

The measurement data is processed to quantitatively determine the kinematics of the landslide using classical computation. Then, qualitative hypotheses regarding average yearly velocities are formulated, typically using a statistical approach [33]. Applying the Twisting Algorithm to slow landslides means analyzing the measurements whilst taking into account not only the values of displacements, but also the density of monitored vertices, the mutual relationships between each measured point, and the interaction between each measured point and the geometrical space where they are located. This approach provides additional information compared to analyzing displacement values alone and improves the description of landslide phenomena, as it enables the creation of a continuous model using only a few discrete points. Additionally, TWT can be used to predict the future positions of GNSS points and the global shape of the landslide at the next time step by incorporating classic Artificial Neural Networks. However, to implement TWT, several parameters must be defined:
  • The number of observational times (T) and the number of GNSS points (P), that are inherent in the data set.
  • The number of substeps (S) in which each time step must be divided (see Figure 1 and Figure 2). The more the substeps in a simulation, the more detailed the time flow of TWT will be.
  • The number of grid cells (L) that represent the entire region where the landslide can occur. The higher the density, the more spatial details will be defined by TWT.
  • A matrix of parameters (F) representing the rigidity/elasticity of each cell of the analyzed area. In this application, it will be assumed that each grid cell has the same elasticity.
The label “Assisi (T8-P13-S10-F25-L2)” indicates that the TWT analysis of the Assisi landslide involved 8 observational time steps (T), monitored by 13 GNSS points (P). Each time step was divided into 10 substeps (S), and each grid cell is square, has a fixed elasticity of 25 (F), and the side length is equal to 2 pixels (L). The parameters T and P are determined by the data set, while variations in the parameters S, L, and F alter the temporal or spatial resolution of the simulation or the strength of deformation in the landslide. However, the overall shape of the landslide remains relatively unchanged. Figure 11 displays the shape of the Assisi landslide as defined by geologists at the end of the observational period and the positions and displacements of the 13 GNSS points, which are the only information available for the TWT analysis.
Consequently, the data set processed by TWT has the following format:
A s s i s i D a t a = x i t , y i t i = 1 P = 13 t = 1 T = 8 .
Figure 12a,b display the results of the TWT, where different grid densities are used. The larger cells indicate areas where the landslide kinematics is faster, while deformed cells indicate areas with a higher concentration of potential energy. Figure 13 uses a color code to present the same study area, with a higher intensity of red indicating a greater degree of deformation caused by the landslide.
Figure 12 detects the local deformations of the area caused by the landslide, indicating the nonuniform directions of the landslide and variations in its kinematics. Figure 13 depicts the overall shape of the landslide, which is similar to the shape of the manually drawn landslide by geologists after a thorough investigation on the site (see Figure 11a).

3.3. The Predictive TWT

The preprocessing of the TWT grids is necessary in order to prepare the data for ANNs learning and subsequent prediction (as outlined in Section 2.3). As the Assisi landslide consists of 13 GNSSs monitored over 8 observation campaigns, TWT will generate 6 grids in this case. Five of these grids will be coded to form the training set for learning with the SVCm ANN (to associate each input time (t) with the position of the GNSSs at time (t + 1)). The last grid (the sixth) will be used as the blind testing set, to predict the position of the GNSSs at step 7 and assess the accuracy of the prediction phase. Table 8 illustrates how the data of a typical grid is coded
Once the SVCm has estimated the new positions of the GNSSs at time step #7, these estimations will be added to the dataset. TWT will then be run again to generate the entire grid at step #7. Two comparisons can be made at this point: (a) a correlation between the real positions and the predicted positions of GNSSs at step #7, and (b) a comparison between the estimated grid and the grid generated by TWT using the real positions of the GNSSs at step #7. A simple linear correlation between the real grid at step #7 and the estimated grid at the same step will be used to measure the accuracy of the prediction.
Figure 14a shows the resulting state generated at the time step #7 by TWT using all available data, while Figure 14b displays the TWT state generated using the predicted positions of the 13 GNSS points at the same time step. The linear correlation between the longitude and latitude of the two maps is 0.999712 for longitude and 0.999747 for latitude, respectively, resulting in an average correlation of 0.9997295.
Figure 15 shows instead the correlation between the real positions of the GNSSs at step #7, and the estimation of the SVCm at the same step.
This application showcases the effectiveness of crowd clustering and twisting theory algorithms in modeling the dynamics of slow-moving landslides. The predicted movement of the landslide in the next time period is found to correspond closely with its actual movement.

4. Application 2: The Corvara Landslide

Corvara in Badia is a town located in the Autonomous Province of Bolzano, and it is the largest town in Val Badia, at the heart of the Dolomites (see map in Figure 4). The town is impacted by a slow-moving landslide, which is causing repeated damage to the SS 244 highway in the Campolongo Pass area, as well as posing a threat to the town of Corvara. To study this landslide, a historical series of 51 GNSS stations monitored over 26 measurement campaigns conducted from 2001 to 2008 by the Province of Bolzano was collected, with an accuracy of less than one centimeter. Figure 16 shows the processing of the 51 GNSS stations using TWT.

4.1. Corvara Landslide: TWT Application

The two objectives of TWT in this case are:
  • To reconstruct the movement of the entire landslide, including its edges, the varying sliding speeds in each region, and the potential subsidence areas;
  • To predict its evolution in the future based on the data collected from 2001 to 2005, in order to anticipate the shape of the landslide in 2008.
Figure 17, Figure 18 and Figure 19 illustrate the reconstruction of the landslide shape after 2008 by a group of geologists [34,35].
The TWT algorithm was applied to 51 points and was trained using the 26 observational monitoring collections (see Figure 16b) with the parameters T26-P51- S10-F10-L2. At the end of training phase TWT reconstructed the entire landslide pattern (Figure 20). Finally, in Figure 21, the reconstruction of the landslide perimeter is compared to the actual perimeter assessed post facto by local experts.

4.2. Corvara Landslide: TWT Prediction

The second objective of this application was to predict the Corvara landslide. A deep ANN (SVCm 24x24) was applied to the first 12 grids generated by TWT (digital grid maps). Then, the trained ANN was tested using 51 GNSS locations from the 13th to the 22nd observational data sets, to assess the accuracy of the ANN predictions.
Table 9 illustrates the construction of the training set. Table 10 displays the test results from step 13 to step 22, while Table 11 compares the results of predicting the GNSS positions in the 22nd observational data set with the actual GNSS positions at that step. The accuracy of the ANN is quite good; so, an attempt was made to reconstruct the entire pattern of the landslide at step 22 using only the predictive estimations of the ANN, as was done previously with the Assisi application. A comparison with the pattern of the TWT map trained with all available data reveals that the two patterns are highly similar.
In order to evaluate the quality of the reconstruction visually, Figure 22 presents a comparison of maps. The maps include the real map with annotations made by geologists, a reconstruction using all available data, a prediction map with the intensity of the landslide represented by colors, and a reconstruction using only the estimations of the ANN, which was trained and tested up to specific steps.

4.3. Corvara Landslide: CCA Application

Figure 23 shows the 51 GNSS monitoring the Corvara landslide on the digital map with their labels. Figure 24 and Figure 25 depict the graphs of cause-and-effect relationships among the GNSS, as identified by CCA at the end of 26 temporal intervals. Figure 24 illustrates the excitatory cause-and-effect relationships among the GNSS, while Figure 25 shows the inhibitory cause effect relations among the same GNSS.
Both in the excitatory graph (Figure 24) and in the inhibitory graph (Figure 25), the main cluster identifies GNSS#51 as the driving force of the entire landslide, as it activates most of the GNSSs around it and expands its effects through GNSS#16 to faraway GNSSs and to the border of GNSSs belonging to minor clusters. GNSS#51 also acts as the main cause of the inhibitory effects of the landslide, through GNSS#50 and GNSS#21 (See Figure 26).
In contrast to the Assisi landslide, we also reported the inhibitory cause-and-effect relationships here. Let us first analyze the graph of excitatory cause-and-effect relationships. We can see that Figure 24 is divided into three subgraphs: the largest one on the left and two smaller ones on the right. A similar division is present in Figure 25, which depicts the inhibitory cause-and-effect relationships. In both larger left subgraphs in Figure 24 and Figure 25, if we examine the direction of the arrows, we can observe that point 51 is the source of all the inhibitory and excitatory relationships in the two graphs. So, on the one hand, it has an important excitatory function in the landslide while, on the other hand, it also acts as a brake. For instance, in the inhibitory cause-and-effect relationship, point 51 is connected with point 3; however, in Figure 24, point 3 does not belong to the left subgraph, but to the right subgraph.

5. Discussion

Classical studies in various countries (see studies conducted in localized areas in Iran [15,16,17], South Korea [36], Saudi Arabia [37], Japan [38], India [39], and China [40,41,42]) often followed a typical approach: collecting expert-identified variables, gathering data from landslides and geomorphological risk situations in order to directly estimate the susceptibility levels of the examined area.
However, this procedure has some limitations and weaknesses: significant time and expertise is needed to select the significant features, preprocessing requires a lot of work, and the historical sample may not be representative enough of the phenomenon, leading to potentially incorrect classification of new geomorphological patterns. Additionally, the entire process lacks transparency in defining the cause and causal chain between the various parts of the territory and the introduced patterns, leading to a lack of clarity in defining the output.
From this perspective, the approach presented here aims to obtain specific answers about the geodynamics of locations by focusing on the singularities of the territory and their interconnections. TWT and CCA are adaptive systems that take into account the individual landslide and the positioned localization sensors, and provide a map of the local deformations that are occurring in real time throughout the territory. Neural networks are trained to predict the position of the sensors in the subsequent steps. These predictions are reintroduced into the two algorithms, which redraw the global shape of the landslide and estimate the cause-and-effect relationships between the sensors.
As a result, the identification of the landslide shape and its progression is explained locally and the various parts of the territory are localized through explicit and local cause-and-effect relationships. This supports analysts to understand, in advance, the path of the entire landslide and which parts of the territory need to be secured to halt the process. Compared to traditional systems or even deep-learning systems, TWT is proposed as a complementary method, as it aims to leverage the peculiarities of different types of predictive models capable of providing non-redundant information and perspectives. Using an ecosystem approach, multiple machine learning and evolutionary algorithms collaborate on the same database to extract a higher quantity of information, which can be advantageous for enhancing landslide susceptibility maps by offering an advanced understanding of on-site applications.
The main point we would like to highlight is that the current potential of applications of artificial intelligence tools and machine learning algorithms is not limited solely to applying a large and complex artificial neural network (ANN) to a problem with a large volume of data. To solve a practical problem, it can be more profitable to consider our machine learning algorithms as a complex ecosystem of algorithms whose cooperative architecture must be found in a new way every time. The more different algorithms we have, the more we can find an efficient solution to a new problem; in artificial intelligence, as in nature, diversity is the greatest asset of any evolution.

6. Conclusions

The presented theory is useful not only for monitoring slowly moving landslides, but it can also model the dynamics and identify the location of where the landslide is most “energetic” and its potential source. The theory, when combined with ANNs (deep and shallow [42]), it can also predict the movement of a landslide with high accuracy. The main objectives of the theory, to find cause-and-effect relationships within a landslide field over time and determine the shape of the landslide, have been accomplished. The model can be applied to large areas of national territory, with the only requirement being the placement of GNSS sensors in the areas of interest. From the comparison between the model estimates and the expert’s posterior measurements, a linear correlation of 0.9997 was obtained at both examined sites (Assisi and Corvara). Further study and additional monitoring instruments may improve the accuracy of the predicted movement of landslides.
However, a more general conclusion can be drawn about TWT. TWT is a new algorithm that expands the typology of algorithms within the Artificial Adaptive Systems family. In more detail, we propose considering TWT as a unique type of Artificial Neural Network (as discussed at the end of Section 2.2).
In this article, we aim to show an architecture in which three different adaptive algorithms (CCA, TWT, and classical ANNs) work together to detect the hidden information within a simple dataset: cause-and-effect relationships, the hidden forms of a landslide, and the prediction of the next forms of such a landslide. This is just a small example of how to approach a new problem more effectively, moving from a basic thought (what is the best algorithm to apply) to a complex architecture of thoughts (which is the most efficient artificial ecosystem capable of simulating the problem in a reasonable way). This research represents only a small step in this direction, but we believe that the methodological approach we propose is clearly outlined here: not a single machine learning algorithm to be applied everywhere, but many different learning algorithms capable of designing a specific artificial ecosystem in which the data can reveal most of its hidden information.

Author Contributions

Conceptualization, P.M.B.; Methodology, P.M.B., W.A.L., M.A.-Z., F.N., M.B., R.P, G.M. and D.D.; Software P.M.B. and G.M.; Validation, P.M.B., G.M. and D.B.; Formal Analysis, P.M.B.; Investigation, P.M.B., M.B., R.P., G.M. and D.B; Resources, D.D. and F.R.; Data Curation, P.M.B., M.B., R.P., G.M. and D.D.; Writing—Original Draft Preparation, P.M.B., W.A.L., M.A.-Z., F.N., M.B., R.P., G.M., D.B and D.D.; Writing—Review & Editing, P.M.B., M.B. and R.P.; Visualization, P.M.B., M.B., R.P., G.M. and D.B.; Supervision, P.M.B., M.B., R.P. and G.M.; Project Administration, P.M.B. and M.B.; Funding, P.M.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data are not publicly available due to legal agreements between the research institution and collaborating parties.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Huang, F.; Zhang, J.; Zhou, C.; Wang, Y.; Huang, J.; Zhu, L. A Deep Learning Algorithm Using a Fully Connected Sparse Autoencoder Neural Network for Landslide Susceptibility Prediction. Landslides 2020, 17, 217–229. [Google Scholar] [CrossRef]
  2. Conforti, M.; Pascale, S.; Robustelli, G.; Sdao, F. Evaluation of Prediction Capability of the Artificial Neural Networks for Mapping Landslide Susceptibility in the Turbolo River Catchment (Northern Calabria, Italy). Catena 2014, 113, 236–250. [Google Scholar] [CrossRef]
  3. Gomez, H.; Kavzoglu, T. Assessment of Shallow Landslide Susceptibility Using Artificial Neural Networks in Jabonosa River Basin, Venezuela. Eng. Geol. 2005, 78, 11–27. [Google Scholar] [CrossRef]
  4. Bui, D.T.; Tsangaratos, P.; Nguyen, V.-T.; Liem, N.V.; Trinh, P.T. Comparing the Prediction Performance of a Deep Learning Neural Network Model with Conventional Machine Learning Models in Landslide Susceptibility Assessment. Catena 2020, 188, 104426. [Google Scholar] [CrossRef]
  5. Dong, V.D.; Abolfazl, J.; Mahmoud, B.; Davood, M.-G.; Qi, C.; Hossein, M.; Tran, V.P.; Hai-Bang, L.; Tien-Thinh, L.; Phan, T.T.; et al. A Spatially Explicit Deep Learning Neural Network Model for the Prediction of Landslide Susceptibility. Catena 2020, 188, 104451. [Google Scholar] [CrossRef]
  6. Rudin, C. Stop Explaining Black Box Machine Learning Models for High Stakes Decisions and Use Interpretable Models Instead. Nat. Mach. Intell. 2019, 1, 206–215. [Google Scholar] [CrossRef] [Green Version]
  7. Gunning, D.; Stefik, M.; Choi, J.; Miller, T.; Stumpf, S.; Yang, G.-Z. XAI-Explainable Artificial Intelligence. Sci. Robot. 2019, 4, eaay7120. [Google Scholar] [CrossRef] [Green Version]
  8. Adadi, A.; Berrada, M. Peeking Inside the Black-Box: A Survey on Explainable Artificial Intelligence (XAI). IEEE Access 2018, 6, 52138–52160. [Google Scholar] [CrossRef]
  9. Wang, H.; Zhang, L.; Yin, K.; Luo, H.; Li, J. Landslide Identification Using Machine Learning. Geosci. Front. 2021, 12, 351–364. [Google Scholar] [CrossRef]
  10. Luti, T.; Segoni, S.; Catani, F.; Munafò, M.; Casagli, N. Integration of Remotely Sensed Soil Sealing Data in Landslide Susceptibility Mapping. Remote Sens. 2020, 12, 1486. [Google Scholar] [CrossRef]
  11. Xiao, T.; Segoni, S.; Chen, L.; Yin, K.; Casagli, N. A Step beyond Landslide Susceptibility Maps: A Simple Method to Investigate and Explain the Different Outcomes Obtained by Different Approaches. Landslides 2020, 17, 627–640. [Google Scholar] [CrossRef] [Green Version]
  12. Cui, Y.; Cheng, D.; Choi, C.E.; Jin, W.; Lei, Y.; Kargel, J.S. The Cost of Rapid and Haphazard Urbanization: Lessons Learned from the Freetown Landslide Disaster. Landslides 2019, 16, 1167–1176. [Google Scholar] [CrossRef] [Green Version]
  13. Froude, M.J.; Petley, D.N. Global Fatal Landslide Occurrence from 2004 to 2016. Nat. Hazards Earth Syst. Sci. 2018, 18, 2161–2181. [Google Scholar] [CrossRef] [Green Version]
  14. Huang, R.; Fan, X. The Landslide Story. Nat. Geosci. 2013, 6, 325–326. [Google Scholar] [CrossRef]
  15. Azarafza, M.; Azarafza, M.; Akgün, H.; Atkinson, P.M.; Derakhshani, R. Deep Learning-Based Landslide Susceptibility Mapping. Sci. Rep. 2021, 11, 24112. [Google Scholar] [CrossRef]
  16. Nikoobakht, S.; Azarafza, M.; Akgün, H.; Derakhshani, R. Landslide Susceptibility Assessment by Using Convolutional Neural Network. Appl. Sci. 2022, 12, 5992. [Google Scholar] [CrossRef]
  17. Nanehkaran, Y.; Mao, Y.; Azarafza, M.; Kockar, M.; Zhu, H.-H. Fuzzy-Based Multiple Decision Method for Landslide Susceptibility and Hazard Assessment: A Case Study of Tabriz, Iran. Geomech. Eng. 2021, 24, 407–418. [Google Scholar] [CrossRef]
  18. Massimi, V.; Asadi-Zeydabady, M.; Buscema, M.; Dominici, D.; Lodwick, W.; Simeoni, L. The Contribution of Artificial Adaptive System to Limit the Influence of Systematic Errors in the Definition of the Kinematic Behavior of an Extremely-Slow Landslide. Eng. Geol. 2016, 203, 30–44. [Google Scholar] [CrossRef]
  19. Buscema, M.; Sacco, P.L.; Grossi, E.; Lodwick, W. Spatiotemporal Mining: A Systematic Approach to Discrete Diffusion Models for Time and Space Extrapolation. In Data Mining Applications Using Artificial Adaptive Systems; Springer: New York, NY, USA, 2013; pp. 231–275. ISBN 978-1-4614-4222-6. [Google Scholar]
  20. Bengio, Y. Learning Deep Architectures for AI. Found. Trends® Mach. Learn. 2009, 2, 1–127. [Google Scholar] [CrossRef]
  21. Buscema, P.M.; Massini, G.; Breda, M.; Lodwick, W.A.; Newman, F.; Asadi-Zeydabadi, M. Artificial Adaptive Systems Using Auto Contractive Maps: Theory, Applications and Extensions, 1st ed.; Springer: Berlin, Germany, 2018; ISBN 978–3-319–75048–4. [Google Scholar]
  22. Kohonen, T. Self-Organized Formation of Topologically Correct Feature Maps. Biol. Cybern. 1982, 43, 59–69. [Google Scholar] [CrossRef]
  23. Kohonen, T. Improved versions of learning vector quantization. In Proceedings of the 1990 IJCNN International Joint Conference on Neural Networks, San Diego, CA, USA, 17–21 June 1990. [Google Scholar] [CrossRef]
  24. Buscema, M.; Catzola, L. AVQ1 Basic and AVQ2 Advanced. In Semeion Report; Semeion: Rome, Italy, 2007. [Google Scholar]
  25. Liu, W.; Wang, Z.; Liu, X.; Zeng, N.; Liu, Y.; Alsaadi, F.E. A Survey of Deep Neural Network Architectures and Their Applications. Neurocomputing 2017, 234, 11–26. [Google Scholar] [CrossRef]
  26. Buscema, M. Recirculation Neural Networks. Subst. Use Misuse 1998, 33, 383–388. [Google Scholar] [CrossRef] [PubMed]
  27. Buscema, M.; Benzi, R. Quakes Prediction Using Highly Non Linear Systems and A Minimal Dataset. In Advanced Networks, Algorithms and Modeling for Earthquake Prediction; River Publishers: Aalborg, Denmark, 2011. [Google Scholar]
  28. Buscema, M.; Massini, G.; Maurelli, G. Artificial Adaptive Systems to Predict the Magnitude of Earthquakes. Boll. Geofis. Teor. Ed Appl. 2015, 56, 227–256. [Google Scholar]
  29. Buscema, P.M.; Grossi, E.; Massini, G.; Breda, M.; Della Torre, F. Computer Aided Diagnosis for Atrial Fibrillation Based on New Artificial Adaptive Systems. Comput. Methods Programs Biomed. 2020, 191, 105401. [Google Scholar] [CrossRef]
  30. Le Cun, Y.; Kanter, I.; Solla, S.A. Eigenvalues of Covariance Matrices: Application to Neural-Network Learning. Phys. Rev. Lett. 1991, 66, 2396–2399. [Google Scholar] [CrossRef] [PubMed]
  31. Le Cun, Y.; Bottou, L.; Bengio, Y.; Haffner, P. Gradient-Based Learning Applied to Document Recognition. Proc. IEEE 1998, 86, 2278–2324. [Google Scholar] [CrossRef] [Green Version]
  32. Bovenga, F.; Nitti, D.O.; Fornaro, G.; Radicioni, F.; Stoppini, A.; Brigante, R. Using C/X-Band SAR Interferometry and GNSS Measurements for the Assisi Landslide Analysis. Int. J. Remote Sens. 2013, 34, 4083–4104. [Google Scholar] [CrossRef]
  33. Raina, R.; Madhavan, A.; Ng, A.Y. Large-Scale Deep Unsupervised Learning Using Graphics Processors. In Proceedings of the Proceedings of the 26th Annual International Conference on Machine Learning, Montreal, QC, Canada, 14 June 2009; ACM: New York, NY, USA, 2009; pp. 873–880. [Google Scholar]
  34. Raiko, T.; Valpola, H.; Lecun, Y. Deep Learning Made Easier by Linear Transformations in Perceptrons. In Proceedings of the Proceedings of the Fifteenth International Conference on Artificial Intelligence and Statistics, Auckland, New Zealand, 21 March 2012; pp. 924–932. [Google Scholar]
  35. Schmidhuber, J. Deep Learning in Neural Networks: An Overview. Neural Netw. 2015, 61, 85–117. [Google Scholar] [CrossRef] [Green Version]
  36. Sameen, M.I.; Pradhan, B.; Lee, S. Application of Convolutional Neural Networks Featuring Bayesian Optimization for Landslide Susceptibility Assessment. Catena 2020, 186, 104249. [Google Scholar] [CrossRef]
  37. Youssef, A.; Pradhan, B.; Dikshit, A.; Al-Katheeri, M.; Matar, S.; Mahdi, A. Landslide Susceptibility Mapping Using CNN-1D and 2D Deep Learning Algorithms: Comparison of Their Performance at Asir Region, KSA. Bull. Eng. Geol. Environ. 2022, 81, 165. [Google Scholar] [CrossRef]
  38. Kikuchi, T.; Sakita, K.; Nishiyama, S.; Takahashi, K. Landslide Susceptibility Mapping Using Automatically Constructed CNN Architectures with Pre-Slide Topographic DEM of Deep-Seated Catastrophic Landslides Caused by Typhoon Talas. Nat. Hazards 2023, in press. [Google Scholar] [CrossRef]
  39. Mandal, K.; Saha, S.; Mandal, S. Applying Deep Learning and Benchmark Machine Learning Algorithms for Landslide Susceptibility Modelling in Rorachu River Basin of Sikkim Himalaya, India. Geosci. Front. 2021, 12, 101203. [Google Scholar] [CrossRef]
  40. Yi, Y.; Zhang, W.; Xu, X.; Zhang, Z.; Wu, X. Evaluation of Neural Network Models for Landslide Susceptibility Assessment. Int. J. Digit. Earth 2022, 15, 934–953. [Google Scholar] [CrossRef]
  41. Chen, Y.; Ming, D.; Ling, X.; Lv, X.; Zhou, C. Landslide Susceptibility Mapping Using Feature Fusion Based CPCNN-ML in Lantau Island, Hong Kong. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2021, 14, 3625–3639. [Google Scholar] [CrossRef]
  42. Suto, J.; Oniga, S. Efficiency Investigation from Shallow to Deep Neural Network Techniques in Human Activity Recognition. Cogn. Syst. Res. 2019, 54, 37–49. [Google Scholar] [CrossRef]
Figure 1. Splitting of the entity trajectory into substeps.
Figure 1. Splitting of the entity trajectory into substeps.
Geosciences 13 00115 g001
Figure 2. The subtrajectory of the entity and the accumulation of potential energy.
Figure 2. The subtrajectory of the entity and the accumulation of potential energy.
Geosciences 13 00115 g002
Figure 3. How to code each GNSS of each TWT grid into a pattern to be learnt by an ANN.
Figure 3. How to code each GNSS of each TWT grid into a pattern to be learnt by an ANN.
Geosciences 13 00115 g003
Figure 4. Map showing the location of the two study areas within Italy: Assisi and Corvara.
Figure 4. Map showing the location of the two study areas within Italy: Assisi and Corvara.
Geosciences 13 00115 g004
Figure 5. Taken from Bovenga et al., 2013 [32]: ”(a) The Assisi landslide and its location with respect to the town center; triangles and dots represent GNSS fiducial and control points, respectively; (b) The Assisi GNSS network and the 1995–2008 displacement vectors”.
Figure 5. Taken from Bovenga et al., 2013 [32]: ”(a) The Assisi landslide and its location with respect to the town center; triangles and dots represent GNSS fiducial and control points, respectively; (b) The Assisi GNSS network and the 1995–2008 displacement vectors”.
Geosciences 13 00115 g005
Figure 6. The GNSS points renamed from P1 to P13.
Figure 6. The GNSS points renamed from P1 to P13.
Geosciences 13 00115 g006
Figure 7. Graph of cause-and-effect relationships after steps 1 and 2. The points with“++” symbol are the more active GNSS and the points with “--“ symbol are the most resilient to the displacement. Cause-and-effect relationship graph after the first two steps of the landslide shows two partially independent clusters. The cluster on the left appears to be the dominant one, with two prominent points, GNSS P7 and GNSS P3. GNSS P7 shows the highest level of activity in this phase, while GNSS P3 displays the highest level of resilience.
Figure 7. Graph of cause-and-effect relationships after steps 1 and 2. The points with“++” symbol are the more active GNSS and the points with “--“ symbol are the most resilient to the displacement. Cause-and-effect relationship graph after the first two steps of the landslide shows two partially independent clusters. The cluster on the left appears to be the dominant one, with two prominent points, GNSS P7 and GNSS P3. GNSS P7 shows the highest level of activity in this phase, while GNSS P3 displays the highest level of resilience.
Geosciences 13 00115 g007
Figure 8. Graph of cause-and-effect relationships after the steps 3 and 4. The points with“++” symbol are the more active GNSS and the points with “--“ symbol are the most resilient to the displacement. The landslide still displays two partially independent clusters, which differ slightly from the previous graph. The cluster on the left continues to be the dominant one, with two key points belonging to it. This time, GNSS (P9) is the most active, while GNSS (P3) remains the most resilient.
Figure 8. Graph of cause-and-effect relationships after the steps 3 and 4. The points with“++” symbol are the more active GNSS and the points with “--“ symbol are the most resilient to the displacement. The landslide still displays two partially independent clusters, which differ slightly from the previous graph. The cluster on the left continues to be the dominant one, with two key points belonging to it. This time, GNSS (P9) is the most active, while GNSS (P3) remains the most resilient.
Geosciences 13 00115 g008
Figure 9. Graph of cause-and-effect relationships after the steps 4 and 5. The points with“++” symbol are the more active GNSS and the points with “--“ symbol are the most resilient to the displacement. The landslide exhibits a singular change of structure, with the newly formed clusters still being unstable.
Figure 9. Graph of cause-and-effect relationships after the steps 4 and 5. The points with“++” symbol are the more active GNSS and the points with “--“ symbol are the most resilient to the displacement. The landslide exhibits a singular change of structure, with the newly formed clusters still being unstable.
Geosciences 13 00115 g009
Figure 10. Graph of cause-and-effect relationships after the steps 5 and 6. The points with“++” symbol are the more active GNSS and the points with “--“ symbol are the most resilient to the displacement. The landslide displays its final configuration regarding the cause-and-effect relations among the GNSS points.
Figure 10. Graph of cause-and-effect relationships after the steps 5 and 6. The points with“++” symbol are the more active GNSS and the points with “--“ symbol are the most resilient to the displacement. The landslide displays its final configuration regarding the cause-and-effect relations among the GNSS points.
Geosciences 13 00115 g010
Figure 11. (a) The shape of the Assisi landslide at the end of the observational period (2008); (b) The positions of the 13 GNSS points at the start of the observational period (1995), serving as the unique input for the TWT.
Figure 11. (a) The shape of the Assisi landslide at the end of the observational period (2008); (b) The positions of the 13 GNSS points at the start of the observational period (1995), serving as the unique input for the TWT.
Geosciences 13 00115 g011
Figure 12. (a) The output of TWT on the Assisi landslide (T8-P13-S10-F10-L2); (b) The output of TWT on the Assisi landslide (T8-P13-S10-F30-L8).
Figure 12. (a) The output of TWT on the Assisi landslide (T8-P13-S10-F10-L2); (b) The output of TWT on the Assisi landslide (T8-P13-S10-F30-L8).
Geosciences 13 00115 g012
Figure 13. The output of TWT on the Assisi landslide (T8-P13-S20-F20-L2) using a color code.
Figure 13. The output of TWT on the Assisi landslide (T8-P13-S20-F20-L2) using a color code.
Geosciences 13 00115 g013
Figure 14. (a) TWT with real data; (b) TWT using predicted data of the last step.
Figure 14. (a) TWT with real data; (b) TWT using predicted data of the last step.
Geosciences 13 00115 g014
Figure 15. Comparison between real latitude and longitude and predicted latitude and longitude of the 13 GNSSs at the step #7.
Figure 15. Comparison between real latitude and longitude and predicted latitude and longitude of the 13 GNSSs at the step #7.
Geosciences 13 00115 g015
Figure 16. The 51 GNSS stations monitored over 26 time units from 2001 to 2008: (a) map that displays the locations of the GNSS stations; (b) digital format that is ready for the application of the TWT.
Figure 16. The 51 GNSS stations monitored over 26 time units from 2001 to 2008: (a) map that displays the locations of the GNSS stations; (b) digital format that is ready for the application of the TWT.
Geosciences 13 00115 g016
Figure 17. Reconstruction of the Corvara landslide after 2008. The size of the arrows indicates speed, while light brown represents inactive landslides and deep brown represents active landslides.
Figure 17. Reconstruction of the Corvara landslide after 2008. The size of the arrows indicates speed, while light brown represents inactive landslides and deep brown represents active landslides.
Geosciences 13 00115 g017
Figure 18. Monitoring system for the Corvara landslide in 2016.
Figure 18. Monitoring system for the Corvara landslide in 2016.
Geosciences 13 00115 g018
Figure 19. The active Corvara landslide up to 2016.
Figure 19. The active Corvara landslide up to 2016.
Geosciences 13 00115 g019
Figure 20. TWT final pattern and (on right) the projection on the geographical map of the two main location where the landslide occurred.
Figure 20. TWT final pattern and (on right) the projection on the geographical map of the two main location where the landslide occurred.
Geosciences 13 00115 g020
Figure 21. Comparison of the final perimeter of the landslide: (a) TWT’s final pattern, (b) the real map.
Figure 21. Comparison of the final perimeter of the landslide: (a) TWT’s final pattern, (b) the real map.
Geosciences 13 00115 g021
Figure 22. (a) Real map annotated by geologists; (b) TWT reconstruction using all available data; (c) Prediction map showing the intensity of the landslide through colors; (d) TWT reconstruction using only the estimations of the ANN trained up to the 12th step and tested up to the 22nd step.
Figure 22. (a) Real map annotated by geologists; (b) TWT reconstruction using all available data; (c) Prediction map showing the intensity of the landslide through colors; (d) TWT reconstruction using only the estimations of the ANN trained up to the 12th step and tested up to the 22nd step.
Geosciences 13 00115 g022
Figure 23. Digital map of the Corvara landslide with the positions of the 51 GNSS labeled.
Figure 23. Digital map of the Corvara landslide with the positions of the 51 GNSS labeled.
Geosciences 13 00115 g023
Figure 24. Excitatory cause-and-effect relationships among the 51 GNSSs.
Figure 24. Excitatory cause-and-effect relationships among the 51 GNSSs.
Geosciences 13 00115 g024
Figure 25. Inhibitory casual relationships among the 51 GNSSs.
Figure 25. Inhibitory casual relationships among the 51 GNSSs.
Geosciences 13 00115 g025
Figure 26. Main excitatory and inhibitory cause-and-effect relationships among the 51 GNSSs.
Figure 26. Main excitatory and inhibitory cause-and-effect relationships among the 51 GNSSs.
Geosciences 13 00115 g026
Table 1. Comparison of advantages and disadvantages between the TWT–CCA method and Classic LSA ANNs.
Table 1. Comparison of advantages and disadvantages between the TWT–CCA method and Classic LSA ANNs.
FeatureTWT–CCAClassic LSA ANNs
MonitoringCan monitor slowly moving landslides.Cannot provide specific information about monitoring.
ModelingCan model the dynamics of a landslide and identify its most energetic location.Can define a global landslide risk map through macro-division into regions.
PredictionCan provide specific information on the evolution of landslides over time.Can estimate the levels of territory risk aimed at land use planning.
Achieving GoalsFinding cause-and-effect relationships and determining the shape of the landslide.Definition of a “global landslide risk map through macro-division into regions”.
ApplicationCan be applied to large areas with the requirement of GNSS sensors only.Needs specific measures regarding topographic, geologic, climatologic, seismic, human activity-related factors.
Table 2. Example of dataset to be processed by CC.
Table 2. Example of dataset to be processed by CC.
ProcessSensor 1Sensor 2Sensor …Sensor M
Lat.Long.Lat.Long.Lat.Long.Lat.Long.
t(1)x1(1)y1(1)x2(1)y2(1)xM(1)yM(1)
t(2)x1(2)y1(2)x2(2)y2(2)xM(2)yM(2)
t(N)x1(N)y1(N)x2(N)y2(N)xM(N)yM(N)
Table 3. Supervised Contractive Map: (a) Assisi Training Architecture; (b) Corvara Training Architecture.
Table 3. Supervised Contractive Map: (a) Assisi Training Architecture; (b) Corvara Training Architecture.
(a) SVCm Architecture (Assisi)(b) SVCm Architecture (Corvara)
Input: 18 Input: 18
Hidden 1: 48LCoef: 0.0100
Function: 24 Contractives
Hidden 1: 24LCoef: 0.0100
Function: 24 Contractives
Hidden 2: 48LCoef: 0.0100
Function: 24 Contractives
Hidden 2: 24LCoef: 0.0100
Function: 24 Contractives
Output: 2LCoef: 0.0100
Function: 2 Contractives
Output: 2LCoef: 0.0100
Function: 2 Contractives
Training Patterns: 612
Testing Patterns: 51
Epochs (Train): 2771
RMSE (Train): 0.0051103
Training Patterns: 65
Testing Patterns: 13
Epochs (Train): 49,572
RMSE (Train): 0.00117701
Table 4. From step 1 to step 2. The points with“++” symbol are the more active GNSS and the points with “--“ symbol are the most resilient to the displacement.
Table 4. From step 1 to step 2. The points with“++” symbol are the more active GNSS and the points with “--“ symbol are the most resilient to the displacement.
Excitations [0,1]
Cause EffectStrenghtCause Module T_1Effect Module T_2Cause Degree T_1Effect Degree T_2
P2P3--0.0152170.2308370.455262247.948959°312.468597°
P40.0430360.2308371247.948959°254.608795°
P7++0.0043370.2308370.322421247.948959°268.315308°
P3--P10.0015680.270310.13408652.326408°51.972771°
P50.0022370.270310.38500252.326408°32.795769°
P4P20.0083670.2167450.215642268.388977°258.407837°
P7++P60.0954220.4928190.869689239.630249°268.394012°
P100.0001440.4928190.42529239.630249°316.548187°
P8P130.060980.4074360.507571192.280960°207.658875°
P9P120.0002440.2094480.411627220.279877°304.003754°
P12P110.0192940.2015050.66613216.253845°239.995087°
P13P80.0951640.6812950.701071209.010361°248.589050°
P90.34110.6812950.193786209.010361°243.434967°
Table 5. GNSS (P13) becomes in this phase the most active point and GNSS (P2) works as the new resilient point of the entire landslide.
Table 5. GNSS (P13) becomes in this phase the most active point and GNSS (P2) works as the new resilient point of the entire landslide.
Excitations [0,1]
Cause EffectStrenghtCause Module T_3Effect Module T_4Cause Degree T_3Effect Degree T_4
P1P3--0.0018540.042970.069856178.106628°232.431427°
P100.000010.042970.087575178.106628°109.903755°
P120.0000050.042970.071665178.106628°123.690071°
P2P40.0122320.0509720.329537257.125000°238.073135°
P6P20.0035440.175830.470611196.898651°257.011871°
P50.0013020.175830.245087196.898651°324.192261°
P70.0188730.175830.550277196.898651°227.771057°
P7P60.0306520.1012520.499697167.035507°243.871765°
P9++P110.0063460.1570260.511014194.664459°193.495728°
P130.0690170.1570261194.664459°206.565063°
P11P1-0.0000010.1687320.05679104.620872°269.641907°
P13P80.0461290.0462230.939036222.510452°223.345932°
P9++0.0329140.0462230.510582222.510452°210.772507°
Table 6. GNSS (P9) becomes in this phase the most active point and GNSS (P3) remains the resilient point of the entire landslide.
Table 6. GNSS (P9) becomes in this phase the most active point and GNSS (P3) remains the resilient point of the entire landslide.
Excitations [0,1]
Cause EffectStrenghtCause Module T_4Effect Module T_5Cause Degree T_4Effect Degree T_5
P1P30.0008820.0430260.15963269.641907°124.707367°
P2--P40.04280.356550.510164257.011871°254.938263°
P4P2--0.223460.2496680.485026238.073135°241.161392°
P5P10.0074020.1856860.128173324.192261°315.850037°
P120.0000040.1856860.099819324.192261°307.116852°
P6P70.0912830.3785860.764968243.871765°230.276962°
P7P60.1053790.4169080.710315227.771057°236.237747°
P8P13++0.2339030.7114430.856142223.345932°221.077576°
P11P100.004480.3871610.126509193.495728°174.143997°
P12P5−0.0000010.0542950.217289123.690071°80.094704°
P13++P80.1909760.7576320.723519206.565063°220.629013°
P90.3117050.7576321206.565063°213.416580°
P110.0147450.7576320.410018206.565063°181.804001°
Table 7. GNSS (13) continues to be the most active point, and GNSS (P3) reverts to its role as the resilient point for the entire landslide.
Table 7. GNSS (13) continues to be the most active point, and GNSS (P3) reverts to its role as the resilient point for the entire landslide.
Excitations [0,1]
Cause EffectStrenghtCause Module T_5Effect Module T_6Cause Degree T_5Effect Degree T_6
P1P3--0.0016520.1281730.152114315.850037°324.061981°
P2P40.0506610.4850260.092368241.161392252.725571°
P4P20.0212110.5101640.255459254.938263348.832031°
P5P10.0108750.2172890.13327480.09470417.364714°
P6P70.1105960.7103150.12913236.237747268.329376°
P7P60.1366120.7649680.144331230.276962°242.480103°
P8P13++0.26680.7235190.221791220.629013°182.501633°
P9P100.01601910.089762213.416580°224.028992°
P120.00455510.140562213.416580°204.407669°
P12P5-0.0000010.0998190.031462307.116852°0.000000°
P13++P80.2255870.8561420.163312221.077576°193.715912°
P90.3405710.8561420.128086221.077576°187.721832°
P110.0163860.8561420.136656221.077576°154.846893°
Table 8. How all the GNSS points (yellow) must be coded with the grid points around them (columns).
Table 8. How all the GNSS points (yellow) must be coded with the grid points around them (columns).
Time StepCoord.Point 1Point 2Point 3Point 4Point 5Point 6Point 7Point 8Point 9Point 10Point 11Point 12Point 13
Time Nx0.004253176.194458239.60281434.95283920.625507175.710831278.025146394.52475511.235687684.584656733.727356814.494507464.439484
y785.632751612.993774610.464661557.697205870.910645862.164673835.528748962.201661098.5428471138.179443962.4678961094.2716061030.632324
x−34.371132141.171097199.122299−5.431697−5.462948141.161331258.783783376.433807492.280426669.680725698.658875787.339905434.339417
y792.45105615.077881615.067688557.165955879.317871850.369507821.381409967.9916991085.6483151143.585815968.048951085.6741941025.897339
x−34.384129141.164474199.124237−5.438146−5.461166141.169312258.800354376.446411492.265625669.694092698.664124787.345581434.34549
y761.693726586.110962586.105591528.199158850.359436821.416809792.437683939.0440061056.6732181114.640137939.0828251056.719604996.943298
x−5.461305170.12471229.89283823.51189423.525455170.128159287.734192405.402313523.063049698.660522727.617371816.303833463.325378
y821.398987645.840332645.877686586.109497908.306274879.337219850.331848996.9602661114.6209721172.566772997.0101931114.6374511056.693359
x0.004253176.194458239.60281434.95283920.625507175.710831278.025146394.52475511.235687684.584656733.727356814.494507464.439484
y785.632751612.993774610.464661557.697205870.910645862.164673835.528748962.201661098.5428471138.179443962.4678961094.2716061030.632324
x−5.401855170.115128229.89659123.51218423.514975170.134537287.76767405.423584523.061218698.653931727.620056816.310181463.323578
y761.671997586.100281586.151428528.18927850.375305821.420898792.444702939.0609131056.6986081114.63855939.0862431056.723267996.960876
x23.515982199.119385258.85473652.47298452.476437199.083801316.723663434.34549552.03418727.618286756.582214845.264038492.265625
y821.415955645.840942645.880127586.110352908.296997879.332581850.361023996.9432981114.6318361172.564453997.0055541114.638551056.673218
x23.537193199.122299258.85736152.48767952.46249199.076218316.743805434.35611552.0354727.611877756.579895845.268799492.268951
y792.463623615.067688615.112557.164734879.322754850.364746821.420837967.9934081085.6727291143.597168968.0477911085.6822511025.90625
x0.004253176.194458239.60281434.95283920.625507175.710831278.025146394.52475511.235687684.584656733.727356814.494507464.439484
y785.632751612.993774610.464661557.697205870.910645862.164673835.528748962.201661098.5428471138.179443962.4678961094.2716061030.632324
Time N + 1x0.006601176.193222239.61157234.94527120.634718175.710144278.024872394.517456511.233154684.593506733.717896814.501038464.426666
y785.635742612.987732610.455139557.669678870.916565862.139954835.519531962.1830441098.5378421138.171143962.4514771094.2618411030.62561
Table 9. Data for the ANN SVCm Training Set.
Table 9. Data for the ANN SVCm Training Set.
Imput Vector = 18
(x-y of each GNSS + the x-y of first 8 neighbours)
Imput Vector = 2
(x-y of the GNSS at the next step)
Number of Patterns
Step 1Step 251
Step 2Step 351
Step 3Step 451
Step 4Step 551
Step 5Step 651
Step 6Step 751
Step 7Step 851
Step 8Step 951
Step 9Step 1051
Step 10Step 1151
Step 11Step 1251
Step 12Step 1351
Number of Training Steps12
Number of Training Patterns612
Table 10. Main Cost Functions: Average Prediction Accuracy of the Trained ANN Tested on GNSS Positions from Step 13th to Step 22nd. Error is the difference between the actual position of the sensors in the next step and the predicted position of each sensor in the next step.
Table 10. Main Cost Functions: Average Prediction Accuracy of the Trained ANN Tested on GNSS Positions from Step 13th to Step 22nd. Error is the difference between the actual position of the sensors in the next step and the predicted position of each sensor in the next step.
ANNRMSEReal ErrorAbsolute ErrorLinear Corr.
D_FF_SVCm(48 × 48)(Step22)0.006309360.0015911811.365991930.99970293
Table 11. The prediction of the ANN regarding the 51 GNSS locations during the 22nd observation data collection effort.
Table 11. The prediction of the ANN regarding the 51 GNSS locations during the 22nd observation data collection effort.
Prediction at the Step 22ndXY
RealANNRealANN
GNSS_1399.974465.66418735.5558738.035
GNSS_2515.66530.31928536.5878536.059
GNSS_3620.589609.18688588.9258586.828
GNSS_4620.508618.08588510.7468510.3
GNSS_5865.304855.35448470.2258471.014
GNSS_6763.245743.90148383.5878385.679
GNSS_7919.8909922.09098261.1788261.914
GNSS_81093.0981091.9298173.6488174.809
GNSS_91238.3161252.5548202.0368201.741
GNSS_101422.4221424.2598099.3668098.508
GNSS_111171.8751175.4428023.998024.207
GNSS_121288.7731288.4947917.837922.718
GNSS_131505.481505.1597948.9887949.073
GNSS_141648.8271642.1847965.8477967.287
GNSS_151842.6551845.0738007.3248009.797
GNSS_161978.8951971.3058082.5238082.572
GNSS_172045.6612026.6377955.0347953.752
GNSS_181973.3141964.6958182.9718184.664
GNSS_191821.7371814.4138281.338283.653
GNSS_201983.4161971.5288459.5168461.534
GNSS_212366.5392362.0088087.8278087.436
GNSS_222542.8932549.1938057.9858057.367
GNSS_232576.4252577.6068210.3258210.667
GNSS_242675.2742677.1648296.9438297.728
GNSS_252932.1992950.7018269.9918270.685
GNSS_262782.0212803.3828219.5538220.042
GNSS_272911.6712922.0298370.898371.251
GNSS_282861.9022852.8338487.4888484.992
GNSS_293255.4073224.3928594.4648592.706
GNSS_303297.2953282.948443.8468443.185
GNSS_313451.0293413.7568004.1918012.435
GNSS_323461.2883421.6797886.7437885.536
GNSS_333268.8183258.6947581.0997591.052
GNSS_342698.6382697.548864.8948886.041
GNSS_352422.9332426.7738989.1948958.019
GNSS_361806.6631793.8168959.018940.632
GNSS_371779.9111778.6158711.9448715.452
GNSS_382336.9382331.2048166.5588164.779
GNSS_391297.5791304.438686.9788686.729
GNSS_401407.3761409.5388584.6398583.33
GNSS_411140.481150.4548464.1038464.272
GNSS_421133.9421139.688804.3588813.357
GNSS_43838.0848829.42698834.5978842.23
GNSS_44829.9848821.14788711.1078712.023
GNSS_45702.3591693.56928735.0028736.819
GNSS_46982.203988.35028563.4838560.64
GNSS_471997.3051984.4078345.6188348.143
GNSS_482599.6362591.5498410.3928410.652
GNSS_492532.7472521.9658450.6828452.008
GNSS_502414.1192401.2718389.0948390.367
GNSS_512324.1492302.2298215.9978218.747
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Buscema, P.M.; Lodwick, W.A.; Asadi-Zeydabadi, M.; Newman, F.; Breda, M.; Petritoli, R.; Massini, G.; Buscema, D.; Dominici, D.; Radicioni, F. Twisting Theory: A New Artificial Adaptive System for Landslide Prediction. Geosciences 2023, 13, 115. https://doi.org/10.3390/geosciences13040115

AMA Style

Buscema PM, Lodwick WA, Asadi-Zeydabadi M, Newman F, Breda M, Petritoli R, Massini G, Buscema D, Dominici D, Radicioni F. Twisting Theory: A New Artificial Adaptive System for Landslide Prediction. Geosciences. 2023; 13(4):115. https://doi.org/10.3390/geosciences13040115

Chicago/Turabian Style

Buscema, Paolo Massimo, Weldon A. Lodwick, Masoud Asadi-Zeydabadi, Francis Newman, Marco Breda, Riccardo Petritoli, Giulia Massini, David Buscema, Donatella Dominici, and Fabio Radicioni. 2023. "Twisting Theory: A New Artificial Adaptive System for Landslide Prediction" Geosciences 13, no. 4: 115. https://doi.org/10.3390/geosciences13040115

APA Style

Buscema, P. M., Lodwick, W. A., Asadi-Zeydabadi, M., Newman, F., Breda, M., Petritoli, R., Massini, G., Buscema, D., Dominici, D., & Radicioni, F. (2023). Twisting Theory: A New Artificial Adaptive System for Landslide Prediction. Geosciences, 13(4), 115. https://doi.org/10.3390/geosciences13040115

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop