CN106169014A - 基于遥感数据的区域地表感热/潜热通量反演方法及*** - Google Patents
基于遥感数据的区域地表感热/潜热通量反演方法及*** Download PDFInfo
- Publication number
- CN106169014A CN106169014A CN201610423128.8A CN201610423128A CN106169014A CN 106169014 A CN106169014 A CN 106169014A CN 201610423128 A CN201610423128 A CN 201610423128A CN 106169014 A CN106169014 A CN 106169014A
- Authority
- CN
- China
- Prior art keywords
- flux
- radiation
- region
- remotely
- inverting
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Testing Or Calibration Of Command Recording Devices (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明公开基于遥感数据的区域地表感热/潜热通量反演方法,包括以下步骤:确定研究区域和遥感数据;根据所述遥感数据制备遥感地表和区域气象参数,其中,所述遥感地表参数包括:归一化植被指数NDVI、植被覆盖度f、反照率albedo、地表比辐射率Emiss、地表温度Ts、叶面积指数LAI;所述区域气象参数包括:空气温度Ta和相对湿度RH;根据所述遥感地表和区域气象参数进行净辐射通量Rn反演;根据所述叶面积指数LAI、空气温度Ta以及净辐射通量Rn进行土壤热通量G反演;根据所述净辐射通量Rn、土壤热通量G、值被覆盖度f和地表温度Ts的理论二维空间以及感热通量估算的温度廓线方程进行辐射‑对流阻抗rae反演;根据所述辐射‑对流阻抗rae实现区域地表感热H/潜热LE通量反演。
Description
技术领域
本发明涉及区域地表感热/潜热通量的反演方法及***。更具体地,涉及基于遥感数据的区域地表感热/潜热通量反演方法及***。
背景技术
区域地表感热/潜热通量的反演是流域水资源管理、农作物灌溉制度制定、农业用水效率评价、区域旱情监测、农业节水等应用领域的基础,同时也是构建土壤-植被-大气连续体能量和水量迁移转换,进行陆-气耦合模拟研究的前提。区域地表感热/潜热通量的反演结果能够为水利、农业、林业专业技术人员提供重要的下垫面状况信息,为辅助决策提供重要支持信息。
遥感作为地面观测的重要手段,为区域地表感热/潜热通量的反演提供了重要的技术手段。现有的基于遥感的地表感热/潜热通量的方法主要分为如下几类:
统计经验模型、地表能量平衡模型、地表温度-植被指数特征空间模型、地表阻抗类模型。统计经验模型将站点观测数据同遥感数据结合,建立地表感热/潜热通量的经验回归关系,统计经验模型依赖于地表观测数据,可移植性差,难以应用于区域地表感热/潜热通量的估算。地表能量平衡模型包括单层模型和双层模型,单层模型将下垫面植被和土壤统一考虑,代表性模型如SEBAL模型和SEBS模型。
由于遥感反演的地表温度和地表空气动力学温度的差异,导致单层模型的应用和推广存在较大的不确定性。双层模型则将下垫面土壤和植被分开考虑,分别计算土壤和植被的感热/潜热通量,代表性模型如TSEB模型。双层模型参数众多,实际操作和应用受到限制。地表温度-植被指数特征空间模型基于遥感反演的地表温度和植被指数构成的二维空间来反演下垫面的感热/潜热通量信息。但特征空间的边界确定存在较大的主观性和不确定性。地表阻抗类模型是将遥感反演的下垫面参数同传统的感热/潜热模型相结合,实现区域尺度感热/潜热通量的估算。水汽阻抗的参数化是地表阻抗类模型的关键,区域尺度水汽阻抗的刻画仍然存在较多不足。
目前,基于遥感的地表感热/潜热通量模型方法仍然具有较大的不确定性,现有的模型方法中,仍然没有一种能够有效的解决无资料区、复杂非均匀下垫面的地表感热/潜热通量估算问题。尤其是无资料区、复杂非均匀下垫面的地表感热/潜热通量缺乏有效的观测手段,是目前研究和实际应用的难点。
因此,需要提供一种基于遥感数据的区域地表感热/潜热通量反演方法和***,以解决上述问题。
发明内容
本发明的目的在于提供基于遥感数据的区域地表感热/潜热通量反演方法及***,在不依赖风速和不需要地表空气动力学阻抗相关参数的模型情况下实现区域地表感热/潜热通量的反演。
为达到上述目的,本发明采用下述技术方案:
基于遥感数据的区域地表感热/潜热通量反演方法,包括以下步骤:
S1:确定研究区域和遥感数据,所述遥感数据为Landsat、ASTER或MODIS遥感数据;
S2:根据所述遥感数据或/和研究区域站点观测数据制备遥感地表参数和区域气象参数,其中,
所述遥感地表参数包括:归一化植被指数NDVI、植被覆盖度f、反照率albedo、地表比辐射率Emiss、地表温度Ts、叶面积指数LAI;
所述区域气象参数包括:空气温度Ta和相对湿度RH;
S3:根据所述遥感地表和区域气象参数进行净辐射通量Rn反演;
S4:根据所述叶面积指数LAI、空气温度Ta以及净辐射通量Rn进行土壤热通量G反演;
S5:根据所述植被覆盖度f和地表温度Ts的理论二维空间以及感热通量H估算的温度廓线方程进行辐射-对流阻抗rae反演;
所述感热通量H估算的温度廓线方程为:
其中,ρ为空气密度;Cp为空气定压比热;Taero为动力学温度;Ta为空气温度;ra为空气动力学阻抗;Ts为地表温度;rex为调和地表温度和空气动力学温度的附加阻抗;rae为辐射-对流阻抗;
基于所述净辐射通量Rn、土壤热通量G、值被覆盖度f和地表温度Ts的理论二维空间,感热通量H估算公式为:
H=(Rn-G)-LE=(Rn-G)-(1-WDI)Ep (2)
其中,LE为潜热通量;WDI为水分亏缺指数;Ep为潜在蒸散;
所述潜在蒸散Ep采用下述P-T公式:
其中,Δ为饱和水汽压曲线斜率,γ为干湿表常数;
所述水分亏缺指数WDI计算公式为:
其中,a为植被覆盖度f和地表温度Ts的理论二维空间干边同估算点的温差;b为植被覆盖度f和地表温度Ts的理论二维空间湿边同估算点的温差;计算公式分别为:
其中,
其中,Ta为空气温度;VPD为水汽亏缺,其由空气相对湿度RH计算得到;rcp为植被最小冠层阻抗;rcx为冠层最大阻抗;
联立上述公式(1)(2)(3)(4)(5)(6)可得辐射-对流阻抗rae方程:
其中,c4、c3、c2、c1以及c0为联立公式(1)(2)(3)(4)(5)(6)所得系数,所述辐射-对流阻抗rae为辐射-对流阻抗rae方程(7)中实数解。
S6:根据所述辐射-对流阻抗rae以及感热通量H估算的温度廓线方程(1) 和感热通量H估算公式(2)进行区域地表感热H/潜热LE通量反演,其中,在极端干燥情况下,所述辐射-对流阻抗rae为所述实数解中最小值。
优选的,S3中净辐射通量Rn计算公式为:
Rn=[Rsd×(1.0-albedo)+Emiss×Rld-Rlu]
其中,Rsd为短波入射辐射,albedo为反照率,Emiss为地表比辐射率,Rld为大气下行辐射,Rlu为地表上行辐射;
所述短波入射辐射Rsd、大气下行辐射Rld和地表上行辐射Rlu计算公式分别为:
Rsd=Gsc×dr×Tsw×cosθ
Rlu=δ×Emiss×Ts 4
其中,Gsc为太阳常数,大小为1367瓦/平方米;dr为日地距离修订系数;Tsw为大气透过率;θ为太阳天顶角;ea为水汽压;Ta为空气温度;δ为斯特凡波尔兹曼常数;Ts为地表温度,Emiss为地表比辐射率;
所述日地距离修订系数dr和所述大气透过率Tsw计算公式为:
其中,nDays为一年中第n天;PI等于3.1415;W为大气含水量,单位为mm;Kt为浊度系数,范围在0到1之间;θ为太阳天顶角;Pa为大气压强。
优选的,S4中,所述土壤热通量G反演过程如下:当所述叶面积指数LAI≥0.5时,所述土壤热通量G估算公式为:
G=[0.05+0.18exp(-0.521×LAI)]×Rn
当所述叶面积指数LAI<0.5时,所述土壤热通量G估算公式为:
G=1.80(Ts-273.15)+0.084×Rn
其中,Ts为地表温度。一种实现上述方法的反演***,包括:初始化模块,遥感地表参数和区域气象参数制备模块、净辐射通量Rn反演模块、土壤热通量G反演模块、辐射-对流阻抗rae反演模块以及区域地表感热H/潜热LE通量反演模块;其中,
所述初始化模块被配置为确定研究区域和遥感数据,所述遥感数据为 Landsat、ASTER或MODIS遥感数据;
所述遥感地表参数和区域气象参数制备模块被配置为根据所述遥感数据制备遥感地表参数和区域气象参数,所述遥感地表参数包括:归一化植被指数NDVI、植被覆盖度f、反照率albedo、地表比辐射率Emiss、地表温度Ts、叶面积指数LAI,所述区域气象参数包括:空气温度Ta和相对湿度RH;
所述净辐射通量Rn反演模块被配置为根据所述遥感地表和区域气象参数确定净辐射通量Rn;
所述土壤热通量G反演模块被配置为根据所述叶面积指数LAI、空气温度Ta以及净辐射通量Rn确定土壤热通量G;
所述辐射-对流阻抗rae反演模块被配置为根据方程(7)确定所述辐射-对流阻抗rae;
所述区域地表感热H/潜热LE通量反演模块被配置为根据所述辐射-对流阻抗rae以及感热通量H估算的温度廓线方程(1)和感热通量H估算公式(2)确定区域地表感热H/潜热LE通量,其中,在极端干燥情况下,所述辐射-对流阻抗rae为所述辐射-对流阻抗rae方程(7)实数解中最小值。
优选的,所述净辐射通量Rn反演模块根据短波入射辐射Rsd、反照率albedo、地表比辐射率Emiss、大气下行辐射Rld以及地表上行辐射Rlu计算净辐射通量Rn,计算公式为:
Rn=[Rsd×(1.0-albedo)+Emiss×Rld-Rlu]
所述短波入射辐射Rsd、大气下行辐射Rld和地表上行辐射Rlu计算公式分别为:
Rsd=Gsc×dr×Tsw×cosθ
Rlu=δ×Emiss×Ts 4
其中,Gsc为太阳常数,大小为1367瓦/平方米;dr为日地距离修订系数;Tsw为大气透过率;θ为太阳天顶角;ea为水汽压;Ta为空气温度;δ为斯特凡波尔兹曼常数;Ts为地表温度,Emiss为地表比辐射率;
所述日地距离修订系数dr和所述大气透过率Tsw计算公式为:
其中,nDays为一年中第n天;PI等于3.1415;W为大气含水量,单位为mm;Kt为浊度系数,范围在0到1之间;θ为太阳天顶角;Pa为大气压强。
优选的,当所述叶面积指数LAI≥0.5时,所述净辐射通量Rn反演模块中土壤热通量G估算公式为:
G=[0.05+0.18exp(-0.521×LAI)]×Rn
当所述叶面积指数LAI<0.5时,所述净辐射通量Rn反演模块中土壤热通量G估算公式为:
G=1.80(Ts-273.15)+0.084×Rn。
本发明的有益效果如下:
提供基于遥感数据的区域地表感热/潜热通量反演遥感方法及***,在不依赖风速的情况下,利用地表感热通量估算的廓线方程和地表温度和植被覆盖度所构成的理论特征空间四个角点方程,可同时实现在无资料区和复杂非均匀下垫面的地表情况下,求解地表感热通量和辐射-对流阻抗具有区域应用的有实用性和可操作性,满足区域水资源管理、农业用水效率评价、区域旱情监测、农业节水等方面的应用领域的需求。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明。
图1示出本发明所述方法流程图。
图2示出植被覆盖度f和地表温度Ts的理论二维空间。
图3示出本发明所述***结构示意图。
图4示出黑河地区地表感热H和地表潜热LE通量估算结果同观测结果的对比。
具体实施方式
为了更清楚地说明本发明,下面结合优选实施例和附图对本发明做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解,下面所具体描述的内容是说明性的而非限制性的,不应以此限制本发明的保护范围。
如图1所示,基于遥感数据的区域地表感热/潜热通量反演遥感方法包括 以下步骤:
S1:根据研究区域所需,下载Landsat、ASTER或MODIS遥感卫星产品数据;其中,Landsat和ASTER满足分辨率为100米的通量反演计算,MODIS遥感产品数据满足分辨率为1km的通量反演计算。
S2:根据所述遥感数据制备遥感地表参数和区域气象参数,所述遥感地表参数包括:归一化植被指数NDVI、植被覆盖度f、反照率albedo、地表比辐射率Emiss、地表温度Ts、叶面积指数LAI;所述区域气象参数包括:空气温度Ta和相对湿度RH。
基于Landsat卫星遥感数据的归一化植被指数NDVI、植被覆盖度f、反照率albedo、地表比辐射率Emiss、地表温度Ts以及叶面积指数LAI可参照Allen,R.G.,M.Tasumi,andR.Trezza(2007),Satellite-based energy balance for mapping evapotranspirationwith internalized calibration(METRIC)-Model,Journal of Irrigation andDrainage Engineering-Asce,133(4),380-394.的文献。
基于ASTER卫星遥感数据的归一化植被指数、植被覆盖度、反照率、比辐射率、地表温度、叶面积指数可参照French,A.N.,et al.(2005),Surface energy fluxes with theAdvanced Spaceborne Thermal Emission and Reflection radiometer(ASTER)at theIowa 2002SMACEX site(USA)(vol 99,pg 55,2005),Remote Sensing of Environment,99(4),471-471.的文献。
基于MODIS卫星遥感数据,归一化植被指数NDVI、植被覆盖度f、反照率albedo、地表比辐射率Emiss、地表温度Ts以及叶面积指数LAI等遥感参数产品可以直接计算获得。
具体的,归一化植被指数NDVI的提取,运用MODIS遥感产品数据第一波段和第二波段的反射率进行计算,计算方法为:
其中,b1和b2分别为MODIS遥感产品数据第一和第二波段的波段反射率。
植被覆盖度f为特定地域植物垂直投影面积与该地域面积之比,使用归一化植被指数NDVI计算如下:
其中,NDVImax和NDVImin分别为研究区内NDVI的最大值和最小值。
反照率albedo的遥感反演,采用Liang,S.L.,2001.Narrowband to broadbandconversions of land surface albedo I Algorithms.Remote Sensing ofEnvironment,76(2):213-238.提出的加权方法:
albedo=0.160×b1+0.291×b2+0.243×b3+0.116×b4+0.112×b5+0.018×b7-0.0015
其中,b1,b2,b3,b4,b5,b7分别为MODIS波段1,2,3,4,5,7的波段反射率。
比辐射率Emiss的反演如下:
对于裸土区,归一化植被指数NDVI小于0.2,
Emiss=0.9825-0.051×b1
对于植被和裸土混合区,归一化植被指数NDVI介于0.2和0.5之间,
Emiss=0.971-0.018×f
对于植被区,归一化植被指数NDVI大于0.5,
Emiss=0.990
对于水体,
Emiss=0.995
地表温度Ts:
Ts=DNTS×0.02
其中,DNTS为MODIS遥感产品数据中MOD11A1产品中的白天地表温度数据层。
叶面积指数LAI:
LAI=Lai1km×0.1
其中,Lai1km为MODIS遥感产品数据中MOD15A2产品中的叶面积指数数据层。
所述区域气象参数中空气温度Ta和相对湿度RH的确定,如果选择基于Landsat遥感产品数据,则需要获得研究区域内站点观测数据的支持,利用研究区域内的站点观测数据进行空间插值得到区域的空气温度Ta和相对湿度RH;如果选择基于MODIS遥感数据产品,则不需要区域内的气象观测数据支持,MODIS遥感卫星的MOD07L2产品提供卫星过境时刻的大气廓线产品,可以从MODIS大气廓线产品中获得近地层的压强Pa、空气温度Ta和相对湿度RH。
具体的,当选择基于MODIS遥感数据产品时,首先要对MOD07L2数据进行处理,采用MOD07L2产品中20个层的空气温度Ta和相对湿度RH中最近地表层的空气温度Ta和相对湿度RH参与计算,对MOD07L2的20个层从 地表层进行循环计算,直到取得有效值为止。
大气压强Pa的计算如下(单位为kPa):
Pa=0.01×SurfacePressure
其中,SurfacePressure为MOD07L2中的压强数据层;
空气温度Ta的计算如下(单位为K):
Ta=(AirTemperature+15000)×0.01
其中,AirTemperature为MOD07L2中的空气温度数据层;
露点温度Tdew的计算如下(单位为℃):
Tdew=(Mositure+15000)×0.01-273.15
其中,Mositure为MOD07L2中的露点温度数据层。
水汽压ea的计算如下(单位为hpa):
大气含水量W的计算如下(单位为mm):
W=0.01×WaterVapor
其中,WaterVapor为MOD07L2中的大气含水量数据层。
S3:根据所述遥感地表和区域气象参数进行净辐射通量Rn反演;
具体的,净辐射通量Rn计算公式为:
Rn=[Rsd×(1.0-albedo)+Emiss×Rld-Rlu]
其中,Rsd为短波入射辐射,albedo为反照率,Emiss为地表比辐射率,Rld为大气下行辐射,Rlu为地表上行辐射;
进一步,所述短波入射辐射Rsd、大气下行辐射Rld和地表上行辐射Rlu计算公式分别为:
Rsd=Gsc×dr×Tsw×cosθ
Rlu=δ×Emiss×Ts 4
其中,Gsc为太阳常数,大小为1367瓦/平方米;dr为日地距离修订系数;Tsw为大气透过率;θ为太阳天顶角;ea为水汽压;Ta为空气温度;δ为斯特凡波尔兹曼常数;Ts为地表温度,Emiss为地表比辐射率。
进一步,所述日地距离修订系数dr和所述大气透过率Tsw计算公式为:
其中,nDays为一年中第n天;PI等于3.1415;W为大气含水量,单位为mm;Kt为浊度系数,范围在0到1之间;θ为太阳天顶角;Pa为大气压强。
S4:根据所述叶面积指数LAI、空气温度Ta以及净辐射通量Rn进行土壤热通量G反演;
具体的,当所述叶面积指数LAI≥0.5时,所述土壤热通量G估算公式为:
G=[0.05+0.18exp(-0.521×LAI)]×Rn
当所述叶面积指数LAI<0.5时,所述土壤热通量G估算公式为:
G=1.80(Ts-273.15)+0.084×Rn
其中,在极端干燥情况下,所述辐射-对流阻抗rae为所述实数解中最小值。
S5:根据所述值被覆盖度f和地表温度Ts的二维空间以及感热通量H估算的温度廓线方程进行辐射-对流阻抗rae反演;
具体的,所述感热通量H估算的温度廓线方程为:
其中,ρ为空气密度;Cp为空气定压比热;Taero为动力学温度;Ta为空气温度;ra为空气动力学阻抗;Ts为地表温度;rex为调和地表温度和空气动力学温度的附加阻抗;rae为辐射-对流阻抗。
基于所述净辐射通量Rn、土壤热通量G、值被覆盖度f和地表温度Ts的二维空间,感热通量H估算公式为:
H=(Rn-G)-LE=(Rn-G)-(1-WDI)Ep (2)
其中,LE为潜热通量;WDI为水分亏缺指数;Ep为潜在蒸散。
所述潜在蒸散Ep采用下述P-T公式:
其中,Δ为饱和水汽压曲线斜率,γ为干湿表常数。
所述水分亏缺指数WDI计算公式为:
其中,如图2所示,a为植被覆盖度f和地表温度Ts二维空间干边同估算点的温差;b为植被覆盖度f和地表温度Ts二维空间湿边同估算点的温差;计算公式分别为:
其中,
联立上述公式(1)(2)(3)(4)(5)(6)可得:
其中,c4、c3、c2、c1以及c0为联立公式(1)(2)(3)(4)(5)(6)所得系数,所述辐射-对流阻抗rae为公式(7)中实数解。
S6:根据所述辐射-对流阻抗rae进行区域地表感热H/潜热LE通量反演;
具体的,根据所述公式(1)(2)估算所述区域地表感热通量H和区域地表感热潜热通量LE;其中,在极端干燥情况下,所述辐射-对流阻抗rae为所述实数解中最小值。
如图3所示,基于遥感数据的区域地表感热/潜热通量的反演***,包括:初始化模块,遥感地表参数和区域气象参数制备模块、净辐射通量Rn反演模块、土壤热通量G反演模块、辐射-对流阻抗rae反演模块以及区域地表感热H/潜热LE通量反演模块。
其中,所述初始化模块被配置为确定研究区域和遥感数据,所述遥感数据为Landsat、ASTER或MODIS遥感数据;
所述遥感地表参数和区域气象参数制备模块被配置为根据所述遥感数据或/和研究区域站点观测数据制备遥感地表参数和区域气象参数;所述遥感地 表参数包括:归一化植被指数NDVI、植被覆盖度f、反照率albedo、地表比辐射率Emiss、地表温度Ts、叶面积指数LAI,所述区域气象参数包括:空气温度Ta和相对湿度RH;
具体的,选择使用Landsat遥感数据时,需要获得研究区的站点观测数据,选择使用MODIS遥感数据时,不需要区域内的气象观测数据支持。
所述净辐射通量Rn反演模块被配置为根据所述遥感地表和区域气象参数确定净辐射通量Rn;
所述土壤热通量G反演模块被配置为根据所述叶面积指数LAI、空气温度Ta以及净辐射通量Rn确定土壤热通量G;
所述辐射-对流阻抗rae反演模块被配置为根据公式(7)确定所述辐射-对流阻抗rae;
所述区域地表感热H/潜热LE通量反演模块被配置为根据所述辐射-对流阻抗rae确定区域地表感热H/潜热LE通量。
实施例
以黑河流域中游为研究区,使用ASTER卫星遥感数据基于本发明方法反演了地表感热H/潜热LE通量,并将反演结果同研究区内的涡度相关仪器测量的感热H/潜热LE通量值进行对比。
具体的,如图4所示,地表感热H/潜热LE通量估算的均方根误差为50W/m2和60W/m2,这一结果和其他模型方法在这一区域的验证结果相近,说明本发明提出的模型方法精度较高,结果合理可信。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。
Claims (6)
1.基于遥感数据的区域地表感热/潜热通量反演方法,其特征在于,包括以下步骤:
S1:确定研究区域和遥感数据,所述遥感数据为Landsat、ASTER或MODIS遥感数据;
S2:根据所述遥感数据或/和研究区域站点观测数据制备遥感地表参数和区域气象参数,其中,
所述遥感地表参数包括:归一化植被指数NDVI、植被覆盖度f、反照率albedo、地表比辐射率Emiss、地表温度Ts、叶面积指数LAI;
所述区域气象参数包括:空气温度Ta和相对湿度RH;
S3:根据所述遥感地表和区域气象参数进行净辐射通量Rn反演;
S4:根据所述叶面积指数LAI、空气温度Ta以及净辐射通量Rn进行土壤热通量G反演;
S5:根据所述感热通量H估算的温度廓线方程和植被覆盖度f和地表温度Ts的理论二维空间进行辐射-对流阻抗rae反演;
所述感热通量H估算的温度廓线方程为:
其中,ρ为空气密度;Cp为空气定压比热;Taero为动力学温度;Ta为空气温度;ra为空气动力学阻抗;Ts为地表温度;rex为调和地表温度和空气动力学温度的附加阻抗;rae为辐射-对流阻抗;
基于植被被覆盖度f和地表温度Ts的理论二维空间,感热通量H估算公式为:
H=(Rn-G)-LE=(Rn-G)-(1-WDI)Ep (2)
其中,LE为潜热通量;WDI为水分亏缺指数;Ep为潜在蒸散;
所述潜在蒸散Ep采用下述P-T公式:
其中,Δ为饱和水汽压曲线斜率,γ为干湿表常数;
所述水分亏缺指数WDI计算公式为:
其中,a为植被覆盖度f和地表温度Ts的理论二维空间干边同估算点的温差;b为植被覆盖度f和地表温度Ts的理论二维空间湿边同估算点的温差;计算公式分别为:
其中,
其中,Ta为空气温度;VPD为水汽亏缺,其由空气相对湿度RH计算得到;rcp为植被最小冠层阻抗;rcx为冠层最大阻抗;
联立上述公式(1)(2)(3)(4)(5)(6)可得辐射-对流阻抗rae方程:
其中,c4、c3、c2、c1以及c0为联立公式(1)(2)(3)(4)(5)(6)所得系数,所述辐射-对流阻抗rae为辐射-对流阻抗rae方程(7)中实数解。
S6:根据所述辐射-对流阻抗rae以及公式(1)和(2)进行区域地表感热H/潜热LE通量反演,其中,在极端干燥情况下,所述辐射-对流阻抗rae为所述实数解中最小值。
2.根据权利要求1所述的方法,其特征在于,S3中净辐射通量Rn计算公式为:
Rn=[Rsd×(1.0-albedo)+Emiss×Rld-Rlu]
其中,Rsd为短波入射辐射,albedo为反照率,Emiss为地表比辐射率,Rld为大气下行辐射,Rlu为地表上行辐射;
所述短波入射辐射Rsd、大气下行辐射Rld和地表上行辐射Rlu计算公式分别为:
Rsd=Gsc×dr×Tsw×cosθ
Rlu=δ×Emiss×Ts 4
其中,Gsc为太阳常数,大小为1367瓦/平方米;dr为日地距离修订系数;Tsw为大气透过率;θ为太阳天顶角;ea为水汽压;Ta为空气温度;δ为斯特凡波尔兹曼常数;Ts为地表温度,Emiss为地表比辐射率;
所述日地距离修订系数dr和所述大气透过率Tsw计算公式为:
其中,nDays为一年中第n天;PI等于3.1415;W为大气含水量,单位为mm;Kt为浊度系数,范围在0到1之间;θ为太阳天顶角;Pa为大气压强。
3.根据权利要求1所述的方法,其特征在于,S4中,所述土壤热通量G反演过程如下:当所述叶面积指数LAI≥0.5时,所述土壤热通量G估算公式为:
G=[0.05+0.18×exp(-0.521×LAI)]×Rn
当所述叶面积指数LAI<0.5时,所述土壤热通量G估算公式为:
G=1.80(Ts-273.15)+0.084×Rn
其中,Ts为地表温度。
4.实现权利要求1-3中任意一项所述方法的反演***,其特征在于,包括:
初始化模块,
遥感地表参数和区域气象参数制备模块,
净辐射通量Rn反演模块,
土壤热通量G反演模块,
辐射-对流阻抗rae反演模块,以及,
区域地表感热H/潜热LE通量反演模块;
其中,所述初始化模块被配置为选择研究区域和遥感数据,所述遥感数据为Landsat、ASTER或MODIS遥感数据;所述遥感地表参数和区域气象参数制备模块被配置为根据所述遥感数据制备遥感地表参数和区域气象参数,所述遥感地表参数包括:归一化植被指数NDVI、植被覆盖度f、反照率albedo、地表比辐射率Emiss、地表温度Ts、叶面积指数LAI,所述区域气象参数包括:空气温度Ta和相对湿度RH;
所述净辐射通量Rn反演模块被配置为根据所述遥感地表和区域气象参数确定净辐射通量Rn;
所述土壤热通量G反演模块被配置为根据所述叶面积指数LAI、空气温度Ta以及净辐射通量Rn确定土壤热通量G;
所述辐射-对流阻抗rae反演模块被配置为根据方程(7)确定所述辐射-对流阻抗rae;
所述区域地表感热H/潜热LE通量反演模块被配置为根据所述辐射-对流阻抗rae以及感热通量H估算的温度廓线方程(1)和感热通量H估算公式(2)确定区域地表感热H/潜热LE通量,其中,在极端干燥情况下,所述辐射-对流阻抗rae为所述辐射-对流阻抗rae方程(7)实数解中最小值。
5.根据权利要求4所述的***,其特征在于,所述净辐射通量Rn反演模块根据短波入射辐射Rsd、反照率albedo、地表比辐射率Emiss、大气下行辐射Rld以及地表上行辐射Rlu计算净辐射通量Rn,计算公式为:
Rn=[Rsd×(1.0-albedo)+Emiss×Rld-Rlu]
所述短波入射辐射Rsd、大气下行辐射Rld和地表上行辐射Rlu计算公式分别为:
Rsd=Gsc×dr×Tsw×cosθ
Rlu=δ×Emiss×Ts 4
其中,Gsc为太阳常数,大小为1367瓦/平方米;dr为日地距离修订系数;Tsw为大气透过率;θ为太阳天顶角;ea为水汽压;Ta为空气温度;δ为斯特凡波尔兹曼常数;Ts为地表温度,Emiss为地表比辐射率;
所述日地距离修订系数dr和所述大气透过率Tsw计算公式为:
其中,nDays为一年中第n天;PI等于3.1415;W为大气含水量,单位为mm;Kt为浊度系数,范围在0到1之间;θ为太阳天顶角;Pa为大气压强。
6.根据权利要求5所述的***,其特征在于,
当所述叶面积指数LAI≥0.5时,所述净辐射通量Rn反演模块中土壤热通量G估算公式为:
G=[0.05+0.18exp(-0.521×LAI)]×Rn
当所述叶面积指数LAI<0.5时,所述净辐射通量Rn反演模块中土壤热通量G估算公式为:
G=1.80(Ts-273.15)+0.084×Rn。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610423128.8A CN106169014B (zh) | 2016-06-15 | 2016-06-15 | 基于遥感数据的区域地表感热/潜热通量反演方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610423128.8A CN106169014B (zh) | 2016-06-15 | 2016-06-15 | 基于遥感数据的区域地表感热/潜热通量反演方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106169014A true CN106169014A (zh) | 2016-11-30 |
CN106169014B CN106169014B (zh) | 2018-06-12 |
Family
ID=57359260
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610423128.8A Expired - Fee Related CN106169014B (zh) | 2016-06-15 | 2016-06-15 | 基于遥感数据的区域地表感热/潜热通量反演方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106169014B (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106771073A (zh) * | 2016-12-28 | 2017-05-31 | 中国科学院地理科学与资源研究所 | 一种基于端元信息模型估算土壤和植被蒸散发的方法 |
CN106770446A (zh) * | 2017-02-15 | 2017-05-31 | 北京佳格天地科技有限公司 | 植被蒸腾计算装置、方法及计算机*** |
CN107247690A (zh) * | 2017-06-09 | 2017-10-13 | 中国科学院寒区旱区环境与工程研究所 | 估算气温的方法及服务终端 |
CN108106676A (zh) * | 2018-02-05 | 2018-06-01 | 中国农业大学 | 一种基于遥感数据的农作物晚霜冻害的监测方法和装置 |
CN108549986A (zh) * | 2018-04-17 | 2018-09-18 | 北京师范大学 | 基于水位梯度的退化湿地生态需水量的测量方法 |
CN108984867A (zh) * | 2018-06-28 | 2018-12-11 | 中国科学院地理科学与资源研究所 | 基于gf-4和modis结合的蒸散发遥感反演方法和*** |
CN109188465A (zh) * | 2018-08-02 | 2019-01-11 | 中国科学院地理科学与资源研究所 | 基于参考像元信息的区域地表蒸散发遥感估算方法 |
CN109271605A (zh) * | 2018-10-12 | 2019-01-25 | 中国科学院地理科学与资源研究所 | 高空间分辨率遥感地表温度数据计算方法与装置 |
CN110263299A (zh) * | 2019-05-31 | 2019-09-20 | 西南大学 | 一种基于遥感的高寒草甸生态***呼吸碳排放估算方法 |
CN110991032A (zh) * | 2019-11-29 | 2020-04-10 | 中南大学 | 一种地表辐射平衡土地利用变化贡献的定量评估方法 |
CN111597692A (zh) * | 2020-04-26 | 2020-08-28 | 北京师范大学 | 一种地表净辐射估算方法、***、电子设备和存储介质 |
CN111881407A (zh) * | 2020-07-30 | 2020-11-03 | 中国科学院地理科学与资源研究所 | 一种基于遥感信息的地表水、热和碳通量耦合估算方法 |
CN111999251A (zh) * | 2020-08-14 | 2020-11-27 | 中国水利水电科学研究院 | 基于热红外遥感的区域植被蒸腾和土壤蒸发反演的遥感模型方法 |
CN112182882A (zh) * | 2020-09-28 | 2021-01-05 | 河海大学 | 一种顾及物侯信息的植被冠层蒸腾反演算法 |
CN114563353A (zh) * | 2022-01-21 | 2022-05-31 | 河海大学 | 一种基于多源卫星遥感数据土壤热通量预测方法 |
CN116737858A (zh) * | 2023-06-20 | 2023-09-12 | 长江水利委员会水文局 | 基于卫星遥感数据反演的林间穿透雨计算方法 |
CN117131313A (zh) * | 2023-10-23 | 2023-11-28 | 沈阳仪表科学研究院有限公司 | 一种中药材土壤墒情参数计算方法、计算机设备及介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101551459A (zh) * | 2008-10-15 | 2009-10-07 | 北京天宏金睛信息技术有限公司 | 基于遥感的区域蒸散量监测方法 |
CN102176002A (zh) * | 2010-12-30 | 2011-09-07 | 中国科学院地理科学与资源研究所 | 基于地表水热通量遥感反演的干旱监测方法及*** |
CN102253184A (zh) * | 2011-06-29 | 2011-11-23 | 南京信息工程大学 | 一种针对干旱半干旱地区的地表蒸散遥感反演方法 |
US20140348201A1 (en) * | 2013-05-24 | 2014-11-27 | The Regents Of The University Of California | Determining sensible heat flux without calibration |
-
2016
- 2016-06-15 CN CN201610423128.8A patent/CN106169014B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101551459A (zh) * | 2008-10-15 | 2009-10-07 | 北京天宏金睛信息技术有限公司 | 基于遥感的区域蒸散量监测方法 |
CN102176002A (zh) * | 2010-12-30 | 2011-09-07 | 中国科学院地理科学与资源研究所 | 基于地表水热通量遥感反演的干旱监测方法及*** |
CN102253184A (zh) * | 2011-06-29 | 2011-11-23 | 南京信息工程大学 | 一种针对干旱半干旱地区的地表蒸散遥感反演方法 |
US20140348201A1 (en) * | 2013-05-24 | 2014-11-27 | The Regents Of The University Of California | Determining sensible heat flux without calibration |
Non-Patent Citations (2)
Title |
---|
杨永民 等: "基于理论参数空间的遥感蒸散模型构建及验证", 《遥感技术与应用》 * |
王开存 等: "利用卫星遥感资料反演感热和潜热通量的研究综述", 《地球科学进展》 * |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106771073B (zh) * | 2016-12-28 | 2018-10-26 | 中国科学院地理科学与资源研究所 | 一种基于端元信息模型估算土壤和植被蒸散发的方法 |
CN106771073A (zh) * | 2016-12-28 | 2017-05-31 | 中国科学院地理科学与资源研究所 | 一种基于端元信息模型估算土壤和植被蒸散发的方法 |
CN106770446A (zh) * | 2017-02-15 | 2017-05-31 | 北京佳格天地科技有限公司 | 植被蒸腾计算装置、方法及计算机*** |
CN107247690A (zh) * | 2017-06-09 | 2017-10-13 | 中国科学院寒区旱区环境与工程研究所 | 估算气温的方法及服务终端 |
CN107247690B (zh) * | 2017-06-09 | 2019-10-01 | 中国科学院寒区旱区环境与工程研究所 | 估算气温的方法及服务终端 |
CN108106676B (zh) * | 2018-02-05 | 2019-06-04 | 中国农业大学 | 一种基于遥感数据的农作物晚霜冻害的监测方法和装置 |
CN108106676A (zh) * | 2018-02-05 | 2018-06-01 | 中国农业大学 | 一种基于遥感数据的农作物晚霜冻害的监测方法和装置 |
CN108549986B (zh) * | 2018-04-17 | 2020-06-23 | 北京师范大学 | 基于水位梯度的退化湿地生态需水量的测量方法 |
CN108549986A (zh) * | 2018-04-17 | 2018-09-18 | 北京师范大学 | 基于水位梯度的退化湿地生态需水量的测量方法 |
CN108984867A (zh) * | 2018-06-28 | 2018-12-11 | 中国科学院地理科学与资源研究所 | 基于gf-4和modis结合的蒸散发遥感反演方法和*** |
CN109188465A (zh) * | 2018-08-02 | 2019-01-11 | 中国科学院地理科学与资源研究所 | 基于参考像元信息的区域地表蒸散发遥感估算方法 |
CN109271605A (zh) * | 2018-10-12 | 2019-01-25 | 中国科学院地理科学与资源研究所 | 高空间分辨率遥感地表温度数据计算方法与装置 |
CN110263299A (zh) * | 2019-05-31 | 2019-09-20 | 西南大学 | 一种基于遥感的高寒草甸生态***呼吸碳排放估算方法 |
CN110991032A (zh) * | 2019-11-29 | 2020-04-10 | 中南大学 | 一种地表辐射平衡土地利用变化贡献的定量评估方法 |
CN110991032B (zh) * | 2019-11-29 | 2022-09-13 | 中南大学 | 一种地表辐射平衡土地利用变化贡献的定量评估方法 |
CN111597692A (zh) * | 2020-04-26 | 2020-08-28 | 北京师范大学 | 一种地表净辐射估算方法、***、电子设备和存储介质 |
CN111597692B (zh) * | 2020-04-26 | 2024-02-27 | 北京师范大学 | 一种地表净辐射估算方法、***、电子设备和存储介质 |
CN111881407A (zh) * | 2020-07-30 | 2020-11-03 | 中国科学院地理科学与资源研究所 | 一种基于遥感信息的地表水、热和碳通量耦合估算方法 |
CN111999251A (zh) * | 2020-08-14 | 2020-11-27 | 中国水利水电科学研究院 | 基于热红外遥感的区域植被蒸腾和土壤蒸发反演的遥感模型方法 |
CN111999251B (zh) * | 2020-08-14 | 2021-03-23 | 中国水利水电科学研究院 | 基于热红外遥感的区域植被蒸腾和土壤蒸发反演的遥感模型方法 |
CN112182882A (zh) * | 2020-09-28 | 2021-01-05 | 河海大学 | 一种顾及物侯信息的植被冠层蒸腾反演算法 |
CN112182882B (zh) * | 2020-09-28 | 2022-10-14 | 河海大学 | 一种顾及物侯信息的植被冠层蒸腾反演方法 |
CN114563353B (zh) * | 2022-01-21 | 2022-11-18 | 河海大学 | 一种基于多源卫星遥感数据土壤热通量预测方法 |
CN114563353A (zh) * | 2022-01-21 | 2022-05-31 | 河海大学 | 一种基于多源卫星遥感数据土壤热通量预测方法 |
CN116737858A (zh) * | 2023-06-20 | 2023-09-12 | 长江水利委员会水文局 | 基于卫星遥感数据反演的林间穿透雨计算方法 |
CN116737858B (zh) * | 2023-06-20 | 2024-03-22 | 长江水利委员会水文局 | 基于卫星遥感数据反演的林间穿透雨计算方法 |
CN117131313A (zh) * | 2023-10-23 | 2023-11-28 | 沈阳仪表科学研究院有限公司 | 一种中药材土壤墒情参数计算方法、计算机设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN106169014B (zh) | 2018-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106169014A (zh) | 基于遥感数据的区域地表感热/潜热通量反演方法及*** | |
Wang et al. | Comparing the diurnal and seasonal variabilities of atmospheric and surface urban heat islands based on the Beijing urban meteorological network | |
Sun et al. | Monitoring surface soil moisture status based on remotely sensed surface temperature and vegetation index information | |
CN101551459B (zh) | 基于遥感的区域蒸散量监测方法 | |
CN102176002B (zh) | 基于地表水热通量遥感反演的干旱监测方法及*** | |
Kato et al. | Estimation of storage heat flux in an urban area using ASTER data | |
Mo et al. | Simulating temporal and spatial variation of evapotranspiration over the Lushi basin | |
Sumner et al. | Utility of Penman–Monteith, Priestley–Taylor, reference evapotranspiration, and pan evaporation methods to estimate pasture evapotranspiration | |
Van den Hurk et al. | A new methodology for assimilation of initial soil moisture fields in weather prediction models using Meteosat and NOAA data | |
Wallace et al. | Soil evaporation from tiger-bush in south-west Niger | |
Jang et al. | Mapping evapotranspiration using MODIS and MM5 four-dimensional data assimilation | |
Gao et al. | A coupled remote sensing and the Surface Energy Balance with Topography Algorithm (SEBTA) to estimate actual evapotranspiration over heterogeneous terrain | |
Gao et al. | Estimating clear-sky regional surface fluxes in the Southern Great Plains Atmospheric Radiation Measurement site with ground measurements and satellite observations | |
Peng et al. | Estimation of evapotranspiration from MODIS TOA radiances in the Poyang Lake basin, China | |
Genanu et al. | Remote sensing based estimation of Evapo-transpiration using selected algorithms: the case of Wonji Shoa sugar cane estate, Ethiopia | |
Yang et al. | Daily actual evapotranspiration estimation of different land use types based on SEBAL model in the agro-pastoral ecotone of northwest China | |
Liu et al. | Actual daily evapotranspiration estimated from MERIS and AATSR data over the Chinese Loess Plateau | |
Bashir et al. | Spatial and Temporal Analysis of Evapotranspiration Using Satellite Remote Sensing Data: A Case Study in the Gezira Scheme, Sudan. | |
Li et al. | The mass and energy exchange of a Tibetan glacier: distributed modeling and climate sensitivity | |
Huo et al. | Temporal and spatial variation of surface evapotranspiration based on remote sensing in Golmud Region, China | |
Almhab et al. | Estimation of Evapotranspiration with Modified SEBAL model using landsat-TM and NOAA-AVHRR images in arid mountains area | |
Ibanez et al. | Estimation of the latent heat flux over full canopy covers from the radiative temperature | |
Xiong et al. | Estimating evapotranspiration using remote sensing in the Haihe basin | |
Ambast et al. | Estimating regional evapotranspiration using remote sensing: Application to sone low level canal system, India | |
Mullens | Evaluation and improvements of the offline CLM4 using ARM data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180612 Termination date: 20190615 |