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.
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. 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:
- 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)
- ga = momentum surface aerodynamic conductance (m s−1)
- δe = vapor pressure deficit (Pa)
- λ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/(m²·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.
- 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 = (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
Therefore , mmHg/K
- Ta = air temperature in kelvins
- Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. (1998). Crop Evapotranspiration—Guidelines for Computing Crop Water Requirements. FAO Irrigation and drainage paper 56. Rome, Italy: Food and Agriculture Organization of the United Nations. ISBN 92-5-104219-5. Retrieved 2007-10-08.
- Shuttleworth, J., Putting the vap' into evaporation http://www.hydrol-earth-syst-sci.net/11/210/2007/hess-11-210-2007.pdf
- Merva, G.E. 1975. Physio-engineering Principles. AVI Publishing Company, Westport, CT.
- Jarvis, P.G. (1976) The interpretation of the variations in leaf water potential and stomatal conductance found in canopies in the field. Phil. Trans. R. Soc. Lond. B. 273, 593-610.
- Neitsch, S.L.; J.G. Arnold; J.R. Kliniry; J.R. Wolliams. 2005. Soil and Water Assessment Tool Theoretical Document; Version 2005. Grassland, Soil and Water Research Laboratory; Agricultural Research Service. and Blackland Research Center; Texas Agricultural Experiment Station. Temple, Texas. http://www.brc.tamus.edu/swat/downloads/doc/swat2005/SWAT%202005%20theory%20final.pdf
- Penman, H.L. (1948): Natural evaporation from open water, bare soil and grass. Proc. Roy. Soc. London A(194), S. 120-145.