-
我国二十大报告中指出,要积极稳妥推进碳达峰和碳中和[1]。土壤碳代表着地球陆地生态系统中最大的碳储量[2],其中约50%的有机碳储存在森林土壤中[3],高达787 Pg [4]。森林土壤有机碳主要分布于土层1 m深度以内,约占陆地土壤碳库的40%[4]。有研究表明,如果100 cm深的土壤碳含量增加4‰,就能抵消当前CO2的净排放,从而实现相对“零排放”[5]。因此,研究森林土壤有机碳特征变化是维持碳平衡和调节碳循环的关键一环。
目前,对森林土壤碳含量变化的研究因研究尺度、指标和角度不同而存在差异。研究尺度主要从国家[6-8]和区域[9-10]的宏观尺度、不同林地类型实地取样分析的局部尺度[11-12]展开;研究指标主要包括有机碳的含量、质量分数、组分、密度,以及碳储量、碳氮比等指标[13-18];研究角度主要从土地利用方式[19]、人工处理措施[20]和林地类型[21]等方面展开对森林土壤碳的研究。研究表明,近几十年的气候变化和人类干扰导致自然生态系统中大量土壤碳流失[22],这已成为土壤动态研究中最具有争议的课题[23]。而森林土壤碳变化主要与造林、树种、轮作[3]、土地利用方式、气候[24]、土壤理化性质[25]及地形[15]等因素有关。Mathias等研究表明人工林代替原生林会对土壤碳储量造成负面影响[3];Li等研究发现不同树种造林后土壤碳储量产生明显差异[24];Cools等研究得出欧洲森林树种不同是影响土壤碳氮比的主要因素的结论[26]。由此可见,土壤有机碳含量在不同林地类型变化情况各异。因此,亟需探明林地类型间土壤有机碳含量的影响因素,这对森林土壤碳库变化预测具有重要意义,也可为提高森林土壤固碳能力和养分有效性提供理论参考[27]。
云南省作为桉树引种省份之一[28],森林资源生物种类繁多、植被类型多样[29]。然而,针对桉树引种区土壤有机碳含量变化的研究较少且存在争议[24,30],其影响因素尚未明确。因此,论文基于大量野外采样数据,选取亚热带地区澜沧县外来树种桉树人工林(Eucalyptus robusta Smith plantation forest)、当地特有种思茅松纯林(Pinus kesiya var. langbianensis pure forest)和天然常绿阔叶林(Evergreen broad–leaf forest)为研究对象,通过单因素和Duncan方差分析,揭示不同林地类型在0~20 cm、20~40 cm及40~60 cm土层土壤有机碳含量分布特征及差异性;在回归分析基础上,利用冗余分析(RDA)探讨3种林地类型在0~60 cm土层深度土壤有机碳与土壤理化指标、环境及地形因子的相关性,揭示影响不同林地类型土壤有机碳的主要因素及影响程度,以期为当地双碳目标的实现和生态效益评估提供科学依据。
-
研究区位于云南省普洱市澜沧拉祜族自治县(99°29′~100°35′ E,22°31′~23°16′ N),地处澜沧江以西,总面积为8 807 km2,年平均气温19.1℃,年均降雨量1 626.5 mm,属亚热带山地季风气候(图1)。境内海拔高差达1 936 m,水系纵横交织,均属澜沧江水系,土壤类型以赤红壤和红壤为主。立体气候明显,干湿季分明。植被类型多样,水平地带性植被属季风常绿阔叶林,垂直地带性植被主要有热带季雨林、热带季风常绿阔叶林、以思茅松为代表的暖热性针叶林、半湿润常绿阔叶林和中山湿性常绿阔叶林。其中桉树人工林主要分布于研究区中部,面积达到32 514.68 hm2,占县区总面积的3.73%[31]。
-
在综合考察植被类型和地形地貌分布的基础上,依据代表性和典型性原则,于2021年3月在澜沧县选取勐朗镇、东回乡、酒井乡、惠民乡、谦六乡和糯福乡6个乡镇共设置22块20m × 20 m标准样地,对境内桉树人工林(EP)、思茅松纯林(PK)和常绿阔叶林(EB)3种主要林地类型进行土壤样品采集(图1、表1)。桉树人工林主要分布于研究区中部,参照林样地点均位于桉树人工林附近,采样点土壤类型均为红壤和赤红壤。采样方法为在每块标准样地随机选取3个样点进行土壤剖面调查。土壤剖面分三层,即0~20 cm、20~40 cm和40~60 cm,首先,用容积100 cm3环刀于每个土层中部取样,当天环刀称量质量并密封保存,用于测定土壤密度和土壤孔隙度。同时,在每个环刀采样底相近位置另取土样装入有盖铝盒以供土壤含水量的测定。除此之外,在各土层中部用四分法另取1 kg左右的土壤混合均匀后装入布袋,0~20 cm土层样品以“S”形布点法采集后同样用四分法混合,以避免人为扰动。将布袋土壤样品带回室内晾干后供土壤养分的测定。采样时间避开了当地的雨季,减少了降水对样品的影响。
林地类型
Woodland
type样地编号
Sample site
number样地数
Number of
sample plots/plot海拔
Altitude/m年均相对湿度
Annual average
relative humidity/%年均温度
Average annual
temperature/
℃年降水量
Annual
rainfall/
mm坡度
Slope/(°)坡向
Slope
direction桉树人工林
Eucalyptus robusta Smith plantation forestEP 7 1 236~1 617 72~79 18~20 75~127 7~23 阳坡
阴坡
半阴坡思茅松纯林
Pinus kesiya var. langbianensis pure forestPK 8 1 279~1 693 72~79 18~19 76~128 3~29 阴坡
半阳坡常绿阔叶林
Evergreen broad–leaf forestEB 7 1 344~1 659 72~79 18~20 75~128 11~24 阳坡
半阳坡
半阴坡Table 1. Basic situation of sample plot
-
土壤测定参照国家标准执行。土壤密度依据LY/T 1215—1999,土壤含水量依据重量法HJ 613—2011,土粒密度依据LY/T 1224—1999。pH值测定参照NY/T 1377—2007,有机质参照NY/T 1121.6—2006;水解性氮参照LY/T 1228—2015中4,有效磷参照NY/T 1121.7—2014,速效钾参照NY/T 889—2004,全氮参照LY/T 1228—2015,全磷采用NY/T 88—1988,全钾参照NY/T 87—1988。所有指标由具备测试资质的公司实验室测定。
土壤有机碳密度计算公式[16]为:
SOCDi = SOCi × BDi × Di
式中:SOCDi和SOCi分别代表第i层土壤有机碳密度/(kg·m−2)和含量/(g·kg−1);BDi为第i层的土壤密度/(g·cm−3);Di为第i层的土层厚度/m,实验设置土层厚度均为0.2 m。
-
所有数据经Excel 2016预处理后,运用SPSS 26.0进行单因素方差分析(one-way ANOVA)及Duncan法显著性检验(α=0.05),采用Canoco 5.0进行前向选择和冗余分析,采用ArcMap 10.5和Origin 2022进行作图。
-
从水平层面看,0~60 cm土层,EB分别与EP和PK土壤有机碳含量具有显著差异性(P<0.05)。3种林地类型土壤有机碳含量变化范围为14.89~22.05 g·kg−1 (表2),有机碳含量表现为EB(22.05 g·kg−1 )> EP(17.72 g·kg−1 )> PK(14.89 g·kg−1 )。其中,EB土壤有机碳含量分别高出EP和PK土壤有机碳含量的24%和48%。3种林地类型土壤有机碳含量空间变异系数介于55.58%~72.97%,均属中等程度变异[32],表现为EP(72.97%)>PK(62.87%)>EB(55.58%),说明EP林下土壤有机碳空间异质性较强。经偏度、峰度计算和K-S检验,P值接近0.05,但在峰度和偏度绝对值分别小于10和3的情况下,3种林地类型土壤有机碳含量基本服从正态分布。
林地类型
Woodland type平均值
Mean value/
( g·kg−1 )标准差
Standard
deviation最大值
Max value/
( g·kg−1 )最小值
Min. value/
( g·kg−1 )变异系数
Variation
Coefficient /%峰度
Kurtosis偏度
SkewnessPK-S EP 17.72b 9.49 25.63 7.35 72.97 0.81 1.03 0.074 PK 14.89b 9.36 41.54 2.86 62.87 0.14 0.83 0.000 EB 22.05a 12.25 55.28 4.29 55.58 0.21 0.88 0.005 注:表中不同小写字母表示不同林地类型土壤有机碳含量差异显著(P<0.05)
Notes: different letters indicate significant differences in soil organic carbon content of varied vegetation types (P<0.05)Table 2. Statistical characteristics of soil organic carbon content in different forest land types
从垂直剖面看,3种林地类型土壤有机碳含量垂直分布特征明显(图2a)。EP在0~20 cm土层分别与20~40 cm和40~60 cm土层土壤有机碳含量具有显著差异性,而PK和EB在3种土层间均具有显著差异性。分土层看,0~20 cm土层,EB分别与EP和PK土壤有机碳含量具有显著差异性,而EP和PK间不具有显著差异性;20~40 cm土层,仅PK与EB土壤有机碳含量具有显著差异性,而EP与这两种林地类型均差异不显著;40~60 cm土层,EP仅与PK土壤有机碳含量差异性显著,而EB与这两种林地类型均差异不显著。3种林地类型土壤有机碳含量在0~20 cm和20~40 cm土层均表现为EB> EP> PK,在40~60 cm土层表现为EP> EB> PK。0~20 cm土层土壤有机碳含量分别占整个采样土壤剖面的40%、50%和53%,占比最高,呈现明显的表层聚集特征,以EB最为显著。随着土壤深度增加,3种林地类型土壤有机碳含量均呈下降趋势,从表层0~20 cm到中层20~40 cm下降速度较快,而中层20~40 cm到深层40~60 cm下降速度变缓,以EB下降最快,其表层土壤有机碳含量约是深层的3.15倍,极差达23.88 g·kg−1 ;而EP下降最慢,表层和深层土壤有机碳含量的极差仅为12.29 g·kg−1 。
-
从水平层面看,0~60 cm土层,3种林地类型土壤有机碳密度整体表现为EB(平均4.12 kg·m−2)>EP(平均4.09 kg·m−2)>PK(平均3.42 kg·m−2),EB有机碳密度均值分别高出EP 0.73%和PK 20.47%(图2b)。从垂直剖面看,EP土壤有机碳密度在0~20 cm分别与20~40 cm和40~60 cm土层具有显著差异性,而20~40 cm和40~60 cm土层之间不具有显著差异性。PK和EB土壤有机碳密度3个土层间均具有显著性差异(P<0.05),表明3种林地类型土壤有机碳密度垂直分布差异较大。分土层看,0~20 cm土层,PK与EB土壤有机碳密度具有显著差异性,而EP与其他两种林地类型差异均不显著;20~40 cm土层,3种林地类型之间土壤有机碳密度差异均不显著;40~60 cm土层,EP与PK土壤有机碳密度差异性显著,而EB与其他两种林地类型差异均不显著(P<0.05),表明同一土层之间各林地类型土壤有机碳密度差异较小。在0~20 cm土层,3种林地类型土壤有机碳密度与有机碳含量表现出一致的表层聚集现象,分别占整个土壤剖面的46%、49%、50%。随着土壤深度增加,3种林地类型土壤有机碳密度均逐渐降低。其中,以EB下降幅度最大,表层约是深层的2.63倍,极差达3.85 kg·m−2;而EP变化幅度最小,极差为2.76 kg·m−2。
-
从水平层面看,0~60 cm土层,3种林地类型土壤C/N整体表现为EP(平均13.82)>EB(平均13.47)>PK(平均12.99),以EP最高,分别高出EB和PK 2.60%和6.39%。而PK土壤C/N最低,表明思茅松林土壤有机质矿化速率高,可为土壤提供氮素的能力优于其他两种林地类型,而土壤碳储量低于其他两种林地类型。从垂直剖面看,PK不同土层间土壤C/N差异显著;EP和EB土壤C/N在40~60 cm土层分别与0~20 cm和20~40 cm土层差异显著,而0~20 cm和20~40 cm土层之间不具有显著差异性,表明PK土壤C/N垂直分异规律较明显。分土层看,0~20 cm和20~40 cm土层,3种林地类型土壤C/N差异均不显著;40~60 cm土层,EP和PK 土壤C/N差异显著,但均与EB不具有显著差异性,表明两种人工林深层土壤C/N差异性较大。随着土壤深度增加,3种林地类型土壤C/N呈降低趋势,变化幅度均表现为0~20 cm到20~40 cm土层下降平缓,20~40 cm到40~60 cm土层下降速度加快。
-
以3种林地类型有机碳(SOC)、有机碳密度(SOCD)及碳氮比(C/N)作为响应变量,以土壤理化指标、环境及地形因子作为解释变量,最终选取了21个影响因子。其中土壤理化指标包括土壤含水量、毛管持水量、土壤密度、总孔隙度、毛管孔隙度、pH、有机质、水解性氮、全氮、速效钾、全钾、速效磷及全磷;环境因子包括年降雨量、年均温度和年均相对湿度;地形因子包括经度、纬度、海拔、坡度和坡向。在进行冗余分析(RDA)之前,首先通过SPSS 26.0共线性检验删除了强烈共线性因子和因子贡献率小于等于0.1的因子,EP、PK、EB土壤碳特征的影响因子最终分别保留12、14和12个,如图3所示。进行前向选择后,3种林地类型所保留的因子分别解释了96.41%、92.82%和97.00%的数据变异,数据可信性较高。RDA结果示意图表示:因子箭头在排序轴上投影的长度反映该因子对响应变量的解释量,投影长度越长,解释率越高,反之则越低;箭头夹角可以反映出该因子与土壤有机碳特征的相关性,夹角小于90°为正相关,反之为负相关,夹角越小相关性越强[33]。
EP的轴1的解释率为77.62%,轴2的解释率为18.79%(图3a)。土壤有机质、全氮、坡向和土壤密度在EP土壤有机碳含量和密度的投影长度最长,说明对EP土壤有机碳含量和密度的解释力较强。其中,土壤密度与有机碳含量和密度箭头连线方向相反,因此与EP土壤有机碳含量和密度呈负相关关系,即随着土壤密度的增大,土壤有机碳含量和密度降低;而土壤有机质、全氮和坡向与土壤有机碳含量和密度的夹角小于90°,呈正相关关系。其中以有机质夹角最小,表明其与EP土壤有机质含量和密度相关性最大;影响因子在土壤C/N上投影最长的为土壤含水量、有机质、速效钾和坡向,表明其对EP土壤C/N解释性较高。其中,C/N与速效钾夹角大于90°呈负相关关系,与土壤含水量、有机质和坡向夹角小于90°呈正相关关系,以坡向夹角最小,表明坡向与EP土壤C/N变化相关性最大。
PK的轴1解释率为70.24%,轴2的解释率为22.58%(图3b)。在PK有机碳含量和密度上投影最长的影响因子依次为土壤毛管持水量、全氮、水解性氮、有效磷、全磷、土壤密度、速效钾,表明以上7个因子对PK土壤有机碳含量和密度解释能力较强。其中,土壤密度和速效钾与PK有机碳含量和密度夹角大于90°,因此,二者与有机碳含量和密度呈负相关关系。而土壤毛管持水量、全氮、水解性氮、有效磷和全磷与有机碳含量和密度夹角小于90°,呈正相关关系。以毛管持水量夹角最小,表明其与PK土壤有机碳含量和密度相关性最大。在土壤C/N上投影最长的解释因子依次为土壤有机质、pH、毛管持水量、土壤密度、坡向和年均温度,说明上述6个因子对PK土壤有机碳含量和密度解释性较强。土壤密度与C/N箭头连线方向夹角大于90°,呈负相关关系;而土壤有机质、pH、毛管持水量、坡向和年均温度与C/N夹角均小于90°,呈正相关关系,以年均温度夹角最小,说明其与土壤C/N相关性最大。
EB的轴1解释率为85.49%,轴2的解释率为11.51%(图3c)。土壤有机质、全氮、水解性氮、全磷、有效磷、总孔隙度、pH和土壤密度在EB土壤有机碳含量和密度方向上的投影最大,说明对EB土壤有机碳含量和密度解释能力较强。其中,土壤pH和土壤密度与EB土壤有机碳含量和密度夹角大于90°,说明二者与土壤有机碳含量和密度呈负相关关系。而土壤有机质、全氮、水解性氮、全磷、有效磷、总孔隙度与有机碳含量和密度夹角小于90°,表明与土壤有机碳含量和密度呈正相关关系,以土壤有机质与有机碳含量、有效磷与土壤有机碳密度夹角最小。土壤有机质、有效磷、pH和海拔在土壤C/N上投影最长,表明4个因子能很好解释土壤C/N变化情况。其中,土壤C/N与pH的夹角大于90°呈负相关关系,与土壤有机质、有效磷和海拔夹角小于90°,呈正相关关系,以海拔与C/N夹角最小,表明两者相关性最大。
总之,土壤有机质、全氮、全磷、水解性氮和土壤密度是3种林地类型土壤有机碳含量和密度的主要影响因子,而环境和地形因子对其产生的影响较小。土壤C/N的主要影响因素为土壤有机质、pH和坡向,环境因子对其影响较小。
Soil Carbon Changes and Its Influencing Factors in Major Forest Types in the Subtropical Area of Yunnan Province
- Received Date: 2023-02-03
- Accepted Date: 2023-10-29
- Available Online: 2024-02-20
Abstract: