(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111491921.9 (22)申请日 2021.12.08 (71)申请人 中国科学院地理科 学与资源研究所 地址 100101 北京市朝阳区大屯路甲1 1号 申请人 中国标准 化研究院 (72)发明人 袁烨城　刘伟　李宝林　王双　 高锡章　蒋育昊　李影　 (74)专利代理 机构 武汉领君知识产权代理事务 所(普通合伙) 42248 代理人 汪俊锋 (51)Int.Cl. G06F 30/20(2020.01) G06F 119/08(2020.01) (54)发明名称 一种基于遥感的陆表蒸散 计算方法 (57)摘要 本发明属于地理学与生态学领域， 涉及一种 基于遥感的陆表蒸散计算方法。 陆表蒸散由植被 蒸腾、 土壤蒸发与拦截蒸发组成。 本发明主要对 植被蒸腾中的温度和土壤蒸发中的土壤水分进 行了优化， 主要包括两部分： 使用物理模型计算 叶片温度代替原方法中环境温度进行植被蒸腾 计算； 使用不同土壤层的土壤含 水量来重新定义 土壤蒸发过程中的土壤水分限制。 本发明的方法 简单实用， 所涉及的数据均可从遥感平台、 中国 气象数据网、 涡旋通量塔以及相关文献免费获 取。 本发明属于机理模型， 外推时无需参数校正， 可以更低的成本快速地对陆表蒸散进行可靠的 估算。 权利要求书2页 说明书13页 附图1页 CN 114417553 A 2022.04.29 CN 114417553 A 1.一种基于遥感的陆表蒸散计算方法， 其特征在于， 实际蒸散(E)与冠层 表面蒸腾(Ec)、 土壤表面蒸发(Es)及拦截蒸发(Ei)有关， 具体 计算方法如下： LE＝L(Ec+Es+Ei) 其中， L表示蒸发潜热(MJ/k g)； E为蒸散(mm)； LE为蒸发潜热与蒸散的乘积， 表示潜热通 量(W/m2)； Ec为冠层表面蒸腾(mm)； Es为土壤表面蒸发(mm)； Ei表示拦截蒸发(mm)； LEc， LEs与 LEi分别表示各部分的潜热通量； αPT为无量纲系数， 一般取值1.26； Δl， Δa与Δs分别表示饱 和水汽压对叶片温度、 气温及土壤温度的曲线斜率(kPa/℃)； γ是干湿表常数(kPa/℃)； Rnc 为冠层表面的净辐射(W/m2)； Rns为土壤表面净辐射(W/m2)； G为土壤热通 量(W/m2)； 由以上公式， 可 得潜在蒸散(PE)与实际蒸散(E)分别为： PE在kc、 ki与ks下转为E， kc、 ki与ks分别表示植被潜在蒸腾到实际蒸腾、 植被表面潜拦截 蒸发到实际拦截蒸发、 土壤拦截潜在蒸发到实际蒸发的限制因子， 计算方法分别为： kc＝(1‑fwet)fgfTfM ks＝fwet+(1‑fwet)fsm ki＝fwet 其中， fwet表示拦截表面占比； fgfTfM分别表示植物表面受到的绿度、 温度以及 水分的胁 迫系数； fsm为土壤表面的水分胁迫系数； 以上五个指标的取值范围都是0 ‑1， 当数值为1时， 表示不存在环境胁迫； 数值越小， 表示环境胁迫越大； 各因子的计算和改进方法如下： (1)温度胁迫系数(fT) 其中， Tl表示叶片温度； Topt表示最适温度； fT的数值范围为0 ‑1。 (2)拦截表面占比(fwet) 其中， Rh表示相对湿度； (3)叶温(Tl) 叶温(Tl)表达式为：权　利　要　求　书 1/2 页 2 CN 114417553 A 2其中， Ta为空气温度； Rnc为落在植被冠层的净辐射； cP为空气定压比热； δb为叶片边界导 度， 其中， μ代 表风速， d与叶片宽度有关； (4)植被绿度(fg) 植被绿度表示 植被蒸腾中， 被绿色植被吸 收的部分， 植被绿度的计算方法为： 其中， fg表示植被绿度； fAPAR表示活体植被吸收的陆表太阳辐射占比； fIPAR表示由植被 拦截的所有太阳辐射占比； (5)植被水分胁迫系数(fM) 植被水分限制的表达式为： 其中， 表示该像元绿色植被吸 收的太阳辐射占比的最大值； (6)土壤水分胁迫系数(fsm) 其中， SM表示土壤水分含量， SMmax与SMmin分别表示土壤水分含量的最大值与最小值； VPD表示水汽压偏差； (7)干湿表常数(γ) 其中， Cp为空气定压比热， 表示在一定气压 下， 单位体积的空气温度每升高1℃所需的能 量； P为大气压， 与海拔有关； λ为蒸发潜热； ε为水汽分子量与干空气分子量之比； (8)太阳辐射： 植被表面(Rnc)和土壤表面(Rns) Rnc＝FVC×Rn Rns＝(1‑FVC)×Rn‑G 其中， FVC表示 植被覆盖度， Rn为植被净辐射， G表示土壤热通 量。权　利　要　求　书 2/2 页 3 CN 114417553 A 3