2012–2016 年三江源垂穗披碱草人工草地碳水热通量观测数据集

贺福全，李奇，陈懂懂，赵亮*

1.中国科学院西北高原生物研究所，西宁 810008

2.中国科学院三江源国家公园研究院，西宁 810008

引 言

三江源被誉为“中华水塔”，是我国重要的水源地。明晰碳、水、热交换对气候变化的响应方式和碳、水循环的机制是三江源区生物多样性维持和保育的核心，长期定位观测是探明碳、水循环的重要手段。三江源地处青藏高原东南部，源区生物多样性极丰富、特有物种多，是高原生物多样性最集中的地区，也是全球气候变化敏感区。三江源区草地面积约27.59×104km2，占全区总面积的71.2%[1-2]。在过去的50 年中，由于气候突变和人为活动增加，青藏高原超过50%草地退化或沙化[3-4]，在极度退化的“黑土型”退化草地上建植人工草地是可行的，也是快速恢复植被的主要途径[5-6]。为了恢复退化草地和发展畜牧业，近年来三江源区人工草地面积不断扩大，目前约为67000 hm2[7]。

目前对青藏高原碳水热通量的研究主要集中在天然草地[8-17]，生态系统的源汇功能被认为主要受生长季的降水量、强度和分配的影响[9,12]。有研究通过比较不同植被类型CO2通量，指出CO2通量是由自下而上的气候因子（温度、降水）和自上而下的生物因子（群落结构、叶面积和放牧）两条途径控制[18]。三江源站自2009 年开始，在三江源东部农牧交错区同德牧场开展垂穗披碱草人工草地碳水热通量的观测，为评估源区人工草地碳源、汇功能和碳储量提供数据支撑。

为了推动三江源乃至全青藏高原草地生态系统物质循环和能量交换等相关研究的发展，三江源站公开发表多年连续观测的碳水热通量及相关常规气象数据，以供更多科研工作者使用数据，更大发挥数据价值。本数据集包含2012–2016 年垂穗披碱草人工草地的常规气象数据（空气温度、空气相对湿度、土壤温度、土壤湿度、光合有效辐射和降水）和碳水热通量数据（净生态系统CO2交换通量、生态系统CO2呼吸通量、总生态系统CO2交换通量、潜热通量、显热通量），形成了30 分钟、日、月和年4 种时间尺度的数据产品。

1 数据采集和处理方法

1.1 数据来源

三江源站（34°21′ 17″ N, 100°28′48″ E；海拔3950 m）位于青海省果洛州玛沁县大武镇。站区为典型的高原大陆性气候，无明显四季之分，年均气温和降水分别为-2.9 ℃和531.0 mm。垂穗披碱草人工草地涡度相关观测位于三江源站东北部约200 km 的海南州同德县同德牧场（35°15′11″N，100°41′57″E，海拔为3305 m）[19]。

垂穗披碱草人工草地涡度相关系统建于2009 年9 月，由30 m×30 m 的围栏保护，碳水热通量的观测高度为2.2 m。该区域属典型高原大陆性气候，辐射强烈，无绝对无霜期，年均温为3.7 ℃，年降水量为427.2 mm，降水主要集中在5–9 月，占全年的87%，气温和降水的峰值均出现在8 月前后，雨热同期，土壤类型为暗棕壤[19]。人工草地于2008 年建成，由翻耕草带和未翻耕草带组成，翻耕草带和未翻耕草带宽均4 m，相间分布，全年不放牧。翻耕草带每5 年重新翻耕（5 年为一个种植周期），草带为单播的垂穗披碱草（Elymus nutans），在6 月初种植，种植当年不收割牧草，其他年份8 月底收割牧草，全年禁牧。未翻耕草带主要物种为兰石草（Lancea tibetica）、紫花苜蓿（Medicagoruthenica）、冷蒿（Artemisia frigida）和垂穗披碱草[19-20]。2011 年6 月种植的垂穗披碱草因仪器故障数据缺测。故2012 年为种植第2 年，2013 年为种植第3 年，2014 年为种植第4 年，2015 年为种植第5 年，2016 年6 月重新翻耕种植[19-20]。由于2011、2017 年涡度相关系统红外气体分析仪工作异常时段较多，数据可用性较低，三江源站没有共享2011 年和2017 年的数据，2017 年之后数据经过进一步处理质控后发布共享。根据多年的观测和文献总结[12,14,20], 结合种植和收割时间，该人工草地生长季确定为4 月中旬至10 月下旬。

1.2 数据采集方法

本数据集中包含的观测数据均是仪器自动观测采集的，测定所用仪器型号、仪器制造商及数据采集器等相关信息详见表1。空气温度、湿度、光合有效辐射、土壤温度、土壤含水量等气象数据的采样频率1 分钟，计算并存储30 分钟的平均数据；CO2/H2O 快速红外气体分析仪及三维超声风速仪的原始数据采样频率为10 Hz，计算并存储30 分钟的平均数据。

表1 三江源站垂穗披碱草人工草地观测仪器信息表Table 1 Table of the information about the observation instruments for the Elymus nutans artificial grassland at Sanjiangyuan Station

1.3 数据处理和产品加工方法

为了保障三江源站垂穗披碱草人工草地通量数据可靠，采用了ChinaFLUX 推荐的通量数据处理的标准方法[21]，具体包括原始数据异常值剔除、超声虚温校正、坐标轴二次旋转、WPL 密度校正、频率损失校正、湍流稳态测试、夜间摩擦风速阈值筛选和通量异常值剔除，以及能量闭合评价。由于仪器故障和数据质量控制导致数据的缺失，需要对缺失数据进行插补。对缺失2 小时以内的气象和通量数据，采用线性内插法进行插补；气象数据缺失数据大于2 小时的利用平均日变化进行插补；通量数据缺失大于2 小时的采用线性回归方程进行拟合、插补，最小插补时间窗口为7 天。对于长时间缺失的CO2通量数据，首先利用有效通量数据和环境因子的非线性回归关系构建方程，然后利用该方程和缺失通量数据对应的环境因子完成数据插补[19-20]。对缺失的白天CO2通量数据采用时间窗口的方程（1）进行插补[22]：

RE(d)日间生态系统呼吸，a是生态系统表观量子效率，Pmax是饱和光合成速率，PPFD是光量子密度通量。

夜间缺失CO2通量及白天RE 通量采用方程（2）插补[23]：

RE(n)为夜间生态系统呼吸，REref为土壤温度为10 ℃时的参考生态系统呼吸速率，Ts为5 cm 深处的土壤温度，Swc为5 cm 深处的土壤含水量 (m3/m3)，a、b、c、d 是适合全年数据的特定地点参数，其中b > 0，d ≤ 0。

CO2通量拆分。涡度相关技术测得的CO2通量是净生态系统CO2交换量（NEE）。利用夜间的有效观测数据和方程（2）可获得RE 与土壤温度(Ts)和土壤含水量(Swc)之间的关系函数，再利用此函数估算生长季白天的生态系统呼吸RE(d)。日RE 为白天呼吸RE(d) 和夜间呼吸RE(n)之和。NEE 减去RE，则可得到生态系统总光合固碳量（GEE）方程（3）：

2 数据样本描述

2.1 数据子集与数据量

本数据集为2012–2016 年三江源站垂穗披碱草人工草地的连续碳水热通量观测数据，分为常规气象数据子集和通量数据子集两类数据文件。每类文件每年均有4 个时间尺度，包括30 分钟、日、月和年尺度，数据集内共计40 个文件，总数据量39.8 MB。

2.2 数据文件示例

以2012 年30 分钟数据文件为例，表2、表3 分别为30 分钟常规气象数据和碳水热通量数据表头说明，所有数据项观测数据均以浮点型数字形式表示。

表2 三江源站垂穗披碱草人工草地常规气象30 分钟数据表说明及观测高度Table 2 Description of conventional half-hourly meteorological data and observation height of the Elymus nutans artificial grassland at Sanjiangyuan Station

表3 三江源站垂穗披碱草人工草地通量30 分钟数据数据表说明Table 3 Description of the half-hourly flux data of the Elymus nutans artificial grassland at Sanjiangyuan Station

3 数据质量控制与评估

能量闭合程度分析是检验通量数据质量的有效手段，虽然能量不闭合现象普遍存在，但大部分被人们广泛接受。能量平衡方程：H + LE = Rn - G。其中，Hs 为显热通量，LE 为潜热通量，Rn 为净辐射，G 为土壤热通量，土壤热通量是两块热通量板平均值。冠层热储量忽略不计。能量闭合方程斜率分别为0.76（2012 年），0.76（2013 年），0.73（2014 年），0.72（2015 年）和0.89（2016 年）[19]。生长季能量平衡方程斜率（0.76～0.9）大于非生长季（0.65～0.71），能量闭合程度在可接受的范围之内 (0.7～0.9)[24-25]。导致能量不闭合的原因非常复杂，可能的原因包括观测系统的采样误差、仪器偏差、湍流通量观测中的高低频损失、平流效应等[26]。

4 数据使用方法和建议

本数据集可用于垂穗披碱草人工草地碳水热生态过程及相关陆面模型开发、验证及对比等分析研究。本数据集采用了ChinaFLUX 推荐的通量数据处理方法[21-22]，与青藏高原的ChinaFLUX 站点对比使用，开展区域尺度的联网分析。需要说明的是由于插值方法不同导致碳水热通量计算结果存在差异，即使年际尺度上通量计算结果相似，也可能在日、月尺度上存在一定的差异，其中30 分钟数据为未插补的原始数据，日、月、年尺度的数据是经过插补计算后的数据，数据使用者根据需要进一步判定，也可利用气象数据进行自行插补使用。

数据作者分工职责

贺福全（1993—），男，青海省大通县人，研究生，工程师，研究方向为高寒草地碳循环。主要承担工作：数据整理和论文撰写。

李奇（1983—），男，甘肃省兰州市人，研究生，高级工程师，研究方向为草地生态学。主要承担工作：数据监测与论文修改。

陈懂懂（1982—），女，山东省东营市人，研究生，高级工程师，研究方向为高寒草地碳循环。主要承担工作：数据分析与论文修改。

赵亮（1974—），男，青海省西宁市人，大学，研究员，研究方向为高寒草地碳循环。主要承担工作：数据管理和统筹。