《Remote Sensing》：Long-Horizon Mining Subsidence Forecasting and Ecological Time-Lag Assessment Using Multi-Source Remote Sensing
编辑推荐：
地下采煤诱发的地表形变具有强非线性和空间异质性特征，这给预警与生态评估带来挑战。尽管基于干涉合成孔径雷达(InSAR)的数据同化模型与基于遥感的生态指标已广泛用于长期监测，仍存在三个未解决的根本局限：(1)形变样本空间分布严重失衡使数据驱动模型偏向均值回归预测
地下采煤诱发的地表形变具有强非线性和空间异质性特征，这给预警与生态评估带来挑战。尽管基于干涉合成孔径雷达(InSAR)的数据同化模型与基于遥感的生态指标已广泛用于长期监测，仍存在三个未解决的根本局限：(1)形变样本空间分布严重失衡使数据驱动模型偏向均值回归预测；(2)递归多步预测累积误差，导致长时序外推不稳定；(3)生态监测中植被韧性引发多年观测滞后，造成光学指标“伪稳定”偏差。为应对上述问题，研究人员提出一种融合多步形变预测与生态时滞分析的统一框架。以大同煤田为研究区，利用2017年3月至2024年12月231景Sentinel-1A影像进行小基线集干涉合成孔径雷达(SBAS-InSAR)形变反演；采用空间分层采样策略提取5894个代表点，重构24步回看与15步前向窗口，系统比较六种预测模型；同时基于Landsat数据反演遥感生态指数(RSEI)进行交叉滞后分析。结果表明：(1)最大形变速率达?276.75 mm/年，累积沉陷超?2000 mm；(2)3步短期预测各模型均稳健，长短期记忆网络(LSTM)最优(RMSE=5.78 mm)，15步极端外推时传统递归模型显著发散(卡尔曼滤波RMSE=45.70 mm)，而N-BEATS借助多输入多输出(MIMO)策略抑制时序误差级联，RMSE=17.98 mm；(3)塌陷核心区生态退化同步(Lag 0)，边缘盆地存在1–2年隐性退化滞后期。该研究为复杂矿区长效主动安全管理和扰动生态系统全周期跟踪修复提供了可靠科学支撑。
研究背景方面，高强度地下采煤引发覆岩运动与地表非线性沉陷，威胁矿山设施安全并重塑水文植被生境。现有InSAR与遥感生态监测面临三方面的瓶颈：一是形变场被大量稳定背景点主导，全局模型优化易陷入“均值回归”而低估核心塌陷趋势；二是主流单步递归多步外推依赖中间态反馈，长时序误差级联放大；三是植被干旱韧性与土壤水分渐进耗竭使光学生态指标出现“伪稳定”时滞偏差。为此，研究人员依托大同煤田（山西朔北，黄土覆盖、双系叠采煤层、半干旱大陆季风）开展研究，旨在突破空间采样偏倚与时序误差传播双重限制，并量化形变—生态响应的像素级时滞机制。论文发表于《Remote Sensing》。
关键技术方法上，研究人员使用231景Sentinel-1A SLC影像（2017.3–2024.12，IW/VV，引入POE精密轨道与GACOS大气延迟校正，以SRTM 30 m DEM辅助，GNSS SXDT站验证）实施SBAS-InSAR反演；形变序列经排除>?30 mm/年稳定点后，按形变梯度施以空间分层异质网格化抽样（网格0.005°至0.0003°），得5894个代表点时序；以24步回看、15步盲测窗口，对比ARIMA、Kalman、RF、XGBoost、LSTM（单步递归）与N-BEATS（MIMO）六模型；生态端选取Landsat 8/9生长季无云影像，由NDVI、湿度(WET)、干度(NDBSI)、地表温度(LST)经主成分分析(PCA)构建RSEI，并在像素尺度以Pearson相关提取Lag∈{0,1,2}年的形变—生态交叉滞后。
研究结果部分，首先“Long-Term Surface Deformation Monitoring Results in the Datong Mining Area”表明大同煤田2017–2024年形变速率场空间极不平衡，外围?2000 mm，GNSS验证STD 5.58 mm、RMSE 8.79 mm。“Results of Multi-Model Multi-Step Forward Subsidence Forecasting”显示3步预测全模型RMSE<7 mm（LSTM 5.78 mm最优），15步时ARIMA/Kalman分别发散至26.55/45.70 mm，N-BEATS以MIMO并行输出将RMSE压至17.98 mm，DM检验p<0.01确认统计显著。“Evolution of Ecological Quality and Cross-Lagged Spatial Patterns in Mining Areas”指出RSEI低值区与沉陷漏斗重合；交叉滞后图中95.17%为背景无响应，显著耦合区占4.83%，呈“中心—边缘”格局：核心Lag 0（1.37%）同步退化，边缘Lag 1（1.60%）与Lag 2（1.86%）环带状滞后响应。
讨论部分归纳：InSAR精度受单站GNSS限制待补水准数据；核心区突变破裂越植被恢复阈值致Lag 0，边缘缓沉叠加地下水降深与根系水分胁迫、半干旱植被生理耐旱形成1–2年隐性退化；N-BEATS的MIMO切断递归反馈链是其长时序鲁棒根因；框架需在多矿区推广验证。
结论翻译浓缩如下：(1)地面形变高度异质，外围?2000 mm。(2)六模型对比证伪递归累积误差可控性——3步内皆RMSE<7 mm；15步时Kalman达45.70 mm，N-BEATS借MIMO抑时序误差级联稳于17.98 mm。(3)交叉滞后显核心—外围格局，中心Lag 0（1.37%）即时退化，外围Lag 1（1.60%）与Lag 2（1.86%）延迟响应，由植被韧性与土壤水渐进耗竭平衡控制。未来研究将框架扩展至多矿区验证泛化性。
