东北中国多年冻土区森林物候–活动层厚度关系的时空记忆与非线性转换

《Ecological Indicators》：Spatiotemporal memory and nonlinear transitions in forest phenology–active layer relationships in the permafrost region of Northeast China

　　

在全球变暖背景下，高纬度多年冻土正经历广泛而快速的退化，其主要表现之一是活动层厚度（ALT）持续加深，深刻重塑地下热水文状况与深根生境，并通过自下而上的方式调节寒区森林获取水分与养分的途径，进而改变其生长周期。森林物候作为追踪生态系统环境响应的敏感指标，直接调控陆地碳、水与能量通量。因此，准确阐明ALT波动与森林物候的关系，对评估寒区森林碳汇潜力及优化全球陆面模式中的冻土–森林参数化方案至关重要。尽管局地观测已报道ALT加深与森林动态的关联，但粗空间分辨率、对简单线性假设的依赖以及长期气候变暖的混杂效应，使空间异质性、时间记忆与非线性形式仍知之甚少。研究人员聚焦东北中国多年冻土区，集成多源数据、机器学习与物理模型，构建2000–2024年30米分辨率的森林物候与ALT数据集；通过去趋势、时滞与分段回归，定量刻画与ALT波动相关的空间分化物候关系及非线性 breakpoint。过去25年区域ALT显著加深（0.13 cm/yr）；去趋势前伴随生长季始期（SOS）提前（?0.17 天/yr）与生长季长度（LOS）延长（+0.16 天/yr）。对ALT与物候时序去趋势后，最优滞后–累积关联呈清晰南北对比：北部不连续多年冻土区SOS与ALT变率负相关，南部边缘带正相关；北部物候与多年累积ALT条件最强相关，南部主要关联当年变率。基于像元–年均值中心化的分段回归识别出指标特异性非线性 breakpoint：SOS与LOS约0.4 cm，EOS约?0.1 cm；跨 breakpoint 后SOS–ALT斜率由?2.68转为+0.24 天/cm，LOS–ALT斜率由+2.78减弱为+0.87 天/cm；更高ALT异常下LOS持续延长越来越多地与EOS推迟相关（+0.93 天/cm）。这些统计关系为未来基于过程模型与碳通量观测的碳固存评估提供了可检验假设。

研究背景方面，高纬度多年冻土退化伴随活动层厚度（ALT，active layer thickness）加深，改变地下热水文与深根生境，自下而上调节寒区森林水分养分获取并改变生长周期；森林物候直接调控陆地碳、水、能量通量。现有研究受限于1公里及更粗遥感产品，无法刻画亚公里尺度ALT空间异质性及物候的局部响应，且多依赖同期线性相关，忽视长期变暖背景下的时间记忆效应与可能的非线性转换。东北中国多年冻土区位于欧亚多年冻土带南缘，自北向南呈不连续（DPR）–稀疏（SPR）–孤岛状（IPR）分布，是检验空间观测尺度效应与非线性关系的理想区域。研究人员由此集成多源遥感、站点观测与机器学习，构建30米分辨率森林物候与ALT长时序数据集，并在去除共享变暖趋势后量化时空记忆与非线性 breakpoint。该研究发表于《Ecological Indicators》。

关键技术方法上，研究人员以2000–2024年30米Landsat合成反射率与光谱指数结合500米GLASS LAI（leaf area index，叶面积指数）纯像元样本，经机器学习检索与STARFM时空融合重建30米/8天LAI时序；以16个气象站0厘米地温计算地面解冻指数（GTI，ground thawing index），耦合机器学习空间化GTI与优化Stefan方程（ALT=E×√GTI，E因子按森林类型–坡向查找表赋值）估算30米ALT；基于稳定森林覆盖掩膜，用动态阈值法从平滑插值LAI提取SOS、EOS、LOS；对ALT与物候像元时序分别去线性趋势，遍历10种滞后–累积配置选最优相关系数，再对均值中心化异常做0.1 cm分箱与两段式分段回归，以AIC判据识别非线性 breakpoint；趋势分析采用Sen斜率与Mann–Kendall检验。

研究结果部分，2.1研究区表明东北北部（116.27°–130.38°E，46.72°–53.42°N）含呼伦贝尔高原与大小兴安岭，气候冷湿，多年冻土由北向南从DPR经SPR至IPR，平均地温?1.5至0℃、平均厚度约1.5米，热稳定性差。2.2数据来源显示研究人员整合Landsat、GLASS LAI、Sentinel-2、MODIS、气象站0厘米地温、野外ALT实测（46点定标、其余独立验证）、GLC_FCS30D植被类型与30米ASTER DEM。2.3.2空间建模表明，九种机器学习比选后以最优模型空间化GTI，再用E因子查找表（如DBF阳坡0.030、阴坡0.015）驱动Stefan方程逐像元算ALT。2.3.3 LAI重建与物候提取表明，STARFM以Landsat为空间基、GLASS为时间轨迹融合出30米/8天LAI，Savitzky–Golay滤波与线性插值至日值，动态阈值（SOS 50%、EOS 45% LAI幅值）提取物候，ENF未解析像元以同类型邻近填充（占比0.001146%）。2.3.4时间记忆与非线性分析表明，去趋势后2003–2024年遍历LAGrid选最优(L,A)；SOS在北部DPR最优关联为负向多年累积ALT，南部IPR为当年ALT正向；分段回归在SOS/LOS异常约+0.4 cm、EOS约?0.1 cm处ΔAIC<?10，确认非线性：SOS–ALT斜率?2.68转+0.24 天/cm，LOS–ALT +2.78转+0.87 天/cm，高ALT异常端LOS延长伴EOS推迟+0.93 天/cm。2.3.5趋势评估表明区域ALT以0.13 cm/yr加深（MK显著），原始序列SOS?0.17 天/yr、LOS+0.16 天/yr。

讨论与结论翻译总结：研究人员指出，粗分辨率会平滑ALT空间细节并抵消相反物候响应，去趋势后南北记忆结构分化揭示土壤热迟滞的累积效应——北部不连续冻土区物候更多记忆多年ALT状态，南部边缘带对当年ALT变率更敏感。分段回归给出的指标特异性 breakpoint 说明适度ALT加深可能通过增加深层水肥可及性促进SOS提前与LOS延长，但跨过临界异常后浅根区水分胁迫与冻融扰根削弱甚至逆转该关联，LOS延长在高ALT端转而依赖EOS推迟。这些统计关联非因果，但为过程模型（如冻土–植被耦合陆面模式）与通量塔观测的碳固存评估提供了可检验假设，尤其提示未来寒区森林碳汇预测需嵌入空间显式ALT、滞后记忆与非线性参数化，而非单一同期线性响应。