火烧严重度塑造亚热带马尾松人工林沿对比性土壤恢复路径的草本与灌木异步恢复

《Trees, Forests and People》：Fire severity shapes asynchronous herbaceous and shrub recovery along contrasting soil-recovery pathways in a subtropical Pinus massoniana plantation

　　

摘要
火后生态系统恢复取决于火烧严重度以及土壤与植被过程的耦合，但草本层和灌木层可能遵循不同的恢复路径。研究人员在华中地区大冶市的一片亚热带马尾松（Pinus massoniana）人工林中，研究了火烧严重度如何塑造表层土壤恢复和林下植被动态。在火后0.5年、1.5年和3.5年，分别对未燃烧参照样地、低严重度、中严重度和高严重度样地进行调查。最终数据集包含20个永久样地，在三个恢复阶段重复调查，共获得60个样地-阶段观测数据。研究人员采用线性混合效应模型（linear mixed-effects model, LMM）、主成分分析（principal component analysis, PCA）、扩展生物量模型（extended biomass model）和分段结构方程模型（piecewise structural equation modeling, piecewise SEM）分析了土壤理化性质、草本和灌木多样性及地上生物量。较高的火烧严重度与早期强烈的土壤位移相关。在火后0.5年，土壤有机碳（soil organic carbon, SOC）从未燃烧样地的37.33 ± 0.89 g kg^-1降至高严重度样地的22.72 ± 0.68 g kg^-1，而容重（bulk density, BD）从1.175 ± 0.013 g cm^-3增加至1.314 ± 0.018 g cm^-3。草本生物量呈现相反响应，从未燃烧样地的108.65 ± 30.94 g m^-2增加至高严重度条件下的442.92 ± 14.71 g m^-2。相比之下，灌木生物量在火后0.5年从720.92 ± 126.53 g/4 m²降至110.31 ± 2.95 g/4 m²，并且在3.5年后高严重度样地仍低于未燃烧参照。土壤主成分1（Soil PC1）解释了土壤性质变异的68.4%，代表了一个综合的土壤恢复梯度，其中土壤含水量（soil water content, SWC）、SOC和总氮（total nitrogen, TN）为正向载荷，而BD、pH和有效磷（available phosphorus, AP）为负向载荷。最终的分段结构方程模型显示出可接受的拟合度（Fisher’s C = 0.208, df = 4, P = 0.995），并揭示了对比性路径：火烧严重度与草本生物量呈正相关（标准化系数 = 0.486），但与灌木生物量呈负相关（-0.355），而土壤恢复与灌木生物量呈正相关（0.311）。这些结果表明，早期草本扩张可能掩盖灌木和土壤的延迟恢复，强调应将火后恢复评估为一种分层特异性的土壤-植被耦合过程，而非单一的绿度或生物量轨迹。

论文解读：火烧严重度调控下亚热带马尾松人工林的土壤-植被异步恢复机制

**一、研究背景与科学问题**

近年来，全球多个森林区域的火烧制度（fire regime）与火烧严重度（fire severity）正在发生显著变化，火后恢复的可预测性随之降低。虽然火是许多生态系统的自然干扰因子，但气候变暖与土地利用变化共同导致了更频繁、更持久或更强烈的火季。在景观尺度上，火烧严重度描述了火对植被、地表可燃物和土壤造成的生态改变程度，其空间分布往往形成未燃烧、低、中、高严重度斑块镶嵌的异质格局。同一火烧景观内，因初始干扰程度不同，即便处于相同区域气候和火后天气条件下，不同严重度斑块的土壤和植被恢复轨迹也可能存在差异。

表层土壤是连接火烧严重度与植被恢复的关键纽带，火通过消耗凋落物、细根和表层有机质改变土壤碳氮库，同时灰分沉降和矿物暴露可短暂提高pH和有效养分。然而，土壤各属性的响应方向并不一致，某一指标可能呈现养分富集而另一指标表现出退化，单一土壤变量难以准确反映再生植被所经历的整体土壤状况。在植被层面，火后恢复常表现出强烈的层特异性：草本植物因冠层开阔、凋落物减少、矿质土壤裸露和木本竞争减弱而快速定居；灌木则依赖地下器官存活、萌蘖能力、种子库和幼苗建成，需要更长时间重建地上结构。然而，既有研究多聚焦于总体植被盖度、光谱绿度或乔木更新，较少在相同的永久样地内明确对比草本层和灌木层，并将两层同时与同一综合土壤恢复梯度相关联。这一区分至关重要，因为草本生物量的快速增加可能仅代表对受干扰空间的早期占据，而非林下结构的真正恢复。因此，识别草本层与灌木层是否沿不同土壤关联路径异步恢复，对于准确评估火后生态系统恢复状态具有重要意义。

亚热带马尾松（Pinus massoniana）人工林是我国南方重要的针叶人工林类型，广泛用于荒山绿化、用材生产和退化红壤区生态恢复。此类人工林林冠结构相对简单、针叶凋落物连续、地表可燃物积累明显，林下植被发育较为完整，是研究层特异性火后恢复的代表性系统。2022年4月湖北省大冶市发生森林火灾，在马尾松人工林景观内形成了火烧严重度异质斑块，为在同一区域气候和林分背景下对比不同严重度样地的恢复过程提供了天然实验条件。

**二、研究设计与技术方法**

研究样地设置在湖北省大冶市（114°31′33″–115°10′38″ E，29°51′16″–30°19′41″ N），该区域属北亚热带季风气候，年均温17.5°C，年降水量约1495.2 mm。火后依据树干炭化高度、灌木层燃烧程度和凋落物消耗量三个现场指标，划分未燃烧参照（CK）、低严重度（LS）、中严重度（MS）和高严重度（HS）四个等级，每等级设置5个永久样地（约25.82 m × 25.82 m），共20个样地。在火后0.5年、1.5年和3.5年三个恢复阶段进行重复调查，形成60个样地-阶段观测记录。研究采用的样地来源为2022年4月大冶市野火形成的自然火烧斑块，属于分层观测比较设计而非完全随机化实验。

关键技术方法包括：采用线性混合效应模型（LMM）检验火烧严重度、恢复阶段及其交互作用对各土壤和植被指标的影响，以样地身份作为随机截距处理重复测量非独立性；采用主成分分析（PCA）对SWC、BD、pH、SOC、TN和AP六个土壤变量进行降维，提取综合土壤恢复梯度（Soil PC1）；采用扩展线性混合效应模型分析土壤恢复轴和多样性对草本及灌木生物量的解释力；采用分段结构方程模型（piecewise SEM）区分火烧严重度对植被生物量的直接与间接路径，并保留样地重复测量结构。

**三、主要研究结果**

**3.1 火烧严重度与恢复阶段对土壤性质的影响**

土壤性质的恢复轨迹受火烧严重度和恢复阶段共同调控，但不同变量的响应强度和持续性存在差异。恢复阶段对全部六个土壤指标均有显著影响，火烧严重度显著影响BD、pH、SOC、TN和AP，对SWC仅边缘显著（F = 3.11, P = 0.056）。火后0.5年，SWC从CK的23.32 ± 0.49%降至HS的18.40 ± 0.62%；BD从1.175 ± 0.013 g cm^-3增至1.314 ± 0.018 g cm^-3；pH从5.04 ± 0.04升至HS的5.50 ± 0.03；AP从12.44 ± 1.13 mg kg^-1增至HS的25.02 ± 0.73 mg kg^-1。SOC随严重度递增而显著降低，从CK的37.33 ± 0.89 g kg^-1降至HS的22.72 ± 0.68 g kg^-1；TN呈类似但较平缓的下降模式。至3.5年，多数土壤变量的等级间差异减弱，但HS的SOC（35.10 ± 0.27 g kg^-1）仍显著低于CK（38.96 ± 0.93 g kg^-1），表明高严重度对土壤碳库的影响可持续至早期恢复阶段之后。

**3.2 草本层恢复**

草本层对火烧严重度呈现强烈正响应。草本Shannon多样性受火烧严重度（F = 70.46, P < 0.001）、恢复阶段（F = 50.89, P < 0.001）及其交互作用（F = 7.30, P < 0.001）的显著影响。火后0.5年，草本Shannon从CK的1.439 ± 0.078增至MS的2.591 ± 0.209和HS的2.754 ± 0.056。草本生物量同样大幅增加，0.5年时从CK的108.65 ± 30.94 g m^-2升至HS的442.92 ± 14.71 g m^-2。虽然火烧样地的草本生物量在1.5年时有所回落，但整个研究期内始终高于CK。相对响应分析证实，所有火烧严重度等级在所有恢复阶段的草本Shannon和生物量均相对CK呈正值，早期阶段的正响应最强。这表明火后草本层的快速扩张代表了火烧后开放条件下对新生境的迅速占据。

**3.3 灌木层恢复**

灌木层恢复模式与草本层形成鲜明对比。灌木Shannon多样性受火烧严重度（F = 6.59, P = 0.004）、恢复阶段（F = 195.66, P < 0.001）及交互作用（F = 3.11, P = 0.016）的显著影响。火后0.5年，HS的灌木Shannon低于CK，至3.5年各等级差异缩小，但HS仍略低于CK。灌木生物量对火烧严重度呈现强烈负响应：0.5年时从CK的720.92 ± 126.53 g/4 m²降至LS的293.62 ± 28.00 g/4 m²、MS的212.00 ± 5.62 g/4 m²和HS的110.31 ± 2.95 g/4 m²。1.5年时LS和MS有所回升，HS仍明显偏低；3.5年时LS和MS基本恢复，但HS仍低于其他等级。相对响应分析显示，灌木生物量在0.5年时相对CK呈强烈负值，HS在3.5年时仍为负响应。灌木Shannon的响应较草本Shannon更弱且变异更大，表明灌木层的恢复是缓慢的结构重建过程。

**3.4 土壤恢复梯度与生物量模型**

PCA提取的Soil PC1解释了土壤变量总变异的68.4%，PC2解释了12.4%。SWC、SOC和TN在PC1上为正向载荷，BD、pH和AP为负向载荷，因此较高的Soil PC1代表土壤含水量较高、碳氮积累较多、容重较低、火致pH和磷素效应较弱的恢复状态。PCA双标图显示，CK样地集中于Soil PC1正向一侧，火烧样地尤其HS偏向负向一侧，证实该梯度捕获了火烧严重度增加导致的多元土壤位移。扩展LMM结果表明，草本生物量与火烧严重度呈正相关（HS系数最大），与Soil PC1呈负相关，与草本Shannon呈正相关；灌木生物量则与火烧严重度呈负相关（HS负效应最强），与Soil PC1呈正相关。这表明同一土壤恢复梯度对草本和灌木生物量具有方向相反的关联效应。

**3.5 分段结构方程模型**

最终分段SEM拟合良好（Fisher’s C = 0.208, df = 4, P = 0.995），各内生变量的R²分别为：土壤恢复0.671、草本Shannon 0.541、灌木Shannon 0.644、草本生物量0.811、灌木生物量0.821。火烧严重度与Soil PC1呈显著负关联（标准化系数 = -0.558, P < 0.001），恢复年份与Soil PC1呈显著正关联（0.606, P < 0.001）。火烧严重度显著正向预测草本Shannon（0.666）和草本生物量（0.486），而Soil PC1对草本生物量呈负向关联（-0.274）。灌木Shannon受恢复年份（0.495）和Soil PC1（0.429）正向预测；灌木生物量与火烧严重度呈负关联（-0.355），与恢复年份（0.238）、Soil PC1（0.311）和灌木Shannon（0.285）呈正关联。此外，灌木Shannon对草本生物量存在显著的跨层负关联（-0.385），表明灌木结构恢复可能对草本扩张产生抑制作用。将火烧严重度作为分类预测变量的敏感性检验结果与有序模型一致，支持了研究结论的稳健性。

**四、讨论与研究结论**

火后表层土壤响应呈现短期对比性强、随后逐步趋同的特征。pH和AP的早期升高最可能源于灰分驱动的养分脉冲，此类效应不稳定，可随淋溶、径流、微生物固持和植被吸收而消退。SOC和TN的恢复更直接依赖有机质输入，高严重度火消耗了凋落物、细根和表层有机质，其恢复需要凋落物归还、根系生产和微生物加工的持续作用。BD的恢复相对滞后，表明土壤物理结构的重建往往慢于部分化学指标的快速变化，化学性质趋同不应被视作土壤完全恢复的证据。

草本与灌木的异步恢复反映了干扰释放与结构重建的过渡过程。火后冠层开阔、凋落物减少和木本竞争减弱为草本快速生长创造了条件，而灌木恢复依赖于幸存个体、萌蘖能力、种子萌发和木质结构逐步积累。本文的SEM结果显示，草本生物量在土壤位移更强的条件下反而更高，而灌木生物量更密切关联于土壤恢复状态和时间推移，证实了两种生活型沿不同土壤关联路径恢复。相对响应分析进一步表明，仅依据草本多样性或生物量会得出火后快速恢复的结论，而灌木生物量揭示了持续的结构亏缺，高严重度斑块尤其明显。因此，火后恢复评估应分离草本层与灌木层，并综合土壤指标（如SOC、TN、SWC和BD）进行判断，绿度或总生物量指标可能高估实际恢复程度。

研究结论：火烧严重度通过塑造不同的土壤-植被耦合路径，驱动了亚热带马尾松人工林林下草本层和灌木层的异步恢复。高严重度火烧导致早期土壤有机碳和总氮显著损耗、容重增加、pH和有效磷短暂升高，同时促进了草本层多样性和生物量的迅速扩张，但造成了灌木生物量的持续亏缺。至火后3.5年，中低严重度样地的土壤性质和灌木生物量已接近未燃烧参照，而高严重度样地仍表现出土壤碳库未完全恢复和灌木生物量偏低的结构性缺陷。整合的土壤恢复梯度（Soil PC1）与草本生物量呈负关联而与灌木生物量呈正关联，表明草本扩张主要受益于火烧造成的开放环境，而灌木恢复更依赖于土壤水分、碳氮积累和物理结构的逐步改善。这一异步恢复模式强调，火后生态系统恢复不应被视为单一的绿度或生物量轨迹，而应理解为层特异性的土壤-植被耦合过程，高严重度火烧区域的监测和管理需同时关注地上结构重建与地下土壤功能恢复。