2005–2015年中国生态系统研究网络农田生态系统生态站表层土壤有机碳密度数据集 中华文库

摘要：耕地有机碳是农田土壤碳库的重要组成部分。中国耕地土壤分布广泛，土壤的区域异质性高，利用方式和种植制度复杂多样，因此研究耕地有机碳及储量的变化规律，摸清耕地土壤的固碳能力，认识耕地固碳与生产力保持的耦合机制，对中国来讲尤其重要。国内很多学者利用全国第二次土壤普查数据、文献报道的数据开展相关研究，但总体上看，系统的、长期的耕地土壤观测与试验数据相对缺乏。本数据集利用中国生态系统研究网络19个涵盖不同农田生态系统的生态站多年系统的表层土壤理化性质实测数据，结合有机碳估算模型，计算出2005–2015年来中国主要农田生态系统生态站的表层土壤有机碳密度，可以为更准确地估算中国农田耕层土壤有机碳储量、研究其空间分布特征提供数据基础。

关键词：表层土壤有机碳密度；农田生态系统；中国生态系统研究网络；长期观测

Abstract: Cultivated land organic carbon is an important component of farmland soil carbon pool. Because of the wide distribution and the high soil heterogeneity of the cultivated land in China, the land use patterns and planting systems are complex and diverse. It is particularly important for China to study the changes of organic carbon and its storage in cultivated land, to understand the carbon sequestration capacity of cultivated land, and to understand the coupling mechanism between carbon sequestration and productivity maintenance of cultivated land. Many Chinese scholars use the data of the second National Soil Census and the data reported in the literature to carry out relevant research, but generally speaking, systematic and long-term observation data of cultivated land soil are relatively scarce. This dataset calculated the topsoil organic carbon density of ecological stations covering main agricultural ecosystems in China from 2005 to 2015 based on the measured data of physical and chemical properties of topsoil from 19 China Ecosystem Research Network (CERN) ecological stations covering different agricultural ecosystems, and combined with the estimation model of organic carbon. It will be the basis of estimating more accurately the soil organic carbon storage and studying its spatial distribution characteristics of farmland in China.

Keywords: topsoil organic carbon density; agricultural ecosystem; Chinese Ecosystem Research Network; long-term observation

土壤有机碳是地球陆地生态系统最重要和活跃的碳库, 同时又是土壤肥力和基础地力的最重要的物质基础。由于受到人类活动的强烈影响，农田土壤碳库的研究一直是碳循环研究的热点，也是全球变化、温室气体减排和粮食安全等问题研究的核心内容之一^[1][2]。耕地有机碳是农田土壤碳库的重要组成部分。在自然因素和农业管理措施的作用下，农田耕层土壤碳库在不断地变化，这种变化不仅改变土壤肥力，而且影响区域乃至全球碳循环^[3][4]。自然资源部《2016中国国土资源公报》显示，中国耕地面积为135 Mhm2左右，约占我国国土面积的14%。由于中国耕地土壤分布广泛，利用方式和种植制度复杂多样，土壤的区域异质性高^[5]，因此，研究耕地有机碳及储量的变化规律，摸清耕地土壤的固碳能力, 认识耕地固碳与生产力保持的耦合机制，对中国来讲尤其重要^[6][2]。

20世纪70年代末到80年代初开展的全国第二次土壤普查，获得了较为详细的全国范围内的土壤状况数据^[7]。国内很多学者以这些数据为基础，利用不同时期的耕地质量监测数据、采样数据、发表文献的数据^[4]，研究国家^[5]、区域^[8]、省域^[9][4]等不同尺度上耕地土壤有机碳时空分布、碳储量及固碳潜力。也有学者利用全国长期农田试验有机碳含量变化的文献资料、国内农田土壤有机碳和土壤有机质研究文献报道的农田耕层有机碳数据进行相关研究^[10][7][3]。但是总体上看，系统的、长期的耕地土壤观测与试验数据相对缺乏^[11]。

中国生态系统研究网络（Chinese Ecosystem Research Network，CERN）由中国科学院于1988年组建，目前有农田、森林、草地、荒漠、湖泊、海湾、沼泽、喀斯特和城市9类生态系统的44个生态站^[12][13]。其中，农田生态系统生态站涵盖了松嫩平原农田生态系统、下辽河平原农田生态系统、华北太行山前平原农田生态系统、黄淮海平原农田生态系统、太湖平原农田生态系统、南方丘岗区农田生态系统、洞庭湖平原农田生态系统、川中丘陵农田生态系统、黄土高原农田生态系统和荒漠绿洲农田生态系统^[14]。CERN制定了统一的大气、水分、土壤和生物要素变化观测指标体系，1998年以来，每个生态站根据指标体系，按规定的时间频度进行观测^[15][12]。农田生态系统土壤观测通过观测生态站所在地区的代表性农田生态系统主要类型土壤的组分、结构和重要生态过程的变化规律，为区域农田土壤质量的变化规律提供长期、系统的观测数据^[16]。土壤有机碳密度是表征土壤质量及陆地生态系统对全球变化贡献大小和衡量土壤中有机碳储量的重要指标^[8]。利用CERN农田生态系统生态站的土壤观测数据，计算各站多年表层土壤有机碳密度，可以为今后的农田土壤碳库研究提供数据支持。

耕作土壤剖面一般分作三层：表土层、心土层和底土层。其中最上部的表土层深度大致0～15（20）cm，是物质转化较快、干湿寒温变化较大的一层^[17]。本数据集以CERN农田生态站表层土壤理化性质实测数据为基础，结合有机碳估算模型，计算表层土壤有机碳密度。

1.1 基础数据

一般认为，土壤有机碳密度是由土壤有机碳含量、砾石（粒径＞2 mm）含量和容重共同决定。有研究认为，中国农业土壤经过长期人为的耕作熟化，耕层中粒径＞2 mm的砾石体积含量可能不是很大，可以忽略不计^[18][5]。因此本数据集主要收集了CERN农田生态站的表层土壤有机质与土壤容重实测数据以及采样地点地理位置、海拔、主要作物、管理方式等信息。数据全部来自“国家生态系统观测研究网络科技资源服务系统”（http://www.cnern.org.cn）中的生态系统长期监测数据集。

CERN在2002–2005年间系统地建立了土壤长期监测规范^[19]。为在国家和区域尺度上保证农田生态系统观测数据的可比性，并为区域农业的可持续发展提供决策依据，CERN农田生态站在综合观测场、辅助观测场采用相同的观测项目和观测方法进行土壤观测，以同时观测不同管理方式下农田生态过程的演化；各站还设置了站区调查点，以获取区域农户管理模式下土壤质量的变化。综合观测场设置在生态站所在区域内最具代表性的农田类型的典型地段，其土壤类型、种植模式与管理制度等具有典型性和代表性；辅助观测场是作为综合观测场的对照和补充，一般设置空白（不施肥）、化肥、化肥配合秸秆处理3种观测类型；站区调查点主要选择耕作、轮作及土壤类型与综合观测场一致或相近的、有代表性的农户地块^[20][16]。按照《CERN陆地生态系统土壤观测规范》，表层土壤有机质的观测频率为2～3年一次，有些生态站根据自身观测能力与条件，提高了观测频率；土壤容重的观测频率则为5年一次。本数据集选取了表层土壤有机质观测时间序列不少于5年的19个生态站综合观测场和辅助观测场的相关数据，时间范围为2005–2015年。19个生态站包含农田生态站14个，荒漠农田复合生态站5个（表1）。

表1 生态站基本信息

各指标测定方法如下：土壤容重全部采用环刀法测定；土壤有机质除海伦站、沈阳站外，均采用重铬酸钾氧化–外加热法测定，沈阳站主要采用元素分析仪测定法，海伦站土壤有机质测定方法具体情况见表2。

表2 海伦站历年土壤有机质测定方法汇总

1.2 数据处理方法

原始数据表中包含了表层和剖面的数据，其中表层土壤数据记录了各个样地不同采样区的数据；同时数据存在采样深度、观测频度等不一致的现象。因此在计算土壤碳密度之前，需要对数据进行预处理，计算出各样地表层土壤的有机质平均含量和平均容重。

1.2.1 土壤有机质数据的处理

首先需要从原始数据表中提取出表层土壤有机质数据。数据的选取分3种情况：（1）单一深度表层土壤有机质，多数生态站表层土壤深度为0～20 cm，有些生态站为0～15 cm；（2）个别生态站表层土壤有机质分为0～10 cm、10～20 cm，需要选取两个土层数据；（3）有些生态站个别年份仅有土壤剖面有机质，因此采用土壤剖面0～10 cm、10～20 cm土壤有机质作为替代值。筛选结果数据进行异常值剔除后，按年、月、作物类型计算各样地表层土壤有机质的算数平均值，得到多年样地尺度的表层土壤有机质数据。

1.2.2 土壤容重数据的处理

土壤容重数据与土壤有机质采用同样的方法处理，得到多年样地尺度的表层土壤容重数据。再根据样地、年、月、采样深度等与有机质数据进行匹配。由于表层土壤容重的观测频度低于有机质，因此多数年份的有机质缺少对应的土壤容重数据。对于缺失的土壤容重，本数据集采取以下方法确定：

（1）以该样地已有的相应深度的容重替代，原则是：2010年之前缺失容重的以2005年容重值替代，2010年与2015年之间缺失容重的以2010年容重值替代；

（2）对无容重数据的土层，根据该样地已有容重数据，建立其与有机质含量的回归关系^[21][22]，从而估算出缺失的容重值，并在备注中注明为“计算容重”。

1.3 土壤有机碳密度计算

土壤有机碳密度是指单位面积一定深度的土层中土壤有机碳的质量。鉴于数据实际情况，本数据集采用了不同的计算方式得到表层土壤有机碳密度。

单一土层（20cm或15cm）土壤有机碳密度的计算公式为^[23][21]：

SOC=0.58×OM×BD×H/100 （1）

式中，SOC为表层土壤有机碳密度（kg/m2）；0.58 为Bemmelen换算系数，该系数乘以有机质含量求得碳含量；OM为有机质含量（g/kg）；BD为土层容重（g/cm3）；H为土层厚度（cm）。

包含2个实测土层的表层土壤有机碳密度计算方法为：按照公式（1）分别计算0～10 cm、10～20 cm土壤有机碳密度，再将2个层次的土壤有机碳密度累加^[24][23][21]即得到0～20 cm的土壤有机碳密度。

本数据集的数据存储格式为Excel文件。文件包含2个数据表单，“土壤有机碳密度”表单存放19个农田生态站2005–2015年的表层土壤有机碳密度数据，共有748条记录；“样地信息”表单为数据涉及样地的描述信息。表3和表4列出2个数据表单所包含的具体字段名称、类型及示例。

表3 土壤有机碳密度表内容

表4 样地信息表内容

3.1 异常值检查

异常数据的判定统计学上有不同的方法，对于排除小样本容量的异常数据，Grubbs检验法目前被认为较其他方法更为有效^[25][26]，有学者应用该方法进行土壤数据异常值的剔除^[27][28]。本数据集采用Grubbs检验法对基础数据进行异常值的判断与剔除，具体方法如下：

对于每个样地的数据，按从小到大的顺序排列后计算平均值以及标准差，然后计算格拉布斯检验统计量；显著性水平取0.05，与格拉布斯表中的相应临界值进行比较，若计算的格拉布斯统计量大于格拉布斯临界值，则认为该数据为异常数据；将其剔除后重新计算余下的数据，直至确认无异常数据为止。

3.2 数据规范化

（1）样地名称与代码的规范化

样地名称以生态站简称加观测场分类表示，辅助观测场还需包含观测类型，如：常熟站综合观测场、常熟站土壤生物辅助观测场（空白）。样地代码采用7位编码，从左至右，第1–3位为站代码，4–5位为观测场分类码，第6–7位为观测场在本类观测场中的序号。站代码的编码规则为：前2位为生态站中文名称的汉语拼音缩写；第3位代表监测的主要生态系统类型，其中A表示农田生态系统，D表示荒漠生态系统。观测场分类码中，“ZH”表示综合观测场，“FZ”表示辅助观测场。例如，常熟站综合观测场的代码为“CSAZH01”、常熟站土壤生物辅助观测场（空白）的代码为“CSAFZ01”。

（2）要素项内容的规范整理

同一样地的采样分区描述、土壤类型、土壤母质、作物等要素项的内容，由于观测人员不同，造成不同年份的表述方式有所差异。因此分样地对这些要素项的描述进行规范化，使多年的记录内容保持一致，如：同一样地的土壤母质，有的年份记为“第四纪红黏土”，有的年份记为“第四纪红色黏土”，统一为“第四纪红黏土”。其中，根据《中国土壤分类与代码》（GB/T 17296-2009）对“土壤类型”要素项内容进行了重新划分，原则上填写土类，有些生态站填写的是亚类，但也予以保留，只是检查其名称是否符合标准。

本数据集以CERN多年系统的表层土壤理化性质实测数据为基础，结合有机碳估算模型，计算出近10年来中国主要农田生态系统系统生态站的表层土壤有机碳密度，具有较强的实际使用价值，可以为更准确地估算中国农田耕层土壤有机碳储量、研究其空间分布特征提供数据基础。

苏文. 2005–2015年中国生态系统研究网络农田生态系统生态站表层土壤有机碳密度数据集[DB/OL]. SDB, 2019. (2019-08-31). DOI: 10.11922/sciencedb.878.