藏北高寒草甸位于青藏高原腹地,平均海拔在4 500 m以上,因独特的地理环境以及生态构造,使其在气候调节、水源涵养、水土保持、生物多样性保护等生态系统服务的维持有着重要影响
[1-2]。近年来,由于人类活动频繁和全球气候变化的共同影响,藏北优良高寒草甸发生了不同程度的退化,退化面积逐年增大,逐步退化为“黑土滩”或裸地,严重威胁到当地的可持续发展
[3-4]。
高寒草甸的退化主要表现为植被退化和土壤退化,其中土壤退化最明显的标志是土壤养分含量降低,而这一现象又直接导致了高寒草甸植被的退化。植被退化也影响着土壤结构和养分含量,退化程度越高,土壤养分流失越明显
[5]。土壤是植物演替的重要驱动力,研究退化草甸土壤理化性质的演变过程,能准确了解高寒草甸生态系统受损的实质
[6]。已有研究表明,全氮和全磷含量随高寒草甸退化程度加重呈下降趋势,表明氮、磷已成为高寒草甸的主要养分限制因子
[7-8]。同时,随着高寒草甸退化程度的加剧,土壤结构也急剧劣化,如土壤体积质量增加,土壤孔隙度以及含水量、养分有效性降低,微生物功能受到抑制等
[9-12],而土壤理化性质的变化导致植被盖度极速下降以及毒草入侵程度扩大,这进一步引发了退化高寒草甸的贫瘠化。尽管国内外学者对高寒草甸退化机理进行了大量研究,但现有研究多聚焦单一退化阶段或孤立指标
[13-14],对退化梯度下土壤理化性质的系统性演变规律仍缺乏深入解析。
本试验以藏北退化高寒草甸土壤为研究对象,研究不同退化阶段高寒草甸土壤理化性质,对土壤化学性质与物理性质关系进行初步分析,探讨高寒草甸退化过程中土壤理化性质影响因素,以期为西藏退化高寒草甸恢复与重建提供理论支撑。
1 材料与方法
1.1 研究区概况
本研究区位于西藏那曲市色尼区那玛切乡(29°54′41″-30°44′49″ N,90°58′32″-92°22′23″ E),海拔高度在4 570 m左右。该地区主要以高寒草甸为主,总面积达到139.1万hm²,可利用的高寒草甸面积为134.62万hm²;高寒草甸退化面积共计49.08万hm²,占草地总面积的35.11%。其中,轻度退化的草地面积为39.46万hm²,占总面积的28.37%;中度退化的草地面积为7.41万hm²,占总面积的5.33%;而重度退化的草地面积为1.97万hm²,占总面积的1.41%。地区气候类型为高原亚寒带半干旱季风性气候,白天气温较高,夜间温度低,温差大。全年四季不分明,主要为冷暖季。大风天气多,风力10~12级。年平均气温-9.2~11.30 ℃,年平均降水量在300.80~592.10 mm,其中6-9月的降水量占全年总降水量的80%以上。土壤类型从表层至深层依次为砂质壤土和砂土,其中表层(0~30 cm)土壤呈灰黑色,属于土壤-植物根系层,为砂质壤土;而30 cm土层以下的土壤主要为粗砂土以及少量碎石。
1.2 样地设置与试验方法
1.2.1 样地设置
综合考虑研究区的植被情况,根据草地退化程度的相关划分评价标准
[15],在研究区内选取均位于东西坡向下坡位的3种不同退化阶段的高寒草甸,分别为轻度退化(lightly degraded,LD)、中度退化(moderately degraded,MD)和重度退化(severely degraded,SD)高寒草甸,计为3个处理,其特征见
表1。
1.2.2 试验方法
分别在轻度退化、中度退化、高度退化高寒草甸内,各设置1个大样方(100 m×100 m);采用对角线取样法,在每个大样方内设置3个0.5 m×0.5 m的小样方。采集土壤前先去除小样方内土壤表层的枯枝落叶,用直径5 cm环刀取0~30 cm土层的原状土壤,重复3次,带回实验室进行土壤体积质量和团聚体、砾石含量的测定。用直径4 cm的土钻在小样方框内随机采集0~30 cm土层土壤样品,重复3次,放入自封袋中及时带回实验室,共计18份土壤样品。随后将18份土壤样品分别分成两部分,一部分新鲜土样于4 ℃保存,用于土壤含水率测定;其余土样风干后,用于测定土壤理化性质。
1.3 土壤理化性质测定
参考《土壤农化分析》
[16]的方法进行土壤化学性质测定,其中有机质含量用K
2CrO
7容量法-外加热法测定,速效氮含量用碱解扩散法测定,速效磷含量测定采用0.5 mol/L NaHCO
3浸取-钼锑抗比色法测定,速效钾含量用乙酸铵提取法测定,全氮含量用凯氏定氮法测定,全磷含量用H
2SO
4-HClO
4联合消解钼锑抗比色法测定,全钾含量用氢氧化钠溶融法测定,pH值采用电位计法测定。参考《土壤理化分析》
[17] 的方法进行土壤物理性质测定,其中土壤含水率采用烘干法测定;土壤体积质量通过环刀法测定。土壤总孔隙度(
Pt)的计算公式为:
,其中
b为土壤体积质量。
将环刀内的土壤过2 mm土壤筛,采用排水法测定砾石含量。使用Mastersizer 2000激光粒度仪测定土壤黏粒(粒径<2 μm)、粉粒(2~<20 μm)、砂粒(200~<2 000 μm)含量,计算土壤粒径分形维数(
D)。用振荡式机械筛分仪(干筛法)测定不同粒径(<0.25,0.25~<0.5,0.5~<1,1~<2,2~<3,3~<5,5~<7,7~<10,≥10 mm)机械稳定性团聚体含量,根据机械稳定性团聚体组成,计算土壤机械稳定性团聚体平均质量直径(MWD)、几何平均直径(GWD)及>0.25 mm团聚体比重(
R0.25);用团粒分析仪(湿筛法)测定不同粒径(<0.25,0.25~<0.5,0.5~<1,1~<2,2~<5,≥5 mm)水稳性团聚体含量,根据水稳性团聚体组成,计算土壤水稳性团聚体MWD、GWD和
R0.25以及团聚体破坏率(PDA)。
D、
R0.25、MWD、GWD和PDA的计算参考文献[
18]进行。
1.4 数据处理
不同退化高寒草甸土壤特征数据采用Excel 2023分析整理,利用SPSS 23.0软件进行单因素方差分析、Pearson相关性分析和主成分分析。用Origin 2018软件绘图。
2 结果与分析
2.1 不同退化阶段高寒草甸的土壤理化性质
表2显示,不同退化阶段高寒草甸的土壤养分含量整体表现为LD>MD>SD,其中LD有机质含量显著高于MD和SD,速效磷、速效氮和速效钾含量均极显著高于MD、SD。不同退化阶段高寒草甸土壤pH值呈现出SD>MD>LD,土壤含水率和总孔隙度均表现为LD>MD>SD,其中总孔隙度在LD与SD和MD之间存在显著差异(
P<0.05);而土壤体积质量和砾石含量则表现为SD>MD>LD,SD和MD的砾石含量以及土壤体积质量与LD有极显著差异(
P<0.01)。
2.2 不同退化阶段高寒草甸土壤粒径和团聚体分布特征
2.2.1 土壤粒径分布特征
土壤粒径分形维数(
D)反映了土壤质地的均匀程度,其值越低,土壤结构稳定性越高
[19]。藏北不同退化高寒草甸土壤粒径分布与分形维数特征见
图1。
图1显示,在不同退化阶段的高寒草甸中,土壤粒径分形维数表现为SD>MD>LD,可知LD的结构稳定性最高。在高寒草甸的不同退化阶段,土壤黏粒含量为9.95%~14.09%,均低于15%,表现为SD>MD>LD;粉粒含量在17.21%~26.11%,表现为SD>MD>LD;而砂粒含量则为59.95%~72.93%,表现为LD>MD>SD。
2.2.2 土壤机械稳定性团聚体分布特征
由
图2可见,在高寒草甸退化过程中,LD土壤机械稳定性团聚体以<0.25 mm为主,其含量为27.98%;0.25~<0.5 mm(14.94%)和0.5~<1 mm团聚体(12.91%)次之,≥10 mm团聚体最少(5.03%)。MD土壤机械稳定性团聚体以<0.25 mm为主,其含量为26.82%;0.25~<0.5 mm(16.79%)和0.5~<1 mm团聚体(14.83%)次之;5~<7 mm团聚体最少(4.98%)。SD土壤机械稳定性团聚体以<0.25 mm为主,其含量为34.94%;0.25~<0.5 mm(14.97%)和0.5~<1 mm团聚体(12.96%)次之;2~<3 mm团聚体最少(1.29%)。退化高寒草甸的土壤机械稳定性团聚体
R0.25表现为MD(73.18%)>LD(72.02%)>SD(65.06%),当高寒草甸由轻度退化到中度退化转变时,
R0.25增加了1.16%;由中度退化到重度退化转变时,
R0.25减少8.12%,可见由中度退化到重度退化是高寒草甸土壤机械稳定性大团聚体向微团聚体转变的主要阶段。
2.2.3 土壤水稳性团聚体分布特征
不同退化阶段高寒草甸各级粒径水稳性团聚体含量的分布特征见
图3。
图3显示,高寒草甸退化过程中,LD土壤水稳性团聚体以<0.25 mm为主,其含量为56.28%;0.5~<1 mm(13.26%)和2~<5 mm团聚体(11.17%)次之;1~<2 mm团聚体最少(6.88%)。MD土壤水稳性团聚体以<0.25 mm为主,其含量为57.12%;0.5~<1 mm(15.27%)和0.25~<0.5 mm团聚体(8.98%)次之;2~<5 mm团聚体最少(4.46%)。SD土壤水稳性团聚体以<0.25 mm为主,其含量为59.86%;0.5~<1 mm(14.75%)和0.25~<0.5 mm团聚体(9.03%)次之;≥5 mm团聚体最少(0.59%)。退化高寒草甸土壤水稳性团聚体比重(
R0.25)表现为LD(43.72%)>MD(42.88%)>SD(40.14%),可知高寒草甸由轻度退化到中度退化转变时,
R0.25减少的幅度较小(0.84%);由中度退化到重度退化转变时,
R0.25减少的幅度较大(2.74%),可知由中度退化到重度退化是高寒草甸土壤水稳性大团聚体向微团聚体转变的主要阶段。
2.2.4 土壤团聚体平均质量直径和几何平均直径
土壤团聚体的平均质量直径(MWD)和几何平均直径(GWD)是衡量土壤团聚体大小分布的关键指标,其值越高,表明土壤团聚体的稳定性越强
[20]。由
表3可见,不同退化阶段高寒草甸土壤机械稳定性团聚体的MWD为1.86~2.25 mm,从LD到MD以及MD到SD,MWD分别下降了5.78%和12.26%; GWD为1.36~1.49 mm,从LD到MD以及MD到SD, GWD分别下降了1.34%和7.48%。水稳定性团聚体MWD为1.02~1.36 mm,从LD到MD以及MD到SD,MWD分别下降了1.47%和23.88%; GWD为0.67~0.76 mm,从LD到MD以及MD到SD,GWD分别下降了2.63%和9.46%。可知随着高寒草甸退化程度的加剧,土壤机械稳定性团聚体的MWD和GWD均高于水稳定性团聚体,且在不同退化阶段中,土壤机械稳定性表现为LD>MD>SD。
2.2.5 团聚体破坏率
不同退化阶段高寒草甸土壤团聚体的破坏率表现为MD(38.02%)>LD(44.14%)>SD(41.21%),由此可知随着高寒草甸退化程度的加剧,团聚体的破坏率总体增大,土壤团聚体稳定性下降,其中MD样地团聚体破坏率最大,表明高寒草甸演替到中度退化阶段时,土壤团聚体稳定性最差。
2.3 不同退化高寒草甸土壤理化性质的相关性分析
不同退化阶段高寒草甸土壤理化性质的相关性分析结果(
表4)表明,土壤总孔隙度与速效钾含量呈极显著正相关(
P<0.01),与土壤有机质、速效氮含量呈显著正相关(
P<0.05);砂粒与速效钾、有机质含量呈极显著正相关(
P<0.01),与全氮、速效磷、速效氮含量呈显著正相关(
P<0.05);土壤黏粒含量与pH值呈显著正相关(
P<0.05);水稳性团聚体几何平均直径与速效氮含量呈显著正相关(
P<0.05)。
2.4 不同退化高寒草甸土壤理化性质的主成分分析
根据3个不同高寒草甸退化梯度进行分组,再基于土壤理化性质共19个指标进行主成分分析,结果见
图4,其中正负箭头方向反映了变量与主成分的正负相关性,线段长度代表变量对主成分的贡献度,而箭头夹角则表明了变量之间的相关性。
图4显示,PC1与PC2的累计贡献率为78.8%,能很好地代表原始数据的大部分信息。在PC1中,速效氮、速效磷、速效钾、有机质等指标载荷分别为0.261,0.252,0.267和0.258,较其他指标高且为正值,可知土壤速效养分对PC1的贡献较大。在PC2中,全钾、黏粒、土壤含水率等的载荷分别为0.389,0.375和0.302,明显高于其他指标,可知土壤全量养分和质地对PC2的贡献较大。PC1贡献率(58.1%)高于PC2(20.7%),说明PC1对原始数据的解释能力更强,因此土壤速效养分含量对退化高寒草甸的影响大。根据LD、MD、SD分布可以看出,LD在第一和第四象限,MD在第一和第三象限,SD在第二和第三象限,可知不同退化阶段高寒草甸样本各项指标存在较大差异。
3 讨 论
退化高寒草甸通过改变植被-土壤反馈机制,导致土壤碳氮流失、pH值升高及土壤元素失衡,形成了恶性循环
[21]。本研究中,随着高寒草甸退化程度的加剧,土壤养分含量也随之降低,这与杨元武等
[22]的研究结果一致,其原因为高寒草甸的退化可能加速了植物残体的分解过程,导致有机养分向无机养分转化,进而引发了土壤养分的流失;此外,由于高寒草甸的退化,地上植被数量锐减,使得地表裸露面积增加,淋溶作用导致土壤养分大量流失。罗亚勇等
[23]发现,在青藏高原东缘的玛曲县,随着高寒草甸退化程度的加剧,土壤速效养分含量呈先增加后降低趋势,这与本研究结果不一致,可能是由于该地区中度退化草甸地下生物量迅速减少,牧草对土壤养分的吸收量也逐渐减少,进而导致了高寒草甸速效养分的增加。
土壤pH值是土壤退化演替的主要标志
[24]。本研究中,随着高寒草甸退化程度的加剧,土壤pH值逐渐增加,可能由于高寒草甸退化导致地表凋落物减少,使地表保水能力减弱,从而使土壤pH升高;此外,藏北地区强烈的阳光辐射,使土壤中的盐分随着水分大量蒸发并向上运输,从而使土壤pH升高,这与孙磊等
[14]的研究结果一致。说明藏北退化高寒草甸有向土壤盐碱化发展的趋势。
土壤物理性质影响着土壤结构、养分有效性以及保水能力,其对维持土壤稳定起着非常重要作用
[25]。本研究中,随着高寒草甸退化程度的加剧,土壤体积质量和砾石含量增加,总孔隙度下降,这与字洪标等
[9]的研究结果一致,其原因为高寒草甸的退化演替导致草地地上生产力降低,而地下根系数量的减少影响了土壤对水分的吸收和利用,使得土壤含水率下降,导致土壤更加紧实,孔隙度降低,体积质量显著增加;另一方面,随着植被覆盖率的减少,土壤蒸发和风蚀、水蚀作用加剧,土壤蓄水能力和通透性也随之下降。
土壤粒径组成决定了土壤养分和水分的转运,高寒草甸退化成为了影响土壤粒径变化的主要因素
[19]。本研究中,高寒草甸土壤退化过程中,土壤黏粒、粉粒含量升高,砂粒含量下降,这与柳妍妍等
[26]研究结果一致。其原因为退化高寒草甸植被盖度降低后,土壤裸露导致其温度升高,土壤呼吸加剧,微生物对有机质的分解加速,土壤有机碳大量流失,团聚体稳定性下降,导致土壤结构被破坏
[27]。此外,土壤粒径分形维数随着高寒草甸退化程度增加而逐渐增大,与宫阿都等
[28]的研究结果一致。高寒草甸植被退化后,土壤裸露在地表,风蚀、水蚀作用加剧,导致土壤结构被破坏,土壤粒径分布不均匀,土壤也随之发生退化,说明高寒草甸土壤退化滞后于草地退化。
土壤团聚体的稳定性是评价土壤结构是否良好的重要指标。本研究中,随着高寒草甸退化程度的加剧,土壤机械稳定性团聚体和水稳定性团聚体的MWD、GMD和
R0.25均呈下降趋势,且退化草地<0.25、0.25~<2 mm团聚体含量始终高于2~<5和≥5 mm团聚体,这与陈帅等
[29]的研究结果一致。可能是因为随着高寒草甸退化程度的加剧,土壤有机质含量逐渐下降,降低了土壤团聚作用,使得大团聚体形成作用减弱,大团聚体不断向小团聚体转化,导致土壤团聚体稳定性逐渐下降。本研究中,土壤机械稳定性团聚体和水稳性团聚体均以<0.25 mm为主,与陈秋捷等
[21]的研究结果并不一致,这可能与高寒草甸的地形、土壤质地及气候环境影响有关。本研究中,退化高寒草甸土壤机械稳定性团聚体的MWD和GWD均大于土壤水稳性团聚,这与江仁涛等
[30]的研究结果一致,说明水蚀作用使大团聚体破碎,平均粒径减小。本研究中,中度退化高寒草甸的土壤团聚体破坏率最大,这与江仁涛等
[30]的研究结果不一致,可能与研究地区退化草地土壤质地存在差异有关。
本研究结果表明,退化高寒草甸土壤总孔隙度与速效钾、有机质、速效氮含量呈显著或极显著正相关,砂粒与速效钾、有机质、全氮、速效磷、速效氮含量呈显著或极显著正相关,土壤黏粒含量与pH值呈显著正相关关系,土壤水稳性团聚体几何平均直径与速效氮呈显著正相关。说明土壤理化性质之间存在复杂的相关性。本研究的主成分分析结果表明,高寒草甸退化主要是由土壤养分和质地共同作用下的协同过程,这与罗亚勇等
[23]的研究结果一致。由本研究结果可知,高寒草甸退化过程中,土壤养分含量急剧下降,同时土壤结构变差,水土保持能力也开始下降,通过人为方式改变土壤养分含量,可以改善退化高寒草甸土壤结构,提高土壤养分,促进高寒草甸草本植物生长,从而改变群落结构和功能
[31]。
4 结 论
随着藏北地区高寒草甸退化程度的加剧,土壤养分含量、土壤含水率、总孔隙度、砂粒含量均有所下降,土壤体积质量、pH值及土壤砾石、黏粒、粉粒含量升高,土壤粒径分形维数增大,土壤团聚体MWD和GMD均呈下降趋势,其中中度退化高寒草甸土壤团聚体的PAD最大,因此中度退化阶段是藏北高寒草甸退化的转折期,同时土壤有机质、速效磷、速效氮和速效钾含量和质地是高寒草甸退化的关键影响因子。
国家重点研发计划项目“青藏高原典型自然保护地生态系统保护恢复及多功能提升技术与示范”(2023YFF1304300)
西藏农牧学院林学博士点(一期)项目(533325001)