On-line version ISSN 1807-1929
Rev. bras. eng. agríc. ambient. vol.16 no.6 Campina Grande June 2012
MANEJO DE SOLO, ÁGUA E PLANTA
Potassium estimation in the soil solution based on electrical conductivity and soil water content
Estimativa de potássio na solução do solo baseada na condutividade elétrica e umidade do solo
Torquato M. de Andrade NetoI, Eugênio F. CoelhoII; José A. do V. SantanaIII; Edvaldo B. Santana JúniorIII; Márcio da S. AlvesIII
IDoutorando em ciências agrárias, Universidade Federal do Recôncavo Baiano (UFRB), Rua Embrapa, s/n, CEP 44380-000, Cruz das Almas, BA. Fone: (75) 3312-8021, E-mail: email@example.com
IIEmbrapa Mandioca e Fruticultura. Rua Embrapa, s/n, CEP 44380-000, Cruz das Almas, BA. Fone: (75) 3312-8021. E-mail: firstname.lastname@example.org
IIIUniversidade Federal do Recôncavo Baiano (UFRB), Campus Cruz das Almas. E-mail:email@example.com; firstname.lastname@example.org; email@example.com
The objective of this work was to evaluate and to validate models for estimating potassium in the soil solution as a function of bulk electrical conductivity (ECw), soil water content (q) and a soil solution electrical conductivity (ECss). Treatments consisted of using three concentrations of injecting solution of potassium chloride (1.0, 2.5 and 4.0 g L-1) which were applied by two trickle irrigation systems (microsprinkler and drip) during the first cycle of the banana crop cv. Terra Maranhão. Results showed that it is feasible to estimate potassium concentration in the soil solution from data of ECss and q obtained by time domain reflectometry (TDR) using an equation that combined a linear and a potential model. The estimated values of potassium concentration were close to the ones measured along the crop cycle under field conditions, with a mean normalized deviation of 10.0%, maximum and minimum deviation of 5.0 and 13.0%, respectively.
Keywords: TDR, trickle irrigation, fertirrigation
O objetivo deste trabalho foi avaliar e validar modelo de estimativa de potássio na solução do solo, como função da condutividade elétrica aparente (CEa), da umidade do solo (q) e da condutividade elétrica da solução do solo (CEss). Os tratamentos consistiram no uso de três concentrações de cloreto de potássio da solução de injeção (1,0, 2,5 e 4,0 g L-1) aplicadas por microaspersão e por gotejamento, durante o primeiro ciclo da cultura da bananeira cultivar Terra Maranhão. Os resultados mostraram que é viável estimar a concentração de K+ na solução do solo a partir de dados de q e CEa, obtidos por meio da técnica da reflectometria no domínio do tempo (TDR) para condições de campo, com uso de equação resultante da combinação de um modelo linear e um potencial. Os valores de K+ estimados se aproximaram dos medidos ao longo do ciclo da cultura da bananeira "Terra" nas condições de campo, com desvio normalizado médio de 10%, desvio máximo e mínimo de 5,0 e 13,0%, respectivamente.
Palavras-chave: TDR, irrigação localizada, fertirrigação
The use of fertilizers by means of irrigation water has been increasing in the irrigation districts due mainly to the increase of trickle systems and the positive effects of fertirrigation (Santana et al., 2006). Fertirrigation is being practiced without technical criteria in Brazilian irrigation districts with consequences in soil chemical properties such as transient salinity and leaching of íons with risk of groundwater contamination (Pinto, 2001). On the other hand, this technique has contributed for optimizing fertilizer use in irrigated agriculture (Oliveira & Villas-Boas, 2008).
An adequate fertirrigation is the one in which nutrients are applied at right moment according to crop needs. The soil chemical analysis demands long time since sampling in the field. The use of the soil solution might be a good alternative that provides the knowledge of ionic state of the soil solution in a short time. The soil solution may be an alternative for obtaining nutrient levels by laboratory analysis or by using kits for fast determination of levelsof these nutrients. Nutrient levels in soil solution may also be estimated by models from bulk electrical conductivity (ECw) or soil solution electrical conductivity (ECss) (Rhoades et al., 1976; Vogeler & Clotheir, 1996; Muñoz-Carpena et al., 2001; Santana et al., 2007).
Time domain reflectometry might be a feasible alternative in order to obtain a better precision of soil ion dynamics as for soil water dynamics (Nobório, 2001; Muñoz-Carpena et al., 2001) since it provides soil water content and bulk electrical conductivity data at different locations in the soil space along time.
Models of Rhoades et al. (1976) and Vogeler & Clothier (1996) are the best for adjusting ECw data as a function of soil water content and ECss (Santana et al., 2007). Muñoz-Carpena et al. (2001) have obtained relations between ECss and different ions in the soil solution like nitrate (Coelho et al., 2005) and potassium (Santana et al., 2007). Also, reasonable fitting has been noticed relating ECss and ion concentration under field conditions (Muñoz-Carpena et al., 2001; Nobório, 2001; Vogeler & Clotheir, 2001) and under controlled environment (Santana et al., 2007; Ritter et al., 2005; Regalado et al., 2005).
This work had as objective to evaluate and validate a model for estimating potassium concentration in the soil solution as a function of bulk electrical conductivity (ECw), soil water content (q) and soil solution electrical conductivity (ECss) under field conditions.
MATERIAL AND METHODS
The work was carried out in the experimental area of Embrapa Cassava & Tropical Fruits, located at Cruz das Almas municipality, Bahia State, (12º 48" S, 39º 06" W, 225 m). The climate is classified as humid to sub humid. The soil is an Yellow Alic and Distrofic Latossol silty textured with 444 g kg-1 total sand, 131 g kg-1 silt, 425 g kg-1 clay, bulk density of 1.55 kg dm-3, 0.23 m3 m-3 soil water content at the superior limit of soil water availability and 0.16 m3 m-3 for inferior limit of soil water availability for the Coastal Tablelands soils of Northeast Brazil. The chemical characteristics at the beginning of the experiment were: pH 6.3; 11.0 mg dm-3 - P; 0.06 cmolc dm-3 - K; 3.4 cmolc dm-3 - Ca + Mg; 0.09 cmolc dm-3 - Na; 1.32 cmolc dm-3 - H + Al; 3.56 cmolc dm-3 - sum of bases; CTC of 4.88 cmolc dm-3; 73% for base saturation.
The work was performed inside an experiment with banana crop, cv. Terra Maranhão in a 0.1 ha area to evaluate the effect of three irrigation water with different potassium chloride concentrations (1.0; 2.5 and 4.0 g L-1) on chemical soil attributes. Two irrigation systems considered were: micro sprinkler with 43.0 L h-1 emitter per four plants and drip with three 4.0 L h-1 emitters per plant in one lateral line per plant row. The amount of water per irrigation event was calculated based upon reference evapotranspiration and crop coefficients as recommended by Allen et al. (1998). Fertirrigation was applied every week and the injection solution was based on Borges & Caldas (2004).
A parametric model was fitted to data of soil water content and electrical conductivity as a function of potassium concentration. At least twenty different values of independent and dependent variables of the model were used in the optimization process. The optimization process was based upon minimization of the sum of square differences between obser-ved and estimated dependent variables (objective function).
Model parameters for estimating electrical conductivity and potassium concentration in the soil solution
Experimental activity was carried inside banana crop where TDR probes and water samplers were installed in the soil at depths of 0.20 and 0.40 m at 0.30 m from plant between two emitters along lateral line (drip system) and between plant and emitter (micro sprinkler). TDR probes were installed at these locations for evaluating bulk electrical conductivity (ECw) and soil water content (q). TDR probes of 0.10 m length were built (Silva et al., 2005) and soil water content was determined by Eq. 1 (Ledieu et al., 1986). Electrical conductivity of soil solution was estimated (Giese & Tiemann, 1975) and their values were corrected to 25 ºC temperature, according to Eq. 2 proposed by Richards (1954).
where: q - soil water content, cm3 cm-3 ECw - soil bulk electrical conductivity
Soil solution samples were collected at the same time of readings of ECw and q, these readings were made every 15 min beginning at the suction application with a vacuum hand pump until soil solution be collected two hours after suction. Soil solution electrical conductivity (ECss) was measured by means of a desk conductivity meter and potassium concentration [:K]: was estimated by using a kit for quick determination (Card Horiba). The average reading of q and ECw taken during soil solution suction in water sampler, ECss and [:K]: data were related by mathematical models. The model of Vogeler & Clothier (1996) was used for estimating ECss as a function of ECw and q according to Eq. 3:
where: ECss - soil solution electrical conductivity ECw - Soil bulk electrical conductivity a, b, c and d - parameters of Vogeler & Clothier (1996) equation
The data of [:K]: and ECss were related by a potential function in order to yield ECss as a function [:K]::
The substitution of Eq. 4 in Eq. 3 resulted in the model:
Statistical evaluation of model performance
The adjustment of the mathematical model and data was accomplished by means of minimization of root square deviation among estimated and observed values. The statistical indices ME (mean errors), RMSE (root mean square error), d (Willmott agreement index) and R2 (goodness of fit) were used for model evaluation. The root mean square error (RMSE) was defined by the equation:
where: n - number of data Oi - measured value Ei - estimated value
The absolute mean error (MEA) was calculated according to the equation:
The normalized error mean (MEN) was also calculated:
Borges et al. (2010) reported that these statistical parameters are good indicators for model efficiency. The units used for the deviations were the same as used in the evaluation of variables in order to facilitate interpretation of results (Legates & McCabe Jr., 1999). Also a regression analysis considering measured and estimated dependent variable with zero intercept was evaluated by the angular coefficient and goodness of fit.
Validation of the model under field conditions
Soil solution samples were collected every 15 days in each plot with three replications by using water samplers installed radially to the micro sprinkle at 0.30 m from plant at depths of 0.20 and 0.40 m during one cycle of banana crop. The water samplers were installed at 0.3 m from plant in the direction of plant-dripper at a fixed distance of 0.15 m from the dripper at two depths (0.20 and 0.40 m) with three replications. The porous cup of the sampler were located in the wetted volume between two drippers.
Water samples were colleted followed by a suction with a vacuum hand pump (-70 kPa) and the suction period lasted for two hours. Afterwards, soil solution samples were taken to laboratory and ECss and [:K]: were determined. Therefore, data of ECw, ECss and [:K]: were registered at an interval of 15 days during a whole banana cycle for the three irrigation waters with different potassium chloride concentrations (1.0; 2.5 and 4.0 g L-1).
Once the parameters for the models for estimating [:K]: as a function of ECw and [:K]: as a function of ECss and q were defined (Eq. 5); the model was applied to the data of ECss, ECw and q which were collected every 15 days along the crop cycle. The efficiency of the model for estimating [:K]: concentration related to the measured concentrations along banana cycle was calculated by using Nash & Suttcliffe (1970) equation.
where: EF - efficiency of the model Oi - potassium measured Pi - potassium estimated
Besides efficiency, RMSE and MEA were evaluated for comparing data of [:K]: observed in the field and estimated by modeling. Statistical analysis was made considering [:K]: observed along crop cycle and estimated at 0.30 m depth for the two irrigation systems by using the t test for comparing means of estimated and measured values.
RESULTS AND DISCUSSION
Soil solution electrical conductivity as a function of bulk electrical conductivity and soil water content
Table 1 illustrates the adjustment of Vogeler & Clothier (1996) model to data of ECss as a function of ECw and q for 1.0, 2.5 and 4.0 g L-1 concentrations of K in water by micro sprinkler and drip irrigation. Results showed that 70, 78 and 93% of ECss variations were due to the variation of ECw and q for 4.0; 2.3 e 1.0 g L-1 concentrations, respectively, in drip system. These results were in agreement to the ones obtained by Silva et al. (2005), they verified that 92% of variation of ECss were due to variations of ECw and q in a similar study. Santana et al. (2006) noticed that the empirical model of Vogeler & Clothier (1996) resulted in best estimate of ECss as a function of ECw and q compared to other models, in which, 82.7% of ECss variations was explained by ECw e q variations.
Vogeler & Clothier (1996) models fitted reasonably to the data of ECss as a function of ECw and q only for 4.0 g L-1 concentration. This result is in agreement to the ones of Santana et al. (2007), who observed that 79% of ECss variation was explained by variation of ECw and q using Vogeler & Clothier (1996) model for a loam soil. The minimum RMSE value and the smaller MEA, which indicate the best fitting of data to Eq. 3 was verified for the irrigation water concentration of 1.0 g L-1 for drip system. The ECss model underestimated the measured values for 1.0 g L-1 concentration in case of drip system and for 4.0 g L-1 in case of micro sprinkler. The model overestimated the measured values with small fitting coefficients in case of micro sprinkler. This result did not agree to the one verified by Santana et al. (2007). These authors found that ECss from the model underestimed the measured values by about 5%.
Potassium concentration as a function of bulk electrical conductivity and soil water content
Table 2 shows model parameters as a result of adjustment of Eq. 5 to potassium data as a function of ECw and q. Also correlation and difference coefficients are shown in the same table for both drip and sprayer systems and considering the irrigation water concentrations. Eq. 5 explained 92.42; 81.52 and 84.29% of variations in [:K]: due to variations of ECw and q for water concentrations respectively of 1.0; 2.5 and 4.0 g L-1, in case of drip system. These coefficients were near to the ones obtained by Santana et al. (2007), who noticed that 81% of variations of [:K]: were explained by variations of ECw and q in a sandy silty soil. The same occurred for sprayer system, i.e., the Eq. 5 explained 87, 84 e 80.1% of variation of [:K]: as a function of ECw and q for concentrations of potassium chloride in irrigation water of 1.0, 2.5 and 4.0 g L-1. The results were in agreement to the ones of Santana et al. (2007), the model explained 81,78% of [:K]: as function of ECw and q in a sandy silty soil.
Treatments that showed the smallest value of RMSE were the ones that presented the largest goodness of fit as observed by Borges et al. (2010), except for the concentration of 2.5 g L-1 in micro sprinkler. Eq. 5 overestimated about 2.5% the measured values of [:K]: for concentrations of 4.0 g L-1 in drip system and about 1.4% the measured values for 2.5 g L-1 concentration in micro sprinkler, however the model underestimated measured [:K]: values in the range of 1.16 to 1.45% for the other irrigation water concentrations in both irrigation systems.
The goodnes of fit for relations Kes= b Kmed (Kes - estimated potassium concentration and Kmed - measured concentration) for all irrigation water concentrations were larger than those obtained by Santana et al. (2007), who worked with different models for estimating [:K]: and obtained 74% of [:K]: variations explained by ECw and q variations. There was statistical difference of deviation modules between the averages of measured and estimated [:K]: by the t test, considering the three concentrations of irrigation water. Only the average deviation (5.68%) for 1.0 g L-1 concentration differed (P < 0.05) fom the others which were close to 12,0% in case of micro sprinkler. The average deviation between [:K]: estimated and measured (5.42%) for 4.0 g L-1 concentration was statistically different from and smaller than the others whose average deviation varied from 10.70 to 13.83% in case of drip irrigation system.
Validation of the model under field conditions
Figure 1 depicts concentrations of [:K]: observed from soil solution and estimated by Eq. 5 along time for concentrations of potassium chloride in irrigation water of 1.0; 2.5 and 4.0 g L-1, in case of drip system (Figure 1A, 1B and 1C) and of micro sprinkler (Figure 1D, 1E and 1F). Eq. 5 presented efficiencies of 0.94, 0.97 and 0.94 for concentrations of 1.0; 2.5 and 4.0 g L-1 in case of sprinkler system. These efficiencies were close to values obtained in case of drip indicating that the model fitted reasonably to field data for both irrigation systems (Table 3). The average variation of normalized errors (MEN) ranged from 5.0 to 13.0% for both irrigation systems.
MEN for 4.0 g L-1 irrigation water concentration differed statistically from the others by t test with a value of 5.4%, in case of drip system. The deviation was about 11% for the other treatments (concentrations). The model efficiency stayed between 0.89 and 0.94 for drip system. The fiited values by Eq. 5 have showed dispersions (RMSE), in the range of 1,43 to 4,78 mg L-1.
1. The model of ECss as a function of ECw and soil water content underestimated measured K íon concentration values for 1.0 g L-1 potassium chloride concentration of irrigation water in drip system and for 4.0 g L-1 in micro sprinkler. The model overstimated measured values with poor fit to data for the other concentrations in case of sprinkler.
2. Eq. 5 showed larger performance coefficients for estimating potassium ion concentration compared to Eq. 3 for estimating ECss.
3. Eq. 5 showed efficiencies ranged from 0.89 to 0.97 for all concentrations of irrigation water potassium chloride, with mean variation of normalized errors of 5 to 13% in both irrigation systems.
Allen, R. G.; Pereira, L. S.; Raes, D.; Smith, M. Crop evapotranspiration: Guidelines for computing crop water requirements. Rome: FAO, 1998, 297p. FAO. Irrigation and Drainage Paper, 56 [ Links ]
Borges, A. L.; Caldas, R. C. Teores de nutrientes nas folhas de bananeira, cv. Pacovan, sob irrigação. Ciência e Agrotecnologia , v.28, p.1099-1106, 2004. [ Links ]
Borges, V. P.; Oliveira, A. S. de.; Coelho Filho, M. A.; Silva, T. S. M.; Pamponet, B. M. Avaliação de modelos de estimativa da radiação solar incidente em Cruz das Almas, Bahia. Revista Brasileira de Engenharia Agrícola e Ambiental, v.14, p.74-80, 2010. [ Links ]
Coelho, E. F.; Vellame, L. M.; Coelho Filho, M. A. Sondas de TDR para estimativa da umidade e condutividade elétrica do solo com uso de multiplexadores. Revista Brasileira de Engenharia Agrícola e Ambiental, v.9, p.475-480, 2005. [ Links ]
Giese, K.; Tiemann, R. Determination of the complex permittivity from the sample time domain reflectometry. Advanced Molex Relaxes Processes, v.7, p.45-49, 1975. [ Links ]
Ledieu, J.; d. Ridder, P.; d. Cleck, P.; Dautrebande, S. A method measuring soil water moisture by time-domain reflectometry. Journal of Hydrology, v.88, p.319-328, 1986. [ Links ]
Legates, D. R.; McCabe Jr., G. J. Evaluating the use of "goodness-of-fit" measures in hydrologic and hydroclimatic model validation. Water Resources Research, v.35, p.233-241, 1999. [ Links ]
Muñoz-Carpena, R.; Regalado, C. M.; Alvarez-Benedí, J.; Socorro, A. R.; Pérez, N. Determinación simultánea mediante TDR del transporte de agua y un soluto salino en el suelo. In: López, J. J.; Quemada, M. Temas de investigación en zona no saturada. 5.ed., Pamplona: Universidade Pública de Navarra, 2001. p.1-7. [ Links ]
Nash, J. E; Sutcliffe, J. V. River flow forecasting through conceptual models part I - A discussion of principles. Journal of Hydrology, v.10, p.282-290, 1970. [ Links ]
Nobório, K. Measurement of soil water content and electrical conductivity by time domain reflectometry: A review. Computers and Electronics in Agriculture, v.31, p. 213-237, 2001. [ Links ]
Oliveira, M. V. A. M de; Villas-Boas, R. L. Uniformidade de distribuição do potássio e do nitrogênio em sistemas de irrigação por gotejamento. Engenharia Agrícola, v.28, p 95-103, 2008. [ Links ]
Pinto, J. M. Fertirrigação em fruticultura irrigada. Revista Irrigação & Tecnologia Moderna, v.49, p.14-23, 2001. [ Links ]
Regalado, C. M.; Ritter, A.; Álvarez-Benedí, J.; Muñoz-Carpena, R. Simplified method to estimate the green-ampt wetting front suction and soil sorptivity with the Philip-dunne falling-head Permeameter. Vadose Zone Journal, v.4, p.291-299, 2005. [ Links ]
Richards, L. A. Diagnosis and improvement of saline and alkali soils. Washington: United States Salinity Laboratory, 1954. 160p. USDA Agriculture Handbook, 60. [ Links ]
Rhoades, J. D.; Raats, P. A.; Prather, R. J. Effects of liquid phase electrical conductivity, water content and surface conductivity on bulk soil electrical conductivity. Soil Science Society of America Journal, v.40, p.651-655, 1976. [ Links ]
Ritter, A, R. Muñoz-Carpena, C. M. Regalado, M. Javaux; Vanclooster, M. Using TDR and Inverse modeling to characterize solute transport in a layered agricultural volcanic Soil. Vadose Zone Journal, v.4, p.300-309. 2005 [ Links ]
Santana, G. da S.; Coelho, E. F.; Silva, T. S. M.; Ramos, M. M. Estimativa da condutividade elétrica da solução do solo a partir do teor de água e da condutividade elétrica aparente do solo. Revista Brasileira de Engenharia Agrícola, v.26, p.672-685, 2006. [ Links ]
Santana, G. S.; Coelho, E. F.; Silva, T. M.; Ramos, M. M. Relação entre potássio na solução do solo, umidade e condutividade elétrica aparente do solo. Revista Brasileira de Engenharia Agrícola e Ambiental, v.11, p.142-151, 2007. [ Links ]
Silva, T. S. M.; Coelho, E. F.; Paz, V. P. S.; Vellame, L. M.; Santana, G. S. Teor de potássio na solução do solo com uso da técnica de reflectometria no domínio do tempo. Revista Irriga, v.10, p.393-402, 2005. [ Links ]
Vogeler, I.; Clothier, B. E. Characterizing water and solute movement by TDR and disk permeametry. Soil Science Society of America Journal, v.60, p.5-12, 1996. [ Links ] Willmott, C. J. On the validation of models. Physical Geography, v.2, p.184-94, 1981. [ Links ]
Protocolo 225.10 – 10/12/2010
Aprovado em 28/03/2012