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