poster - International Conference of Agricultural Engineering

poster - International Conference of Agricultural Engineering poster - International Conference of Agricultural Engineering

cigr.ageng2012.org
from cigr.ageng2012.org More from this publisher
10.11.2014 Views

Monteith equation to climatic and parametric factors in a semi-arid climate for a reference grass surface, grain sorghum and sweet sorghum in Italy. For grass, available energy and aerodynamic resistance played a major role. For sweet sorghum, the model was most sensitive to vapor pressure deficit. For grain sorghum under water stress, the most sensitive term was canopy resistance. Irmak et al. (2006) calculated the sensitivity coefficient of the standardized daily ASCE-Penman-Monteith equation in different climates of United States. Recently, Ali et al. (2009) calculated the sensitivity coefficient of the FAO Penman-Monteith equation under the environment of a semi-humid sub-tropic region of Bangladesh. The results showed that the ET o estimates are most sensitive to maximum temperature, relative humidity, sunshine duration, wind speed and minimum temperature, respectively. Estevez et al. (2009) calculated the sensitivity coefficients of the standardized ASCE-Penman-Monteith equation from 87 automatic weather stations in Spain. The results showed a large degree of daily and seasonal variability, especially for temperature and relative humidity. Also, the sensitivity of ET o to the same climatic variables showed significant differences among locations. Sensitivity of the daily standardized ASCE equation to meteorological variables has not yet been studied in Iran. Thus, the objective of this study was to investigate the sensitivity of the standardized daily ASCE ET o equation to climatic variables under different subclasses of arid and semiarid climatic conditions in Kerman province, southeast of Iran and to derive sensitivity coefficients for each one. 2. Materials and Methods 2.1. Study area and weather data source The study was carried out in Kerman province (Southeast of Iran), located between the latitudes 27° and 30° N and the longitudes 55° and 58° W. Site elevations range from 400 to 2775 m above mean sea level. All agricultural productions in this area are irrigation-based. In this region irrigation water resources are supplied mostly from groundwater and slightly spring and Qanat. Surface irrigation is the most popular method of irrigation in this area, however, frequent droughts have led to switch over to pressurized irrigation systems to improve water use efficiency and prevent depletion of groundwater resources. Meteorological data used in this analysis were obtained over 4 stations in Kerman province from the Kerman Meteorological Department. These stations are located in different subclasses of arid and semiarid climates of the province based on extended-De Martonne classification (Khalili, 1997). Daily weather data were used in this study. Table 1 lists the annual average weather data of 4 meteorological stations in Kerman province, including the corresponding site elevations and coordinates. Average annual rainfall ranged from 61.3 mm at Bam to 261.6 mm at Baft; the range for the annual air temperature was 14.8-25.0 °C; 31- 43 %, for relative humidity; and 0.5-1.6 m s-1 for wind speed at 2 m. TABLE 1: Summary of weather station sites characteristics used in the study 2.2. ASCE standardized reference evapotranspiration equation The purpose of the standardized reference ET equation and calculation procedures is to bring commonality to the calculation of reference ET and to provide a standardized basis for determining or transferring crop coefficients for agricultural and landscape use (ASCE- EWRI, 2005). For the standardization, the ASCE-Penman-Monteith method was applied for 3

two types of reference surfaces representing clipped grass and alfalfa, and for daily or hourly time step. This equation is simplified to a reduced form of the ASCE–PM. The equation is expressed as: 900 0.408 ( Rn G) U 2 ( es ea ) ET T 273 o (1) (1 0.34U 2 ) where ET o is standardized grass reference evapotranspiration (mm day-1), R n is calculated net radiation at the crop surface (MJ m -2 d -1 ), G is soil heat flux density at the soil surface (MJ m -2 d -1 ), T is mean daily air temperature at 2 m height (°C), U 2 is mean daily wind speed at 2 m height (m s-1), e s is saturation vapor pressure (kPa), calculated for daily time steps as the average of saturation vapor pressure at maximum and minimum air temperature, ea is mean actual vapor pressure (kPa), Δ is slope of the saturation vapor pressure-temperature curve (kPa °C -1 ), γ is psychrometric constant (kPa °C -1 ), Units for the 0.408 coefficient are m 2 mm MJ -1 . All the calculations of daily values (R n , e s , e a ) and other parameters were made by ASCE-EWRI (2005). 2.3. Sensitivity analyses and sensitivity coefficients In ecological and hydrometeorological studies (e.g., McCuen, 1974; Saxton, 1975; Beven, 1979; Anderton et al., 2002) a number of sensitivity coefficients have been defined depending on the purpose of the analyses. For example, Saxton (1975) mathematically differentiated the equation under investigation to derive equations for the rate of change of the independent variable with respect to each dependent variable. Smajstrla et al. (1987) defined the sensitivity coefficient as the slope of the curve of ET o versus the climatic variable being studied. Slopes computed in this manner represented the rates of change in ET o with respect to change in the climatic variable. For multi-variable models (e.g., the Penman- Monteith equation), different variables have different dimensions and different ranges of values, which makes it difficult. In general, a model of evapotranspiration can be written as: ET f p p , p ,..., (2) o 1, 2 3 p N where i p is input data variable and N is the number of parameters and input data variables. The error in ET o ( ETo ) that results from errors in the pi can then be expressed as: ETo ETo f p1 p1, p2 p2 ,..., pN p N (3) Expanding equation 3 in Taylor series and ignoring second-order terms and above, leads to: ETo ETo ETo ETo p1 p2 ... pN (4) p p p 1 ET where the differential p 2 i o N is the absolute sensitivity of the estimation of p i , and pi is the individual error associated with p i (Beven, 1979). Following McCuen (1974) and Beven (1979), the partial derivative is transformed into a nondimensional form to display the sensitivity for the variables: ETo / ETo ETo pi S lim / P i pi pi pi ETo (5) Δp 0 i where S P is the sensitivity coefficient represents the fraction of the change in p i i that is transmitted through to the estimate of ET o . However, the relative coefficients are sensitive to the values of p and ET o . p in this study represents air temperature, solar radiation, relative i i humidity and wind speed. Basically, a positive/negative sensitivity coefficient of a variable 4

two types <strong>of</strong> reference surfaces representing clipped grass and alfalfa, and for daily or hourly<br />

time step. This equation is simplified to a reduced form <strong>of</strong> the ASCE–PM. The equation is<br />

expressed as:<br />

900<br />

0.408 (<br />

Rn<br />

G)<br />

U<br />

2<br />

( es<br />

ea<br />

)<br />

ET<br />

T 273<br />

o<br />

<br />

<br />

(1)<br />

(1 0.34U<br />

2<br />

)<br />

where ET o is standardized grass reference evapotranspiration (mm day-1), R n is calculated<br />

net radiation at the crop surface (MJ m -2 d -1 ), G is soil heat flux density at the soil surface<br />

(MJ m -2 d -1 ), T is mean daily air temperature at 2 m height (°C), U 2 is mean daily wind speed<br />

at 2 m height (m s-1), e s is saturation vapor pressure (kPa), calculated for daily time steps as<br />

the average <strong>of</strong> saturation vapor pressure at maximum and minimum air temperature, ea is<br />

mean actual vapor pressure (kPa), Δ is slope <strong>of</strong> the saturation vapor pressure-temperature<br />

curve (kPa °C -1 ), γ is psychrometric constant (kPa °C -1 ), Units for the 0.408 coefficient are m 2<br />

mm MJ -1 . All the calculations <strong>of</strong> daily values (R n , e s , e a ) and other parameters were made by<br />

ASCE-EWRI (2005).<br />

2.3. Sensitivity analyses and sensitivity coefficients<br />

In ecological and hydrometeorological studies (e.g., McCuen, 1974; Saxton, 1975; Beven,<br />

1979; Anderton et al., 2002) a number <strong>of</strong> sensitivity coefficients have been defined<br />

depending on the purpose <strong>of</strong> the analyses. For example, Saxton (1975) mathematically<br />

differentiated the equation under investigation to derive equations for the rate <strong>of</strong> change <strong>of</strong><br />

the independent variable with respect to each dependent variable. Smajstrla et al. (1987)<br />

defined the sensitivity coefficient as the slope <strong>of</strong> the curve <strong>of</strong> ET o versus the climatic variable<br />

being studied. Slopes computed in this manner represented the rates <strong>of</strong> change in ET o with<br />

respect to change in the climatic variable. For multi-variable models (e.g., the Penman-<br />

Monteith equation), different variables have different dimensions and different ranges <strong>of</strong><br />

values, which makes it difficult. In general, a model <strong>of</strong> evapotranspiration can be written as:<br />

ET f p p , p ,...,<br />

(2)<br />

<br />

o 1, 2 3<br />

p N<br />

where<br />

i<br />

<br />

p is input data variable and N is the number <strong>of</strong> parameters and input data variables.<br />

The error in ET o ( ETo<br />

) that results from errors in the pi<br />

can then be expressed as:<br />

ETo<br />

ETo<br />

f p1 p1, p2<br />

p2<br />

,..., pN<br />

p<br />

N<br />

<br />

(3)<br />

Expanding equation 3 in Taylor series and ignoring second-order terms and above, leads to:<br />

ETo<br />

ETo<br />

ETo<br />

ETo<br />

p1<br />

p2<br />

... pN<br />

(4)<br />

p<br />

p<br />

p<br />

1<br />

ET<br />

where the differential<br />

p<br />

2<br />

i<br />

o<br />

N<br />

is the absolute sensitivity <strong>of</strong> the estimation <strong>of</strong><br />

p<br />

i<br />

, and<br />

pi<br />

is the<br />

individual error associated with p<br />

i (Beven, 1979).<br />

Following McCuen (1974) and Beven (1979), the partial derivative is transformed into a nondimensional<br />

form to display the sensitivity for the variables:<br />

ETo<br />

/ ETo<br />

ETo<br />

pi<br />

S lim<br />

<br />

<br />

/<br />

<br />

P i<br />

pi<br />

pi<br />

pi<br />

ETo<br />

(5)<br />

Δp 0<br />

i<br />

<br />

where S<br />

P is the sensitivity coefficient represents the fraction <strong>of</strong> the change in p<br />

i<br />

i that is<br />

transmitted through to the estimate <strong>of</strong> ET o . However, the relative coefficients are sensitive to<br />

the values <strong>of</strong> p and ET o . p in this study represents air temperature, solar radiation, relative<br />

i<br />

i<br />

humidity and wind speed. Basically, a positive/negative sensitivity coefficient <strong>of</strong> a variable<br />

4

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!