不确定性量化FAO-56 参考蒸散发模型 · FAO-56 Reference ET Model — Uncertainty Quantification

01

模型定义

标准化的参考作物（草/苜蓿）蒸散发计算模型，用于灌溉调度与作物系数法。

02

公式

ET0 由 Penman-Monteith 在固定参考面（rs=70 s/m, h=0.12 m）参数下计算。

公式中各量的定义与单位见"输入 / 输出"。

03

输入与输出

项	内容
输入	气温、湿度、风速、太阳辐射（或日照时数）、海拔；量化随机/系统误差支持年累计与决策。
输出	参考蒸散发 ET0（mm/d）

04

参考依据

Allen et al. (1998) FAO-56

05

成熟度与边界

成熟度

Established

领域

通量与涡度 · Flux & Eddy Covariance

适用边界ET0×Kc 得到作物 ET；通量塔实测 ET 可验证 Kc。传播误差；报告随机与系统分量。

06

机理与过程链

通量观测以涡度相关法为主：三维超声风速仪与气体分析仪以 10 Hz 采样，通过协方差计算 CO₂/H₂O/热量通量，需做坐标旋转、频率响应修正、WPL 密度修正与夜间 u* 过滤。

通量的代表性由 footprint 决定：90% 贡献区通常在上风向 50~200 m（视高度与大气稳定度），footprint 越界（下垫面不均一、地形起伏）时通量不能代表目标生态系统。能量闭合率 EBR=(LE+H)/(Rn−G) 常为 0.7~0.9，缺口来自储存项、平流与测量偏差。

07

核心假设与成立条件

假设不成立时，模型偏差通常最先体现在趋势而非绝对值上。

假设	工程含义与成立条件
稳态与水平均一	30 min 内湍流统计平稳且下垫面均一
充分混合	u* 超过阈值（森林 0.1~0.2 m/s、草地 0.05 m/s），夜间弱湍流需过滤
谱相似	惯性子区谱斜率符合 −5/3，频率响应损失经修正 <10%
传感器线性	气体分析仪在量程内线性响应，标气定期校准

08

误差来源与失效模式

出现下列表现时应先复核假设与输入，而不是调整输出。

失效模式	表现与判据
弱湍流夜间	u* 低于阈值时通量系统性低估（夜间泄流），需 u* 过滤与补缺
footprint 越界	风向偏转至道路/水体/建筑，通量混入异质源
降雨污染	雨滴与露水污染光路，通量出现尖峰或负值
传感器漂移	气体分析仪零点漂移使月尺度通量偏差 5%~15%

09

输入数据条件与精度要求

本模型最小可用数据集：气温、湿度、风速、太阳辐射（或日照时数）、海拔

数据项	采样与精度要求
三维风速	10 Hz，超声风速仪，路径长度 ≤15 cm，每半年检查零点
气体浓度	10 Hz，CO₂/H₂O 闭路或开路，每月标气校准
辐射与土壤热通量	1 min/30 min，用于能量闭合核对
湍流统计	30 min 窗口，输出 u*、稳定性、footprint 与 QC 标志

输入

气温、湿度、风速、太阳辐射（或日照时数）、海拔；量化随机/系统误差支持年累计与决策。

输出

参考蒸散发 ET0（mm/d）

10

部署、标定与运行要点

围绕不确定性量化FAO-56 参考蒸散发模型的运行前检查清单。

环节	要点
采样与计算	10 Hz 原始数据全存，30 min 通量计算，每段输出 0~3 级 QC 标志
标定	每月标气校准（零气/跨度气），每季检查光路清洁
夜间过滤	u* 阈值按站点湍流统计确定，过滤率 20%~40% 属正常
能量核对	与辐射、土壤热通量同步观测，月度核对闭合率

11

验证指标与参考阈值

阈值为工程常用参考区间，具体项目按合同与规范取值。

指标	参考阈值
能量闭合率	0.7~0.9（森林常 0.7~0.8）
有效数据覆盖率	>60%（年尺度）
谱检验	惯性子区斜率 −5/3 ± 0.2
月尺度 NEE 不确定度	±10%~15%（与生物量/夜间池法互校）

12

常见问答

夜间通量为什么常为负或缺失？
弱湍流下涡度相关低估夜间呼吸，需 u* 过滤并用夜间池法或模型补缺；缺失率 30%~50% 属正常。

能量不闭合说明数据不能用吗？
闭合率 0.7~0.9 是涡度相关的已知局限（储存项与平流），不代表数据失效；长期趋势与年总量仍可用，须注明闭合率。

footprint 怎么判断？
按风向、稳定度与高度计算 90% 贡献区，越出均一下垫面即标记为异质源时段，不参与生态系统尺度统计。

14

与相邻方法的取舍

同一目标量常有多种建模路线，取舍取决于数据条件与精度要求。

方法	适用	局限
涡度相关法	直接、连续、生态系统尺度	夜间弱湍流失效，能量不闭合
廓线/梯度法	适合弱湍流与分层	需多高度观测与相似理论假设
静态箱法	单点精度高、可分源	空间代表性差，难以连续