R语言（“REddyProc”包）对涡动小时通量数据处理的方法介绍3387
发表时间：2020-06-04 14:41 目录 1. 后处理方法介绍 1.1 Ustar阈值判断（主要针对夜间NEE） 1.2 数据插补 1.2.1 查表法插补（LUT法） 1.2.2 平均日变化曲线法（MDC法） 1.2.3 样本边缘分布采样法（MDS法） 1.3 数据拆分 2. REddyProc包处理数据格式介绍 2.1 输入需要处理数据的格式 2.2 输出处理完毕数据的格式 3. REddyProc包的R代码介绍 3.1 准备—R程序包安装、运行、目标数据导入和调整 3.2 数据后处理 3.2.1 Ustar阈值计算 3.2.2 数据插补 3.2.3 NEE拆分插补 3.2.4整合处理结果并输出数据 涡动通量数据处理分为在线处理（online-processing）和后处理（post-processing）。其中在线处理针对高频通量数据（e.g.10Hz data）通过一系列标准方法进行计算，最后得到带有质量评价的低频通量数据（e.g.half-hour data），后处理主要包括Ustar阈值估计、数据插补和碳通量（NEE）拆分（植被总生产力GPP和呼吸消耗Re）及其结果的可视化表达。 当夜间大气湍流运动较弱时，摩擦风速u∗降低，涡动相关系统测量碳通量NEE时会出现低估的现象，数据漂移值增多。通常需要判断出u∗阈值，剔除这些低于u∗阈值的NEE；对缺失的数据进行插补，有利于得到完整的时间序列并得到更长时间尺度（月或年）下的均值；NEE通过主流的模型方法进行拆分，以便进一步了解研究区NEE两大组分：（1）生态系统总生产力（或总初级生产力）（2）生态系统呼吸。REddyProc 程序包通过R语言平台实现了以上三个方面的数据后处理，以及对其计算结果实现基本可视化功能。 1. 后处理方法介绍 数据后处理所使用的通量数据是已经过异常值剔除后的数据，NEE拆分或可插补的数据包括碳通量（NEE，umolm-2s-1）、感热通量（sensible heat flux (H) Wm-2）、潜热通量（latent heat flux (LE) Wm-2）、摩擦风速（friction velocity (u∗) ms-1）、入射短波辐射（global radiation (Rg) Wm-2）, 空气或土壤温度（air or soil temperature (Tair, Tsoil) ℃）和水汽压饱和差（vapor pressure deficit (VPD) hPa）或相对湿度（relative humidity (RH) %）。其中u∗、Rg、VPD、Tair和RH是NEE滤除、插补和拆分默认使用数据。 数据后处理主要流程包括（图 1）：  ●确定和滤除湍流发展较弱的时期下的NEE（计算u∗阈值）。 ●插补缺失的小时数据。 ●拆分碳通量小时数据,得到GPP和Reco。 图1 数据后处理流程，以某一站点数据为例 （Wutzler et al.(2018)） 1.1 Ustar阈值判断（主要针对夜间NEE） 仪器所在高度处可以测量到下垫面全部碳通量（无平流损失），对应的最小u∗称为u∗阈值，u∗阈值通常出现在夜间（Rg <10 Wm-2）。由于下垫面粗糙度在不同时期（季节）发生变化，导致u∗阈值会产生季节变化。（the u∗ threshold is the minimum u∗ above which respiration reaches aplateau. This threshold is specific for each season of a site year.）。 当前REddyProc 包计算Ustar阈值方法主要有移动点法（the moving point method，MPT）和断点检测法（the breakpoint detection method，CPT），其中MPT较常用。 1.2 数据插补 Ustar阈值滤除NEE后，会有更多的NEE缺失数据，需要插补。 1.2.1 查表法插补（LUT法） 在REddyProc包的查表法中（look-up table (LUT)），所有通量数据以特定的时间窗口内的相似气象条件为依据进行分类并计算平均值，最后得到可供参照的速查表。缺失的数据可利用同时间序列中已知的气象数据与速查表匹配，对应的通量数据即为所缺失的数据。 1.2.2 平均日变化曲线法（MDC法） 该方法可在其他气象数据缺失条件下进行通量数据插补。假设植物晚上只进行呼吸作用，白天发生光合和呼吸作用，且NEE具有较为规律的日变化特征。则缺失的数据可根据临近天同时刻（或前后一小时）已知的通量数据进行插补（mean diurnal course (MDC)）。 1.2.3 样本边缘分布采样法（MDS法） 边缘分布采样法（marginal distribution sampling (MDS)）结合了以上LUT和MDC两种方法，根据通量数据与气象因子之间的关系（covariation）以及通量数据在时间上的自相关进行插补。MDS可针对较大缺失范围的NEE和LE数据插补，该方法目前最受欢迎。 利用Rg, Tair和 VPD三种气象数据，（1）如果三个气象数据皆未缺失, 使用LUT 方法，三个气象因子默认边际条件（default margins）为50 Wm−2, 2.5 ◦C和5.0 hPa；（2）Tair 或VPD 缺失, 则只利用 Rg；(3) 如果三种气象数据都缺失，使用 MDC方法。另外，很多站点没有Rg的观测数据，可用光合有效辐射par代替，并设置par的边际条件（可尝试使用100-200 μmol m-2 s-1） 1.3 数据拆分 NEE、Reco（↑）和GPP（↓）三者关系为NEE = Reco– GPP，当前NEE拆分为Reco 和GPP主要方法有利用夜间NEE数据拆分和利用白天NEE数据拆分两种。当前夜间NEE数据拆分方法最常用。 夜间NEE数据拆分方法是假设植被呼吸Reco只与Tair变化有关，且夜间植被只进行呼吸作用，因此可以通过夜间NEE对Tair的响应变化曲线推出白天植被的呼吸Reco变化，最后根据以上关系式求出植被总生产力GPP。 白天NEE数据拆分方法是将白天NEE和总辐射的关系假设为Rg和VPD对GPP的影响以及Tair对Reco的影响的综合。 2. REddyProc包处理数据格式介绍 本节图片来源： https://www.bgc-jena.mpg.de/bgi/index.php/Services/REddyProcWebDataFormat 注意虽然REddyProc包是基于该网页在线工具所开发的，但是二者的算法还有一些区别，详情参见Wutzler et al.(2018)。 2.1 输入需要处理数据的格式 输入数据格式如图2所示，输入文件类型为“文本文件（制表符分隔）（*.txt）”  图2 数据输入类型及格式 2.2 输出处理完毕数据的格式 输出的数据主要包括数据插补结果（图 3），u∗阈值估计结果（图 4）和NEE拆分为GPP和Reco的结果（图 5）。  图3 数据插补数据结果格式  图4 Ustar阈值数据结果格式   图5 NEE数据拆分结果格式 3. REddyProc包的R代码介绍 3.1 准备—R程序包安装、运行、目标数据导入和调整
3.2 数据后处理 按照Ustar阈值计算，数据插补和NEE拆分三个流程分别进行处理。 3.2.1 Ustar阈值计算   3.2.2 数据插补   3.2.3 NEE拆分    3.2.4 整合处理结果并输出数据
使用过程中，如您遇到相关问题，敬请与我们联系 电话：010-51290427、13910124057 网址：www.servirst.com 邮箱：quan.yang@servirst.com ning.zheng@servirst.com 或扫描下方公众号二维码，进行留言。 微信号 : 北京萨维福特科技有限公司长按识别二维码关注  |
