2007–2013年栾城典型灌溉农田生态系统日通量数据集

中国生态系统研究网络CERN专题 I 区论文（评审中） • 版本 ZH1

Daily fluxes dataset of the typical irrigated agro-ecosystem in Luancheng (2007 – 2013)

张玉翠，姜寒冰，张传伟，沈彦俊

Zhang Yucui, Jiang Hanbing, Zhang Chuanwei, Shen Yanjun

DOI: 10.11922/csdata.2020.0002.zh
PID: 21.86101.1/csdata.2020.0002.zh
: 10.11922/sciencedb.939

： 2020 - 01 - 14

： 2020 - 01 - 17

443 3 0

摘要&关键词

摘要：基于微气象学原理的涡度相关技术已成为陆地生态系统与大气间CO₂和水热通量交换的标准观测手段之一。本数据集为中国科学院栾城农业生态系统试验站自涡度相关系统架设完成至今积累的全部通量数据，时间跨度从2007年10月至2013年9月，覆盖6个连续作物年，观测样地为冬小麦–夏玉米一年两熟轮作制农田，属于华北地区典型的潮褐土高产农业生态类型，具有较强的代表性。本数据集的联网观测、质量控制和处理存储过程均严格遵守ChinaFLUX数据管理技术体系，数据可靠性高，可为华北地区典型农田的水平衡、适水种植制度调整、农田生态系统碳水循环过程以及作物模型模拟等相关研究提供坚实的数据支撑。

关键词：栾城；涡度相关系统；灌溉农田生态系统；水–热–碳日通量；冬小麦–夏玉米

Abstract & Keywords

Abstract: Based on micrometeorology, the Eddy Covariance System (EC) has become one of the standard observation methods for the exchange of CO₂ and water heat flux between terrestrial ecosystems and the atmosphere. The dataset includes all the fluxes data accumulated in Luancheng Agro-ecosystem Experimental Station of the Chinese Academy of Sciences from the EC System, which spans from October 2007 to September 2013 and covers 6 consecutive crop growth years. And the fluxes consist of evapotranspiration(ET), precipitation(P), net radiation (Rn), sensible heat Flux (Hs) and net ecosystem exchange of CO₂(NEE). The observation object is winter wheat-summer maize rotation cropping system. It can represent the general crop growth system in the North China Plain (NCP). The network observation, quality control and storage process of this dataset strictly abide by ChinaFLUX data management technology system to make sure the reliability. The EC flux data provide important support for the study of crop field water balance and typical cropping system adjustment in the North China Plain, the mechanism of ecosystem carbon-water coupling cycles and crop growth model simulation.

Keywords: Luancheng; eddy covariance system; irrigated agro-ecosystem; water-heat-carbon daily fluxes; winter wheat-summer maize cropping system

数据库（集）基本信息简介

数据库（集）名称	2007–2013年栾城典型灌溉农田生态系统日通量数据集
数据作者	张玉翠、姜寒冰、张传伟、沈彦俊
数据通信作者	沈彦俊（yjshen@sjziam.ac.cn）
数据时间范围	2007年10月1日至2013年9月30日
地理区域	观测站点：中国科学院栾城农业生态系统试验站(37°53′N，114°41′E)；代表区域：太行山山前平原（东经112°30′–119°30′，北纬34°46′–40°25′）。
数据量	295 KB
数据格式	*.xlsx
数据服务系统网址	http://www.sciencedb.cn/dataSet/handle/939
基金项目	国家生态系统观测研究共享服务平台专题服务；国家自然科学基金项目（40871021、31400375、41471027、31870422、41930865）；国家重点研发计划课题（2016YFC0401403）；中国科学院青年创新促进会项目（2017138）。
数据库（集）组成	本数据集由栾城试验站通量观测数据组成，含1个Excel文件，包括5个表单，分别为农田蒸散量、降水量、净辐射量、感热通量、碳通量数据记录。每个表单分别包括2192条日尺度数据和6条生育季尺度数据。

Dataset Profile

Title	Daily fluxes dataset of the typical irrigated agro-ecosystem in Luancheng (2007 – 2013)
Data authors	Zhang Yucui, Jiang Hanbing, Zhang Chuanwei, Shen Yanjun
Data corresponding author	Shen Yanjun (yjshen@sjziam.ac.cn)
Time range	From October 1, 2007 to September 30, 2013
Geographical scope	Luancheng Agro-ecosystem Experimental Station, CAS (37°53′N，114°41′E); geographical scope represented: The Taihang piedmont plain of North China (112°30′ E – 119°30′E, 34°46′ – 40°25′).
Data volume	295 KB
Data format	*.xlsx
Data service system	<http://www.sciencedb.cn/dataSet/handle/939>
Source(s) of funding	Special Services of National Ecosystem Observation and Research Sharing Service Platform; National Natural Science Found of China (40871021, 31400375, 41471027, 31870422, 41930865); National key research and development program (2016YFC0401403); The Youth Innovation Promotion Association of Chinese Academy of Sciences (2017138).
Dataset composition	The dataset consists of the flux and meteorological measurements of Luancheng Agro-ecosystem Experimental Station. This Excel file includes five sheets: Evapotranspiration (ET), Precipitation (P), Net Radiation (Rn), Sensible Heat Flux (Hs), Net ecosystem exchange of CO₂(NEE), respectively. Each sheet consists of 2192 daily data and 6 growth seasonal scale data.

引言

陆地生态系统碳水循环是地球化学循环的核心过程。研究陆地生态系统碳水循环过程对于人类调节生物圈–水圈–大气圈相互作用，维持全球生态系统的物质运输与能量循环具有重要的指导意义，也是人为适应和干预全球变暖，缓解水资源供需矛盾的战略需求^[1-2]。经过长期的理论发展和技术改革，基于微气象学原理的涡度相关技术已成为陆地生态系统与大气间的CO₂和水热通量交换的标准观测手段之一^[2]。涡度相关系统的观测原理就是将空气流动视为无数个涡流，通过物理量的脉动与垂直风速脉动的协方差计算湍流通量^[3]。目前，涡度相关系统已被广泛应用于陆地生态系统碳水循环过程研究，长期连续的高频率联网观测积累的通量数据可为碳水循环机理、碳源碳汇的时空分布特征以及生物地球化学循环模型的优化等相关研究提供有力的数据支撑^[4-5]。

2001年，中国陆地生态系统通量观测研究网络（ChinaFLUX）依托中国生态系统研究网络（CERN）成立，经过十几年的发展，ChinaFLUX目前已建设完成100余个通量观测站，覆盖森林、草地、农田等多种类型的陆地生态系统，形成了全国尺度的碳水通量观测研究体系，实现了我国通量观测事业从起步到国际前沿水平的跨越式发展^{[6,7,8 ]}。作为ChinaFLUX重要成员之一的中国科学院栾城农业生态系统试验站（以下简称“栾城试验站”）始建于1981年，隶属于中国科学院遗传与发育生物学研究所农业资源研究中心（简称“中科院农业资源研究中心”）。1989年加入中国生态系统研究网络（CERN），成为第一批基本站。2005年成为国家生态系统观测研究网络（CNERN）台站，同年加入中国通量观测研究联盟，成为ChinaFLUX农田生态系统的重要观测站点之一。

栾城试验站涡度相关系统长期连续的碳水通量联网观测开始于2007年，本数据集包含了栾城试验站2007–2013年的碳水通量观测数据，选取了其中较为重要的组成部分，包括农田蒸散量、降雨量、净辐射量、感热通量和碳通量，整理形成了日尺度和生育季节尺度数据产品，并以此数据论文形式对数据集进行系统性描述，以便于读者溯源定位，理解数据集内涵，推动数据开放共享和规范化使用。

1 数据采集和处理方法

1.1 数据来源及采集方法

基于ChinaFLUX的顶层设计，经过观测塔选址、观测仪器选型和野外观测系统安装与调试等技术工作，中科院栾城农业生态系统试验站于2007年11月6日正式开始长期连续的碳水通量联网观测。本数据集为2007年10月1日至2013年9月30日的农田碳水通量数据，其中2007年10月1日至11月5日的数据由历年相应通量项平均值并结合栾城试验站气象数据、土壤含水量和参考蒸散等插补得到。

栾城试验站全称为中国科学院栾城农业生态系统试验站（台站代码：LCA），该站地处太行山山前平原中部(37°53′N、114°41′E，海拔50.1 m)，暖温带半湿润季风气候，年平均气温12.9℃，多年平均降雨量490 mm（2007–2013年），无霜期200天，全年平均日照时数2608小时，≥0℃的积温为4710℃，雨热同期，适宜农业发展。土壤以潮褐土为主，质地为壤土，随土层由浅到深依次为粉砂壤、壤土和粘壤。农业生产以华北地区种植面积最为广泛的冬小麦–夏玉米一年两熟轮作制农田为主。冬小麦10月初播种，次年6月中旬收获；夏玉米6月中旬播种，9月底收获。其中冬小麦冠层高度约为1.0 m，夏玉米冠层高度在2.0 m左右。

结合上述植被类型、冠层高度以及通量观测下垫面要求，栾城站架设了通量观测塔，并将碳水分析仪安装高度设置为3.5 m。涡度相关系统主要由开路式CO₂/H₂O红外分析仪和三维超声风速计组成，同时辅有辐射分量和常规气象要素（降雨和空气温湿度等）的同步观测，观测设备型号及制造商等信息见表1。数据观测和存储除人工观测降水量外由系统自动完成，通常情况下每0.5–1个月定期维护一次，每两年校准一次。

表1 各观测系统主要测定要素、关键传感器型号、制造商信息及安装高度

观测系统	观测要素	传感器和分析仪	制造商	安装高度（m）
气象要素	自动降水量	HMP155	CAMPBELL (USA)	3.5
气象要素	人工降水量	20 cm直径雨量筒+1000 ml 量筒	-
CO₂和水热通量	三维超声风速	CSAT3	CAMPBELL (USA)
	CO₂/H₂O密度	LICOR7500	LI-COR (USA)
	净辐射量	CNR-1	Kipp&Zonen (Netherlands)
数据采集与通讯	常规气象要素	CR1000	CAMPBELL (USA)
数据采集与通讯	碳水通量要素	CR5000	CAMPBELL (USA)

通量观测塔下垫面为充分灌溉管理下的典型冬小麦–夏玉米一年两熟轮作制农田生态系统，下垫面均匀且地势平坦，风浪区足够大（fetch＞350 m）。观测系统利用CR5000型数据采集器以10 Hz的频率进行原始数据采样和存储，连续自动监测生态系统的碳水通量，每30 min输出一组平均值。常规气象要素的观测记录输出频度同样为30 min，由相应的数据采集系统自动完成数据获取和存储。通过U盘下载原始数据并存储备份，然后根据ChinaFLUX数据管理要求进行后续的质量控制和标准化处理。

1.2 数据处理方法

通量数据的原始观测数据统一遵循ChinaFLUX技术体系进行标准化的质量控制和处理，处理流程如图1。为规范和便于读者对数据的使用，本数据集将10 Hz的原始数据处理计算形成了日尺度和生育季节尺度数据产品^[9]。

图1 通量数据质量控制与处理流程

数据质量控制过程主要包括原始数据分析、超声虚温校正、坐标轴旋转、WPL校正、频率损失校正、野点剔除、稳态测试、湍流积分特性、夜间摩擦风速阈值筛选以及能量闭合度评价^{[10,11,12,13,14,15 ]}。野点剔除指原始数据中的异常值剔除，剔除方法包括绝对值法和相对值法^[16]。绝对值法即依据观测指标理论上可能出现的数值范围设置上下限，将明显超出正常范围的数值直接剔除；相对值法即计算数据的平均值和方差，以4倍方差作为检验标准，剔除偏离过大的数据。坐标轴旋转也称为倾斜校正，其目的是保证超声风速仪与地表平行，减少由于仪器倾斜导致观测的垂直风速中包含水平风速分量而造成的通量数据误差。本研究采用中倾斜校正，主要包括三次坐标轴旋转（Tiple Coordinate Rotation，TR）、WPL（Webb-Pearman-Leuning）校正、频率响应校正。其中三次坐标轴旋转可使x轴平行于平均风速方向，侧向平均风速与垂直平均风速及两者的协方差为零；WPL校正针对由于水热通量输送引起气体密度变化而导致的潜热通量和感热通量观测误差；频率响应校正针对仪器的系统误差，即由于仪器对高频和低频信号响应存在的缺陷。根据平均值检验法（Average Values Test，AVT）对夜间CO₂通量数据进行校正^[17]。本数据集前期主要利用编写的程序代码（2007–2011年）以及软件EddyPro（LI-COR公司提供，2011–2013年）对通量数据进行了以上质量控制和处理，2017年后期利用ChinaFLUX统一程序代码对原始数据又进行了一次处理，由于所选计算方法基本一致，数据一致性达99.5%。针对缺失数据的插补，短时段（小于2 h）的缺失数据采用线性内插法，较长时段（大于2 h）的缺失数据利用平均昼夜变化法（mean diurnal variation，MDV）进行插补^[18]。对于蒸散量（ET）连续多日（7日以上）的数据缺失时段或MDV插值后仍存在的异常值，我们根据蒸散发生的机制对数据处理方法进行了改进，利用土壤含水量与ET/ET₀的定量关系进行重新插补，这也是我们对数据插补方法的一点创新。具体方法如下：

①利用联合国粮农组织（food and agricultural organization，FAO）推荐的Penman-Montieth 公式计算参考作物腾发量ET₀^[19] ：

\({ET}_{0}=\frac{0.408∆\left({R}_{n}-G\right)+900U*\gamma *VPD/\left({T}_{a}+273.3\right)}{∆+\gamma \left(1+0.34U\right)}\) （1）

式中：ET₀ 为潜在蒸散发(mm)，R_n 为作物表层净辐射(MJ/(m²·day))，G为土壤热通量(MJ/(m²·day))，T_a 日平均气温(^oC)，U为2 m高度处风速(m/s)，VPD为饱和水汽压差(kPa)，Δ为饱和水汽压曲线斜率(kPa/^oC)，γ为干湿表常数(kPa/^oC)。

②筛选出连续7日以上的数据缺失时段或异常值时段。

③筛选出存在缺值的日期及其前后2–3天的数据。

④从筛选结果中选择与缺值处土壤含水量相近的日期（日期相同，年份不同），将各ET/ET₀ 取平均值作为缺值日的ET/ET₀ 。其中土壤含水量数据由IH-II型中子仪和自动传感器共同测定，观测深度为0–2 m。

⑤将平均所得的ET/ET₀ 乘以当天对应的ET₀ ，得到插补后的ET。

2 数据样本描述

2.1 数据集命名规则及数据量

本数据集为栾城站2007年10月1日至2013年9月30日的农田生态系统日尺度碳水通量观测数据，共计1个EXCEL数据文件，大小为295 KB。文件命名为“栾城典型灌溉农田生态系统2007–2013年日通量数据集”，其中包括5个表单，分别命名为“蒸散量”“降水量”“净辐射量”“感热通量”和“碳通量”。每个表单分别包括2192条日尺度数据和6条生育季尺度数据，其中降雨量数据包括了同时期的人工观测和自动观测数据。

2.2 数据文件示例

表2为数据表头示例，详细介绍了数据项含义、数据类型、计量单位等信息。

表2 通量观测数据表头说明

数据项	数据类型	计量单位	数据项说明	示例
Date	日期	——	年/月/日	2007/10/1
Year	日期	——	年-次年（表示生育年）	2007–2008
ET	数字	mm	日尺度蒸散量	1.3
ET-wheat	数字	mm	冬小麦生育季蒸散量	400.2
ET-maize	数字	mm	夏玉米生育季蒸散量	299.3
ET-Year	数字	mm	生育年尺度蒸散量	699.5
P-M	数字	mm	人工观测日尺度降水量	6.2
P-A	数字	mm	自动观测日尺度降水量	6.2
P-wheat	数字	mm	冬小麦生育季降水量	240.2
P-maize	数字	mm	夏玉米生育季降水量	357.3
P-Year	数字	mm	生育年尺度降水量	597.5
Rn	数字	MJ/m²/day	日尺度净辐射量	5.9
Rn-wheat	数字	MJ/m²	冬小麦季净辐射量	1232.2
Rn-maize	数字	MJ/m²	夏玉米季净辐射量	1003.0
Rn-Year	数字	MJ/m²/year	生育年尺度净辐射量	2235.2
Hs	数字	MJ/m²/day	日尺度感热通量	2.4
Hs-wheat	数字	MJ/m²	冬小麦季感热通量	270.7
Hs-maize	数字	MJ/m²	夏玉米季感热通量	202.6
Hs-Year	数字	MJ/m²/year	生育年尺度感热通量	473.3
NEE	数字	gC/m²/day	日尺度CO₂通量	0.9
NEE-wheat	数字	gC/m²	冬小麦季CO₂通量	−309.2
NEE-maize	数字	gC/m²	夏玉米季CO₂通量	−393.1
NEE-Year	数字	gC/m²/year	生育年尺度CO₂通量	−702.3

数据表头说明：

（1）Year表示生育年，10月1日至次年9月30日，即从小麦播种到玉米收获的一个轮作周期。

（2）P（Precipitation）表示降水量（mm），表示从天空降落到地面上的液态或固态（经融化后）水，未经蒸发、渗透、流失，而在水平面上积聚的深度。

（3）ET（Evapotranspiration）表示农田蒸散量（mm），指单位时间单位面积农田土壤蒸发量和植物蒸腾量的总和。

（3）Rn（Net Radiation）表示净辐射量（MJ/m²），是反映太阳辐射强弱的物理量。

（4）Hs（Sensible Heat Flux）：表示感热通量（MJ/m²），指由于温度变化引起的大气与下垫面之间发生的湍流形式的热交换。

（5）NEE（Net ecosystem exchange of CO₂）表示CO₂通量（gC/m²）。NEE为正值表示农田生态系统向大气释放CO₂；NEE为负值表示农田生态系统从大气吸收CO₂。

3 数据质量控制和评估

目前全球通量观测研究领域普遍使用的数据质量评价体系主要包括一贯性检验和完全湍流假设的检验。本数据集基于此对数据质量进行了系统评价。功率谱和协方差谱检验结果表明，三维风速、CO₂浓度和H₂O浓度的功率谱变化模态符合惯性副区－2/3 斜率理论值，表明涡度相关系统设备响应正常，不存在系统性的相移或失真^[20]。结合湍流稳定性检验和积分统计特性检验对数据进行总体质量评价和等级划分，结果显示观测样本符合湍流的方差相似性规律，数据质量较高^[21]。能量闭合度是通量数据质量评价的重要参考标准之一，Wilson等通过统计分析了全球通量观测网络（FLUXNET）22个站点的50年通量数据，发现所有植被类型和气候条件下的观测均存在不同程度的能量不闭合现象，平均不闭合度达20%^[22]。为反映生态系统通量实际变化情况，本数据集中通量观测数据未进行强制能量闭合处理。但是借助于改进的ET（LE）插补方法，能量闭合度可达91%。其中小麦季观测结果的能量闭合度可达95%，玉米季可达88%，远高于全球通量观测网络台站汇交数据能量闭合程度的平均水平，数据质量较好。此外，为保证降水数据质量，本数据集不仅包含了人工观测降水数据，也收录了同时期的自动观测降水数据。

本数据集中，日尺度蒸散量和碳通量有效观测数据为90%和93%；人工降水量有效观测数据为100%，自动观测为96%；净辐射量和感热通量有效观测数据约为95%。数据缺失的原因主要分为两类：一是观测仪器运行故障导致的数据缺失，包括供电故障、设备维护等；二是数据处理过程导致的数据缺失，比如夜间和降水天气数据的筛选等。

4 数据使用方法和建议

本数据集为栾城试验站自通量观测塔架设完成至今积累的全部通量数据，其观测、处理和质量控制与评估均采用国际通用方法，并根据自身条件和观测情况进行了改进，可靠性高。数据跨度从2007年10月至2013年9月，覆盖6个连续作物年，观测样地属于华北地区典型的潮褐土高产农业生态类型，具有较强的代表性。本数据集可为农业节水理论与水文水资源研究等相关领域的发展提供坚实的数据支撑，包括农田耗水特征与蒸散结构分离^[23-24]，土壤水利用层次与地下水补给途径解析^[23,25-26]，农田生态系统碳氮水循环过程^[27-28]，农田水平衡和适水种植制度调整^{[16-17,29-30]}，农业生产力与水资源可持续利用程度评价^[31]，作物模型的验证与改进^[29,30,32]等。此外，本数据集可与其他观测台站的数据进行综合集成，服务于陆地生态系统对全球变化的响应、区域尺度陆地生态系统物质循环与能量流动过程以及生态系统管理政策制定等相关领域。

目前通量观测数据仍普遍存在不同程度的能量不闭合现象，数据处理和质量控制也存在一定的不确定性。因此，本数据集在使用过程中需注意以下几个方面：

（1）随着ChinaFLUX数据管理技术体系的更新和完善，不同年份汇交的通量数据处理方法不完全相同，因此本数据集可能与早期发表的数据存在一定差异（误差允许范围内）。

（2）本数据集的数据管理和质量控制方法详细信息可参考Zhang等^{[16, 26]}、Shen等^[17]、于贵瑞等^[33]、Wen等^[20]和Yu等^[6]已发表的文献。

（3）本数据集各项观测数据受仪器运行状态或数据处理过程影响，存在不同程度的数据缺失。数据文件中对插补数据进行了标注，建议优先采用未插补数据，以减少不确定性。

（4）针对降水数据，人工观测可能会在降水量较大的连续降水日无法及时观测导致数据偏小，而自动观测可能对小雨0.2 mm的降水难以捕捉，有时也会由于仪器问题出现缺测，这是目前降水观测普遍存在的问题，所以建议根据研究实际情况选择所需降雨数据。

致谢

衷心感谢CERN综合中心和ChinaFLUX数据资源管理人员以及栾城站观测人员的支持！感谢中科院地理所张雷明、陈智老师在数据质量控制和处理程序方面给予的指导！感谢中科院农业资源研究中心王玉英、程一松、闵雷雷、齐永青和郭英博士在数据收集、处理与监管方面给予的帮助！

[1]

GU F, ZHANG Y, HUANG M, et al. Climate-driven uncertainties in modeling terrestrial ecosystem net primary productivity in China[J]. Agricultural & Forest Meteorology, 2017, 246: 123-132.