环境科学  2025, Vol. 46 Issue (6): 3846-3855   PDF    
滇东南岩溶与非岩溶区土壤有机碳空间分异及其影响因素
涂纯1,2,3, 毛诗画4, 孔祥胜1,2,3, 罗为群1,2,3, 覃星铭5, 李发东6,7, 张妍4     
1. 中国地质科学院岩溶地质研究所,桂林 541004;
2. 广西平果喀斯特生态系统国家野外科学观测研究站,百色平果喀斯特生态系统广西野外科学观测研究站,平果 531406;
3. 自然资源部、广西岩溶动力学重点实验室,联合国教科文组织国际岩溶研究中心,桂林 541004;
4. 西北大学城市与环境学院,西安 710127;
5. 南宁师范大学北部湾环境演变与资源利用教育部重点实验室,南宁 530001;
6. 中国科学院地理科学与资源研究所,北京 100101;
7. 中国科学院大学资源与环境学院,北京 100190
摘要: 岩溶与非岩溶区土壤有机碳(SOC)差异影响县域尺度土壤碳汇的精确评估,是岩溶区土壤碳研究的热点与难点. 以滇东南典型岩溶区云南省文山市为例,采用1km×1km网格采样法,调查研究了文山市0~20 cm SOC含量、密度及储量在岩溶与非岩溶区的空间差异,探讨了不同土地利用方式、石漠化、海拔及土壤理化性质对表层SOC的影响. 结果表明,文山市0~20 cm土层SOC储量为1 261.90万t,呈西部高、中北部低的分布格局;其中非岩溶区ω(SOC)和SOC密度平均值分别为20.88 g·kg-1和4.93 kg·m-2,分别比岩溶区高27.94%和23.55%(P < 0.01). 不同土地利用方式下ω(SOC)差异显著:灌木林地(24.58 g·kg-1) > 草地(20.92 g·kg-1) > 有林地(20.52 g·kg-1) > 园地(19.89 g·kg-1) > 其他林地(18.87 g·kg-1) > 旱地(16.16 g·kg-1) > 水田(15.58 g·kg-1),其中非岩溶区灌木林地、草地和其他林地SOC高出岩溶区35.95%~66.52%,表明非岩溶区丰富的植物凋落物输入有助于SOC累积. 文山市石漠化主要分布在北部和中东部地区,与ω(SOC)低值区(8.0~15.0 g·kg-1)分布较为一致,轻度、中度和重度石漠化区ω(SOC)为12.98~15.55 g·kg-1,显著低于无石漠化区和潜在石漠化区(17.34~17.82 g·kg-1P < 0.01). 地理探测分析表明,土壤总氮(TN)含量对SOC空间分异的单因子解释力超过0.70,其次为海拔、土壤容重和C/P,三者对SOC的单因子解释力在0.31~0.52;而TN与C/N、TP与C/P交互作用对SOC的解释力在0.73~0.95;表明土壤氮、磷等养分状况对SOC有重要影响. 由于石漠化的影响,文山市岩溶区ω(TN)为1.54~2.10 g·kg-1,低于岩溶区11.47%~27.25%,导致上述影响因子及其交互作用对SOC的解释力低于非岩溶区. 研究可为岩溶区县域尺度的土壤碳汇评估提供科学参考.
关键词: 岩溶区      非岩溶区      土壤有机碳(SOC)      土地利用方式      石漠化      影响因素     
Spatial Partitioning and Influencing Factors of Soil Organic Carbon in Karst and Non-Karst Regions of Southeast Yunnan Province
TU Chun1,2,3 , MAO Shi-hua4 , KONG Xiang-sheng1,2,3 , LUO Wei-qun1,2,3 , QIN Xing-ming5 , LI Fa-dong6,7 , ZHANG Yan4     
1. Institute of Karst Geology, Chinese Academy of Geological Sciences, Guilin 541004, China;
2. Pingguo Guangxi, Karst Ecosystem, National Observation and Research Station/Pingguo Baise, Karst Ecosystem, Guangxi Observation and Research Station, Pingguo 531406, China;
3. Key Laboratory of Karst Dynamics Ministry of Natural Resources, Guangxi/International Research Center on Karst under the Auspices of United Nations Educational, Scientific and Cultural Organization, Guilin 541004, China;
4. College of Urban and Environmental Sciences, Northwest University, Xi'an 710127, China;
5. Environmental Evolution and Resources Utilization in Beibu Gulf under Ministry of Education, Nanning Normal University, Nanning 530001, China;
6. Institution of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China;
7. College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100190, China
Abstract: The difference in soil organic carbon (SOC) between karst and non-karst areas affects the accurate assessment of soil carbon sinks at the county scale, making it a hot and difficult topic in soil carbon research in karst areas. Taking Wenshan County, a typical karst area in southeast Yunnan Province, as an example, a 1 km×1 km grid sampling method was used to conduct the spatial differences in SOC content, density, and storage in the 0-20 cm soil layer between karst and non-karst areas. The impacts of different land use patterns, karst rocky desertification, elevation, and soil physicochemical properties on surface SOC were also explored. The results showed that the SOC storage in the 0-20 cm soil layer was 12.619 million tons, with a distribution pattern of high in the west and low in the central and northern regions in Wenshan County. The mean values of ω(SOC) and SOC density in non-karst areas were 20.88 g·kg-1 and 4.93 kg·m-2, respectively, which were 27.94% and 23.55% higher than those in karst areas (P < 0.01). Significant differences in ω(SOC) were observed under different land use types, with shrubland (24.58 g·kg-1) > grassland (20.92 g·kg-1) > forestland (20.52 g·kg-1) > orchard (19.89 g·kg-1) > other woodland (18.87 g·kg-1) > dryland (16.16 g·kg-1) > paddy field (15.58 g·kg-1). The SOC in shrubland, grassland, and other woodland in non-karst areas was 35.95%-66.52% higher than that in karst areas, indicating that abundant plant litter input in non-karst areas contributed to SOC accumulation. Karst rocky desertification was mainly distributed in the northern, central, and eastern regions of Wenshan County, which was consistent with the distribution of low SOC areas (8.0-15.0 g·kg-1). The ω(SOC) in mild, moderate, and severe karst rocky desertification areas ranged from 12.98 to 15.55 g·kg-1, significantly lower than that in non-rocky desertification and potential rocky desertification areas (17.34-17.82 g·kg-1, P < 0.01). Geographical detection analysis showed that the single-factor explanatory power of soil total nitrogen (TN) content for SOC spatial differentiation exceeded 0.70, followed by elevation, soil bulk density, and C/P, with single-factor explanatory powers ranging from 0.31 to 0.52. The interaction effects of TN and C/N, as well as TP and C/P, had explanatory powers of 0.73 to 0.95 for SOC, indicating that soil nutrients such as nitrogen and phosphorus had important impacts on SOC. Due to the influence of karst rocky desertification, the ω(TN) in the karst area of Wenshan County was 1.54-2.10 g·kg-1, which was 11.47%-27.25% lower than that in non-karst areas, resulting in lower explanatory powers of the aforementioned influencing factors and their interactions for SOC in karst areas compared to those in non-karst areas. This study provides a scientific reference for assessing soil carbon sinks at the county scale in karst regions.
Key words: karst area      non-karst area      soil organic carbon(SOC)      land use types      rocky desertification      influence factors     

土壤有机碳(soil organic carbon,SOC)是表征土壤质量的重要指标,同时也是陆地生态系统最大的碳库,全球大约有1 500 Pg碳(1 Pg=1015 g)以有机碳形式存在于土壤中[1],其变化对全球气候和碳循环具有重要影响. 我国西南地区裸露型岩溶面积约53.26万km2,占其国土面积的27.36%,是我国石漠化最为严重、生态环境极为脆弱的地区[2],深入认识西南岩溶地区SOC特征,对于评估该地区土壤碳汇潜力和土壤健康状况,提升我国陆地生态系统碳汇能力具有重要意义[3].

在岩溶地区,强烈的岩溶作用塑造了地表-地下二元水文结构和地表复杂多样的微地貌形态,导致其土壤环境具有基岩出露,土层稀薄和分布不连续等基本特征[245],加上气候、海拔、土地利用和人类活动等诸多因素影响,相比非岩溶区,岩溶地区SOC分布具有更为强烈的空间异质性,增加了区域碳汇评估的难度[3]. 目前,省域尺度[67]、流域尺度[389]和样地尺度[1011]的SOC研究已广泛报道,谷佳慧等[12]以图幅调查发现,云南广南表层SOC含量高于非岩溶区;肖筱怡等[13]和方芳等[14]在广西桂林的研究认为,岩溶区水稻田、柑橘园和玉米地的SOC含量均高于非岩溶区. Zhu等[15]对岩溶和非岩溶区人工造林地的对比分析发现,植树造林前期非岩溶区土壤碳含量高于岩溶区,但随着年限延长,岩溶区具有更为持久的碳固持能力. 土地利用和石漠化影响下的SOC空间分异研究较为普遍[16 ~ 19];不同土地利用方式上,岩溶和非岩溶地区均表现为林地、灌木地和草地SOC比耕地富集,主要是林草地有更多的凋落物和根系输入土壤,为SOC累积提供了充足的外源碳;另外,石漠化过程也会显著影响区域SOC含量及分布,陈飞[8]认为,石漠化等级越低,植被生长状态越好,SOC越高;反之,石漠化加剧将促使SOC丧失. 刘方等[20]认为,岩溶区植被遭受严重破坏时,基岩大面积裸露,土壤数量和有机碳等养分逐步流失,随着石漠化程度的加剧,SOC下降显著;而王霖娇等[16]和盛茂银等[21]研究发现,石漠化区SOC含量的演变并不随石漠化程度的加剧而一直退化,随着石漠化等级不断增加,裸露岩石对大气沉降养分和岩溶产物的聚集效应逐渐凸显,重度石漠化环境下SOC反而得到明显改善.

综上,目前的研究虽然解释了岩溶与非岩溶区、不同土地类型和石漠化环境下SOC的空间差异,但研究多集中在样地或流域尺度,受调查样本的限制,区域SOC评估的精确性和准确性还有待验证,尤其是以县域尺度的高密度调查研究尚未涉及. 因此,本研究以滇东南石漠化相对严重的文山市为例,按照1 km×1 km网格布点采样,探讨岩溶与非岩溶区SOC含量、密度和储量的差异,明确县域尺度SOC的空间分异及土地利用方式、石漠化和土壤理化性质对SOC的影响,以期为岩溶地区土壤碳库的精确评估和制定土壤碳汇管理措施提供科学依据.

1 材料与方法 1.1 研究区概况

研究区位于云南省东南部文山壮族苗族自治州文山市(图 1),地理位置:103°43′~104°27′E,23°06′~23°44′N. 全市总面积2 967.17 km2,本研究使用1∶20万区域地质图提取的裸露型碳酸盐数据,确定文山市岩溶区面积为2 126.31 km2,占国土面积的71.66%,属于滇东南地区典型岩溶县;非岩溶区面积840.86 km2,占国土面积的28.34%,主要分布在文山市西部地区. 文山市地势总体西高东低,最高海拔2 974 m,最低海拔627 m,相对高差2 347 m,岩溶地貌广泛分布,属滇东南典型的岩溶峰丛洼地地貌. 根据第二次全国土壤普查数据,文山市分布有棕壤、黄棕壤、黄壤、红壤、赤红壤、水稻土、石灰土和紫色土等8个土类. 区内年平均气温18.4 ℃,年均降雨量1 187.8 mm,处于中亚热带季风气候区.

图 1 研究区地理位置示意 Fig. 1 Geographical location of the study area

1.2 样品采集与分析

根据《多目标区域地球化学调查规范(1∶250 000)》(DZ/T 0258-2014)[22],按照1 km×1 km网格采集0~20 cm表层土壤样品,采样密度为1个点·km-2;为减少随机误差的影响,根据现场地形地貌条件,在主采样点附近100 m范围内等量采取2~4个分样混合后,四分法留取1.00~1.50 kg土壤样品,并记录主采样点的经纬度坐标、海拔以及周围土地利用类型和石漠化状况等信息;2016年7月至2017年2月共采集土壤样品2 921个. 样品采集后进行自然风干,碾碎过20目孔筛,装瓶在阴凉环境下保存,将4 km2大格内样品等量组合成一个分析样,共计716个. 样品测试由安徽省地质实验研究所完成,采用重铬酸钾-外加热法测定SOC含量,蒸馏水浸提-酸度计测定土壤pH值,凯氏氮容量法测定土壤总氮(total nitrogen,TN)含量,X射线荧光光谱法测定土壤总磷(total phosphorus,TP)含量.

1.3 石漠化遥感解译数据来源及处理

本研究使用的遥感影像来源于Google earth engine云平台,获取的Landsat8 OLI影像时间为2017年12月,影像空间分辨率30 m,云量小于1%;影像通过辐射定标、大气校正、几何校正和图像增强等预处理. 石漠化信息提取采用人机交互解译法,根据遥感影像的光谱、纹理、色调及形状等特征,在岩溶地区建立石漠化遥感解译判读标志;并通过波段运算提取文山市归一化植被指数和归一化岩石指数,进一步计算植被覆盖率和基岩裸露率[20],综合后形成石漠化提取信息(表 1);之后结合前期野外调查数据进行石漠化解译数据修正,结合地形、坡度、坡向和土地利用类型分析,重新选取不同等级石漠化影像特征标志进行解译,综合解译准确率达90.91%.

表 1 文山市石漠化等级划分及遥感影像判读标志 Table 1 Karst rocky desertification grade and interpretation signs of remote sensing image in Wenshan County

目前石漠化等级存在众多划分方法,主要依据基岩裸露、植被和土被覆盖条件、坡度、土壤厚度等因素进行石漠化等级划分[23 ~ 25]. 本研究结合实际情况,选择植被覆盖率、基岩裸露率作为研究区石漠化分级指标,将文山市石漠化划分为无石漠化、潜在石漠化、轻度石漠化、中度石漠化和重度石漠化这5个等级,其中石漠化面积为轻度、中度和重度石漠化面积之和[26].

1.4 数据处理与分析

SOC密度计算公式如下[22]

(1)

式中,SOCD为SOC密度(kg·m-2),D为土壤深度(cm),ρ为土壤容重(g·cm-3),C为SOC含量(%),T为土壤中直径大于2 mm的砾石占比(%),本调查中土壤的砾石粒径均小于2 mm,因此T取值为0[27]. 因调查缺失土壤容重数据,参考已有研究[28],土壤容重使用以下公式计算:

(2)
(3)

SOC储量计算公式如下:

(4)

式中,SCR为SOC储量(t),SOCDi为采样点i有机碳密度(kg·m-2),Ai为采样点所在网格面积(m2).

采用地理探测器中的单因子探测和交互作用探测,分析土地利用、石漠化、海拔、土壤容重、pH、TN、TP、土壤碳氮比(C/N)、土壤碳磷比(C/P)和土壤氮磷比(N/P)等影响因素对SOC空间分异的影响. 在因子探测时通过q值大小来度量单个影响因素对研究区SOC空间分异的驱动力度,其取值范围为0~1,q值越大表明影响因素对SOC空间分异的影响越大. 交互作用探测用于识别两个影响因素对SOC的解释力强弱,可分为非线性减弱、单因子非线性较弱、双因子增强、独立和非线性增强这5类[29].

1.5 数据分析与制图

采用Excel 2010、SPSS 13.0、ArcGIS 10.6和Origin 2022软件对数据进行处理和作图. SOC含量、密度和储量数据在岩溶区与非岩溶区比较时,使用T检验法检验其差异显著性(P < 0.05);不同土地利用方式和不同等级石漠化区SOC比较时,使用多重比较LSD法检验其差异显著性(P < 0.05). 使用Pearson分析SOC和海拔及土壤理化性质的相关性,运用ArcGIS 10.6的克里金插值法绘制SOC含量、密度及储量的空间分布图,使用Origin 2022绘制影响因子交互作用图.

2 结果与分析 2.1 文山市表层SOC总体分布特征

统计显示,文山市非岩溶区0~20 cm土层中的ω(SOC)和SOC密度均大于岩溶区,其平均值分别为20.88 g·kg-1和4.93 kg·m-2,高出岩溶区27.94%和23.55%(表 2). 空间分布显示(图 2),SOC含量、SOC密度和储量均呈现西部高、中北部低的特征;在薄竹镇南部、坝心乡西部和小街镇北部等文山国家自然保护区,ω(SOC)、SOC密度和储量分别在40.0~65.0 g·kg-1、9.5~13.5 kg·m-2和35.0万~55.0万t,显著高于红甸乡、卧龙街道和马塘镇等农业活动相对集中的地区. 文山市表层SOC储量为1 261.90万t,其中岩溶区因其面积占全市国土空间的71.66%,其储量为858.39万t,高出非岩溶区1.13倍;但从SOC密度平均值来看,岩溶区仅为4.10 kg·m-2,比非岩溶区低7.35%;其中在小街镇等森林覆盖面积较大的地区,SOC密度高达6.38 kg·m-2,虽然该区非岩溶区面积小于岩溶区,但SOC储量高于非岩溶区,表明文山市非岩溶区表层土壤具有更大的固碳能力.

表 2 文山市不同乡镇0~20 cm土层SOC含量、SOC密度和SOC储量统计 Table 2 Statistics of SOC content, SOC density and SOC storage in 0-20 cm soil layer in different towns of Wenshan County

图 2 文山市0~20 cm土层SOC含量、SOC密度和SOC储量空间分布 Fig. 2 Spatial distribution of SOC content, SOC density, and SOC storage in 0-20 cm soil layer in Wenshan County

2.2 不同土地利用方式下SOC差异

不同土地利用方式下ω(SOC)存在明显差异,总体表现为:灌木林地(24.58 g·kg-1) > 草地(20.92 g·kg-1) > 有林地(20.52 g·kg-1) > 园地(19.89 g·kg-1) > 其他林地(18.87 g·kg-1) > 旱地(16.16 g·kg-1) > 水田(15.58 g·kg-1). 其中灌木林地、草地和有林地SOC含量显著大于旱地和水田(26.98%~57.78%,P < 0.01,表 3),表明受人类活动干扰较少的土地利用类型区具有更高的SOC含量. 非岩溶区各土地类型中的SOC含量均高于岩溶区,其中在灌木林地、其他林地和草地,非岩溶区SOC含量高出岩溶区35.95%,而在水田、旱地和园地,非岩溶区SOC含量高出岩溶区10.50%~16.26%,说明非岩溶区相对岩溶区SOC的高累积现象更多体现在凋落物输入较大的林地和草地等土地类型中.

表 3 文山市不同土地利用方式下0~20 cm土层SOC含量统计1) Table 3 Statistics of SOC content in 0-20 cm soil layer under different land use types in Wenshan County

2.3 不同等级石漠化区SOC差异

文山市2017年石漠化面积为837.00 km2,其中轻度、中度和重度石漠化面积分别为639.64、176.30和21.06 km2,轻度石漠化面积占石漠化总面积的76.42%,主要分布在文山市北部和中东部地区(图 3). 对比分析SOC和石漠化的空间分布(图 2图 3),发现文山市ω(SOC)在8.00~15.00 g·kg-1的地区与石漠化分布区域较好重叠. 进一步分析不同等级石漠化区的SOC含量差异(表 4),显示轻度、中度和重度石漠化区ω(SOC)平均值分别为15.55、15.19和12.98 g·kg-1,显著低于无石漠化区(17.82 g·kg-1P < 0.01)和潜在石漠化区(17.34 g·kg-1P < 0.01),其中重度石漠化区SOC含量相比无石漠化区和潜在石漠化区分别显著降低37.29%和33.59%(P < 0.01),表明文山市石漠化程度加重将加速表层SOC损失.

图 3 文山市石漠化空间分布 Fig. 3 Spatial distribution of rocky desertification in Wenshan County

表 4 文山市石漠化区不同土地利用方式下0~20 cm土层SOC含量统计1) Table 4 Statistics of SOC content in 0-20 cm soil layer under different land use types in rocky desertification area in Wenshan County

随着石漠化严重程度减轻,SOC含量也逐步增加,这种增加量在灌木林地表现最为明显. 研究显示在重度石漠化区灌木林地的ω(SOC)仅为14.84 g·kg-1表 4),随着石漠化问题逐步向中度、轻度、潜在和无石漠化改善时,SOC含量分别增加了6.13%、24.39%、52.90%和69.54%,表明文山市灌木林地的恢复对于改善石漠化问题和提升表层SOC水平具有重要作用.

2.4 SOC空间分异影响因素分析

单因子探测结果显示,TN、海拔、土壤容重、C/P、石漠化和N/P是解释SOC的主要影响因子,其中TN对SOC空间分异的解释力(Q值)超过0.70,其次为海拔、土壤容重和C/P,其对有机碳的空间分异解释力超过0.30(表 5). 交互作用探测结果表明,任意两个影响因子协同对SOC空间分异的解释能力明显高于单因子的作用,呈现双因子增强或非线性增强的特征(图 4);其中TN与各因子交互后对SOC的解释力均在0.70以上;其次TP和C/P、N/P的交互作用对其解释力也较强,表明土壤氮、磷等养分状况对SOC有重要影响. 相比土地利用和石漠化等人为影响因子,海拔作为外部环境影响因子对SOC的空间分异具有更强的解释能力,表明地理环境在驱动SOC形成和分布的重要性. Pearson相关分析显示,SOC与TN、TP、C/N、C/P、N/P及海拔呈显著正相关,而与pH、土壤容重呈显著负相关(表 6). 无论单因子探测、交互因子探测或者相关性分析,岩溶区各影响因子对SOC的解释能力均低于非岩溶区,表明多因素协同研究对于分析复杂岩溶环境下的SOC空间变异十分必要.

表 5 文山市0~20 cm土层SOC空间分异的单因子探测分析1) Table 5 Single factor analysis of SOC spatial partitioning in 0-20 cm soil layer in Wenshan County

X1~X10依次表示土地利用、石漠化、海拔、土壤容重、pH、TN、TP、C/N、C/P和N/P;*表示双因子增强,**表示非线性增强 图 4 文山市0~20 cm土层SOC影响因子交互作用探测结果 Fig. 4 Interaction detection results of SOC influencing factors in 0-20 cm soil layer in Wenshan County

表 6 文山市0~20 cm土层SOC与海拔及土壤理化性质的相关性分析1) Table 6 Correlation analysis of SOC in 0-20 cm soil layer with elevation and soil physicochemical properties in Wenshan County

3 讨论 3.1 文山市表层SOC空间分异特征

通过对文山市进行网格化的高密度土样采样,研究发现文山市岩溶区ω(SOC)平均值为16.32 g·kg-1,低于非岩溶区(20.88 g·kg-1),与之前在西南地区的研究相反;比如钟聪等[7]和方芳等[14]研究认为广西岩溶区表层ω(SOC)在17.23~42.52 g·kg-1,高于非岩溶区(13.88~25.99 g·kg-1). Zhu等[15]通过野外调查和资料收集,发现西南岩溶区ω(SOC)平均值在29.15 g·kg-1,显著高于非岩溶区(18.40 g·kg-1). 另外,黄先飞等[30]在贵州(15.44~47.41 g·kg-1),邵晗等[31]在云南(16.22~44.23 g·kg-1)的研究均认为岩溶区具有较高的SOC含量.

本研究结果与之前研究不一致,可能是文山市处于云南石漠化较为严重的地区[32 ~ 34],超过30%的采样点位于轻度以上石漠化区. 有研究认为[2035],石漠化地区地表植被遭到破坏,外源有机碳输入减少;同时,植物的减少致使土壤表层丧失保护作用,导致其受侵蚀风险加强,促进了表层土壤有机质和养分流失. 而在非岩溶区,土壤呈连续分布,植被生长相对较好,地表动植物残体相对聚集,可为土壤提供相对充足的碳源[35]. 从空间分布来看(图 2),文山市SOC低值区主要分布在红甸乡、德厚镇、马塘镇和卧龙街道等文山市中北部地区,该区域也是石漠化严重区域;而SOC高值区位于文山市西部等非岩溶区,土层深厚且乔灌植被生长茂盛,土壤中外源碳输入高于北部和中东部岩溶区,显示了文山市非岩溶区具有良好的土壤固碳能力. 研究证实了县域尺度高密度采样对于精确评估岩溶与岩溶区SOC库水平的重要性.

3.2 土地利用方式和石漠化条件下SOC的空间差异

不同土地利用方式下的植被类型差异决定了归还土壤的凋落物和根系分泌物的数量和质量,进而影响SOC含量[3637]. 本研究表明,无论在岩溶区还是非岩溶区,灌木林地、有林地和草地表层SOC含量均大于水田和旱地. 这与黄先飞等[30]、白义鑫等[36]和邵晗等[31]研究的结果一致. 林灌草地受人为干扰少,生态系统稳定,尤其是林地和灌木,每年有大量的凋落物富集在表土层,利于SOC累积[36];而且林地根系生物量较大,根系活动及分泌物有助于土壤团聚体形成,增强了SOC的物理保护作用[11]. 相反,水田和旱地中的农作物会定期收割,地表凋落物数量有限,同时传统的耕作方式破坏了土壤的团聚体结构,影响SOC的物理稳定性,进而促进土壤微生物对其分解利用,致使耕地SOC相对减少[38].

本研究显示无石漠化和潜在石漠化区SOC含量高于轻度、中度和重度石漠化区,体现了石漠化发展过程与土壤退化的一致性关系,这与王纳伟等[35]、魏兴琥等[39]和王兴富等[40]研究的结果一致. 石漠化地区由于植被破坏、土壤侵蚀严重、土层稀薄及养分含量低等不利因素的共同作用,很大程度上限制了SOC的积累和稳定. 相反,石漠化程度较轻的地区,植被生存环境相对较好,物种多样性更加丰富,SOC含量会高于石漠化严重地区[41]. 研究显示随着石漠化程度降低,SOC含量也逐步增加,尤其是灌木林地SOC增加最为显著,表明灌木林地恢复对于改善石漠化区SOC水平有重要作用.

3.3 文山市SOC空间分异的影响因素

本研究显示在单因子探测中,TN对SOC的解释力最为突出,且两者显著正相关(表 5表 6),表明土壤中碳与氮元素存在强烈的耦合关系. 外源氮输入能促使植物叶片增大,增加叶片有机碳并降低凋落物碳氮比,进而促进凋落物分解和土壤腐殖质形成[4243]. 另外,土壤中氮素富集可以抑制微生物对SOC矿化和土壤呼吸,利于SOC含量提升[44]. 海拔和土壤容重对SOC也具有较强的单因子解释力,其中海拔与SOC呈显著负相关性,土壤容重与SOC呈显著负相关. 海拔对土壤所处区域的水热条件具有调控作用,海拔越高,土壤温度越低、湿度越大,有机质的矿化速率和微生物活性变慢,利于有机质累积[2745]. 土壤容重影响土壤透气、持水和入渗能力,容重较低的土壤疏松多孔,分布较多大团聚体,利于土壤微生物活动和植物根系发育,因此SOC含量较高[4647].

从交互作用探测结果来看,TN与其他影响因子协同对SOC产生双因子增强或非线性增强影响,尤其是TN与C/N的交互作用对SOC解释力超过0.90. C/N可在一定程度上影响微生物活性,常用来指示SOC的分解或累积,较低C/N有利于微生物对有机质分解,反之微生物在分解有机质过程中存在氮限制,从而有利于SOC保存[4849]. 另外,TP与C/P交互作用对SOC的解释力也较强;在亚热带和热带森林的研究认为[5051],土壤氮素相对富集,而磷是限制微生物生物量的主要元素,磷添加可减少微生物从土壤有机质中获取磷,并能改善土壤微生物群落结构和真菌/细菌比,进而减少微生物对SOC的分解. 较低的C/P有助于微生物在分解有机质过程中释放有效磷进入土壤,促进植物生长;而C/P较高,微生物会出现磷限制,进而与植物竞争土壤磷,抑制植物生长.

本研究显示,岩溶区各影响因子对SOC的解释能力均低于非岩溶区,可能原因是文山市作为石漠化较为严重的地区,土壤环境的空间异质性较大,土壤养分流失也较快. 研究发现石漠化作为影响因子对SOC有显著影响,其与TN交互作用对SOC的解释力达到0.75,而岩溶区ω(TN)为1.54~2.10 g·kg-1,分别低于非岩溶区11.47%~27.25%,同时也低于之前在西南岩溶区的研究(4.71~13.80 g·kg-1[1552],这表明文山市石漠化区土壤可能存在氮限制,这将进一步影响植物生长和微生物活动,进而限制SOC形成. 另外,海拔作为外部环境因子,其在岩溶区对SOC的解释力低于非岩溶区,其原因是文山市非岩溶区主要分布在西部,采样点平均海拔1 642.47 m,而岩溶区分布在北部和中东部,采样点平均海拔1 536.78 m. 在高海拔区,温度较低,枯枝落叶和腐殖质不易被微生物分解利用,且受人类活动的干预较少,土壤中的有机质能最大限度保留下来;而在低海拔区,水热条件不仅有利于微生物对有机碳的分解利用,而且植物的群落结构和物种丰富度也不如非岩溶区,加上低海拔区人为扰动频繁,进一步加速了岩溶区SOC的分解. 因此,文山市应在高海拔的西部非岩溶区,维持林地和灌木地的保护,稳定现有SOC储量;而在北部和中东部低海拔岩溶区,应积极开展石漠化生态修复,实施坡耕地退耕还林、封山育林和植树造林等工程措施,优化土地利用结构;同时在水田、旱地实施秸秆还田、免耕或少耕等管理措施,提升SOC库水平.

4 结论

文山市0~20 cm表层SOC、密度和储量均呈现西部高、中北部低的分布格局,其中非岩溶区SOC含量显著高于岩溶区;不同土地利用方式下SOC含量差异显著,在岩溶和非岩溶区均表现为灌木林地 > 草地 > 有林地 > 园地 > 其他林地 > 旱地 > 水田;文山市石漠化主要分布在北部和中东部地区,与SOC低值区分布具有一致性,随着石漠化程度加剧,表层SOC逐步下降,体现了石漠化发展过程与土壤退化的一致性关系. 地理探测器分析表明,TN是SOC的空间分异主要影响因子,其次为海拔、土壤容重和C/P;同时,TN与C/N、TP与C/P交互作用对SOC的解释能力较强;然而由于石漠化的影响,文山市岩溶区土壤可能存在氮限制,这导致上述影响因子对SOC的解释能力低于非岩溶区.

参考文献
[1] Lal R. Soil carbon sequestration impacts on global climate change and food security[J]. Science, 2004, 304(5677): 1623-1627. DOI:10.1126/science.1097396
[2] 曹建华, 袁道先, 章程, 等. 受地质条件制约的中国西南岩溶生态系统[J]. 地球与环境, 2004, 32(1): 1-8.
Cao J H, Yuan D X, Zhang C, et al. Karst ecosystem constrained by geological conditions in Southwest China[J]. Earth and Environment, 2004, 32(1): 1-8. DOI:10.3969/j.issn.1672-9250.2004.01.001
[3] 张珍明, 周运超, 田潇, 等. 喀斯特小流域土壤有机碳空间异质性及储量估算方法[J]. 生态学报, 2017, 37(22): 7647-7659.
Zhang Z M, Zhou Y C, Tian X, et al. Study on spatial heterogeneity and reserve estimation of soil organic carbon in a small Karst catchment[J]. Acta Ecologica Sinica, 2017, 37(22): 7647-7659.
[4] 张殿发, 王世杰, 周德全, 等. 贵州省喀斯特地区土地石漠化的内动力作用机制[J]. 水土保持通报, 2001, 21(4): 1-5.
Zhang D F, Wang S J, Zhou D Q, et al. Intrinsic driving mechanism of land rocky desertification in Karst regions of Guizhou Province[J]. Bulletin of Soil and Water Conservation, 2001, 21(4): 1-5. DOI:10.3969/j.issn.1000-288X.2001.04.001
[5] 张信宝, 王世杰, 曹建华, 等. 西南喀斯特山地水土流失特点及有关石漠化的几个科学问题[J]. 中国岩溶, 2010, 29(3): 274-279.
Zhang X B, Wang S J, Cao J H, et al. Characteristics of water loss and soil erosion and some scientific problems on Karst rocky desertification in Southwest China karst area[J]. Carsologica Sinica, 2010, 29(3): 274-279. DOI:10.3969/j.issn.1001-4810.2010.03.009
[6] 薛梦雨. 云南省表层土壤有机碳时空分异特征及影响因素分析[D]. 昆明: 云南大学, 2017.
Xue M Y. Spatiotemporal variations and environmental controls of topsoil organic carbon in Yunnan Province[D]. Kunming: Yunnan University, 2017.
[7] 钟聪, 李小洁, 何园燕, 等. 广西土壤有机质空间变异特征及其影响因素研究[J]. 地理科学, 2020, 40(3): 478-485.
Zhong C, Li X J, He Y Y, et al. Spatial variation of soil organic matter and its influencing factors in Guangxi, China[J]. Scientia Geographica Sinica, 2020, 40(3): 478-485.
[8] 陈飞. 喀斯特区域土壤有机碳的空间分布特征及其影响因素[D]. 贵阳: 贵州大学, 2022.
Chen F. Spatial distribution characteristics and influencing factors of soil organic carbon in Karst Area[D]. Guiyang: Guizhou University, 2022.
[9] 申楷慧, 魏识广, 李林, 等. 漓江流域喀斯特森林土壤有机碳空间分布格局及其驱动因子[J]. 环境科学, 2024, 45(1): 323-334.
Shen K H, Wei S G, Li L, et al. Spatial distribution patterns of soil organic carbon in Karst forests of the Lijiang River Basin and its driving factors[J]. Environmental Science, 2024, 45(1): 323-334.
[10] 蓝家程, 王俊贤, 王莎莎, 等. 喀斯特石漠化治理措施对土壤颗粒有机碳与团聚体有机碳的影响[J]. 中国岩溶, 2022, 41(5): 773-783.
Lan J C, Wang J X, Wang S S, et al. Impact of controlling Karst rocky desertification on soil particulate organic carbon and aggregate-associated organic carbon[J]. Carsologica Sinica, 2022, 41(5): 773-783.
[11] 周业恒, 肖孔操, 夏世斌, 等. 广西喀斯特峰丛洼地造林对土壤有机碳、氮及其组分的影响[J]. 农业现代化研究, 2024, 45(2): 316-325.
Zhou Y H, Xiao K C, Xia S B, et al. Effects of afforestation on soil carbon and nitrogen fractions in the Karst-cluster depression region of Guangxi[J]. Research of Agricultural Modernization, 2024, 45(2): 316-325.
[12] 谷佳慧, 杨奇勇, 蒋忠诚, 等. 广南县幅岩溶区与非岩溶区土壤碳氮磷生态化学计量比空间变异分析[J]. 中国岩溶, 2018, 37(5): 761-769.
Gu J H, Yang Q Y, Jiang Z C, et al. Spatial variation analysis of soil carbon, nitrogen and phosphorus eco-stoichiometric ratios in Karst and Non-Karst areas of Guangnan county, Yunnan, China[J]. Carsologica Sinica, 2018, 37(5): 761-769.
[13] 肖筱怡, 靳振江, 冷蒙, 等. 岩溶与非岩溶区水稻土团聚体细菌群落结构和功能类群比较[J]. 环境科学, 2022, 43(7): 3865-3875.
Xiao X Y, Jin Z J, Leng M, et al. Comparison of bacterial community structure and functional groups of paddy soil aggregates between Karst and Non-Karst areas[J]. Environmental Science, 2022, 43(7): 3865-3875.
[14] 方芳, 靳振江, 李强, 等. 岩溶区与非岩溶区土壤有机碳、养分及特征元素对比[J]. 桂林理工大学学报, 2016, 36(3): 550-556.
Fang F, Jin Z J, Li Q, et al. Comparison of soil organic carbon nutrients and characteristic elements between Karst and Non-Karst area[J]. Journal of Guilin University of Technology, 2016, 36(3): 550-556. DOI:10.3969/j.issn.1674-9057.2016.03.021
[15] Zhu X C, Ma M G, Tateno R, et al. Effects of vegetation restoration on soil carbon dynamics in Karst and Non-Karst regions in Southwest China: a synthesis of multi-source data[J]. Plant and Soil, 2022, 475(1): 45-59.
[16] 王霖娇, 盛茂银, 李瑞. 中国南方喀斯特石漠化演替过程中土壤有机碳的响应及其影响因素分析[J]. 生态科学, 2016, 35(1): 47-55.
Wang L J, Sheng M Y, Li R. Response of soil organic carbon to rocky desertification succession in South China Karst and the analysis of its influence factors[J]. Ecological Science, 2016, 35(1): 47-55.
[17] 苏同庆, 崔婷婷, 张建兵, 等. 土地利用方式对广西平果喀斯特土壤碳氮磷全量与易利用组分的影响[J]. 中国岩溶, 2023, 42(2): 311-320.
Su T Q, Cui T T, Zhang J B, et al. Effect of land utilization patterns on total and easy-to-use components of soil carbon, nitrogen and phosphorus in the Karst area of Pingguo, Guangxi[J]. Carsologica Sinica, 2023, 42(2): 311-320.
[18] Chen H S, Zhang W, Wang K L, et al. Soil organic carbon and total nitrogen as affected by land use types in Karst and Non-Karst areas of Northwest Guangxi, China[J]. Journal of the Science of Food and Agriculture, 2012, 92(5): 1086-1093. DOI:10.1002/jsfa.4591
[19] 龙启霞, 蓝家程, 姜勇祥. 生态恢复对石漠化地区土壤有机碳累积特征及其机制的影响[J]. 生态学报, 2022, 42(18): 7390-7402.
Long Q X, Lan J C, Jiang Y X. Effects of soil organic carbon fractions on aggregates under ecological restoration in rocky desertification region[J]. Acta Ecologica Sinica, 2022, 42(18): 7390-7402.
[20] 刘方, 王世杰, 刘元生, 等. 喀斯特石漠化过程土壤质量变化及生态环境影响评价[J]. 生态学报, 2005, 25(3): 639-644.
Liu F, Wang S J, Liu Y S, et al. Changes of soil quality in the process of Karst rocky desertification and evaluation of impact on ecological environment[J]. Acta Ecologica Sinica, 2005, 25(3): 639-644. DOI:10.3321/j.issn:1000-0933.2005.03.035
[21] 盛茂银, 刘洋, 熊康宁. 中国南方喀斯特石漠化演替过程中土壤理化性质的响应[J]. 生态学报, 2013, 33(19): 6303-6313.
Sheng M Y, Liu Y, Xiong K N. Response of soil physical-chemical properties to rocky desertification succession in South China Karst[J]. Acta Ecologica Sinica, 2013, 33(19): 6303-6313.
[22] DZ/T 0258-2014, 多目标区域地球化学调查规范(1∶250000)[S].
DZ/T 0258-2014, Specification of multi-purpose regional geochemical survey (1∶250000)[S].
[23] 李瑞玲, 王世杰, 周德全, 等. 贵州岩溶地区岩性与土地石漠化的相关分析[J]. 地理学报, 2003, 58(2): 314-320.
Li R L, Wang S J, Zhou D Q, et al. The correlation between rock desertification and lithology in Karst area of Guizhou[J]. Acta Geographica Sinica, 2003, 58(2): 314-320. DOI:10.3321/j.issn:0375-5444.2003.02.019
[24] 张信宝, 王世杰, 白晓永, 等. 贵州石漠化空间分布与喀斯特地貌、岩性、降水和人口密度的关系[J]. 地球与环境, 2013, 41(1): 1-6.
Zhang X B, Wang S J, Bai X Y, et al. Relationships between the spatial distribution of Karst land desertification and geomorphology, lithology, precipitation, and population density in Guizhou Province[J]. Earth and Environment, 2013, 41(1): 1-6.
[25] 王正雄, 蒋勇军, 张远嘱, 等. 基于GIS与地理探测器的岩溶槽谷石漠化空间分布及驱动因素分析[J]. 地理学报, 2019, 74(5): 1025-1039.
Wang Z X, Jiang Y J, Zhang Y Z, et al. Spatial distribution and driving factors of Karst rocky desertification based on GIS and geodetectors[J]. Acta Geographica Sinica, 2019, 74(5): 1025-1039.
[26] 蒋忠诚, 罗为群, 童立强, 等. 21世纪西南岩溶石漠化演变特点及影响因素[J]. 中国岩溶, 2016, 35(5): 461-468.
Jiang Z C, Luo W Q, Tong L Q, et al. Evolution features of rocky desertification and influence factors in Karst areas of Southwest China in the 21st century[J]. Carsologica Sinica, 2016, 35(5): 461-468.
[27] 楚梦玮, 侯巍楹, 净婷菲, 等. 云南省深层土壤有机碳的估算及其影响因素分析[J]. 云南农业大学学报(自然科学), 2021, 36(3): 516-523.
Chu M W, Hou W Y, Jing T F, et al. Estimation of deep soil organic carbon in Yunnan province and its influencing factors[J]. Journal of Yunnan Agricultural University (Natural Science), 2021, 36(3): 516-523.
[28] Wu H B, Guo Z T, Peng C H. Distribution and storage of soil organic carbon in China[J]. Global Biogeochemical Cycles, 2003, 17(2). DOI:10.1029/2001GB001844
[29] 王劲峰, 徐成东. 地理探测器: 原理与展望[J]. 地理学报, 2017, 72(1): 116-134.
Wang J F, Xu C D. Geodetector: principle and prospective[J]. Acta Geographica Sinica, 2017, 72(1): 116-134.
[30] 黄先飞, 周运超, 张珍明. 喀斯特石漠化区不同土地利用方式下土壤有机碳分布特征[J]. 水土保持学报, 2017, 31(5): 215-221.
Huang X F, Zhou Y C, Zhang Z M. Distribution characteristics of soil organic carbon under different land uses in a Karst rocky desertification area[J]. Journal of Soil and Water Conservation, 2017, 31(5): 215-221.
[31] 邵晗, 王虎, 王妍, 等. 岩溶石漠化地区不同利用方式对土壤肥力和重金属质量分数的影响[J]. 浙江农林大学学报, 2022, 39(3): 635-643.
Shao H, Wang H, Wang Y, et al. Effects of different land use modes on soil fertility and heavy metal contents in Karst rocky desertification area[J]. Journal of Zhejiang A & F University, 2022, 39(3): 635-643.
[32] 陈棋, 张超, 田湘云, 等. 云南省2000—2020年石漠化时空演变分析[J]. 浙江农林大学学报, 2023, 40(2): 417-426.
Chen Q, Zhang C, Tian X Y, et al. Spatio-temporal evolution process of rocky desertification in Yunnan Province from 2000 to 2020[J]. Journal of Zhejiang A & F University, 2023, 40(2): 417-426.
[33] 罗旭玲, 王世杰, 白晓永, 等. 西南喀斯特地区石漠化时空演变过程分析[J]. 生态学报, 2021, 41(2): 680-693.
Luo X L, Wang S J, Bai X Y, et al. Analysis on the spatio-temporal evolution process of rocky desertification in Southwest Karst area[J]. Acta Ecologica Sinica, 2021, 41(2): 680-693.
[34] 郭丽琴, 赵志芳, 代启学, 等. 基于RS和GIS文山州石漠化时空演变特征及成因研究[J]. 国土资源遥感, 2017, 29(S1): 106-113.
Guo L Q, Zhao Z F, Dai Q X, et al. Temporal and spatial evolution and genesis of rocky desertification based on RS and GIS in Wenshan Prefecture[J]. Remote Sensing for Land & Resources, 2017, 29(S1): 106-113. DOI:10.6046/gtzyyg.2017.s1.18
[35] 王纳伟, 吴立潮, 钟杰, 等. 湘中石漠化区土壤有机碳和全氮的变化特征[J]. 南方农业学报, 2013, 44(4): 616-620.
Wang N W, Wu L C, Zhong J, et al. Study on soil organic carbon and total nitrogen variation in Rocky Desertification area of central Hunan[J]. Journal of Southern Agriculture, 2013, 44(4): 616-620. DOI:10.3969/j:issn.2095-1191.2013.4.616
[36] 白义鑫, 盛茂银, 胡琪娟, 等. 西南喀斯特石漠化环境下土地利用变化对土壤有机碳及其组分的影响[J]. 应用生态学报, 2020, 31(5): 1607-1616.
Bai Y X, Sheng M Y, Hu Q J, et al. Effects of land use change on soil organic carbon and its components in Karst Rocky desertification of Southwest China[J]. Chinese Journal of Applied Ecology, 2020, 31(5): 1607-1616.
[37] 徐广平, 李艳琼, 沈育伊, 等. 桂林会仙喀斯特湿地水位梯度下不同植物群落土壤有机碳及其组分特征[J]. 环境科学, 2019, 40(3): 1491-1503.
Xu G P, Li Y Q, Shen Y Y, et al. Soil organic carbon distribution and components in different plant communities along a water table gradient in the Huixian Karst Wetland in Guilin[J]. Environmental Science, 2019, 40(3): 1491-1503.
[38] 严毅萍, 曹建华, 杨慧, 等. 岩溶区不同土地利用方式对土壤有机碳碳库及周转时间的影响[J]. 水土保持学报, 2012, 26(2): 144-149.
Yan Y P, Cao J H, Yang H, et al. The impact of different soil types on soil organic carbon pool and turnover in Karst Area[J]. Journal of Soil and Water Conservation, 2012, 26(2): 144-149.
[39] 魏兴琥, 李森, 罗红波, 等. 粤北石漠化过程土壤与植被变化及其相关性研究[J]. 地理科学, 2008, 28(5): 662-666.
Wei X H, Li S, Luo H B, et al. Changes and correlation of soil and vegetation in process of rock desertification in Northern Guangdong Province[J]. Scientia Geographica Sinica, 2008, 28(5): 662-666. DOI:10.3969/j.issn.1000-0690.2008.05.012
[40] 王兴富, 黄先飞, 胡继伟, 等. 喀斯特石漠化过程中小生境及岩性的演替对土壤有机碳的影响[J]. 水土保持学报, 2020, 34(4): 295-303.
Wang X F, Huang X F, Hu J W, et al. Relationship among soil organic carbon and small environment and lithology in the rocky desertification process in different karst landforms[J]. Journal of Soil and Water Conservation, 2020, 34(4): 295-303.
[41] Qi D, Wieneke X, Zhou X, et al. Succession of plant community composition and leaf functional traits in responding to Karst Rocky desertification in the Wushan County in Chongqing, China[J]. Community Ecology, 2017, 18(2): 157-168. DOI:10.1556/168.2017.18.2.5
[42] 李成, 王让会, 李兆哲, 等. 中国典型农田土壤有机碳密度的空间分异及影响因素[J]. 环境科学, 2021, 42(5): 2432-2439.
Li C, Wang R H, Li Z Z, et al. Spatial differentiation of soil organic carbon density and influencing factors in typical croplands of China[J]. Environmental Science, 2021, 42(5): 2432-2439.
[43] 刘德燕, 宋长春, 黄靖宇. 沼泽湿地植物光合特性及固"碳"潜势对外源氮输入的响应[J]. 环境科学学报, 2008, 28(2): 305-312.
Liu D Y, Song C C, Huang J Y. Effects of exogenous nitrogen supply on photosynthetic characteristics and carbon sequestration potential of calamagrostis angustifolia in freshwater marshes[J]. Acta Scientiae Circumstantiae, 2008, 28(2): 305-312. DOI:10.3321/j.issn:0253-2468.2008.02.013
[44] Wang X J, Wang J P, Xu M G, et al. Carbon accumulation in arid croplands of Northwest China: pedogenic carbonate exceeding organic carbon[J]. Scientific Reports, 2015, 5(1). DOI:10.1038/srep11439
[45] 尚斌, 邹焱, 徐宜民, 等. 贵州中部山区植烟土壤有机质含量与海拔和成土母质之间的关系[J]. 土壤, 2014, 46(3): 446-451.
Shang B, Zou Y, Xu Y M, et al. Relationship between SOM contents of tobacco fields and elevation and parent materials in central region of Guizhou Province[J]. Soils, 2014, 46(3): 446-451.
[46] 张穗粒, 盛茂银, 王霖娇, 等. 西南喀斯特长期植被修复对土壤有机碳组分的影响[J]. 生态学报, 2023, 43(20): 8476-8492.
Zhang S L, Sheng M Y, Wang L J, et al. Effects of long term vegetation restorations on soil organic carbon fractions in the Karst Rocky desertification ecosystem, Southwest China[J]. Acta Ecologica Sinica, 2023, 43(20): 8476-8492.
[47] 郑维熙, 周忠发, 朱粲粲, 等. 岩溶区土壤理化特征与土壤有机碳含量关系研究[J]. 环境科学与技术, 2021, 44(7): 91-99.
Zheng W X, Zhou Z F, Zhu C C, et al. Study on the relationship between soil physical and chemical characteristics and soil organic carbon content in Karst area[J]. Environmental Science & Technology, 2021, 44(7): 91-99.
[48] Li Y F, Yang R, Hu P L, et al. Lower sensitivity of soil carbon and nitrogen to regional temperature change in Karst forests than in Non-Karst forests[J]. Forests, 2023, 14(2). DOI:10.3390/F14020355
[49] 曹小玉, 李际平, 闫文德. 不同龄组杉木林土壤有机碳与氮磷钾分布特征及耦合关系[J]. 土壤通报, 2014, 45(5): 1137-1143.
Cao X Y, Li J P, Yan W D. Variation of contents of soil organic carbon, N, P and K and their coupling relationships in different age groups of Chinese fir plantations[J]. Chinese Journal of Soil Science, 2014, 45(5): 1137-1143.
[50] Fang X M, Zhang X L, Chen F S, et al. Phosphorus addition alters the response of soil organic carbon decomposition to nitrogen deposition in a subtropical forest[J]. Soil Biology and Biochemistry, 2019, 133: 119-128. DOI:10.1016/j.soilbio.2019.03.005
[51] Six J, Frey S D, Thiet R K, et al. Bacterial and fungal contributions to carbon sequestration in agroecosystems[J]. Soil Science Society of America Journal, 2006, 70(2): 555-569. DOI:10.2136/sssaj2004.0347
[52] 武亚楠, 喻理飞, 张丽敏, 等. 喀斯特高原区植被恢复过程中土壤碳特征及其影响因素[J]. 生态环境学报, 2020, 29(10): 1935-1942.
Wu Y N, Yu L F, Zhang L M, et al. Characteristics and influencing factors of soil carbon pool during vegetation restoration in Karst Plateau[J]. Ecology and Environmental Sciences, 2020, 29(10): 1935-1942.