马尾松(
Pinus massoniana Lamb.)是中国南方的主要常绿针叶乔木,具有适应性强,耐干旱、瘠薄等优点,在中国林业生产及人工林生态系统中发挥着重要作用
[1]。川南马尾松林分布广、面积大,但森林质量不高
[2],林分单位面积蓄积仅59.04 m
3/hm
2,远低于中国乔木树种的平均值(89.79 m
3/hm
2)
[3],主要原因是过度砍伐、粗放经营或无抚育经营
[4]。马尾松分布于中国16个省区,包括陕西、河南、江西、浙江、福建、四川、贵州、广东和云南等地
[5],是中国分布最广的松属树种。马尾松经济价值高,是中国南方林区主要的用材树种之一
[6],其木材广泛应用于建筑、矿柱、坑木和制浆造纸等,在中国的森林资源总量中占有举足轻重的位置
[7-8]。
土壤线虫是地球上数量最多的后生动物,具有迁移能力弱、世代周期短、功能类群丰富、占据食物网关键链、对土地利用等变化扰动响应敏感等特点
[9],是最常用的土壤质量与功能指示生物之一。一般情况下,将土壤线虫划分为四大营养类群:食细菌线虫(bacterivores,Ba)、食真菌线虫(fungivores,Fu)、植物寄生线虫(plant parasites,Pp)、捕食杂食线虫(predatorsomnivores,Op)
[10]。
土壤线虫不仅能维持土壤生态体系的稳定,还能加快土壤中的物质循环及能量流动
[11]。作为土壤指示动物,其群落特征可以反映出土壤的健康状况及土壤生态系统的稳定性。土壤线虫丰度、多样性等特征都会随着土壤环境的变化而表现出不同的变化趋势,其生态指数能够反映受扰动后或不同生态系统中土壤线虫群落结构和功能的变化,指示土壤有机物降解路径及食物网结构变化等特征
[9]。研究表明,随着种植年限增加,植物寄生线虫将引起连作障碍,对农业可持续发展产生严重的影响。种植年限增加会使土壤的主要分解途径由细菌转化为真菌,且对土壤肥力、理化性质、酶和微生物群落产生负面影响,从而使土壤食物网遭到破坏,土壤健康状况恶化
[12]。因此,本研究选取川南11,19,26,32年生4个林龄马尾松人工林为对象,分析不同林龄林下土壤线虫群落特征及其群落特征与环境因子的关系,明确不同种植年限的马尾松人工林土壤的健康状况,从而对人工林进行更有针对性的管理,实现川南马尾松人工林的可持续发展。
1 材料与方法
1.1 研究区概况
研究区位于长江上游四川省宜宾市高县来复镇森林经营所毛顛坳工区(28.589°-28.603° N,104.558°-104.573° E),地处川云公路中段,地貌以低山丘陵为主。土壤是由砂岩发育而成的山地黄壤,土壤发育完善,层次过渡明显,土层厚度约50 cm。研究区海拔457~533 m,属于亚热带湿润季风气候,年平均气温18.1 ℃,5-10月为雨季,雨量充沛,雨热同季,年降水量1 021 mm。
1.2 样品采集与测定方法
1.2.1 土壤样品采集及土壤线虫分离鉴定
综合考虑地貌和林场不同地形情况,在林场选择11,19,26,32年生4个不同林龄马尾松人工林为研究对象,4个林龄样地基本情况如
表1所示。于2025年1月在4个林龄样地中各设置3个20 cm×20 cm样方,小心去除表层枯枝落叶,按“S”形五点随机取样法确定采样点,小心挖出根系,轻轻抖落大块不含根系的土壤,用无菌刷子采集根部附近4 mm附着的土壤,随机选取100 g后均分为两份,保存于保温箱中,一部分用于测定土壤线虫,另一部分用于测定土壤理化性质。
称取50 g土样,用Baermann湿漏斗法
[13]分离土样中的土壤线虫,分离48 h后,用规格为10 mL的离心管收集分离液,静置约24 h后,去除上清液,保留2 mL液体,水浴锅中60 ℃恒温水浴3 min,再加入6 mL的TAF固定液固定线虫,静置6 h后,小心去除上清液,保留约2 mL液体。在显微镜(OLYMPUSCX23)下鉴定前100条线虫,对所采集到的土壤线虫鉴定至属,随后统计所有线虫数量,依据土壤含水量换算成100 g干土中的线虫数量(条/hg)。
1.2.2 土壤理化性质测定
参照中华人民共和国林业行业标准分别测定土壤含水量、酸碱度以及硝态氮、铵态氮、全磷、有机碳和全氮含量。土壤含水量采用环刀法测定(LY/T 1213-1999);土壤pH采用电位法(LY/T 1239-1999)测定;土壤硝态氮含量采用紫外分光光度法(LY/T 1230-1999)测定;土壤铵态氮含量采用浸提-靛酚蓝比色法(LY/T 1231-1999)测定;土壤全磷含量采用钼锑抗比色法(LY/T 1270-1999)测定;土壤有机碳含量采用重铬酸钾氧化-外加热法(LY/T 1237-1999)测定;土壤全氮含量采用半微量凯式定氮法(LY/T 1228-1999)测定。
1.2.3 土壤线虫指数计算
依据土壤线虫食性可将其分为食细菌线虫(Ba)、食真菌线虫(Fu)、植物寄生线虫(Pp)、捕食杂食线虫(Op)4个营养类群
[14-15];依据土壤线虫r策略到k策略的生活史策略将其分为5个类群,分别赋予功能团值(c-p值),用于评估其在生态系统中的功能
[16]。根据线虫类群及c-p值,计算不同林龄马尾松人工林土壤线虫香农-威纳多样性指数(Shannon-Wiener diversity index,
H')、均匀度指数(evenness index,
J')、优势度指数(Simpson index)、丰富度指数(generic richness index)、营养多样性指数(trophic diversity index,TD)、通路比值(nematode channel ratio,NCR)、自由线虫成熟度指数(maturity index,MI)、植物寄生线虫成熟度指数(plant parasite index,PPI)、线虫富集指数(enrichment index,EI)和结构指数(structure index,SI)
[17]。
1.3 数据分析
使用Excel 2020软件处理原始数据并计算土壤线虫多样性指数和生态功能指数,采用SPSS 26.0对数据进行单因素方差分析(ANOVA)和Duncan’s多重比较,分析不同林龄马尾松样地之间土壤线虫群落的差异。土壤线虫指标与环境因子间的关系采用Pearson’s相关性分析;采用主成分分析(principal component analysis,PCA)对不同样地中的土壤线虫功能类群进行排序。采用Origin 2024绘图。
2 结果与分析
2.1 不同林龄马尾松人工林土壤线虫的群落组成
本研究共分离得到土壤线虫8 351条,隶属于8目30科62属,其中威尔斯属(Wilsonema)在4个林龄中均为优势属,在11,19,26和32年生马尾松人工林总个体数中分别占34.09%,16.72%,11.37%和48.49%。
11年生马尾松人工林土壤中分离得到线虫133条,个体平均密度为115条/hg。其中,优势属为威尔斯属(Wilsonema),占总个体数的34.09%;常见属为穿咽属(Nygolaimus)、滑刃属(Aphelenchoides)、野外垫刃属(Aglenchus)等13属,占总个体数的56.83%;稀有属为克拉尔属(Clarkus)、垫刃属(Tylenchus)、头垫刃属(Cephalenchus)等12属,占个体总数的9.08%。
19年生马尾松人工林土壤中共分离得到线虫963条,个体平均密度为816条/hg。其中,优势属为潜根属(Hirschmanniella)、威尔斯属(Wilsonema)和垫刃属(Tylenchus),占总个体数的47.00%;常见属为针属(Paratylenchus)、盆咽属(Panagrolaimus)、裸矛属(Psilenchus)等17属,占总个体数的46.00%;稀有属为真滑刃属(Aphelenchus)、膜皮属(Diphtherophora)等19属,占总个体数的7.00%。
26年生马尾松人工林土壤中共分离得到线虫1 304条,个体平均密度为1 339条/hg。其中,优势属为针属(Paratylenchus)、潜根属(Hirschmanniella)、威尔斯属(Wilsonema),占总个体数的52.33%;常见属为裸矛属(Psilenchus)、扭沟属(Plectonchus)、肾状属(Rotylenchulus)等14属,占总个体数的40.67%;稀有属为三孔属(Tripyla)、中矛线属(Mesodorylaimus)、拟丽突属(Acrobeloides)等16属,占总个体数的7.00%。
32年生马尾松人工林土壤中共分离得到线虫587条,平均个体密度为512条/hg。其中,优势属为威尔斯属(Wilsonema),占总个数的48.33%;常见属为潜根属(Hirschmanniella)、单宫属(Monhystera)、绕线属(Plectus)等12属,占总个数的41.00%;稀有属为三等齿属(Pelodera)、头垫刃属(Cephalenchus)、真矛线属(Eudorylaimus)等27属,占总个体数的10.51%。
随着马尾松人工林林龄的增加,土壤线虫数量呈先增加后降低的变化趋势。在11年生马尾松人工林中土壤线虫总数量最低,随后在26年生马尾松人工林中土壤线虫数量达最高值,但与19年生马尾松人工林差异并不显著,而在32年生马尾松人工林中线虫数量显著下降(
图1)。
2.2 不同林龄马尾松人工林土壤线虫群落生态指数及区系分析
图2显示,马尾松人工林土壤线虫群落的香农-威纳多样性指数、优势度指数和丰富度指数均随着林龄的增加呈现出先增加后下降的变化趋势,但各林龄间均无显著性差异。
其中,香农‐威纳多样性指数和丰富度指数分别在19和26年生马尾松人工林中最高,分别在32和11年生马尾松人工林中最低;土壤线虫群落优势度指数在32年生马尾松人工林中最低,其他林龄之间相近,表明长期的种植使林地土壤中可利用的有机物和养分减少。同时,林地土壤线虫的均匀度指数随着马尾松人工林种植年限的增加而逐渐降低,但各林龄间均无显著性差异;林地土壤线虫营养多样性在32年生林地中最高,且显著高于11年生和19年生马尾松人工林,但其余3个林地间均无显著差异。所有样地的通路指数在0.80~0.94,且有随着林龄的增加而增大的趋势,但在各林龄林地间均无显著差异。
不同林龄马尾松人工林土壤线虫区系分析结果(
图3)表明,11,19,26年生样地的土壤线虫区系均位于象限B,表明土壤受干扰程度低,食物网成熟,该区域拥有较高的物种丰富性和更健康的生态系统,能够更好地应对环境变化;32年生样地则处于象限D,表明样地土壤受干扰程度最高,食物网退化。
2.3 不同林龄马尾松人工林土壤线虫群落的营养结构特征
根据土壤线虫各功能类群的c-p值计算得出的自由线虫成熟度指数(MI)、植物寄生线虫成熟度指数(PPI)、PPI/MI如
图4所示。
图4显示,PPI和PPI/MI随马尾松人工林林龄的增加而增大,均以32年生马尾松人工林最高,但仅PPI在不同林龄间有显著差异(
P<0.05),表现为32和11年林龄均与其他林龄存在显著差异,而19与26年林龄间无显著差异;MI的变化趋势则相反,在11年生马尾松人工林中最大,但在不同林龄间均无显著差异(
P>0.05)。
本研究结果表明植物寄生线虫(Pp)在4个林龄马尾松人工林土壤线虫中占比均最高,且随着林龄的增加而同步增加,分别占11,19,26,32年生马尾松人工林土壤线虫的61.61%,71.85%,75.27%,77.22%(
图5)。各林龄土壤中食细菌线虫(Ba)和食真菌线虫(Fu)占比随着林龄的增加而减少。Ba土壤线虫在4个林龄中都极少,Op土壤线虫在11和32年生马尾松人工林中占比更高,分别占各林龄土壤线虫的19.83%和15.86%。
不同林龄马尾松人工林中,c-p2类群土壤线虫在4个林龄中表现出绝对优势,分别占11,19,26,32年林龄土壤线虫的70.71%,62.17%,61.20%,71.07%;c-p3类群土壤线虫随着林龄增加呈现先增加后下降的趋势,在19和26年生马尾松人工林中占比更高;c-p4类群土壤线虫变化则与c-p3类群土壤线虫相反,在11和32年生马尾松人工林中占比更高;c-p5、c-p4、c-p1类群土壤线虫在4个林龄中的占比都较少。
2.4 不同林龄马尾松人工林土壤理化性质特征
研究区内土壤理化性质如
表2所示。
表2显示,随林龄增加,马尾松人工林土壤理化性质在26年生林龄时出现明显“拐点”:土壤含水量呈先增加后下降的趋势,并在26年林龄时最大,且显著高于11,32年生人工林(
P<0.05);26年生马尾松人工林土壤pH显著低于11年生;土壤硝态氮、铵态氮含量呈先升后降趋势,在26年生人工林时达峰值,且显著高于其他3个林龄的马尾松人工林;土壤全磷、有机碳、全氮含量均在26年生人工林时骤升至最高,均显著高于32年生马尾松人工林,表明长期经营后土壤肥力可能退化,26年生是马尾松人工林土壤养分与水肥条件最丰沛的临界期。
2.5 马尾松人工林土壤线虫群落与环境因子的相互关系
以不同林龄马尾松人工林土壤线虫营养类群和c-p类群数量为响应变量,以土壤理化性质为解释变量进行主成分分析。结果(
图6)表明,前两轴分别解释了土壤线虫类群的59.64%和15.92%。Ba、Fu、Op线虫营养类群与土壤有机碳、土壤全氮含量呈正相关,Pp线虫营养类群与土壤含水量以及土壤硝态氮、铵态氮、全磷、有机碳含量呈正相关。c-p1、c-p3线虫类群与土壤含水量、铵态氮、全磷含量呈正相关,c-p2线虫类群与土壤有机碳含量呈正相关,c-p4、c-p5土壤线虫类群数量与土壤全氮含量呈正相关。
进一步由土壤线虫群落指标与环境因子相关性分析可知,土壤线虫群落与环境因子密切相关(
图7)。其中,马尾松人工林土壤线虫数量与土壤含水量、硝态氮、有机碳含量呈极显著正相关,与土壤铵态氮、全磷和全氮含量呈显著正相关,与土壤pH呈现显著负相关。土壤线虫类群数量与土壤含水量及土壤硝态氮、铵态氮、全磷、有机碳、全氮含量呈正相关,其中与土壤硝态氮和有机碳含量显著正相关,与土壤pH呈现极显著负相关。土壤线虫香农-威纳多样性指数、均匀度指数、优势度指数与土壤含水量及土壤硝态氮、有机碳、全氮含量呈正相关,与土壤pH均呈负相关,均匀度指数、优势度指数还与土壤全磷和铵态氮含量呈负相关。土壤线虫丰富度指数与土壤pH呈极显著负相关,与土壤有机碳含量呈显著正相关。土壤线虫营养多样性指数、通路比值、自由线虫成熟度指数、植物寄生线虫成熟度指数与土壤pH呈正相关,营养多样性指数、通路比值、植物寄生线虫成熟度指数与土壤pH、有机碳含量和全氮含量呈负相关。
3 讨论
3.1 不同林龄马尾松人工林土壤线虫群落组成及多样性特征
土壤线虫作为指示生物
[18-19],其丰度、多样性等特征都会随着土壤环境的变化而呈现出不同的变化趋势,土壤线虫的生态指数能反映土壤中可利用资源或不同生态系统中土壤线虫群落结构和功能的变化,指示土壤有机物降解路径及食物网结构变化等特征
[12]。
本研究结果表明,林龄对马尾松人工林土壤线虫的个体密度有显著影响,土壤线虫密度随着林龄的增加呈现出先上升后下降的变化趋势,个体密度以26年生马尾松人工林最高,显著高于11和32年生马尾松人工林,该结果与王巧莉等
[20]、王楠等
[21]研究结果一致,表明26年生马尾松人工林土壤可利用有机物等种类更丰富,含量更高,能为多种多样的土壤线虫提供更好的生存条件。本试验共分离得到土壤线虫62属,优势属为威尔斯属、潜根属和针属,4个林龄土壤中线虫的优势种存在差异,这可能与不同林龄马尾松人工林凋落物种类和土壤食物网结构有关
[22-23]。
本研究中,不同林龄马尾松人工林土壤线虫群落的香农-威纳多样性指数、优势度指数和丰富度指数均随着林龄增大呈现出先上升后下降的变化趋势,这与马尾松人工林在发育过程中造成的栖息地异质性有关
[24]。在林龄中期(19和26年)群落物种多样性最高,而至群落演替后期(32年),受优势种影响,群落物种多样性降低。此外,在马尾松林龄中早期,林下植物受强光照、充足水分和较高土壤温度的影响,提高了凋落物分解速率
[25],同时因马尾松树冠尚未完全展开,林分中的喜光灌木和草本迅速萌发,人工林中的植物产生大量有机物,为土壤线虫提供了丰富的食物源,土壤线虫群落的香农-威纳多样性指数、优势度指数和丰富度指数上升。随着林龄的增大,通路比值逐渐增加,表明土壤生态系统的逐渐退化和相对不稳定性
[26]。本试验中,所有林龄样地的通路指数在0.81~0.94,表明细菌是不同林龄马尾松人工林土壤有机质的主要分解者,这与刑树文等
[27]研究结果一致。32年生马尾松人工林中土壤线虫通路比值最大,表明随着林龄的增长,土壤中真菌的生物量或活性相对增加,而细菌的生物量或活性相对减少,通路比值的增加,表明随着马尾松人工林种植年限的增加,土壤食物网中存在更多的营养级和更复杂的相互作用。均匀度指数是反映土壤线虫群落组成结构特征的一个重要指标,较为稳定的群落一般均匀度指数都较高,本试验中土壤线虫均匀度指数随林龄的增加而减小,表明在32年生马尾松人工林中,土壤线虫群落最不稳定。
3.2 不同林龄马尾松人工林土壤线虫群落的营养结构特征
土壤线虫群落的营养结构特征可以反映土壤生态系统的稳定性、土壤食物网结构与功能和土壤肥力与养分状况
[28]。植物寄生线虫成熟度指数(PPI)可作为评价扰动后土壤恢复状况的指标之一。PPI值的大小是由植物组织对植物寄生线虫的承载能力决定的。PPI值越大,表明植被受植物寄生线虫的侵害越严重
[29]。本试验中,马尾松人工林土壤线虫的PPI、PPI/MI随林龄的增加而增大,MI的变化趋势则相反。PPI、PPI/MI以32年生马尾松人工林最高,结果表明随着林龄的增加,土壤线虫中植物寄生线虫也随之增加,种植年限的增加使土壤稳定性减弱,线虫群落多以扩繁速度快的机会主义者为主,该种类的土壤线虫能够寄生于植物的各种组织,使植物发育不足,且能快速适应环境变化并大量繁殖。MI减小也反映了微生物生物量的减少和微生物群落结构的变化,土壤的受干扰程度加大
[30]。这可能与各林龄阶段树木的生长和养分积累及森林的生态系统稳定性有关。26年生的树木生长状况良好,根系发达,为寄生线虫提供了优良的栖息和繁殖环境。该时期的马尾松人工林土壤中有更丰富的有机质和养分,吸引了更多的寄生线虫,促进其繁殖。26年生人工林的生态系统相对稳定,植物寄生性线虫与树木之间的相互作用处于动态平衡。32年生马尾松人工林土壤受干扰程度最大,人工林中树木长时间的竞争导致土壤养分持续被消耗,影响土壤肥力,进而影响土壤微生物活动,加剧了生态系统的脆弱性;另外,32年生人工林中可能存在较多的人为活动,如施肥、病虫害防治等,也可能使土壤质量下降,增加其受干扰程度。
植物寄生线虫在4个林龄马尾松人工林土壤线虫占比最高,且随着林龄的增加占比也相应增加,分别占11,19,26和32年林龄土壤线虫的61.61%,71.85%,75.27%,77.22%;食细菌线虫在4个林龄中占比极少。11和19年生马尾松人工林的植物寄生线虫占比减少,食细菌线虫占比增大,表明11和19年生马尾松人工林土壤线虫营养结构更健康。4个林龄较高的植物寄生线虫占比和极低的食细菌线虫占比,反映出该地区的土壤食物网中以细菌为食物来源的营养通道不畅通,食物网结构趋于简化,生态系统功能受损
[31]。这可能与马尾松人工林种植树种较单一有关。
大量研究表明,较高c-p值的线虫具有较长的生命周期、较低的繁殖率和较强的竞争力,如c-p4和c-p5类群,这类线虫在稳定的生态系统中常表现出优势
[32-33]。较低c-p值的线虫有较短的生命周期、较高的繁殖率和较强的适应性,如c-p1和c-p2类群,在受到干扰的生态系统中表现出优势。本研究的马尾松人工林中c-p2类群土壤线虫在4个林龄中均表现出绝对优势,分别占11,19,26和32年林龄土壤线虫的70.71%,62.17%,61.20%,71.07%,而c-p4、c-p5类群土壤线虫在4个林龄中的占比都较少,表明4个林龄的马尾松人工林都受到不同程度的干扰,生态系统较不稳定
[34],但19和26年生马尾松人工林c-p3类群土壤线虫占比更大,说明两者在4个林龄中相对稳定,受干扰程度更小。
3.3 马尾松人工林土壤线虫群落与环境因子的相互关系
土壤水分和养分含量是影响土壤线虫群落的重要环境因素。本研究中,土壤含水量显著影响了马尾松人工林中土壤线虫的个体数量,这是因为大多数线虫主要依附在土壤孔隙的水膜中,因此线虫密度对土壤含水量变化的响应十分灵敏。在26年生马尾松人工林土壤含水量显著高于11,32年生马尾松人工林,而26年生的土壤线虫个体数量也显著高于其他两个林龄,该结果与吴雨薇等
[35]的研究结果相似。同时,土壤含水量对土壤线虫的生态指数和功能类群也有显著影响。本研究发现,捕食/杂食类线虫数量与土壤含水量呈显著正相关。在含水量更高的林龄土壤中拥有更多的捕食/杂食类土壤线虫,这是因为较高营养级的捕食/杂食类线虫对土壤含水量的变化响应更加灵敏。
森林中土壤铵态氮、硝态氮、有机碳和全氮含量会显著影响土壤线虫的生态指数和功能团。大量研究也表明土壤碳氮资源的变化通过影响微生物群落和植物生长
[36],从而能够改变不同营养类群线虫的食物资源状况,进而对线虫生态指数和功能团产生影响
[12]。本研究中,马尾松人工林土壤食细菌线虫和食真菌线虫在总线虫数中的占比随各林龄土壤中铵态氮含量的增加而减少,该结果与张爱林等
[37]的研究结果一致。过量的铵态氮会使土壤酸化,土壤酸化会抑制土壤微生物的生长和有机质的转化,改变了土壤中的营养物质平衡,进而影响线虫的食物来源和生存环境
[38],使得某些线虫的食物来源减少或质量下降,同时铵态氮的增加还可能改变土壤的物理结构,增加土壤的紧实度,影响线虫的活动和扩散,导致线虫的均匀度下降,反而使优势度增加。本研究发现,c-p1和c-p3线虫类群数量与土壤含水量、铵态氮含量和全磷含量呈显著正相关,26年生马尾松人工林中的土壤含水量、铵态氮含量和全磷含量都显著高于其他3个林龄,表明这类线虫对土壤肥力和水分条件的变化较为敏感
[39];c-p2线虫类群与土壤有机碳量呈正相关,说明土壤有机碳含量的增加有助于c-p2线虫类群的生长和繁殖。
4 结论
不同林龄马尾松木人工林根际土壤线虫群落的研究表明,林龄对土壤线虫个体密度有显著影响,土壤线虫个体密度随林龄的增加而呈现先增加后减少的趋势。林龄会影响林地土壤线虫的功能团和营养类群,这与不同林龄的土壤理化性质有密切的关系,土壤含水量、pH值以及硝态氮、铵态氮含量和有机碳含量均会影响土壤线虫的群落组成。长期种植单一树种会导致人工林生态稳定性下降并影响土壤线虫的生长和分布。
四川省自然科学基金面上项目(2022NSFSC0084)
四川省科技计划项目(2019YJ0427)
国家自然科学基金青年科学基金项目(31400457)