1、第 45 卷 第 1 期2023 年 2 月Vol.45，No.1Feb.，2023冰川冻土JOURNAL OF GLACIOLOGY AND GEOCRYOLOGY长江源区布曲流域土壤有机碳分布特征及其影响因素张超飞1，2，马素萍1，何晓波3，4，汪少勇3，5（1.中国科学院 西北生态环境资源研究院 甘肃省油气资源研究重点实验室，甘肃 兰州 730000；2.中国科学院大学 地球与行星科学学院，北京 100049；3.中国科学院 西北生态环境资源研究院 冰冻圈科学国家重点实验室，甘肃 兰州 730000；4.中国科学院 西北生态环境资源研究院 内陆河流域生态水文重点实验室，甘肃 兰州 730
2、000；5.中国科学院大学 资源与环境学院，北京 100049）摘要：基于长江源区布曲流域23个采样点的138个土壤样品，分析了土壤有机碳的分布特征，并探讨1 m以内土壤有机碳含量的影响因素。结果表明，长江源区布曲流域010 cm、1020 cm、2030 cm、3040 cm、4050 cm 和 50100 cm 的 6层土壤有机碳含量的平均值标准差分别为（10.234.84）g kg-1、（10.185.19）g kg-1、（9.345.20）g kg-1、（9.044.41）g kg-1、（8.014.74）g kg-1、（9.404.67）g kg-1，其随土壤深度增加而降低（R2=0
3、.511）。研究区土壤有机碳含量与海拔并非线性相关（P0.05），而是在4 7005 100 m的海拔范围内，随海拔的升高而增加，然后在5 0005 100 m海拔处达到最大值后降低。土壤有机碳含量与pH值呈现极显著的负相关（P0.01），与碳氮比呈现极显著的正相关（P0.01），与体积含水量呈现显著的正相关关系（P0.05），表明长江源区布曲流域1 m以内土壤有机碳含量的主要相关因素是土壤的pH值、碳氮比和体积含水量。研究结果可以为长江源区布曲流域土壤碳循环的研究提供基础数据和科学依据。关键词：长江源区；布曲流域；多年冻土；有机碳；影响因素中图分类号：P593 文献标志码：A 文章编号：10
4、00-0240（2023）01-0233-100 引言 多年冻土是冰冻圈的重要组成部分，其是指温度低于0 时至少连续保存两年的岩土层1。多年冻土对气候变化的响应极为敏感2。在全球气候变暖的大背景下，目前多年冻土正在发生显著、快速和广泛的退化，主要表现有多年冻土面积减少和活动层厚度增加3，导致多年冻土区生物地球化学循环发生改变4，进而产生了一系列生态、水文和环境问题5-6。多年冻土退化会导致多年冻土区土壤中的大量有机碳分解释放，其通常以温室气体（如CO2和 CH4）的形式进入大气圈7，或以溶解性有机碳（dissolved organic carbon，DOC）的形式进入水圈8，进而对气候变暖产生
5、强烈的正反馈9。青藏高原是全球中低纬度高海拔冰冻圈分布的典型区，被称为世界的“第三极”10，其多年冻土主要分布在4 0005 500 m的海拔范围内，多年冻土实际分布面积约为115.02104 km2 11。目前青藏高原土壤在研究方法、尺度和深度上具有较大差异，这造成了对青藏高原多年冻土区土壤有机碳储量估算结果存在较大差异12。Mu等13评估青藏高原多年冻土区 025 m 深度土壤有机碳储量为（13277）Pg，Zhao等14评估青藏高原多年冻土区2 m以DOI：10.7522/j.issn.1000-0240.2022.0007ZHANG Chaofei，MA Suping，HE Xiaob
6、o，et al.Distribution characteristics and influencing factors of soil organic carbon in Buqu catchment in the source region of the Yangtze River，China J.Journal of Glaciology and Geocryology，2023，45（1）：233-242.张超飞，马素萍，何晓波，等.长江源区布曲流域土壤有机碳分布特征及其影响因素 J.冰川冻土，2023，45（1）：233-242.收稿日期：2021-09-09；修订日期：2021-1
7、1-21基金项目：三江源国家公园联合研究专项（LHZX-2020-11-1；LHZX-2020-10-1）；中国水科院科研专项（HY110145B0012021）；冰冻圈科学国家重点实验室自主课题项目（SKLCS-ZZ-2022）资助作者简介：张超飞，硕士研究生，主要从事多年冻土碳循环研究.E-mail：zhang_通信作者：马素萍，副研究员，主要从事沉积地球化学研究.E-mail：45 卷冰川冻土内的土壤有机碳储量为17.07 Pg，Ding等15评估青藏高原多年冻土区 3 m 以内的土壤有机碳储量为34.238.9 Pg，Wang等16评估青藏高原多年冻土区土壤总有机碳储量为 35.786
8、9.02 Pg，Wang 等17评估青藏高原多年冻土区3 m以内的土壤有机碳储量为21.69 Pg。从这些结果可以看出，青藏高原多年 冻 土 区 土 壤 碳 储 量 巨 大，但 是 不 确 定 范 围较大18。长江源区位于青藏高原腹地，是对全球气候变暖响应最为敏感的区域之一19，布曲流域是长江南源当曲流域内一条典型的多年冻土流域，受全球气候变暖的影响十分明显，由于其高海拔的地理位置和恶劣的气候条件，样品和数据获取难度较大，关于布曲流域土壤有机碳的研究相对匮乏。本研究通过在布曲流域实地挖取多年冻土剖面，进行土壤样品的采集，进而分析土壤有机碳的分布特征，在一定程度上弥补了布曲流域土壤有机碳数据的不
9、足，为布曲流域多年冻土区土壤有机碳储量的估算和模拟提供基础数据。此外，探讨了布曲流域土壤有机碳的影响因素，以期为认识区域上多年冻土区碳循环及其对气候变暖的响应和反馈提供科学依据。1 研究区概况 研究区（9124 9220 E，3249 3338 N）位于长江源区的布曲流域（图1），该流域位于重要气 候 分 界 线 的 唐 古 拉 山 脉 中 段，流 域 面 积 为 4 504.72 km2，海拔处于4 7036 081 m之间，平均海拔为5 152 m。气候寒冷干燥，空气稀薄，属于典型的 高 原 亚 寒 带 半 干 旱 气 候，多 年 平 均 气 温 为 -6.0，多年平均降水为443 mm。
10、流域内多年冻土广泛发育，其分布海拔范围为4 7005 300 m，土壤类型以沼泽土与高山草甸土为主。植被类型以高寒草甸（alpine meadow，AM）为主，AM种类组成较为丰富，主要是以莎草科为主的蒿草类和苔草类，如高山蒿草（Kobresia pygmae）、藏蒿草（Kobresia tibetica）、矮蒿草（Kobresia humilis）、青藏苔草（Carex moorcroftii）等，植被覆盖度为30%98%，生长最盛时高度为510 cm，根系深且发达20。2 材料与方法 2.1样品采集2020年08月中旬，在长江源区布曲流域进行土壤样品的采集（图1），参考唐古拉站气象、水文站
11、点的分布，按照海拔的变化，共取23个采样点。采用五点混合取样法采集土壤样品，在每个样点取 010 cm、1020 cm、2030 cm、3040 cm、4050 cm和50100 cm的6层土壤样品，3次重复。首先随机选取10 m10 m的样区，将每个样区平均分为九格，然后在中间一格以及任意周边两格采集 3份土壤样品，然后将同层位的3份土壤样品混合均匀成一份总的土壤样品，编号后进行装袋。使用手持GPS记录经度（longitude，LON）、纬度（latitude，LAT）和高程（altitude，ALT）。土壤样品在阴凉处风干，并清除根系、石砾，然后在中国科学院沙漠与沙漠化重点实验室进行土壤理
12、化性质的测定。实验室测定土壤样品的质地、pH 值、有机碳（soil organic carbon，SOC）、无机碳（total inorganic carbon，TIC）、全氮（total nitrogen，TN）、全磷（total phosphorus，TP）、全钾（total kalium，TK）和阳离子交换量（cation exchange capacity，CEC），同时用 200 cm3环刀取原状土用于测量土壤的容重（bulk density，BD）和体积含水量（volumetric water content，VWC）。2.2样品测定土壤SOC的测定方法为重铬酸钾-硫酸外加热容量法
13、，土壤样品中的有机碳和过量的重铬酸钾-硫酸溶液发生氧化还原反应，多余的重铬酸钾用硫酸亚铁标准液进行滴定，通过消耗的重铬酸钾量来计算土壤样品中的有机碳含量。TIC 的测定方法为气量法，TN 的测定方法为凯氏定氮法，碳氮比（C/N）通过SOC与TN的比值计算获得。土壤TP的测定采用碱熔-钼锑抗比色法，TK的测定方法为碱熔-火焰光度法，CEC的测定方法为乙酸铵法，pH值的测定方法为溶液-电极法。土壤质地的测定方法为筛分-沉降法，BD和VWC分别采用环刀法和图1研究区及采样点位置示意Fig.1Location of the study area and sampling sites2341 期张超飞等
14、：长江源区布曲流域土壤有机碳分布特征及其影响因素烘干法测定21。2.3数据处理年均气温（mean annual air temperature，MAT）数 据 和 年 均 降 水 量（mean annual precipitation，MAP）数据分别来自于国家青藏高原科学数据中心（https：/ km分辨率逐月平均气温数据集（http：/ km分辨率逐月降水量数据集（http：/ ArcGIS 10.5 计算并提取每个采样点的MAT和MAP数据（表1）。使用Origin 2018软件制图。使用SPSS 25对土壤有机碳含量和密度以及其他理化性质进行描述性分析和单因素方差分析，对有机碳含量和气
15、候因子进行线性回归分析，对有机碳含量和其他土壤理化性质进行皮尔逊相关分析。使用公式（1）计算有机碳存量15。SOCS=i=1nTi BDi SOCCi(1-Ci)100（1）式中：SOCS为有机碳存量（kg m-2）；Ti为第i层土壤的厚度（cm）；BDi为第 i 层土壤的容重（g cm-3）；SOCCi为第i层土壤有机碳含量（g kg-1）；Ci为第i层土壤中粒径大于2 mm的砾石所占的比率（%）。3 实验结果 3.1土壤有机碳分布特征3.1.1土壤有机碳在剖面中的分布特征在长江源区布曲流域 010 cm、1020 cm、2030 cm、3040 cm、4050 cm 和 50100 cm
16、的 6 个土层中，土壤有机碳含量的平均值标准差分别为（10.234.84）g kg-1、（10.185.19）g kg-1、（9.345.20）g kg-1、（9.044.41）g kg-1、（8.014.74）g kg-1、（9.404.67）g kg-1（图 2）。土壤有机碳含量的平均值虽然有着随土壤深度增加而降低的趋势（R2=0.511），但在剖面中的分布并没有显著差异（P0.05）。研究区 6层土壤有机碳存量的平均值 标准差分别为（0.760.38）kg m-2、（0.760.50）kg m-2、（0.600.41）kg m-2、（0.570.30）kg m-2、（0.530.37）kg m-2、（3.172.34）kg m-2，其 中50100 cm的土层因为厚度的增加而具有较大的土壤有机碳存量。在研究区 050 cm 内的 5 个土层中，土壤有机碳存量有随深度增加而降低的趋势（R2=0.893），而且 1020 cm、4050 cm 两个土层的有机碳存量呈现出显著的差异（P0.05）。在 4 7005 100 m的海拔范围内，研究区土壤的有机碳含量随着海拔的升高而增加，然后
