-
森林是陆地生态系统的重要组分部分,土壤有机碳(SOC)储量远高于植被和大气总碳,是全球生态系统碳循环中的活跃碳库,对维持全球碳平衡发挥着不可忽视的作用[1]。全球森林土壤有机碳的微小变化亦会对大气二氧化碳(CO2)变化产生实质影响[2]。根据有机碳的稳定性和分解速率差异可将土壤有机碳分为惰性有机碳、缓效性有机碳和活性有机碳[3]。土壤活性有机碳如易氧化碳(ROC)、颗粒有机碳(POC)等,具有活跃性较高、易被微生物氧化和分解、移动快等特点[4],能显著影响土壤中化学物质的溶解、吸附、吸收、解吸、迁移等,对植物养分供应有最直接的作用[5]。活性有机碳库中,微生物生物量碳(MBC)等对土壤管理、森林经营措施的响应更敏感,是土壤质量的重要指示指标[6]。森林生态系统复杂性高,不同森林类型及同一类型中的不同树种组成、不同演替阶段的土壤有机碳及其活性组分特征亦有明显差异,同时也存在同一林分土壤有机碳随土层深度增加有效性降低的规律[7]。值得注意的是土壤有机碳各组分间能够相互转化,不同组分间又共同调控土壤碳库平衡、碳循环过程,影响和指示森林土壤质量[8]。因此,厘清同一优势树种下不同演替阶段林分或管理措施对森林土壤有机碳活性组分分布积累的影响及其内在机制,对于准确评估土壤质量及碳库动态变化具有重要意义。
松材线虫Bursaphelenchus xylophilus病是一种主要危害马尾松Pinus massoniana林分的毁灭性森林病害。中国感染该病的松林面积已超过180万hm2[9]。马尾松林为南方低山丘陵区群落演替的主要先锋林分,浙江省的松林面积约占全省森林总面积的54.87%。2023年浙江省秋季普查发现:已遭受松材线虫病危害的松林面积约28.72万hm2,病死木达153.2万株[10−11]。当前松材线虫病病害林分的主要清理方式为卫生伐疫木,采伐后出现的林窗空地则让其自然恢复形成马尾松次生林[12]。由于对不同病害程度马尾松林分采伐强度不一致,所以不同采伐强度下自然恢复后的植被类型、组成及林分格局会有很大差异[13]。究其原因为不同采伐强度会导致不同大小的林木生长空间,缓解被伐疫木周围活立木对光照和养分的竞争[14],且自然恢复还使林分中的物种组成发生变化,改善了土壤的理化性质和林分微环境[15]。MA等[16]和DON等[17]研究表明:森林土壤中的碳循环对林业经营管理措施极为敏感,植被恢复能有效提高退化生态系统的土壤有机碳质量分数。然而,由于不同植被恢复措施下的林分物种种类、地上地下生物量、凋落物量均差异显著,因此其对土壤有机碳库的影响亦显著不同[18],且适度采伐可提升森林的土壤有机碳质量分数[16]。因此,系统研究马尾松松材线虫病林分土壤活性有机碳对采伐强度、植被自然恢复的响应规律,有利于深入了解采伐及植被自然修复对土壤活性有机碳的影响以及与林分间的互作机制。
土壤活性有机碳主要由林分中鲜有机质与腐殖质间的过渡性物质分解而成,林分结构、土壤理化性质以及微环境的变化均会对土壤有机碳的空间变异产生影响,增加它的不稳定性,从而快速表征土壤的碳库变化。王晓荣等[19]发现:抚育间伐1 a后,马尾松林土壤活性有机碳空间异质性增加,总有机碳与微生物生物量碳显著降低。翟凯燕等[20]则发现:中度间伐9 a后的马尾松林土壤总有机碳、颗粒有机碳、易氧化有机碳质量分数提升,而重度间伐林分表现为下降。另有学者研究发现:不同植被恢复措施作用后土壤活性有机碳库出现以下3种情况:显著增加[21]、无显著变化[22]、有所减少[23]。可见,采伐调节林分密度对土壤有机碳的影响以及与林分间的互作机制[24],以及植被恢复后不同演替阶段的森林土壤有机碳变化规律须进一步探索[25]。本研究以中度采伐、重度采伐等2种强度采伐后自然恢复5、15 a的马尾松松材线虫病次生林为研究对象,探讨其土壤有机碳及其活性组分的累积变化特征,旨在深入理解森林经营管理对马尾松次生林土壤有机碳及其活性组分的变异特性,探明采伐后自然恢复对土壤活性有机碳的长期效应,准确评估和调控松材线虫病危害后的马尾松次生林土壤碳库情况,为遭受松材线虫病危害的马尾松林经营管理提供科学依据。
-
研究区位于亚热带地区的浙江省杭州市余杭区与临安区内(29°56′14″~30°34′17″N,118°50′22″~120°20′21″E)。该地属亚热带季风气候,四季分明,海拔为90~243 m,年平均气温为15.4~18.4 ℃,平均相对湿度为70%~81%,年降水量为1 454.0~1 882.0 mm,年无霜期为240.0 d,年日照时数为1 352.0~1 765.0 h。
研究林分为马尾松健康林及遭受松材线虫病危害后的马尾松次生林。林分优势层以马尾松为主,另有木荷Schima superba、杉木Cunninghamia lanceolata、青冈Quercus glauca、朴树Celtis sinensis、樟树Camphora officinarum、漆树Toxicodendron vernicifluum等树种。亚优势层主要有檵木Loropetalum chinense、枫香Liquidambar formosana等。灌木层有赤楠Syzygium buxifolium、山矾Symplocos sumuntia、栀子Gardenia jasminoides、小蜡Ligustrum sinense、柃木Eurya japonica、总序桂Phillyrea latifolia、映山红Rhododendron simsii等。草本层有芒萁Dicranopteris pedata、狗脊Woodwardia japonica、红盖鳞毛蕨Dryopteris erythrosora、阔鳞鳞毛蕨Dryopteris championii、蹄盖蕨Athyrium filix-femina、菝葜Smilax china、欧洲蕨Pteridium aquilinum、青绿薹草Carex breviculmis、蓝眼庭菖蒲Sisyrinchium angustifolium、马银花Rhododendron ovatum、土茯苓Smilax glabra、六月雪Serissa japonica、络石Trachelospermum jasminoides、金樱子Rosa laevigata等。林分土壤类型为红壤,土层深度为45~55 cm。
-
采用空间替代时间法,于2023年春季选取地形、海拔、土壤类型、坡向、坡位等立地条件相近的马尾松林为研究对象。设置5个处理:未遭受松材线虫病危害且未采伐的马尾松林为对照(ck);多年累计采伐株数强度为30%~40%(中度采伐),经自然恢复5 a的马尾松次生林(ML5);多年累计采伐株数强度为60%~80%(重度采伐),经自然恢复5 a的马尾松次生林(HL5);多年累计采伐株数强度为30%~40%(中度采伐),经自然恢复15 a的马尾松次生林(ML15);多年累计采伐株数强度为60%~80%(重度采伐),经自然恢复15 a的马尾松次生林(HL15)。上述林分均为20世纪80年代荒山飞播造林培育而成。被采伐的林木均为疫木。采伐后均未对林地开展补植,自然恢复期间林木均未遭受病虫危害,且无其他经营管理措施。
2023年3—5月,进行林地踏查后在每个林分设置3个625 m2(25 m×25 m)的标准样地,共15块样地,进行样地植物种类调查(表1),并于2023年6月进行土壤物理性质调查及样品采集。具体方法:采用对角线法在每个样地的对角线上设置3个土壤采样点,挖取1 m宽的土壤剖面,按照0~10、10~20、20~40 cm的分层方式进行环刀和土壤样品采样,将土壤样品装入自封袋中混匀后,放于低温保温箱中带回实验室测理化性质。共挖取45个土壤剖面,采集135个环刀及135份土壤样品。环刀用于测定土壤容重、孔隙度、含水率等物理性质。土壤样品剔除石粒、根系、凋落物后过2 mm孔径土筛,将样品分为2份,一份新鲜土样置于4 ℃冰箱中用于测定微生物生物量碳(MBN)、水溶性有机碳(WSOC)、硝态氮(${\mathrm{NO}}_3^- $-N)、铵态氮(NH4 +-N),另一份置于通风处自然风干后用于测定土壤有机碳(SOC)、颗粒有机碳(POC)、易氧化碳(ROC)。
表 1 不同采伐强度及自然恢复年限下林分基本概况
Table 1. Basic overview of forests under different logging intensities and natural restoration years
处理 采伐强度 马∶阔数量比 恢复年限/a 主要植物 乔木层 灌草层 ck 未采伐 - - 马尾松、檵木、青冈等 蹄盖蕨、菝葜、欧洲蕨、青绿薹草、蓝眼庭菖蒲等 ML5 中度 7∶3 5 马尾松、木荷、朴树等 芒萁、柃木、山矾、赤楠、马银花等 HL5 重度 2∶8 5 木荷、青冈、马尾松等 赤楠、马银花、山矾、土茯苓、总序桂等 ML15 中度 6∶4 15 马尾松、青冈、樟树、漆树等 栀子、总序桂、映山红、六月雪、络石、金樱子等 HL15 重度 6∶4 15 枫香、漆树、樟树、檵木等 山矾、金樱子、总序桂、菝葜、阔鳞鳞毛蕨等 说明:马指马尾松,阔指阔叶树。−表示无此项。林龄均为40 a。 -
土壤容重、pH、含水率、孔隙度等参照《土壤农化分析》中的方法测定[26];土壤有机碳质量分数采用K2Cr2O7-H2SO4消煮法测定[27];土壤易氧化碳质量分数采用高锰酸钾氧化法测定[28];土壤水溶性有机碳质量分数采用水浸提法测定[29];土壤微生物生物量碳质量分数采用氯仿熏蒸浸提法测定[30];土壤硝态氮质量分数采用紫外分光光度法测定[31];土壤铵态氮质量分数采用靛酚兰比色法测定[32];土壤颗粒有机碳质量分数采用六偏磷酸钠分散法测定[33]。
-
采用Excel 2016、SPSS 27.0、Origin 2021进行数据统计与分析。正态分布检验后,采用双因素方差分析法分析采伐强度与恢复年限对土壤有机碳及其活性组分的影响,并采用最小显著差异法(LSD)进行显著性检验。采用 Pearson 相关分析法分析有机碳及其组分与土壤理化性质的关系(P<0.05)。
-
由图1可知:不同采伐强度与恢复年限的马尾松次生林分土壤有机碳质量分数随土层加深而降低,且各处理0~10 cm土层土壤有机碳质量分数均高于20~40 cm土层(P<0.05)。各林分不同土层土壤有机碳质量分数方面,0~10与20~40 cm的2个土层中HL5相差最小,为9.75 g·kg−1,ML5相差最大,为13.86 g·kg−1,各采伐林分均高于ck的7.32 g·kg−1。且各土层ML15显著高于HL5 (P<0.05)。
图 1 不同采伐强度及自然恢复年限下各土层中土壤有机碳变化
Figure 1. Changes of soil SOC content in different soil layers under different logging intensities and natural restoration years
不同采伐强度林分的土壤有机碳质量分数方面,与ck相比,采伐后各林分不同土层的土壤有机碳质量分数均有所增加,其中ML5、ML15、HL15显著高于ck (P<0.05),且同一恢复年限下中度采伐的林分(ML5、ML15)均高于重度采伐林分(HL5、HL15)。
采伐后不同恢复年限的林分土壤有机碳质量分数方面,中度采伐与重度采伐林分各土层土壤有机碳质量分数均随着恢复年限的增加而上升。0~10 cm土层的ML15土壤有机碳质量分数最高,为(33.53±7.86) g·kg−1,ML5、ML15 、HL15均显著高于ck (P<0.05)。20~40 cm土层的ML15土壤有机碳质量分数显著高于ML5 (P<0.05),重度采伐的HL15显著大于ck。
-
各处理土壤水溶性有机碳质量分数随土层加深而降低(图2A),ML15、HL15林分0~10 cm土层的水溶性有机碳质量分数均显著高于20~40 cm (P<0.05)。同一土层深度下,ML15的水溶性有机碳质量分数最高,显著大于ck、HL5、HL15 (P<0.05),0~10、10~20、20~40 cm各土层的水溶性有机碳质量分数均为ck最低,分别为(163.09±2.33)、(145.57±8.92)、(136.55±9.54) mg·kg−1。不同采伐强度林分的土壤水溶性有机碳质量分数方面,采伐后各林分不同土层均有所增加,0~20 cm土层,ML15、HL15分别显著高于ML5、HL5及ck (P<0.05),且ML5、ML15、HL15均显著高于ck (P <0.05)。采伐后不同恢复年限的林分土壤水溶性有机碳质量分数方面,中度采伐与重度采伐林分各土层均随着恢复年限的增加而上升,最高为0~10 cm土层下ML15,为(417.96±26.13) mg·kg−1,是ck的2.56倍。
图 2 不同采伐强度及自然恢复年限下各土层土壤水溶性有机碳、微生物生物量碳、颗粒有机碳、易氧化碳变化
Figure 2. Changes of soil ROC, POC, MBC, and WSOC content in different soil layers under different logging intensities and natural restoration years
各处理土壤微生物生物量碳随土层加深而降低(图2B),且各林分0~10 cm土层的微生物生物量碳质量分数均显著高于20~40 cm (P<0.05),ML15的各土层微生物生物量碳质量分数差异显著(P<0.05)。各处理不同土层的微生物生物量碳质量分数方面,0~10与20~40 cm 土层中 HL5相差最小,为229.13 mg·kg−1,ML5相差最大,为298.99 mg·kg−1,各采伐林分均高于ck的176.18 mg·kg−1,且ML15显著高于HL5 (P<0.05)。不同采伐强度林分的土壤微生物生物量碳质量分数方面,采伐后各林分不同土层均高于ck,同一恢复年限下各土层的中度采伐林分(ML5、ML15)均高于重度采伐林分(HL5、HL15)。采伐后不同恢复年限的林分土壤微生物生物量碳质量分数方面,恢复年限为15 a的2种不同采伐林分(ML15、HL15)均分别显著高于ML5、HL5 (P<0.05),且中度采伐与重度采伐林分各土层微生物生物量碳质量分数均随着恢复年限的增加而上升,最高为0~10 cm 土层下ML15,为(895.87±158.71) mg·kg−1,其中ML15、HL15均显著高于ck (P<0.05)。
各处理土壤颗粒有机碳质量分数随土层加深而降低(图2C),其中ML5、HL5、ML15、HL15的0~10 cm土层的颗粒有机碳质量分数显著高于20~40 cm (P<0.05)。各林分不同土层颗粒有机碳质量分数方面,ML15均为最高,ck均为最低,ML15显著高于ck、ML5、HL5 (P<0.05),为0~10、10~20、20~40 cm相同土层下ck的3.0、7.4、11.8倍。不同采伐强度林分的土壤颗粒有机碳质量分数方面,采伐后各林分不同土层均有所增加,其中ML15、HL15显著高于ck (P<0.05),且同一恢复年限下重度采伐林分(HL5、HL15)均低于中度采伐林分(ML5、ML15),ML15与HL15差异显著(P<0.05)。采伐后不同恢复年限的林分土壤颗粒有机碳质量分数方面,中度采伐与重度采伐林分各土层颗粒有机碳质量分数均随着恢复年限的增加而上升,最高为0~10 cm 土层下ML15,为(11.74±0.26) g·kg−1,15 a恢复期的林分(ML15、HL15)分别显著高于5 a恢复期的林分(ML5、HL5,P<0.05),且ML15、HL15的颗粒有机碳质量分数至少为0~10、10~20、20~40 cm相同土层下ML5、HL5林分的2.60倍以上。
如图2D所示:各林分土壤易氧化碳质量分数随土层加深而降低,其中HL5、ML15的0~10 cm土层易氧化碳质量分数显著高于20~40 cm (P<0.05), ML15的各土层间差异显著(P<0.05)。各林分不同土层的易氧化碳质量分数方面,ML15的0~10 cm最高,为(4.34±0.03) g·kg−1,且10~20、20~40 cm土层ML15均显著高于其他林分(P<0.05)。不同采伐强度林分的土壤易氧化碳质量分数方面,中度采伐与重度采伐林分各土层均随着恢复年限的增加而上升,其中ML5、ML15、HL15显著高于ck (P<0.05),且同一恢复年限下重度采伐林分(HL5、HL15)均低于中度采伐的林分(ML5、ML15)。采伐后不同恢复年限的林分土壤易氧化碳质量分数方面,各采伐强度的林分土壤易氧化碳质量分数随恢复年限延长而显著增加(P<0.05)。
-
双因素方差检验结果(表2)表明:采伐强度及植被自然恢复年限是土壤有机碳及其活性组分质量分数的显著影响因素(P<0.05),其中,采伐强度对土壤活性有机碳组分质量分数的影响达极显著(P<0.001),植被自然恢复年限则对土壤微生物生物量、水溶性有机碳、易氧化碳、颗粒有机碳有极显著影响(P<0.001)。而采伐强度与植被自然恢复年限的交互作用仅对土壤颗粒有机碳有显著影响(P<0.05)。
表 2 采伐及自然恢复年限对土壤有机碳及其活性组分影响的双因素方差分析
Table 2. Two-way ANOVA of soil organic carbon and its active components under logging and natural restoration
影响因子 采伐强度 恢复年限 采伐强度×恢复年限 F P F P F P 有机碳 11.440 <0.050 6.247 <0.050 0.647 0.426 微生物生物量碳 14.838 <0.001 29.273 <0.001 1.100 0.301 水溶性有机碳 133.725 <0.001 54.664 <0.001 0.223 0.640 易氧化碳 15.781 <0.001 9.498 <0.001 0.346 0.560 颗粒有机碳 18.922 <0.001 76.350 <0.001 4.111 <0.05 -
如表3所示:各林分土壤pH为4.23~4.47,容重为0.96~1.40 g·cm−3,土壤含水率为12.05%~20.00%,土壤孔隙度41.05~52.40 mg·kg−1。硝态氮与铵态氮质量分数呈现为随着土层加深而降低的趋势,且随植被恢复年限的延长而增加,最高均为0~10 cm土层ML15林分。采伐强度对土壤pH、容重、硝态氮有极显著影响(P<0.01),对土壤含水量、铵态氮有显著性影响(P<0.05);恢复年限仅对硝态氮有极显著影响(P<0.01),对土壤孔隙度、铵态氮有显著性影响(P<0.05);采伐强度及恢复年限的交互作用对pH和土壤容重有显著性影响(P<0.05)。
表 3 不同采伐强度及自然恢复年限下林分土壤理化因子及双因素方差分析
Table 3. Soil physical and chemical factors and two-way ANOVA analysis in forest stands under different logging intensities and natural restoration years
处理 土层/cm pH 容重/(g·cm−3) 含水量/% 孔隙度/% 硝态氮/(mg·kg−1) 铵态氮/(mg·kg−1) ck 0~10 4.31±0.04 Aa 1.33±0.20 Aa 12.05±5.84 Aa 42.92±1.36 Ab 2.21±0.21 Ac 2.89±0.29 Ac 10~20 4.30±0.06 Ab 1.40±0.18 Aa 14.63±5.27 Aa 42.53±1.25 Ab 1.52±0.23 Bc 2.33±0.78 Ad 20~40 4.27±0.04 Aa 1.35±0.09 Aa 17.84±4.80 Aa 41.05±2.62 Ab 1.17±0.06 Bb 2.08±0.68 Aa ML5 0~10 4.32±0.06 Aa 0.96±0.10 Bc 20.00±1.53 Aa 52.40±2.72 Aa 3.84±0.54 Ab 5.52±1.41 Abc 10~20 4.33±0.04 Ab 1.05±0.08 Bc 15.33±0.67 Ba 51.16±2.46 Aa 2.79±0.38 Bb 4.58±0.66 Abc 20~40 4.24±0.03 Aa 1.15±0.05 Ab 16.00±0.58 Ba 49.28±0.19 Aa 2.05±0.40 Bab 4.02±0.47 Aa HL5 0~10 4.24±0.01 Aa 1.02±0.12 Ac 12.33±0.33 Aa 45.43±2.93 Aab 2.54±0.31 Ac 4.33±0.50 Abc 10~20 4.34±0.03 Ab 1.10±0.12 Ac 15.67±1.45 Aa 48.65±1.36 Aab 2.14±0.48 ABbc 3.39±1.17 Acd 20~40 4.30±0.09 Aab 1.16±0.08 Aab 10.67±2.60 Aa 43.95±2.49 Aab 1.57±0.25 Bb 2.46±0.99 Aa ML15 0~10 4.32±0.06 Aa 1.29±0.12 Aab 19.03±1.74 Aa 49.21±2.30 Aab 6.68±0.46 Aa 13.71±2.61 Aa 10~20 4.31±0.09 Ab 1.29±0.13 Aab 18.01±1.51 Aa 47.48±1.15 Aab 4.87±0.41 ABa 6.33±0.93 Ba 20~40 4.23±0.07 Aa 1.36±0.04 Aa 18.06±1.32 Aa 48.25±2.20 Aab 3.72±0.93 Ba 4.02±1.88 Ba HL15 0~10 4.31±0.17 Aa 1.07±0.06 Abc 19.11±2.21 Aa 49.92±0.75 Aab 4.53±0.47 Ab 6.27±0.66 Ab 10~20 4.23±0.03 Aa 1.25±0.04 Aab 18.38±0.79 Aa 46.30±2.44 Aab 3.01±0.57 ABb 5.46±0.68 ABab 20~40 4.47±0.09 Aa 1.20±0.19 Aab 17.83±0.48 Aa 45.35±2.98 Aab 2.10±0.50 Bab 3.77±0.71 Ba L - ** ** * - ** * Y - - - - * ** * L×Y - * * - - - - 说明:不同大写字母表示相同采伐强度与自然恢复年限下不同土层土壤理化性质差异显著(P<0.05);不同小写字母表示不同采伐强度与自然恢复年限下同一土层土壤理化性质差异显著 (P<0.05)。L表示采伐强度,Y表示恢复年限。*表示0.05水平下差异显著,**表示0.01水平下差异极显著,-表示无显著影响。 -
相关分析可知(图3):土壤有机碳及其活性组分之间均呈显著正相关关系(P<0.05),其中,土壤易氧化碳与与土壤微生物生物量碳相关系数最高,达0.87。土壤容重与土壤有机碳显著负相关(P<0.05),土壤含水率与土壤微生物生物量碳、土壤水溶性碳、土壤颗粒有机碳间表现出显著正相关(P<0.05);土壤硝态氮、土壤铵态氮与土壤有机碳及其活性组分间亦为显著正相关(P<0.05)。
Effects of natural vegetation restoration after logging on soil organic carbon and its active components in Pinus massoniana secondary forests
-
摘要:
目的 研究遭受松材线虫Bursaphelenchus xylophilus病危害的马尾松Pinus massoniana林土壤有机碳及其活性组分分布特征,揭示采伐后自然恢复对马尾松次生林土壤碳汇的影响机制,为有机碳库管理提供科学支撑。 方法 2023年3—6月在浙江省杭州市余杭区与临安区,采用空间替代时间调查法,以未遭受危害的马尾松林为对照(ck),设置4种植被恢复类型:中度采伐恢复5 a (ML5)、中度采伐恢复15 a (ML15)、重度采伐恢复5 a (HL5)和重度采伐恢复15 a (HL15),比较不同采伐强度采伐后自然恢复下马尾松次生林各土层(0~10、10~20、20~40 cm)土壤有机碳、易氧化碳、颗粒有机碳、微生物生物量碳、水溶性有机碳的变化特征,分析其与土壤理化性质的关系,探究其机制。 结果 ①各处理不同土层土壤有机碳及其活性组分质量分数从大到小依次均为土壤有机碳、易氧化碳、颗粒有机碳、微生物生物量碳、水溶性有机碳,且随着土层加深其质量分数相应降低, ML15的0~10 cm土层中土壤有机碳质量分数最高,为33.53 g·kg−1, ck的20~40 cm土层中水溶性有机碳质量分数最低,为136.55 mg·kg−1。②相同采伐强度下的马尾松次生林土壤各活性有机碳质量分数均表现为随着其自然恢复年限增加而提高,其中各土层间颗粒有机碳质量分数变化最大,15 a的各土层高于5 a各土层的2.6倍以上;③ML15的0~10 cm土层中各活性有机碳组分质量分数均表现为最高,且均显著高于ck (P<0.05);④相关性分析表明:土壤中各活性有机碳组分之间及其与硝态氮、铵态氮均呈显著正相关关系(P<0.05)。 结论 中度采伐对马尾松次生林土壤有机碳及其活性组分具有更显著的促进作用,随自然恢复时间越长越有利于土壤活性有机碳组分质量分数增加,适度采伐马尾松纯林并让其自然恢复可促进土壤碳库积累。图3表3参46 Abstract:Objective The objective is to study the distribution characteristics of soil organic carbon (SOC) and its active components in Pinus massoniana forest affected by pine wilt disease, so as to reveal the influence mechanism of natural vegetation restoration after logging on soil carbon sink in P. massoniana secondary forest and provide scientific support for organic carbon pool management. Method From March to June 2023, in Yuhang District and Lin’an District of Hangzhou City, the undamaged P. massoniana forest was used as the control (ck), four vegetation restoration types were set up by using the space substitution time survey method: moderate logging restoration for 5 a (ML5), moderate logging restoration for 15 a (ML15), heavy logging restoration for 5 a (HL5) and heavy logging restoration for 15 a (HL15) to compare the variation characteristics of SOC, readily oxidizable carbon (ROC), particulate organic carbon (POC), microbial biomass carbon (MBC) and water-soluble organic carbon (WSOC) in different soil layers ( 0−10, 10−20, 20−40 cm ) of P. massoniana secondary forest under natural recovery after different logging intensities, and their relationship with soil physicochemical properties was analyzed to explore their mechanisms. Result (1) The content of SOC and its active components in each soil layer ranking from high to low were SOC, ROC, POC, MBC and WSOC, and the content decreased accordingly with the deepening of soil layer. The highest SOC content was 33.53 g·kg−1 in 0−10 cm soil layer of ML15, and the lowest WSOC content was 136.55 mg·kg−1 in 20−40 cm soil layer of ck. (2) Under the same logging intensity, the content of active organic carbon in soil of P. massoniana secondary forest increased with the increase of natural recovery years, among which the content of POC in each soil layer changed the most, and the content of POC in each soil layer of 15 a was over 2.6 times higher than that in each soil layer of 5 a. (3) The content of active organic carbon component in 0−10 cm soil layer of ML15 was the highest, and was significantly higher than that in ck (P<0.05 ). (4) The results of correlation analysis showed that there were significant positive correlations among active organic carbon components in soil and between active organic carbon components and ${\mathrm{NO}}_3^- $-N and NH4 + -N ( P<0.05 ). Conclusion Moderate logging has a significant promoting effect on SOC and its active components in P. massoniana secondary forest. The longer the natural recovery time of vegetation, the more favorable it is to increase the content of soil active organic carbon components. Moderate logging of P. massoniana pure forests and their natural recovery can promote soil carbon pool accumulation. [Ch, 3 fig. 3 tab. 46 ref.] -
表 1 不同采伐强度及自然恢复年限下林分基本概况
Table 1. Basic overview of forests under different logging intensities and natural restoration years
处理 采伐强度 马∶阔数量比 恢复年限/a 主要植物 乔木层 灌草层 ck 未采伐 - - 马尾松、檵木、青冈等 蹄盖蕨、菝葜、欧洲蕨、青绿薹草、蓝眼庭菖蒲等 ML5 中度 7∶3 5 马尾松、木荷、朴树等 芒萁、柃木、山矾、赤楠、马银花等 HL5 重度 2∶8 5 木荷、青冈、马尾松等 赤楠、马银花、山矾、土茯苓、总序桂等 ML15 中度 6∶4 15 马尾松、青冈、樟树、漆树等 栀子、总序桂、映山红、六月雪、络石、金樱子等 HL15 重度 6∶4 15 枫香、漆树、樟树、檵木等 山矾、金樱子、总序桂、菝葜、阔鳞鳞毛蕨等 说明:马指马尾松,阔指阔叶树。−表示无此项。林龄均为40 a。 表 2 采伐及自然恢复年限对土壤有机碳及其活性组分影响的双因素方差分析
Table 2. Two-way ANOVA of soil organic carbon and its active components under logging and natural restoration
影响因子 采伐强度 恢复年限 采伐强度×恢复年限 F P F P F P 有机碳 11.440 <0.050 6.247 <0.050 0.647 0.426 微生物生物量碳 14.838 <0.001 29.273 <0.001 1.100 0.301 水溶性有机碳 133.725 <0.001 54.664 <0.001 0.223 0.640 易氧化碳 15.781 <0.001 9.498 <0.001 0.346 0.560 颗粒有机碳 18.922 <0.001 76.350 <0.001 4.111 <0.05 表 3 不同采伐强度及自然恢复年限下林分土壤理化因子及双因素方差分析
Table 3. Soil physical and chemical factors and two-way ANOVA analysis in forest stands under different logging intensities and natural restoration years
处理 土层/cm pH 容重/(g·cm−3) 含水量/% 孔隙度/% 硝态氮/(mg·kg−1) 铵态氮/(mg·kg−1) ck 0~10 4.31±0.04 Aa 1.33±0.20 Aa 12.05±5.84 Aa 42.92±1.36 Ab 2.21±0.21 Ac 2.89±0.29 Ac 10~20 4.30±0.06 Ab 1.40±0.18 Aa 14.63±5.27 Aa 42.53±1.25 Ab 1.52±0.23 Bc 2.33±0.78 Ad 20~40 4.27±0.04 Aa 1.35±0.09 Aa 17.84±4.80 Aa 41.05±2.62 Ab 1.17±0.06 Bb 2.08±0.68 Aa ML5 0~10 4.32±0.06 Aa 0.96±0.10 Bc 20.00±1.53 Aa 52.40±2.72 Aa 3.84±0.54 Ab 5.52±1.41 Abc 10~20 4.33±0.04 Ab 1.05±0.08 Bc 15.33±0.67 Ba 51.16±2.46 Aa 2.79±0.38 Bb 4.58±0.66 Abc 20~40 4.24±0.03 Aa 1.15±0.05 Ab 16.00±0.58 Ba 49.28±0.19 Aa 2.05±0.40 Bab 4.02±0.47 Aa HL5 0~10 4.24±0.01 Aa 1.02±0.12 Ac 12.33±0.33 Aa 45.43±2.93 Aab 2.54±0.31 Ac 4.33±0.50 Abc 10~20 4.34±0.03 Ab 1.10±0.12 Ac 15.67±1.45 Aa 48.65±1.36 Aab 2.14±0.48 ABbc 3.39±1.17 Acd 20~40 4.30±0.09 Aab 1.16±0.08 Aab 10.67±2.60 Aa 43.95±2.49 Aab 1.57±0.25 Bb 2.46±0.99 Aa ML15 0~10 4.32±0.06 Aa 1.29±0.12 Aab 19.03±1.74 Aa 49.21±2.30 Aab 6.68±0.46 Aa 13.71±2.61 Aa 10~20 4.31±0.09 Ab 1.29±0.13 Aab 18.01±1.51 Aa 47.48±1.15 Aab 4.87±0.41 ABa 6.33±0.93 Ba 20~40 4.23±0.07 Aa 1.36±0.04 Aa 18.06±1.32 Aa 48.25±2.20 Aab 3.72±0.93 Ba 4.02±1.88 Ba HL15 0~10 4.31±0.17 Aa 1.07±0.06 Abc 19.11±2.21 Aa 49.92±0.75 Aab 4.53±0.47 Ab 6.27±0.66 Ab 10~20 4.23±0.03 Aa 1.25±0.04 Aab 18.38±0.79 Aa 46.30±2.44 Aab 3.01±0.57 ABb 5.46±0.68 ABab 20~40 4.47±0.09 Aa 1.20±0.19 Aab 17.83±0.48 Aa 45.35±2.98 Aab 2.10±0.50 Bab 3.77±0.71 Ba L - ** ** * - ** * Y - - - - * ** * L×Y - * * - - - - 说明:不同大写字母表示相同采伐强度与自然恢复年限下不同土层土壤理化性质差异显著(P<0.05);不同小写字母表示不同采伐强度与自然恢复年限下同一土层土壤理化性质差异显著 (P<0.05)。L表示采伐强度,Y表示恢复年限。*表示0.05水平下差异显著,**表示0.01水平下差异极显著,-表示无显著影响。 -
[1] LAL R. Soil carbon sequestration impacts on global climate change and food security[J]. Science, 2004, 304(5677): 1623 − 1627. [2] 刘世荣, 王晖, 栾军伟. 中国森林土壤碳储量与土壤碳过程研究进展[J]. 生态学报, 2011, 31(19): 5437 − 5448. LIU Shirong, WANG Hui, LUAN Junwei. A review of research progress and future prospective of forest soil carbon stock and soil carbon process in China[J]. Acta Ecologica Sinica, 2011, 31(19): 5437 − 5448. [3] 武燕, 黄青, 刘讯, 等. 西南喀斯特地区马尾松人工林林龄对土壤理化性质的影响[J]. 南京林业大学学报(自然科学版), 2023, 48(3): 99 − 107. WU Yan, HUANG Qing, LIU Xun, et al. Effects of Pinus massoniana plantation age on soil physical and chemical properties in karst area in southwest China[J]. Journal of Nanjing Forestry University (Natural Sciences Edition), 2023, 48(3): 99 − 107. [4] BONGIORNO G, BÜNEMANN E K, OGUEJIOFOR C U, et al. Sensitivity of labile carbon fractions to tillage and organic matter management and their potential as comprehensive soil quality indicators across pedoclimatic conditions in Europe[J]. Ecological Indicators, 2019, 99: 38 − 50. [5] 柳敏, 宇万太, 姜子绍, 等. 土壤活性有机碳[J]. 生态学杂志, 2006, 25(11): 1412 − 1417. LIU Min, YU Wantai, JIANG Zishao, et al. A research review onsoil active organic carbon[J]. Chinese Journal of Ecology, 2006, 25(11): 1412 − 1417. [6] YUAN Guoyin, HUAN Weiwei, SONG Hang, et al. Effects of straw incorporation and potassium fertilizer on crop yields, soil organic carbon, and active carbon in the rice-wheat system [J/OL]. Soil and Tillage Research, 2021, 209 : 104958[2024-02-26]. doi: 10.1016/j.still.2021.104958. [7] 张金硕, 李素艳, 孙向阳, 等. 山东省不同植被类型土壤有机碳及其组分分布特征[J]. 土壤, 2024, 56(2): 350 − 357. ZHANG Jinshuo, LI Suyan, SUN Xiangyang, et al. Characteristics of soil organic carbon and its components under different vegetation types in Shandong Province[J]. Soils, 2024, 56(2): 350 − 357. [8] 窦艳星, 侯琳, 马红红, 等. 间伐对松栎混交林土壤活性有机碳的影响[J]. 中南林业科技大学学报, 2015, 35(5): 64 − 69. DOU Yanxing, HOU Lin, MA Honghong, et al. Effects of forest thinning on soil labile organic carbon in a pine-oak mixed forest[J]. Journal of Central South University of Forestry & Technology, 2015, 35(5): 64 − 69. [9] 邓莉丽, 刘青华, 周志春, 等. 抗松材线虫病马尾松种质资源遗传多样性分析及核心种质构建[J]. 浙江农林大学学报, 2024, 41(1): 67 − 78. DENG Lili, LIU Qinghua, ZHOU Zhichun, et al. Genetic diversity analysis and core collection of pinewood nematodiasis-resistant Pinus massoniana germplasm resources[J]. Journal of Zhejiang A&F University, 2024, 41(1): 67 − 78. [10] 浙江省林业局. 浙江省林业局林业概况[EB/OL]. 2023-07-05[2024-02-28]. http://lyj.zj.gov.cn/col/col1275946/. Department of Forestry of Zhejiang Province. Forestry Overview of Zhejiang Forestry Bureau [EB/OL]. 2023-07-05[2024-02-28]. http://lyj.zj.gov.cn/col/col1275946/. [11] 浙江省林业局. 中国绿色时报: “数字森防”助力浙江高质量防控松材线虫病 [EB/OL]. 2024-01-16[2024-02-28]. http://lyj.zj.gov.cn/art/2024/1/16/art_1275964_59065639.html. Department of Forestry of Zhejiang Province. China Green Times: “Digital Forest Defense” Helps Zhejiang to Prevent and Control Pine Wilt Disease with High Quality [EB/OL]. 2024-01-16[2024-02-28]. http://lyj.zj.gov.cn/art/2024/1/16/art_1275964_59065639.html. [12] 吴敏娟, 尤誉杰, 张晓红, 等. 不同干扰模式对受害马尾松人工纯林林分结构的影响[J]. 应用生态学报, 2019, 30(1): 58 − 66. WU Minjuan, YOU Yujie, ZHANG Xiaohong, et al. Effects of different disturbance patterns on stand structure of infected pure Pinus massoniana plantation[J]. Chinese Journal of Applied Ecology, 2019, 30(1): 58 − 66. [13] 张华锋, 陈思宇, 刘刚, 等. 松材线虫病疫木卫生伐对马尾松纯林林分结构的影响[J]. 浙江农林大学学报, 2020, 37(4): 745 − 751. ZHANG Huafeng, CHEN Siyu, LIU Gang, et al. Effects of sanitation cutting pine wilt diseased trees on the stand structure of pure Pinus massoniana plantation[J]. Journal of Zhejiang A&F University, 2020, 37(4): 745 − 751. [14] JANDL R, LINDNER M, VESTERDAL L, et al. How strongly can forest management influence soil carbon sequestration?[J]. Geoderma, 2007, 137(3/4): 253 − 268. [15] PARIONA W, FREDERICKSEN T S, LICONA J C. Natural regeneration and liberation of timber species in logging gaps in two Bolivian tropical forests[J]. Forest Ecology and Management, 2003, 181(3): 313 − 322. [16] MA Junyong, KANG Fengfeng, CHENG Xiaoqin, et al. Moderate thinning increases soil organic carbon in Larix principis-rupprechtii (Pinaceae) plantations[J]. Geoderma, 2018, 329: 118 − 128. [17] DON A, SCHUMACHER J, FREIBAUER A. Impact of tropical land-use change on soil organic carbon stocks: a meta-analysis[J]. Global Change Biology, 2011, 17(4): 1658 − 1670. [18] ZHANG Yaohua, XU Xianli, LI Zhenwei, et al. Effects of vegetation restoration on soil quality in degraded karst landscapes of southwest China[J]. Science of the Total Environment, 2019, 650(2): 2657 − 2665. [19] 王晓荣, 雷蕾, 曾立雄, 等. 抚育间伐对马尾松林土壤活性有机碳的短期影响[J]. 生态学杂志, 2021, 40(4): 1049 − 1061. WANG Xiaorong, LEI Lei, ZENG Lixiong, et al. Short-term effects of tending thinning on soil labile organic carbon in Pinus massoniana stands[J]. Chinese Journal of Ecology, 2021, 40(4): 1049 − 1061. [20] 翟凯燕, 马婷瑶, 金雪梅, 等. 间伐对马尾松人工林土壤活性有机碳的影响[J]. 生态学杂志, 2017, 36(3): 609 − 615. ZHAI Kaiyan, MA Tingyao, JIN Xuemei, et al. Effects of thinning intensity on soil active organic carbon in Pinus massoniana plantation[J]. Chinese Journal of Ecology, 2017, 36(3): 609 − 615. [21] 赵元, 张伟, 胡培雷, 等. 桂西北喀斯特峰丛洼地不同植被恢复方式下土壤有机碳组分变化特征[J]. 生态学报, 2021, 41(21): 8535 − 8544. ZHAO Yuan, ZHANG Wei, HU Peilei, et al. Responses of soil organic carbon fractions to different vegetation restoration in at typical karst depression[J]. Acta Ecologica Sinica, 2021, 41(21): 8535 − 8544. [22] 李文杰, 张祯皎, 赵雅萍, 等. 刺槐林恢复过程中土壤微生物碳降解酶的变化及与碳库组分的关系[J]. 环境科学, 2022, 43(2): 1050 − 1058. LI Wenjie, ZHANG Zhenjiao, ZHAO Yaping, et al. Changes in soil microbial carbon-degrading enzymes and their relationships with carbon pool components during the restoration process of Robinia pseudoacacia[J]. Environmental Science, 2022, 43(2): 1050 − 1058. [23] 苏静, 赵世伟, 马继东, 等. 宁南黄土丘陵区不同人工植被对土壤碳库的影响[J]. 水土保持研究, 2005, 12(3): 50 − 52, 179. SU Jing, ZHAO Shiwei, MA Jidong, et al. Influence of man-made vegetation on carbon pool in southern Ningxia region in Loess Plateau[J]. Research of Soil and Water Conservation, 2005, 12(3): 50 − 52, 179. [24] DUVENECK M, SCHELLER R, WHITE M. Effects of alternative forest management on biomass and species diversity in the face of climate change in the northern Great Lakes region (USA)[J]. Canadian Journal of Forest Research, 2014, 44(7): 700 − 710. [25] 王伟, 李占斌, 李鹏, 等. 生态建设对坡面土壤有机碳分布的影响[J]. 水土保持研究, 2020, 27(2): 35 − 41. WANG Wei, LI Zhanbin, LI Peng, et al. Effect of ecological construction on soil organic carbon distributionon slope land[J]. Research of Soil and Water Conservation, 2020, 27(2): 35 − 41. [26] 鲍士旦. 土壤农化分析[M]. 北京: 中国农业出版社, 2005. BAO Shidan. Soil Agrochemical Analysis [M]. Beijing: China Agricultural Publishing House, 2005. [27] 李鹏, 刘晓君, 刘苑秋, 等. 红壤侵蚀区不同植被恢复模式土壤碳储量特征及其影响因素[J]. 浙江农林大学学报, 2024, 41(1): 12 − 21. LI Peng, LIU Xiaojun, LIU Yuanqiu, et al. Characteristics and influencing factors of soil carbon stocks in different vegetation restoration models in red soil erosion areas[J]. Journal of Zhejiang A&F University, 2024, 41(1): 12 − 21. [28] BLAIR G J, LEFROY R D B, LISLE L. Soil carbon fractions based on their degree of oxidation, and the development of a carbon management index for agricultural systems[J]. Australian Journal of Agricultural Research, 1995, 46(7): 1459 − 1466. [29] 李伟成, 盛海燕, 蒋跃平, 等. 基塘系统不同竹林土壤CO2通量特征及其影响因子[J]. 林业科学, 2018, 54(8): 13 − 22. LI Weicheng, SHENG Haiyan, JIANG Yueping, et al. Soil CO2 flux and its influence factors of different bamboo plantations in the dike-pond ecosystem[J]. Scientia Silvae Sinicae, 2018, 54(8): 13 − 22. [30] 胡婵娟, 刘国华, 吴雅琼. 土壤微生物生物量及多样性测定方法评述[J]. 生态环境学报, 2011, 20(6/7): 1161 − 1167. HU Chanjuan, LIU Guohua, WU Yaqiong. A review of soil microbial biomass and diversity measurements[J]. Ecology and Environmental Sciences, 2011, 20(6/7): 1161 − 1167. [31] 国家环境保护总局. 水质 硝酸盐氮的测定 紫外分光光度法: HJ/T 346—2007 [J]. 北京: 中国标准出版社, 2007. State Environmental Protection Administration. Water Quality−Determination of Nitrate-nitrogen−ultraviolet Spectrophotometry: HJ/T 346−2007 [J]. Beijing: Standards Press of China, 2007. [32] 陈希, 陈伏生, 叶素琼, 等. 丘陵红壤茶园根际氮磷转化对不同强度酸雨胁迫的响应[J]. 应用生态学报, 2015, 26(1): 1 − 8. CHEN Xi, CHEN Fusheng, YE Suqiong, et al. Responses of rhizosphere nitrogen and phosphorus transformations to different acid rain intensities in a hilly red soil tea plantation[J]. Chinese Journal of Applied Ecology, 2015, 26(1): 1 − 8. [33] GARTEN, C T, POST, HANSON P J, et al. Forest soil carbon inventories and dynamics along an elevation gradient in the southern Appalachian Mountains[J]. Biogeochemistry, 1999, 45(2): 115 − 145. [34] 王昌亮, 王庆成, 张程, 等. 间伐强度对落叶松人工林土壤有机碳的影响[J]. 森林工程, 2015, 31(1): 12 − 16. WANG Changliang, WANG Qingcheng, ZHANG Cheng, et al. Effect of thinning intensity on soil organic carbon in Larix gmelinii plantation[J]. Forest Engineering, 2015, 31(1): 12 − 16. [35] 袁知洋, 邓邦良, 张学玲, 等. 武功山草甸植被小群落土壤活性有机碳与土壤养分的典型相关分析[J]. 中南林业科技大学学报, 2016, 36(2): 84 − 90. YUAN Zhiyang, DENG Bangliang, ZHANG Xueling, et al. Canonical correlation analysis between soil active organic carbon and soil nutrient in different vegetation small communities of Wugong Mountain alpine meadow[J]. Journal of Central South University of Forestry & Technology, 2016, 36(2): 84 − 90. [36] 林鑫宇, 惠昊, 王亚茹, 等. 不同林分类型下土壤活性有机碳含量和分布特征[J]. 安徽农业大学学报, 2021, 48(3): 437 − 443. LIN Xinyu, HUI Hao, WANG Yaru, et al. The content and distribution characteristics of soil active organic carbon under different stand types[J]. Journal of Anhui Agricultural University, 2021, 48(3): 437 − 443. [37] GREGORICH E G, CARTER M R, ANGERS D A, et al. Towards a minimum data set to assess soil organic matter quality in agricultural soils[J]. Canadian Journal of Soil Science, 1994, 74(4): 367 − 385. [38] 尚瑶, 傅民杰, 孙宇贺, 等. 温带阔叶林土壤有机碳及其颗粒组分空间分布特征[J]. 水土保持学报, 2014, 28(5): 176 − 181, 301. SHANG Yao, FU Minjie, SUN Yuhe, et al. Spatial distribution of soil organic carbon and particulate organic carbon in temperate broad-leaved forest in northeast of China[J]. Journal of Soil and Water Conservation, 2014, 28(5): 176 − 181, 301. [39] 韩海荣. 华北落叶松人工林根际与非根际土壤活性有机碳及酶活性对密度调控的响应[D]. 北京: 北京林业大学, 2020. HAN Hairong. Response of Active Organic Carbon and Enzyme Activities to Density Adjustment in Rhizosphere and Non-rhizosphere Soil of Larix principis-rupprechtii Plantation [D]. Beijing: Beijing Forestry University, 2020. [40] GEISEN S, HU Shenran, CRUZ T E E, et al. Protists as catalyzers of microbial litter breakdown and carbon cycling at different temperature regimes[J]. The ISME Journal, 2021, 15(2): 618 − 621. [41] PIAO H C, HONG Y T, YUAN Z Y. Seasonal changes of microbial biomass carbon related to climatic factors in soils from karst areas of southwest China[J]. Biology and Fertility of Soils, 2000, 30(4): 294 − 297. [42] 王斐, 马锐豪, 夏开, 等. 森林转换对土壤活性有机碳组分的影响[J]. 水土保持研究, 2023, 30(1): 233 − 240. WANG Fei, MA Ruihao, XIA Kai, et al. Response of soil labile organic carbon fractions to forest conversions[J]. Journal of Soil and Water Conservation, 2023, 30(1): 233 − 240. [43] 朱浩宇, 王子芳, 陆畅, 等. 缙云山5种植被下土壤活性有机碳及碳库变化特征[J]. 土壤, 2021, 53(2): 354 − 360. ZHU Haoyu, WANG Zifang, LU Chang, et al. Variation characteristics of soil active organic carbon and carbon pools under five vegetation types in Jinyun Mountain[J]. Soil, 2021, 53(2): 354 − 360. [44] 徐雅洁, 郭月峰, 姚云峰, 等. 不同林分配置对土壤水分物理性质的影响[J]. 四川农业大学学报, 2021, 39(3): 370 − 377. XU Yajie, GUO Yuefeng, YAO Yunfeng, et al. Effect of different forest-stand configuration on soil-moisture physical properties[J]. Journal of Sichuan Agricultural University, 2021, 39(3): 370 − 377. [45] 辜翔. 中亚热带植被恢复对土壤有机碳库积累及其稳定性的影响[D]. 长沙: 中南林业科技大学, 2019. GU Xiang. The Effect of Vegetation Restoration on the Accumulation and Stability of Soil Organic Carbon Pool in the Mid-Subtropical of China [D]. Changsha: Central South University of Forestry & Technology, 2019. [46] 袁星明, 朱宁华, 郭耆, 等. 南亚热带不同人工林对土壤理化性质的影响及土壤质量评价[J]. 林业科学研究, 2022, 35(3): 112 − 122. YUAN Xingming, ZHU Ninghua, GUO Qi, et al. Effects of different plantations on soil physical and chemical properties and soil quality evaluation in south subtropical zone[J]. Forest Research, 2022, 35(3): 112 − 122. -
链接本文:
https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.20240264