Penman–Monteith equation

The Penman–Monteith equation approximates net evapotranspiration (ET) from meteorological data, as a replacement for direct measurement of evapotranspiration. The equation is widely used, and was derived by the United Nations Food and Agriculture Organization for modeling reference evapotranspiration ET0.

Significance
Evapotranspiration contributions are very significant in a watershed's water balance, yet are often not emphasized in results because the precision of this component is often weak relative to more directly measured phenomena, e.g. rain and stream flow. In addition to weather uncertainties, the Penman–Monteith equation is sensitive to vegetation specific parameters, e.g. stomatal resistance or conductance.

Various forms of crop coefficients (Kc) account for differences between specific vegetation modeled and a reference evapotranspiration (RET or ET0) standard. Stress coefficients (Ks) account for reductions in ET due to environmental stress (e.g. soil saturation reduces root-zone O2, low soil moisture induces wilt, air pollution effects, and salinity). Models of native vegetation cannot assume crop management to avoid recurring stress.

Equation
Per Monteith's Evaporation and Environment, the equation is:
 * $$ \overset{\text{Energy flux rate}}{\lambda_v E=\frac{\Delta (R_n-G) + \rho_a c_p \left( \delta e \right) g_a }

{\Delta + \gamma \left ( 1 + g_a / g_s \right)}} ~ \iff ~ \overset{\text{Volume flux rate}}{ET=\frac{\Delta (R_n-G) + \rho_a c_p \left( \delta e \right) g_a } { \left( \Delta + \gamma \left ( 1 + g_a / g_s \right) \right) L_v }} $$


 * λv = Latent heat of vaporization. Energy required per unit mass of water vaporized. (J g−1)
 * Lv = Volumetric latent heat of vaporization. Energy required per water volume vaporized. (Lv = 2453 MJ m−3)


 * E = Mass water evapotranspiration rate (g s−1 m−2)
 * ET = Water volume evapotranspired (mm s−1)


 * &Delta; = Rate of change of saturation specific humidity with air temperature. (Pa K−1)
 * Rn = Net irradiance (W m−2), the external source of energy flux
 * G = Ground heat flux (W m−2), usually difficult to measure
 * cp = Specific heat capacity of air (J kg−1 K−1)
 * ρa = dry air density (kg m−3)
 * δe = vapor pressure deficit (Pa)
 * ga = Conductivity of air, atmospheric conductance (m s−1)
 * gs = Conductivity of stoma, surface or stomatal conductance (m s−1)
 * γ = Psychrometric constant (&gamma; ≈ 66 Pa K−1)

Note: Often resistances are used rather than conductivities.
 * $$ g_a = \tfrac{1}{ r_a} ~ ~ \And ~ ~ g_s = \tfrac{1}{ r_s} = \tfrac{1}{ r_c}$$

where rc refers to the resistance to flux from a vegetation canopy to the extent of some defined boundary layer.

The atmospheric conductance ga accounts for aerodynamic effects like the zero plane displacement height and the roughness length of the surface. The stomatal conductance gs accounts for effect of leaf density (Leaf Area Index), water stress and concentration in the air, that is to say to plant reaction to external factors. Different models exist to link the stomatal conductance to these vegetation characteristics, like the ones from P.G. Jarvis (1976) or Jacobs et al. (1996).

Accuracy
While the Penman-Monteith method is widely considered accurate for practical purposes and is recommended by the Food and Agriculture Organization of the United Nations, errors when compared to direct measurement or other techniques can range from -9 to 40%.

FAO 56 Penman-Monteith equation
To avoid the inherent complexity of determining stomatal and atmospheric conductance, the Food and Agriculture Organization proposed in 1998 a simplified equation for the reference evapotranspiration ET0. It is defined as the evapotranpiration for "[an] hypothetical reference crop with an assumed crop height of 0.12 m, a fixed surface resistance of 70 s m-1 and an albedo of 0.23." This reference surface is defined to represent "an extensive surface of green grass of uniform height, actively growing, completely shading the ground and with adequate water". The corresponding equation is:


 * $$ ET_o = \frac{0.408 \Delta (R_n-G) + \frac{900}{T} \gamma u_2 \delta e }{\Delta + \gamma (1 + 0.34 u_2)}

$$


 * ET0 = Reference evapotranspiration, Water volume evapotranspired (mm day−1)
 * &Delta; = Rate of change of saturation specific humidity with air temperature. (Pa K−1)
 * Rn = Net irradiance (MJ m−2 day−1), the external source of energy flux
 * G = Ground heat flux (MJ m−2 day−1), usually equivalent to zero on a day
 * T = Air temperature at 2m (K)
 * u2 = Wind speed at 2m height (m/s)
 * δe = vapor pressure deficit (kPa)
 * γ = Psychrometric constant (&gamma; ≈ 66 Pa K−1)

N.B.: The coefficient 0.408 and 900 are not unitless but account for the conversion from energy values to equivalent water depths: radiation [mm day−1] = 0.408 radiation [MJ m−2 day−1].

This reference evapotranspiration ET0 can then be used to evaluate the evapotranspiration rate ET from unstressed plant through crop coefficients Kc: ET = Kc * ET0.

Variations
The standard methods of the American Society of Civil Engineers modify the standard Penman–Monteith equation for use with an hourly time step. The SWAT model is one of many GIS-integrated hydrologic models estimating ET using Penman–Monteith equations.

Priestley–Taylor
The Priestley–Taylor equation was developed as a substitute to the Penman–Monteith equation to remove dependence on observations. For Priestley–Taylor, only radiation (irradiance) observations are required. This is done by removing the aerodynamic terms from the Penman–Monteith equation and adding an empirically derived constant factor, $$\alpha$$.

The underlying concept behind the Priestley–Taylor model is that an air mass moving above a vegetated area with abundant water would become saturated with water. In these conditions, the actual evapotranspiration would match the Penman rate of reference evapotranspiration. However, observations revealed that actual evaporation was 1.26 times greater than reference evaporation, and therefore the equation for actual evaporation was found by taking reference evapotranspiration and multiplying it by $$\alpha$$. The assumption here is for vegetation with an abundant water supply (i.e. the plants have low moisture stress). Areas like arid regions with high moisture stress are estimated to have higher $$\alpha$$ values.

The assumption that an air mass moving over a vegetated surface with abundant water saturates has been questioned later. The lowest and turbulent part of the atmosphere, the atmospheric boundary layer, is not a closed box, but constantly brings in dry air from higher up in the atmosphere towards the surface. As water evaporates more easily into a dry atmosphere, evapotranspiration is enhanced. This explains the larger than unity value of the Priestley-Taylor parameter $$\alpha$$. The proper equilibrium of the system has been derived and involves the characteristics of the interface of the atmospheric boundary layer and the overlying free atmosphere.

History
The equation is named after Howard Penman and John Monteith. Penman first published his equation in 1948 and Monteith revised it in 1965.