Abstract：Currently, research on slow-moving slopes is hindered by a lack of integration between the temporal and spatial dimensions in analyzing movement characteristics and their underlying causes. Additionally, the potential instability of the No. 3 slope at Wangjiapo in Bailuyuan, an important residential and tourist area, poses significant risks of casualties and economic losses. To investigate the movement mechanism of this slow slope, we developed a collaborative analytical framework consisting of “temporal InSAR monitoring, multi-scale signal decomposition, three-dimensional numerical simulation, and field verification.” Initially, we extracted the annual average velocity field and time sequence displacement of the landslide mass movement by comprehensively utilizing SBAS-InSAR and StaMPS-InSAR technologies. Furthermore, we decomposed the high-precision InSAR time series displacement into periodic and trend components using wavelet decomposition. The trend component was fitted using a multi-segment nonlinear fitting method to calculate the temporal movement velocity of the landslide. We conducted Pearson correlation analysis between deep displacement measurements, GNSS data, and the wavelet-decomposed trend displacement. The relationship between rainfall and the displacement changes of the periodic components was quantified through linear regression combined with correlation analysis. Subsequently, we employed three MESH finite difference numerical models with varying degrees of freedom to assess the stability of the landslide and simulate its motion characteristics using the intensity reduction method. Finally, we performed field verification based on the velocity field derived from time-series InSAR, the wavelet-decomposed time-series displacement, and the numerical simulation results. The findings indicate that: (1) StaMPS-InSAR demonstrates higher accuracy compared to SBAS-InSAR. The InSAR velocity field reveals that the upper portion of the landslide moves more slowly, while the middle and lower sections exhibit faster velocities. (2) The Pearson correlation coefficients between the wavelet-decomposed InSAR trend displacement and the deep displacement measurements and GNSS data were 0.96 and 0.95, respectively. (3) The calculated time series velocity, following multi-segment nonlinear fitting of the trend displacement, shows spatiotemporal variations in the landslide’s movement, with fewer acceleration events in the upper part compared to the middle and lower sections. This variation in velocity changes correlates with differences in soil strength at different locations. Notably, there is no significant linear correlation between the periodic displacement changes and rainfall. The safety factor obtained from numerical simulations was 0.967, indicating that the lower part of the slope is a stress concentration area, consistent with findings from InSAR and on-site investigations. This analytical framework has achieved a multi-dimensional spatiotemporal analysis of landslide movement through the organic coupling and progressive verification of multiple technical methods, providing a scientific basis for addressing the Wangjiapo No. 3 landslide and offering valuable insights for studying the movement mechanisms of similar slow-moving landslides.
王洪明，师 芸. 联合时序InSAR技术与数值模拟的王家坡三号滑坡体慢速运动机制解析[J]. 岩石力学与工程学报, 2025, 44(12): 3367-3385.
WANG Hongming，SHI Yun. Analysis of slow motion mechanism of Wangjiapo No. 3 landslide based on InSAR and numerical simulation. , 2025, 44(12): 3367-3385.
[1] XU Y，LU Z，LESHCHINSKY B. Kinematics of irrigation-induced landslides in a Washington Desert：impacts of basal geometry[J]. Journal of Geophysical Research：Earth Surface，2022，127(2)：006355.
[2] CRUDEN D M，VARNES D J. landslide types and processes，transportation research board，us national academy of sciences，special report[J]. Transp Res Board，1996，247(1)：36–57.
[3] LACROIX P，HANDWERGER A L，BIEVRE G. Life and death of slow-moving landslides[J]. Nature Reviews Earth and Environment，2020，1(8)：404–419.
[4] MANSOUR M F，MORGENSTERN N R，MARTIN C D. Expected damage from displacement of slow-moving slides[J]. Landslides，2011，8(1)：117–131.
[5] ALONSO E，ZERVOS A，PINYOL N M. Thermo-poro-mechanical analysis of landslides：from creeping behaviour to catastrophic failure[J]. Géotechnique，2016，66(3)：202–219.
[6] INTRIERI E，RASPINI F，FUMAGALL A，et al. The Maoxian landslide as seen from space：detecting precursors of failure with Sentinel-1 data[J]. Landslides，2018，15(1)：123–133.
[7] HANDWERGER A L，FIELDING E J，HUANG M H，et al. Widespread initiation，reactivation，and acceleration of landslides in the northern California Coast Ranges due to extreme rainfall[J]. Journal of Geophysical Research：Earth Surface，2019，124(7)： 1 782–1 797.
[8] FERRETTI A，PRATI C，ROCCA F. Permanent scatterers in SAR interferometry[J]. IEEE Transactions on geoscience and remote sensing，2002，39(1)：8–20.
[9] BERARDINO P，FORNARO G，LANARI R，et al. A new algorithm for surface deformation monitoring based on small baseline differential SAR interferograms[J]. IEEE Transactions on Geoscience and Remote Sensing，2002，40(11)：2 375–2 383.
[10] ZHANG L，LU Z，DING X，et al. Mapping ground surface deformation using temporarily coherent point SAR interferometry：Application to Los Angeles Basin[J]. Remote Sensing of Environment，2012，117(1)：429–439.
[11] HOPPER A. A multi-temporal InSAR method incorporating both persistent scatterer and small baseline approaches[J]. Geophysical research letters，2008，35(16)：302–307
[12] FERRETTI A，FUMAGALLI A，NOVALI F，et al. A new algorithm for processing interferometric data-stacks：SqueeSAR[J]. IEEE transactions on geoscience and remote sensing，2011，49(9)：3 460–3 470.
[13] BERTI M，CORSINI A，FRANCESCHINI S，et al. Automated classification of persistent scatterers interferometry time series[J]. Natural Hazards and Earth System Sciences，2013，13(8)：1 945– 1 958.
[14] LIU Y，QIU H，YANG D，et al. Deformation responses of landslides to seasonal rainfall based on InSAR and wavelet analysis[J]. Landslides，2022，19(1)：1–12.
[15] URGILEZ V A，HANDWERGER A L，BAKKER M，et al. A new method to detect changes in displacement rates of slow-moving landslides using InSAR time series[J]. Landslides，2022，19(9)： 2 233–2 247.
[16] 罗军尧，朱国金，冯业林，等. 基于SBAS-In SAR技术的澜沧江哑贡倾倒体变形演化特征[J]. 岩石力学与工程学报，2024，43(12)：3 018–3 031.(LUO Junyao，ZHU Guojin，FENG Yelin，et al. Deformation evolution characteristics of Lantang River dummy tribute dump based on SBAS-In SAR technique[J]. Chinese Journal of Rock Mechanics and Engineering，2019，43(12)：3 018–3 031.(in Chinese))
[17] ZHU Y，YAO X，YAO L，et al. Identifying the mechanism of toppling deformation by InSAR：A case study in Xiluodu reservoir，Jinsha river[J]. Landslides，2022，19(10)：2 311–2 327.
[18] 赵超英，刘晓杰，张 勤，等. 甘肃黑方台黄土滑坡InSAR识别、监测与失稳模式研究[J]. 武汉大学学报：信息科学版，2019，44(7)：996–1 007.(ZHAO Chaoying，LIU Xiaojie，ZHANG Qin，et al. Study on InSAR recognition，monitoring and instability model of Heifangtai Loess landslide in Gansu province[J]. Journal of Wuhan University：Information Science，2019，44(7)：996–1007(in Chinese))
[19] DHATT G，LEFRANCOIS E，TOUZOT G. Finite element method[M]. [S. l]：John Wiley and Sons，2012：256–259.
[20] MOCZO P，KRISTEK J，HALADA L. The finite-difference method for seismologists[J]. An introduction，2004，161(1)：3–12.
[21] PANDA S D，KUMAR S，PRADHAN S P，et al. Effect of groundwater table fluctuation on slope instability：a comprehensive 3D simulation approach for Kotropi landslide，India[J]. Landslides，2023，20(3)：663–682.
[22] MA P，CUI Y，WANG W，et al. Coupling InSAR and numerical modeling for characterizing landslide movements under complex loads in urbanized hillslopes[J]. Landslides，2021，18(5)：1 611–1 623.
[23] 凌 晴，张 勤，张 静，等. 融合工程地质资料与GNSS高精度监测信息的黑方台党川黄土滑坡稳定性研究[J]. 测绘学报，2022，51(10)：2 226–2 238.(LING Qing，ZHANG Qin，ZHANG Jing，et al. Study on the stability of loess landslide in Taidangchuan，Heilongjiang Province by integrating engineering geological data and GNSS high-precision monitoring information[J]. Journal of Geodesy and Cartography，2019，51(10)：2 226–2 238.(in Chinese))
[24] FERNANDEZ-MERODO J A，GARCIA-DAVALILLO J C，HERRERA G，et al. 2D viscoplastic finite element modelling of slow landslides：the Portalet case study (Spain)[J]. Landslides，2014，11(1)：29–42.
[25] BUR G，FERNANDEZ-MERODO J A，GARCIA-DAVALILLO J C，et al. Site scale modeling of slow-moving landslides，a 3D viscoplastic finite element modeling approach[J]. Landslides，2018，15(2)：257–272.
[26] HOPPER A，ZEBKER H A. Phase unwrapping in three dimensions with application to InSAR time series[J]. Journal of the Optical Society of America A，2007，24(9)：2 737–2 747.
[27] ZHOU C，CAO Y，GAN L，et al. A novel framework for landslide displacement prediction using MT-InSAR and machine learning techniques[J]. Engineering Geology，2024，334(2)：107497.
[28] LI H，XU Q，HE Y，et al. Modeling and predicting reservoir landslide displacement with deep belief network and EWMA control charts：a case study in Three Gorges Reservoir[J]. Landslides，2020，17(3)：693–707.
[29] PAEK J M，LEE B C，CHAE S W，et al. Surface reconstruction from FE mesh model[J]. Journal of Computational Design and Engineering，2019，6(2)：197–208.
[30] ZIENKIEWUCZ O C，TAYLOR R L，ZHU J Z. The finite element method：its basis and fundamentals[M]. [S. l.]：Elsevier，2005：361–342.
[31] SANSOSTI E，BERARDINO P，MANUNTA M，et al. Geometrical SAR image registration[J]. IEEE Transactions on Geoscience and Remote Sensing，2006，44(10)：2 861–2 870.
[32] HANSSEN R F. Radar interferometry：data interpretation and error analysis[M]. [S. l.]：Springer Science and Business Media，2001：85–89.
[33] QIN Y X，PERISSIN D，BAI J. Investigations on the coregistration of Sentinel-1 TOPS with the conventional cross-correlation technique[J]. Remote Sensing，2018，10(9)：1 405.
[34] SEYMOUR M S，CUMMING I G. Maximum likelihood estimation for SAR interferometry[C]// IEEE International Geoscience and Remote Sensing Symposium. [S. l.]：[s. n.]，1994：2 272–2 275.
[35] TONG X P，SCHMIDT D. Active movement of the Cascade landslide complex in Washington from a coherence-based InSAR time series method[J]. Remote Sensing of Environment，2016，186(1)：405–415.
[36] HOOPER A，BEKAERT D，SPAANS K，et al. Recent advances in SAR interferometry time series analysis for measuring crustal deformation[J]. Tectonophysics，2012，514(5)：1–13.
[37] NOTTI D，HERRERA G，BIANCHINI S，et al. A methodology for improving landslide PSI data analysis[J]. International Journal of Remote Sensing，2014，35(6)：2 186–2 214.
2025, Vol. 44 
