Potential evapotranspiration

Thornthwaite equation


The Thornthwaite equation is a method developed by Thornthwaite (1948) based on an empirical approach in order to estimate potential evapotransipiration. Evapotranspiration is originally destinated to be computed on a monthly basis but daily estimations are possible as well (cf. Thornthwaite & Mather (1957).


The potential evapotranspiration according to the Thornthwaite (1948) formulation \(PET_{Thorn} \)[mm/day] is calculated as follows (the equation for \(T>26\) is an approximation of the original table proposed by Willmott et al. (1985) :

\[PET_{Thorn}= \begin{cases} 0, & \mbox{if } T < 0\\ {16}\cdot\frac {N}{360}\cdot\bigg(\dfrac{10\cdot T}{I}\bigg)^a, & \mbox{if }0 \leqslant T \leqslant 26 \\ \frac {N}{360}\cdot\big(-415.85+30.5332.24 \cdot T - 0.43 \cdot T^2\big), & \mbox{if } T > 26\\ \end{cases} \]

where \(N\) is the duration of sunlight in hours, varying with season and latitude (cf. calculation of daylight hours after the FAO formulation, or Table 8 in Thornthwaite & Mather (1957), \(T\) average daily air temperature [°C], and \(I\)a heat index calculated as follows:

\[I=\displaystyle \sum_{Jan}^{Dec} \bigg(\dfrac{max[0,T_m]}{5}\bigg)^{1.514}\]

where \(T_m\) is monthly mean temperature [°C].

The exponent \(a\) is calculated as follows:


An improved version of the daily Thornthwaite equation was proposed by Camargo et al. (1999) and Pereira & Pruitt (2004). Instead of average daily air temperature \(T\) [°C], Camargo et al. (1999) suggested to use an "effective temperature" \(T_{ef}\) [°C] given by:


where \(k=0.72\), \(T_{max}\) is maximum daily temperature [°C], and \(T_{min}\) minimum daily temperature [°C].

Pereira & Pruitt (2004) corrected this equation by including a day-night ratio:


\[\mbox{with }\dfrac{T_{max}+T_{min}}{2}\leqslant T_{ef}^* \leqslant T_{max}\]

where \(T_{ef}^*\) [°C] is the corrected \(T_{ef}\), and \(N\) the daylight hours.


Thornthwaite (1948)
Thornthwaite & Mather (1957)
Willmott et al. (1985)
Camargo et al. (1999)
Pereira & Pruitt (2004)

Penman equation


Penman (1948) developed an equation based on a more theoretical approach than Thornthwaite (1948) in order to estimate potential evaporation from open water. His equation has also been widely used to estimate evapotranspirations from vegetation covers. Shuttleworth (1993) reformulated the Penman equation in metric units.


The potential evapotranspiration after Penman (1948) \(PET_{Pen}\) [mm/day] is calculated as follows (Shuttleworth 1993):


where \(\Delta\) is the the slope of the saturation vapor pressure curve [kPa/°C], \(\gamma\) the psychrometric constant [kPa/°C], \(R_n\) the net radiation [MJ·m-2·d-1], \(\lambda\) the latent heat of vaporization, \(U\) the wind speed [m/s], and \(\Delta{e}\) the vapor pressure deficit [kPa].

The slope of the saturation vapor pressure \(\Delta\) is calculated as follows (Allen et al. 1998):


where \(T\) is air temperature [°C].

The latent heat of vaporization \(\lambda\) [MJ/kg] is calculated as follows (Shuttleworth 1993):


The psychrometric constant \(\gamma\) [kPa/°C] is calculated as follows (Shuttleworth 1993):


where \(p_{atm}\) is the atmospheric pressure [kPa] and calculated as follows (Allen et al. 1998):


where \(z\) is elevation above sea level [m].


Penman (1948)
Shuttleworth (1993)
Allen et al. (1998)


Variable Description Unit
\(T\) air temperature °C
\(T_{dew}\) dew point temperature °C
\(H\) air humidity %
\(P\) rainfall mm
\(U\) windspeed m/s
\(w\) days since last rain
(or rain above threshold)
\(rr\) days with consecutive rain d
\(\Delta t\) time increment d
\(\Delta{e}\) vapor pressure deficit kPa
\(e_s\) saturation vapor pressure kPa
\(e_a\) actual vapor pressure kPa
\(p_{atm}\) atmospheric pressure kPa
\( PET\) potential evapotranspiration mm/d
\(r\) soil water reserve mm
\(r_s\) surface water reserve mm
\(EMC\) equilibrium moisture content %
\(DF\) drought factor -
\(N\) daylight hours hr
\(D\) weighted 24-hr average moisture condition hr
\(\omega\) sunset hour angle rad
\(\delta\) solar declination rad
\(\varphi\) latitude rad
\(Cc\) cloud cover Okta
\(J\) day of the year (1..365/366) -
\(I\) heat index -
\(R_n\) net radiation MJ⋅m-2⋅d-1
\(R_a\) daily extraterrestrial radiation MJ⋅m-2⋅d-1
\(R_s\) solar radiation MJ⋅m-2⋅d-1
\(R_{so}\) clear-sky solar radiation MJ⋅m-2⋅d-1
\(R_{ns}\) net shortwave radiation MJ⋅m-2⋅d-1
\(R_{nl}\) net longwave radiation MJ⋅m-2⋅d-1
\(\lambda\) latent heat of vaporization MJ/kg
\(z\) elevation m a.s.l.
\(d_r\) inverse relative distance Earth-Sun -
\(\alpha\) albedo or canopy reflection coefficient -
\(\Delta\) slope of the saturation vapor pressure curve kPa/°C
\(Cc\) cloud cover eights
\(ROS\) rate of spread m/h
\(RSF\) rate of spread factor -
\(WF\) wind factor -
\(WRF\) water reserve factor -
\(FH\) false relative humidity -
\(FAF\) fuel availability factor -
\(PC\) phenological coefficient -

Suffix Description
\(-\) mean / daily value
\(_{max}\) maximum value
\(_{min}\) minimum value
\(_{12}\) value at 12:00
\(_{13}\) value at 13:00
\(_{15}\) value at 15:00
\(_{m}\) montly value
\(_{y}\) yearly value
\(_{f/a}\) value at fuel-atmosphere interface
\(_{dur}\) duration
\(_{soil}\) value at soil level

Constant Description
\(e\) Euler's number
\(\gamma\)psychrometric constant
\(G_{SC}\)solar constant
\(\sigma\)Stefan-Bolzmann constant