Stalk yield of sugarcane cultivars under different water regimes by subsurface drip irrigation

Produtividade de colmos de cultivares de cana-de-açúcar submetidas a diferentes regimes hídricos por gotejamento subsuperficial

ABSTRACT

This study aimed to evaluate the stalk yield of ten sugarcane cultivars (RB962962, RB931011, RB931530, RB98710, RB92579, RB867515, RB863129, SP791011, RB72545 and VAT90212) subjected to water deficit, full irrigation and water surplus by subsurface drip irrigation, during three cropping seasons (2011-2014). The experiment was conducted at the Experimental Field of Embrapa Meio-Norte, Teresina, Piauí State, Brazil, in dystrophic Red Yellow Argisol. The cultivars RB962962 (162.3 Mg ha^-1) and RB867515 (158.5 Mg ha^-1) have better stalk yield compared with other cultivars in all water regimes and cropping seasons.

Key words:
Saccharum officinarum L.; irrigation management; drip irrigation; cultivar selection

RESUMO

Objetivou-se, neste estudo, avaliar a produtividade de colmos de dez cultivares de cana-de-açúcar (RB962962, RB 931011, RB931530, RB98710, RB92579, RB867515, RB863129, SP791011, RB72454 e VAT90212) submetida aos regimes hídricos de deficiência, irrigação plena e excedente hídrico, por gotejamento subsuperficial, durante três ciclos de cultivo (2011-2014). O ensaio foi conduzido no Campo Experimental da Embrapa Meio-Norte, Teresina, PI, em Argissolo Vermelho-Amarelo Distrófico. As cultivares RB962962 (162,3 Mg ha^-1) e RB867515 (158,5 Mg ha^-1) apresentam melhor produtividade de colmos em comparação às demais cultivares, em todos os regimes hídricos e ciclos de cultivo.

Palavras-chave:
Saccharum officinarum L.; manejo de irrigação; fertirrigação; seleção de cultivares

Introduction

In the Northeast region, sugarcane agribusiness is already well consolidated in the states of Alagoas, Paraíba and Pernambuco. However, the crop has occupied areas in non-traditional regions, especially in the states of Tocantins, Maranhão and Piauí, due to their characteristics of soil, topography, climate (solar radiation) and favorable logistic conditions of exportation (MAPA, 2006). However, these potential areas are situated in regions with semi-arid to dry sub-humid climate, with irregular rainfall (800 to 1200 mm per year) (Bastos & Andrade Júnior, 2014) and that do not meet the water demand of the crop during its cycle, requiring the use of irrigation to achieve satisfactory yield levels.

Among the pressurized irrigation methods used in sugarcane cultivation, sprinkler irrigation and the localized methods stand out, such as subsurface drip irrigation (SDI). In the Northeast region, SDI has been widely used because of the greater adaptation to the characteristics of the sugarcane cultivation (Ferreira Júnior et al., 2014) and better water application efficiency (Gil et al., 2008).

Studies conducted in different production environments have demonstrated increment in the stalk yield of sugarcane genotypes in response to irrigation (Carvalho et al., 2009; Gava et al., 2011; Oliveira et al., 2011; Macêdo et al., 2012; Silva et al., 2014). However, for the genotype-environment interaction, it is crucial to conduct tests aiming at the expression of the productive potential of the sugarcane cultivars in different water regimes and production environments (Silva et al., 2008).

Thus, this study aimed to evaluate the stalk yield of ten sugarcane cultivars subjected to different water regimes by subsurface drip irrigation under the edaphoclimatic conditions of Teresina-PI, Brazil, during three cropping seasons.

Material and Methods

The experiments were carried out at the Experimental Field of Embrapa Meio-Norte, in Teresina-PI, Brazil (05° 05’ S; 42° 48’ W and altitude of 74.4 m). The soil of the area is classified as dystrophic Red Yellow Argisol (Melo et al., 2014) and its chemical and physico-hydraulic characteristics (Donagema et al., 2011) are presented in Table 1. The local climate is dry sub-humid, megathermal, with moderate water surplus in the summer and 32.2% of the potential evapotranspiration concentrated from September through November. The extract of the climatological water balance (historic series from 1980 to 2013) records an 8-month period of water deficit, extending from May to December (Bastos & Andrade Júnior, 2014).

Table 1
Chemical and physico-hydraulic characterization of the soil of the experimental area

• ¹ Moisture at field capacity measured at -10 kPa
• ² Moisture at permanent wilting point measured at -1500 kPa
• ³ AWC – available water capacity
• ⁴ Easily available water (50% of AWC)

The genotypes were evaluated during three crop cycles (plant cane, 1^st ratoon and 2^nd ratoon), during the cropping seasons of 2011/2012, 2012/2013 and 2013/2014 (Table 2). Ten genotypes were evaluated (RB962962, RB 931011, RB931530, RB98710, RB92579, RB867515, RB863129, SP791011, RB72454 and VAT90212), subjected to three water regimes (WRs), during the stage of full development of the crop: water deficit (WD), full irrigation (FI) and water surplus (WS). The terms WD, FI and WS refer to the partial, full and excess meeting of the sugarcane water demand, in the stage of full development, because the crop coefficient recommended in this phase by FAO is equal to 1.0 (Silva et al., 2012). The WRs were imposed with the application of different irrigation depths, based on the reference evapotranspiration (ETo): WD - replacement of 50% ETo, FI - replacement of 100% ETo and WS - replacement of 150% ETo. The different WRs were always applied on Mondays, Wednesday and Fridays, replacing the ETo accumulated in the period between two consecutive irrigations.

Table 2
Dates of beginning and end of planting/emergence, application of the different water regimes, harvest and total crop cycle in each cropping season

• ¹ Cycle (d), number of days between the end of planting/emergence and the end of harvest

Before and after the period of imposition of the different WRs (DWR), a uniform irrigation depth equivalent to 100% ETo was applied in all treatments. During the rainy season, irrigations were applied only when soil water storage reached 50% of available water capacity (AWC). Irrigation was suspended 30 days before harvest. The reference evapotranspiration (ETo) was calculated through the Penman- Monteith method (Allen et al., 1998), with daily data of air temperature, relative air humidity, wind speed and global solar radiation, obtained by an automatic agrometeorological station close to the experiment.

A subsurface drip irrigation system was used, with polyethylene drip lines and emitters with nominal flow rate of 2 L h^-1, spaced by 0.5 m. The laterals were spaced by 2 m and buried at depth of 0.25 m. Each drip line was buried in the center of two sugarcane rows, spaced by 0.5 m. Water meters were installed to control and measure the water volume applied in each regime.

The adopted experimental design was randomized blocks, with four replicates, and treatments arranged in split plots (WRs in the plots and cultivars in the subplots). The plots consisted of one 10-m-long double row of each cultivar. The evaluated area of the plot comprehended the same 10 m in the central length of the plot.

Soil tillage consisted of one plowing and two cross harrowings. Two months before planting, soil liming was performed through the application of 1 Mg ha^-1 of dolomitic limestone. Planting was manually performed. The seedstalks were cut, maintaining three buds per set, which were distributed in the planting furrows, to achieve 15 budS m^-1.

The fertilizers (macro- and micronutrients), in the three crop cycles, were applied totally through fertigation, with weekly frequency for macronutrients and monthly frequency for micronutrients, based on the following recommendation: 120 kg ha^-1 of N (urea), 120 kg ha^-1 of P₂O₅ (MAP), 100 kg ha^-1 of K₂O (potassium chloride) and a formulation of micronutrients composed of 6.5 kg ha^-1 of boric acid, 7 kg ha^-1 of zinc oxide, 6 kg ha^-1 of copper oxide, 11 kg ha^-1 of manganese oxide and 1 kg ha^-1 of sodium molybdate (Andrade Júnior et al., 2012).

At harvest, the stalks of the evaluated area of each experimental plot were weighed using a PR30‑3000 dynamometer, with precision of 500 g, to determine stalk yield (tons of cane per hectare - TCH, in Mg ha‑1). The daily records of total rainfall were obtained from an automatic agrometeorological station close to the experiment. The effective rainfall was estimated by the methodology of the Soil Conservation Service of the United States (USDA), available in the software Cropwat v. 8.0 (FAO, 2009).

Stalk yield data were subjected to the joint analysis of variance and compared by the Scott-Knott test, at 0.05 probability level. The joint analysis of the data is possible, because the relationship between the highest and lowest residual mean squares of the individual tests was not higher than seven (Gomes, 2000).

Results and Discussion

The total rainfalls (P) were equal to 1,237.8 mm (2011-2012), 1,166.9 mm (2012-2013) and 1,431.7 mm (2013-2014). Since there was not an accentuated variation in the number of rainy days (NRD) in the three cropping seasons (from 98 to 117 days), the rainfall efficiency (Ef) was very similar, with values of 58.5% (2011-2012), 56.4% (2012-2013) and 51.8% (2013-2014) (Figure 1). In the cropping season of 2011-2012, during the period of imposition of the different water regimes (DWRs), there was an effective rainfall (Pef) of 187.8 mm, which particularly damaged the application of the water deficit (WD) regime for the evaluated varieties. In the 1^st and 2^nd ratoons, the effective rainfall contributed equally to the expression of the productive potential and biometry of the evaluated cultivars, and the probable detected differences were due to the greater contribution of the irrigation applied notably in the period prior to the rainy season.

Figure 1
Daily rainfall recorded during the three crop cycles

The total water depth (TWD) (effective rainfall + irrigation), in the evaluated WRs, did not vary in an accentuated way in the three cropping seasons (Table 3). In the WD regime, the total water depth applied oscillated from 1,156.5 mm (2013-2014) to 1,248.2 mm (2011-2012), while in the WS regime this oscillation was from 1,462.1 mm (2011-2012) to 1,579.1 mm (2012-2013) (Table 3). In percent terms, the mean TWD applied in the WS regime (1,510 mm) was only 24.4% higher than that applied in the WD regime (1,213.7 mm). These TWD values were lower than those recommended for 12-month-old sugarcane (1,500 to 2,000 mm) (Silva et al., 2014), a reflex of the better water application efficiency of the subsurface drip irrigation (Oliveira et al., 2011).

Table 3
Irrigation depths (ID) applied in response to the imposed water regimes (WR) and effective rainfall (Pef) recorded during the cultivation cycle of the sugarcane cultivars in the three cropping seasons (2011-2014)

• DWR, different water regimes; Pef, effective rainfall; I, irrigation; WR, water regime; WD, water deficit; FI, full irrigation; WS, water surplus

On the other hand, the TWD applied during the period of DWR imposition varied substantially between the WRs applied in the three cropping seasons. In 2011-2012, the contribution of irrigation to the TWD applied in the three WRs was lower than those in the other cropping seasons, due to the greater contribution of Pef (187.8 mm), which particularly damaged the imposition of the WR regime. Because of that, the lowest irrigation depths were recorded in this cropping season, which varied from 138.7 mm (WD) to 352.6 mm (WS). However, in 2012-2013 and 2013-2014, the contribution of Pef was very small (9.1 mm), which increased the contribution of the irrigation depth (ID), in all WRs (Figure 1). The values of ID in 2012-2013 were higher than those in 2013-2014, due to the superiority of the daily ETo (6.1 mm) and the longer period of DWR imposition (87 days) (Table 2).

In the period after DWR application, the ID applied in 2012-2013 was 330.1 mm, superior to those applied in 2011-2012 (184.3 mm) and 2013-2014 (171.6 mm), precisely to complement the greater ETo–Pef difference of 363.4 mm occurred in this period (Figure 1). In addition, the rainfall concentrated in only 98 days, with an effectivity of 56.4% (Figure 1). This combined effect of the variation of Pef and ID, especially during DWR application, promoted high oscillation in soil water storage between the WRs, affecting the expression of the productive potential of the evaluated cultivars.

The joint analysis of the TCH data revealed that there was no significant interaction (p > 0.05) between the periods, water regimes and cultivars, indicating that the interannual variability of the imposed WRs was not able to promote alteration in the productive performance of the evaluated cultivars (Gava et al., 2011; Oliveira et al., 2011). The TCH of the cultivars showed significant interaction (p < 0.05) of periods versus water regimes (Figure 2A), periods versus cultivars (Figure 2B) and water regimes versus cultivars (Figure 2C), a reflex of the different application of the WRs in each cropping season.

Figure 2
Stalk yield (TCH) of the evaluated cultivars in response to the water regimes (WR) in the three crop cycles. (A) interaction WR versus PP; (B) interaction PP versus C; (C) interaction WR versus C. PP: Planting period; C: Cultivar; WD: water deficit, FI: full irrigation, WS: water surplus

For all cultivars, the means of TCH in the cropping season of 2011-2012 were superior to those of 2012-2013 and 2013-2014, a consequence of the better water availability in the soil in this season, especially due to the occurrence of rainfalls that totaled 187.8 mm during the period of DWR imposition (Figure 1). On the other hand, in the cropping season of 2012-2013, the means of TCH were lower because of the strong water restriction in the soil, particularly under WD and FI regimes, during the period of DWR imposition (Figure 2A). The mean values of TCH under the water regimes of FI (142.5 Mg ha^-1) and WS (154.9 Mg ha^-1) were higher than those of WD (125.1 Mg ha^-1) in all crop cycles evaluated (Figure 2A), reinforcing the importance of adequate water supplementation to sugarcane, to achieve high stalk yields (Silva et al., 2014).

The cultivars with best performances were RB867515 (210.9 Mg ha^-1) and RB931011 (200.0 Mg ha^-1), which produced 24.7 and 18.2% more stalks in comparison to the mean yield of the experiment in 2011-2012 (169.2 Mg ha^-1). However, in 2012-2013, the most critical cropping season with respect to the water availability in the soil, the cultivar RB962962 (155.6 Mg ha^-1) produced 33.7 and 10.5% more stalks than the mean TCH of 2012-2013 (116.3 Mg ha^-1) and the overall mean of the experiment (140.8 Mg ha^-1) (Figure 2B).

The cultivars VAT90212, RB92579 and RB72454 form a group with intermediate behavior regarding stalk production, with mean TCH of three crop cycles of 151.3, 148.0 and 145.5 Mg ha^-1, respectively, which are 7.4, 5.1 and 3.3% higher than the overall mean of TCH of the experiment (Figure 2B). Oliveira et al. (2011) also found superior performance of TCH for the cultivars RB92579 (255.6 Mg ha^-1) and RB72454 (191.6 Mg ha^-1), in cultivation under conventional big gun sprinkler irrigation (total water depth of 1,396.6 mm), plant-cane cycle and full irrigation, under the conditions of soil (dystrophic Yellow Argisol) and climate (rainy tropical) of Carpina-PE.

On the other hand, the cultivars SP791011 (122.2 Mg ha^-1), RB863129 (125.9 Mg ha^-1), RB931530 (126.1 Mg ha^-1) and RB98710 (129.1 Mg ha^-1) showed the lowest performance regarding stalk production and the TCH values were 13.3, 10.6, 10.4 and 8.3% lower in relation to the overall mean of the experiment (140.8 Mg ha^-1), respectively (Figure 2B). The overall performance of these cultivars was substantially affected by the low yields in TCH under the condition of water deficit in the soil. Lower performances in TCH for the cultivars SP791011 (155.0 Mg ha^-1) and RB863129 (181.5 Mg ha^-1) were also reported by Oliveira et al. (2011), in cultivation irrigated by conventional big gun sprinklers (total water depth of 1,396.6 mm), plant-cane cycle, under the conditions of soil (dystrophic Yellow Argisol) and climate (rainy tropical) of Carpina-PE. The values of TCH of these cultivars were superior to those of the present study, because they were not subjected to periods of water deficit in the soil.

The cultivars RB962962 (162.3 Mg ha^-1) and RB867515 (158.5 Mg ha^-1) showed higher productive performance in relation to the others, in all cropping seasons and WRs evaluated (Figures 2B and 2C). Under conditions of both low and high water availability in the soil, the cultivars RB867515 (156.7 Mg ha^-1, under WD, and 175.6 Mg ha^-1, under WS) and RB962962 (141.7 Mg ha^-1, under WD, and 174.1 Mg ha^-1, under WS) showed similar productive behavior (Figure 2C), which is a desirable characteristic, considering the high rainfall instability of the Northeast region. Silva et al. (2014) also obtained superior performance in terms of TCH for the cultivar RB867515, in the cycles of plant cane (140.9 Mg ha^-1) and 1^st ratoon (118.8 Mg ha^-1), under subsurface drip irrigation (total water depths of 1,880 mm, in plant cane and 1,714 mm, in 1^st ratoon), under the conditions of soil (eutrophic Red Latosol) and climate of Jaú-SP. On the other hand, the cultivars RB931530 (103.7 Mg ha^-1), SP791011 (107.3 Mg ha^-1), RB931011 (111.7 Mg ha^-1) and RB98710 (111.4 Mg ha^-1) should not be recommended for cultivation in the region of Teresina-PI, because of the low capacity of stalk production, notably under the condition of reduced water availability in the soil. Sugarcane plants subjected to water deficit exhibited reduced daytime CO₂ assimilation, due to the decrease in stomatal conductance (diffusive factor), reduction in the activity of PEPcase and RuBisCO (biochemical factors) and lower photochemical activity (Machado et al., 2013), with negative effects on stalk production.

Increments in stalk yield under regimes of FI and WS can be attributed to the adequate water availability during the growth period, which promotes greater availability of nutrients to the sugarcane root system (Oliveira et al., 2011). Under this condition, there is an increase in tillering and elongation of the stalks due to the large transport of carbohydrates from the roots to the shoots (Inman-Bamber & Smith, 2005; Singh et al., 2007).

Conclusion

The cultivars RB962962 (162.3 Mg ha^-1) and RB867515 (158.5 Mg ha^-1) showed better performance of stalk yield in comparison to the others, in all water regimes, in the three cropping seasons.

Literature Cited

• Allen, R. G.; Pereira, L. S.; Raes, D.; Smith, M. Crop evapotranspiration: Guidelines for computing crop water requirements. Rome: FAO, 1998. 300p. Irrigation and Drainage Paper, 56
• Andrade Júnior, A. S.; Bastos, E. A.; Ribeiro, V. Q.; Duarte, J. A. L.; Braga, D. L.; Noleto, D. H. Níveis de água, nitrogênio e potássio por gotejamento subsuperficial em cana-de-açúcar. Pesquisa Agropecuária Brasileira, v.47, p.78-84, 2012. https://doi.org/10.1590/S0100-204X2012000100011
  » https://doi.org/10.1590/S0100-204X2012000100011
• Bastos, E. A.; Andrade Júnior, A. S. Boletim agrometeorológico de 2013 para o município de Teresina, Piauí. Teresina: Embrapa Meio-Norte, 2014. 38p. Documentos 228
• Carvalho, C. M. de; Azevedo, H. M. de; Dantas Neto, J.; Farias, C. H. de A.; Silva, C. T. S.; Gomes Filho, R. R. Rendimento de açúcar e álcool da cana-de-açúcar submetida a diferentes níveis de irrigação. Revista Brasileira de Ciências Agrárias, v.4, p.72-77, 2009. https://doi.org/10.5039/agraria.v4i1a12
  » https://doi.org/10.5039/agraria.v4i1a12
• Donagema, G. K.; Campos, D. V. B.; Calderano, S. B.; Teixeira, W. G.; Viana, J. H. M. Manual de métodos de análise de solos. Rio de Janeiro: Embrapa Solos, 2011. 230p.
• FAO - Food and Agriculture Organization of the United Nations. Cropwat 8.0. Rome: FAO, 2009. &lt;http://www.fao.org/nr/water/infores_databases_ cropwat.html&gt;. 15 Jul. 2012.
  » http://www.fao.org/nr/water/infores_databases_ cropwat.html
• Ferreira Júnior, R. A.; Souza, J. L.; Escobedo, J. F.; Teodoro, I.; Lyra, G. B.; Araújo Neto, R. A. Cana-de-açúcar com irrigação por gotejamento em dois espaçamentos entrelinhas de plantio. Revista Brasileira de Engenharia Agrícola e Ambiental, v.18, p.798-804, 2014. https://doi.org/10.1590/1807-1929/agriambi.v18n08p798-804
  » https://doi.org/10.1590/1807-1929/agriambi.v18n08p798-804
• Gava, G. J. C.; Silva, M. A.; Silva, R. C.; Jerônimo, E. M.; Cruz, J. C. S.; Kölln, O. T. Produtividade de três cultivares de cana-de-açúcar sob manejos de sequeiro e irrigado por gotejamento. Revista Brasileira de Engenharia Agrícola e Ambiental, v.15, p.250-255, 2011. https://doi.org/10.1590/S1415-43662011000300005
  » https://doi.org/10.1590/S1415-43662011000300005
• Gil, M.; Sinobas, L. R.; Juana, L.; Sanchez, R.; Losada, A. Emitter discharge variability of subsurface drip irrigation in uniform soils: Effect on water-application uniformity. Irrigation Science, v.26, p.451-458, 2008. https://doi.org/10.1007/s00271-008-0116-1
  » https://doi.org/10.1007/s00271-008-0116-1
• Gomes, F. P. Curso de estatística experimental. 14.ed. Piracicaba: Nobel, 2000. 477p.
• Inman-Bamber, N. G.; Smith, D. M. Water relations in sugarcane and response to water deficits. Field Crops Research, v.92, p.185-202, 2005. https://doi.org/10.1016/j.fcr.2005.01.023
  » https://doi.org/10.1016/j.fcr.2005.01.023
• Macêdo, G. A. R.; Costa, E. L. da; Viana, M. C. M.; Ferreira, J. J.; Pires, J. F.; Freire, F. M. Características agronômicas e químicas das variedades de cana-de-açúcar RB83-5486 e RB86-7515 sob irrigação e sequeiro. Revista Brasileira de Engenharia Agrícola e Ambiental, v.16, p.599-603, 2012. https://doi.org/10.1590/S1415-43662012000600002
  » https://doi.org/10.1590/S1415-43662012000600002
• Machado, D. F. S. P.; Lagoa, A. M. M. A.; Ribeiro, R. V.; Marchiori, P. E. R.; Machado, R. S.; Machado, E. C. Baixa temperatura noturna e deficiência hídrica na fotossíntese de cana-de-açúcar. Pesquisa Agropecuária Brasileira, v.48, p.487-495, 2013. https://doi.org/10.1590/S0100-204X2013000500004
  » https://doi.org/10.1590/S0100-204X2013000500004
• MAPA - Ministério da Agricultura, Pecuária e Abastecimento. Plano nacional de agroenergia 2006-2011. 2.ed. Brasília: Embrapa Informação Tecnológica, 2006. 110p.
• Melo, F. B.; Andrade Júnior, A. S.; Pessoa, B. L. O. Levantamento, zoneamento e mapeamento pedológico detalhado da área experimental da Embrapa Meio-Norte em Teresina, PI. Teresina: Embrapa Meio-Norte, 2014. 47p. Documentos 231
• Oliveira, E. C. A.; Freire, F. J.; Oliveira, A. C.; Simões Neto, D. E.; Rocha, A. T.; Carvalho, L. A. Produtividade, eficiência de uso da água e qualidade tecnológica de cana-de-açúcar submetida a diferentes regimes hídricos. Pesquisa Agropecuária Brasileira, v.46, p.617-625, 2011. https://doi.org/10.1590/S0100-204X2011000600007
  » https://doi.org/10.1590/S0100-204X2011000600007
• Silva, M. A.; Arantes, M. T.; Rhein, A. F. de L.; Gava, G. J. C.; Kolln, O. T. Potencial produtivo da cana-de-açúcar sob irrigação por gotejamento em função de variedades e ciclos. Revista Brasileira de Engenharia Agrícola e Ambiental, v.18, p.241-249, 2014. https://doi.org/10.1590/S1415-43662014000300001
  » https://doi.org/10.1590/S1415-43662014000300001
• Silva, M. de A.; Silva, J. A. G. da; Enciso, J.; Sharma, V.; Jifon, J. Yield components as indicators of drought tolerance of sugarcane. Scientia Agricola, v.65, p.620-627, 2008. https://doi.org/10.1590/S0103-90162008000600008
  » https://doi.org/10.1590/S0103-90162008000600008
• Silva, T. G. F. da; Moura, M. S. B. de; Zolnier, S.; Soares, J. M.; Vieira, V. J. de S.; Gomes Júnior, W. F. Requerimento hídrico e coeficiente de cultura da cana-de-açúcar irrigada no semiárido brasileiro. Revista Brasileira de Engenharia Agrícola e Ambiental, v.16, p.64-71, 2012. http://dx.doi.org/10.1590/S1415-43662012000100009
  » http://dx.doi.org/10.1590/S1415-43662012000100009
• Singh, P. N.; Shukla, S. K.; Bhatnagar, V. K. Optimizing soil moisture regime to increase water use efficiency of sugarcane (Saccharum spp. Hybrid complex) in subtropical India. Agricultural Water Management, v.90, p.95-100, 2007. https://doi.org/10.1016/j.agwat.2007.02.008
  » https://doi.org/10.1016/j.agwat.2007.02.008

Publication Dates

• Publication in this collection
  Mar 2017

History

• Received
  04 Apr 2016
• Accepted
  25 Nov 2016