三江源冻土区浅层土壤有机碳及其氧化稳定性特征
1.
2.
3.
Organic carbon and its oxidation stability characteristics of shallow soil in frozen soil area of Three River Source Region
1.
2.
3.
通讯作者:
收稿日期: 2021-11-30 修回日期: 2022-08-23
基金资助: |
|
Received: 2021-11-30 Revised: 2022-08-23
关键词:
Keywords:
本文引用格式
孙丰豪, 张志春, 蒋平安, 杨凯, 周少龙.
SUN Fenghao, ZHANG Zhichun, JIANG Ping’an, YANG Kai, ZHOU Shaolong.
0 引言
冻土是温度低于0 ℃并出现冻结现象的土壤(或地表疏松层)。活动层则是随季节转换而发生周期性冻融的冻土类型,其年平均气温高于0 ℃,受土壤性质影响较大[1-2]。三江源气候属青藏高原气候系统,年平均气温低,无明显四季之分,浅层土壤(0~40 cm)是随冷暖季节转换而发生冻融的活动层。冻土区浅层土壤有机碳的数量与质量关系到季节转换时土壤与大气之间的热量交换,影响着冻土的发育、保护和退化。土壤有机碳氧化稳定性能够反映成土条件及过程[3],是衡量土壤有机碳抗氧化能力的重要指标。氧化稳定性系数(Kos)是土壤有机碳氧化稳定性的量化指标,其数值越大土壤有机碳就越稳定。土壤有机碳氧化稳定性的变化主要发生在易氧化分解的部分,不同氧化程度有机碳的占比能直观地反映出土壤有机碳分解的难易程度[4]。根据土壤有机碳氧化的难易程度可划分为四个不同组分(Fraction 1~4):高氧化活性有机碳(F1)、中氧化活性有机碳(F2)、低氧化活性有机碳(F3)和难氧化活性有机碳(F4)[5]。各组分的多寡共同决定着土壤有机碳氧化稳定性的大小。
三江源平均海拔在4 000 m以上,冻土分布广泛,自西向东土壤类型逐渐由草甸土转变为草原土,土壤有机碳含量较高[6],特殊的地理位置及气候条件使得冻土区浅层土壤有机碳对环境变化极为敏感[7-9]。复杂的大尺度区域研究,样点稀疏且高寒地区恶劣的自然环境又加大了工作难度,使得土壤样品采集困难且容易出现误差,传统土壤制图难以对有机碳特征做出精准反映[10]。数字土壤制图基于环境协变量来定量土壤景观模型和栅格数据,能克服传统土壤调查制图仅利用斑块的局限性,多种因素综合考虑,在一定程度上可弱化外界条件影响并减少误差,较准确反映出土壤性质的变化,已在青藏高原农耕区广泛应用于碳、氮分布等研究[11-16]。本文采用随机森林模型对三江源冻土区浅层土壤有机碳氧化稳定性分布特征进行解析,探索其分布规律和对环境的响应。
1 研究区概况
2011—2019年间全区年均气温为0.9 ℃,年均降水量475.50 mm,年均日照时数为2 622.62 h,气温低、降水少但日照充足。年均气温、年均降水量呈自东南向西北逐渐降低的特征,年均日照时数呈自东南向西北逐渐升高的特征(图1)。
图1
图1
2011—2019年三江源地区水热、日照条件
Fig. 1
Water, heat and sunshine conditions in Three River Source Region from 2011 to 2019: average annual air temperature (a), precipitation (b) and sunshine hours (c)
2 材料与方法
2.1 样品采集与测定
采样点布设以中国1∶1 000万冻土区划及类型图、1∶100万土壤图、青海省土壤普查资料、青海省行政区划图和青海省交通图为依据,在能代表当地及植被类型的平坦地区布设采样点采样,样点共86个,间距在30~100 km之间(图2)。
图2
图2
三江源地区土壤采样点空间分布
Fig. 2
Spatial distribution of soil sampling points in Three River Source Region
确定采样点后,在能够代表本地植被特征且相对平坦区域按100 m×100 m的样方布置样点,按“X”形5点采样,采集0~40 cm或至冰碛物的土壤。样品捡去石块等杂物后,用四分法保留1 kg布袋封装带回实验室阴凉处自然风干。将风干土再次混匀,以四分法分出150~200 g,仔细拣出枯枝落叶、砾石、植物根系等杂物后,研磨过0.25 mm筛待测。
对于土壤有机碳氧化稳定性的测定,普遍采用强酸加不同浓度氧化剂分别在一定温度下氧化土壤,由此来获得有机碳不同的氧化量,常用的氧化剂有K2Cr2O7和KMnO4。本文采用不同浓度梯度的硫酸-重铬酸钾溶液在不同温度下氧化土壤(对同一样品分别称取5份相同的量)获取相应条件下土壤有机碳的氧化量,确定各级氧化活性组别的绝对量,据此计算各组别的相对含量和氧化稳定性系数。用Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ表示不同的氧化强度,序号越高氧化能力越弱[19]。不同的氧化强度条件如下:
Ⅰ 0.40 mol 1/6 K2Cr2O7同1∶1 H2SO4,170~180 ℃条件下于油浴煮沸5 min;
Ⅱ 0.27 mol 1/6 K2Cr2O7同1∶2 H2SO4,140~150 ℃条件下于油浴煮沸5 min;
Ⅲ 0.20 mol 1/6 K2Cr2O7同1∶3 H2SO4,130~140 ℃条件下于油浴煮沸5 min;
Ⅳ 0.13 mol 1/6 K2Cr2O7同1∶5 H2SO4,120~130 ℃条件下于油浴煮沸5 min;
Ⅴ 0.073 mol 1/6 K2Cr2O7同1∶10 H2SO4,110~120 ℃条件下于油浴煮沸 5 min。
2.2 制图参数
2.3 结果评估
用goof函数进行交叉验证验证得出决定系数(R2)、一致性系数(CC)、平均误差(ME)和均方根误差(RMSE)作为模型预测结果评估指标,公式[24]为
式中:N为建模样点的个数;
3 结果与分析
3.1 样品数据统计
采样结果表明,三江源冻土区TOC在5.27~157.35 g·kg-1之间,均值为(33.00±17.51) g·kg-1。不同氧化程度有机碳含量差异显著(P<0.05),均值F4最高,F3均值最低,分别为18.15 g·kg-1、3.88 g·kg-1,氧化稳定性系数(Kos)在0.44~2.92之间,均值为1.41(表1)。
表1 采样点土壤有机碳(TOC、F1~F4、Kos)
Table 1
土壤有机碳 | 均值 | 最小值 | 最大值 | 方差 | 偏度 | 峰度 |
---|---|---|---|---|---|---|
TOC | 33.00 | 5.27 | 157.35 | 337.38 | -0.64 | 0.13 |
F1 | 5.12 | 0.31 | 18.77 | 8.92 | -0.36 | 0.36 |
F2 | 6.35 | 0.70 | 21.36 | 15.11 | -1.12 | 0.17 |
F3 | 3.88 | 0.15 | 13.68 | 4.48 | -0.51 | 0.55 |
F4 | 18.15 | 1.25 | 105.31 | 102.94 | -0.51 | 0.23 |
Kos | 1.41 | 0.44 | 2.92 | 0.14 | 1.38 | 1.25 |
3.2 环境变量相关性分析
表2 采样点土壤有机碳(TOC、F1~F4、Kos)与环境变量的相关系数
Table 2
环境变量 | 土壤有机碳 | |||||
---|---|---|---|---|---|---|
TOC | F1 | F2 | F3 | F4 | Kos | |
年均降水量 | 0.51** | 0.44** | 0.46** | 0.35** | 0.47** | -0.15* |
年均日照时数 | -0.50** | -0.43** | -0.48** | -0.35** | -0.45** | 0.15* |
年均气温 | 0.36** | 0.29** | 0.43** | 0.17* | 0.33** | -0.16* |
海拔 | -0.31** | -0.22** | -0.31** | -0.21** | -0.29** | 0.08 |
3.3 制图结果评估
预测结果中决定系数(R2)均在0.54以上,模型对上述六种有机碳属性有54%的解释度。一致性系数(CC)均大于0.58,预测值与实际值一致性较高。除TOC预测结果均方根误差(RMSE)为14.02、F4预测结果均方根误差(RMSE)为9.18之外,其余预测结果均方根误差皆小于2.90。随机森林模型能较好地对土壤有机碳进行预测,并对其变异性有较为精准的反映(表3)。
表3 土壤有机碳(TOC、F1~F4、Kos)预测值交叉验证结果
Table 3
土壤有机碳 | 一致性系数(CC) | 平均误差(ME) | 均方根误差(RMSE) | 决定系数(R2) |
---|---|---|---|---|
TOC | 0.79 | -0.26 | 14.02 | 0.69 |
F1 | 0.78 | -0.05 | 2.43 | 0.67 |
F2 | 0.81 | 0.01 | 2.90 | 0.73 |
F3 | 0.78 | -0.07 | 1.84 | 0.72 |
F4 | 0.76 | -0.18 | 9.18 | 0.66 |
Kos | 0.58 | -0.01 | 0.57 | 0.54 |
3.4 土壤有机碳制图结果
TOC经度地带性特征明显,高值区分布在三江源东部(100°00′~102°23′ E)和中部(96°30′~97°30′ E),西部(90°00′~95°00′ E)低值区分布广泛。
F1相对含量(F1/TOC,下同)在9.75%~24.03%之间,高值区分布在东北一隅(34°50′~36°10′ N,100°00′~101°27′ E)并向西递减至低值区(32°30′~33°25′ N,98°20′~101°27′ E),中部(95°00′~97°30′ E)F1占比适中,在12.93%~19.27%之间,西部(90°00′~95°00′ E)多为低值区分布。
F2相对含量在5.34%~30.62%之间,同F1相近,高值区分布在东北一隅,但向西F2占比降幅度小于F1,中部(95°00′~97°30′ E)F2占比适中,西部偏低,自东向西逐渐递减。
F3相对含量在4.15%~30.21%之间,与F1、F2不同,西部(90°0′~94°40′ E)F3高值区广泛分布,中部(94°40′~98°30′ E)低值区广泛分布,中值区在东部(98°30′~101°27′ E)广泛分布且东北一隅(34°50′~36°10′ N,100°00′~101°27′ E)F1、F2的高值区反而是F3的中值区。
F4相对含量在34.71%~65.19%之间,同F1、F2相近,仅在东北一隅存在小范围高值区,依次向西递减,低值区分布广泛。
Kos在0.77~2.60之间,东北一隅(32°30′~35°25′ N,100°20′~101°27′ E)和西北(34°30′~35°25′ N,92°30′~94°47′ E)存在部分高值区,低值区偏向南部分布。北高南低,地带性特征并不突出(图3)。
图3
图3
三江源土壤有机碳(TOC、F1/TOC~F4/TOC、Kos)空间分布
Fig. 3
Spatial distribution of soil organic carbon (TOC, F1/TOC~F4/TOC, Kos) in Three River Source Region
3.5 环境变量贡献率分析
节点纯度(IncNodePurity)是随机森林用于回归过程中残差的平方和,定量地反映了环境因子在预测过程中所占权重,其数值越大贡献率越高,越小则反之。
图4
图4
土壤有机碳(TOC、F1~F4、Kos)预测结果的节点纯度
Fig. 4
Nodal purity of soil organic carbon (TOC, F1~F4, Kos) prediction results
4 讨论
不同氧化难易程度组分影响着土壤有机碳的理化性质,其与土壤有机碳含量的占比反映了土壤有机碳的组成和状态[33-35]。与高寒地区低温少雨特性有关,三江源冻土区浅层土壤难氧化活性有机碳组分占比(F4/TOC)最高,而低氧化活性有机碳组分占比(F3/TOC)最低,前者对土壤有机碳的氧化稳定性的影响最大,后者最小。全区自东南向西北随着海拔升高、气候的变化,土壤类型逐渐由草甸土向草原土过渡,植被覆盖度也随之降低,同时海拔等因素也对局部小气候造成不可忽视的影响。在种种因素共同作用下,不同氧化难易程组分有着不同的分布特征。F1/TOC、F2/TOC、F4/TOC自东向西递减地带性特征明显,TOC阶梯性递减,中部地区偶有升高;水热光照条件良好且草甸土广泛分布的东北一隅不仅是TOC的高值区,也是F1/TOC、F2/TOC和F4/TOC较高的区域。气候寒冷干燥、植被覆盖度低且以草原土为主导的西部地区TOC、F1/TOC、F2/TOC和F4/TOC都较低,F3/TOC却较高。中部TOC高值区有机碳组分分布表现出了较大差异(F3相对含量更低),其原因可能是中部高寒草甸、草原镶嵌分布,植被种类相对丰富,加快了低氧化活性有机碳的转化,也正是这种各组分间的分布差异导致了Kos和TOC在大尺度范围内变异性的产生。
作为具有较为综合性质的土壤有机碳氧化稳定性系数(Kos),数字制图结果显示全区Kos处于较高水平。三江源低温少雨的气候特征导致土壤有机碳活性较低、分解较慢,同时也意味着冷、干的环境有利于冻土区浅层土壤有机碳的长期赋存。
5 结论
基于实验数据与2011—2019年气候资料,采用随机森林模型对三江源冻土区浅层土壤有机碳氧化稳定性进行了多要素数字化制图,得出以下结论:
(1)三江源冻土区浅层土壤有机碳受降水和日照时数影响较大,温度对其作用有限。
(2)三江源冻土区浅层土壤有机碳不同氧化难易程度组分分布各异,除低氧化活性有机碳组分(F3)外,其余组分随水热条件的地带性变化而变化。
(3)三江源冻土区浅层土壤有机碳氧化稳定性系数总体呈北高南低的特征,其地带性特征表征力度不及各不同氧化难易程度组分。
参考文献
Study on water vapor conditions of precipitation on the Tibetan Plateau during rainy season
[J].
青藏高原雨季降水的水汽条件研究
[J].
The space-time variation analyses of air temperature over the Qinghai-Xizang Plateau in warm and cold seasons
[J].
青藏高原暖季与冷季气温的时空演变分析
[J].
Labile soil organic matter as influenced by cropping practices in an arid environment
[J].
Storage of organic carbon in aggregate and density fractions of silty soils under different types of land use
[J].
Oxidizable organic carbon fractions and soil quality changes in an Oxic Paleustalf under different pasture leys
[J].
Vegetation coverage change and its response to climate change in Three-River Headwaters Region
[J].
三江源地区植被覆盖度变化及其与气候因子的关系
[J].
Distribution characteristics and influencing factors of soil organic carbon in alpine ecosystems on Tibetan Plateau transect
[J].
青藏高原样带高寒生态系统土壤有机碳分布及其影响因子
[J].
Changes in climate, glaciers and permafrost during 1979—2018 in Tibetan Autonomous Region and their potential impacts on renewable energy
[J].
1979—2018年西藏自治区气候与冰川冻土变化及其对可再生能源的潜在影响
[J].
Soil organic carbon storage in typical regions of China
[J].
中国典型地区土壤碳储量研究
[J].
Soil information system of China (SISChina) and its application
[J].
中国土壤信息系统(SISChina)及其应用基础研究
[J].
Review of digital soil mapping in soil nutrients
[J].
数字土壤制图在土壤养分方面的研究综述
[J].
Operationalizing soil security by making soil data useful: digital soil mapping, assessment and return-on-investment
[J].
Chapter 1 spatial soil information systems and spatial soil inference systems: perspectives for digital soil mapping
[J].
Spatio-temporal variation of soil total nitrogen in Huangshui River basin and its affecting factors in the past 30 years
[J].
近30年湟水流域土壤全氮时空变异及影响因素
[J].
Comparison on soil depth prediction among different spatial interpolation methods: a case study in the Three-River Headwaters Region of Qinghai Province
[J].
土壤厚度的空间插值方法比较: 以青海三江源地区为例
[J].
Research frame of vulnerability and adaptation for the cryosphere and its changes
[J].
冰冻圈及其变化的脆弱性与适应研究体系
[J].
Research on the establishment of national “northwest Sichuan national ecological protection and ecological economic development comprehensive pilot area”: based on the establishment of “Qinghai Sanjiangyuan national ecological protection comprehensive pilot area” and the investigation and analysis of desertification control in northwest Sichuan
[J].
创建国家级“四川西北国家生态保护与生态经济发展综合实验区”的研究: 基于“青海三江源国家生态保护综合实验区”的建立及四川西北地区沙化治理的调研分析
[J].
Biodiversity decreasing of Three-River Source Area and countermeasures for control
[J].
三江源生物多样性减少及防治对策
[J].
A critical examination of a rapid method for determining organic carbon in soils: effect of variations in digestion conditions and of inorganic soil constituents
[J].
Soil readily oxidizable carbon and carbon pool management index in spruce plantation (Picea asperata) with different thinning intensity in western Sichuan
[J].
间伐强度对川西亚高山人工云杉林土壤易氧化碳及碳库管理指数的影响
[J].
Soil organic carbon pools and productivity relationships for a 34 year old rice-wheat-jute agroecosystem under different fertilizer treatments
[J].
Research on urban fire risk evaluation based on the random forest algorithm
[J].
基于随机森林算法的城市火灾风险评估研究
[J].
Robust and accurate shape model matching using random forest regression-voting
[J].
Performance evaluation of classification algorithms by k-fold and leave-one-out cross validation
[J].
Digital mapping of topsoil attributes in Qinghai Province
[J].
青海省表层土壤属性数字制图
[J].
The method of selecting an integrated development territory for the high-rise unique constructions
[J].
Hyper-scale digital soil mapping and soil formation analysis
[J].
The review and outlook of digital soil mapping
[J].
数字土壤制图研究综述与展望
[J].
Prediction of soil total nitrogen for rubber plantation at regional scale based on environmental variables and random forest approach
[J].
基于多源环境变量和随机森林的橡胶园土壤全氮含量预测
[J].
Construction of membership functions for predictive soil mapping under fuzzy logic
[J].
Study on the oxidative stability of soil humus
[J].
土壤腐殖质氧化稳定性的研究
[J].
Soil organic carbon characteristics of different alpine grasslands in Qinghai-Tibet Plateau
[J].
青藏高原不同类型草地土壤有机碳特征研究
[J].
Advances in researches on soil organic carbon storages, affecting factors and its environmental effects
[J].
土壤有机碳储量、影响因素及其环境效应的研究进展
[J].
Characteristics of SOC and ROC under different land use types at arid valley in the upper reaches of Minjiang River
[J].
岷江上游干旱河谷不同土地利用类型的土壤有机碳和易氧化态碳特征
[J].
Distribution of oxidizable organic C fractions in soils under cacao agroforestry systems in Southern Bahia, Brazil
[J].
Comparative analysis of the soil organic carbon oxidation stability in Altai Mountain, Xinjiang
[J].
阿尔泰山区土壤有机碳氧化稳定性的初步比较分析
[J].
/
〈 |
|
〉 |
