关于潜热的理解Latent heat flux

以前只知道可以根据能量平衡的方法计算蒸发蒸腾量,“余项法”就是把潜热做为剩余项计算出。也知道,涡度相关计或者通量塔也可以测量潜热,由潜热就可以计算蒸发蒸腾量。那么潜热的单位是什么,和ET之间如何转换?

潜热的单位一般为W/m2,表示单位面积上的功率,它与时间没有关系,一般每1小时或半小时测量一次,它每天的值就是每天测量n次的平均值。

  • Energy flux–Latent heat flux, Hl (W/m2)

蒸发就比较容易理解,

  • Water flux–Evaporation rate, E (mm/day)

它们之间的关系为:

Hl= ρlv E

ρ= 1000 kg/m3

lv = 2.5 MJ/kg

那么, 1mm/day应该相当于多少W/m2

1000(kg/m3)*2.5*106(J/kg)*1mm/day*(1/86400)(day/s)*(1/1000)(mm/m)

=28.94W/m2

当然也会受温度的一些影响:

Temp Lv Density Conversion
0 2501000 999.9 28.94
10 2477300 999.7 28.66
20 2453600 998.2 28.35
30 2429900 995.7 28.00
40 2406200 992.2 27.63
  • Energy Balance of Earth

http://www.uwsp.edu/geo/faculty/ritter/geog101/textbook/energy/radiation_balance.html

energy-balance-of-earth

可以编写下列程序将测量的LE转换为每天的ET值:

%%
% time range
datedaily=[datenum('2010-11-15'):datenum('2009-11-9')]';

% get rid of the lower accuracy data
LE(LE<-1000)=NaN;
Tem(Tem<-1000)=NaN;

% daily data, because meausure 48 times each day;
LE_daily=zeros(372,1);
Tem_daily=zeros(372,1);

for i=1:372
 LE_daily(i)=nanmean(LE((i*48-47:i*48),1));
 Tem_daily(i)=nanmean(Tem((i*48-47:i*48),1));
end
% interpolate conversion parameter 
% http://www.uwsp.edu/geo/faculty/ritter/geog101/textbook/energy/radiation_balance.html
Temp=[0;10;20;30;40]
Conversion=[28.9400000000000;28.6600000000000;28.3500000000000;28;27.6300000000000]
Conversion1=interp1(Temp,Conversion,Tem_daily,'linear')
ET_daily=LE_daily./Conversion1;

Leave a Reply

Your email address will not be published. Required fields are marked *