Code and Data for: Mechanically consistent inversion of landslide slip surface geometry by integrating time-series InSAR observations and numerical modeling

=============================================================================Dataset for: Mechanically consistent inversion of landslide slip surface geometry by integrating time-series InSAR observations and numerical modeling============================================================================= # OverviewThis repository contains the source code, numerical model geometries, InSAR deformation inputs, and validation datasets required to reproduce the results presented in the associated manuscript. To facilitate peer review and future research, the files are systematically organized into two primary case studies: the Bajiaocun landslide (a data-sparse scenario utilizing single-orbit radar) and the Xishancun landslide (a data-rich scenario utilizing dual-orbit radar for methodological comparison). # File Structure and Detailed Descriptions -----------------------------------------------------------------------------1. Bajiaocun Landslide Case (Data-Sparse Scenario & Physical Validation)-----------------------------------------------------------------------------* Bajiaocun_InSAR_Deformation_Map_ASC_Track.csv Contains the processed single-orbit (Ascending) Line-of-Sight (LOS) deformation velocity field covering the landslide area. * profile_2-2_.dxf The AutoCAD DXF file containing the 2D topographical and geological boundary coordinates of profile 2-2', which was imported into FLAC3D to construct the numerical mesh. * Bajiaocun_Dataset_Inversion_Fitting.csv Records the spatial coordinates, observed InSAR displacements, simulated model displacements, and the final fitting residuals (in mm) for the monitoring points along the profile. * Bajiaocun_Dataset_Borehole_Validation.csv Contains the independent geotechnical borehole records. It explicitly lists the inverted slip surface depths and the actual measured depths at specific borehole locations (e.g., ZK2, ZK7, etc.) used for absolute physical validation. -----------------------------------------------------------------------------2. Xishancun Landslide Case (Methodological Comparison & Dual-Orbit Integration)-----------------------------------------------------------------------------* Xishancun_Profile_AA'_InSAR_Data.csv Contains the extracted dual-orbit (Ascending and Descending) InSAR deformation velocities, along with their spatial coordinates, distributed along the continuous A-A' profile. * Xishancun_Optimized_Control_Points_AA'.csv Contains the final optimized spatial coordinates (X_Ctrl, Z_Ctrl) of the Akima spline control points, which mathematically define the inverted basal slip surface along the A-A' profile. * Xishancun_Profile_1_1'_InSAR_Data.csv Contains the extracted dual-orbit InSAR deformation velocities distributed along the 1-1' engineering geological profile. * Xishancun_Optimized_Control_Points_11'.csv Contains the final optimized spatial coordinates (X_Ctrl, Z_Ctrl) of the Akima spline control points determining the inverted slip surface along the 1-1' profile. -----------------------------------------------------------------------------3. Source Code-----------------------------------------------------------------------------* main_inversion.py The core Python script executing the coupled geomechanical inversion framework. It drives the Particle Swarm Optimization (PSO) heuristic algorithm, iteratively updates the Akima spline control points, and interacts with the finite difference numerical modeling environment to minimize kinematic residuals. # Usage NotesThe script (“main_inversion.py”) is designed to be executed within the embedded Python environment of Itasca FLAC3D (Version 7.00 or higher). The “.dxf”geometry file can be directly imported into FLAC3D for mesh generation. # Data Sources & Acknowledgments- SAR imagery provided by the European Space Agency (ESA), accessed via ASF DAAC.- DEM data provided by NASA, accessed via USGS EarthExplorer.- Bajiaocun borehole logs and geological profiles provided by the Hunan Provincial Survey and Design Institute Co., Ltd.- Xishancun comparative dual-orbit InSAR datasets and profile configurations generously provided by Shen et al. (2025). =============================================================================对应研究数据集：融合时序合成孔径雷达干涉测量（Interferometric Synthetic Aperture Radar, InSAR）观测与数值模拟的滑坡滑面几何力学一致性反演 ## 概况 本仓库包含复现相关论文手稿中结果所需的源代码、数值模型几何文件、InSAR形变输入数据与验证数据集。为便于同行评审与后续研究，文件按两个核心案例系统组织：八角村滑坡（数据稀缺场景，采用单轨道雷达数据）与西山村滑坡（数据丰富场景，采用双轨道雷达数据用于方法学对比）。 ## 文件结构与详细说明 -----------------------------------------------------------------------------1. 八角村滑坡案例（数据稀缺场景与物理验证）----------------------------------------------------------------------------- * Bajiaocun_InSAR_Deformation_Map_ASC_Track.csv：包含经处理的单轨道（升轨）视线向（Line-of-Sight, LOS）形变速度场，覆盖滑坡区域。 * profile_2-2_.dxf：AutoCAD DXF格式文件，包含2-2'剖面的二维地形与地质边界坐标，该文件被导入FLAC3D以构建数值网格。 * Bajiaocun_Dataset_Inversion_Fitting.csv：记录该剖面监测点的空间坐标、实测InSAR位移、模拟模型位移与最终拟合残差（单位：毫米）。 * Bajiaocun_Dataset_Borehole_Validation.csv：包含独立的岩土工程钻孔记录，明确列出了反演得到的滑面深度与特定钻孔位置（如ZK2、ZK7等）的实际测量深度，用于绝对物理验证。 -----------------------------------------------------------------------------2. 西山村滑坡案例（方法学对比与双轨道数据融合）----------------------------------------------------------------------------- * Xishancun_Profile_AA'_InSAR_Data.csv：包含沿连续A-A'剖面分布的双轨道（升轨与降轨）InSAR形变速度数据及其空间坐标。 * Xishancun_Optimized_Control_Points_AA'.csv：包含Akima样条控制点的最终优化空间坐标（X_Ctrl, Z_Ctrl），该控制点通过数学方式定义了A-A'剖面上反演得到的基底滑面。 * Xishancun_Profile_1_1'_InSAR_Data.csv：包含沿1-1'工程地质剖面分布的双轨道InSAR形变速度数据。 * Xishancun_Optimized_Control_Points_11'.csv：包含Akima样条控制点的最终优化空间坐标（X_Ctrl, Z_Ctrl），该控制点定义了1-1'剖面上反演得到的滑面。 -----------------------------------------------------------------------------3. 源代码----------------------------------------------------------------------------- * main_inversion.py：核心Python脚本，实现耦合岩土力学反演框架。该脚本驱动粒子群优化（Particle Swarm Optimization, PSO）启发式算法，迭代更新Akima样条控制点，并与有限差分数值模拟环境交互以最小化运动学残差。 ## 使用说明 本脚本（main_inversion.py）需在Itasca FLAC3D（7.00及更高版本）的内嵌Python环境中运行。.dxf格式几何文件可直接导入FLAC3D以生成网格。 ## 数据来源与致谢 - SAR影像由欧洲空间局（European Space Agency, ESA）提供，通过ASF DAAC获取。 - DEM数据由美国国家航空航天局（National Aeronautics and Space Administration, NASA）提供，通过USGS EarthExplorer获取。 - 八角村钻孔日志与地质剖面由湖南省勘测设计研究院有限公司提供。 - 西山村双轨道InSAR对比数据集与剖面配置由Shen等人（2025）慷慨提供。