Evaluation of Crop Water Stress Index, Canopy Temperature and Grain Yield of Five Iranian Wheat Cultivars Under Late Season Drought Stress

Document Type : Research Paper

Authors

1 Shiraz University, Shiraz, Iran

2 Department of Plant Production and Genetics, Shiraz University, Shiraz, Iran

Abstract

Abstract
In order to evaluate crop water stress index (CWSI) and canopy temperature of wheat cultivars under terminal drought stress, a field experiment was conducted at the Agricultural Research Station of Shiraz University, Shiraz, during 2009 growing season. Five wheat cultivars including Shiraz, Bahar, Pishtaz, Sistan and Yavaros and four levels of water regime including well watering [Irrigation according to 100% field capacity (FC)], excess watering (125% FC), and mild (75% FC) and severe drought (50% FC) stress were used in a split plot design experiment with three replicates. Results showed that Yavaros and Shiraz cultivars with 7.36 and 6.81°C had the highest canopy-air temperature differences (Tc-Ta), respectively, while in Bahar this difference was 3.9°C. In all cultivars, slope (a) and intercept (b) of lower base line equation between Tc-Ta and vapour pressure deficit (VPD) were increased significantly due to more limitation in water and increasing VPD. Yavaros and Shiraz cultivars with higher a value were found to be more sensitive to increasing VPD. Shiraz and Yavaros cultivars with 0.73 and 0.71 had the highest seasonal mean CWSI, respectively, while CWSI in Bahar, Pishtaz and Sistan ranged from 0.61 to 0.64 under severe drought. A negative relationship was found between CWSI and amount of water supply and net photosynthesis of flag leaf. Maximum grain yield was obtained in Shiraz and Yavaros under well and excess watering and CWSI in these cultivars ranged from 0.31 to 0.36, whereas by decreasing water supply and increasing CWSI, grain yield in these cultivars decreased significantly. Bahar, Pishtaz and Sistan cultivars with lower Tc-Ta, water supply and CWSI had better performance than Shiraz and Yavaros cultivars, especially when exposed to water stress conditions. The role of these traits should be further investigated as potential indirect selection criteria for grain yield of wheat cultivars in semi-arid conditions.
 
 

Keywords


Introduction

Wheat is an important cereal crop and is adapted to a wide range of climatic conditions (Emam 2007). The success of sustained wheat production in arid and semi-arid regions of the world depends entirely on water availability (Alderfarsi and Nielsen 2001). Efficient use of water in the Mediterranean region is becoming an important issue due to increasing irrigation water requirements as well as environmental sustainability (Emekli et al. 2007).

Canopy temperature is a part of the canopy energy balance. As solar radiation is absorbed by leaves, leaf temperatures increase (Panda et al. 2003). Leaf cooling takes place as some of the thermal energy drives transpirational water loss. Under water deficit conditions, stomata close in response to loss of turgor pressure, causing a lowering of transpiration rate and hence, an increase in canopy temperature (Kramer 1983). This is the basis for the use of canopy temperature to determine plant water status (Jackson et al. 1981).

The crop water stress index (CWSI) calculation is based on three main environmental variables: plant canopy temperature (Tc), air temperature (Ta) and atmospheric vapor pressure deficiency (VPD). All these three variables have much influence on water used by plants (Braunworth 1989). An infrared thermometer measures the surface temperature of a crop canopy without making direct physical contact (Howell et al. 1986).

Idso et al. (1981) defined CWSI based on the empirical linear relationship between midday Tc-Ta and VPD under high net radiation and well watered conditions. The CWSI has been used to quantify water status in the field based on canopy temperature (Yuan et al. 2004; Emekli et al. 2007) and irrigation scheduling of wheat in many places (Alves and Pereira 2000; Alderfarsi and Nielsen 2001; Orta et al. 2004). Feng et al. (2009) asserted that wheat cultivars with low canopy temperature could maintain superiority to cultivars with high canopy temperature and low canopy temperature in wheat could be used as an index to evaluate physiological capacities of wheat under drought stress and also as a useful marker in wheat breeding for drought tolerance. Al-Faraj et al. (2001) reported that Tc-Ta was increased with a decrease in soil water content for tall fescue (Festuca arundinacea Schreb.). They suggested that CWSI could be used for irrigation timing in turfgrass.  Furthermore, Jalali-Farahani et al. (1993) showed that changes in CWSI depended on the applied irrigation volume.

Little research has been done to quantify the CWSI of Iranian wheat cultivars especially in south of the country, where water stress in wheat is pervasive and frequent during grain filling period. The objectives of the present study were: (I) to develop a baseline equation which can be used to calculate CWSI for monitoring water status of Iranian wheat cultivars and (II) to evaluate the relationship of CWSI with amount of water supply, net photosynthesis rate and yield of wheat cultivars under late season drought stress.

 

 

Materials and Methods

A field experiment was conducted during November 2009 - June 2010 at the Agricultural Research Station of Shiraz University, Shiraz, Iran, for establishing the crop water stress index of wheat crop. The research station is located at a latitude of 29°44' N, a longitude of 52°37' E and an altitude of 1810 m. Ten-day averages of some meteorological data measured daily in the study area during April to June 2010  are shown in Table 1. The research area has Mediterranean climate with hot and dry summers and cool and rainy winters. Five wheat cultivars including Shiraz, Bahar, Pishtaz and Sistan (as bread wheat) and Yavaros (as durum wheat) were arranged in sub-plots and four levels of water regime including well watering [Irrigation according to 100% field capacity (FC)], excess watering (125% FC) and mild (75% FC) and severe drought (50% FC) stress in main plots of a split plot experiment with three replicates.

On November 7th 2009, commercial wheat seeds were sown in rows 30 cm apart with 250 plants/m2 in plots of 2×5 m. Before planting, 60 kg P/ha, as super phosphate, and 60 kg N/ha, as urea, were applied. Another 60 kg N/ha was added at the end of tillering stage. The soil water status was monitored in each plot by gravimetric method at 30 cm intervals down to 120 cm and irrigation regimes were applied at booting stage until late season (i.e. physiological ripening). The amount of water supply was measured by time-volume technique according to Grimes et al. (1987) and is presented in Table 2 for each cultivar under different irrigation regimes.

 

 

 

 

 

 

 

Table 1. Ten-day means of climatic data measured daily at experimental site

 

Wind speed (m/s)

Relative humidity (%)

Mean evaporation (mm)

Temperature (°C)

Month

 

 

 

 

April

1.15

57.2

6.1

12.5

1-10

1.23

65.7

4.4

13.3

11-20

1.16

56.7

5.6

15.2

21-30

 

 

 

 

May

1.24

58.8

5.3

15.7

1-10

1.08

56.6

7.7

19.6

11-20

0.84

49.4

9.8

20.7

21-31

 

 

 

 

June

0.82

46.9

8.9

22.5

1-10

0.68

38.5

19.5

24.1

11-20

1.02

36.1

11.6

23.8

21-30

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Table 2. Total water applied (mm) in each irrigation regime and wheat cultivar

 

Irrigation regime

(according to field capacity)

Wheat cultivar

Water applied (mm)

Severe drought (50% FC)

Shiraz

302.7

 

Bahar

296.0

 

Pishtaz

292.7

 

Sistan

283.2

 

Yavaros

 

350.2

Mild drought (75% FC)

Shiraz

454.1

 

Bahar

444.0

 

Pishtaz

439.1

 

Sistan

424.7

 

Yavaros

 

525.2

Well watered (100% FC)

Shiraz

605.4

 

Bahar

592.0

 

Pishtaz

585.4

 

Sistan

566.3

 

Yavaros

 

700.3

Excess watered (125% FC)

Shiraz

756.8

 

Bahar

740.0

 

Pishtaz

731.8

 

Sistan

707.9

 

Yavaros

875.4

               FC: Field Capacity

 

 

 

 

 

 

 

 

To determine CWSI of wheat cultivars, an infrared thermometer (Kyoritsu Electronic Instrument, Model 5500, Japan) was used and the canopy temperature was measured (3, 6 and 9 days after each irrigation) from 4 April to 21 June 2010 (151- 233 days after planting). To ensure collection of accurate data, the infrared thermometer was held with a horizontal angle of 45° during measurements. Temperature measurement was done when there was no cloud. According to Idso et al. (1981), midday canopy temperature is the best indicator to detect the crop water stress. The measurements were carried out from four directions (east, west, north and south) in each experimental plot.

Simultaneously, air temperature and relative humidity were recorded using thermohygrograph (Lambrecht, Model 252, Germany) and psychrometer (Lambrecht, Model 1030, Germany) as basis for calculating vapour pressure deficit (VPD) (Monteith and Unsworth 1990).  VPD was computed from standard psychrometer equation (Allen et al. 1998). Then, CWSI values were alculated using the empirical method of Idso et al. (1981). The relationship between canopy-air temperature differences (Tc-Ta) and VPD were plotted under stressed and non- stressed conditions (Figure 1). In this graph, the non-stressed baseline for each wheat cultivar was determined from the data collected three days after irrigation in excess watered treatment between 08:00 and 17:00 h with 30-min intervals. The Idso’s empirical non-water-stressed baseline can be expressed as Equation (1):

Tc Ta = aVPD +b                                      (1)

where Tc–Ta is the measured canopy and air temperature differences for non- stressed treatment (°C) and VPD is vapour pressure deficit (kPa) and a (slope) and b (intercept) are the linear regression coefficients of TcTa on VPD. The upper baseline was determined using the average Tc–Ta values measured at 13:00, 14:00 and 15:00 h before each irrigation. Using the upper and lower limit estimates, a CWSI can be defined by the following Equation (2) (Idso et al. 1981):

CWSI=[(Tc-Ta)m-(Tc-Ta)ll]/[(Tc-Ta)ul-(Tc-Ta)ll]                                                   

(2)           

where (Tc–Ta)m, (Tc–Ta)ll and (Tc–Ta)ul are the measured canopy and air temperature differences at the moment and the lower and upper limit values (°C), respectively.

 Furthermore, net photosynthesis rate (Pn) was measured from the flag leaf of main shoot of each plant using a portable photosynthesis system (IGRA model LCA4-ADC, Hoddeson, UK) after flowering to maturity. Grain yield measured from centre of 1 m2 final harvest area in each plot. The data were analyzed using SAS (2003) software and means were compared by Duncan’s multiple range test at 0.05 probability level.

 

Results and Discussion

Canopy temperature

Variation in Tc-Ta observed among wheat cultivars, was significant at 5% probability level. Yavaros and Shiraz cultivars with 7.36 and 6.81 °C had the higher canopy temperature differences, while in Bahar this difference reached to 3.9 °C (Table 3 and Figure 1). In a similar study, Feng et al. (2009) concluded that canopy temperature could be considered as a consistent character for each wheat genotype. They declared the difference in canopy temperature between low temperature wheat cultivars and high temperature cultivars could be observed mainly during the grain filling period, a key period for wheat to form grain. Xiaoyan 6 could be considered as a low canopy temperature wheat genotype (LTW), whereas Yanshi 9, NR 9405, 9430 as high canopy temperature wheat genotypes (HTW). In the present study, differences in LTW among Bahar, Pishtaz and Sistan and in HTW between Shiraz and Yavaros were significant (Table 3).

 

 

 

 

Table 3. Comparison of the upper limits values of canopy and air temperature difference (Tc -Ta)ul  and slopes and intercepts for lower limit [(Tc -Ta)ll = a VPD + b] of five wheat cultivars

 

 

 

Wheat cultivars

 

 

 

Shiraz

Bahar

Pishtaz

Sistan

Yavaros

Tc-Ta

6.81a

3.95b

4.12b

4.43b

7.36a

Slope (a)

-1.41a

-0.96bc

-0.93c

1.0b

-1.45a

Intercept (b)

2.89a

1.40b

1.03c

0.89c

3.01a

(Tc - Ta)ul = canopy-air temperature difference at the lower limit values (no stress); (Tc -Ta)ll = canopy-air temperature difference at the upper limit values (stress). Means in each row by the same letters are not significantly different at 5%

 

probability level using Duncan’s multiple range test.

 

 

 

 

 

Determination of lower base line

For all wheat cultivars slope (a) and intercept (b) were compared in Table 3. In all cultivars, a and b of lower base line equation between Tc-Ta and VPD were significantly increased due to more limitation in water and increasing VPD (Figure 1). Our result was in agreement with Orta et al. (2004) who declared that Tc-Ta measured above a crop was negatively related to the atmospheric VPD (R2=0.88). The value of a varied from -1.45 in Yavaros to -0.93 in Pishtaz (Table3). It appeared that Yavaros and Shiraz cultivars with higher a value were more sensitive to increasing VPD (Table 3 and Figure 1). On the other hand, in these cultivars difference between upper base line (under stress) and lower base line (non-stress) was more than Bahar, Pishtaz and Sistan cultivars (Figure 1). Also, the value of b ranged from 3.01 to 0.89 and was significantly different among wheat cultivars (Table3) In India, Gontia and Tiwari (2008) showed that the lower baseline equation obtained for wheat crop was (Tc-Ta)ll= -1.1141(VPD)+1.0827 during flowering to maturity and a and b parameters in the following equation were close to parameter of Bahar, Pishtaz and Sistan (Figure 1). On the other hand, in Shiraz and Yavaros cultivars the value of a was very close to that reported by Alderfarsi

 

    

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Figure 1.  Stressed and non-stressed baselines for calculation of CWSI in five wheat cultivars.

VPD = vapour pressure deficit

 

 

and Nielsen (2001) for winter wheat in Colorado [(Tc-Ta)ll= -1.5VPD+0.41], however, b in this equation was smaller than that for Shiraz and Yavaros cultivars of our study. This might be attributed to higher temperature in our experimental site, i.e. Shiraz, from April to June (Table 1), compared to Colorado.  Indeed, many researchers pointed out that type of cultivar and local conditions could influence the baseline equation causing differences in slopes and intercepts (Alves and Pereira 2000; Panda et al. 2003; Yuan et al. 2004).

 

CWSI variation

In all cultivars and irrigation regimes, high amount of variation (0.19 to 0.83) was observed from April to June for monthly CWSI and increased by progressing drought stress late in the season (Table 4). Garrot et al. (1994) reported that in durum wheat (CV. Aldura) mean CWSI varied from 0.11 under well watered to 0.82 under severe drought conditions. Type of wheat cultivar and amount of water applied had noticeable effect on seasonal mean CWSI. Shiraz and Yavaros cultivars with 0.73 and 0.71 had the highest seasonal mean CWSI, respectively, while in Bahar, Pishtaz and Sistan cultivars, it ranged from 0.61 to 0.64 under severe drought and CWSI variation of these cultivars was less than Shiraz and Yavaros (Table 4). Gontia and Tiwari (2008) reported that the maximum CWSI of 0.52, 0.58, 0.68 and 0.89 were foundunder irrigation according to 100, 60, 40 and 20% of FC, respectively. Also, similar to our study, they reported genotypic variation among wheat cultivars for CWSI so that Bezostaya and Sandy cultivars had the highest and the lowest CWSI, respectively.

 

Relationship between CWSI and water applied

The value of CWSI was negatively correlated with water supply in wheat (Figure 2a). Linear regression showed that with decreasing water supply under stress CWSI was increased and the slope of linear regression from 525 to 283 mm water supply was more than that of well and excess watered conditions (R2=0.80). Stokcle and Dugas (1992) reported that as plants closed their stomata due to water shortage, and hence stomatal conductivity, heat flux, transpiration and the cooling effects of evaporation were decreased, the canopy temperature  and CWSI were increased, compared with well watering conditions. In the present study, Shiraz and Yavaros cultivars used more water (Table 2) and had more CWSI (Table 4) when subjected to mild and severe drought stress conditions.

 

Relationship between CWSI and net
photos-ynthesis

A negative relationship was found between CWSI and net photosynthesis (Pn) of flag leaf (R2=0.77), especially at CWSI values greater than 0.51 (Figure 2b). These results confirmed those of Yuan et al. (2004) who reported that with increasing CWSI after anthesis, Pn decreased sharply. Also, Feng et al. (2001) declared that severe drought, during grain filling period, by closing the stomata to avoid leaf transpiration, decreased Pn and increased CWSI significantly.

 

Relationship between CWSI and grain yield

Maximum grain yield was obtained in Shiraz and Yavaros cultivars under well and excess watering and CWSI in these cultivars ranged from 0.31 to 0.36, whereas by lowering water supply and increasing CWSI, grain yield in these cultivars decreased significantly (Table 4). Our results were somehow similar to those of Garrot et al. (1994) who reported that the highest grain yield (606 g/m2) was achieved at CWSI levels between 0.37 and 0.3.  These results illustrate the value of using CWSI as an indicator of crop water status and many researchers suggest that CWSI could be used to measure crop water status, improve irrigation  scheduling  and  obtain  optimum  grain

 


Table 4.  Effect of irrigation regimes on monthly and seasonal mean CWSI values and grain

yield  of wheat cultivars.

Irrigation regime

(according to filed capacity)

Wheat cultivar

 

Month

 

Seasonal mean CWSI

Grain yield

(g/m2)

 

 

April

May

June

 

 

Severe drought (50% FC)

Shiraz

0.62

0.74

0.83

0.73a

401.2 l

 

Bahar

0.41

0.63

0.79

0.61b

435.1k

 

Pishtaz

0.43

0.66

0.80

0.63b

451.1j

 

Sistan

0.47

0.65

0.80

0.64b

461.3i

 

Yavaros

0.55

0.76

0.82

0.71a

403.3l

 

 

 

 

 

 

 

Mild drought (75% FC)

Shiraz

0.46

0.67

0.70

0.61b

510.3g

 

Bahar

0.33

0.56

0.64

0.51c

501.3h

 

Pishtaz

0.38

0.59

0.59

0.52c

580.1f

 

Sistan

0.44

0.54

0.64

0.54c

592.1e

 

Yavaros

0.49

0.62

0.66

0.59b

500.4h

 

 

 

 

 

 

 

Well watered (100% FC)

Shiraz

0.23

0.38

0.47

0.36d

771.4a

 

Bahar

0.24

0.29

0.37

0.30ef

621.1c

 

Pishtaz

0.27

0.30

0.45

0.34de

602.7d

 

Sistan

0.29

0.31

0.39

0.33def

707.8b

 

Yavaros

0.29

0.33

0.43

0.35de

771. 6a

 

 

 

 

 

 

 

Excess watered (125% FC)

Shiraz

0.21

0.33

0.39

0.31def

772.1a

 

Bahar

0.19

0.35

0.36

0.30ef

622.3c

 

Pishtaz

0.23

0.29

0.32

0.28f

604.4d

 

Sistan

0.19

0.31

0.40

0.30ef

708.6b

 

Yavaros

0.22

0.34

0.43

0.33def

773.8a

                 

FC: Field Capacity. Means in each column by the same letters are not significantly different at 5% probability level using Duncan’s multiple range test.

 

 

 

yield especially under water shortage conditions (Gardner et al. 1992; Alderfarsi and Nielsen 2001; Emekli et al. 2007).    

Fischer et al. (1998) also studied the association among wheat yield progress with higher photosynthetic rate and cooler canopies and concluded that measuring photosynthetic rate, canopy temperature and CWSI, should be further investigated as potential indirect selection criteria for grain yield. In our study it was found that generally, Bahar, Pishtaz and Sistan had lower Tc-Ta and CWSI and higher grain yield compared to Shiraz and Yavaros cultivars under drought stress (Table 4 and Figure1).

The grain yield was correlated with mean seasonal CWSI values (Figure 2c) by the following polynomial Equation (3):

Y= -997.19 (CWSI)2+259.74(CWSI)+707.29    (3)                                           

where Y is grain yield (g/m2). As was shown in Figure 2c, the seasonal mean CWSI was related to wheat grain yield negatively (R2=0.86). This equation could be used for yield prediction under different water status in wheat. Predicting yield response to crop water stress is important in developing strategies and decision-making by researchers and farmers for irrigation scheduling under limited water conditions (Gardner et al. 1992; Yuan et al. 2004; Orta et al. 2004).

 

 

 

 

 

 

Figure 2. Relationships of CWSI with water supply (a), net photosynthesis (b) and grain yield (c) of wheat

 

 

 

 

 

 

 

 

 

 

 

 

Conclusion

Crop canopy temperature reflects the interactions among plants, soil and atmosphere. The application of canopy–air temperature difference was appropriate for crop water stress determination as it is non-destructive, non-contact, reliable, provides considerably precise estimation and represents actual crop water demand. The CWSI can be estimated using semi-empirical approach with observations of Tc-Ta and VPD. A negative relationship was observed between CWSI and net photosynthesis of flag leaf and water applied under different irrigation regimes. The seasonal mean CWSI was related to kernel yield of wheat, negatively and a polynomial equation (Equation 3) can be used to predict the yield potential. Indeed more CWSI could lead to less grain yield due to more water limitation. Bahar, Pishtaz and Sistan, with lower Tc-Ta, water supply, and CWSI had better performances than Shiraz and Yavaros, especially when exposed to water stress. Determination of canopy temperature and CWSI should be further investigated as potential indirect selection criteria for grain yield of wheat cultivars under semi arid conditions

Alderfasi AA and NielsenDC, 2001. Use of crop water stress index for monitoring water status and scheduling irrigation in wheat. Agricultural Water Management 47: 69–75.
Al-Faraj A, Meyer GE and Horst GL, 2001. A crop water stress index for tall fescue (Fetusca arundinacea Schreb.) irrigation decision-making: a traditional method. Commercial Agriculture 31: 107–124.
Allen RG, Pereira LS, Raes D and Smith M, 1998. Crop evapotranspiration. FAO Irrigation and Drainage Paper 56. FAO, Rome.
Alves I and Pereira LS, 2000. Non-water-stressed baselines for irrigation scheduling with infrared thermometers: a new approach. Irrigation Science 19: 101–106.
Braunworth  WS, 1989. The possible use of the crop water stress index as an indicator of evapotranspiration deficits and yield reductions in sweet corn. Journal of American Society of Horticulture Science 114: 542–546.
Emam Y, 2007. Cereal Production. 3rd Edition. ShirazUniversity Press, Shiraz, 190 pp.
Emekli Y, Bastug R, Buyuktas D and EmekliNY, 2007. Evaluation of a crop water stress index for irrigation scheduling of bermudagrass. Agricultural Water Management 90: 205–212.
Feng BL, Wang CF and Miao F, 2001. Leaf gas exchange character of low canopy temperature wheat in drought conditions. Journal of Triticale Crop 21: 48–51.
Feng BL, Yu H, Hu Y, Gao X, Gao J, Gao D and Zhang S, 2009. The physiological characteristics of the low canopy temperature wheat (Triticum aestivum L.) genotypes under simulated drought condition. Acta Physiology Plantarum 31: 1229–1235.
Fischer RA, Rees D, Sayre KD, Lu Z, Condon AG and Saavendra AL, 1998. Wheat yield progress associated with higher stomatal conductance and photosynthetic rate, and cooler canopies. Crop Science 36: 1467–1475.
Gardner  BR, NielsenDC and Shock CC, 1992. Infrared thermometry and the crop water stress index. II. Sampling procedures and interpretation. Journal of Production Agriculture 5: 466–475.
Garrot DJ, Ottman, DD, Fangmeier DD and Hunman SH, 1994. Quantifying wheat water stress with the crop water stress index to schedule irrigations. Agronomy Journal 86: 195-199.
Gontia NK and Tiwari KN, 2008. Development of crop water stress index of wheat crop for scheduling irrigation using infrared thermometry. Agricultural Water Management 95: 1144–1152.
Grimes DW, Yamada H and Hughes SW, 1987. Climate-normalized cotton leaf water potentials for irrigation scheduling. Agricultural Water Management 12: 293-304.
Howell TA, Musick JT and Tolk JA, 1986. Canopy temperature of irrigated winter wheat. Transition ASAE 29: 1692–1699.
Idso SB, Jackson RD, Pinter JR, Reginato, RJ and Hatfield JL, 1981. Normalizing the stress-degree-day parameter for environmental variability. Agriculture Meteorology 24: 45–55.
Jackson RD, Idso RB, Reginato RJ and Pinter PJ, 1981. Canopy temperature as a crop water stress indicator. Water Resource 17: 1133–1138.
Jalali-Farahani HR, Slack DC, Kopec DM, and Matthias AD, 1993. Crop water-stress index models for bermudagrass. Agronomy Journal 85: 1210–1217.
Kramer PJ, 1983. Water Relations of Plants. Academic Press, New York.
Li L, Nielsen DC, Yua Q, Mac L and Ahuja LR, 2008. Evaluating the crop water stress index and its correlation with latent heat and CO2 fluxes over winter wheat and maize in the North China plain. Agricultual Water Management 7: 1146-1155.
Monteith JL and Unsworth MH, 1990. Principles of Environmental Physics. Edward Arnold, London, pp. 243.
Orta AH, Baser I, Sehirali S, Erdem T and Erdem Y, 2004. Use of infrared thermometry for developing baseline equations and scheduling irrigation in wheat. Cereal Research 32: 363–370.
Panda, RK, BeheraSK and Kashyap PS, 2003. Effective management of irrigation water for wheat under stressed conditions. Agricultural Water Management 63: 37–56.
SAS, 2003. SAS for windows. V. 9.1. SAS Inst, Cary, USA.
StokcleCO and DugasWA, 1992. Evaluating canopy temperature-based indices for irrigation scheduling. Irrigation Science 13: 31–37.
Yuan G, Luo Y, Sun X and Tang D, 2004. Evaluation of a crop water stress index for detecting water stress in winter wheat in theNorth China Plain. Agricultural Water Management 64: 29–40.