Penman equation explained

The Penman equation describes evaporation (E) from an open water surface, and was developed by Howard Penman in 1948. Penman's equation requires daily mean temperature, wind speed, air pressure, and solar radiation to predict E. Simpler Hydrometeorological equations continue to be used where obtaining such data is impractical, to give comparable results within specific contexts, e.g. humid vs arid climates.

Details

Numerous variations of the Penman equation are used to estimate evaporation from water, and land. Specifically the Penman–Monteith equation refines weather based potential evapotranspiration (PET) estimates of vegetated land areas.[1] It is widely regarded as one of the most accurate models, in terms of estimates.

The original equation was developed by Howard Penman at the Rothamsted Experimental Station, Harpenden, UK.

The equation for evaporation given by Penman is:

Emass=

mRn+\rhoacp\left(\deltae\right)ga
λv\left(m+\gamma\right)

where:

m = Slope of the saturation vapor pressure curve (Pa K−1)

Rn = Net irradiance (W m−2)

ρa = density of air (kg m−3)

cp = heat capacity of air (J kg−1 K−1)

δe = vapor pressure deficit (Pa)

ga = momentum surface aerodynamic conductance (m s−1)

λv = latent heat of vaporization (J kg−1)

γ = psychrometric constant (Pa K−1)

which (if the SI units in parentheses are used) will give the evaporation Emass in units of kg/(m2·s), kilograms of water evaporated every second for each square meter of area.

Remove λ to obviate that this is fundamentally an energy balance. Replace λv with L to get familiar precipitation units ETvol, where Lv=λvρwater. This has units of m/s, or more commonly mm/day, because it is flux m3/s per m2=m/s.

This equation assumes a daily time step so that net heat exchange with the ground is insignificant, and a unit area surrounded by similar open water or vegetation so that net heat & vapor exchange with the surrounding area cancels out. Some times people replace Rn with and A for total net available energy when a situation warrants account of additional heat fluxes.

Temperature, wind speed, relative humidity impact the values of m, g, cp, ρ, and δe.

Shuttleworth (1993)

In 1993, W.Jim Shuttleworth modified and adapted the Penman equation to use SI, which made calculating evaporation simpler.[2] The resultant equation is:

Emass=

mRn+\gamma*6.43\left(1+0.536*U2\right)\deltae
λv\left(m+\gamma\right)

where:

Emass = Evaporation rate (mm day−1)

m = Slope of the saturation vapor pressure curve (kPa K−1)

Rn = Net irradiance (MJ m−2 day−1)

γ = psychrometric constant =

0.0016286*PkPa
λv
(kPa K−1)

U2 = wind speed (m s−1)

δe = vapor pressure deficit (kPa)

λv = latent heat of vaporization (MJ kg−1)

Note: this formula implicitly includes the division of the numerator by the density of water (1000 kg m−3) to obtain evaporation in units of mm d−1

Some useful relationships

δe = (es - ea) = (1 – relative humidity) es

es = saturated vapor pressure of air, as is found inside plant stoma.

ea = vapor pressure of free flowing air.

es, mmHg = exp(21.07-5336/Ta), approximation by Merva, 1975[3] Therefore

m=\Delta=

des
dTa

=

5336
2
T
a
\left(21.07-
5336
Ta
\right)
e
, mmHg/K

Ta = air temperature in kelvins

See also

References

Notes and References

  1. Book: Allen, R.G. . Pereira, L.S. . Raes, D. . Smith, M. . Crop Evapotranspiration—Guidelines for Computing Crop Water Requirements . 2007-10-08 . FAO Irrigation and drainage paper 56 . 1998 . Food and Agriculture Organization of the United Nations . Rome, Italy . 92-5-104219-5 .
  2. Shuttleworth, J., Putting the vap' into evaporation http://www.hydrol-earth-syst-sci.net/11/210/2007/hess-11-210-2007.pdf
  3. Merva, G.E. 1975. Physio-engineering Principles. AVI Publishing Company, Westport, CT.