多源动态分析揭示喀斯特景观人类活动强度与生态系统韧性的尺度依赖、非线性阈值及季节反转耦合

《Environmental and Sustainability Indicators》：Linking Human Activities and Ecosystem Resilience in Karst Landscapes: Nonlinear Responses Revealed by Multi-source Dynamic Analysis

　　

喀斯特生态系统对人类与气候扰动高度敏感。空间静态指标难以刻画景观韧性的多尺度时空动态。为环境管理提供有效的可持续性指标，该研究量化云南喀斯特地区生态系统韧性（Ecosystem Resilience, ER）与人类活动强度（Human Activity Intensity, HAI）之间尺度依赖的非线性耦合。研究人员利用2014—2024年多源数据，构建抵抗-适应-恢复（Resistance-Adaptation-Recovery, RAR）框架作为生态可持续性的动态指标。采用池化广义相加模型（Pooled Generalized Additive Model, Pooled GAM）将地形约束与HAI解耦，得到净化人类活动强度（Purified Human Activity Intensity, pHAI），以准确测度人为压力。研究人员利用Spearman相关、双变量Moran's I本地空间关联指标（Bivariate Moran's I Local Indicators of Spatial Association, LISA）、极端梯度提升（eXtreme Gradient Boosting, XGBoost）、沙普利加性解释（SHapley Additive exPlanations, SHAP）与增量Moran's I识别尺度交互与非线性驱动因子。结果表明：（1）ER与pHAI呈异步时空轨迹。（2）Spearman相关与LISA识别10级子流域为最佳交互尺度。（3）pHAI影响贡献从约16%（旱季）升至约28%（湿季）。（4）出现季节反转：旱季ER仅与pHAI强度弱相关，但对累计趋势/总量敏感，而湿季ER易受强度峰值影响。（5）ER基线与跨季韧性关联最强，超过约0.54（旱季）和约0.51（湿季）阈值后出现自维持状态。这些发现建立了稳健的指标体系，用于评估人类压力增加下的系统韧性，为脆弱山地环境动态、基于节奏的空间治理提供定量依据。

研究背景与问题：喀斯特生态系统对人类和气候扰动高度敏感，但空间静态指标难以捕捉景观韧性的多尺度时空动态。人类活动从外部扰动转向内生驱动，而ER作为扰动下维持关键结构与功能的能力，其空间显式、尺度依赖的HAI响应仍难以量化。既有HAI-ER耦合分析多依赖线性或空间平稳回归，难以识别非线性阈值、抵抗饱和与尺度依赖响应异质性；喀斯特高度破碎地形又会使地形约束被误归因于人类影响，产生统计内生性。云南喀斯特连续碳酸盐岩覆盖高、干湿季分明，为检验去地形约束后的真实非线性机制提供实验场。该研究构建RAR动态框架、Pooled GAM去地形约束获得pHAI，并结合多尺度相关、双变量Moran's I LISA、XGBoost-SHAP与增量Moran's I，揭示pHAI-ER非线性驱动与季节响应异质性，为脆弱山地适应性管理提供依据。

关键技术方法概括：数据主要覆盖2014—2024年，整合MODIS遥感生态、ERA5-Land气候、WorldPop人口、夜间灯光、OSM道路及文化遗产生成等多源数据，并回溯气象至1981年、遥感至2001年建立基线；研究队列为云南连续碳酸盐岩覆盖率超过27%的9个州（市）。ER以遥感生态指数（Remote Sensing EcoLogical Index, RSEI）为基础，按RAR构建干/湿季动态指标；HAI按显式/隐式维度经指标间相关性重要性准则（Criteria Importance Through Intercriteria Correlation, CRITIC）赋权，再用Pooled GAM回归海拔与坡度取残差得pHAI；采用Sen斜率与Mann-Kendall检验、贝叶斯优化XGBoost与SHAP、增量Moran's I进行趋势、非线性和尺度诊断。

研究结果：
3.1 HAI与ER评估结果：基于最大信息—最小不稳定性准则采用13年统一窗口计算动态ER；ER的Moran's I旱季0.54—0.68、湿季0.52—0.60，L3—L4约占60%，旱季低韧性区比例更高且更波动；Pooled GAM显示原始HAI与pHAI的R²=0.896，pHAI相对HAI平均降低23.7%，高值区由昆明扩散，L1约40—50%。
3.2 时空动态与错配：重心分析显示pHAI持续南移约10 km，旱季ER重心2019—2024年向西北移动>11 km；标准差椭圆主轴近东西向，pHAI轴偏约26°；ER呈“高位波动—相对稳定—震后缓慢恢复”三阶段并总体略降，旱季标准差0.017高于湿季0.009；pHAI呈升—降两阶段，二者明显异步。
3.3 多尺度效应与时空相关：3.3.1 空间尺度依赖：Spearman相关在网格、行政区、子流域比较显示网格4 km与8 km局部峰值，子流域层级相关稳定，10级子流域（平均约140 km²）被选为过程尺度；双变量Moran's I LISA在8 km网格与10级子流域干湿季聚类动态高度一致。3.3.2 时间趋势一致性：Sen斜率与MK检验显示pHAI略增、两季ER下降且旱季更明显，8 km网格与10级子流域的Sen斜率Moran's I差异小；旱季ER趋势周边多升、昆明核心显著降，湿季呈南升北降纬度梯度；ML、MM、MH组合占pHAI-ER交互>90%，双变量散点R²<0.1提示非线性阈值而非简单剂量效应。
3.4 共线性控制与全局重要性：3.4.1 因子相关与分布：因子间相关<0.9，干湿季ER Sen斜率仅中度相关（r=0.40），pHAI标准差与均值r=0.87，干湿季ER均值r=0.89。3.4.2 模型性能与全局重要性：干季模型需更深树结构与较低学习率，湿季模型结构较浅、学习率较高；SHAP显示旱季SRI趋势、Tmin趋势、DEM、Zssrd重要，人类活动因子贡献16.31%、生态基线6.04%；湿季ER均值、SRI量级、DEM重要，人类活动因子贡献28.28%、生态基线19.73%；排除控制变量后，旱季ER基线>pHAI均值>pHAI波动，湿季ER基线>pHAI Sen趋势>pHAI总和。
3.5 SHAP偏依赖与空间聚集：3.5.1 局部依赖分析：ER基线与ER稳定正相关，阈值约0.54（旱）与0.51（湿）以上进入自维持；旱季pHAI强度呈倒U形，极高强度>489.00 person/km²转正，波动STD>3.52为正，但pHAI总和>18×10³ person转负、趋势>0.07迅速转负；湿季趋势>0.04稳定为正、总和>19×10³为正，但强度>141.12 person/km²受抑制、波动高时迅速转负。总体上旱季受益空间集聚和短期方差但牺牲累计总量与加速趋势，湿季受益适度增长与规模效应但敏感强度峰值与脉冲波动。3.5.2 驱动带宽与情景异质性：模型残差Moran's I旱季0.004、湿季0.101，增量Moran's I将驱动尺度分为流域（<80 km）、局地（80—200 km）、半全球（200—400 km）、全球（>400 km）；控制变量呈大尺度平滑背景，pHAI与ER基线呈季节弹性。湿季ER基线为全球尺度441.43 km，干季收缩至局地155.80 km；pHAI强度旱季半全球207.73 km、湿季流域25.97 km；pHAI趋势旱季流域77.90 km、湿季局地181.76 km；pHAI波动旱季流域64.92 km、湿季流域12.98 km，显示湿季“广防御—局部风险”、干季“防御收缩—风险扩张”。

讨论部分总结：该研究认为，RAR动态框架与Sen斜率揭示区域ER略降，说明韧性提升不能简单等同于植被覆盖或生产力增加；旱季信号与岩性控制及西南森林韧性下降报道一致。喀斯特薄土层与地下管道造成快速水文漏失，放大季节干旱与人类用水叠加的复合扰动。低ER区集中于城市扩展区与石漠化带，高ER区位于低pHAI边缘；重心与椭圆表明pHAI与ER并非同步线性权衡，长期轨迹呈“人类扩散—风险聚焦”的季节锁定，符合多稳态与阈值滞后框架。子流域是pHAI-ER交互的功能尺度，网格受可塑性面积单元问题（Modifiable Areal Unit Problem, MAUP）影响，子流域更能代表水文汇流过程；Pooled GAM显示地形主导人类活动分布，净化链平均降低HAI约23.7%，体现地形—社会双约束。pHAI贡献旱季约16%、湿季约28%，干季对趋势与总量敏感，湿季对强度与波动敏感；ER基线跨季关联最稳定并存在自维持阈值。管理启示强调ER基线为长期韧性结构基础，湿季自维持阈值较低或利于基于自然的措施，干季为约束窗口，治理重点应从年活动规模转向季节节奏，水密集型产业宜湿季增总量但削峰值。局限包括年HAI未充分捕捉年内脉冲、滞后项未显式估计、气候异常时间窗与ER参考期不完全对齐、时间自相关需谨慎、机制阈值尚待云南外检验，观察设计不能区分人类影响、土地利用遗留、政策干预与内在敏感性差异，关联不等同因果，pHAI幅度降低23.7%与地形或净化链整体的关系仍需澄清；未来应以多尺度面板模型分离独立效应并检验快慢变量滞后反馈。

研究结论翻译：该研究基于多源数据与双季节情景建立“格局—驱动—响应”框架，系统揭示云南喀斯特地区pHAI与ER的非线性耦合结构。主要发现包括：（1）pHAI较原始HAI约低23.7%，这反映整个净化链的综合作用，地形亦起主导作用；旱季ER重心位移明显大于湿季，表明系统生态波动集中于缺水阶段。（2）ER呈“高位波动—相对稳定—震后缓慢恢复”循环并总体略降，pHAI呈升—降趋势，二者轨迹明显异步与错配。（3）10级子流域被识别为刻画人地交互的最佳空间单元，验证喀斯特人为扰动通过水文连通而非行政或网格边界传播，即合适尺度为过程尺度而非格局尺度。（4）pHAI贡献从约16%（旱季）升至约28%（湿季），确认活跃季物理放大；敏感机制反转：旱季系统对pHAI累计总量与趋势敏感，湿季对强度峰值与波动脆弱。（5）存在清晰ER基线阈值（约0.54旱季/0.51湿季），高于该水平与生态系统自组织及可持续性相关，但该结论在云南喀斯特之外的可迁移性仍需检验。（6）影响尺度随季节显著变化：控制变量呈广尺度关联，ER基线在干季收缩至局地尺度，人类强度风险足迹扩张至半全球尺度，形成结构性脆弱。总体而言，喀斯特区域ER呈空间结构化、季节调制的双重韧性架构：跨季持续性与高基线生态区固有抵抗能力密切相关，年内波动反映对比性季节耐受机制与人类活动空间强度、速率及累计足迹的非线性交互。该研究建立“时空格局—尺度依赖过程—非线性阈值响应”框架，为区域景观治理从静态面积控制转向适应性、空间显式、节奏敏感的人地交互管理提供定量基础。