PenmanMonteith equation
Calculation procedures with missing data
Pan evaporation method
This chapter demonstrates how the crop reference evapotranspiration (ET_{o}) is determined either from meteorological data or from pan evaporation.
The FAO PenmanMonteith method is maintained as the sole standard method for the computation of ET_{o} from meteorological data. The method itself is introduced in Chapter 2, and the computation of all data required for the calculation of ET_{o} is discussed in Chapter 3. This chapter presents guidelines to calculate ET_{o} with different time steps, ranging from hours to months, and with missing climatic data. The ET_{o} calculation can be done by hand with the help of a calculation sheet, or by means of a computer.
ET_{o} can also be estimated from the evaporation loss from a water surface. The procedure to obtain ET_{o} from pan evaporation and the coefficients for different types of pans are presented in this chapter.
Calculation procedure
ET_{o} calculated with different time steps
From the original PenmanMonteith equation and the equations of the aerodynamic and canopy resistance, the FAO PenmanMonteith equation has been derived in Chapter 2:
_{} (6)
where
ET_{o} reference evapotranspiration [mm day^{1}],
R_{n} net radiation at the crop surface [MJ m^{2} day^{1}],
G soil heat flux density [MJ m^{2} day^{1}],
T air temperature at 2 m height [°C],
u_{2} wind speed at 2 m height [m s^{1}],
e_{s} saturation vapour pressure [kPa],
e_{a} actual vapour pressure [kPa],
e_{s}  e_{a} saturation vapour pressure deficit [kPa],
D slope vapour pressure curve [kPa °C^{1}],
g psychrometric constant [kPa °C^{1}].
The FAO PenmanMonteith equation determines the evapotranspiration from the hypothetical grass reference surface and provides a standard to which evapotranspiration in different periods of the year or in other regions can be compared and to which the evapotranspiration from other crops can be related.
Calculation sheet
ET_{o} can be estimated by means of the calculation sheet presented in Box 11. The calculation sheet refers to tables in Annex II for the determination of some of the climatic parameters. The calculation procedure consists of the following steps:
1. Derivation of some climatic parameters from the daily maximum (T_{max}) and minimum (T_{min}) air temperature, altitude (z) and mean wind speed (u_{2}).2. Calculation of the vapour pressure deficit (e_{s}  e_{a}). The saturation vapour pressure (e_{s}) is derived from T_{max} and T_{min}, while the actual vapour pressure (e_{a}) can be derived from the dewpoint temperature (T_{dew}), from maximum (RH_{max}) and minimum (RH_{min}) relative humidity, from the maximum (RH_{max}), or from mean relative humidity (RH_{mean}).
3. Determination of the net radiation (R_{n}) as the difference between the net shortwave radiation (R_{ns}) and the net longwave radiation (R_{nl}). In the calculation sheet, the effect of soil heat flux (G) is ignored for daily calculations as the magnitude of the flux in this case is relatively small. The net radiation, expressed in MJ m^{2} day^{1}, is converted to mm/day (equivalent evaporation) in the FAO PenmanMonteith equation by using 0.408 as the conversion factor within the equation.
4. ET_{o} is obtained by combining the results of the previous steps.
Examples 17 and 20 present typical examples using the calculation sheet.
Computerized calculations
Calculations of the reference crop evapotranspiration ET_{o} are often computerized. The calculation procedures of all data required for the calculation of ET_{o} by means of the FAO PenmanMonteith equation are presented in Chapter 3. Typical sequences in which the calculations can be executed are given in the calculation sheets. The procedures presented in Boxes 7 (vapour pressure deficit), 9 (extraterrestrial radiation and daylight hours), 10 (net radiation) and 11 (ET_{o}) can be used when developing a spreadsheet or computer program to calculate ET_{o}.
Many software packages already use the FAO PenmanMonteith equation to assess the reference evapotranspiration. As an example, the output of CROPWAT, the FAO software for irrigation scheduling, is presented in Figure 18.
The selection of the time step with which ET_{o} is calculated depends on the purpose of the calculation, the accuracy required and the time step of the climatic data available.
BOX 11. Calculation sheet for ET_{o} (FAO PenmanMonteith) using meteorological tables of Annex 2  
Parameters  
T_{max} 

°C 
 
T_{min} 

°C 
T_{mean}  (T_{max} + T_{min})/2 

°C 
T_{mean} 

°C 
D (Table 2.4 of Annex 2) 

kPa/°C 
Altitude 

m 
g (Table 2.2 of Annex 2) 

kPa/°C 
u_{2} 

m/s 
(1 + 0.34 u_{2}) 


D /[D + g (1 + 0.34 u_{2})] 

 
g /[D + g (1 + 0.34 u_{2})] 

 
[900/(T_{mean} + 273)] u_{2} 

 
Vapour pressure deficit  
T_{max} 

°C 
e°(T_{max}) (Table 2.3) 

kPa 
T_{min} 

°C 
e°(T_{min}) (Table 2.3) 

kPa 
Saturation vapour pressure e_{s} = [(e°(T_{max}) + e°(T_{min})]/2 

kPa  
e_{a} derived from dewpoint temperature:  
T_{dew} 

°C 
e_{a} = e°(T_{dew}) (Table 2.3) 

kPa 
OR e_{a} derived from maximum and minimum relative humidity:  
RH_{max} 

% 
e°(T_{min}) RH_{max}/100 

kPa 
RH_{min} 

% 
e°(T_{max}) RH_{min}/100 

kPa 

e_{a}: (average) 

kPa  
OR e_{a} derived from maximum relative humidity: (recommended if there are errors in RH_{min})  
RH_{max} 

% 
e_{a} = e°(T_{min}) RH_{max}/100 

kPa 
OR e_{a} derived from mean relative humidity: (less recommended due to nonlinearities)  
RH_{mean} 

% 
e_{a} = e_{s} RH_{mean}/100 

kPa 
Vapour pressure deficit (e_{s}  e_{a}) 

kPa  
Radiation  
Latitude 

° 



Day 


R_{a} (Table 2.6) 

MJ m^{2} d^{1} 
Month 


N (Table 2.7) 

hours 
n 

hours 
n/N 


If no R_{s} data available: R_{s} = (0.25 + 0.50 n/N) R_{a} 

MJ m^{2} d^{1}  
R_{so} = [0.75 + 2 (Altitude)/100000] R_{a} 

MJ m^{2} d^{1}  
R_{s}/R_{so} 

 
R_{ns} = 0.77 R_{s} 

MJ m^{2} d^{1}  
T_{max} 


_{} (Table 2.8) 

MJ m^{2} d^{1} 
T_{min} 


_{} (Table 2.8) 

MJ m^{2} d^{1} 
_{} 

MJ m^{2} d^{1}  
e_{a} 

kPa 
(0.340.14 Ö e_{a}) 


R_{s}/R_{so} 


(1.35 R_{s}/R_{so}  0.35) 


_{} 

 
R_{n} = R_{ns}  R_{nl} 

 
T_{month} 

°C 
G_{day} (assume) 
0 

T_{month1} 

°C 
G_{month} = 0.14 (T_{month}  T_{month1}) 


R_{n}  G 

MJ m^{2} d^{1}  
0.408 (R_{n}  G) 

mm/day  
Grass reference evapotranspiration  
_{} 

mm/day  
_{} 

mm/day  
_{} 

mm/day 
FIGURE 18. ET_{o} computed by CROPWAT
MONTHLY REFERENCE EVAPOTRANSPIRATION PENMAN MONTEITH  
Mateostation: CABINDA 
Country: Angola  
Altitude: 20 m. 
Coordinates: 5.33 
South 
12.11 East  
Month 
MinTemp 
MaxTemp 
Humid. 
Wind 
Sunshine 
Radiation 
EToPenMon 

°C 
°C 
% 
km/day 
Hours 
MJ/m^{2}/day 
mm/day 
January 
22.8 
29.6 
81 
78 
4.0 
15.7 
3.4 
February 
22.7 
30.3 
82 
69 
4.6 
16.9 
3.7 
March 
23.0 
30.6 
80 
78 
5.1 
17.4 
3.8 
April 
23.0 
30.2 
82 
69 
5.0 
16.4 
3.5 
May 
22.0 
28.6 
84 
69 
3.8 
13.5 
2.9 
June 
19.2 
26.5 
81 
69 
3.3 
12.2 
2.6 
July 
17.6 
25.1 
78 
78 
3.2 
12.3 
2.6 
August 
18.6 
25.3 
78 
78 
2.6 
12.4 
2.6 
September 
20.5 
26.5 
78 
104 
2.0 
12.4 
2.8 
October 
22.5 
28.0 
79 
130 
2.2 
12.9 
3.1 
November 
23.0 
28.7 
80 
104 
3.2 
14.4 
3.3 
December 
23.0 
29.1 
82 
95 
3.8 
15.2 
3.4 
Year 
21.5 
28.2 
80 
85 
3.6 
14.3 
3.1 
CROPWAT 7.0 Climate file: C:\PROFP~1\CROPWAT\CROPWAT\CLI\CABINDA.PEN 03/07/98
Tenday or monthly time step
Notwithstanding the nonlinearity in the PenmanMonteith equation and some weather parameter methods, mean tenday or monthly weather data can be used to compute the mean tenday or monthly values for the reference evapotranspiration. The value of the reference evapotranspiration calculated with mean monthly weather data is indeed very similar to the average of the daily ET_{o} values calculated with daily average weather data for that month.
The meteorological data consist of:
· Air temperature: tenday or monthly average daily maximum (T_{max}) and average daily minimum temperature (T_{min}).· Air humidity: tenday or monthly average of the daily actual vapour pressure (e_{a}) derived from psychrometric, dewpoint or relative humidity data.
· Wind speed: tenday or monthly average of daily wind speed data measured at 2 m height (u_{2}).
· Radiation: tenday or monthly average of daily net radiation (R_{n}) computed from the mean tenday or monthly measured shortwave radiation or from actual duration of daily sunshine hours (n). The extraterrestrial radiation (R_{a}) and daylight hours (N) for a specific day of the month can be computed using Equations 21 and 34 or can be selected from Tables 2.5 and 2.6 in Annex 2.
When the soil is warming (spring) or cooling (autumn), the soil heat flux (G) for monthly periods may become significant relative to the mean monthly R_{n}. In these cases G cannot be ignored and its value should be determined from the mean monthly air temperatures of the previous and next month. Chapter 3 outlines the calculation procedure (Equations 43 and 44).
EXAMPLE 17. Determination of ET_{o} with mean monthly data
Given the monthly average climatic data of April of Bangkok (Thailand) located at 13°44'N and at an elevation of 2 m: 

 
Monthly average daily maximum temperature (T_{max}) = 
34.8 
°C 
 
Monthly average daily minimum temperature (T_{min}) = 
25.6 
°C 
 
Monthly average daily vapour pressure (e_{a}) = 
2.85 
kPa 
Measured at 2 m 
Monthly average daily wind speed (u_{2}) = 
2 
m/s 
 
Monthly average sunshine duration (n) = 
8.5 
hours/day 
For April 
Mean monthly average temperature (T_{month, i}) = 
30.2 
°C 
For March 
Mean monthly average temperature (T_{month, i1}) = 
29.2 
°C 
Determination according to outline of Box 11 (calculation sheet ET_{o}) 

Parameters 



 
T_{mean} = [(T_{max} = 34.8) + (T_{min} = 25.6)]/2 = 
30.2 
°C 
From Table 2.4 or Eq.13: 
D = 
0.246 
kPa/°C 
From Table 2.1 and Table 2.2 or Eq. 7 and Eq. 8: 
Altitude = 
2 
m 
P= 
101.3 
kPa 

g = 
0.0674 
kPa/°C 

 
(1 + 0.34 u_{2}) = 
1.68 
 
 
D /[D + g (1 + 0.34u_{2})] = 0.246/[(0.246 + 0.0674 (1.68)] = 
0.685 
 
 
g /[D + g (1 + 0.34u_{2})] = 0.0667/[0.246 + 0.0674 (1.68)] = 
0.188 
 
 
900/(T_{mean} + 273) u_{2} = 
5.94 
 
Vapour pressure deficit 

From Table 2.3 or 
T_{max} = 
34.8 
°C 
Eq. 11: 
e°(T_{max}) = 
5.56 
kPa 
From Table 2.3 or 
T_{min} = 
25.6 
°C 
Eq. 11: 
e°(T_{min}) = 
3.28 
kPa 
 
e_{s} = (5.56 + 3.28)/2 = 
4.42 
kPa 
Given 
e_{a} = 
2.85 
kPa 
 
Vapour pressure deficit (e_{s}  e_{a}) = (4.42  2.85) = 
1.57 
kPa 
Radiation (for month = April) 

From Table 2.6 or 2.5 or Eq. 21: 
J = (for 15 April) 
105 
 
Latitude = 13°44'N = (13 + 44/60) = 
13.73 
°N 

R_{a} = 
38.06 
MJ m^{2} day^{1} 

N (Table 2.7 or Eq. 34): 
Daylength N = 
12.31 
hours 
 
n/N = (8.5/12.31) = 
0.69 
 
 
R_{s} = [0.25 + 0.50 (0.69)] 38.06 = 
22.65 
MJ m^{2} day^{1} 
 
R_{so} = (0.75 + 2 (2)/100000) 38.06 = 
28.54 
MJ m^{2} day^{1} 
 
R_{s}/R_{so} = (22.65/28.54) = 
0.79 
 
 
R_{ns} = 0.77 (22.65) = 
17.44 
MJ m^{2} day^{1} 
From Table 2.8: 
T_{max} = 
34.8 
°C 
_{} 
44.10 
MJ m^{2} day^{1} 

From Table 2.8: 
T_{min} = 
25.6 
°C 
_{} 
39.06 
MJ m^{2} day^{1} 

 
_{} 
41.58 
MJ m^{2} day^{1} 
For: 
e_{a} = 
2.85 
kPa 
Then: 
(0.340.14Ö e_{a}) = 
0.10 
 
For: 
R_{s}/R_{so} = 
0.79 
 
Then: 
(1.35 R_{s}/R_{so}0.35) = 
0.72 
 
 
R_{nl} = 41.58 (0.10)0.72 = 
3.11 
MJ m^{2} day^{1} 
 
R_{n} = (17.443.11) = 
14.33 
MJ m^{2} day^{1} 
 
G =0.14 (30.229.2) = 
0.14 
MJ m^{2} day^{1} 
 
(R_{n}  G) = (14.330.14) = 
14.19 
MJ m^{2} day^{1} 
 
0.408 (R_{n}  G) = 
5.79 
mm/day 
Grass reference evapotranspiration 

 
0.408 (R_{n}  G) D /[D + g (1+0.34 u_{2})] = 


 
(5.79) 0.685 = 
3.97 
mm/day 
 
900 u_{2}/(T + 273) (e_{s}  e_{a}) g /[D + g (1+0.34 u_{2})] = 5.94(1.57)0.188 = 
1.75 
mm/day 
 
ET_{o} = (3.97+1.75) = 
5.72 
mm/day 
The grass reference evapotranspiration is 5.7 mm/day. 
Daily time step
Calculation of ET_{o} with the PenmanMonteith equation on 24hour time scales will generally provide accurate results. The required meteorological data consist of:
· Air temperature: maximum (T_{max}) and minimum (T_{min}) daily air temperatures.· Air humidity: mean daily actual vapour pressure (e_{a}) derived from psychrometric, dewpoint temperature or relative humidity data.
· Wind speed: daily average for 24 hours of wind speed measured at 2 m height (u_{2}).
· Radiation: net radiation (R_{n}) measured or computed from solar and longwave radiation or from the actual duration of sunshine (n). The extraterrestrial radiation (R_{a}) and daylight hours (N) for a specific day of the month should be computed using Equations 21 and 34. As the magnitude of daily soil heat flux (G) beneath the reference grass surface is relatively small, it may be ignored for 24hour time steps.
EXAMPLE 18. Determination of ET_{o} with daily data
Given the meteorological data as measured on 6 July in Uccle (Brussels, Belgium) located at 50°48'N and at 100 m above sea level: 

 
Maximum air temperature (T_{max}) = 
21.5 
°C 
 
Minimum air temperature (T_{min}) = 
12.3 
°C 
 
Maximum relative humidity (RH_{max}) = 
84 
% 
 
Minimum relative humidity (RH_{min}) = 
63 
% 
 
Wind speed measured at 10 m height = 
10 
km/h 
 
Actual hours of sunshine (n) = 
9.25 
hours 
Conversion of wind speed 

At 10 m height 
Wind speed = 10 km/h or u_{z} = 
2.78 
m/s 
From Eq. 47, with z = 10 m: 
At standard height, u_{2} = 0.748 (2.78) = 
2,078 
m/s 
Parameters 

From Eq. 7, for: 
altitude = 
100 
m 
 
P = 
100.1 
kPa 
 
T_{mean} = (21.5 + 12.3)/2 = 
16.9 
°C 
From Eq. 13, for: 
T_{mean} = 
16.9 
°C 
D = 
0.122 
kPa/°C 

From Eq. 8, for: 
P = 
100.1 
kPa 
g = 
0.0666 
kPa/°C 

 
(1 + 0.34 u_{2}) = 
1.71 
 
 
D /[D + g (1+0.34 u_{2})] = 0.122/[(0.122 + 0.0666 (1.71)] = 
0.518 
 
 
g /[D + g (1+0.34 u_{2})] = 0.0666/[0.122 + 0.0666 (1.71)] = 
0.282 
 
 
900/(T_{mean} + 273) u_{2} = 
6.450 
 
Vapour pressure deficit 

From Eq. 11, for: 
T_{max} = 
21.5 
°C 
Then: 
e°(T_{max}) = 
2.564 
kPa 
From Eq. 11, for: 
T_{min} = 
12.3 
°C 
Then: 
e°(T_{min}) = 
1.431 
kPa 
 
e_{s} = (2.564 + 1.431) = 
1.997 
kPa 
Given relative 
RH_{max} = 
84 
% 
humidity data 
RH_{min} = 
63 
% 
From Eq. 17: 
e_{a} = [1.431 (0.84) + 2.564 (0.63)]/2 = 
1.409 
kPa 
 
Vapour pressure deficit (e_{s}  e_{a}) = (1.9971.409) = 
0.589 
kPa 
Radiation 

From Table 2.5: 
Month 7, Day = 6 


J = 
187 
 

From Eq. 21: 
Latitude = 50°48'N = 
50.80 
°N 
J = 
187 
 

R_{a} = 
41.09 
MJ m^{2} day^{1} 

From Eq. 34: 
Latitude = 50°48'N = 
50.80 
°N 
J = 
187 
 

N = 16.1 
16.1 
hours 

n/N = 9.25/16.3 = 
0.57 
 

From Eq. 35 
R_{s} = [0.25 + 0.50 (0.57)] 41.09 
22.07 
MJ m^{2} day^{1} 
From Eq. 37 
R_{so} = (0.75 +2(100)/100000) 41.09 = 
30.90 
MJ m^{2} day^{1} 
 
R_{s}/R_{so} = 
0.71 
 
From Eq. 38 
R_{ns} = 0.77 (22.07) = 
17.00 
MJ m^{2} day^{1} 
For: 
T_{max} = 
21.5 
°C 
Then: 
T_{max, K} = 21.5+273.16 = 
294.7 
K 
_{} 
36.96 
MJ m^{2} day^{1} 

For: 
T_{min} = 
12.3 
°C 
Then: 
T_{min, K} = 12.3 +273.16 = 
285.5 
K 
_{} 
32.56 
MJ m^{2} day^{1} 

_{} 
34.76 
MJ m^{2} day^{1} 

 
(0.34  0.14Ö e_{a}) = 
0.17 
 
 
(1.35(R_{s}/R_{so})  0.35) = 
0.61 
 
From Eq. 39 
R_{nl} = 34.76(0.17)0.61 = 
3.71 
MJ m^{2} day^{1} 
From Eq. 40 
R_{n} = (17.00  3.71) = 
13.28 
MJ m^{2} day^{1} 
From Eq. 42 
G = 
0 
MJ m^{2} day^{1} 
 
(R_{n}  G) = (13.28  0) = 
13.28 
MJ m^{2} day^{1} 
 
0.408 (R_{n}  G) = 
5.42 
mm/day 
Grass reference evapotranspiration 

 
0.408 (R_{n}  G) D /[D + g (1 + 0.34 u_{2})] = 
2.81 
mm/day 
 
900/(T+273) u_{2} (e_{s}  e_{a}) g /[D + g (1+0.34 u_{2})] = 
1.07 
mm/day 
 
ET_{o} (Eq. 6) = 2.81 + 1.07 = 3.88 » 
3.9 
mm/day 
The grass reference evapotranspiration is 3.9 mm/day. 
Hourly time step
In areas where substantial changes in wind speed, dewpoint or cloudiness occur during the day, calculation of the ET_{o} equation using hourly time steps is generally better than using 24hour calculation time steps. Such weather changes can cause 24hour means to misrepresent evaporative power of the environment during parts of the day and may introduce error into the calculations. However, under most conditions, application of the FAO PenmanMonteith equation with 24hour data produces accurate results.
With the advent of electronic, automated weather stations, weather data are increasingly reported for hourly or shorter periods. Therefore, in situations where calculations are computerized, the FAO PenmanMonteith equation can be applied on an hourly basis with good results. When applying the FAO PenmanMonteith equation on an hourly or shorter time scale, the equation and some of the procedures for calculating meteorological data should be adjusted for the smaller time step. The FAO PenmanMonteith equation for hourly time steps is:
_{} (53)
where
ET_{o} reference evapotranspiration [mm hour^{1}],
R_{n} net radiation at the grass surface [MJ m^{2} hour^{1}] (Equation 40),
G soil heat flux density [MJ m^{2} hour^{1}] (Equations 45 and 46),
T_{hr} mean hourly air temperature [°C],
D saturation slope vapour pressure curve at T_{hr} [kPa °C^{1}] (Equation 13),
g psychrometric constant [kPa °C^{1}] (Equation 8),
e°(T_{hr}) saturation vapour pressure at air temperature T_{hr} [kPa] (Equation 11),
e_{a} average hourly actual vapour pressure [kPa] (Equation 54),
u_{2} average hourly wind speed [m s^{1}].
Given relative humidity measurements, the actual vapour pressure is determined as:
_{} (54)
where
e_{a} average hourly actual vapour pressure [kPa],
e°(T_{hr}) saturation vapour pressure at air temperature T_{hr} [kPa] (Equation 11),
RH_{hr} average hourly relative humidity [%].
The net radiation is the difference between the net shortwave radiation (R_{ns}) and the net longwave radiation (R_{nl}) at the hourly time steps. Consequently:
· If R_{ns} and R_{nl} need to be calculated, the extraterrestrial radiation value (R_{a}) for the hourly period (Equation 28) should be used.· In the computation of R_{nl} by means of Equation 39, _{} is replaced by _{} and the StefanBoltzman constant becomes:
s = (4.903/24) 10^{9} = 2.043 10^{10} MJ m^{2} hour^{1}.
Since the ratio R_{s}/R_{so} is used to represent cloud cover, when calculating R_{nl} for hourly periods during the nighttime, the ratio R_{s}/R_{so} can be set equal to the R_{s}/R_{so} calculated for a time period occurring 23 hours before sunset, before the sun angle becomes small. This will generally serve as a good approximation of cloudiness occurring during the subsequent nighttime. The hourly period that is 2 to 3 hours before sunset can be identified during computation of R_{a} as the period where w, calculated from Equation 31, is within the range (w _{s}  0.79) £ w £ (w _{s}  0.52), where w _{s} is calculated using Equation 25. As a more approximate alternative, one can assume R_{s}/R_{so} = 0.4 to 0.6 during nighttime periods in humid and subhumid climates and R_{s}/R_{so} = 0.7 to 0.8 in arid and semiarid climates. A value of R_{s}/R_{so} = 0.3 presumes total cloud cover.
Soil heat flux is important for hourly calculations. Equations 45 and 46 can be used to derive G for the hourly periods.
The required meteorological data consist of:
· Air temperature: mean hourly temperature (T_{hr}).
· Air humidity: average hourly relative humidity (RH_{hr}).
· Wind speed: average hourly wind speed data measured at 2 m height (u_{2}).
· Radiation: total hourly solar (R_{s}) or net radiation (R_{n}).
Because of the need for standardization, the constants in Equation 53 presume a constant surface resistance (r_{s}) of 70 s/m during all periods. This constant resistance may cause some underprediction of hourly ET_{o} during some daytime periods when actual r_{s} may be somewhat lower. The constant resistance may cause some overprediction of hourly ET_{o} during evening periods when actual r_{s} may be somewhat higher. However, when the calculations of hourly ET_{o} from Equation 53 are summed over 24 hour periods to produce an equivalent 24hour ET_{o}, the hourly differences tend to compensate one another and the results are generally equivalent to calculations of ET_{o} made on a 24hour time step. Precise estimates of ET_{o} for specific hourly periods may require the use of aerodynamic stability functions and functions for modifying the value for r_{s} based on levels of radiation, humidity and temperature. Application of these functions are not normally required when hourly calculations are to be summed to 24hour totals. Therefore, these functions are not described here.
EXAMPLE 19. Determination of ET_{o} with hourly data
Given mean average hourly data between 02.00 and 03.00 hours and 14.00 and 15.00 hours on 1 October in N'Diaye (Senegal) at 16°13'N and 16°15'W and 8 m above sea level. In the absence of calibrated coefficients, indicative values for a_{s} and b_{s} (Eq. 35 Angstrom formula) and for the coefficients of the net longwave radiation (Eq. 39) are used. 

Measured climatic data 
02.0003.00 h 
14.0015.00h 
Units 

T_{hr}: mean hourly temperature = 
28 
38 
°C 

RH_{hr}: mean hourly relative humidity = 
90 
52 
% 

u_{2}: mean hourly wind speed = 
1.9 
3.3 
m/s 

R_{s}: total solar radiation = 
 
2.450 
MJ m^{2} hour^{1} 

Parameters 

From Eq. 13 
D = 
0.220 
0.358 
kPa °C^{1} 
From Eq. 8 
g = 
0.0673 
0.0673 
kPa °C^{1} 
Vapour pressure deficit 

From Eq. 11 
e°(T) = 
3.780 
6.625 
kPa 
From Eq. 54 
e_{a} = 
3.402 
3.445 
kPa 
 
e_{s}  e_{a} = 
0.378 
3.180 
kPa 
Extraterrestrial radiation 
02.0003.00 h 
14.0015.00h 
Units 

From Table 2.5 



for 1 October: 
J = 274 
 

From Eq. 22: 
j = p /180 (16.22) = 0.2830 
rad 

From Eq. 23: 
d_{r} = 1.0001 
 

From Eq. 24: 
d =  0.0753 
rad 

From Eq. 33: 
b = 3.3315 
 

From Eq. 32: 
S_{c} = 0.1889 
hour 

 
L_{z} = 15 
degrees 

 
L_{m} = 16.25 
degrees 

 
t = 
2.5 
14.5 
hour 
From Eq. 31: 
w = 
2.46 
0.682 
rad 
 
t_{l} = 
1 
1 
hour 
From Eq. 29: 
w _{1} = 
 
0.5512 
rad 
From Eq. 30: 
w _{2} = 
 
0.8130 
rad 
From Eq. 28: 
R_{a} = 
0 
3.543 
MJ m^{2} hour^{1} 
Radiation 

Given 
R_{s} = 
0 
2.450 
MJ m^{2} hour^{1} 
From Eq. 37: 
R_{so} = 
0 
2.658 
MJ m^{2} hour^{1} 
From Eq. 38: 
R_{ns} = 
0 
1.887 
MJ m^{2} hour^{1} 
 
_{} 
1.681 
1.915 
MJ m^{2} hour^{1} 
 
(0.340.14 Ö e_{a}) = 
0.082 
0.080 
 
 
R_{s}/R_{so} = 
0.8 _{ (assumed)} 
0.922 
 
 
(1.35 R_{s}/R_{so}  0.35) = 
0.730 
0.894 
 
From Eq. 39: 
R_{nl} = 
0.100 
0.137 
MJ m^{2} hour^{1} 
From Eq. 40: 
R_{n} = 
0.100 
1.749 
MJ m^{2} hour^{1} 
From Eq. 46, 45: 
G = 
0.050 
0.175 
MJ m^{2} hour^{1} 
 
(R_{n}  G) = 
0.050 
1.574 
MJ m^{2} hour^{1} 
 
0.408(R_{n}  G) = 
0.020 
0.642 
mm/hour 
Grass reference evapotranspiration 

 
0.408(R_{n}  G) 



 
D /[D + g (1+0.34 u_{2})] = 
0.01 
0.46 
mm/hour 
37/(T + 273) u_{2} (e_{s}  e_{a}) 




From Eq. 53: 
g /[D + g (1 +0.34 u_{2})] = 
0.01 
0.17 
mm/hour 
ET_{o} = 
0.00 
0.63 
mm/hour 

The grass reference evapotranspiration is 0.00 mm/hour between 02.00 and 03.00 hours and 0.63 mm/hour between 14.00 and 15.00 hours. 
The meteorological data, required to estimate ET_{o} by means of the FAO PenmanMonteith equation, consist of air temperature, air humidity, wind speed and radiation. Where some of these data are missing or cannot be calculated, it is strongly recommended that the user estimate the missing climatic data with one of the procedures presented in Chapter 3 and that the FAO PenmanMonteith method be used for the calculation of ET_{o}. The use of an alternative ET_{o} calculation procedure, requiring only limited meteorological parameters, is less recommended.
Example 20 illustrates the estimation of monthly ET_{o} with the FAO PenmanMonteith for a data set containing only maximum and minimum air temperature. The procedures given in Chapter 3 to estimate missing humidity, radiation and wind speed data should be validated by comparing ET_{o} calculated with full and with limited data sets for weather stations in the region with complete data sets.
EXAMPLE 20. Determination of ET_{o} with missing data
Given the monthly average daily maximum and average daily minimum air temperature of July from a station near Lyon, France (45°43'N, altitude 200 m). No other climatic data were recorded. 

 
Monthly average daily maximum temperature (T_{max}) = 
26.6 
°C 
 
Monthly average daily minimum temperature (T_{min}) = 
14.8 
°C 
Determination according to Box 11 (calculation sheet ET_{o}) 

Estimation of wind speed: 

Parameters: 

 
T_{mean} = (26.6+14.8)/2 = 
20.7 
°C 
From Table 2.4 or 
T_{mean} = 
20.7 
°C 
Eq. 13: 
D = 
0.150 
kPa/°C 
From Table 2.2 or 
Altitude = 
200 
m 
Eq. 8: 
g = 
0.066 
kPa/°C 
 
(1 + 0.34 u_{2}) = (1 + 0.34 (2)) = 
1.68 
 
 
D /[D + g (1 + 0.34 u_{2})] = 0.150/[(0.150 + 0.066(1.68)] = 
0.576 
 
 
g /[D + g (1 + 0.34 u_{2})] = 0.0658/[0.150 + 0.066(1.68)] = 
0.252 
 
 
900/(T_{mean} + 273) u_{2} = 
6.13 
 
Estimation of humidity data: 

Assume (Eq. 48): 
T_{dew} » T_{min} = 
14.8 
°C 
Consequently(Eq. 14 or Table 2.3) for: 
T_{dew} = 
14.8 
°C 
Then 
e_{a} = 
1.68 
kPa 
From Table 2.3 or Eq. 11, for: 
T_{max} = 
26.6 
°C 
Then: 
e°(T_{max}) = 
3.48 
kPa 
From Table 2.3 or Eq. 11, for: 
T_{min} = 
14.8 
°C 
Then: 
e°(T_{min}) = 
1.68 
kPa 
 
e_{s} = (3.48+1.68)/2 = 
2.58 
kPa 
 
(e_{s}  e_{a}) = (2.581.68) = 
0.90 
kPa 
This corresponds with: 

 
RH_{max} = 100e_{a}/e°(T_{min}) = 
100 
% 
 
RH_{min} = 100 e_{a}/e°(T_{max}) = 100 (1.68/3.48) = 
48 
% 
 
RH_{mean} = (RH_{max} + RH_{min})/2 = 
74 
% 
Estimation of radiation data: 

R_{s} can be derived from the difference between T_{max} and T_{min}: 

From Eq. 50 
R_{s} = 0.16 Ö (26.614.8) R_{a} 
 
MJ m^{2} day^{1} 
 
R_{s} = 0.55 R_{a} 
 
MJ m^{2} day^{1} 
Table 2.6 or Eq. 21, for: 
For Day 15, Month = July, J = 
196 
 
Latitude = 45°43'N = 
45.72 
°N 

Then: 
R_{a} = 
40.55 
MJ m^{2} day^{1} 
 
R_{s} = 0.55 R_{a} = 0.55 (40.55) = 
22.29 
MJ m^{2} day^{1} 
 
R_{so} = (0.75 + 2 (200)/100000) 40.55 = 
30.58 
MJ m^{2} day^{1} 
 
R_{s}/R_{so} = 
0.73 
 
 
R_{ns} = 0.77 (22.29) = 
17.16 
MJ m^{2} day^{1} 
Table 2.8, for: 
T_{max} = 26.6°C 
26.6 
°C 
T_{max, K} = 26.6+273.16 = 
299.76 
K 

Then: 
_{} 
39.59 
MJ m^{2} day^{1} 
Table 2.8, for: 
T_{min} = 14.8°C 
14.8 
°C 
T_{min, K} = 14.8 +273.16 
287.96 
K 

Then: 
_{} 
33.71 
MJ m^{2} day^{1} 
 
_{} 
36.65 
MJ m^{2} day^{1} 
For: 
e_{a} = 
1.68 
kPa 
Then: 
(0.34  0.14Ö e_{a}) = 
0.16 
 
For: 
R_{s}/R_{so} = 
0.73 
 
Then: 
(1.35 R_{s}/R_{so}  0.35) = 
0.63 
 
 
R_{nl} = 36.65 (0.16) 0.63 = 
3.68 
MJ m^{2} day^{1} 
 
R_{n} = (17.163.68) = 
13.48 
MJ m^{2} day^{1} 
Assume: 
G = 
0 
MJ m^{2} day^{1} 
 
(R_{n}  G) = (13.480) = 
13.48 
MJ m^{2} day^{1} 
 
0.408 (R_{n}  G) = 
5.50 
mm/day 
Grass reference evapotranspiration: 

 
0.408 (R_{n}  G) D /[D + g (1+0.34 u_{2})] = 
3.17 
mm/day 
 
900/(T + 273) u_{2} (e_{s}  e_{a}) g /[D + g (1+0.34 u_{2})] = 
1.39 
mm/day 
 
ET_{o} = (3.17 + 1.39) = 
4.56 
mm/day 
The estimated grass reference evapotranspiration is 4.6 mm/day. If instead of 2 m/s, the wind speed is estimated as 1 or 3 m/s, ET_{o} would have been 7% lower (4.2 mm/day) or 6% higher (4.8 mm/day) respectively. In comparison, the Hargreaves equation (Equation 52) predicts ET_{o} = 5.0 mm/day 
The evaporation rate from pans filled with water is easily obtained. In the absence of rain, the amount of water evaporated during a period (mm/day) corresponds with the decrease in water depth in that period. Pans provide a measurement of the integrated effect of radiation, wind, temperature and humidity on the evaporation from an open water surface. Although the pan responds in a similar fashion to the same climatic factors affecting crop transpiration, several factors produce significant differences in loss of water from a water surface and from a cropped surface. Reflection of solar radiation from water in the shallow pan might be different from the assumed 23% for the grass reference surface. Storage of heat within the pan can be appreciable and may cause significant evaporation during the night while most crops transpire only during the daytime. There are also differences in turbulence, temperature and humidity of the air immediately above the respective surfaces. Heat transfer through the sides of the pan occurs and affects the energy balance.
Notwithstanding the difference between panevaporation and the evapotranspiration of cropped surfaces, the use of pans to predict ET_{o} for periods of 10 days or longer may be warranted. The pan evaporation is related to the reference evapotranspiration by an empirically derived pan coefficient:
ET_{o} = K_{p} E_{pan} (55)
where
ET_{o} reference evapotranspiration [mm/day],
K_{p} pan coefficient [],
E_{pan} pan evaporation [mm/day].
Pan types and environment
Different types of pans exist. Descriptions of Class A and Colorado sunken pans are given in Boxes 12 and 13. As the colour, size, and position of the pan have a significant influence on the measured results, the pan coefficients are pan specific.
In selecting the appropriate pan coefficient, not only the pan type, but also the ground cover in the station, its surroundings as well as the general wind and humidity conditions, should be checked. The siting of the pan and the pan environment also influence the results. This is particularly so where the pan is placed in fallow rather than cropped fields. Two cases are commonly considered: Case A where the pan is sited on a short green (grass) cover and surrounded by fallow soil; and Case B where the pan is sited on fallow soil and surrounded by a green crop (Figure 19).
FIGURE 19. Two cases of evaporation pan siting and their environment
Pan coefficients
Depending on the type of pan and the size and state of the upwind buffer zone (fetch), pan coefficients will differ. The larger the upwind buffer zone, the more the air moving over the pan will be in equilibrium with the buffer zone. At equilibrium with a large fetch, the air contains more water vapour and less heat in Case A than in Case B. Pan coefficients for the Class A pan and for the Colorado sunken pan for different ground cover, fetch and climatic conditions are presented in Tables 5 and 6. Regression equations derived from the tables are presented in Table 7. Where measured data from other types of sunken pans are available, such data should first be related to Colorado sunken pan data or to the FAO PenmanMonteith equation to develop K_{p}. Ratios between evaporation from sunken pans and from the Colorado sunken pan for different climatic conditions and pan environment are given in Table 8.
Where data are missing, wind speed can be estimated by taking a global value of u_{2} = 2 m s^{1} or as indicated in Table 4 (page 63). RH_{mean} can be approximated from air temperature as RH_{mean} = 50 e°(T_{min})/e°(T_{max}) + 50.
Adjustments
Under some conditions not accounted for in the tables, the presented K_{p} coefficients may need some adjustment. This is the case in areas with no agricultural development, or where the pans are enclosed by tall crops. Not maintaining the standard colour of the pan or installing screens can affect the pan readings and will require some adjustment of the pan coefficient.
In areas with no agricultural development and extensive areas of bare soils (large fetch, Case B), as found under desert or semidesert conditions, the listed values for K_{p} given for arid, windy areas may need to be reduced by up to 20%; for areas with moderate levels of wind, temperature and relative humidity, the listed values may need to be reduced by 510%; no or little reduction in K_{p} is needed in humid, cool conditions.
Where pans are placed in a small enclosure but surrounded by tall crops, for example 2.5 m high maize, the listed pan coefficients will need to be increased by up to 30% for dry windy climates whereas only a 510% increase is required for calm, humid conditions.
Painting the pans may affect the pan evaporation. The pan coefficients presented apply to galvanized pans annually painted with aluminium and to stainless steel pans. Little difference in E_{pan} will occur where the inside and outside surfaces of the pan are painted white. An increase in E_{pan} of up to 10% may occur when they are painted black. The material from which the pan is made may account for variations of only a few percent.
The level at which the water is maintained in the pan is important; resulting errors may be up to 15% when water levels in the Class A pan fall 10 cm below the accepted standard of between 5 and 7.5 cm below the rim. Screens mounted over pans will reduce E_{pan} by up to 10%. In an attempt to avoid pans being used by birds for drinking, pans filled to the rim with water can be placed near the Class A pan; birds may prefer to use the fully filled pan. The evaporation pan should be placed in a large, secure, wire enclosure to prevent animals from entering and drinking. The turbidity of the water in the pan usually does not affect E_{pan} by more than 5%. The overall variation in E_{pan} is not constant with time because of ageing, deterioration and repainting.
TABLE 5. Pan coefficients (K_{p}) for Class A pan for different pan siting and environment and different levels of mean relative humidity and wind speed (FAO Irrigation and Drainage Paper No. 24)
Class A pan 
Case A: Pan placed in short green cropped area 
Case B: Pan placed in dry fallow area 

RH mean (%) ® 

low < 40 
medium 40  70 
high > 70 

low < 40 
medium 40  70 
high > 70 
Wind speed (m s^{1}) 
Windward side distance of green crop (m) 



Windward side distance of dry fallow (m) 



Light 
1 
.55 
.65 
.75 
1 
.7 
.8 
.85 
< 2 
10 
.65 
.75 
.85 
10 
.6 
.7 
.8 
100 
.7 
.8 
.85 
100 
.55 
.65 
.75 

1000 
.75 
.85 
.85 
1000 
.5 
.6 
.7 

Moderate 
1 
.5 
.6 
.65 
1 
.65 
.75. 
.8 
25 
10 
.6 
.7 
.75 
10 
.55 
.65 
.7 
100 
.65 
.75 
.8 
100 
.5 
.6 
.65 

1000 
.7 
.8 
.8 
1000 
.45 
.55 
.6. 

Strong 
1 
.45 
.5 
.6 
1 
.6 
.65 
.7 
58 
10 
.55 
.6 
.65 
10 
.5 
.55 
.65 
100 
.6 
.65 
.7 
100 
.45 
.5 
.6 

1000 
.65 
.7 
.75 
1000 
.4 
.45 
.55 

Very strong 
1 
.4 
.45 
.5 
1 
.5 
.6 
.65 
> 8 
10 
.45 
.55 
.6 
10 
.45 
.5 
.55 
100 
.5 
.6 
.65 
100 
.4 
.45 
.5 

1000 
.55 
.6 
.65 
1000 
.35 
.4 
.45 
TABLE 6. Pan coefficients (K_{p}) for Colorado sunken pan for different pan siting and environment and different levels of mean relative humidity and wind speed (FAO Irrigation and Drainage Paper No. 24)
Sunken Colorado 
Case A: Pan placed in short green cropped area 
Case B: Pan placed in dry fallow area (1) 

RH mean (%) ® 

low < 40 
medium 40  70 
high > 70 

low < 40 
medium 40  70 
high > 70 
Wind speed (m s^{1}) 
Windward side distance of green crop (m) 



Windward side distance of dry fallow (m) 



Light 
1 
.75 
.75 
.8 
1 
1.1 
1.1 
1.1 
< 2 
10 
1.0 
1.0 
1.0 
10 
.85 
.85 
.85 
³ 100 
1.1 
1.1 
1.1 
100 
.75 
.75 
.8 





1000 
.7 
.7 
.75 

Moderate 
1 
.65 
.7 
.7 
1 
.95 
.95 
.95 
25 
10 
.85 
.85 
.9 
10 
.75 
.75 
.75 
³ 100 
.95 
.95 
.95 
100 
.65 
.65 
.7 





1000 
.6 
.6 
.65 

Strong 
1 
.55 
.6 
.65 
1 
.8 
.8 
.8 
58 
10 
.75 
.75 
.75 
10 
.65 
.65 
.65 
³ 100 
.8 
.8 
.8 
100 
.55 
.6 
.65 





1000 
.5 
.55 
.6 

Very strong 
1 
.5 
.55 
.6 
1 
.7 
.75 
.75 
> 8 
10 
.65 
.7 
.7 
10 
,55 
.6 
.65 
³ 100 
.7 
.75 
.75 
100 
.5 
.55 
.6 





1000 
.45 
.5 
.55 
(1) For extensive areas of barefallow soils and no agricultural development, reduce K_{pan} by 20% under hot, windy conditions; by 510% for moderate wind, temperature and humidity conditions.
TABLE 7. Pan coefficients (K_{p}): regression equations derived from Tables 5 and 6
Class A pan with green fetch 
K_{p} = 0.108  0.0286 u_{2} + 0.0422 ln(FET) +0.1434 ln(RH_{mean})  0.000631 [ln(FET)]^{2} ln(RH_{mean}) 
Class A pan with dry fetch 
K_{p} = 0.61 + 0.00341 RH_{mean}  0.000162 u_{2} RH_{mean}  0.00000959 u_{2} FET + 0.00327 u_{2} ln(FET)  0.00289 u_{2} ln(86.4 u_{2})  0.0106 ln(86.4 u_{2})ln(FET) + 0.00063 [ln(FET)]^{2}ln(86.4 u_{2}) 
Colorado sunken pan with green 
K_{p} = 0.87 + 0.119 ln(FET)0.0157[ln(86.4 u_{2})]^{2} In(RH_{mean})  0.000053 ln(86.4u_{2})ln(FET)RH_{mean} 
Colorado sunken pan with dry fetch 
K_{p} = 1.145  0.080 u_{2} + 0.000903(u_{2})^{2}ln(RH_{mean})  0.0964 ln(FET) + 0.0031 u_{2} In(FET) + 0.0015[ln(FET)]^{2}ln(RH_{mean}) 
Coefficients and parameters 
K_{p} pan coefficient [] 
Range for variables 
1 m £ FET £ 1000 m (these limits must be observed) 
Recommendations
The above considerations and adjustments indicate that the use of tables or the corresponding equations may not be sufficient to consider all local environmental factors influencing K_{p} and that local adjustment may be required. To do so, an appropriate calibration of E_{pan} against ET_{o} computed with the PenmanMonteith method is recommended.
It is recommended that the pan should be installed inside a short green cropped area with a size of a square of at least 15 by 15 m. The pan should not be installed in the centre but at a distance of at least 10 m from the green crop edge in the general upwind direction.
Where observations of wind speed and relative humidity, required for the computation of K_{p}, are not available at the site, estimates of the weather variables from a nearby station have to be utilized. It is then recommended that these variables be averaged for the computation period and that E_{pan} be averaged for the same period.
Equation 1 in Table 7 yields K_{p} = 0.83 for data in Example 21 as shown in Example 22.
TABLE 8. Ratios between the evaporation from sunken pans and a Colorado sunken pan for different climatic conditions and environments (FAO Irrigation and Drainage Paper No. 24)

Ratio E_{pan} mentioned and E_{pan} Colorado  
Climate 
Humidtemperate climate 
Arid to semiarid (dry season)  
Ground cover surrounding pan (50 m or more) 
Short green cover 
Dry fallow 
Short green cover 
Dry fallow  

Pan area (m^{2}) 



 
GGI 20 diameter 5 m, depth 2 m (former Soviet Union) 
19.6 
1.0 
1.1 
1.05 
1.25  
Sunken pan diameter 12 ft (3.66 m) depth 3.3 ft (Israel) 
10.5 












BPI diameter 6 ft (1.83 m), depth 2 ft (0.61 m) (USA) 
2.6 












Kenya pan diameter 4 ft (1.22 m) depth 14 in (0.356 m) 
1.2 












Australian pan diameter 3 ft (0.91 m) depth 3 ft (0.91 m) 
0.7 



1.0 



1.0  
Symmons pan 6 ft^{2} (0.56 m^{2}) depth 2 ft (0.61 m) 
0.6 












Aslyng pan 0.33 m^{2}, depth 1 m (Denmark) 
0.3 






1.0 


 
GGI 3000 diameter 0.618 cm, depth 6080 cm (former Soviet Union) 
0.3 












Sunken pan diameter 50 cm, depth 25 cm (Netherlands) 
0.2 
1.0 
0.95 
1.0 
0.95 
EXAMPLE 21. Determination of ET_{o} from pan evaporation using tables
Given the daily evaporation data for the first week of July for a Class A pan installed in a green area surrounded by short irrigated field crops: 8.2, 7.5, 7.6, 6.8, 7.6, 8.9 and 8.5 mm/day. In that period the mean wind speed is 1.9 m/s and the daily mean relative humidity is 73%. Determine the 7day average reference evapotranspiration.  
Pan is installed on a green surface: Case A  
Pan is surrounded by irrigated crops: 
fetch_{max} = 
1000 
m 
Wind speed is light: 
u < 
2 
m/s 
Relative humidity is high: 
RH_{mean} > 
70 
% 
From Table 5 (for above conditions): 
K_{p} = 
0.85 
 
 
E_{pan} = (8.2 + 7.5 + 7.6 + 6.8 + 7.6 + 8.9 + 8.5)/7 = 
7.9 
mm/day 
From Eq. 55: 
ET_{o} = 0.85 (7.9) = 
6.7 
mm/day 
The 7day average of the crop reference evapotranspiration is 6.7 mm/day 
BOX 12. Description of Class A pan The Class A Evaporation pan is circular, 120.7 cm in diameter and 25 cm deep. It is made of galvanized iron (22 gauge) or Monel metal (0.8 mm). The pan is mounted on a wooden open frame platform which is 15 cm above ground level. The soil is built up to within 5 cm of the bottom of the pan. The pan must be level. It is filled with water to 5 cm below the rim, and the water level should not be allowed to drop to more than 7.5 cm below the rim. The water should be regularly renewed, at least weekly, to eliminate extreme turbidity. The pan, if galvanized, is painted annually with aluminium paint. Screens over the pan are not a standard requirement and should preferably not be used. Pans should be protected by fences to keep animals from drinking. The site should preferably be under grass, 20 by 20 m, open on all sides to permit free circulation of the air. It is preferable that stations be located in the centre or on the leeward side of large cropped fields. Pan readings are taken daily in the early morning at the same time that precipitation is measured. Measurements are made in a stilling well that is situated in the pan near one edge. The stilling well is a metal cylinder of about 10 cm in diameter and some 20 cm deep with a small hole at the bottom.

BOX 13. Description of Colorado sunken pan The Colorado sunken pan is 92 cm (3 ft) square and 46 cm (18 in) deep, made of 3 mm thick iron, placed in the ground with the rim 5 cm (2 in) above the soil level. Also, the dimensions 1 m square and 0.5 m deep are frequently used. The pan is painted with black tar paint. The water level is maintained at or slightly below ground level, i.e., 57.5 cm below the rim. Measurements are taken similarly to those for the Class A pan. Siting and environment requirements are also similar to those for the Class A pan. Sunken Colorado pans are sometimes preferred in crop water requirements studies, as these pans give a better direct estimation of the reference evapotranspiration than does the Class A pan. The disadvantage is that maintenance is more difficult and leaks are not visible.

EXAMPLE 22. Determination of ET_{o} from pan evaporation using equations
Given the 7day average evaporation measurement from Example 21, estimate the ET_{o} for the two types of pans and two types of fetch conditions represented by equations in Table 7. Assume that fetch is 1000 m in both fetch cases (green and dry). 


F_{etch} = 

1000 
m 
u_{2} = 
1.9 
m/s 

RH_{mean} = 
73 
% 

Class A pan with green fetch 
K_{p} = 0.108  0.0286 u_{2} + 0.0422 ln(FET) + 0.1434 ln(RH_{mean})  0.000631 [ln(FET)]2 ln(RH_{mean}) 
K_{p} = 0.108  0.0286 (1.9) + 0.0422 ln(1000) +0.1434 ln(73)  0.000631 [ln(1000)]2 ln(73) 
0.83 
 
E_{pan} = 

7.9 
mm/day 

ET_{o} = 
ET_{o} = 0.83 (7.9) 
6.6 
mm/day 

Class A pan with dry fetch 
K_{p} = 0.61 + 0.00341 RH_{mean}  0.000162 u_{2} RH_{mean}  0.00000959 u_{2} FET + 0.00327 u_{2} ln(FET)  0.00289 u_{2} ln(86.4 u_{2})  0.0106 ln(86.4u_{2})ln(FET) + 0.00063 [ln(FET)]^{2}ln(86.4 u_{2}) 
K_{p} = 0.61 + 0.00341 (73)  0.000162 (1.9) (73)  0.00000959 (1.9)(1000) + 0.00327 (1.9) ln(1000)  0.00289 (1.9) ln(86.4(1.9))  0.0106 ln(86.4(1.9))ln(1000) + 0.00063 [ln(1000)]^{2} ln(86.4 (1.9)) 
0.61 
 
E_{pan} = 

7.9 
mm/day 

ET_{o} = 
ET_{o} = 0.61 (7.9) 
4.8 
mm/day 

Colorado sunken pan with green fetch 
K_{p} = 0.87 + 0.119 ln(FET)  0.0157[ln(86.4 u_{2})]^{2}  0.0019 [ln(FET)]^{2}ln(86.4 u_{2}) + 0.013 ln(86.4 u_{2}) ln(RH_{mean})  0.000053 ln(86.4 u_{2})ln(FET) RH_{mean} 
K_{p} = 0.87 + 0.119 ln(100)  0.0157 [ln(86.4(1.9))]^{2}  0.0019 [ln(1000)]2 ln(86.4 (1.9)) + 0.013 ln(86.4(1.9)) ln(73)  0.000053 ln(86.4 (1.9)) ln(1000) (73) 
0.97 
 
E_{pan} = 

7.9 
mm/day 

ET_{o} = 
ET_{o} = 0.97(7.9) 
7.7 
mm/day 

Colorado sunken pan with dry fetch 
K_{p} = 1.145  0.080 u_{2} +0.000903(u_{2})^{2}ln(RH_{mean})  0.0964 ln(FET) + 0.0031 u_{2} ln(FET) + 0.0015 [ln(FET)]^{2} ln(RH_{mean}) 
K_{p} = 1.145  0.080(1.9) + 0.000903(1.9)^{2}ln(73)  0.0964 ln(1000) + 0.0031 (1.9) ln(1000) +0.0015 [ln(1000)]^{2}ln(73) 
0.69 
 
E_{pan} = 

7.9 
mm/day 

ET_{o} = 
ET_{o} = 0.69 (7.9) 
5.4 
mm/day 

The 7day average of the crop reference evapotranspiration for the four pan/fetch conditions is 6.6, 4.8, 7.7, and 5.4 mm/day 