加拿大落基山脉东坡高山树线生态交错带对气候变暖的响应时空动态

《Climate》：Spatiotemporal Dynamics of the Alpine Treeline Ecotone in Response to Climate Warming Across the Eastern Slopes of the Canadian Rocky Mountains Behnia Hooshyarkhah, Dan L. Johnson, Locke Spencer, Hardeep S. Ryait and Amir Chegoonian

　　

山岳生态系统易受气候变化影响，高山树线生态交错带（Alpine Treeline Ecotone, ATE）是气候驱动环境变化的重要响应指标之一。本研究基于1984至2023年加拿大落基山脉东坡（Eastern Slopes of the Canadian Rocky Mountains, ESCR）的长期ATE时空动态，旨在评估区域气候变暖是否在不同坡向和流域尺度上影响ATE的分布范围与海拔高度。研究人员将多年代际Landsat影像、ERA5-Land气温数据及地形变量整合于谷歌地球引擎（Google Earth Engine, GEE）框架中，利用高山树线生态交错带指数（Alpine Treeline Ecotone Index, ATEI）（一种用于捕捉过渡植被带的概率方法）对ATE进行制图。采用非参数统计、相关分析以及基于流域和坡向的比较评估时间趋势。结果显示，2023年ESCR的ATE总面积比1984年大约大13.3%，但其范围和海拔的演变呈现非单调特征，线性趋势分析未在整个研究期内检测到统计学上显著的增加或减少趋势。ATE的海拔和扩展表现出显著的空间异质性，在北坡和西北坡以及特定流域内变化较大。相比之下，夏季（7–9月）气温显著升高（+2.84 °C），超过全球陆地区域的变暖速率，且植被绿度（Normalized Difference Vegetation Index, NDVI）与气温呈强且统计学显著的正相关关系。上述发现表明，尽管气候变暖显著提高了植被生产力，但海拔方向的ATE动态仍保持空间异质性，且时间上与夏季气温趋势不同步。

**论文解读：高山树线生态交错带对气候变暖的时空响应——以加拿大落基山脉东坡为例**

**研究背景与科学问题**
高山树线生态交错带（ATE）是亚高山郁闭林与高山苔原之间的过渡地带，对气候变化极为敏感。尽管全球研究显示超过一半的ATE站点在过去一个世纪发生上移，但ATE响应存在强烈的空间异质性、非线性特征和时间滞后性。加拿大落基山脉东坡（ESCR）是重要的水文和生态区域，然而先前研究多局限于局地样线或短期观测，缺乏多流域、多坡向、数十年的区域尺度综合分析。此外，传统阈值法难以刻画ATE的渐变特性。因此，本研究旨在回答以下问题：1984–2023年ATE面积如何变化？海拔和分布在坡向与流域间呈现何种系统性差异？夏季气温变化与ATE范围、海拔及植被绿度有何关联？概率性ATEI方法相较阈值法能否更好地捕捉渐变动态？该论文发表于《Climate》期刊。

**主要技术方法**
研究人员整合了Landsat 5 TM和Landsat 8 OLI的30 m地表反射率数据（1984–2023），基于最佳可用像素法构建9个多时相合成期（每期3–5年）。采用SRTM数字高程模型提取地形数据，ERA5-Land再分析数据获取夏季（7–9月）月均2 m气温。利用谷歌地球引擎（GEE）计算归一化差异植被指数（NDVI）最大值，并构建高山树线生态交错带指数（ATEI），该指数包含NDVI-海拔梯度分量、中等NDVI分量和空间协变分量，通过二项逻辑回归生成概率图。采用分层随机抽样选取386个验证像素，分类总体精度81.30%。统计方法包括Mann–Kendall检验、Spearman相关、LOESS平滑及HAC标准误回归。

**研究结果**
*3.1 ATEI随时间分布变化*：基于ATEI≥0.7阈值，ATE面积从1984年约1494 km²增至2023年约1693 km²，净增13.32%。但Mann–Kendall检验显示无显著单调趋势（p>0.05），表明扩展呈非单调、间歇性特征。
*3.2 不同坡向海拔变化*：各坡向平均ATE海拔相对一致（2067–2160 m），北坡海拔始终最低。1984–2023年北坡净上升约+40 m，西北坡+27 m，东北坡+18 m，但各坡向均无显著单调趋势。2000–2007年期间所有坡向出现短暂下降，2007年后快速回升。
*3.3 不同流域海拔变化*：流域尺度平均上移23.41 m，但单样本t检验不显著（p=0.108）。Bow流域上移最大（+65.79 m），Red Deer流域几乎不变（?1.30 m），Oldman、North Saskatchewan和Athabasca流域分别上升+15.02、+13.20、+24.32 m。
*3.4 流域与坡向交互对ATE扩展的影响*：ATE扩展面积在阿萨巴斯卡流域最大（264.00 km²），其次北萨斯喀彻温（258.55 km²）和Bow（246.25 km²）。按流域面积归一化后，Bow流域扩展强度最高（0.0277 km²/km²）。北坡和西南坡平均扩展较大，但坡向效应受流域地形强烈调制。
*3.5 ATE变化与7–9月等温线的关系*：LOESS趋势显示1984–2023年夏季平均气温上升约+2.1 °C（每十年+0.48 °C），Mann–Kendall检验显著（τ=0.313，p=0.0045）。但ATE面积与气温的Spearman相关系数（ρ=0.50）在考虑时间依赖后不显著（p=0.941），一阶差分和滞后分析也未发现显著同步关系。ATE海拔与气温呈中等相关（ρ=0.427）但未达显著（p=0.252）。区域升温幅度（+2.84 °C）超过全球陆地平均异常变化（+1.46 °C）。
*3.6 NDVI与夏季平均气温的关系*：NDVI与气温的2年移动平均趋势均拟合良好（R²=0.69和0.60），Spearman相关系数ρ=0.71（p=0.0003），表明气温升高与植被绿度显著正相关。

**讨论与结论**
讨论部分指出：ATE面积净增但无单调趋势，符合ATE非线性、滞后响应的全球观察；坡向差异受局地微气候调节，但无显著趋势佐证非热力因子的约束；流域间异质性凸显区域地貌和干扰史的作用；ATE范围与海拔对升温的响应不同步，体现了生态惰性；NDVI显著上升说明植被生产力已率先响应，可作为气候变化的早期指示。
研究结论翻译：本研究利用Landsat时间序列数据和概率性ATEI方法，调查了1984–2023年ESCR的ATE多年代动力学。通过整合植被绿度梯度、中间植被状态及与海拔的空间协变性，该途径实现了ATE作为过渡带而非固定边界的空间显式制图。结果揭示了ATE在过往四十年间净向上迁移，伴随显著空间异质性：在较冷、日照较少坡向变化更强，其他地形设置下变化较温和；时间上呈非单调性，短期收缩穿插于长期扩展趋势中。方法上，ATEI等基于概率和梯度的框架通过捕捉位移发生前的渐变结构变化，相比传统阈值法更具优势。研究强调多年代际、区域尺度分析对于理解复杂山区ATE行为的价值。局限性包括验证主要依赖近期高分辨率影像、未直接量化干扰过程，以及多年合成影像无法解析年际超前-滞后动态。