Ciencia e investigación agraria
versión ISSN 0718-1620
Cienc. Inv. Agr. vol.38 no.2 Santiago ago. 2011
Cien. Inv. Agr. 38(2):243-251. 2011
Effects of physicochemical soil properties of five agricultural soils on herbicide soil adsorption and leaching
Relación entre la adsorción y la lixiviación de
Claudio Alister, Manuel Araya, and Marcelo Kogan
Escuela de Ciencias Agrícolas, Universidad de Viña del Mar. Av. Agua Santa 7055. Viña del Mar, Chile.
Once pesticides reach the soil, there are several factors that affect their soil behavior. To identify the principal soil and herbicide properties that control their adsorption and leaching, a study of five Chilean agricultural and forest soils was performed. Simazine, diuron, terbuthylazine and MCPA were applied to the top of 45-cm tall by 12-cm diameter disturbed soil columns, filled with either an Andisol, Ultisol, Entisol or one of two Inceptisol soils. After herbicide applications, each lysimeter received 24 mm of simulated rain every 24 hours for five days. Once water percolation stopped, the lysimeters were divided into five sections and herbicide concentrations were quantified using High Pressure Liquid Chromatography. Relationships between soil physicochemical properties, herbicide sorption (adsorption and desorption) and herbicide leaching were determined. All herbicides exhibited the least depth reached from the Andisol soil (10 cm) and the highest from the Ultisol soil (45 cm). The principal soil property that affected herbicide adsorption was the soil organic carbon content, specifically the fulvic acid-humins fraction. Soil leaching was related to the inverse of soil adsorption (1/Kd), cation exchange capacity, humic substances content and herbicide pKa. These results suggest that it is possible to develop simple quantitative models to predict the soil-leaching properties of pesticides.
Key words: Diuron, fulvic acids, humic substances, MCPA, triazines.
Una vez que un plaguicida llega al suelo comienzan a actuar una serie de procesos que van controlando la permanencia de este compuesto en el sitio de aplicación. En base a esto se realizó un estudio en cinco suelos Chilenos para determinar el efecto de las propiedades físico-químicas del suelo y de los herbicidas sobre la lixiviación de esos compuestos. Simazina, diuron, terbutilazina y MCPA fueron aplicados sobre columnas de suelo disturbado de 45 cm de altura y 12 cm de diámetro llenas con suelos Andisol, Ultisol, Entisol y dos Inceptisoles. Una vez aplicados los herbicidas se simuló una lluvia de 24 mm cada 24 horas por 5 días. Cuando la percolación del agua terminó cada columna fue dividida en cinco secciones a las cuales se les determinó, la concentración de los herbicidas mediante HPLC-DAD. Todos los herbicidas mostraron su menor lixiviación en el suelo Andisol (10 cm) y la mayor en el suelo Ultisol (45 cm). El contenido de carbono orgánico fue la principal propiedad del suelo relacionada con la adsorción de los herbicidas, específicamente la fracción ácido fúlvico + humina. La lixiviación los cuatro herbicidas estudiados se relacionó con el inverso del coeficiente de adsorción (1/Kd), capacidad de intercambio catiónico, fracción húmica del suelo y el pka de los herbicidas. Estos resultados permiten pensar en la posibilidad de desarrollar un índice cuantitativo simple para predecir la lixiviación en condiciones productivas.
Palabras clave: Ácido fúlvico, diuron, MCPA, substancias húmicas, triazinas.
When a pesticide reaches the soil, a series of processes that control the persistence of this compound in the site of application occur. The degradation, volatilization, lixiviation and adsorption of the compound to the soil are considered among the main processes that are associated with the dynamics of a pesticide in soil. It is generally accepted that the leaching of a compound depends mainly on its degradation and its sorption to the soil. Therefore, most of the simple indexes used to evaluate the risk of pesticide leaching include degradation and sorption as the main factors (Laskowski et al. 1982; Gustafson, 1989; Spadotto, 2002), as well as some other soil parameters (e.g., density, organic carbon, soil humidity) in the case of more complex indexes (Rao et al, 1985; Meeks and Dean, 1990). However, according to some authors, adsorption and desorption are the phenomena that regulate the speed and magnitude of leaching, and a compound should not be affected by other processes while it is adsorbed to the humic-argillic system (Charnay et al, 2005; Papiernik et al., 2006). Thus, the incorporation of organic soil amendments increase the adsorption of some pesticides, decreasing leaching and increasing persistence. Nevertheless, as pluviometry is varied through simulated rain, persistence is modified, but not leaching, which indicates that adsorption itself could be a better signal of the deep movement of a pesticide compared with other processes, such as degradation (Fouque-Brouard and Fournier, 1996; Muller et al, 2003; Alister et al., 2005; Kogan et al., 2007).
A study based on this information was carried out to determine the mobility of four herbicides in five different soils and the relationships of the physicochemical properties of the soils and the herbicides with the sorption and lixiviation of the herbicides.
MATERIALS AND METHODS
Five soils with different physicochemical properties were selected (Table 1). Analyses of the principal physicochemical properties (CE, pH, cation exchange capacity, organic carbon and texture) were made according to Kalra and Maynar (1991). The methodology proposed by Rocha et al. (1998) and IHSS (2009) were used to determine the humic substances (HS), humic acids (HA) and fulvic acids + humins (FAH) fractions (Table 1).
Disturbed soil columns
Three 50-cm long by 12-cm diameter PVC columns were filled to a depth of 45 cm with each soil, which were sieved previously to <2 mm. Conic base caps made of PVC were installed on these columns and then, filled with quartz sand before soil filling and hermetically fitted to the column. A glass wool cord was placed intermingled in the quartz sand and in direct contact with the soil of each column, allowing the water movement inside the column of disturbed soil to flow similarly as obtained naturally in the soil (Boll et al, 1992).
Once the columns were installed, a micro-aspersion system was arranged in the upper part of the columns and used to create a saturated soil environment. The selected herbicides were applied 48 hours after water percolation reached the bottom of the columns at doses equivalent to those indicated in Table 2. The herbicides were applied with a micropipette over the column surface using a concentric distribution of 50µL drops, for a total volume of 1 mL of herbicide solution applied to each column. A daily rain of 24 mm was simulated for 5 days 24 hours after the application of herbicides, over a period of 4 minutes per event.
Two days after the water percolation stopped in the columns, each column was opened longitudinally by making two lateral cuts, exposing the soil profile, which was separated into five sections: 0-10, 10-20, 20-30, 30-40 and under 40 cm. The herbicides were quantified by HPLC for each soil section (different depths).
Quantification of herbicide
For extractions, 20 g of soil were shaked at 300 rpm for 90 min with 40 mL of methanol. Subsequently, the extracts were filtered (paper Whatman No. 1) and dried in a rotary evaporator. The herbicide residues were resuspended in 1 mL of methanol and transferred to a 1.5 mL glass vial to be analyzed using high-performance liquid chromatography (HPLC; Hitachi Elite LaChrom L-2300 model) with a Diode Array Detector (DAD; Hitachi Elite LaChrom L-2450 model) fitted for reading wavelengths between 190 and 400 nm. The HPLC was equipped with a LiChrospher® 100 RP-18 5µm column of 125 mm long. The specific analytical conditions for each herbicide are listed in Table 3.
Determination of adsorption coefficients
The adsorption coefficient (Kd) and the desorp-tion percentage were determined for each herbicide for each soil. A total of 6 mL of a 0.01 M CaCl2 solution with an herbicide concentration equivalent to 50% of the solubility of each herbicide (Table 2) were put into 15-mL capped polypropylene centrifuge tubes with 3 g of each soil. These suspensions were shaked for 8 hours at 180 rpm at 21±1 °C in darkness.
Finally, the tubes were centrifuged for 15 min at 5,000 rpm, and 1 mL of each supernatant was filtered through a 0.45µm fiberglass filter and used for quantification by HPLC, according to the conditions described above. The amount of adsorbed herbicide was calculated as the difference between the concentration of herbicide determined in the control tubes (herbicide solution without soil) and the remnant quantified in the solution after centrifugation. Equation 1 was used to determine the Kd, where Kd is the adsorption coefficient (L kg-1), Cs is the amount of adsorbed herbicide (mg kg-1), and Ce is the concentration of herbicide in solution after agitation (mg L-1).
Cs= Kd* Ce 
After the supernatant sample was taken to determine the Kd, the remaining supernatant was eliminated from each tube (4.0 mL) and the same volume was replaced with a 0.01 M CaCl2 solution (without herbicide). These new suspensions were agitated for 4 hours following the procedure indicated above. This desorption step was performed three times, and the amount of desorbed herbicide was determined based on the total sum of the herbicide quantified in each wash, according to the methodology described before. Controls were included for the different herbicide solutions during the whole process, and no losses due to hydrolysis or herbicide adsorption to the polypropylene tubes were detected.
The program SAS® was used to analyze the results. Tests of multiple correlation (PROC CORR) were made among the soil parameters and the values of adsorption and desorption. The procedure of multiple linear regression (PROC REG) was used to quantify the effects of the physicochemical properties of the soil and the herbicides on leaching (the depth reached by an herbicide and the percentage that leached to below 10 cm) using the step-wise procedure, with a variable input and output significance of P<0.1.
RESULTS AND DISCUSSION
The adsorption of herbicides and their desorption percentages were variable, depending on the soil type (Table 4). The highest Kd values and lowest percentages of desorption occurred in the Andisol soil. In contrast, the lowest Kd values and the highest desorption percentages were determined for the Ultisol soil. Among the remaining soils studied, the Entisol soil had the second highest adsorptions for the four herbicides, with a texture composed with 83% of sand (Table 1).
The behaviors of most herbicides are generally influenced by the content of organic carbon in soil (Torrents and Jayasundera, 1997; Haberhauer et al, 2002; Hua Guo et al, 2008). The results of this work show that the effect of soil organic carbon (OC) on the adsorption of these herbicides might be explained by the composition of the humic substances (HS) present in this organic fraction. In general, the Kd values obtained for all of the herbicides in the Entisol soil, which had an OC content of 0.9 %, were only slightly higher than those for the two Inceptisol soils, with OC contents over 2%. This may have been due to the similar HS contents of these three soils (Table 1 and 4). An OC-adsorption relationship was observed when correlations between the coefficients of herbicide adsorption and the physicochemical properties of the soils were examined (Table 5), which is consistent with previous reports for simazine, diuron and terbuthylazine (Bollag and Meyers, 1992; Chefetz et al, 2004; Cao et al, 2008). Moreover, this work suggests that this relationship is specifically related to the fulvic acids + humins (FAH) content present in the humic fraction of the soil (Table 5).
The adsorption of MCPA, unlike the other three herbicides, had a high correlation with the pH and the cation exchange capacity (CEC) of the soils (Table 5). The pH values of the Ultisol and Andisol soils (Table 1) are close to the pKa of MCPA (Table 2). Thus, the increases in the adsorption of MCPA in these soils could be associated with their higher CECs due to a higher amount of this herbicide molecule being present in cationic form in these soils (Table 4). In the case of the Ultisol soil, the higher CEC is derived from the clay content, and in the Andisol soil, it is associated with the FAH content, which is highly related to the adsorption of herbicides like MCPA (Iglesias et al, 2009) (Table 5).
In general, the desorption process did not have a highly significant relationship with any parameter, except for the clay content (Table 5). The contrast between the high correlation of the desorption percentages of herbicides with the soil clay content, and not with the OC content, and the high correlation between the adsorption of these herbicides and the soil OC content (especially the HS content), but not the soil clay content, indicates the presence of a certain degree of hysteresis in the adsorption-desorption process. Therefore, the adsorption mechanism of the herbicides in these soils might differ from the desorption mechanism (Weber and Weiling, 1998). However, it is important to consider that the estimation of the desorption process is associated with a high degree of uncertainty because the adsorption-desorption process requires long periods of time to equilibrate and, in some cases, may take weeks or months to approach a real balance (Ball and Roberts, 1991; Weiling and Weber, 1998). Subsequently, the times for which desorption measurements were taken in this study might have been too short to achieve a realistic balance.
The leaching of herbicides in the studied soils depended on the type of soil and pesticide, as evident in Figure 1. Thus, none of the herbicides moved further than 10 cm deep in the Andisol soil; in contrast, all of the herbicides reached the maximum column depth in the Ultisol soil (45 cm). In the Entisol soil, only MCPA reached the 45-cm depth. For the two Inceptisol soils, different patterns for sima-zine, diuron and terbuthylazine leaching were observed. Only MCPA showed a similar behavior in both soils.
|Figure 1. Leaching of four herbicides in disturbed soil columns. A) Andisol soil, B) Entisol soil, C) Inceptisol soil, D) Inceptisol soil and E) Ultisol soil.|
The physicochemical properties of the soil (pH, CEC, OC, HS, HA, FAH, clay, silt and sand), the physicochemical properties of the herbicides (solubility, pKa, and log(Kow)) and specific soil-herbicide parameters (Kd and % desorption) that had the best correlations with the leaching of the different herbicide compounds were selected, and multiple regression equations were determined based on the selected parameters. The results of these regressions show that the leaching of the slightly ionic herbicides (simazine, terbuthylazine) and the non-ionic one (diuron) could be mainly explained by the inverse value of Kd (1/Kd), for both the depth reached and the amount mobilized to depths greater than 10 cm (Table 6). On the other hand, the Kd and HS content of the soil best explained the soil depth reached and the amount of herbicide leached beyond 10 cm for the more ionic herbicide MCPA (Table 5). When considering all of the herbicides, the main parameters determined by regression to explain the depth of movement and the amount of leaching below 10 cm included 1/Kd, CEC, HS and pKa. The inclusion of the herbicide pKa in the regression might occur because three of four herbicides exhibit an ionic behavior (Table 2).
According to the results, the herbicide-soil relationship, measured as the adsorption of herbicide to the soil, might be the main factor determining leaching (depth and mobilized mass) of herbicides. This is coherent with Crosby (1998), who has suggested that there should be a direct relationship between adsorption and the distance of pesticide movement in soil. Likewise, Kogan et al. (2007) found that simple indexes that include adsorption as the main parameter in their equation have better fits with the movement of simazine, diuron, oxifluorfen, pendi-methalin and flumioxazin in the field than the indexes that include degradation, solubility and some hydraulic properties of the soil in addition to adsorption.
According to the results obtained in this study with disturbed soil columns, the phenomenon of herbicide leaching was primarily related to the adsorption of these compounds in the different soils, mainly due to their contents of ful-vic acids + humins in the organic fraction of the soils, followed by parameters such as the cation exchange capacity, the content of humic substances in the organic fraction and the acid dissociation constant (pKa).
Although these results may not be entirely extrapolatable to field conditions, in which the soil has different physicochemical conditions that would not be included in this study (especially in terms of structure), they help to determine the focus of the study of these relationships in natural conditions (undisturbed) to develop simple quantitative models to estimate pesticide leaching that consider only the adsorption of a pesticide to a soil and some other simple parameters.
The authors wish to thank FONDECYT (Chilean Fund for Science and Technology) for funding Project Number N° 11085003.
Alister, C., R. López, and M. Kogan. 2005. Simazine dynamics in a vineyard soil at Casablanca valley, Chile. Pest Manag. Sci. 61: 1083-1088. [ Links ]
Ball, W., and R. Roberts. 1991. Long term sorption of halogenated organic chemicals by aquifer materials: Part 2. Intraparticle diffusion. Environ. Sci. Technol. 25: 1237-1249. [ Links ]
Boll, J., J. Selker, and T. Steenhuis. 1992. Fiberglass wicks for sampling water and solutes in the vadose zone. Soil Sci Soc Am Journal 56: 701-707. [ Links ]
Charnay, M., S. Tuis, Y. Coquet, and E. Barriuso. 2005. Spatial variability in 14C-herbicide degradation in surface and subsurface soils. Pest Manag. Sci. 61: 845-855. [ Links ]
Chefetz, B., Y. Bilkis, and T. Polubesova. 2004. Sorption-desorption behavior of triazines and phenylurea herbicides in Kishon river sediments. Water Research 38: 4383-4394. [ Links ]
Crosby, D. 1998. Environmental Toxicology and Chemistry. Chapter 4: Environmental Transport. p. 51-67. Oxford University Press, New York, USA. [ Links ]
Cao, J., H. Guo, M. Zhu, L. Jiang, and H. Yang. 2008. Effects of SOM, surfactant and pH on the sorption-desorption and mobility of prometryne in soils. Chemosphere 70: 2127-2134. [ Links ]
Fouque-Brouard, C., and J. Fournier. 1996 Adsorption-desorption and leaching of phenylurea herbicides on soils. Talanta 43: 1793-1802. [ Links ]
Gustafson, D. 1989. Groundwater ubiquity score: a simple method for assessing pesticide leachability. Env. Toxicol. Chem. 8: 339-357. [ Links ]
Haberhauer, G, B. Temmel, and M. Gerzabek. 2002. Influence of dissolve humic substances on the leaching of MCPA in a soil column experiment. Chemosphere 46: 495-499. [ Links ]
Hua Guo, J., H. Mei Zhu, L. Jiang, and H. Yang. 2008. Effects of SOM, surfactant and pH on the sorption-desorption and mobility of prometryne in soils. Chemosphere 70: 2127-2134. [ Links ]
Iglesias, A., R. López, D. Gondar, J. Antelo, S. Fiol, and F. Arce. 2009. Effect of pH and ionic strength on the binding paraquat and MCPA by soil fulvic and humic acids. Chemosphere 76: 107-113. [ Links ]
IHSS (International Humic Substances Society). 2009. International Humic Substances Society. Available online at: ihttp://ihss.gatech.edu/ihss2/soilhafa.html (Website accessed: October 20th, 2009). [ Links ]
Kalra, Y., and D. Maynard. 1991. Methods manual for forest soil and plant analysis. Information Report NOR-X-319, Forest Canada Northwest Region, Northern Forestry Center. Edmonton, Alberta, Canada. [ Links ]
Kogan, M., S. Rojas, P. Gómez, F. Suárez, J.F. Muñoz, and C. Alister. 2007. Evaluation of six pesticides leaching indexes using field data of herbicides application in Casablanca Valley. Water Sci. Tech. 56: 169-178. [ Links ]
Laskowski, D., C. Goring, P. Mc Call, and R. Swann. 1982. Terrestrial environment. In: Environmental Risk Analysis for Chemicals. R.A. Conway (ed.). Van Nostrand Reinhold Co. New York, USA pp. 198-240. [ Links ]
Meeks Y., and J. Dean. 1990. Evaluating groundwater vulnerability to pesticides. Journal Water Resource Planning Manag. 116: 693-707. [ Links ]
Muller, K., R. Smith, T. James, P. Holland, and A. Rahman. 2003. Spatial variability of atrazine dissipation in an allophanic soil. Pest Manag. Sci. 59: 893-903. [ Links ]
Rao, P., A. Hornsby, and R. Jessup. 1985. Indices for ranking the potential for pesticide contamination of groundwater. Proceedings of the soil and Crop Science Society of Florida 44: 1-8. [ Links ]
Papiernik, S., W. Koskinen, P. Cox-Rice, S. Clay, N. Werdin-Pfisterer, and K. Norberg. 2006. Sorption-desorption of imidacloprid and its metabolites in soil and vadose zone materials. J. Agric. Food Chem. 54: 8163-8170. [ Links ]
Rocha, J, A. Henrique, and M. Furlan. 1998. An alternative methodology for extraction of humic substances from organic soils. J. Braz. Chem. Soc. 9: 51-56. [ Links ]
Spadotto, C. 2002. Screening method for assessing pesticide leaching potential. Pesticidas. Revista Ecotoxicologia e Meio Ambiente 12: 69-78. [ Links ]
Torrents, A., and S. Jayasundera. 1997. The sorption of noionic pesticides onto clays and the influence of natural organic carbon. Chemospehere 35: 1549-1565. [ Links ]
Weiling, H., and W. Weber. 1998. Hysteresis in the sorption and desorption of hydrophobic organic contaminants by soils and sediments. 1. A comparative analysis of experimental protocols. J. Cont. Hydrol. 31: 129-148. [ Links ]
Weber, W., and H. Weiling. 1998. Hysteresis in the sorption and desorption of hydrophobic organic contaminants by soils and sediments. 2. Effects of soil organic matter heterogeneity. J. Cont. Hydrol. 31: 149-165. [ Links ]
Received October 29, 2010. Accepted April 28, 2011.
Corresponding author: firstname.lastname@example.org