张法伟，李红琴，赵亮，张雷明，陈智，祝景彬，徐世晓，杨永胜，赵新全，于贵瑞，李英年＊

1.中国科学院西北高原生物研究所高原生物适应与进化重点实验室，西宁810008

2.洛阳师范学院生命科学学院，洛阳 471934

3.中国科学院地理科学与资源研究所生态系统网络观测与模拟重点实验室，北京100101

引 言

青藏高原是中国“青藏高原生态屏障”的空间载体，被誉为“中华水塔”。高寒湿地是青藏高原的重要植被类型[1]，面积约为13.3×106hm2，在碳素固持、水源涵养及生物多样性维持方面具有重要意义[2]。长期的科学观测数据是认知高寒湿地碳、水、热交换格局和内在生态过程及关键生态功能评估的基础。涡度相关法是目前唯一能直接测定大气与下界面的碳水热通量的标准方法，被国内及国际通量观测组织广泛采用[3]。青海海北高寒草地生态系统国家野外科学观测研究站（简称海北站）是国内较早利用涡度相关技术开展碳水热通量观测的台站之一，自2000 年陆续开展了高寒矮嵩草（Kobresia humilis）草甸、高寒金露梅（Potentilla fruticosa）灌丛和高寒帕米尔苔草（Carexpamirensis）湿地的碳水热通量观测，积累了大量的原始观测数据，填补了青藏高原通量观测研究的空白，取得了丰硕的研究成果，为理解青藏高原高寒陆地生态系统碳水功能及区域可持续发展提供了翔实的数据支撑和重要贡献[4-14]。

为了促进青藏高原生态系统碳水热交换过程的科学认知及生态功能的精准评估，体现数据共享理念和挖掘数据潜在价值，海北站拟免费公开发表连续观测的高寒湿地碳水热通量及常规气象数据。由于2010 年高寒湿地涡度相关系统的快速红外气体分析仪工作异常时段较多，导致该年度的通量有效观测数据比例偏低，尤其是净生态系统CO2交换量和显热通量分别仅为8.2%和10.0%，数据可利用较低，海北站没有共享该年度通量数据。因此碳水热通量数据子集和常规气象数据子集的年份分别为2004-2009 年和2004-2010 年。本数据集包含碳水热通量数据子集（净生态系统CO2交换通量、生态系统CO2呼吸通量、总生态系统CO2交换通量、潜热通量、显热通量）和常规气象数据子集（空气温度、空气相对湿度、水汽压、风速、风向、土壤温度、总辐射、光合有效辐射和降水等观测指标），时间尺度包括半小时、日、月和年。

1 数据采集和处理方法

1.1 数据来源

海北站（37°37′N、101°19′E，3200 m）位于青海省海北州门源县种马场风匣口，多年平均气温和降水分别为−1.7℃和580 mm，站区为典型的高原大陆性气候，暖季温暖多雨而短暂，冷季寒冷干燥而漫长。高寒湿地涡度相关观测系统位于海北站东部约2 km 的高寒帕米尔苔草湿地生态系统内（37°36′N、101°19′E，3200 m），始建于2003 年10 月，是中国通量观测研究网络（ChinaFLUX）的首批通量观测站点之一。涡度相关系统的观测区地势平坦，可以满足通量观测“风浪区”的要求。研究区植被生长茂盛，分布均匀，外貌整齐，平均总盖度约98%，冠层高度为25-50 cm。植被种类组成较为贫乏，约有25 种植物，多以寒冷湿生、多年生地下芽、具有发达通气组织的草本植物莎草科、毛茛科为主。湿地中央以帕米尔苔草为建群种，边缘以藏嵩草（K.tibetica）为建群种。湿地群落的主要优势种为帕米尔苔草、藏嵩草，次优势种有华扁穗草（Blysmus sinocompressus），黑褐苔草（C.atrofusca）等。在边缘地带还有大量的星状风毛菊（Saussurea stilla）、青藏苔草（C.moorcropt）、乳白香青（Anaphalis lectea）、假龙胆（Gentianella limprichtii）、天山报春（Primula nutans Georgi）和斑唇马先蒿（Pedicularis longiflora Rudolph var）等[11]。土壤类型为有机寒冻潜育土（Organic Cryic Gleysols）。研究区为冬季牧场，放牧强度约为每公顷4.0 羊单位，放牧时间约为10 月下旬至翌年5月上旬[12]。湿地植被群落2003-2005 年平均最大地上生物量（干物质）和最大叶面积指数分别为8月下旬的331.3 g⋅m-2和7 月下旬的3.9 m2⋅m-2[15]。

1.2 数据采集方法

高寒湿地碳水热通量数据集的数据全部为仪器自动观测并存储在相应数据采集器中，具体的观测仪器及型号、仪器制造商及数据采集器等相关信息详见表1。碳水热通量原始数据的采样频率为10 Hz，计算并存储其30 min 的平均通量数据于CR5000 中；常规气象要素的采样频率为1 min，计算并存储30 min 的平均数据于CR23X 中。

表1 海北高寒湿地的观测项目所用分析仪相关信息

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

为了保障数据质量及数据的可对比度，本数据集从数据观测、数据采集、质量控制、处理方法等流程严格遵循ChinaFLUX 制定的标准化数据处理和质量控制技术体系[16-17]。具体内容如下：

数据质量控制：碳水热通量的数据质量控制主要包括原始数据异常值剔除、超声虚温校正、坐标轴二次旋转、WPL 密度校正、频率损失校正、冠层储存项校正、湍流稳态测试；半小时通量数据的质量控制包括夜间摩擦风速阈值筛选、异常值剔除和降水时刻数据剔除以及能量闭合评价。用单点密度估算的方法表明冠层存储项通量分别约占CO2通量、H2O 通量和显热通量的1.1%、0.09%和2.2%，因此本数据集的通量数据没有进行储存项校正。常规气象数据的质量控制主要为阈值判别和异常值剔除。

缺失数据插补：由于通量观测仪器自身的稳定性和异常天气的影响及其他不可控因素，野外观测数据的缺失十分普遍，缺失数据的插补十分必要[18]。数据插补根据数据子集类别和缺失数据长短分别采用不同的方法进行插补。其中，短时间（小于2 小时）缺失的通量和常规气象数据，采用线性内插法进行插补；长时间缺失的气象数据，基于海北站的自动气象站同期的观测资料（降水数据除外）进行插补；如未能完成插补，则利用平均日变化法完成相应数据的插补。

对于长时间缺失的水热通量数据，采用水热通量与净辐射的线性回归关系进行拟合、插补。对于长时间缺失的碳通量数据，采用有效通量数据和环境因子的非线性经验方程进行拟合、插补。其中夜间缺失数据利用CO2通量有效数据与5 cm 土壤温度的Arrhenius 方程插补。白天缺失数据分两种情况进行插补，在植被生长季中利用CO2通量数据与光合有效辐射的直角双曲线方程插补。在非生长季中，首先采用夜间CO2通量数据和5 cm 土壤温度的形成Arrhenius 方程，然后结合缺失数据对应的土壤温度完成白天通量的插补。碳水热缺失通量数据的最小插补时间窗口为7 天。

CO2通量数据拆分：基于边际分布采样法的原理将观测的净生态系统CO2通量数据拆分为生态系统CO2呼吸通量和总生态系统CO2交换通量。首先，依据夜间CO2通量数据与土壤温度拟合的Arrhenius 方程估算白天的生态系统CO2呼吸通量（生态系统暗呼吸通量）；其次，利用插补完成的白天CO2通量数据和估算的同时刻白天生态系统CO2呼吸通量，求和得到总生态系统CO2交换通量（生态系统初级生产力）[17]。生态系统CO2呼吸通量为生态系统夜间呼吸通量和白天呼吸通量之和。

2 数据样本描述

2.1 数据子集与数据量

本数据集为高寒湿地碳水热通量观测数据，分别为2004-2009 年通量数据子集和2004-2010 年常规气象数据子集两类数据文件，每类数据文件均有半小时、日、月和年4 个时间尺度。因此，高寒湿地碳水热通量数据集内含52 个文件，总数据量34 MB。

2.2 数据文件示例

以2004 年数据文件为例，表2、表3 分别为海北站高寒湿地碳水热通量和常规气象数据表头说明，所有数据项观测数据均以浮点型数字形式表示。

表2 海北站高寒湿地通量数据表说明

注：该表为30 min 通量数据的数据表。日尺度数据表中无时、分等2 列数据项，通量数据为30 min 通量累计值，NEE、RE 和GEE 的单位为g C m-2 d-1,LE 和Hs 的单位为MW m-2；月尺度数据表中无日、时、分等3 列数据项，通量数据为日通量累计值，NEE、RE 和GEE 的单位为g C m-2 month-1,LE 和Hs 的单位为MW m-2；年尺度数据表中没有月、日、时、分等4 列数据项，通量数据为月通量累计值，NEE、RE 和GEE 的单位为g C m-2 year-1,LE 和Hs 的单位为MW m-2。

表3 海北高寒湿地常规气象数据表说明及指标观测高度

注：日尺度数据表中无时、分、秒3 列数据项，月尺度数据表中无日、时、分、秒4 列数据项，年尺度数据表中没有月、日、时、分、秒5 列数据项。

3 数据质量控制和评估

2004-2009 年，经过30 min 数据质量控制和异常值剔除后，高寒湿地数据集净生态系统CO2交换量、潜热通量、显热通量的平均有效观测数据比例分别为42.4%、52.7%和55.0%，其中净生态系统CO2交换通量有效观测数据比例的最低年份和最高年份分别为2007年的34.1%和2004年的48.3%，两者之差为14.2%，而潜热通量和显热通量有效观测数据比例的最低年份与最高年份分别为2009 年的44.2%和47.8%与2008 年的63.9%和67.4%，两者相差约为10%（表4）。

表4 海北高寒湿地半小时碳水热通量的有效观测数据比例

4 数据使用方法和建议

本数据集在国家科技资源共享服务平台的国家生态科学数据中心（http://www.cnern.org.cn/data/initDRsearch?cid=SYC_A02）发布。用户注册、登录系统后，在数据资源栏搜索“海北站湿地”，即可下载半小时、日、月和年等时间尺度的碳水热通量及常规气象数据。也可登录Science Data Bank（http://www.sciencedb.cn/dataSet/handle/1010）访问和下载。本数据集可以用于高寒湿地碳水收支评价和相关过程模型开发、验证及对比等研究。另外，需要说明的是由于异常值剔除标准和插值方法及计算平台等因素的差异，会对碳水热通量计算结果产生较大影响[18]，本数据集的数据可能和相关研究人员的计算结果并不完全相同[7,9,12]，因此相关研究人员在数据集使用和结果认知上需要谨慎。

致 谢

感谢瓦金龙长期以来在数据采集方面的贡献。

数据作者分工职责

张法伟（1981—），男，高级工程师，研究方向为高寒草地物质循环和能量交换。主要承担工作：数据监测和论文撰写。

李红琴（1981—），女，副研究员，研究方向为高寒草地碳循环。主要承担工作：数据分析与论文修改。

赵亮（1974—），男，研究员，研究方向为全球变化和鸟类生态学研究。主要承担工作：数据整理与质量控制。

张雷明（1974—），男，副研究员，研究方向为生态系统碳水循环与全球变化。主要承担工作：数据处理和质量控制。

陈智（1981—），女，助理研究员，研究方向为生态系统碳通量时空格局。主要承担工作：数据质量分析。

祝景彬（1989—），男，博士研究生，研究方向为高寒草地碳循环。主要承担工作：数据整理和论文修改。

徐世晓（1973—），男，研究员，研究方向为草地畜牧生态学。主要承担工作：数据整理。

杨永胜（1987—），男，副研究员，研究方向为高寒草地水循环。主要承担工作：数据整理。

赵新全（1959—），男，研究员，研究方向为草地生态学。主要承担工作：数据管理和统筹。

于贵瑞（1959—），男，研究员，研究方向为生态系统与全球变化。主要承担工作：数据管理和统筹。

李英年（1962—），男，研究员，研究方向为全球变化生态学。主要承担工作：数据分析和质量控制。