非线性气候-生态系统功能相互作用决定了祁连山地区地上森林生物量的动态变化：遥感和可解释机器学习的启示

《Ecological Indicators》：Nonlinear climate–ecosystem functional interactions shape aboveground forest biomass dynamics in the Qilian Mountains: Insights from remote sensing and explainable machine learning

森林地上生物量（AGB）时空动态及其对气候因素的响应研究——以祁连山天祝县为例 通过多源遥感数据与机器学习模型融合，揭示天祝县森林AGB的年际变化规律（2009-2023）及对降水、温度的非线性响应机制。研究表明，AGB总体呈缓慢上升趋势（年均增速0.0246 Mg/ha），存在显著空间异质性，西南高海拔区AGB较高（均值80 Mg/ha），东北部因人类干扰和地形限制较低。非线性分析表明，AGB在年降水330-460毫米、温度-2至2摄氏度范围内达到峰值，单因素线性相关弱（降水r=0.045，温度r=0.009），但交互作用（T×P）主导空间变异（SHAP值最高）。滞后效应分析显示，降水1-2年滞后效应显著（R2=0.794），而温度效应较弱且周期性波动。该研究为高寒山地森林碳汇评估和管理提供量化依据。

　　
祁连山天祝藏族自治县森林地上生物量（AGB）的时空动态与气候响应机制研究

### 研究背景与科学问题
祁连山作为青藏高原东北缘的重要生态屏障，其森林生态系统在碳汇功能与气候变化响应方面具有典型性。随着全球气候变暖和极端天气事件频发，揭示森林AGB的时空演变规律及其对气候因素的响应机制，对指导区域生态修复、碳汇监测和适应性管理具有重要科学价值。本研究聚焦天祝县这一冷-旱山地森林生态系统，通过整合多源遥感数据与多维生态参数，系统解析AGB的动态特征及气候驱动机制。

### 研究方法与技术路线
1. **数据集成与预处理**
- 采用Landsat 7/8时序影像（30m分辨率）提取光谱反射率，结合MODIS LAI/FPAR（1000m分辨率）构建植被指数体系。
- 通过Google Earth Engine平台完成影像融合、云掩膜及标准化处理，消除地形阴影和大气干扰。
- 整合气象数据（1901-2024年月尺度）、数字高程模型（DEM）及人类活动 footprint指数，形成包含6大维度42项指标的数据库。

2. **机器学习模型构建**
- 采用LASSO回归进行特征筛选，通过方差膨胀因子（VIF）控制多重共线性，最终保留13项核心变量（包括NDVI、FPAR、DEM等）。
- 构建多元回归（MLR）、随机森林（RF）和极端梯度提升（XGBoost）三种模型，通过交叉验证（10折）评估模型泛化能力。
- 选择XGBoost+LASSO组合模型作为最优方案，因其平衡了预测精度（R2=0.68）与计算效率（训练时间1.48秒/次），并通过SHAP框架实现特征解释。

3. **气候响应解析方法**
- **线性关联分析**：采用偏相关系数评估温度、降水与AGB的单变量关系，发现降水呈微弱正相关（r=0.045），温度无显著独立效应。
- **非线性阈值分析**：基于局部加权回归（LOWESS）拟合单因子响应曲线，揭示降水存在330-460mm的“最优湿润窗口”，而温度-2℃至2℃区间促进生长。
- **交互效应量化**：引入温度-降水交互项（Tem_pre_interaction），通过SHAP值分解显示其贡献度（1.85）超越单一因子（降水1.72、温度1.38），证实水热协同驱动机制。
- **滞后效应分析**：采用最小二乘法评估1-10年滞后气候因子影响，发现降水滞后效应显著（R2=0.794），短期1-2年正反馈为主，3-4年负反馈占优，呈现“先促后抑”的波动模式。

### 关键研究发现
1. **时空动态特征**
- 2009-2023年AGB均值增长率为0.0246 Mg·ha?1·yr?1，西南高海拔区增速显著（年增0.05 Mg·ha?1），东北部因人类干扰（过度放牧、道路建设）出现局部衰退。
- 空间分布呈现“西南高-东北低”格局，峰值区域集中在海拔3200-4500m的云杉林带，最低值区（<80 Mg·ha?1）多分布于海拔<2500m的次生林与荒漠过渡带。

2. **非线性气候响应机制**
- **降水效应**：在年降水280-460mm区间AGB增幅达峰值（SHAP值+2.4），超过460mm后转为抑制（SHAP值-1.2），揭示“湿优旱劣”阈值规律。
- **温度效应**：-2℃至2℃区间呈现正响应（SHAP值+1.8），低温区（<-2℃）负贡献率达93%，高温区（>2℃）因蒸腾加剧导致抑制效应（SHAP值-1.5）。
- **交互作用**：温度-降水协同效应主导空间异质性，中温（-0.67℃-1.34℃）+中湿（412-452mm）区域AGB达区域均值最高（82 Mg·ha?1）。

3. **空间异质性分异规律**
- 按温度（低/中/高）与降水（低/中/高）三维度划分为9种气候区：
- **高湿-高温区**（Tem>1.34℃+Pre>48.80mm）：AGB均值80.5 Mg·ha?1，但高蒸发量导致负交互效应（SHAP值-2.2）。
- **冷湿过渡带**（Tem<-0.67℃+Pre>451.91mm）：云杉林通过深根系（>2m）和窄叶结构（叶面积指数LAI 4.2-5.8）维持稳定AGB（82 Mg·ha?1）。
- **暖旱胁迫区**（Tem>1.34℃+Pre<412.24mm）：杨树人工林因蒸腾耗水（年蒸发量>2000mm）出现生长抑制（SHAP值-1.8）。

4. **时间滞后效应特征**
- 降水滞后效应呈现“短促正、中长负”规律：1年滞后促进（β=2.07）、3年滞后抑制（β=-0.17）、7年滞后再促进（β=0.14），揭示土壤水分“记忆”周期约3-4年。
- 温度滞后效应表现为“双峰响应”：2年滞后正效应（β=6.68）与7年滞后正效应（β=6.60）对应植被生理恢复期，而3年滞后负效应（β=-7.03）反映干旱胁迫累积。

### 理论创新与生态启示
1. **水热协同驱动模型**
揭示冷-旱区森林AGB存在“双阈值”调控机制：
- **水分阈值**：年降水需维持在330-460mm区间，低于280mm或高于550mm均导致抑制。
- **温度阈值**：最适生长温度为-2℃至2℃，极端低温（<-2℃）或高温（>6℃）均抑制光合作用。

2. **生态适应机制解析**
- **种-境协同**：云杉（Picea crassifolia）通过窄冠幅（叶面积指数LAI 4.2-5.8）和深根系（平均根系深度1.2m）适应低温环境，其优势种地位强化了西南区AGB稳定性。
- **管理-生态互馈**：人工造林区（如2017年划定的国家公园核心区）因均匀立地条件，表现出比自然林更高的抗逆性（SHAP值+1.5）。

3. **碳汇管理优化路径**
- **分区治理**：在西南湿润-中温区（占县域40%面积）优先实施林分 densification（目标AGB>100 Mg·ha?1），东北干旱区（30%面积）应加强保水措施（如土壤改良、滴灌）。
- **灾害预警**：建立降水-温度交互预警模型，当连续2年降水>500mm或温度>6℃时触发生态风险预警。
- **碳汇增强**：通过营造云杉-冷杉混交林（目标胸径>25cm），可提升单位面积固碳量（+18%）。

### 方法论贡献
1. **多尺度特征融合技术**
整合30m光学影像（地形、植被结构）、1km MODIS LAI（冠层动态）和1000m气象再分析数据（气候记忆），突破单一数据源分辨率限制。

2. **SHAP-交互增强解释框架**
开发“特征重要性-交互敏感度”双维度分析模型，量化：
- 单因子贡献率：降水（41%）、地形（29%）、人类活动（9%）
- 交互因子贡献率：水热协同效应占比达38%

3. **动态阈值校准算法**
通过蒙特卡洛交叉验证（1000次重采样），确定气候阈值置信区间：
- 降水最优区间（330-460mm）：95%置信区间[280-500mm]
- 温度最适区间（-2℃-2℃）：90%置信区间[-5℃-5℃]

### 应用前景
1. **碳汇监测**：构建10km×10km网格化AGB动态数据库，精度达15 Mg·ha?1误差范围，可支撑区域碳汇年评估（分辨率1km）。
2. **生态修复**：基于气候-土壤-植被耦合模型，制定分区修复方案：
- 湿润区（>460mm降水）：推广云杉+高山松混交林（固碳量+22%）
- 干旱区（<280mm降水）：实施保水措施+根系固氮菌剂（土壤持水量+35%）
3. **气候适应设计**：提出“冷区保水-暖区控温”管理范式，在低温胁迫区（海拔>3500m）优先配置耐寒树种（如 Betula platyphylla），在高温区（<2500m）增设蒸腾抑制剂。

### 局限与展望
1. **数据局限性**：30m分辨率限制了对个体林分差异（如胸径>30cm的云杉与针叶林）的监测精度，建议补充LiDAR点云数据（1m分辨率）。
2. **模型泛化性**：需在黄土高原其他冷-旱区（如武威、张掖）验证模型稳定性，计划开展多站点联合建模。
3. **机制深化研究**：建议结合土壤同位素（δ1?O）分析水分来源，利用放射性碳测年（14C）验证模型预测精度。

本研究为《联合国气候变化框架公约》第28次缔约方大会（COP28）提交的“山地生态系统碳汇提升技术指南”提供了关键科学依据，特别在干旱-半干旱区碳汇管理方案设计中具有范式意义。