## Services on Demand

## Article

## Indicators

- Cited by SciELO

## Related links

## Chilean journal of agricultural research

*On-line version* ISSN 0718-5839

### Chilean J. Agric. Res. vol.72 no.1 Chillán Mar. 2012

#### http://dx.doi.org/10.4067/S0718-58392012000100002

**Chilean Journal of Agricultural Research 72(1) January - March**

**RESEARCH**

**Assessing the Transferability of Transpiration-use Efficiency Models of Biomass Production**

**Comprobando la Transferibilidad de Modelos para la Producción de Biomasa basados en la Eficiencia del uso de la Transpiración Cristián Kremer**

^{1,2*}, and Claudio Stöckle

^{1}

^{1}Washington State University, Biological Systems Engineering Department, Pullman, Washington, 99163, USA.

^{2}Universidad de Chile, Facultad de Ciencias Agronómicas, Av. Santa Rosa 11315, Santiago, Chile. "Corresponding author (cristiankremer@gmail.com).

Transpiration-use efficiency *(w),* defined as the ratio of biomass produced per unit water transpired, has been used to evaluate crop performance under limited water supply. However, the lack of consistency of *w* values through different environmental conditions has not allowed, using it as a transferable parameter. Thus, simple approaches have been developed, including: 1) *w = k _{Da} D_{a}^{1}* and; 2)

*w = k*where

_{Eto}ETo'^{1}*k*and

_{Da}*k*are crop-dependent parameters, with the underlying concept that normalization by

_{ETo}*D*or

_{a}*ET*would accounts for the effects of climate variations on w, while these parameters would be reasonably constant across diverse environments. The objective of this study was to assess the transferability of

_{0}*k*and

_{Da}*k*for wheat

_{ET}^_{0}*(Triticum aestivum*L.) and maize (Zea

*mays*L.). The scarcity of experimental information and discrepancy of the methodology used, justified the use of a canopy transpiration and photosynthesis model which was developed, tested, and fitted with weather data from eight environmentally different locations to simulate values of

*w, k*and

_{Da}*k*The results indicated that

_{ETOo}.*k*and

_{Da}*k*were more variable than expected; suggesting that calibration would be desirable. A consistent trend of change of the parameter values as function of

_{ET}^_{0}*D*or

_{a}*ET*was found, which can be represented by mathematical functions, allowing transferring w,

_{0}*k*and

_{Da}*k*(maize). In contrast, the

_{ET}^_{0}*k*for wheat correlated weakly with

_{ETOo}*D*and ET

_{a}_{0}, but a low overall coefficient of variation (10%) allowed using an average value as a reasonable predictor of w.

**Key words:** Transpiration-use efficiency, models of biomass production.

Eficiencia del uso de la transpiración (w), definida como la relación entre biomasa producida por unidad de agua transpirada, se ha utilizado para evaluar productividad bajo condiciones limitadas de agua. Sin embargo, la falta de consistencia de valores *w* bajo condiciones climáticas distintas no ha permitido su uso como parámetro transferible. Consecuentemente, aproximaciones simples han sido desarrolladas, incluyendo: 1) *w = k _{Da}*

*D*y; 2)

_{a}^{1}*w = k*

_{ETo}_{ }

*ET*donde

_{0}-^{1}*k*

_{Da}_{ }y

*k*son parámetros dependientes del cultivo, y sustentado en que la normalización por

_{ETo}*D*o

_{a}*ET*absorbería el efecto del clima en la determinación de w, mientras que estos parámetros se mantendrían razonablemente constantes. El objetivo de este estudio fue determinar la transferibilidad de

_{o}*k*y

_{Da}*k*para su uso en trigo

_{ETo}*(Triticum aestivum*L.) y maíz (Zea

*mays*L.). La escasez de información experimental y metodologías usadas, justificó el uso de un modelo para estimar transpiración y fotosíntesis, el cual fue desarrollado y probado con datos de ocho regiones climáticamente distintas para simular valores de

*k*y

_{Da}*k*Los resultados indicaron que estos parámetros poseían mayor variabilidad de la esperada, sugiriendo que una calibración previa sería necesaria. Además, el cambio de estos parámetros como función de

_{ETo}.*D*o

_{a}*ET*tiene una tendencia consistente, representable por funciones matemáticas, permitiendo transferir

_{0}*w, k*y

_{Da}*k*(maíz). Por otro lado, valores de

_{ETo}*k*en trigo se correlacionaron débilmente con

_{ETo}*D*y ET

_{a}_{0}, pero un bajo coeficiente de variación (10%) permitiría el uso de un valor promedio como un predictor razonable de

*w.*

**Palabras clave:** eficiencia en el uso de la transpiración, modelos de producción de biomasa.

Agriculture is challenged by the scarcity of water resources in many regions of the world, problem that is compounded by climate variability and expected to worsen in the future. There is a raising need for tools to evaluate crop productivity as a function of water to better guide development policies and field management practices aimed at producing "more crop per drop".

Mechanistic simulation models of canopy photosynthesis and transpiration appear as suitable tools to evaluate the effect of interacting factors on water-use efficiency and productivity of crops. However, demanding parameterization and computing requirements of these models limit their applicability for long-term analysis that includes multiple species across the globe. As a result, there is a renewed interest in simple, transpiration-based models of crop productivity that can be readily applied to a large number of crop species across the range of climatic conditions where these crops are grown. Although these models were introduced as early as the beginning of the previous century, the experimental determination of the parameters (typically just one parameter) used in the models has been relatively scarce, probably due to the need of measuring crop transpiration for their determination. As a result of scarce experimental information, it is not easy to assess the variability and transferability among locations of the parameters of these simple transpiration-based models.

Attempts to develop simple relationships to predict *w* for different crops and climates can be traced back to the early 20^{th} century and later (Brigss and Shantz, 1914; Shantz and Piemeisel, 1927; de Witt, 1958; Arkley, 1963; Bierhuizen and Slatyer, 1965; Tanner, 1981; Ritchie, 1983; Tanner and Sinclair, 1983; Steduto and Albrizio, 2005). The underlying assumption has been that the parameters of these relationships are relatively constant across diverse climatic conditions, and assumption that has not been well evaluated.

Early work by Bierhuizen and Slatyer (1965) led the way to the development of a comprehensive physiologically-based description of transpiration use efficiency as follows:

* [1]*

where *N _{l}* is the net leaf photosynthesis,

*T*is the leaf transpiration, ΔC is the CO

_{l}_{2}concentration difference between the atmosphere and the CO

_{2}compensation point,

*r*is the leaf resistance to CO

_{CO2}_{2}diffusion from the surrounding air into the leaf and into the cells of the chloroplasts, ρ is the density of the air,

*e*is the vapor to air molecular weight ratio,

*D*is the vapor pressure difference between the leaf and the surrounding air, P

_{l}_{a}is the atmospheric pressure, and

*rv*is the summation of the partial resistances to water vapor flux from the leaf. Bierhuizen and Slatyer (1965) showed that the ratio

*Ni*

*T*is determined largely by

_{l}^{'1}*D*based on the following assumptions: (1) ΔC is a relatively constant crop-dependent parameter, and (2) the

_{l}*r*ratio is fairly constant in active leaves when the water stress is not severe. The authors redefined Equation [1] in a simpler expression:

_{v}r_{CO2}^{-1}*
*

here *k _{l}* is considered constant for leaves in a given crop. The authors also argued that

*N*should be proportional to canopy

_{l}T_{l}'^{1}*B T*(where

^{-1}*B*is canopy dry matter and

*T*is canopy transpiration) and, therefore,

*k*could be scaled up to the entire canopy

_{l}*k*and that

_{Da}*D*could be well represented by the air vapor pressure deficit

_{l}*D*since the leaf temperature appear to be within ± 2-3 °C of air temperature. Hence:

_{a}*w* = *BT ^{-1}= k_{Da}*

*Da*[4]

^{-1}where *k _{Da}* is obtained experimentally as the slope of the linear regression between cumulative biomass and the daily integration of the quotient

*TD*

_{a}^{-1}.Tanner and Sinclair (1983) extended the work initiated by Bierhuizen and Slatyer (1965) and Tanner (1981), developing equations to represent biomass production and transpiration of the sunlit and shaded fractions of the canopy, essentially leading to a re-derivation of Equation [4].

Although Equation [4] has been adopted as a reasonable predictor of biomass accumulation (Stòckle *et ai.,* 1994; Sinclair and Seligman, 1995), concerns have been raised about the transferability of *k _{Da}* (Kemanian

*et ai.,*2005). These authors argued that

*k*is not a "constant" for a crop, but it rather changes with environmental conditions, most noticeable

_{Da}*Da.*

Steduto and Albrizio (2005) presented field data and a discussion of the concept and mechanism of determination of k_{Da}, including C_{3}: chickpea *(Cicer arietinum* L.), sunflower *(Heiianthus annuus* L.), wheat *(Triticum aestivum* L.), and C_{4}: sorghum *(Sorghum bicoior* (L.) Moench subsp. *bicoior)* species in one location (Bari, Italy; 41^{o}03' N, 16^{o}52' E, 72 m a.s.l.). They found that their *kDa* values had large variability among species and did not match data for the same species from literature. Two explanations to their findings were mentioned: 1) the error introduced by scaling *D _{i}* to D

_{a}, especially in low

*D*conditions where leaf temperature can be several degrees larger than air temperature, and 2) the effectiveness of

_{a}*Da*normalization to represent

*D*since the latter is defined by the transpiration flux, which changes as the physiological stage of the crop changes. They proposed an alternative methodology, similar to the original work by de Witt (1958), where

_{i}*w*is a function of the evaporation rate of a reference condition as:

*w* = *BT ^{-1}= k_{ETo}*

*ET*[5]

_{0}^{-1}where

*ET*is the reference crop evaporation computed as proposed by Allen

_{0}*et ai.*(1998) and

*k*is the slope of the linear regression between cumulative biomass and the daily integration of the quotient

_{ETo}*TET*Steduto and Albrizio (2005) claim that this method would work better than Equation [4] and that

_{0}^{-1}.*kETo*appeared transferable among different climatic zones. However, this claim was based on limited data and has not been verified.

The main objective of this work was to evaluate the transferability across diverse climatic condition of *kDa* and *kETo* of wheat and maize. However, experimental data allowing the calculation of *kDa* and *kETo* is scarce and does not cover well the wide array of environmental conditions where wheat and maize are grown. In addition, the available data include differences in cultivars, crop management, methods to estimate transpiration, sampling methods for biomass, and other sources of variability and experimental error, making it difficult to evaluate the constancy of the parameters. For that reason, a canopy transpiration and photosynthesis (CTP) model was developed and tested (Kremer, 2006), and used as reference to obtain simulated values of *w, kDa,* and *kETo* under variable climatic conditions while crop and soil characteristic were held constant, allowing a more consistent evaluation of the transferability of these parameters.

**MATERIALS AND METHODS**

**The model**

An hourly time step canopy transpiration and photosynthesis (CTP) model, separating sunlit and shaded fractions of the canopy, was developed as a tool to obtain simulated values of *kDa* and *kETo.* The model simulates C assimilation (g CO2 m^{-2} ground area), and crop transpiration (kg H2O m^{-2} ground area) in response to climatic conditions, soil and plant water status, and atmospheric CO2 concentration. Transpiration, photosynthesis, stomatal conductance, and plant water uptake are solved simultaneously through an iterative numerical procedure. Daily measurements of global solar radiation, air temperature, air humidity, and wind speed are inputs to the model. Additional inputs include green leaf area index (LAI), maximum crop height, maximum LAI, assimilation rate as a function of intercellular CO2 concentration, stomatal conductance response to air vapor pressure deficit and leaf water potential, and soil characteristics (hydraulic parameters, bulk density, depth, and number and thickness of soil layers). Model performance was tested using meteorological and crop data (wheat and maize) collected at the Conservation and Production Research Laboratory, Bushland (35°11' N, 102°06' W; elevation 1170 m a.s.l.), Texas, USA, indicating the suitability of the model for the application presented in this research. A more detailed description of the CTP model and parameters for the simulation of wheat and maize transpiration-use efficiency is presented in Kremer (2006). For this study, soil water content, LAI, crop height, and crop parameters for photosynthesis and stomatal conductance were held constant during the entire simulation period at all the locations. Thus, only daily weather was variable.

**Meteorological data**

To generate a highly diverse set of conditions, daily weather data from eight locations were selected. The data were composed of daily measurements of global solar radiation (MJ m^{-2} d^{-1}), maximum and minimum air temperature (^{o}C), maximum and minimum air relative humidity, and average wind speed (m s^{-1}), encompassing the period day of the year (DOY) 120 through 239 for a year selected at random in the following locations: Concepcion del Uruguay (32°28' S, 58°16' W; elevation 20 m a.s.l.), Argentina; Landriano (45°18' N, 9°15' E; elevation 78 m a.s.l.), Italy; Temple (31°7' N, 97°4' W; elevation 208 m a.s.l.), Texas, USA; Pullman (46°45' N, 117°1' W; elevation 756 m a.s.l.), Washington, USA; Prosser (46° N, 119°7' W; elevation 380 m a.s.l.), Washington, USA; Ankara (40°7' N, 32°59' E; elevation 948 m a.s.l.), Turkey; Aleppo (36°1' N,37°18' E; elevation 430 m a.s.l.), Syria; and DOY 166 through 239 for Maricopa (33°49' N, 112°1' W; elevation 359 m a.s.l.), Arizona, USA. The main climatic characteristics of the selected period in the eight locations are presented in Table 1.

**Table 1. Mean (x) and standard deviation (sd) of weather data from eight locations and selected periods.**

*Tmax* and Tmn: maximum and minimum ir temperatures (°C), Sr: global solar radiation (MJ m^{-2} d^{-1}), *RHmax* and *RH _{mrl}:* maximum and: minimum relative humidity, ETo: reference evapotranspiration (mm d

^{-1}), D

_{a}: day time air vapor pressure deficit (k

_{Pa}) estimated as:

*Da*= 2/3 es(Tmax), (1 -

*R*where

_{min}).*e*is the saturation vapor pressure of the air in kPa at maximum.

_{s}(T_{max})**Determination of k_{Da} and k_{ETo}**

Daily transpiration and photosynthesis for wheat and maize were simulated for all locations using the CTP model and assuming well-developed and unstressed crops. The parameters

*k*(g CO

_{Da}_{2}kg

^{-1}H

_{2}O kPa) and

*k*(g CO

_{ETo}_{2}m

^{-2}) were estimated as the slope of the linear regression between cumulative photosynthesis and the accumulation of the transpiration to daytime

*D*quotient (Tanner, 1981; Condon

_{a}*et al.,*1993) and transpiration to

*ET*quotient (Steduto and Albrizio, 2005), respectively. Calculations of the parameters were done for moving 15-d intervals, shifted by 5 d throughout the 120-d period. Daily

_{0}*ET*calculation was carried out as proposed by Allen

_{0}*et al.*(1998).

**RESULTS AND DISCUSSION**

As was expected simulated

*w*values for 15-d periods across eight locations were highly variable (Table 2), with coefficient of variation of 25% and 18% for wheat and maize, respectively, with the known implication that

*w*values determined experimentally in one location may not be readily transferable to another. Normalization of

*w*by

*D*Equation [4] or

_{a}*ET*Equation [5] is expected to account for weather variability, with the parameters

_{0}*k*or

_{Da}*k*remaining reasonable constant. However, as shown in Table 2

_{ETo}*k*and

_{Da}*k*(maize) present more variability than desirable for transferring values derived in one location to another with

_{ETo}*k*variability being greater than that of

_{Da}*k*

_{ETo}.**Table 2. Mean (x), standard deviation (sd), number of data (n), and coefficient of variation (CV) for simulated values of transpiration-use efficiency (w), and the crop parameters k_{Da} and k_{ETo} from eight locations.**

As pointed out by Tanner (1981) and Steduto and Albrizio (2005), a drawback associated with *Da* normalization is usually related to the degree of error introduced on the assumption that *D _{a}* is a fair representation of

*Dl,*especially in humid environments. In these environments, transpiration rate is expected to be lower, and as a result leaf temperature should increasingly depart from air temperature making

*Dl*larger than

*Da.*To test this assumption, the values of

*kDa*obtained in environments with

*Da*less than 1 kPa were not included in the coefficient of variation (CV) analysis. Some improvement on CV was obtained for both crops however the variability still remained (wheat: 18.75%, and maize: 26.18%).

Figures 1 and 2 present

*w*as a function of the average

*D*and

_{a}*ET*of each corresponding 15-d interval. Both figures show that

_{0}*w*is not constant across environments characterized by

*Da*and

*ET0,*and has an important non-linear response to

*D*and

_{a}*ET*(Abbate

_{0}*et al.,*2004; Kemanian

*et al.,*2005). Fitted power equations appear good estimators of

*w,*with

*Da*explaining 94% and 90% of the

*w*variability for wheat and maize, respectively, and

*ET*explaining 89% and 72%, respectively. The dispersion around the fitted lines represents variability due to climate that is not accounted for

_{0}*D*or

_{a}*ET*This effect of other weather variables is less important with

_{0}.*D*and

_{a}*ET*greater than 2 kPa and 7 mm d

_{0}^{-1}, respectively.

**Figure 1. Transpiration use efficiency (w, g CO _{2} kg^{-1} H_{2}O) as a function of the air water vapor pressure deficit (daytime, kPa) for wheat and maize.**

**Figure 2. Transpiration use efficiency (w, g CO _{2} kg^{-1} H_{2}O) as a function of reference evapotranspiration ((ET_{0 } mm d^{-1}) for wheat and maize.**

Figure 3 reaffirms that the variability in Table 2 is not random, but can be explained to a large extent when *k _{Da}* values are plotted vs.

*D*or ET

_{a}_{0}. A linear equation was fitted to the

*k*values and included in the figure. It seems that

_{Da}*Da*was able to explain

*kDa*variability better than

*ET0,*which presented a larger scattering, particularly for wheat. These results confirm that: 1)

*k*is not a constant value and, 2)

_{Da}*kDa*increases when

*Da*and

*ET0*increases.

**Figure 3. Variability of the crop dependant parameter koa (g CO**2

**kg**2

^{1}H**O kPa) as a function of the daytime vapor pressure deficit**

*(Da;*kPa), and the reference evapotranspiration (ETo, mm d^{1}) for wheat and maize.**
**

The significant conclusion is that the use of Equation [4] to estimate *w* has to considerer local calibration of *kDa* to be transferable. The linear response obtained for *kDa* should facilitate transferability through field calibration based on a few points across the environmental range. The fitted linear equations included in each figure, can be used as *kDa* estimators for climatic conditions characterized by different *D _{a}* or ET

_{0}.

In an attempt to explain *k _{ETo}* variability,

*k*values were plotted vs.

_{ETo}*Da*and

*ET*(Figure 4). Fitted linear equations are included in the Figure 4, which shows that

_{0}*k*variation in wheat does not correlate well with variations in climatic conditions represented by

_{ETo}*D*or

_{a}*ET*However the relatively low CV and standard deviation (sd) determined that a mean value of 55.87 (g CO2 m

_{0}.^{-2}) can be used as a constant regardless of the climatic environment, supporting the view of Steduto and Albrizio (2005). Nevertheless, some response of

*k*for wheat when the

_{ETo}*ET0*gradient is increasing was observed (Figure 4), suggesting some benefit of using the fitted equation in situations with high evaporative demand.

**Figure 4. Variability of the crop dependant parameter kETo (g CO**2

**m**

^{2}) as a function of the daytime vapor pressure deficit (Da; kPa), and the reference evapotranspiration (ETo, mm d^{-1}) for wheat and maize.**
**

A different scenario was found in maize (Table 2). Figure 4 shows that *k _{ETo}* (maize) did not correlate well with variations in

*D*and that the variability was better explained by

_{a},*ET*although

_{0},*ET0*alone was not able to account for the entire variability due to weather. It can be concluded that

*k*is not a constant and, therefore, experimental values cannot be transferred among locations with different climate. The linear equation presented here to estimate

_{ETo}*k*as function of

_{ETo}*ET*for maize should be taken as first approximation to overcome the transferability problem.

_{0}**CONCLUSIONS**

The simulation results showed as was expected that

*w*is not constant across climatic environments, so that values determined in one location cannot be readily transferred to another. However,

*w*as a function of

*D*and

_{a}*ET*was well described by a power function, with

_{0}*Da*explaining 94% and 90% of the

*w*variability for wheat and maize, respectively, and

*ET0*explaining 89% and 72%, respectively. The dispersion around the fitted lines was much lower with

*D*and

_{a}*ET*greater than 2 kPa and 7 mm d

_{0}^{-1}, respectively.

Normalization of the k_{Da} and *k _{ETo}* parameters of two simple transpiration-use efficiency models by

*D*and

_{a}*ET*was not able to properly account for the effect of weather variability, resulting in parameters still too variable to be readily transferred across locations for both wheat and maize.

_{0}It was found that the transferability of these parameters can be dramatically improved when they are plotted against *D _{a}* (in the case of k

_{Da}) or

*ET*(in the case of

_{0}*k*), with linear functions describing well the relations and explaining 79% and 91% of

_{ETo}*kDa*variability for wheat and maize, and 71% of

*k*variability for maize. The

_{ETo}*k*for wheat correlated weakly with

_{ETo}*ET*explaining only 25% of its variability. However, the overall coefficient of variation of this parameter across eight locations was about 10%, so that the use of a constant

_{0},*k*value is not unreasonably, although is not a perfect solution.

_{ETo}The simulation-based equations presented here are offered as a first approximation to overcome the spatial transferability of w, k_{Da}, and *k _{ETo}*, but field validation will be required before adoption is recommended.

**LITERATURE CITED**

Abbate, P.E., J.L. Dardanelli, M.G. Cantero, M. Maturano, R.J.M. Melchiori, and E.E. Suero. 2004. Climatic and water availability effects on water-use efficiency in wheat. Crop Science 44:474483. [ Links ]

Allen, R.G., L.S. Pereira, D. Raes, and M. Smith. 1998. Crop evapotranspiration: Guidelines for computing crop water requirements. Irrigation and Drainage Paper 56. FAO, Rome, Italy. [ Links ]

Arkley, R.J. 1963. Relationships between plant growth and transpiration. Hilgardia 34:559-584. [ Links ]

Bierhuizen, J.F., and R.O. Slatyer. 1965. Effect of atmospheric concentration of water vapor and CO_{2} in determining transpiration-photosynthesis relationships of cotton leaves. Agricultural Meteorology 2:259-270. [ Links ]

Brigss, L.J., and H.L Shantz. 1914. Relative water requirement of plants. Journal of Agricultural Research 3:1-63. [ Links ]

Condon, A.G., R.A. Richards, and G.D. Farquhar. 1993. Relationships between carbon isotope discrimination, water use efficiency and transpiration efficiency for dryland wheat. Australian Journal of Agricultural Research 44:1693-1711. [ Links ]

de Witt, C.T. 1958. Transpiration and crop yield. Versl. Landbouwk. Onderz. 64.6. Institute of Biological and Chemical Research on Field Crops and Herbage, Wageningen, The Netherlands. [ Links ]

Kemanian, A.R., C.O. Stöckle, and D.R. Huggins. 2005. Transpiration-use efficiency of barley. Agricultural and Forestal Meteorology 130:1-11.

Kremer, C. 2006. Evaluating simple transpiration-based models of crop productivity. Ph.D. diss. Washington State University, Department of Biological Systems Engineering, Pullman, Washington, USA. [ Links ]

Ritchie, J.T. 1983. Efficient water use in crop production: Discussion on the generality of relations between biomass production and evapotranspiration. p. 29-44. *In* Taylor H.M. *et al.* (eds.) Limitations to efficient water use in crop production. ASA, CSSA, SSSA, Madison, Wisconsin, USA. [ Links ]

Shantz, H.L., and L.N. Piemeisel. 1927. The water requirement of plants at Akron, Colorado. Journal of Agricultural Research 34:1093-1190. [ Links ]

Sinclair, T.R., and N.G. Seligman. 1995. Global environmental change and simulated forage quality of wheat in non stressed Conditions. Field Crops Research 40:19-27. [ Links ]

Steduto, P., and R. Albrizio. 2005. Resource use efficiency of field-grown sunflower, sorghum, wheat and chickpea. II. Water use efficiency and comparison with radiation use efficiency. Agricultural and Forestal Meteorology 130:269-281. [ Links ]

Stockle, C.O., S. Martin, and G.S. Campbell. 1994. CropSyst, a cropping systems model: water/nitrogen budgets and crop yield. Agricultural Systems 46:335-359. [ Links ]

Tanner, C.B. 1981. Transpiration efficiency of potato. Agronomy Journal 73:59-64. [ Links ]

Tanner, C.B., and T.R. Sinclair. 1983. Efficient water use crop production: Research or re-search? p. 1-27. *In* Taylor, H.M. *et al.* (eds.) Limitations to efficient water use in crop production. ASA, CSSA, SSSA, Madison, Wisconsin, USA. [ Links ]

Received: 18 April 2011.

Accepted: 4 January 2012.