Abstract. A new method has been proposed for estimating sensible heat flux from singlelevel measurement of air temperature. When turbulent transfer of heat in the lower atmosphere over a homogeneous surface is modeled by a one-dimensional diffusion equation with a constant diffusivity, heat flux can be expressed as a weighted average (half-order derivative) of the time history of air temperature. This formula provides an approximate solution of the diffusion equation where the (eddy) diffusivity characterizing the turbulent flow is not constant. Eddy diffusivity has been formulated based on MoninObukhov similarity theory with Businger-Dyer stability functions or determined by an empirical equation. The knowledge of surface parameters including friction velocity sometimes is needed to apply this method. The method was tested against observations collected during two field experiments, FIFE and ABRACOS. The close agreement between the estimated sensible heat flux and observations suggests that this novel approach is a potentially powerful tool in evaluating the energy balance at the land surface.
IntroductionThe fundamentals of the commonly used approaches for estimating sensible heat flux have remained unchanged for a long time. Direct measurement of sensible heat flux is possible with a fast-response turbulence instrument that can sense the rapid fluctuations of temperature and velocity. Application of this method, known as the eddy correlation method, is limited by its stringent instrumental requirements. Hence many indirect methods have been developed to estimate sensible heat flux using observations of meteorological variables measured in situ or remotely. They include the Bowen ratio method, the bulk transfer method, the gradient method, the profile method, thermodynamic energy equation We assume the following initial profile of T in deriving the desired solution rl,=0 = z > 0.This condition implies that the model starts when the sensible heat flux is zero. It turns out that this specific condition does not affect the generality of the solution.The full solution of T as a function of z and t requires an additional boundary condition at the surface. The surface is understood as the actual height above which (1) is valid. An-2281