文献参数: 冯金元, 蒋志成, 蒋志仁, 等.祁连山中段青海云杉林土壤有机碳分布特征及影响因素[J].水土保持通报,2025,45(3):261-267. Citation:Feng Jinyuan, Jiang Zhicheng, Jiang Zhiren, et al. Characteristics and influencing factors of soil organic carbon at Picea crassifolia forests in middle part of Qilian Mountains[J]. Bulletin of Soil and Water Conservation,2025,45(3):261-267.
森林生态系统作为陆地生态系统的重要组成部分,其土壤有机碳储量占到全球碳库的2/3以上,是陆地生物碳库的3倍、大气碳库的2倍
[1],其中森林土壤碳库约占全球土壤有机碳库的70%以上,即使微小的变化也会引起大气二氧化碳浓度的改变,造成全球气候波动
[2]。因此,在当前“双碳”背景下,深入研究森林土壤有机碳分布特征和影响因素作用机制对维持碳循环和碳平衡具有重要意义。目前,国内外学者在这方面已经开展了大量研究,主要从土壤因素、植物因素、气候因素、地形因素等方面来阐明土壤有机碳的分布特征及影响因素。如Chow等
[3]认为,土壤微生物的数量受控于土壤含水量,当土壤含水量过高时会限制微生物的数量,进而增强了残余微生物的活性,促进土壤有机碳的矿化。赵青等
[4]的研究表明土壤因子和林分密度等条件对土壤有机碳的影响要远大于植物群落结构特征对土壤有机碳的影响,海拔的不同会影响土壤有机碳的转换进而直接影响土壤有机碳的分布。Guo等
[5]对美国本土研究发现土壤有机碳含量随海拔增加而减少,而程浩等
[6]对武夷山研究发现土壤有机碳含量随海拔增加而增加。由此可见,不同的区域不同的植被类型,即使是同一植被类型土壤有机碳含量随海拔梯度的变化具有各自的变化规律和不同的影响因素。
祁连山地处中国生态脆弱带,自然条件复杂,水热条件差异明显
[7],其森林土壤碳储量具有很大的空间异质性,碳循环过程抗干扰能力也更加脆弱,如果碳平衡被打破,将难以恢复,这对高寒山地的森林系统将造成极大的影响
[8]。不少学者对祁连山高寒山地的碳平衡、碳循环等方面已有许多研究。如曾立雄等
[9]的研究发现随着土层深度的增加,土壤环境条件与土壤微生物活性在逐渐降低,进而影响不同土层土壤有机碳含量。马剑等
[10]对祁连山土壤有机碳研究中发现,土壤有机碳呈现明显的“表聚性”,主要存储在0—40 cm土层,并且随着土层深度增加呈现为先增后减的趋势。朱猛等
[11]对祁连山草地土壤有机碳特征研究中发现,地形通过改变植被格局和水热条件影响凋落物输入与土壤有机质的分解,间接影响土壤有机碳的积累。虽然对祁连山土壤有机碳的研究已有很多,但侧重点各有不同
[12-14],而且较少关注对土壤有机碳含量影响较大的土壤和林分因子,以及这些因子对高寒山地森林土壤有机碳的影响程度如何。因此,本研究选择分布在祁连山大野口流域不同海拔的青海云杉(
Picea crassifolia)林土壤为研究对象,分析青海云杉林不同土层土壤有机碳随海拔变化形成的差异特征,分析土壤有机碳与土壤pH值、容重、温度、质量含水率、林分郁闭度和林分密度之间的相关关系,进而分析土壤及林分因子对不同海拔土壤有机碳的影响程度。研究成果有助于进一步了解土壤有机碳对环境变化的响应机制,对维持地区土壤碳平衡具有重要意义。
1 材料与方法
1.1 研究区概况
研究区位于祁连山北麓中段祁连山国家公园内的大野口流域,地理坐标为100°17′0″—100°18′30″E,38°32′0″—38°33′30″N,流域面积约73.32 km
2,海拔2 590~4 645 m
[15],年均温度5.4 ℃,日、年平均日较差高于10 ℃,无霜期约100 d,年均降水量约416 mm,主要集中在5—9月,年潜在蒸发量1 011.3 mm,年日照时间1 562.6 h,年均相对湿度为60%
[16]。流域属于典型的大陆性高寒山地气候,地形复杂,海拔高差较大,植被和土壤类型沿海拔呈明显的垂直分布,植被以建群树种青海云杉为主,为纯林,占乔木林的69.67%
[17],土壤类型为山地灰褐土,腐殖质积累与淋溶作用较强,无灰化和黏化现象,肥力较高
[18]。
1.2 样品采集
于2021年7—8月,在祁连山大野口流域选择青海云杉林设置样地,海拔梯度间隔为200 m,每个海拔段设置6个重复样地,样地规格为20 m×20 m,将每个样地分成16个5 m×5 m的小样方,进行每木检尺,测定每个小样方内胸径大于等于1 cm所有单株的胸径、树高和冠幅。采用随机抽样法在每个样地中选取3个小样方,于小样方中心使用鱼眼镜头测定郁闭度。在每个样方的左上角、中心、右下角挖取0.4 m(深)×1.0 m(长)×0.5 m(宽)的土坑,并在去除表面枯落物后,在距地面深10,20,30和40 cm处取土壤剖面样品,使用容积为100 cm
3的环刀重复取土3次,再将一个样方同一土层的土壤均匀混合,土壤温度和土壤水分通过Campbell cr310自动记录仪采集,土壤容重采用环刀法测定,土壤pH值采用浸提电极法测定,土壤有机碳采用重铬酸钾法测定
[19]。样地情况详见
表1。
1.3 数据处理
土壤有机碳含量数据的正态性采用Shapiro-Wilk法检验
[20],各层土壤有机碳在海拔间的差异性采用非参数检验法Kruskal-Wallis法分析,多重比较采用Dunn’s检验法检验。各层土壤有机碳与土壤及林分因子使用Pearson相关性分析,再进行线性拟合。使用置换多因素方差分析,分析土壤及林分因子对不同海拔各层土壤有机碳差异的影响。
2 结果与分析
2.1 不同海拔土壤有机碳分布特征
祁连山中段青海云杉林土壤有机碳分布特征如
图1所示。所有土层的土壤有机碳含量随海拔增加均呈现出先减后增的趋势,海拔2 700 m土壤有机碳含量平均最高,为90.12~127.63 g/kg;海拔3 100 m土壤有机碳含量平均最低,为62.17~70.50 g/kg。各海拔土壤有机碳含量随土层深度增加呈逐渐减小趋势,其变化幅度也随着土层深度的增加逐渐平缓,10 cm土层土壤有机碳含量最高,为70.50~127.63 g/kg,各海拔间变化幅度最大;40 cm土层土壤有机碳含量最低,为62.17~90.12 g/kg,各海拔间变化幅度最小。非参数检验结果显示,海拔2 700 m与海拔3 100 m的土壤有机碳在不同土层都存在显著差异,随着土层深度增加,不同海拔土壤有机碳含量差异性在逐渐减小。
2.2 土壤有机碳与土壤及林分因子之间的关系
青海云杉林不同土层土壤有机碳与土壤及林分因子相关关系详见
表2。10 cm土层土壤有机碳与pH值、容重呈显著负相关,相关系数分别为-0.67和-0.91,与温度、郁闭度、林分密度呈现显著正相关,相关系数分别为0.53,0.21,0.45。20 cm土层土壤有机碳与pH值、容重呈现显著负相关,相关系数分别为-0.76和-0.92,与温度、郁闭度、林分密度呈现显著正相关,相关系数分别为0.65,0.30,0.62。30 cm土层土壤有机碳与pH值、容重呈现显著负相关,相关系数分别为-0.72和-0.79,与温度、郁闭度、林分密度呈现显著正相关,相关系数分别为0.49,0.30,0.65。40 cm土层土壤有机碳与pH值、容重呈现显著负相关,相关系数为-0.70和-0.74,与温度、郁闭度、林分密度呈现显著正相关,相关系数分别为0.42,0.32,0.64。
土壤有机碳与土壤及林分因子的线性回归分析结果如
图2所示。土壤有机碳与容重的线性拟合关系最好,决定系数为0.81(
p<0.01),其次是pH值和土壤温度,决定系数分别为0.57(
p<0.01)和0.32(
p<0.01),土壤有机碳与土壤质量含水率线性拟合关系最差,决定系数为0.00(
p>0.05),表明容重大小对土壤有机碳含量影响比土壤pH值和温度大,而土壤质量含水率对土壤有机碳含量几乎没有影响。尽管郁闭度和林分密度与土壤有机碳存在较明显的相关关系,但线性拟合关系均不显著(
p>0.05)。
2.3 土壤及林分因子对不同海拔土壤有机碳差异的影响
各层土壤有机碳与土壤及林分因子之间的差异性分析结果如
表3所示。由
表3可知,10 cm土层土壤有机碳在海拔间形成的差异与土壤pH值和容重有关,决定系数分别为0.44和0.40,其中pH值和容重是造成10 cm土层土壤有机碳差异的主要因素,容重影响较小;20 cm土层土壤有机碳在海拔间形成的差异与pH值、容重和林分密度有关,决定系数分别为0.46,0.29和0.01,其中pH值是造成20 cm土层土壤有机碳差异的主要因素,容重的影响有所下降,林分密度影响较小,影响有所上升;30和40 cm土层土壤有机碳在海拔间形成的差异与pH值、容重和林分密度有关,决定系数分别为0.50,0.20,0.03,0.56,0.16和0.05,其中pH值仍是造成30和40 cm土层土壤有机碳差异的主要因素,容重的影响有所下降,林分密度影响最小。
3 讨 论
3.1 不同海拔高度土壤有机碳分布特征
本研究结果表明,各土层土壤有机碳随海拔高度呈现为先降后升的变化趋势,表层土壤有机碳变化幅度最大,深层土壤有机碳变化幅度最小,表明输入表层土壤的有机质较多,在较为复杂的外界环境作用下变化幅度较大;输入深层土壤的有机质较少,外界环境的影响因子也随之减少,因此变化幅度较小。各海拔土壤有机碳表现为随土层深度增加而减少,最大值出现在海拔2 700 m,最小值出现在海拔3 100 m,这可能与海拔梯度选取有关。分析认为,土壤有机碳的固存释放与植物输入土壤的有机质类型和数量、微生物数量、土壤水热等条件密切相关
[21],不同海拔易形成差异性较大的植被类型和生物量
[22]。海拔2 700 m水热条件较好,有利于林下灌草等植物生长,植被输入土壤的有机质丰富,微生物活性较高,因此土壤有机碳输入多转化快。随着海拔升高,土壤温度逐渐降低,降雨量增加,林下灌草等植物种间竞争加剧,输入土壤的有机质逐渐减少,在海拔3 100 m土壤有机碳输入少、转化快,形成了较低含量的土壤有机碳。
3.2 土壤因子与林分因子对土壤有机碳的影响
本研究表明,土壤pH值、容重、温度与土壤有机碳的相关关系随土层深度增加表现出先升后降趋势,林分密度和郁闭度表现为随土层深度增加与深层土壤有机碳有较强相关性。线性拟合结果表明影响土壤有机碳的土壤及林分因子作用表现为:容重>pH值>温度>林分密度>郁闭度>质量含水率,林分密度、郁闭度和质量含水率的相关性作用不显著。分析认为,土壤表层凋落物丰富,是动植物生长活动的主要场所,较小的土壤容重具有良好的通气性,有利于土壤微生物分解土壤有机质
[23],随着土层深度的增加,土壤孔隙度减小,向下输入有机质的过程受到制约
[24],同时土壤通气性也变差,使土壤微生物在分解有机质过程中因二氧化碳浓度升高氧气含量降低从而降低了微生物活性
[25]。土壤pH值主要通过影响矿质元素的溶解与迁移,影响微生物生存环境以及参与有机碳循环的酶活性从而制约土壤有机碳含量
[26]。土壤温度仅影响表层土壤有机碳,可能是因为光照提升了表层土壤温度,有利于提高表层土壤微生物的分解效率,而深层土壤由于常年处于较低温度,变化幅度也较低,对深层土壤影响较小。
土壤pH值、容重是形成所有土层土壤有机碳差异的主要影响因子,林分密度次之。同时,随着土层深度增加,在各因子对土壤有机碳差异的影响中pH值的比例在逐渐增加,容重的影响在逐渐减小,表明容重虽然与土壤有机碳密切相关,但pH值才是土壤有机碳在海拔间形成差异的主要影响因子。分析认为,土壤有机碳固存与释放的主要参与者是土壤微生物和土壤酶
[27],土壤pH值和容重都会影响土壤微生物的活性,但其影响作用和途径有所差异,容重主要通过影响通气透水性间接影响土壤微生物活性,而pH值的改变对土壤微生物和土壤酶参与有机碳循环过程的影响更为直接和敏感。林分密度主要从输入土壤有机质的种类和质量方面影响土壤有机碳的分布和特征,深层土壤在水热条件稳定的基础上才进一步显现出林分密度的微弱影响,由于青海云杉林分结构单一,林下其他灌草植被分布稀疏,因此对土壤有机碳的差异影响较小。
4 结 论
祁连山青海云杉林40 cm土层土壤有机碳随海拔增加均呈现为先降后升的趋势,海拔高度2 700 m土壤有机碳含量平均最高,海拔高度3 100 m土壤有机碳含量平均最低。土壤有机碳在不同土层都存在显著差异,随着土层深度增加,不同海拔高度土壤有机碳含量差异性在逐渐减小。土壤有机碳与土壤及林分因子相关关系表现为:容重>pH值>温度>林分密度>郁闭度>质量含水率,土壤容重、pH值和温度与土壤有机碳呈极显著相关(p<0.01),其他因子的相关性不显著(p>0.05)。pH值是形成海拔间土壤有机碳差异的主要影响因子,容重的影响有所下降,林分密度影响最小。
甘肃省陇原青年创新创业人才(个人)项目“开发祁连山国家公园生态旅游资源在乡村振兴中的作用研究”(2023LQGR12)