Hydrus结果的T_Level.out文件中,在没有降水和灌溉的情况下,地表的潜在通量+根系的潜在吸水速率(rTop+rRoot),应该就是ET0,但是他们之间差了0.24倍。这个0.24是从哪里出来的呢? 查阅技术手册可知p37-38,潜在通量和根系的潜在吸水速率计算的都是瞬时值。 Variations in potential evaporation and transpiration during the day can be generated with HYDRUS-1D using the assumptions that hourly values between 0-6 a.m. and 18-24 p.m. represent 1% of the total daily value and that a sinusoidal shape is followed during the rest of the day [Fayer, 2000], i.e.,
假设上午0-6点到下午18-24点之间的小时值代表全天总日值的1%,并且在一天的其余时间遵循正弦形状[Fayer, 2000],可以用HYDRUS-1D生成白天潜在蒸发和蒸腾的变化。
where Tp- is the daily value of potential transpiration (or evaporation).
Matlab 程序如下:
% 日值转化小时值
% 假如日蒸发蒸腾量为5mm。
V_daily=5/24; % 每天每小时的平均值,mm
V_hourly=zeros(1,24); % 预设每小时的值,mm/h
for i=1:24;
if i>0.736*24 || i<0.264*24;
V_hourly(i)=0.24*V_daily;
else
V_hourly(i)=2.75*V_daily.*sin(2*pi*(i-1)/24-pi/2);
end
end
做图如下:
每时刻的值如表所示:
时刻 ET/mm
0.05mm/小时5mm/day的1%。0.05*24=1.2mm/days是0时刻的值,然而每天的ET为5mm。1.2/5=0.24,这个0.24就是这么产生的,因此,应该用累加的值相减来计算日通量更方便一些(diff)。