Unsupervised Machine Learning for Dynamic Slope Stability Classification: A Comparative Evaluation of PCA-K-Means, SOM, and Hybrid Algorithms Using InSAR Time-Series Data

Interpreting complex, non-linear Interferometric Synthetic Aperture Radar (InSAR) displacement time-series data for infrastructure risk assessment remains a significant geotechnical challenge. This is particularly evident in regions with established road and railway infrastructures, where the primary objective is monitoring the entire network to ensure safety and operational continuity. Because landslide displacement is a highly complex process affected by a combination of internal geological conditions and external triggers, time-series data inherently encode non-linear trends and periodic fluctuations. To address this, a data-driven framework utilizing a sliding-window transformation to engineer temporal-kinematic features is proposed, providing a broader framework for the contextualization of slope stability assessment from a network perspective. This is paired with Principal Component Analysis (PCA) for dimensionality reduction and evaluated across four unsupervised architectures: K-means, Self-Organising Maps (SOMs), Hybrid SOM-K-means, and PCA-K-means. The comparative evaluation reveals that the PCA-K-means pipeline performed best, offering a highly efficient and scalable workflow. The analysis revealed that the optimized PCA-K-means architecture successfully captured 79.20% of the kinematic variance across the first two principal components. Furthermore, it achieved a robust Between-Cluster-to-Total-Sum-of-Squares (BCSS/TSS) ratio of 71.70%, an optimal Silhouette Score of 0.320, and a low Quantisation Error (QE) of 0.90, demonstrating superior spatial separation and geometric accuracy compared to traditional heuristic methods. When cross-validated against static topographic susceptibility models, the dynamic kinematic clusters exhibited a 23% spatial convergence at the polar bounds of risk, successfully grounding the algorithm’s predictions in physically verified geomorphological features. Relying on the statistical volatility of displacements, this optimal model successfully partitioned the data into five distinct geotechnical risk classes, ranging from stable (Class A) to extreme risk (Class E). The results demonstrate that the developed dynamic framework provides a highly reliable, actionable tool for proactive, large-scale slope stability and infrastructure risk assessment.