Effect of prescribed burning on soil anti-erodibility and anti-scourability of Pinus yunnanensis forest

研究区位于滇中高原玉溪市新平县南部的照壁山，海拔为1 990~2 050 m，是典型的高原山地丘陵地貌。代表性土壤为红壤，是云南省亚热带北部与亚热带南部的过渡地带。多年平均降水量为1 050 mm，最高气温为32.0 ℃，最低气温为1.0 ℃，年平均气温为15.0 ℃。

在照壁山上对该区域内林分类型、地质地貌及自然条件进行系统调查，选择立地条件相似，平均树龄为30 a且长势较好的云南松林地为研究对象。在防火通道2侧对称布设计划烧除与未烧除样地。林区内灌木主要有木荷Schima superba、杜鹃Rhododendron simsii等，草本主要有紫茎泽兰Ageratina adenophora、贯众Cyrtomium fortunei、毛茛Ranunculus japonicus、白茅Imperata cylindrica等。从20世纪90年代末期起，该区域制订了相对完善的计划烧除制度，并于每年1—2月开展计划烧除工作，至今20余年从未间断。2021年1月实施了该年度森林防火计划烧除工作，计划烧除火焰高度1.5 m，属低强度火烧。在计划烧除和未烧除2块大样地内分别设置20 m×20 m标准样地3块，于2021年7月1—4日采集土壤样品。

按0~10、10~20、20~40和40~60 cm土层分层采样，测定土壤团聚体相关指标，每个样品设置3个重复且均大于500 g。土壤团聚体采用干筛和湿筛法分级，按照粒径＞5.00、5.00~2.00、2.00~1.00、1.00~0.50、0.50~0.25和＜0.25 mm的分级标准分别称量，并计算各粒级质量百分比，土壤抗蚀指数采用静水崩解法测定。

平均质量直径(D_MW)、平均几何直径(D_MG)^[19]计算公式如下：

其中：w_i为第i粒径团聚体组分的质量占土壤总质量的质量百分比(%)；x_i为第i粒径团聚体组分的平均直径(mm)。

在计划烧除和未烧除2块大样地内设置的3块标准样地中，分别设置3个50 cm×50 cm的样方，收集样方中地表凋落物分解层和半分解层，带回实验室即刻称量。采用室内浸泡法，测定凋落物自然含水率、蓄积量以及持水性能^[20]等指标。

土壤抗冲性采用原状土冲刷水槽法测定^[21]。在收集完凋落物的样方中，用自制20 cm×10 cm×5 cm取样器取原状土样，取完后用保鲜膜和泡沫袋将土样包裹，避免运输过程中对原状土样的破坏。实验前将土样浸水24 h至土壤样品饱和，沥去多余重力水后置于1.2 m×0.1 m冲刷槽内。设计冲刷坡度为20°，冲刷流量为2 L·min⁻¹，冲刷时间固定为10 min，隔1 min取1次浑水样，共取10次。试验结束后，将浑水样静置过滤，并烘干称泥沙质量。土壤抗冲系数计算公式如下：C_AS=(Qt)/M。其中：C_AS为土壤抗冲系数(L·g⁻¹)；M为冲走的土样质量(g)；t为冲刷时间(min)；Q为冲刷流量(L·min⁻¹)。

将冲刷试验结束后取样器内剩余土样置于0.25 mm土筛中反复冲洗，洗出根系后于65 ℃条件下烘干，称量，用于计算根长密度[D_R=N/(LWH)]和根系生物量[B_R=M_r/(LWH)]^[22]。其中：N为取样器中根系总长度(cm)，L、W、H分别为取样器的长、宽、高(cm)，M_r为根系干质量(g)。

采用Excel和SPSS 21.0 进行数据统计分析(试验数据均采用上限排除法进行统计)，并用 Origin 9.0 绘制图形。

土壤团聚体质量百分比差异直接影响土壤孔隙度的大小，导致土壤结构和功能发生变化。计划烧除后地表凋落物和土壤环境受到影响，团聚体组分特征发生变化。如图1所示：干筛后各土层土壤团聚体均以粒径＞5.00 mm团聚体占比最大，粒径0.25~0.50 mm团聚体占比最小；随土层深度增加，粒径＞5.00 mm团聚体所占比例减小，粒径＜0.25 mm团聚体占比增大。计划烧除后，0~40 cm土层粒径＞5.00、2.00~5.00 mm团聚体比例增加，粒径＜0.25 mm团聚体比例减小；其中0~10 cm土层变化较明显。未烧除样地0~10 cm土层土壤粒径＞5.00、2.00~5.00和＜0.25 mm团聚体含量分别为24.73%、14.70%和18.13%，计划烧除样地为28.62%、15.95%和13.58%。

Figure 1.  Distribution characteristics of soil aggregates

图1显示：湿筛后各土层土壤团聚体均以粒径＜0.25 mm占比最大。随着土层深度增加，粒径＞5.00和2.00~5.00 mm团聚体所占比例减小，粒径0.25~0.50和＜0.25 mm团聚体占比增大。未烧除样地0~10、10~20和20~40 cm土层土壤粒径＞2.00 mm土壤团聚体含量分别为29.17%、18.57%和16.01%，计划烧除样地为19.79%、16.80%和12.82%。表明计划烧除后林地内土壤0~40 cm土层粒径＞2.00 mm土壤机械稳定性团聚体含量增加，水稳性团聚体含量减少，40~60 cm土层各粒级团聚体含量无明显变化。

在团聚体组分特征的基础上，平均质量直径和平均几何直径也是评价土壤团聚体大小和结构稳定性的重要指标。一般来说平均质量直径和平均几何直径越大，土体结构稳定性越好^[23]。土壤团聚体稳定性特征见表1。

表1表明：随着土层深度增加，平均质量直径和平均几何直径逐渐减小。干筛处理下，计划烧除使0~40 cm土层平均质量直径和0~10 cm土层平均几何直径显著增大(P＜0.05)，40~60 cm土层平均质量直径和平均几何直径无显著变化。而湿筛处理下，计划烧除使0~40 cm土层平均质量直径和0~10 cm土层平均几何直径显著减小(P＜0.05)，40~60 cm土层平均质量直径和平均几何直径无显著差异，说明计划烧除使土壤团聚度降低，土体结构遇水稳定性减弱。

团聚体结构破坏率是表征土壤团聚体稳定性的重要指标。团聚体结构破坏率越小，土壤结构越稳定，抗侵蚀能力越好。土壤团聚体结构破坏率及抗蚀指数见表2。

表2表明：计划烧除后0~10 cm土层团聚体结构破坏率显著增大(P＜0.05)，40~60 cm土层团聚体结构破坏率无显著差异。采用静水崩解法，能较为直观地分析土壤抗水蚀的能力大小。火烧后，0~20 cm土层抗蚀指数显著减小(P＜0.05)，随土层深度的增加，土粒抵御静水崩解的能力有减弱的趋势，但在40~60 cm土层基本无影响。火烧后林内植被覆盖度降低，使土壤性质和功能发生变化，土壤团聚体遇水更容易崩散，表层团聚体结构破坏率显著增大。可见计划烧除使土壤表层抗蚀性减小，随着土层深度的增加这种影响逐渐变小。

林下凋落物的覆盖能很大程度地减缓降雨对土壤的冲刷。计划烧除后林下草本植物以及地表凋落物全部化为灰烬，灌木生长受到影响，地表裸露面积增加。本研究对多年连续计划烧除后地表凋落物蓄积量及持水性能进行测定，结果见表3。

由表3可以看出：经多年连续计划烧除后，林内地表蓄积量仅为3.82 t·hm⁻²，未进行计划烧除的对照样地地表蓄积量为9.35 t·hm⁻²，是计划烧除样地的2.45倍。凋落物组分变化导致凋落物层持水性能也随之改变，未烧除样地最大持水率、最大拦蓄率、有效拦蓄率均高于计划烧除样地，分别是计划烧除样地的1.47、1.90和2.15倍；分层来看，2种样地类型半分解层蓄积量均低于未分解层，但半分解层最大持水率、最大拦蓄率、有效拦蓄率均高于未分解层。表明计划烧除后凋落物蓄积量显著减少，持水能力下降。

在土壤表层，植物根系分布较广，主要来源于灌木和草本植物。计划烧除后土壤根长密度及生物量见表4。表4表明：在表层土壤根系分布中，计划烧除和未烧除样地＜0.50 mm径级的根生物量及根长密度均大于＞2.00、0.50~2.00 mm径级，均以＜0.50 mm径级根系为主。计划烧除后林地＜0.50 mm根系占总生物量百分比显著降低(P＜0.05)。可见低强度地表火对细根影响较大，一些较粗根系能抵御低强度火烧或在火烧后逐渐恢复生长。

细根在土壤缠绕固结方面发挥重要作用。林火干扰后大量细根烧毁，土壤抗冲性能变化。对烧除前后泥沙流失量及抗冲系数进行配对t检验结果见表5。由表5可见：随着冲刷时间的延长，计划烧除和未烧除样地含沙量总体上呈下降趋势；冲刷过程前3 min，计划烧除样地产沙量较大，计划烧除样地产沙量显著大于未烧除样地(P＜0.05)；3 min后的冲刷产沙量均处于低水平平稳状态，无显著差异；烧除前后抗冲系数也有显著差异(P＜0.05)，未烧除样地抗冲系数是计划烧除样地的4.09倍，表明多年连续计划火烧使土壤抗冲性能显著降低。

林下凋落物的覆盖度和土壤内植物根系分布状况与土壤抗冲性能密切相关，土壤抗冲性与凋落物及根系参数相关性见表6。由表6可知：土壤抗冲系数与凋落物和根系指标均呈正相关关系，其中土壤抗冲系数与凋落物蓄积量、最大持水率、最大拦蓄率、有效拦蓄率、根长密度、根系生物量呈极显著正相关(P＜0.01)，与自然含水率呈显著正相关 (P＜0.05)。因此，林内凋落物蓄积量越大，水分含量越高，越能有效减少地表径流的产生，减轻降雨对土壤的直接侵蚀。同时植被根系在土壤中交错穿插，够显著提髙土壤的抗冲性能，减少水土流失。长期计划烧除使凋落物无法蓄积，根系损毁，是导致抗冲性能大幅下降的主要原因。

土壤团聚体是土壤的基本结构单位，对维持土体结构稳定有重要作用，是评价土壤抗蚀性的指标。林火干扰后，燃烧剩余灰烬堵塞土壤孔隙，从而增加土壤大团聚体。本研究显示：土壤表层粒径＞2.00 mm机械稳定性团聚体增加，这与前人研究结果一致^[24]。地表覆盖物的拦蓄作用以及植物根系的胶结作用能有效减缓降雨雨滴击溅造成的土壤侵蚀，提高水稳性团聚体含量，增强土体结构稳定性，减少土壤侵蚀。火烧后地表裸露，不利于大团聚体的形成，土壤抗蚀性能受到影响。不同火烧强度对土壤抗蚀性的影响不相同。大部分研究认为：林火干扰对多种林型的影响趋势基本一致，火烧后表层凋落物和植物细根被烧毁，土壤黏粒、粉粒结构减少，同时土壤微生物活动受到抑制，有机质含量减少，不利于土壤团粒结构的形成。随着火烧强度增加，土壤的抗蚀性下降^[25−26]。也有研究认为：只有高强度的火烧才会使土壤团聚体破坏严重，中低强度火干扰对土壤团聚体影响不大，还能使土壤团聚体周围形成一层疏水薄膜，减缓土壤的润湿速率，提高土壤团聚体的抗水蚀能力^[27]。本研究显示：在进行周期性计划烧除后，虽是低强度的地表火，但土壤表层水稳性团聚体减少，团聚体结构破坏率增大，随着土层深度的增加，影响逐渐减弱，这可能与火烧的频率有关。已有报道显示：低强度火烧形成的隔水层在2 a后消失^[28]。在本研究中长期的林火干扰可能不利于土壤疏水层的形成，并且使土壤的结构和功能发生变化，再加上表层土壤失去凋落物的拦蓄作用，降雨发生时土壤更易遭到侵蚀。

土壤抗冲性变化是多种因素共同作用的结果。林下凋落物对土壤抗冲性的影响分为直接和间接2个部分：直接影响是凋落物层通过对降雨的截持来减轻降雨动能对土壤的侵蚀，间接影响是凋落物在分解过程中，微生物活性增强，产生大量有机质及胶结物质增强土体稳定性，通过改良土壤质地来减缓水流对土壤的冲刷作用^[29−30]。众多研究表明：土壤凋落物覆盖下土壤容重和土壤孔隙度降低^[31]，土壤水文性质得到有效改良，输沙能力随凋落物蓄积量增加而递减^[32]，且凋落物半分解层的水土保持与水源涵养功能优于未分解层。本研究显示：周期性计划烧除将林地内凋落物烧毁，短期内无法大量蓄积，同时火烧后土壤微生物活动受到抑制，不利于凋落物分解及养分归还，凋落物层稀松，导致计划烧除样地地表大面积裸露，难以抵御降雨和水流的直接冲击与破坏，抗冲系数显著增大。这与陈璐等^[33]对大兴安岭火烧迹地土壤抗冲性研究结果基本一致。根系通过对土壤的覆盖与固持，能有效防止土粒崩散，增强土壤结构稳定性，改善土壤理化性质，抑制水土流失^[34−35]，对土壤抗冲性有重要的积极作用。研究表明：林火干扰后，浅层土壤中的细根系燃烧殆尽，部分粗根系也因林火高温难以存活，无法起到固结土壤的作用，特别是重度火烧迹地由于林地裸露和雨水冲击，土壤形成抗水层，表面板结，孔隙度和渗透率降低，土壤结构严重退化，土壤物理性质受到长期影响^[36]。也有研究认为：死亡的根系腐化会产生大量的胶结物质，促进土壤团聚体的形成，增强土壤的抗侵蚀能力^[37]。本研究发现：与高强度森林火灾不同，计划烧除只将地表凋落物、草本植物及部分灌木烧毁，对云南松生长影响不大。低强度火烧后细根和部分粗根系被烧毁，但仍有一部分存活根系和新生根系，对土壤仍然起着一定的固持作用。本研究是在2021年年度计划烧除结束半年后采样，部分灌木和草本植物恢复生长，在一定程度上减轻了火烧对土壤抗冲性的影响。但长期林火干扰使林内灌草覆盖度远低于未烧除样地，地表裸露面积增加，土壤中根系含量减少，根系间网络串联、根土黏结和生物化学作用减弱，微生物对烧死根系分解速率减缓，可能使其无法产生胶结物质。同时，计划烧除导致降雨发生时雨水与土壤直接接触，土壤抗水蚀能力下降，大量表层土壤易被冲走，存在土壤侵蚀风险。

在云南松林地内开展周期性计划烧除后，大量燃烧剩余灰烬堵塞土壤孔隙，使土壤变得紧实。干筛处理下，计划烧除样地粒径＞2.00 mm团聚体高于未烧除样地，团聚体平均质量直径、团聚体平均几何直径也显著增大；湿筛处理下，计划烧除样地粒径＜0.25 mm团聚体较未烧除样地有所增加，团聚体平均质量直径、团聚体平均几何直径、土壤抗蚀指数均显著低于未烧除样地，土壤团聚体结构破坏率增大，土壤抗水蚀性能减弱。随着土层深度的增加，火干扰对土壤抗蚀性的影响逐渐减弱，对深层土壤基本无影响。

高频率、低强度的地表火对林地内凋落物和植被根系有一定的影响，地表凋落物持水性能减弱，土壤中根系减少，土壤失去凋落物的覆盖和植被根系的固持作用，土壤抗冲性能显著低于未烧除样地。