以前只知道可以根据能量平衡的方法计算蒸发蒸腾量,“余项法”就是把潜热做为剩余项计算出。也知道,涡度相关计或者通量塔也可以测量潜热,由潜热就可以计算蒸发蒸腾量。那么潜热的单位是什么,和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
可以编写下列程序将测量的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;
