华体会体育(中国)hth·官方网站

LICA United Technology Limited

服务热线: 13910499761 010-51292601
企业邮箱
应用支持 Application Support
News 应用支持

R语言(“REddyProc”包)对涡动小时通量数据处理的方法介绍

日期: 2024-04-26
浏览次数: 251
目录
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):
R语言(“REddyProc”包)对涡动小时通量数据处理的方法介绍
●确定和滤除湍流发展较弱的时期下的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)”
R语言(“REddyProc”包)对涡动小时通量数据处理的方法介绍
图2 数据输入类型及格式
2.2
输出处理完毕数据的格式
输出的数据主要包括数据插补结果(图 3),u∗阈值估计结果(图 4)和NEE拆分为GPP和Reco的结果(图 5)。
R语言(“REddyProc”包)对涡动小时通量数据处理的方法介绍
图3 数据插补数据结果格式
R语言(“REddyProc”包)对涡动小时通量数据处理的方法介绍
图4 Ustar阈值数据结果格式
R语言(“REddyProc”包)对涡动小时通量数据处理的方法介绍
R语言(“REddyProc”包)对涡动小时通量数据处理的方法介绍
图5 NEE数据拆分结果格式
3. REddyProc包的R代码介绍
白色字为代码,“###”后仅为代码介绍的文本,无其他功能。“#”为跳过无需运行的代码。
3.1
准备—R程序包安装、运行、目标数据导入和调整
R语言(“REddyProc”包)对涡动小时通量数据处理的方法介绍
3.2
数据后处理
按照Ustar阈值计算,数据插补和NEE拆分三个流程分别进行处理。
3.2.1 Ustar阈值计算
R语言(“REddyProc”包)对涡动小时通量数据处理的方法介绍
R语言(“REddyProc”包)对涡动小时通量数据处理的方法介绍
3.2.2 数据插补
R语言(“REddyProc”包)对涡动小时通量数据处理的方法介绍
R语言(“REddyProc”包)对涡动小时通量数据处理的方法介绍
3.2.3 NEE拆分
R语言(“REddyProc”包)对涡动小时通量数据处理的方法介绍
R语言(“REddyProc”包)对涡动小时通量数据处理的方法介绍
R语言(“REddyProc”包)对涡动小时通量数据处理的方法介绍
3.2.4 整合处理结果并输出数据
R语言(“REddyProc”包)对涡动小时通量数据处理的方法介绍


News / 相关新闻 More
2025 - 08 - 20
每年,全球有超过 70亿只雄性雏鸡在孵化后被淘汰,仅因“性别不符”。这种传统做法在动物福利、资源浪费和伦理层面引发广泛关注。有没有可能——在鸡蛋孵化前,就知道性别?美国伊利诺伊大学研究团队近日在权威期刊《Food Control》发表成果,提出了一种基于高光谱成像与机器学习的非破坏性性别识别技术,成功在孵化前对鸡蛋进行精准“性别筛选”,为蛋禽行业带来颠覆性突破。实验方法简述本研究由美国伊利诺伊大学厄本那-香槟分校团队开展,实验在校内家禽科研平台进行。研究共采集了264枚白壳鸡蛋,鸡种为 White Leghorn,母鸡年龄约29周。所有鸡蛋在清洁、编号后,于12.8°C、相对湿度78%的条件下冷藏过夜,以保证实验一致性。次日,鸡蛋在室温下静置30分钟恢复温度后,使用Resonon Pika L 型高光谱成像仪进行图像采集。系统配置如下:光谱范围:374–1015nm光谱分辨率:2...
2025 - 08 - 18
研究背景:被忽视的“冬季脉冲”释放一氧化二氮(N2O)作为一种高效温室气体,其单位质量对全球变暖的影响是二氧化碳的近300倍。特别是在寒冷草原地区,每年春季的“冻融期”会爆发剧烈的N2O“热时刻”,而这些短暂却强烈的排放事件,往往被全球温室气体模型忽略。研究目标:揭秘雪下土壤N2O的爆发机制中国科学院植物研究所的研究团队,首次结合“原位高频自动通量监测”与“区域土壤柱模拟实验”,系统揭示了加深的冬季积雪如何显著放大草地土壤N2O排放,并进一步明确了水分与微生物联动机制在这一过程中的核心作用。图1. 研究区概况实验一:原位自动监测,捕捉全年N2O变化趋势地点:内蒙古草原生态系统研究站(IMGERS)方法:在天然草地中布设雪围栏制造“深雪处理”,并使用SF-3500系列多通道土壤气体通量测量系统(华体会体育(中国)hth·官方网站),配合高精度激光光腔分析仪,实现全年不间断、每日高频率N2O通量监测。亮...
2025 - 08 - 11
在中国海沿岸的红树林中,有一种名叫桐花树 (Aegiceras corniculatum) 的植物,静静地记录着环境的蛛丝马迹。近日,广东海洋大学与中科院广州地化所、加州理工学院等单位合作,结合叶蜡正构烷烃和脂肪酸的δ2H分析及环境参数监测,首次构建了盐度-代谢调控-同位素分馏的关联模型,不仅修正了红树林同位素分馏理论模型,还为重建热带沿海古环境提供了全新视角。图1.研究区概况· 研究地点:湛江红树林国家级自然保护区植物叶片中的蜡质化合物,能“记住”它们吸收水分的氢同位素特征(δ²H),而这一特征受降水、盐度等因素影响。团队选取了湛江高桥、九洲与营仔三条河口,分别在旱季(2021年11月)与雨季(2022年6月),采集了桐花树的叶片、木质部、叶水、水体和沉积物孔隙水等样品,系统分析了不同季节和盐度梯度下的氢氧同位素变化。在本研究中,科研人员采用LI-2100 全自动真空...
2025 - 08 - 11
在应对气候变化的全球行动中,“土壤碳汇”这一关键词越来越频繁地出现在科学家们的研究中。在全球气候治理日益紧迫的今天,如何提升土壤碳汇能力,成为科学界关注的热点。近日,一项发表在国际期刊《Land Degradation & Development》的最新研究,聚焦非洲萨赫勒地区,通过实地采样与建模分析,揭示了不同放牧强度对草原土壤有机碳储量的影响机制。 图1. 研究区域和实验地点的位置。 研究背景 非洲萨赫勒地区是全球典型的半干旱草原生态系统,牧业在当地生计中占据重要地位。然而,关于放牧对土壤碳储量的具体影响,尤其是在不同强度下的机制,一直缺乏系统的实地数据。研究方法§采用多点样方布设,覆盖不同放牧压力区域§采集0–30cm深度的土壤进行碳氮分析§结合ASDLabSpec4地物光谱仪进行快速土壤属性预测§使用结构方程模型(SEM)解析变量间直...
关闭窗口】【打印
Copyright ©2018-2023 华体会体育(中国)hth·官方网站
犀牛云提供企业云服务

华体会体育(中国)hth·官方网站

地址:北京市海淀区安宁庄东路18号光华创业园5号楼(生产研发)
          光华创业园科研楼四层
电话:13910499761 13910124070  010-51292601
传真:010-82899770-8014
邮箱:info@li-ca.com
邮编:100085

 



 


 


  • 您的姓名:
  • *
  • 公司名称:
  • *
  • 地址:
  • *
  • 电话:
  • *
  • 传真:
  • *
  • 电子邮箱:
  • *
  • 邮政编码:
  • *
  • 留言主题:
  • *
  • 详细说明:
  • *
在线留言
关注我们
  • 官方微信
  • 官方手机端
友情链接:
X
1

QQ设置

3

SKYPE 设置

4

阿里旺旺设置

等待加载动态数据...

等待加载动态数据...

5

电话号码管理

  • 010-51292601
6

二维码管理

等待加载动态数据...

等待加载动态数据...

展开