Skip to main content
Geosciences LibreTexts

12.7: Reference evapotranspiration

  • Page ID
    38828
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

    ( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\id}{\mathrm{id}}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\kernel}{\mathrm{null}\,}\)

    \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\)

    \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\)

    \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    \( \newcommand{\vectorA}[1]{\vec{#1}}      % arrow\)

    \( \newcommand{\vectorAt}[1]{\vec{\text{#1}}}      % arrow\)

    \( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vectorC}[1]{\textbf{#1}} \)

    \( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

    \( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

    \( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

    \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)

    One of the most widespread uses of the surface energy balance concept is for determining evapotranspiration. The actual evapotranspiration rate for a given time for a particular location on the land surface can be substantially influenced by site-specific management practices. Thus, often it is useful to estimate the potential or reference evapotranspiration(ET0)which is the evapotranspiration rate expected for a standard well-watered surface under the prevailing atmospheric conditions. The daily or monthly value of ET0 can be determined for relatively large regions, and those values can be adjusted to estimate actual ET for a specific location if adequate information is available [17]. Methods for estimating ET0 or ET commonly are based on the surface energy balance.

    12.7.1 Hargreaves method

    One of the simplest approaches for estimating ET0, the Hargreaves method,was developed to facilitate irrigation management and water planning. The initial work began in California, and the importance of having a simple approach with low data requirements was solidified when the work expanded to Haiti [18]. Building on earlier work in which ET0 was calculated as a linear function of incoming solar radiation (Rs), Hargreaves found that the product of Rs and air temperature (T) explained much of the variance in measured ET0. The problem was that measurements of Rs are often unavailable. Therefore, Hargreaves and colleagues developed a way to estimate Rs from knowledge of extraterrestrial radiation (Ra, solar radiation at the top of Earth’s atmosphere for a specific place and time) and air temperature[19].Since Ra can be accurately estimated based on site latitude and day of year and air temperature measurements are widely available, the Hargreaves method can be used virtually anywhere. The Hargreaves equation is

    formula Equation 12-5.png

    (Equation 12-5)

    where Ra is the cumulative extraterrestrial radiation for the calculation period converted to units of mm of water, T is the mean temperature (°C) for the calculation period, Tmax is mean daily maximum temperature for the calculation period, and Tmin is mean daily minimum temperature for the calculation period. This equation is designed for use at 5-day or larger calculation periods and is frequently used with monthly calculation periods. The Hargreaves method is a good choice for ET0 estimation when the only meteorological data available are air temperatures. For example, the Hargreaves method was used to estimate a 45 -yr time series of ET0 for a large watershed stretching across southwest Oklahoma and part of west Texas, USA ( Fig. 12-7) [20]. Compiling long time series like this would be more difficult, or in some cases impossible, for more data intensive ET0 estimation methods.

    figure of the annual reference evapotranspiration
    Figure 12-7. Annual reference evapotranspiration (ET0) totals from 1970-2014 for the North Fork of the Red River watershed which lies in southwest Oklahoma and the Texas panhandle. Monthly ET0 values were calculated using the Hargreaves method and summed to get the annual totals. Adapted from Krueger et al. (2017)

    12.7.2 Penman-Monteith method

    A more precise approach was for estimating ET was proposed by Penman, who in 1948 combined the equations for the surface energy balance and for turbulent transport of water vapor away from a saturated evaporating surface[21]. In 1965, Monteith expanded the Penman approach to surfaces where the vapor pressure was less than the saturated vapor pressure [22].The Penman-Monteith equation is

    fraction Equation 12-6.png

    (Equation 12-6)

    where L is the latent heat of vaporization for water (J kg- 1), Rn is the net radiation (W m-2), G is the soil heat flux (W m -2), es-ea is the vapor pressure deficit of the air (kPa), ρa is the density of the air (kg m-3), cp is the specific heat capacity of the air at constant pressure (J kg-1°C-1), Δ is the slope of the relationship between saturated vapor pressure and temperature (kPa °C-1), y is called the psychrometric constant (kPa °C-1), rs is the surface resistance to vapor transport (s m-1), and ra is the aerodynamic resistance to vapor transport (s m-1). These two resistance terms provide a simplified but effective representation of water vapor transport processes near the land surface (Fig. 12-8).

    figure of conceptual representation of near surface
    Figure 12-8. Conceptual representation of near surface water vapor transport as controlled by the surface and aerodynamic resistances. Reproduced from Allen et al. (1998)

    A panel of experts from the United Nations Food and Agricultural Organization(FAO) has defined the reference surface as, “A hypothetical reference crop with an assumed height of 0.12 m, a fixed surface resistance of 70 s m-1, and an albedo of 0.23.” [17]. Based on this definition, and employing some approximations for the resistance terms, the FAO Penman-Monteith equation for ET0 is

    fraction Equation 12-7.png

    (Equation 12-7)

    where ET0 is the reference evapotranspiration (mm d-1), Rn and G are in MJ m -2 d-1, T is the mean daily air temperature at 2-m height (°C), and u 2 is the mean daily wind speed at 2 m (m s -1). Once ET 0 is determined, site - and day - specific “crop coefficients” (Kc) can be estimated to calculate actual ET from ET 0, i.e. ET = ET0 x Kc. With judicious use, the FAO Penman- Monteith equation has proven to provide reliable ET0 estimates for use in crops, grasslands, and forests around the world. For example, the FAO Penman - Monteith equation was used to estimate daily ET0 during a 3- yr bioenergy crop experiment in Oklahoma (Fig.12-9). The results show that actual ET can be substantially less than ET0 during dry periods (e.g. summers of 2011 and 2012) but can exceed ET 0 during wet periods (e.g. summer of 2013) [23].

    figure of reference evapotranspiration
    Figure 12-9. Reference evapotranspiration (shaded region) and actual evapotranspiration under switchgrass (SWG), biomass sorghum (BMS), and mixed perennial grasses (MXG) at Stillwater, Oklahoma. Reproduced from Yimam et al. (2014).

    This page titled 12.7: Reference evapotranspiration is shared under a CC BY 4.0 license and was authored, remixed, and/or curated by Tyson Oschner.