作为集水库、粮库、钱库、碳库于一身的大宝库,森林不仅提供多重生态服务功能,还在缓解气候变化、促进社会经济可持续发展等方面发挥着战略性与基础性的作用
[1]。森林是地球上最大的碳汇,碳储存密度比其他植被更高
[2-3]。鉴于不同森林类型的生态服务功能存在差异,如何针对其碳汇功能开展协同与权衡,对于森林碳库的可持续发展和管理至关重要
[1-2,4]。为明确森林主要生态功能、促进分类经营管理、提升林业效能和科学管理,对不同森林类型土壤碳库进行研究具有重要的意义
[5-7]。
土壤碳库是陆地生态系统中最大的碳库
[8-10]。其构成不仅影响自身的累积与分解过程,还直接作用于全球陆地生态系统的碳平衡,同时影响着土壤质量的变化
[11-12]。土壤碳库是指储存在土壤中的碳,由有机碳和无机碳组成
[13-15]。土壤有机碳的形成抵消了二氧化碳的排放,也促进了土壤结构的形成并提高其稳定性
[16-17],土壤碳库是碳汇最重要的功能结构
[18],土壤结构中不同有机碳库形成过程及影响因素成为全球气候变化研究中的热点之一
[19]。影响土壤碳库的因素可分为气候条件、地形条件、土壤条件以及植被条件。当林分类型改变时,会同步影响土壤温度、水分含量、养分循环以及微生物活动等指标
[20]。因此,研究不同林分类型土壤碳库组成变化及稳定性,对提高其土壤质量有重要价值。
高海拔森林对气候变化更敏感,其土壤碳循环呈现显著特征。研究显示,高海拔区有机碳含量普遍较同纬度低海拔区高,微生物驱动的分解速率呈现随海拔升高指数下降趋势
[21-22],贵州威宁高海拔地区平均海拔2 200 m,土壤碳库相关研究主要集中在草海湿地重金属分布
[23]、土地利用对土壤理化性质影响
[24],土壤酶活性及微生物生物量碳
[25],草海湿地周丛生物碳计量
[26]等方面,而对森林土壤碳库及其组分的研究还存在一定不足。自1985年启动人工造林以来,威宁高海拔地区森林类型主要以云南松针叶纯林、云南松-麻栎针阔混交林以及麻栎-桤木阔叶混交林为主。因此,本研究以该区域的典型的针阔混交林、针叶纯林和常绿阔叶混交林3种典型林分为研究对象,探究不同林分土壤碳库及其稳定性的变化特征,为贵州高海拔地区森林经营管理提供理论依据。
1 研究区概况
研究区位于贵州省西北部典型高海拔地区的威宁县(2 200m)国有沙子坡林场黑石工厂区,样地设在海拔1 800~1 900 m的阔叶混交林(KYH)、针阔混交林(ZKH)和针叶纯林(ZYC)人工林中。研究区内常见优势乔木有云南松(
Pinus yunnanensis)、麻栎(
Quercus acutissima)、白栎(
Quercus fabri)、慈竹(
Bambusa emeiensis);优势灌木有河滩冬青(
Ilex metabaptista)、山胡椒(
Lindera glauca)、中华绣线菊(
Spiraea chinensis),优势草本有荩草(
Arthraxon hispidus)、沿阶草(
Ophiopogon bodinieri)等,详细样地信息见
表1。
2 材料与方法
2.1 样地设置与土壤采样
依据不同经营措施,本研究从森林类型角度出发,评估不同森林类型土壤的物理与生态特征的变化。以相同林龄的针叶林、针阔混交林、阔叶混交林为研究对象,在3种不同森林类型的人工林内分别设置大小为25.8 m×25.8 m的样地,每种类型样地设置3个重复,共9个样地(
表1)。采用五点取样法,每个样地选取具有代表性的采样点,首先去除土壤表面凋落物和腐殖质层,对0~20 cm、20~40 cm、40~60 cm、60~80 cm深度土壤进行取样,且均在晴天进行。将同一土层取样点取得的土壤样品混合,并去除可见的植物细根等杂质,保持4 ℃恒温运回实验室,待自然风干后,测定不同森林类型土壤理化性质等指标。林分类型以树种蓄积量比例确定,针叶纯林划分标准为单个针叶树种蓄积≥90%,阔叶纯林划分标准为单个阔叶树种蓄积≥90%,针阔混交林划分标准为针叶树种或阔叶树种总蓄积占35%~65%,阔叶混交林划分标准为阔叶树种总蓄积>65%(
表1)。
2.2 测试方法
土壤全碳采用TOC元素分析仪测定;土壤pH采用电位法测定;参考焦欢等
[27]分离重组有机碳,研磨过筛后采用重铬酸钾外加热法测定有机碳含量
[28],再用差减法计算轻组有机碳含量。根据粒径大小将土壤区分为颗粒有机碳和矿物结合态有机碳,采用超声波分离法分离砂砾和黏粉粒,再用重铬酸钾-硫酸外加热法分别测定有机碳含量
[29]。土壤可溶性有机碳(DOC)用K
2SO
4溶液提取法测定,土壤有机碳含量(SOC)采用《土壤农化分析》
[30]方法测定,土壤微生物量碳(MBC)采用氯仿熏蒸法
[31]测定,易氧化有机碳(EOC)采用高锰酸钾氧化法测定。
2.3 数据处理
本文以阔叶混交林土壤作为对照,计算土壤碳库质量指标
[32]:
碳库活度(L)=土壤活性有机碳含量/土壤非活性有机碳含量;
碳库活度指数(IL)=样品的碳库活度/参考土壤的碳库活度;
碳库指数(ICP)=样品总碳含量/参考土壤总碳含量;
碳库管理指数(ICM)=ICP×IL×100%。
使用WPS2025对数据进行保存和查看,利用R语言(v4.5.1)tidyverse包进行数据的整理和分析
[33],使用stats包对土壤理化因素和碳库质量指数进行非参数组间差异的Kruskal-Wallis检验,
p值矫正方法为Benjamini-Hochberg,采用ggplot2包绘图
[34]。利用rstatix包anova_test函数进行方差分析,利用lavaan包sem函数构建结构方程模型,利用fitMeasures函数计算结构方程模型参数
[35],利用semPlot包semPaths函数绘制结构方程模型示意图
[36]。
3 结果与分析
3.1 林型和土层对土壤有机碳组分及土壤碳库质量的影响
以森林类型为随机效应因子,以土层为固定效应因子,分别对比土壤有机碳、易氧化有机碳、惰性有机碳、微生物生物量碳和可溶性有机碳的差异。结果表明,森林类型和土层的交互作用对土壤有机碳、易氧化有机碳、惰性有机碳、微生物生物量碳有显著影响。其中,森林类型、土层对易氧化有机碳的单独效应均存在显著性差异,而土层对土壤有机碳、惰性有机碳、微生物生物量碳的单独效应存在显著性差异(
表2),这暗示了土壤有机碳及组分的变异主要来源于土层差异而非森林类型的不同,森林类型对可溶性有机碳的影响存在显著性差异,同时,这种差异还受到土层的影响。此外,森林类型和土层的交互作用对CMI、CPI、LI和L指数的影响不显著,二者的单独效应影响亦不显著(
表2)。
进一步地,本文探究了土层对不同林型土壤有机碳组分含量的影响,结果表明,对于针叶纯林、针阔混交林和阔叶混交林,土层分别单独可以解释约80%、87%和61%的有机碳含量差异且影响显著(
表3),对于阔叶混交林、针阔混交林以及针叶纯林,土层分别可以单独解释约97%、92%和97%的微生物生物量碳含量差异,单独解释约93%、92%和86%的易氧化有机碳含量差异且影响显著(
表3)。综上,土层是土壤有机碳及其组分含量差异主要来源。此外,土层对土壤碳库质量指数变异的单独解释率均较低。
当考虑森林类型对土壤有机碳及其组分的影响时,对于易氧化有机碳,林型分别可以单独解释A层64%、B层65%的含量差异;林型可单独解释92%的可溶性有机碳含量差异;此外,林型对其他指标(有机碳组分、土壤碳库质量指数)的各个土层的单独解释率均较低(
表4)。
3.2 不同林型间土壤有机碳组分及碳库质量变化特征
多重配对比较分析表明,对于针阔混交林以及针叶纯林,A层的总有机碳含量、惰性有机碳含量和易氧化有机碳均显著高于D层(
图1)。在不同森林类型中,易氧化有机碳在不同森林类型中均表现为A层>B层>D层,而在针叶纯林中,C层的易氧化有机碳含量和B、D层间差异不显著,但在阔叶混交林、针阔混交林中,C层易氧化有机碳和B层之间差异均显著,这表明阔叶混交林和针阔混交林中的易氧化有机碳向下周转的速率高于针叶纯林,尤其在40 cm土深以后(
图1)。在不同森林类型中,微生物生物量碳在不同林型中均表现为A层>B层>C层>D层,且A层显著高于B层,B层显著高于D层,阔叶混交林和针叶纯林B层显著高于C层,针阔混交林B层和C层之间差异不显著,这表明阔叶混交林和针叶纯林土壤微生物生物量碳向下周转的速率高于针阔混交林,尤其在20 cm土层以后。碳库质量指数是评估碳库稳定性的常用方法之一,本研究中,森林类型和土层交互作用对碳库质量指数无显著影响,表明森林类型和土层交互作用对碳库稳定性影响不显著。
3.3 不同林型间土壤碳库质量变化特征
针叶树种是本气候区的先锋树种,阔叶混交林是顶极群落类型,林型从针叶纯林到针阔混交林再到阔叶混交林,其实质是主要优势乔木树种的有序替代,即演替阶段的转变。我们定义针叶纯林为演替阶段1,针阔混交林为演替阶段2,阔叶混交林为演替阶段3,以演替阶段为自变量、碳库质量指数为因变量构建回归模型。结果显示,随着正向演替,碳库管理指数显著增加(
P<0.05,
R2=0.18),表现为阔叶混交林(88.73±6.79)>针阔混交林(77.89±9.22)>针叶纯林(59.03±7.41)。表明随着演替过程的进行,碳库稳定性逐渐提升,有机碳更容易被微生物降解和植物利用(
图2)。此外,随着进展演替,不同森林类型间的碳库活度、碳库活度指数和碳库指数均无显著性变化。
3.4 土壤碳库质量影响因子分析
土壤易氧化有机碳是土壤碳库中周转速率较快的组分,是土壤碳库变化的关键动力,土壤易氧化有机碳的积累(EOC)可能是土壤碳库稳定性维持的重要因素。选取土壤碳库中的活跃组分(易氧化有机碳、可溶性有机碳、微生物生物量碳)与土壤碳库管理指数、土壤pH、土壤有机碳(SOC)建立路径分析(卡方5.76,卡方自由度3,
P>0.05,拟合指数0.97,近似均方根误差0.07,见
图3)。随着林型从针叶纯林→针阔混交林→阔叶混交林的转变,碳库管理指数(CMI)逐渐增加,碳库稳定性逐渐提升。在此过程中,土壤有机碳(SOC)和微生物生物量碳(MBC)对碳库管理指数具有显著的负效应,易氧化有机碳(EOC)对土壤碳库管理指数具有显著的正效应,可溶性有机碳(DOC)对土壤碳库管理指数具有正效应,但影响不显著。土壤有机碳和土壤pH对可溶性有机碳具有显著正效应。土壤有机碳和易氧化有机碳对微生物生物量碳具有显著的正效应。
4 讨论
本研究中,不同林型土壤有机碳及其组分、碳库质量指数均无显著性差异,这可能是因为碳输入和输出较为均衡导致的。在广东佛山,不同林分类型凋落物储量表现为针阔混交林>针叶纯林
[37];在安徽大别山,阔叶混交林和针阔混交林的凋落物分解速率高于针叶纯林
[38]。通常来说,阔叶混交林和针阔混交林的凋落物储量以及凋落物分解速率要高于针叶纯林,相应地,土壤碳元素的输入量也有类似变化。然而,贵州高海拔地区阔叶混交林、针阔混交林和针叶纯林土壤养分元素含量无显著性差异,这表明贵州高海拔地区不同森林类型土壤养分元素的输入与分解速率整体平衡,导致土壤养分元素在不同森林类型间无显著性差异。此外,研究区同属同一气候带,地理距离相对较近,该区域温度、降水等气候条件高度相似,可能削弱了植被类型对有机碳分解与积累的差异化影响。
同一林型中,不同土层之间有机碳及其组分差异显著且土层对其具有较好的解释率,推测是因为受到雨水的自然渗透,土壤有机碳及其组分向下周转的结果。然而,不同林分类型土壤碳向下周转的速率存在差异。阔叶混交林和针阔混交林中的易氧化有机碳向下周转的速率高于针叶纯林,尤其在40 cm土深以后。这主要与阔叶树种较高质量的凋落物输入、较深的根系分布以及活跃的微生物群落活动有关
[39],这些因素共同促进易氧化有机碳沿剖面的向下迁移与矿化
[40-41]。此外,阔叶混交林和针叶纯林土壤微生物生物量碳向下周转的速率高于针阔混交林,尤其在20 cm土层以后,这可能源于前者碳源输入的持续性与后者土壤环境异质性较大所致
[42]。
易氧化有机碳是引起土壤碳库容量变化的因子之一,也是土壤有机碳中最为敏感和周转速率最快的组分,可用于表征土壤碳库(如植物群落环境和土壤环境)早期变化,同时也作为土壤肥力变化的指标之一
[43-44]。本研究中,阔叶混交林和针阔混交林土壤易氧化有机碳显著高于针叶纯林,暗示演替后期的土壤易氧化有机碳对植物和土壤环境具有更好的预测性,碳库存储和流转系统更加完善,这和已有研究类似
[45-47]。土壤碳库管理指数表明土壤碳库在稳定性方面有所提高,在针叶纯林→针阔混交林→阔叶混交林演替过程中,土壤碳库管理指数逐渐增加,暗示阔叶混交林土壤碳库稳定性较针阔混交林和针叶纯林高。回归分析结果表明,演替阶段与碳库管理指数之间存在显著的线性关系。然而,方差分析未能检出不同分组均值之间的显著差异。这可能是由于方差分析在将连续自变量离散化时损失了趋势信息,同时组内变异较大、统计功效不足,从而掩盖了整体的线性趋势。因此,结果提示演替阶段对碳库管理指数的作用主要体现为连续梯度效应,而非离散水平之间的显著均值差异。
5 结论
土层是贵州高海拔地区主要森林生态系统土壤有机碳及其组分含量分异的因素而非森林类型。演替后期的林型土壤易氧化有机碳向下周转速率较快。无论是森林类型还是土层,对碳库质量指数的影响均不显著,但随演替阶段推进呈上升趋势,演替后期的阔叶混交林的碳库稳定性最高。土壤易氧化有机碳的积累可能是影响不同森林土壤碳库稳定性的重要因素。
贵州省科技厅计划项目(黔科合基础-ZK[2022]一般241)
2024年森林资源管理补助资金项目(贵州省林业碳汇监测技术研究)
2025年森林资源管理补助资金项目(贵州森林固碳能力及碳汇监测体系建设与研究)