《Ecological Indicators》：Asynchronous convergence of annual satellite-derived forest productivity and within-year production allocation after fire: Counterfactual evidence from 7197 forest fires in China
编辑推荐：
野火可在燃烧后多年内持续影响森林生产力。准确表征火灾后生产力的变化及其恢复轨迹，对于理解生态系统对森林干扰的响应以及开展大尺度生态监测至关重要。卫星遥感能够在广阔空间范围和长时间尺度上提供连续的植被观测，已被广泛用于火灾后恢复评估。已有研究报告的恢复时间差
野火可在燃烧后多年内持续影响森林生产力。准确表征火灾后生产力的变化及其恢复轨迹，对于理解生态系统对森林干扰的响应以及开展大尺度生态监测至关重要。卫星遥感能够在广阔空间范围和长时间尺度上提供连续的植被观测，已被广泛用于火灾后恢复评估。已有研究报告的恢复时间差异显著。全球分析表明，大多数火烧迹地植被的生产力可在相对较短的时间内接近火灾前水平，而加利福尼亚州部分大型火灾的光合作用效应可持续超过二十年。中国东北地区的森林也表现出受造林历史、多年冻土退化和地形影响的漫长恢复轨迹。遥感指标代表不同的植被属性，其恢复轨迹未必一致；例如，光谱信号可能早于森林结构或生产力相关信号恢复。综合来看，这些研究表明大尺度火灾后恢复并不遵循单一的时间尺度，其评估取决于生态背景、所选指标和评估设计。估算恢复轨迹的方法也在不断进步。以受干扰森林的火灾前状况作为固定基线可能混淆火灾效应与气候和植被的同步变化，而以同期未燃烧森林作为时变反事实参考，能更恰当地估计火灾引起的偏离及其后续轨迹。恢复时间的评估之外，恢复债务（recovery debt）的概念将关注扩展至恢复至参考状态之前损失的生态系统功能，已被用于量化火灾后的长期生产力亏缺。这些进展改进了恢复参考和跨年累积变化的表征方式。然而，大尺度研究仍主要采用年度或生长季生产力来表征恢复，此类轨迹捕捉了生产力在连续年份间的广泛变化，但几乎不提供单个年份内部生产如何分配的信息。年度总初级生产力（gross primary productivity, GPP）整合了年内不同时段的生产，其年度值取决于生产持续时间和强度，因此相似的年度GPP可源自不同的年内生产分配。年内生产的时间安排也塑造了碳吸收的季节动态，不同阶段生产力的变化可进一步影响大尺度陆地碳循环的季节性。生产阶段之间亦相互关联，早期生产的变化可通过后续水分利用影响后期生产，使各阶段的增益与损失相互抵消。火灾可能进一步改变这种年内生产结构。已有研究表明，火灾驱动的植被组成变化可使北方森林中增强的光合作用集中在初夏和中夏，并改变季节性碳交换；中国西南地区火灾引起的物候时间变化也可持续数年。因此，即便年度生产力已接近同期未燃烧参考水平，其年内生产分配仍可能存在差异。仅基于年度总量的恢复评估可能忽视生产力恢复的年内结构及其关联的生态过程。年内生产分配是否在年度生产力再次接近同期参考水平后也趋于收敛，因此需要在大尺度火灾后恢复评估中进一步澄清。在此背景下，研究人员基于2005至2016年中国7197场森林火灾，将燃烧森林与同期未燃烧森林进行匹配，利用500米分辨率、8天时间步长的MOD17和PML-V2.2a数据重建了至火灾后第8年的反事实GPP轨迹。研究聚焦两个问题：第一，当卫星遥感的年度生产力接近同期未燃烧反事实水平时，年内生产分配是否同步收敛？第二，若二者不同步收敛，年内差异在不同时间定义和GPP产品间是否表现出一致的方向和时间位置？由于卫星GPP估算受输入数据、参数化和模型结构的影响，两种产品并行分析以区分跨产品一致的年内特征和对产品选择敏感的效应幅度。进一步分析用于表征年内差异的时间结构，评估可能导致该模式的因素，并探索其随火灾严重度和火灾前森林类型的变化。通过将年度生产总量和年内生产分配与同一反事实参考在共同时间支撑下进行比较，本研究旨在阐明年度聚合在多大程度上捕获了火灾后生产力恢复的时间信息，并为大尺度恢复评估中的时间尺度选择和多产品结果解读提供依据。
**论文解读**
**一、研究背景与问题提出**
野火对森林生产力的影响可持续燃烧后多年。准确表征火灾后生产力的变化及其恢复轨迹，对于理解生态系统对森林干扰的响应以及开展大尺度生态监测具有关键意义。卫星遥感凭借其广覆盖、长时序的连续植被观测能力，已被广泛用于火灾后恢复评估。然而，已有研究报告的恢复时间差异显著：全球尺度分析显示多数火烧迹地植被生产力可在较短时间内接近火灾前水平，而加利福尼亚州部分大型火灾的光合效应可持续超过二十年；中国东北森林也表现出受造林历史、多年冻土退化和地形影响的漫长恢复轨迹。不同遥感指标代表不同的植被属性，其恢复轨迹未必一致，光谱信号可能早于森林结构或生产力相关信号恢复。这些研究共同表明，大尺度火灾后恢复并不遵循单一时间尺度，其评估取决于生态背景、所选指标和评估设计。
在方法层面，以火灾前状况作为固定基线可能混淆火灾效应与气候和植被的同步变化；而采用同期未燃烧森林作为时变反事实参考（counterfactual reference），能更恰当地估计火灾引起的偏离及其后续轨迹。此外，恢复债务（recovery debt）概念将关注扩展至恢复至参考状态前损失的生态系统功能，已被用于量化火灾后的长期生产力亏缺。然而，当前大尺度研究仍主要采用年度或生长季生产力来表征恢复——此类轨迹捕捉了生产力在连续年份间的广泛变化，但几乎不提供单个年份内部生产如何分配的信息。年度GPP整合了年内不同时段的生产，相似的年度GPP可源自不同的年内生产分配。年内生产的时间安排塑造碳吸收的季节动态，早期生产的变化可通过后续水分利用影响后期生产，各阶段之间相互关联。火灾可能进一步改变这种年内生产结构：火灾驱动的植被组成变化可使光合作用集中到初夏和中夏并改变季节性碳交换，火灾引起的物候时间变化也可持续数年。因此，即便年度生产力已接近同期未燃烧参考水平，其年内生产分配仍可能存在差异。仅基于年度总量的恢复评估可能忽视生产力恢复的年内结构及其关联的生态过程。年内生产分配是否在年度生产力再次接近同期参考水平后也趋于收敛，成为亟待澄清的科学问题。
**二、研究设计与关键技术方法**
针对上述问题，研究人员利用2005至2016年中国7197场森林火灾开展了反事实分析。燃烧像素来自500米分辨率的MODIS MCD64A1火烧迹地产品，火灾前森林覆盖度（forest fraction）由30米分辨率中国土地覆盖数据集（CLCD）聚合至MODIS像元网格，仅保留火灾前森林覆盖度不低于0.70的像元，并要求火灾前三年内无燃烧记录且火灾年至2024年无后续燃烧。燃烧像元间距不超过1.5公里且燃烧日期相差不超过7天的被归并为单一火灾事件。火灾年份限定为2005至2016年，使得每个事件至少拥有三个完整火灾前年份和八个完整火灾后年份。
反事实对照设计是本研究的关键方法。每个燃烧像元从候选搜索半径内的未燃烧森林像元中匹配对照，要求候选像元在2002至2024年间无MCD64A1燃烧记录且满足相同的森林覆盖度阈值。匹配使用了13个火灾前协变量，涵盖地形（海拔、坡度、北向度、东向度）、植被状态（MODIS树冠覆盖度、火灾前EVI、LAI及GPP水平和趋势、森林年龄）和30米森林覆盖度，并优先匹配相同火灾前森林类型。最终7197场火灾事件匹配至21035个事件-对照对，平均每事件2.92个对照。在9个候选设计（森林覆盖度阈值0.60/0.70/0.80与搜索半径2/5/10公里组合）中，基于火灾前协变量平衡（标准化均值差异|SMD|<0.10）选择0.70/5公里设计，选择过程未使用火灾后结果。
GPP数据采用两种结构不同的卫星遥感产品进行并行分析：MOD17基于光能利用率（light-use-efficiency）模型，PML-V2.2a基于Penman-Monteith-Leuning框架联合估算蒸散发和GPP。两种产品空间分辨率均为500米、时间分辨率均为8天，使用2002至2024年观测数据。GLASS-GPP仅用于年度尺度比较和反事实设计敏感性分析。火灾严重度以MODIS dNBR（差分归一化燃烧比）表征。年内生产阶段依据各事件未燃烧对照的火灾前平均生产曲线定义：以累积生产达到年度总量25%和75%的日期为界，将生产年划分为早期（early）、核心（core）和晚期（late）三个阶段。核心分析采用预设的年内生产分配对比指标D（早期与核心阶段相对效应之差），以火灾后第6至8年平均值作为后期汇总。
**三、主要研究结果**
**3.1 年度生产力接近反事实水平但月度效应不均**
所有三种卫星遥感GPP产品在火灾后均表现出年度偏离同期反事实水平的显著减小。火灾后第1年，MOD17、PML-V2.2a和GLASS-GPP的年度相对效应分别为?6.60%、?4.83%和?8.05%。至火灾后第6至8年，MOD17平均年度相对效应A为?1.13%（95% CI，?1.50%至?0.80%），PML-V2.2a为+0.59%（?0.15%至+1.31%），GLASS-GPP为?2.19%（?2.85%至?1.62%），三者的组均值均落入±5%的操作比较范围内。月度分解显示两种产品均存在明显的年内变异：MOD17在火灾后第8年各月效应均为负值，最大负效应出现在3月和4月；PML-V2.2a在1至4月效应为负，6至11月转为正值。两种产品中，负效应均集中在年初，PML-V2.2a还表现出年中至年末的正效应。
**3.2 年内生产分配对比缩小但未与年度效应同步收敛**
预设的年内生产分配对比D在火灾后第1至8年总体趋向零值：MOD17的D从?8.51个百分点缩小至?1.90个百分点，PML-V2.2a从?8.19缩小至?2.19个百分点，绝对幅度分别减小77.7%和73.2%，大部分变化发生在前三年。然而至火灾后第6至8年，MOD17平均D为?1.98个百分点（?2.59至?1.44），PML-V2.2a为?2.45个百分点（?3.30至?1.78），早期阶段效应在两种产品中均低于核心阶段。该方向在四种预设时间定义下均一致：所有8个产品×定义的D估计值均为负且95%置信区间完全位于零以下。跨年分析进一步显示，在根据火灾后第7年年度效应是否落入±2.5%、±5%或±10%范围筛选的事件中，第8年D始终保持负值；连续回归分析中，当拟合的第7年年度效应为零时，第8年D仍为负值（MOD17为?1.77个百分点，PML-V2.2a为?2.14个百分点）。
**3.3 年内分配对比超过设计内部变异并伴随早期阶段时间和幅度差异**
对照-对照随机化（control-to-control randomization）显示，MOD17和PML-V2.2a的国家均值经验零分布中心分别为?0.276和?0.254个百分点，标准差分别为0.153和0.146个百分点，而观测D值为?1.98和?2.45个百分点，分别偏离零分布中心11.2和15.0个标准差（P=0.0002）。生产窗口的时间变化集中在起始端：火灾后第6至8年，q10（累积生产达10%的日期）延迟约0.87天（MOD17）和1.05天（PML-V2.2a），而q90（累积生产达90%的日期）变化较小。将未燃烧对照曲线整体延迟一天可产生观测D的31%至51%，对应的局部线性等效位移约2.0至3.2天。同时，拟合的早期阶段幅度低于核心阶段幅度，早期减去核心的幅度对比在MOD17为?0.0141、PML-V2.2a为?0.0184。8天数据尚无法可靠区分小幅时间位移与阶段特异性幅度差异的相对贡献。PDSI（帕尔默干旱严重度指数）与核心阶段效应无显著关联；严格筛选仅保留主要反演算法像元后D仍为负值；火灾前三年未检测到与火灾后D同方向的趋势。探索性分层分析中，D随火灾严重度增加趋于更负，火灾前常绿森林的D均值低于落叶森林。
**3.4 群体层面年内模式跨产品一致而事件层面评估和绝对幅度差异显著**
MOD17和PML-V2.2a在群体均值分解中显示相同方向模式：火灾后第6至8年，峰值速率通道均为正值（MOD17为+0.071，PML-V2.2a为+0.059 kg C m
?2 yr
?1），有效持续时间通道均为负值（?0.090和?0.049 kg C m
?2 yr
?1），与早期生产阶段相对效应更低的模式方向一致。然而事件层面一致性明显较弱：两种产品对同一事件主导通道一致的仅占51.8%–53.8%；更换GPP产品导致39.3%事件效应符号反转，更换时间定义为32.0%，更换反事实设计为27.1%。固定事件队列代表性像素的8年累积GPP差异在MOD17（?0.470 Tg C）、PML-V2.2a（?0.170 Tg C）和GLASS-GPP（?0.796 Tg C）之间最大相差4.7倍。总体而言，两种GPP产品在群体层面的年内模式方向和时间位置上的一致性最强，而事件层面评估和GPP绝对幅度对产品选择和分析设置更为敏感。
**四、讨论与结论**
本研究表明，年度卫星遥感生产力和年内生产分配在火灾后均逐步向同期未燃烧反事实水平收敛，但二者的收敛轨迹并不同步。从火灾后第1年至第8年，年内分配对比缩小超过70%，表明年内生产分配本身也随时间显著恢复；然而至火灾后第6至8年，年度效应已较小，早期生产阶段相对于核心阶段仍存在约2个百分点的负分配对比。更重要的是，在跨年连续分析中，当年度效应拟合为零时，后续年份的分配对比仍为负值。因此，年度GPP接近反事实水平表明年度总生产已接近参考状态，但这并不意味年内生产分配已同步收敛。该差异在四种时间定义和跨产品配对上方向一致，且观测值明显偏离对照-对照随机化的经验零分布，表明该差异不依赖特定年度窗口、阶段边界或单一GPP产品。
年内分配差异的特征表现为生产窗口起始端约1天的小幅时间延迟和阶段间幅度差异，但8天数据无法可靠分离二者贡献。火灾后植被组成变化为阶段间生产强度差异提供了合理解释：常绿和落叶组分对年度生产的贡献时间不同，火灾后功能群相对贡献的变化可改变年内生产分配。本研究中D随火灾严重度增加而更负、火灾前常绿森林中D更低的模式，在方向上与早期贡献组分受到更强影响并通过其他组分补偿恢复的时间重组一致。然而由于缺乏物种和功能型直接观测，该机制有待未来实地数据验证。水分利用的季节性结转作为另一种可能的补偿机制，在PDSI分析中未获得跨产品支持。
跨产品比较表明，MOD17和PML-V2.2a虽然月度轨迹不完全相同，但在群体层面早期阶段效应低于核心阶段、最大负月度效应集中在年初的方向一致性，以及峰值速率-有效持续时间分解的相同方向模式，降低了年内模式仅源于单一产品特定公式的可能性。但两种产品均依赖遥感和气象信息，其一致性应理解为跨产品可重复性而非对实际生态系统GPP的独立验证。产品间绝对幅度差异大（8年累积差异最大相差4.7倍），事件层面效应符号反转比例高，表明绝对火灾后GPP变化比群体层面年内模式更依赖产品选择。年内分配对比对应于年度对照GPP的约0.6%–0.7%，这一适中的年度尺度幅度表明，在年度总量接近反事实后，可检测的差异主要体现在年内生产的分布方式上。
研究的主要局限包括：两种GPP产品均为卫星遥感估算而非实地观测，产品间部分共享遥感和气象输入使其一致性不能替代独立验证；500米空间分辨率下燃烧像元内不可避免地混合存活树木、死亡植被、再生植被和裸地；8天时间分辨率限制了对年内过程的精确归因，无法可靠区分小幅起始端时间位移与早期阶段幅度降低；对单个事件的推断不确定性较大，严重度和森林类型对比仍属探索性质。
研究结论如下：基于中国7197场森林火灾的反事实轨迹分析，年度卫星遥感GPP与年内生产分配提供了相关但不同步的火灾后生产力恢复信号。至火灾后第6至8年，年度GPP已接近反事实水平，而早期阶段生产相对于核心阶段仍保持较低水平；该年内对比自火灾后第1年以来已大幅缩小，但在不同时间定义下及当前期年度效应拟合为零时仍保持负值。该模式伴随生产窗口起始端的小幅时间差异和较低的早期阶段生产幅度。跨产品比较显示群体层面模式的方向和时间位置的可重复性强于GPP绝对幅度和事件层面评估。年度生产总量表征了整合的年度状态，而年内生产分配提供了年内生产组织方式的补充信息，综合考虑两者可为大尺度火灾后生产力恢复提供更完整的表征，并为多产品评估结果的解读提供更清晰的基础。
