地表蒸散发遥感产品比较与分析
李佳,辛晓洲,彭志晴,李小军

Remote Sensing Products of Terrestrial Evapotranspiration: Comparison and Outlook
Jia Li,Xiaozhou Xin,Zhiqing Peng,Xiaojun Li
表2 各产品潜热通量及显热通量计算公式对比
Table 2 The comparison of algorithms of latent heat flux and sensible heat flux
蒸散发产品算法参数化方案参考文献
MODIS-MOD16λE=λEwet_c+λEtrans+λEsoilλE为潜热通量,λEwet_cλEtransλEsoil分别表示湿冠层、干冠层及土壤的蒸散通量[36-37]
SSEBopETa=ETf×kET0

ETa为实际蒸散,ETf为蒸发比例,k为系数,ET0为草地的参考蒸散值

ETf=Th-TsTh-Tc=Th-TsdT=ρaCpRnrah(Th-Ts)

k=1.25

[35,38]
ET_PT-JPLAET=?ETs+ETc+ETi

ETc=(1-fwet)fgfTfMαΔΔ+γRnc 植被的潜热通量

ETs=fwet+fSM1-fwetαΔΔ+γ(Rn-G) 土壤的潜热通量

ETi=fwetαΔΔ+γRnc冠层截留的潜热通量

[33]
ET_ALEXIλE=λEs+λEcλEC=αPTfgss+γRn,c,?λEC为植被部分的潜热通量[39]
LSA-SAF

LEi=ζiLEi

H=?ζiHi

LEiHi分别为斑块的潜热通量和显热通量

LEi=Lvρarai+rciqsatTsk,i-qaTa

LEi=Lvρa(rai+rci)[qsatTsk,i-qa(Ta)]

[40]
GLEAME=?Ep×S+I-βI

Ep为潜在蒸散,S为土壤胁迫系数,I为降雨截留蒸散,β为截留系数取0.07

λEp=αΔΔ+ψ(Rn-G)

[34,41]
BESS

aλEj2+bλEj+c=0

λEsoil=ss+γRnsoil-Gsoil×RHD/1000

a、b、c为相应的系数,j分别表示叶片的光照和阴影部分

λEsoil表示土壤部分的潜热通量,s为饱和水汽压随气温变化的速率,γ为湿度常数,RnsoilGsoil分别表示土壤表面的净辐射和土壤热通量,RH为相对湿度,D为饱和水汽压差。

[35,42]
ETWatchLE=ΔRn-G+ρCpraes(Ta-ea)Δ+r(1+rsra)LE为潜热通量,ETWatch先由地表能量平衡模型获取晴好日蒸散发结果,将空气动力学阻抗输入至P-M模型获取连续日蒸散发[43-44]
ETMonitorET=?CpPMc+CsPMs

CcCs分别为冠层和土壤的表面阻抗系数,PMcPMs分别表示植被蒸腾和土壤水分蒸发的变量

PMc=ΔRn-G0+(ρCpVPD-ΔracRns-G0)/(raa+rac))(Δ+γ(1+rss/(raa+rac)))λ

PMs=ΔRn-G0+(ρCpVPD-ΔracRnc/(raa+ras))(Δ+γ(1+rss/(raa+ras)))λ

[45-46]