Rusteau-1D Help
Theory and equations

Penman–Monteith combination equation

With meteorological data, the potential evapotranspiration is calculated with the Penman–Monteith combination equation Monteith (1981); Monteith and Unsworth (1990) in the form of FAO Irrigation and Drainage Paper 56 Allen et al. (1998):

ETp = (1/λ) [ Δ (Rn − G) + ρa cp (es − ea) / ra ] / [ Δ + γ (1 + rc/ra) ]
SymbolMeaning
ETppotential evapotranspiration [mm/d], converted to the project's units
λlatent heat of vaporization, λ = 2.501 − 0.002361 T [MJ/kg]
Δslope of the saturation vapor pressure curve [kPa/°C]
Rnnet radiation at the surface [MJ/m²/d]
Gsoil heat flux [MJ/m²/d], neglected for daily values
ρa cpvolumetric heat capacity of the air
es − eavapor pressure deficit [kPa]
ra, rcaerodynamic and canopy resistance [s/m]
γpsychrometric constant [kPa/°C], γ = 0.0016286 P/λ, with the atmospheric pressure P = 101.3 [(293 − 0.0065 z)/293]5.253 at the altitude z [m]

Vapor pressure

The saturation vapor pressure is the mean of e°(Tmax) and e°(Tmin), with e°(T) = 0.6108 exp[17.27 T/(T + 237.3)], and Δ the mean of the slopes 4098 e°(T)/(T + 237.3)² at Tmax and Tmin. The actual vapor pressure ea follows from the relative humidity or is given directly Allen et al. (1998).

Aerodynamic and canopy resistance

ra = ln[(zm − d)/zom] ln[(zh − d)/zoh] / (k² uz)
rc = rl / LAIactive = 200 / LAI

with the measurement heights of wind zm and humidity zh, the zero-plane displacement d = 0.667 hc, the roughness lengths zom = 0.123 hc and zoh = 0.1 zom of a crop of height hc, the von Kármán constant k = 0.41 and the wind speed uz; rl = 100 s/m and LAIactive = 0.5 LAI Allen et al. (1998). The LAI comes from the crop data or from the crop height (clipped grass: 0.24 hc with hc in cm; alfalfa: 5.5 + 1.5 ln hc).

Radiation

The net radiation is the difference of the net short-wave and the net long-wave radiation:

Rn = (1 − a) Rs − Rnl
Rs = (as + bs n/N) Ra
Rnl = σ [(Tmax,K4 + Tmin,K4)/2] (al + bl √ea) f

a is the albedo, Rs the solar radiation (measured, or from the relative sunshine duration n/N after Ångström (1924); as = 0.25, bs = 0.5), Ra the extraterrestrial radiation from the latitude and the day of the year, N the day length, σ the Stefan–Boltzmann constant, al + bl √ea the net emissivity (0.34, −0.14 Allen et al. (1998)) and f the cloudiness factor: f = a1 n/N + b1 from the sunshine duration (0.9, 0.1 Doorenbos and Pruitt (1977)), or f = ac Rs/Rso + bc from the ratio of the measured to the clear-sky radiation (1.35, −0.35 Allen et al. (1998)).

Partitioning

ETp is divided into potential transpiration and evaporation with the soil cover fraction SCF = 1 − exp(−k LAI) (partitioning) Ritchie (1972); the canopy interception is subtracted from the transpiration (interception).

See also: Hargreaves equation · surface energy balance · page Meteorological data

Literature cited in this topic

  • Allen, R. G., L. S. Pereira, D. Raes, and M. Smith (1998). Crop Evapotranspiration – Guidelines for Computing Crop Water Requirements. FAO Irrigation and Drainage Paper 56, Food and Agriculture Organization of the United Nations, Rome, 300 pp. https://www.fao.org/4/x0490e/x0490e00.htm
  • Doorenbos, J., and W. O. Pruitt (1977). Guidelines for Predicting Crop Water Requirements. FAO Irrigation and Drainage Paper 24, Food and Agriculture Organization of the United Nations, Rome, 144 pp.
  • Monteith, J. L. (1981). Evaporation and surface temperature. Quarterly Journal of the Royal Meteorological Society 107(451): 1–27. doi:10.1002/qj.49710745102
  • Monteith, J. L., and M. H. Unsworth (1990). Principles of Environmental Physics. 2nd ed., Edward Arnold, London, 291 pp.
  • Ritchie, J. T. (1972). Model for predicting evaporation from a row crop with incomplete cover. Water Resources Research 8(5): 1204–1213. doi:10.1029/WR008i005p01204
  • Ångström, A. (1924). Solar and terrestrial radiation. Quarterly Journal of the Royal Meteorological Society 50(210): 121–126. doi:10.1002/qj.49705021008

All sources: References.