基于土壤有机碳含量的黑土层厚度预测及影响因素分析
Organic carbon content-baesd prediction and influencing factors of black soil layer thicknesses
通讯作者: 宋运红(1983-),女,高级工程师,从事生态地质和土地质量调查研究综合研究工作。Email:yun-hong408@163.com
第一作者:
责任编辑: 蒋实
收稿日期: 2023-10-11 修回日期: 2024-01-8
基金资助: |
|
Received: 2023-10-11 Revised: 2024-01-8
黑土层厚度是黑土的一项基本属性,是衡量土壤肥力和侵蚀程度的重要指标,其空间预测的研究对支撑我国黑土地保护工程、保障粮食安全具有重要意义。本文参考土壤系统分类中黑土层的诊断特征,将有机碳含量高于成土母质6×10-3作为黑土层的判定标准,并利用有机碳在土壤垂向剖面上的指数分布规律,推算出黑土层厚度的计算公式。基于多目标区域地球化学调查获取的62 896个表层土壤和15 687个深层土壤的有机碳测试数据,对松辽平原黑土层厚度进行了详尽的空间预测分析,并探讨了黑土层厚度与土壤类型和气候因子之间的关系。 结果显示,松辽平原黑土层厚度在0~165 cm之间,中位数为23.33 cm。黑土层空间分布呈现出显著的非均质性,整体呈西南薄、东北厚的分布特点。沼泽土和泥炭土的黑土层平均厚度最大,在60~80 cm之间,其次为黑土,平均厚度为56 cm,白浆土和草甸土的黑土层平均厚度在40~50 cm之间。黑土层厚度的空间分布与气候条件关系紧密,主要表现为与温度呈显著的负相关,与降雨量呈正相关。同时,研究发现年均温0 ℃是影响黑土厚度的一个重要温度阈值,当年均温高于0 ℃时,黑土层平均厚度在80 cm以上,并且不再随温度发生变化。随着气候变暖,年均温0 ℃等温线的南移可能对黑土层厚度产生重要影响。
关键词:
Black soil layer thicknesses, anessential attribute of black soil,serves as a significant indicator for measuring the fertility and erosion degree of black soil. Their spatial prediction holds critical significance for supporting China's black land conservation and ensuring food security.Considering the diagnostic characteristics of black soil layers in soil system classification, this study regarded soil layers with organic carbon content higher than 6×10-3 of soil parent materials as black soil layers.Moreover, it derived the calculation formula for the thicknesses of black soil layers relying on the exponential distribution pattern of organic carbon in the vertical soil profile. Based on the 62 896 topsoil and 15 687 deepsoil organic carbon data obtained from the multi-purpose regional geochemical survey, this study conducted detailed spatial prediction of the thicknesses of black soil layers in the Songliao Plain and analyzed their relationship with soil types and climate factors. Key findings are as follows:(1) The thicknesses of black soil layers in the Songliao Plain range from 0 to 165 cm, with a median of 23.33 cm;(2) The spatial distribution of black soil layers exhibits significant heterogeneity, characterized by thin southwestern and thick northeastern portions;(3) The black soil layers of swampy soil and peat soil manifest the largest average thicknesses between 60 and 80 cm, followed by those of typical black soil (average thickness: 56 cm) and those of albic soil and meadow soil (average thickness: 40~50 cm);(4) The spatial distributions of the thicknesses of black soil layers are closely associated with climatic conditions, primarily showing a significant negative correlation with temperature and a positive correlation with rainfall;(5) The mean annual temperature of 0 ℃ is a significant temperature threshold for the development of thick black soil layers.Above this temperature, the average thickness of black soil layers exceeds 80 cm and no longer changes with temperature. With global warming, the southward shift of this 0 ℃ is otherm may significantly influence the thicknesses of black soil layers.
Keywords:
本文引用格式
刘凯, 戴慧敏, 刘国栋, 梁帅, 魏明辉, 杨泽, 宋运红.
LIU Kai, DAI Hui-Min, LIU Guo-Dong, LIANG Shuai, WEI Ming-Hui, YANG Ze, SONG Yun-Hong.
0 引言
黑土层厚度是黑土地的一项基本属性,也是衡量土壤肥力和侵蚀程度的重要指标[4-5]。目前关于土壤厚度制图的研究较多,而黑土层厚度的大尺度研究较少。土壤厚度图多根据实地调查的土壤厚度数据和环境变量进行模型预测,例如Wang等[6]利用地理统计学、地理信息系统和第二次全国土壤调查的1 627个土壤数据库,研究了中国土壤厚度的空间分布特征。Malone等[7]通过3个大型观测数据集和空间环境变量创建了澳大利亚国家土壤厚度图。也有研究人员利用探地雷达对土层厚度进行预测[8],但这种方法仅限于局部地区,难以推广到区域尺度。在黑土层时空变化方面,刘凯等[9]于2017年对第二次土壤普查的61个土壤剖面进行了重新调查,发现典型黑土区黑土层厚度在近40年内平均减少了12 cm,并认为气温升高和侵蚀作用是导致黑土层减薄的主要原因,但其研究范围仅限于典型黑土区,并且样品数量有限。
在土壤学中,目前并没有黑土层的明确定义。2022年,我国颁布的《中华人民共和国黑土地保护法》中,将黑土地定义为具有黑色或者暗黑色腐殖质表土层,性状好、肥力高的耕地。在土壤系统分类中,暗沃表层与黑土层的定义最为接近,可以作为黑土层判定标准的重要参考。暗沃表层的诊断特征包括厚度、颜色、有机碳含量、盐基饱和度、土壤结构等5个方面[10],其中关键之一是要求暗沃表层的有机碳含量≥6×10-3,同时,为了消除局部地区成土母质有机碳含量偏高的影响,还规定暗沃表层的有机碳含量应比成土母质高出6×1
本文利用多目标区域地球化学调查获取的松辽平原表层和深层土壤样品有机碳数据,根据黑土层有机碳含量的划分标准,尝试建立黑土层厚度的推算公式,并对松辽平原黑土层厚度的空间分布特征和影响因素进行分析。
1 研究区范围
图1
图1
松辽平原数据范围及气候条件
Fig.1
Data range and climatic conditions of the Songliao Plain
2 数据和方法
2.1 数据来源
土壤有机碳数据来源于多目标区域地球化学调查,采样时间为2003~2018年。表层土壤(0~20 cm)采集密度为1个样品/km2,深层土壤(150~200 cm)采集密度为1个样品/4 km2。样品在室温下风干后,用木槌压碎并通过20目尼龙筛,具体细节可参考文献[12]。表层土壤按4 km2大格、深层土壤按16 km2大格对格子内样品等质量组合成1个分析样品。每件组合样品都进行了有机碳、常量元素、微量元素等54种指标的测试分析,其中有机碳测试方法为重铬酸钾容量法,样品分析过程中进行了内部和外部控制,以确保测试质量[13]。在本研究中,共使用了62 896个表层土壤测试结果和15 687个深层土壤的测试数据,数据分布范围见图1a。
研究使用的气候数据(年均温和年降雨量)来自中国科学院资源环境科学与数据中心的中国气象背景数据集(
2.2 黑土层厚度计算方法
在平原区,土壤有机碳含量均表现为从表层至深层递减,且表层递减速率快,深层逐渐减慢的特点。多数研究表明,有机碳含量在垂向的变化通常符合指数函数模型,即:
式中:TOC为土壤有机碳含量;a、b为系数。该函数必定要通过本数据库中表层土壤(平均深度为10 cm)和深层土壤(平均深度是180 cm)两个点,因此将其代入指数函数:
经整理后,可计算出a和b:
将式(4)和式(5)代入式(1)中,并将黑土层的有机碳含量标准TOC标设定为TOC标=TOC深+0.6×10-2,计算所得的土壤深度即为黑土层厚度。经整理得到以下公式:
式中:TOC表为表层土壤有机碳含量,TOC深为深层土壤有机碳含量,单位均为%,黑土层厚度单位为cm。
为了验证公式的合理性,选取4种情况进行比较。在图2中,4条颜色曲线为不同情况下计算的有机碳含量与深度的指数分布曲线。A表示表层土壤有机碳TOC表=20×10-3,深层土壤有机碳TOC深=4×10-3,黑土层有机碳含量划定标准TOC标为10×10-3,计算的黑土层厚度为60.7 cm;B表示表层土壤有机碳TOC表=15×10-3,深层土壤有机碳TOC深=4×10-3,黑土层有机碳含量划定标准为10×10-3,计算的黑土层厚度为37.1 cm;C表示表层土壤有机碳TOC表=15×10-3,深层土壤有机碳TOC深=2×10-3,黑土层有机碳含量划定标准为8×10-3,计算的黑土层厚度为47.3 cm;D表示表层土壤有机碳TOC表=7×10-3,深层土壤有机碳TOC深=2×10-3,黑土层有机碳含量划定标准为8×10-3,计算的黑土层厚度为5.6 cm。A、B比较说明,成土母质的有机碳含量相同时,表层有机碳含量越高,黑土层越厚;B、C比较说明,表层有机碳含量相同时,成土母质有机碳含量越低,黑土层越厚。以上结果与黑土层判定标准和发育机理一致,说明公式是可靠的。而D的情况表明,当表层土壤有机碳含量与成土母质有机碳含量差值小于6×10-3时,计算的黑土层厚度会小于10 cm,我们将这种情况认定为不存在黑土层。
图2
图2
不同情况下土壤有机碳含量在垂向上的指数分布曲线
Fig.2
Exponential distribution curves of soil organic carbon content in the vertical direction under different conditions
具体操作如下:在ArcGIS中,将每个表层土壤点与对应位置的深层土壤点的有机碳属性进行连接,使每个表层土壤点具有TOC表和TOC深两个指标,代入公式(6),即可获得每个表层土壤点的黑土层厚度值。通过计算共获得松辽平原黑土层厚度值62 896个,随后在ArcGIS中利用克里金插值绘制黑土层厚度的空间分布图。
为判别黑土层厚度的影响因素,将黑土层厚度与气候因子进行相关分析和制图,相关操作在Orgin 2021中完成。
3 结果
3.1 表层和深层土壤有机碳含量空间分布
图3
图3
松辽平原表层土壤(a)和深层土壤(b)有机碳含量空间分布
Fig.3
Spatial distribution of organic carbon content in surface soils (a) and deep soils (b) of the Songliao Plain
3.2 黑土层厚度空间分布
图4
图4
黑土层厚度计算结果
a—黑土层厚度空间分布;b—黑土层厚度统计直方图;c—表层土壤有机碳含量与黑土层厚度散点图
Fig.4
Calculation results of black soil layer thickness
a—spatial distribution of black soil layer thickness; b—statistical histogram of black soil layer thickness; c—scatter plot of surface soil organic carbon content versus black soil layer thickness
图5
图5
不同类型土壤黑土层平均厚度统计
Fig.5
Statistical of the average thickness of the black soil layer in different types of soils
3.3 黑土层厚度与气候的关系
图6a、b显示了黑土层厚度与年均温和年降水量之间的关系。黑土层厚度与年均温之间整体上表现为明显的负相关关系,特别是当年均温在1~6 ℃之间时,黑土层厚度随着气温的升高而线性减小(R2=0.99,P≤0.01);而当年均温高于6 ℃或者低于1 ℃时,黑土层厚度较稳定,与温度关系不明显。黑土层厚度与年降雨量的关系相对复杂,在年降雨量350~530 mm内,黑土层厚度随年降雨量升高而增加;在年降雨量530 mm左右时,黑土层平均厚度达到58.37 cm;而随着降雨量继续增加,黑土层厚度呈波动性变化。
图6
图6
黑土层厚度与年均温和降雨量散点图
Fig.6
Scatter plot of black soil layer thickness and annual mean temperature and precipitation
由于温度和降水之间往往具有一定相关性,为了更直观地反映二者对黑土层厚度的影响,根据年均温和年降雨量的区间对所有样品进行整合。图6c显示,在同样的年降雨量范围内,黑土层厚度均随着年均温降低而增加。在年均温<0 ℃时,黑土层厚度均较大,在60 cm以上,此时黑土层厚度与年降雨量并没有明显关系,而当年均温>0 ℃后,黑土层厚度随年降雨量增加而呈一定的增加趋势。
4 讨论
4.1 黑土层厚度计算公式的合理性
本研究采用的黑土层计算公式有两个前提,一是选取指数函数作为有机碳含量在土壤垂向上的分布函数,二是将1.8 m深度的土壤有机碳含量认定为成土母质的有机碳含量。下面对这两个前提的合理性进行分析。
首先,众多学者对土壤有机碳含量的垂向分布规律开展过研究。不同地区土壤有机碳垂向特征并不一致,但在东北平原广泛分布的黑土、草甸土、黑钙土中,已有多项研究表明有机碳含量均呈现出相似的规律性,即由表层至深层含量递减,且表层递减速率快,深层递减速率减慢,直至有机碳含量稳定,整体上呈指数分布的特征[15-16]。同时,郝翔翔等[17]对黑土区农田、林地、裸地、草地等土壤有机碳含量的对比结果显示,这种指数分布规律并不以土地利用方式的不同而变化。奚小环等[18]对中国14个省市414条土壤剖面数据进行分析发现,有机碳分布状态存在地区性和沉积类型之间的差异,但垂向上仍以指数分布函数为主。在此基础上,众多关于土壤有机碳库的计算都是基于有机碳的指数模型而进行的[19]。因此,本文利用有机碳含量的指数分布模型来推算不同深度的有机碳含量是可行的。
4.2 黑土层厚度的空间分布影响因素
不同土壤类型的黑土层厚度具有一定差异性,沼泽土、黑土、草甸土的黑土层较厚,而风沙土、碱土、棕壤等土壤的黑土层偏薄,这与不同土壤的成土条件和成土过程有关。气候作为影响土壤形成的主要因素,在大尺度上的作用尤为重要[24],从图1和图4a中明显看出,黑土层厚度与气候条件具有较好的空间吻合度。图6a显示,温度对黑土层厚度具有显著的控制作用,是影响黑土层厚度空间分布的最主要因素。整体上看,年均温越高,黑土层厚度呈线性减少,这也是北部黑土层较厚,南部黑土层较薄的主要原因。图6c表明,在固定的温度范围内,降雨量会对黑土层厚度产生一定影响,即降雨量越高,黑土层发育的厚度越大,因此西部地区黑土层厚度普遍比东部要薄。
此外,本次研究中发现,当年均温低于0 ℃时,黑土层普遍较厚(平均值大于80 cm),并且不再随温度的继续降低而变化,这可能代表着年均温0 ℃是限制黑土层发育的一个温度阈值[32]。Yu等[33]对美国黑土区的一项研究也发现了同样的规律,即当年均温高于0 ℃时,土壤有机碳含量呈“陡崖式”下降,这主要反映了冻结条件对微生物分解活性具有重要的控制作用,从而影响了黑土层的发育程度。因此,年均温0 ℃是有利于厚层黑土发育的一个重要条件,而当年均温升高超过0 ℃时,土壤有机碳分解速率将显著加快,并可能导致黑土层厚度的迅速减薄。东北地区是我国气候变化最明显的地区之一,年均温呈现显著上升趋势,升温速率约0.38 ℃/10 a[34],明显高于华北、华南、西北等地区[35]。因此,随着气候变暖,年均温0 ℃等温线可能南移,这对黑土层厚度的未来变化和生态效应将是一个不容忽视的潜在影响因素,应予以重点关注。
5 结论
1)本文以土壤系统分类中暗沃土层的有机碳含量诊断特征,基于土壤地球化学调查获取的表层和深层土壤有机碳数据,利用有机碳在垂向上的指数分布函数,创建了黑土层厚度的计算公式,并预测了松辽平原黑土层的空间分布。
2)温度和降雨量是影响黑土层厚度空间分布的重要因素,黑土层厚度随着年均温的降低和降雨量的升高而变大。
3)年均温0 ℃是影响黑土厚层的重要温度阈值。随着东北地区气温升高,年均温0 ℃等温线的南迁可能对黑土层厚度产生重要影响。
参考文献
松辽平原典型黑土—古土壤剖面AMS14C年龄首次报道
[J].
The first report of the AMS14C age of mollisol-paleosol profile of Songliao Plain
[J].
东北地区黑土退化地球化学指示与退化强度
[J].
Black soil degradation and intensity in Northeast China:Geochemical indication
[J].
35年来东北松辽平原耕地土壤全氮时空变化
[J].
Spatio-temporal variation of total N content in farmland soil of Songliao Plain in Northeast China during the past 35 years
[J].
Quantitative assessment of soil productivity and predicted impacts of water erosion in the black soil region of northeastern China
[J].
Maximum humus horizon thickness as a criterion for identifying standard soils in the Crimean Plain
[J].
Characteristics of spatial variability of soil thickness in China
[J].
Improvements to the Australian national soil thickness map using an integrated data mining approach
[J].
GPR diagnostics of chernozem humus horizon thickness
[J].
东北黑土区黑土层厚度的时空变化
[J].
Spatiotemporal variation of black soil layer thickness in black soil region of Northeast China
[J].
Keys to soil taxonomy (12th ed.)
[R].
National multi-purpose regional geochemical survey in China
[J].
A complete set of analytical schemes and nanlytical data monitoring systems for determinations of 54 components in multi-purpose geochemical mapping
[J].
黑龙江省主要类型土壤剖面养分分布研究
[J].
Study on nutrient distribution of soil profiles of main soil type in Heilongjiang Province
[J].
不同利用方式对黑钙土腐殖质组分剖面分布特征的影响
[J].
Profile distribution characteristics of humus components in chernozem soil under different land use types
[J].
土地利用方式对黑土剖面有机碳分布及碳储量的影响
[J].以典型黑土区29年长期定位试验处理下的土壤为对象,研究了农田、裸地、自然草地和落叶松林地4种土地利用方式下土壤剖面(0~200 cm)有机碳及碳储量的分布特征.结果表明: 不同土地利用方式下表层(0~10 cm)土壤有机碳含量差异最大,表现为草地>农田>林地>裸地.农田10~120 cm各土层有机碳含量均低于草地、林地和裸地.与农田相比,自然草地对土壤有机碳提升作用明显,其0~60 cm各层土壤有机碳含量均显著高于农田;裸地表层(0~10 cm)土壤有机碳含量显著低于农田;落叶松林地0~20 cm有机碳含量与农田相比无明显变化,但其20~140 cm土层有机碳含量均高于农田.土壤剖面有机碳含量与pH值、容重、粉粒和粘粒含量呈显著负相关,与全氮和砂粒含量呈显著正相关.农田0~200 cm剖面有机碳储量显著低于其他3种利用方式,分别比草地、裸地和林地低13.6%、11.4%和10.9%.农田黑土在增加碳储量及改善环境方面具有很大潜力.
Profile distribution and storage of soil organic carbon in a black soil as affected by land use types
[J].
基于多目标区域地球化学调查的中国土壤碳储量计算方法研究
[J].土壤碳储量问题在碳循环和全球变化领域具有重要意义,长期以来这一基本参数的计算受到数据来源的制约。全国多目标区域地球化学调查采用双层网格化采样和分析,取得了大量高密度和高精度土壤地球化学数据,为土壤碳库的高精度计算提供了数据基础。文中重点探讨利用这些数据计算土壤碳储量的方法。首先提出针对多目标区域地球化学调查数据的“单位土壤碳量(USCA)”概念,用以代表调查数据基本面积单元(4 km2)的碳储量,并作为区域和全国尺度土壤碳储量计算的基本单位。在收集分析14个省市414条的土壤剖面数据的基础上,发现土壤有机碳(TOC)的垂直分布与无机碳和其他元素不同,存在指数分布规律,运用面积积分运算方法发现利用直线模型计算土壤有机碳库的误差(+20%~+100%)远大于指数模型的误差(-10%~+20%)。因此,文中建议 “有机碳单位土壤碳量(USCATOC)”需使用指数模型拟合后积分求算,而 “无机碳单位土壤碳量(USCATIC)”则使用直线模型, “全碳单位土壤碳量(USCATC)”采用两者加和计算。文中还分析了其他元素的垂直分布特征,并提出氮储量计算需采用与有机碳类似的方法,其他51种元素(氧化物)储量采用与无机碳类似的方案,按照直线模型计算单位土壤元素量和单位土壤氧化物量(USEA)。
Calculation techniques for soil carbon storage of China based on multi-purpose geochemical survey
[J].
松辽平原土壤碳库变化及其原因分析
[J].
Change in soil carbon pool in Songliao Plain and its cause analysis
[J].
Profile stock of soil organic carbon and distribution in croplands of Northeast China
[J].
Multiple exchange processes on mineral surfaces control the transport of dissolved organic matter through soil profiles
[J].
Dissolved carbon leaching from soil is a crucial component of the net ecosystem carbon balance
[J].
黑土开垦后黑土层厚度的变化
[J].
The thickness changes of ah horizon after the phaeozems cultivated
[J].
The age distribution of global soil carbon inferred from radiocarbon measurements
[J].
Soil organic carbon pools controlled by climate and geochemistry in tropical volcanic regions
[J].
Soil carbon storage controlled by interactions between geochemistry and climate
[J].
The contentious nature of soil organic matter
[J].
Global subsoil organic carbon turnover times dominantly controlled by soil properties rather than climate
[J].
DOI:10.1038/s41467-019-11597-9
PMID:31417092
[本文引用: 1]
Soil organic carbon (SOC) in the subsoil below 0.3 m accounts for the majority of total SOC and may be as sensitive to climate change as topsoil SOC. Here we map global SOC turnover times (τ) in the subsoil layer at 1 km resolution using observational databases. Global mean τ is estimated to be [Formula: see text] yr (mean with 95% confidence interval), and deserts and tundra show the shortest ([Formula: see text] yr) and longest ([Formula: see text] yr) τ respectively. Across the globe, mean τ ranges from 9 (the 5% quantile) to 6332 years (the 95% quantile). Temperature is the most important factor negatively affecting τ, but the overall effect of climate (including temperature and precipitation) is secondary compared with the overall effect of assessed soil properties (e.g., soil texture and pH). The high-resolution mapping of τ and the quantification of its controls provide a benchmark for diagnosing subsoil SOC dynamics under climate change.
Temperature sensitivity of soil carbon decomposition and feedbacks to climate change
[J].
Soil organic carbon stocks and dynamics in a mollisol region:A 1980s-2010s study
[J].
Persistence of soil organic matter as an ecosystem property
[J].
Climatic and geochemical controls on soil carbon at the continental scale:Interactions and thresholds
[J].
1961-2005年东北地区气温和降水变化趋势
[J].
Characteristics of temperature and precipitation in Northeastern China from 1961 to 2005
[J].
/
〈 |
|
〉 |
