Revista Mexicana Ciencias Agrícolas   volume 12   number 2   February 15 - March 31, 2021

DOI: https://doi.org/10.29312/remexca.v12i2.2425

Article

Mobility of atrazine in two types of soil in the state of Puebla

Ana Gabriela Rodríguez-Camacho1

Alvaro Sampieri2

José Víctor Tamariz-Flores1, 3

Ricardo Dario Peña-Moreno4

Sonia Emilia Silva-Gómez1, 5

Rosalía Castelán-Vega1, 3§

1Postgraduate in Environmental Sciences-Meritorious Autonomous University of Puebla. IC6 building (ground floor), University City, Puebla, Mexico. ZC. 72592. (gab.roca@hotmail.com). 2Faculty of Chemical Engineering-Meritorious Autonomous University of Puebla. Av. San Claudio and 18 Sur, University City, Puebla. ZC. 72592. (alvaro.sampieri@correo.buap.mx). 3Department of Research in Agricultural Sciences-Meritorious Autonomous University of Puebla. Av. 14 sur 6301 IC1 building, University City, Puebla. ZC. 72570. (jose.tamariz@correo.buap.mx). 4Chemistry Center-Meritorious Autonomous University of Puebla. IC8 Building, University City, Av. San Claudio Esq. 24 Sur, Col. San Manuel, Puebla. ZC. 72570. (ricardo.pena@correo.buap.mx). 5Department of Sustainable Development-Meritorious Autonomous University of Puebla. Av. 14 sur 6301 IC2 building, Col. San Manuel, Puebla. ZC. 72570. (soemsigo@gmail.com).

Abstract

Surface and groundwater contamination by herbicides such as atrazine occurs by leaching through the soil profile. The objective of the study was to evaluate the mobility of atrazine in plots dedicated to maize cultivation, each with a different type of soil located in the municipalities of San Nicolás Buenos Aires (Regosol) and Los Reyes de Juárez (Calcisol), Puebla, Mexico. Displacement experiments were conducted on soil packed columns (8.5 cm in diameter and 30 cm in length). An aqueous concentration of 6.25 g L-1 of atrazine (Gesaprim®) was applied to each column. First, a pore volume of a soil restorative solution (without atrazine) was used and subsequently, 30 volumes of irrigation water pore containing atrazine were applied. The leachates and three depths of the soil columns (0-2 cm, 2-10 cm and 10-20 cm, respectively) were analyzed. Pearson’s correlation results showed that organic material, texture and cation exchange capacity were the physicochemical characteristics of the soil that affect the mobility of atrazine. Atrazine had a delay coefficient of 1.08 in Regosol and 1.03 in Calcisol. The transport of the herbicide was nine times faster in the Regosol compared to Calcisol. Atrazine was retained in a greater amount in the first two centimeters of soil columns. The degraded mass percentage of atrazine was 94% for Regosol and 69% for Calcisol.

Keywords: atrazine mobility, Calcisol, corn crop, Regosol.

Reception date: December 2020

Acceptance date: March 2021

Introduction

Farmers protect their crops with agrochemicals such as herbicides to fight weeds (Carrero and Planes, 2008). They are organic substances, which can present one or more elements such as S, C, N, Cl, P and that when contact with the environment present transformations at the physical, chemical and biological level  (Josep and Figueras, 1998). That is, they generate complex systems to study, due to adsorption and absorption phenomena in soils, in cultives, rhizosphere, interactions with microorganisms, small vertebrates and invertebrates.

For example, some elements are volatilized towards the troposphere, either by chemical or microbial degradation, whose metabolites are not necessarily innocuous to the biosphere (López-Geta et al., 1992). One of the most commonly used herbicides in Mexican agriculture is atrazine, C8H14ClN5  (CAS: 1912-24-9, 6-chlorine-N-ethyl-N-isopropyl-1,3,5-triazine-2,4-diamine, which is characterized by act on a systemic level and being broad spectrum (Kim et al., 2011). It is used for the emerging control of broadleaf weeds and pasture, acts at root (rhizosphere) and accumulates in the leaves preventing the photosynthesis of plants (Marquez et al., 2009).

Farmers generally apply an atrazine concentration of 0.1- 4 kg ha-1 year-1, in crops of sorghum maize, pineapple, sugar cane, wheat and various types of pastures (Fernández  et al., 2001; Villada-Canela, 2006). When atrazine is used to increase yields in maize production, there is a negative impact on soil fertility and water quality in streams and aquifers (Hernández and Hansen, 2011; Lerch  et al., 2011).  It is común, that the atrazine is applied that by direct spray on the area to be treated, a significant amount is deposited on the surface of the soil, having low sorption and high potential for percolation, which with the effect of rains is quickly transported to the surface waters and into the soil matrix.

Studies have identified that the constant use of atrazine generates, over time, resistance in weeds (decreasing their efficiency), so farmers increase the dose and periodicity of application of the herbicide. Consequently, it causes ‘intermittent and cumulative cycles’ of atrazine, severely increasing contamination of soil, surface water and aquifers (Abdelhafid et al., 2000; Azevedo et al., 2000; Graymore et al., 2001; Shukla et al., 2003; Hernández and Hansen, 2011; Lerch et al., 2011; Mudhoo and Garg, 2011).

The mobility of atrazine in the soil depends on the texture, organic matter content, soil drainage and amount of water applied, either by rain or irrigation. For example, clay soils, high in organic matter, tend to retain atrazine on the surface and their presence in drained water is limited (Hang et al., 2010), while, in sandy soils, it tends to infiltrate and be transported with drained water (Montoya et al., 2006). Heavy rains and the application of excessive irrigation, increases the risk of migration of the herbicide to bodies of water by processes of runoff and infiltration  (Javanshir et al., 2012).

In addition, the adsorption and degradation processes of atrazine are the main natural attenuation mechanisms that control the migration of this herbicide in water and soil. Other studies indicate that adsorption of atrazine depends on soil texture and composition, pH and the applied amount of herbicide (Hang and Sereno, 2002). González-Marquez and Hansen (2009) evaluated the effect of soil parameters on natural attenuation of atrazine, and they noted that when the organic matter content is higher, its adsorption is higher. Adsorbed atrazine is less available to participate in other dissipative processes, such as biodegradation and transport to bodies of groundwater and surface water.

The objective of this research was to evaluate the mobility of atrazine in packed columns of two types of soil that are dedicated to corn crop, the soils are located in the municipalities of San Nicolás Buenos Aires and Los Reyes de Juárez in the state of Puebla, Mexico. According to data from the Secretariat of Agriculture, Livestock, Rural Development, Fisheries and Food (SAGARPA) in 2018, a maize harvest was recorded in the municipalities of 29 951 t and 3 000 t, respectively (SIAP, 2019). It is important to contextualize that the agronomic management of atrazine in the area is a research model to understand what agricultural practices favor the migration of the herbicide and thereby generate recommendations for optimal management.

Materials and methods

Description of the study area and selection of sampling units

The area of study of this research was carried out in two municipalities of Puebla, Mexico (Figure 1). The municipality of San Nicolás Buenos Aires is located between coordinates 19° 08’ 36” and 19 20’ 00” north latitude and 97° 28’ 36” and 97° 34’ 54” west longitude of the Greenwich meridian (Figure 1A). Its approximate extension is 195.19 km2. It is located in the morphological region of the plains of San Juan at an average altitude of 2 380 meters above sea level. It has a predominantly temperate-subhumid climate with rains in summer (C(w1)) in the southern part of the municipality and semi-dry-temperate with rains in the northern part, without rain during the winter (BS1kw). Its main economic activity is agriculture of temporary.

Figure 1. A) Location of the municipalities analyzed. Experimental plots; B) municipality San Nicolás Buenos Aires; and C) municipality Los Reyes de Juárez.

The predominant soil is regosol, although Feozems, Fluvisols and Lithosols can also be identified (INEGI, 2017). The second sampling municipality corresponds to Los Reyes de Juárez located between the geographical coordinates 18° 55’ 36” and 19° 01’ 06” north latitude and the meridians 97° 47’ 54” and 97° 52’ 06” of western longitude (Figure 1C). It has an extension of 30.55 km2. It is located in the Tepeaca Valley, the plain extends to the center of the poblana plateau and is characterized by its eminently limestone soil and marble deposits. The predominant soil types are Feozem and Calcisol (INEGI, 2017). The relief is flat, presenting a slight decline in a north-south direction. The predominant climate is temperate subhumid with rains in summer and very dry during the winter (INEGI, 2017). Its main economic activities are horticulture and agriculture of temporary.

In these municipalities the farmers use atrazine for weed control in corn crop, apply three pre-emergency irrigations two or three weeks before planting, during the months of march or April, incorporate spray atrazine in the first 3 to 5 cm deep, sometimes mixing it with fertilizers and other herbicides, such as carbofuram, without having an established dose. They do not use protective equipment or have controlled management for the preparation of solutions with herbicides. In addition, they rotate crops on an annual basis, sowing pumpkin, bean, carrot or broccoli (SENASICA, 2019).

Soil sampling

Two soil types were selected, one Regosol in San Nicolás Buenos Aires and the other Calcisol in Los Reyes de Juárez. The sampling units in each municipality corresponded to three agricultural plots planted with maize and without application of atrazine. Soil samples were taken in November 2018. From the first 30 cm deep, the soil sample was obtained following the methodology known as ‘five golds’ obtaining five subsamples of each plot.

The subsamples were placed in properly labeled plastic bags, indicating the sampling unit, they were subsequently transferred to the soil laboratory of the Department of Agricultural Science Research of the Meritorious Autonomous University of Puebla (BUAP), for their physicochemical characterization, according to the standardized methods of the Official Mexican Standard NOM-SEMARNAT-021-2000 (Official Mexican Standard, 2002). Given the homogeneity in the results obtained from the individual characterizations of the subsamples of each of the three plots of each municipality, the subsamples were mixed to generate a composite sample of each plot for two sampling units.

Characterization of soils without atrazine

For each soil the following parameters were determined: density, pH, electrical conductivity, texture, cation exchange capacity (CIC), concentrations of organic material, calcium and magnesium. The main minerals present in each soil were characterized by X-ray diffraction using a Bruker D 8m equipment with CuKα radiation, with scan ranges of 0.5° to 5° (2θ) with a step size of 0.025°. The possible presence of atrazine in the sampled soils was analyzed by performing methanol extractions from the soil. To do this, 5 g of each soil was used, which were placed in centrifugal tubes and added 25 mL of a 70/30 v/v methanol/water solution corresponding to a weight/volume ratio extractant-soil 5:1; subsequently, they were placed in an orbital agitator at 180 rpm for 2 h. From the supernatant, 5 ml were passed through 0.45 mm pore diameter nylon filters and stored in amber bottles for further analysis.

The concentration of atrazine residues was determined by a high-efficiency liquid chromatograph (HPLC) coupled to a dynamic equipment consisting of a model 1 525 high-pressure binary pump system, a 2 707 model self-sampler, a UV-Vis detector model 2 489 and an xSelect CSH analytical column (4.6 x 150 mm) of 3.5 µm particle size, all of them of the Waters brand. Acetonitrile HPLC grade was used as a mobile phase in isocratic mode with an extraction flow of 1 ml min-1 with a wavelength of 220 nm and samples were taken automatically every 5 min, the calibration curve was performed for the verification of this method. The chromatographic data obtained were stored in a data acquisition system using Waters Corp’s Empower 3.0 software. This methodology was based on the proposal by Salazar et al. (2018).

Columns packed with soil and study of atrazine mobility

The two types of soil were packed in random consecutive columns of transparent polychloride vinyl (PVC) approximately 8.5 cm internal diameter by 30 cm in height. Each column was packed to a height of 20 cm and its weight and density were determined. Field capacity was estimated using the Bodman and Mahmud formula, soil porosity, infiltration and leaching speed.

At the beginning of the experiment each column was added a restorative solution (0.001M KCl, 0.006M CaCl2, 0.001M MgCl2 and 0.018M NaCl) in continuous flow using a drip irrigation system to simulate the pH, electrical conductivity and ion force conditions of the soil solution (Duwig et al., 2009).

Once pore saturation had been reached, the flow of distilled water was restarted again steadily (2.62 ± 0.3 ml min-1). A 6.25 g L-1 atrazine (Gesaprim®) solution was applied by spraying at the top of the column on days 1, 11 and 21. To do this, a solution volume of 250 ml was used for the columns of Regosol and 200 mL for those of Calcisol. Drip irrigation was used every other day only with distilled water. The experiment was conducted during the months of august and September 2019, having an average temperature of 25 °C.

Because the inflow of water and outflow of the leachates during the experiment was different for each soil column, during the 30 days of the experiment, the percentage of leached water was calculated relative to the amount of volume applied. This compares the leaching process in each of the columns in a more homogeneous way.

Leachate recovery

The leachates were obtained daily from the packed columns for 30 days, the volume of the leached samples was recorded, 5 ml of each sample was filtered and preserved in amber bottles at 4 °C for further analysis. The quantification of atrazine was performed by HPLC, in the same way as in the analysis of methanol extractions (Salazar et al., 2018) .

Atrazine transport parameters in columns packed with soil

At the end of the 30 days of irrigation, the soil was unpacked from the columns previously treated with atrazine, each column was divided into three sections to the following depths: 0-2 cm, 2-10 cm, and 10-20 cm. Each section of the soil was homogeneously mixed, it was weighed and dried on a stove at 60 °C for 8 h, subsequently the extraction of atrazine was performed following the method proposed above. Finally, the delay coefficient was determined.

This coefficient provides a first approximation of the leaching potential of atrazine through the soil, taking into account the adsorption and distribution between the solid and liquid phases, is calculated by the following equation:  (1). Where: R= delay coefficient (dimensionless), ρ: apparent density (g cm-3), η: porosity (dimensionless), Kd: distribution coefficient (cm3 g-1) (Eweis et al., 1999). The mass balance for atrazine was estimated with the following equation (Raymundo et al., 2011):  (2). Where: MA is the applied atrazine mass, ML is the total leached mass, MR is the mass of atrazine retained in the soil of the column and MD is the degraded atrazine mass. By difference, with equation (2) the amount of degraded atrazine is obtained.

A constant rate of degradation (χ) was assumed with the following equation (Raymundo et al., 2011): . Where: t is the total time of the duration of the experiment. The half-life time of atrazine was estimated in both soil types with the following equation (Raymundo et al., 2011):

Statistical analysis of soil characterization results, field capacity, porosity, infiltration rate and infiltration flow rates were performed by a student t-test for independent samples with significance level p≤ 0.05. The data obtained for the retention of atrazine and its transport parameters were statistically analyzed using a variance analysis (Anova) and the study of differences between means per Tukey test. The relationship between the properties of the soils studied and the mobility of atrazine was established through Pearson's linear correlation method. In all cases a significance value of p≤ 0.05 was considered, all analyses were performed in R statistical program version 3.1.3.

Results and discussion

Characterization of soils without atrazine

Table 1 shows the results obtained on the characterization of the soils evaluated. The regosol soil of the municipality of San Nicolás Buenos Aires is moderately acidic (5.8). It does not have high salinity (Castellanos, 2000)  as it shows a relatively low electrical conductivity value (0.142 dS m-1). This soil is of a loam-sandy texture class with 80% sands, with medium concentration organic matter (1.95%). This is because in this region the soil is under constant labor and is cultivated devoid of plant cover that could provide a greater amount of organic material. Cation exchange capacity (CIC) is low (8 cmol+ kg-1), therefore, the soil’s ability to

retain cations and nutrients for maize cultivation (Castellanos, 2000) is limited. Concentrations of interchangeable Ca2+ and Mg2+ ions are medium (4.9 meq 100 g-1 soil) and low (1.16 meq 100 g-1 soil), respectively.

The Calcisol-type soil in the municipality of Reyes de Juárez, presented a pH close to neutrality (7.1); its electrical conductivity is very low (0.25 dS m-1), so they are not considered saline (Castellanos, 2000). This type of soil has a loam textural class with predominance of sands (47%) and clays (36%). The concentration of organic material is 1.84 %, value considered as medium (Castellanos, 2000). Its CIC is medium (18.9 Cmol kg-1)  (Official Mexican Standard, 2002),due to the presence and nature of the clays. It has high concentrations of interchangeable ions of Ca2+ and Mg2+, 9.7 and 3.2 meq 100 g-1 soil, respectively (Castellanos, 2000).

Table 1. Physicochemical properties of the evaluated soils.

Means with different letter are statistically different between soil types (p≤ 0.05).

As for mineralogy, the Regosol presented an X-ray diffractogram (Figure 2) indicating the presence of albite and quartz, as reported by Lebgue-Keleng et al. (2014) for Regosol from the state of Chihuahua. These crystalline minerals do not provide significant contributions of nutrients to the soil  (Sadeghian et al., 2012). Regarding Calcisol, the X-ray diffractogram indicated that it is rich in CaCO3, (Sadeghian et al., 2012), although the presence of albita is also observed.

Weber (1995)  reported that atrazine adsorption is pH-independent when it has values greater than 4 and Koskinen and Clay (1997) reported that the maximum pH value affecting atrazine adsorption is 6. Therefore, in the soils studied there is no influence of pH on the adsorption of atrazine because the pH of the soils used are 5.8 and 7.1. Organic material concentrations are relatively low in both soils without significant differences (p≤ 0.05). According to the studies of Wan Ting et al. (2005), this property can negatively affect the adsorption of atrazine in Regosol soil, while for Calcisol, which has significant difference (p≤ 0.05) and high clay content, atrazine adsorption would be favored, but biodegradation processes would be limited.

Figure 2. X-ray diffractograms of soil samples at large angles (Ca calcite, Al albite, Q quartz).

The presence of atrazine was not observed before starting the treatments in the sampled soils when performing high-efficiency liquid chromatography studies of the extracts from each of the soils, because the peak corresponding to atrazine was not observed.

Atrazine leaching in soil columns

The value obtained for field capacity is consistent with those reported by Saxton and Rawls (2006) for the evaluated soils, Regosol and Calcisol. Porosity presented statistically significant difference (p≤ 0.05) for soil types, being for Regosol of 32% and for Calcisol of 50%, this is a variable that determines the infiltration and leaching processes that influence the transport of water and agrochemicals (Horowitz and Walling, 2005). The volume of water in the pores for Calcisol (254 cm3) and for Regosol (156 cm3) indicated the relationship with porosity, which in both soils was favorable for crop development (Volverás-Mambuscay et al., 2016) . These results can be explained from the perspective of agronomic management of plots and the temporal effect of farm implements, which improve soil properties such as porosity and moisture retention capacity, in soils such as Calcisol, other studies report similar results in tropical agricultural systems (Sustaita et al., 2000)

Infiltration flow and leaching flow, Table 2, present significant statistical difference (p≤ 0.05) for soil types presenting higher values for Regosol, due to the formation of a greater number of preferential flow routes, characteristic of soils with a higher percentage of sand, which capture some of the surface flow and transport it into the soil, facilitating the mobility of water and contaminants (Shakesby, 2000).

Table 2. Leaching characteristics of packaged columns.

m= soil mass; CC= field capacity; ϕ= porosity; υ= volume of water in pores; QI= infiltration flow; QL= leaching flow. Means with different letter are statistically different between soil types (p≤ 0.05).

Elution curves (Figure 3) show the concentration of atrazine in leachates in parts per million (ppm) versus the daily pore volume of 200 ml for Regosol and 180 ml for Calcisol. After 5 volumes of pore, Calcisol presented the concentration of atrazine accumulated in leachates of 99.3 ppm. For the case of the regosol of San Nicolás Buenos Aires, a cumulative concentration of atrazine was obtained in leachates of 114.2 ppm after five volumes of pore, being statistically different (p≤ 0.05).

It is important to note that, even though the atrazine solution was also added on days 11 and 21, the presence of herbicide was not observed in the leachates of both soils. Soils with a low organic matter content have been reported to be more susceptible to leaching (DeSutter  et al., 2003) . However, the leaching flow rate and infiltration flow rate is 9 and 21 times, respectively larger for the Regosol soil, resulting in a higher leaching speed.

Figure 3.  Relative elution curves for the soils studied.

Comparing the atrazine concentration results obtained in this study with those reported by Ndongo et al. (2000) it can be mentioned that the high concentration of atrazine in leachate in the first days is related to the presence of preferential paths in soil pores. The amount of leachate volume decreased 50% after 21 pore volumes for Regosol and 22 pore volumes for Calcisol, due to pore saturation and moisture accumulation within columns (Figure 4).

Figure 4. Percentage of leached irrigation water in relation to the amount of pore volume applied.

Salazar-Ledesma et al. (2018) assessed the adsorption-desorption capacity of atrazine in three soil types under different agricultural management, concluding that the elution of the herbicide depends on the content of organic matter, its chemical composition and the agricultural management of each soil, therefore, having a similar organic matter content, both types of soil contain a similar concentration of atrazine in their leachates (Salazar et al., 2018).

Atrazine retained in soil columns

Table 3 reports the values of atrazine per gram of soil retained at different depths of the columns of each soil type. It is observed that the herbicide is at higher concentration in the first two centimeters of the column, 2.4 and 2.3 mg of atrazine g-1 soil for Regosol and Calcisol, respectively, without presenting significant difference (p≤ 0.05). For Regosol, atrazine had a lower retention inside the column, because sandy soils facilitate leaching, have no retention and have high permeability (Sadeghian and Arias, 2017). In this soil, the function of atrazine in weed control is guaranteed, because it remains on the surface longer. For Calcisol there is migration to a deeper surface of the column (1.8 mg atrazine g-1 soil) in this case, there is a high risk of contamination of surface and groundwater.

Table 3. Atrazine retained at different depths of the columns.

Means with different letter are statistically different between soil types (p≤ 0.05).

Atrazine transport parameters

Table 4 presents the values of the distribution coefficient (Kd) of atrazine, without presenting statistical difference (p≤ 0.05), 1.99 for Regosol and 1.88 for Calcisol. The value of the delay coefficient (R) indicated that the atrazine is mobile in both soil types (Bernard et al., 2005), presenting statistical difference (p≤ 0.05),  being the Regosol the largest R= 1.08 ±0.18 indicator of faster leaching processes in this type of soils and groundwater contamination (Gutiérrez et al., 2007). The importance of calculating the delay coefficient (R) of atrazine in the columns is due to the advection and adsorption effects, considering moisture and porosity in a soil type. Therefore, the conditions of the hydraulic behavior of the soil would be similar to those used in real conditions (Eweis et al.,  1999).

Finally, hydraulic conductivity was significantly higher (p≤ 0.05) for Regosol (1.91 ml min-1) than for Calcisol (0.09 ml min-1), indicating a greater ability of Regosol to transmit water to the soil profile. In addition, the percentage of degraded mass was 94% for Regosol and 69% for Calcisol. The rate of degradation was significantly higher (p≤ 0.05) in Regosol (0.097 days) with respect to Calcisol (0.039 days), therefore, atrazine persists longer in Calcisol, with half-life time of 17.77 days.

Table 4. Atrazine transport parameters.

Kd= foc*Koc; Foc= fraction of organic carbon available in soil; Koc= organic carbon partition coefficient. For the atrazine Koc= 102.32 ml g-1; χ= constant rate of degradation; t0.5= half-life time. Means with different letter are statistically different between soil types (p≤ 0.05).

The results of the correlation analysis between atrazine mobility and soil properties are presented in Table 5. The correlation coefficients obtained indicate that the factors that affect the mobility of atrazine in the soil are the content of clays, organic matter, texture and cationic exchange capacity, which is consistent with Eweis et al. (1999). Soils with higher clay and organic matter content have greater potential to retain herbicides, so their availability decreases in the soil solution to be absorbed by plants.

Table 5. Pearson correlation matrix between atrazine mobility and soil properties.

Conclusions

The characterization of Regosol and Calcisol showed low organic material concentrations (1.95% and 1.84%, respectively), the concentration of clays and cationic exchange capacity were significantly different between the two soils, being higher in Regosol with values of 36% and 18.9 cmol+ kg-1. These properties showed high correlation (>0.8) with the mobility of atrazine in soil columns, so the delay coefficient was 1.08 in Regosol and 1.03 in Calcisol. Because the hydraulic conductivity atrazine was higher in Regosol (1.91 ml min-1) is that it was transported nine times faster than Calcisol (0.09 ml min-1). Atrazine was maintained in greater quantity in the first two centimeters of soil columns, 2.4 and 2.3 mg atrazine g-1 soil for Regosol and Calcisol and the percentage of degraded mass of atrazine was 94% for regosol and 69% for Calcisol.

Cited literature

Abdelhafid, R.; Houot, S. and Barriuso, E. 2000. How increasing availabilities of carbon and nitrogen affect atrazine behaviour in soils. Biology and Fertility of Soils. 30(4):333-340.

Azevedo, A. S.; Kanwarbl, R. S. and Pereira, S. 2000. Atrazine transport in irrigated heavy-and coarse-textured soils, part i: field studies. J. Agric. Eng. Res. 76(2):165-174.

Bernard, H.; Chabalier, P. F.; Chopart, J. L.; Legube, B. and Vauclin, M. 2005. Assessment of herbicide leaching risk in two tropical soils of Reunion Island (France). J. Environ. Quality. 34(2):534-543.

Carrero, J. y Planes, S. 2008. Plagas del campo (Mundi-Prensa (Ed.). 224-230 pp. 

Castellanos, J. Z.; Bueno, U. J. X. y Aguilar, S. A. 2000. Manual de interpretación de análisis de suelos y aguas. Instituto de Capacitación para la Productividad Agrícola. 2da. edición. 198-212 pp.

DeSutter, T. H.; Clay, S. A. and Clay, D. E. 2003. Atrazine sorption and desorption as affected by aggregate size, particle size, and soil type. Weed Sci. 51(3):456-462.

Duwig, C.; Prado-Pano, B. L.; Hidalgo-Moreno, C. I; Gavi-Reyes, F.; Raymundo, R. E.; Nikolskii, G. I. y Figueroa, S. B. 2009. Transporte de atrazina en un Andosol y un Vertisol de México. 34(5):330-338.

Eweis, J. B.; Ergas, S. J.; Chang, D. P. y Schroeder, E. D. 1999. Principios de biorrecuperación, (Ed. Mc Graw-Hill Interamericana. 327 p.

Fernández, M.; González, E.; Villafranca, M. and Flores, F. 2001. Mobility of atrazine from alginate-bentonite controlled release formulations in layered soil. Chemosphere. 43(3):347-353.

Graymore, M.; Stagnitti, F. and Allinson, G. 2001. Impacts of atrazine in aquatic ecosystems. Environ. Inter. 26(7-8):483-495.

Hang, S.; Andriulo, A.; Sasal, C.; Nassetta, M. M.; Portela, S. and Cañas, A. I. 2010. Integral study of atrazine behavior in field lysimeters in Argentinean humid pampas soils. Chilean J. Agric. Res. 70(1):104-112.

Hang, S. y Sereno, R. 2002. Adsorción de atrazina y su relación con las características sedimentológicas y el desarrollo del perfil de dos suelos de la provincia de Córdoba. Rev. Investigaciones Agropec. 31(3):73-87.

Hernández, A. A. y Hansen, A. M. 2011. Uso de plaguicidas en dos zonas agrícolas de México y evaluación de la contaminación de agua y sedimentos. Rev. Inter. Contaminación Amb. 27(2):115-127.

INEGI. 2017. Anuario Estadístico del Estado de Puebla.

Javanshir, A. M.; Hadi, M.; Abdolazim, B.; Saeed-Boroomand, N. and Zahra-Nazari, K. 2012. Atrazine transport and distribution in field soils and comparison of the predictions made by leaching estimation and chemistry model-pesticide (LEACHP) model. Sci. Res. Essays. 7(40):3348-3356.

Josep, M. y Figueras, L. 1998. Acuíferos de Cataluña. 127-133 pp.

Kim, S. H.; Fan, M.; Prasher, S. O.; Patel, R. M. and Hussain, S. A. 2011. Fate and transport of atrazine in a sandy soil in the presence of antibiotics in poultry manures. Agric. Water Management. 98(4):653-660.

Lebgue-Keleng, T.; Soto-Cruz, R. A.; Melgoza-Castillos, A. y Morales-Nieto, C. 2014. Mineralogía del suelo y su relación con las especies de cactáceas en el sureste de Chihuahua, México. Ecología Aplicada. 13(2):205-209.

Lerch, R. N.; Sadler, E. J.; Sudduth, K. A.; Baffaut, C. and Kitchen, N. R. 2011. Herbicide transport in goodwater creek experimental watershed: I. long-term research on Atrazine. J. Am. Water Res. Association. 47(2):209-223.

López-Geta, J. A.; Martínez-Navarrete, C.; Moreno-Merino, L. y Navarrete Martínez, P. 1992. Las aguas subterráneas y los plaguicidas. Instituto Geológico y Minero de España. 149 p.

González-Márquez, L. C. y Hansen, A. M. 2009. Adsorción y mineralización de atrazina y relación con parámetros de suelos del DR 063 Guasave, Sinaloa. Rev. Mex. Cienc. Geol. 26(3):587-599.

Gutiérrez, H. M.; Barba, L. E. and Materón, H. 2007. Mobility of pesticides carbofurane and imidacloprid in Typic Humitropept soil. Agron. Colomb. 25(1):160-167.

Montoya, J.; Costa, J.; Liedi, R.; Bedmar, F. and Daniel, P. 2006. Effects of soil type and tillage practice on atrazine transport through intact soil cores. Geoderma. 137(1-2):161-173.

Mudhoo, A. and Garg, V. K. 2011. Sorption, transport and transformation of Atrazine in soils, minerals and composts: a review. Pedosphere. 21(1):11-25.

Ndongo, B.; Leroux, G. and Fortín, J. 2000. Transport du linuron, d´el imidaclopride et du bromure au travers de colonnes de sol et de lysimètres drainast. Biotechnology Agronomic Society. 4(1):33-40.

Norma Oficial Mexicana. (2002). NOM-021-SEMARNAT-2000 Especificaciones de fertilidad, salinidad y clasificación de suelos, estudio, muestreo y análisis. Diario Oficial de La Federación. 85 p.

Raymundo, R. E.; Prado, P. B. L.; Nikolskii, G. I. y Mendoza, H. J. H. R. 2011. Sorción, degradación y lixiviación de imazapyr en un suelo volcánico de México. Tecnología y Ciencias del Agua. 2(3):61-75.

Sadeghian, K. S. y Zapata, H. R. D. 2012. Propiedades relacionadas con la adsorción de cationes intercambiables en algunos suelos de la zona cafetera de Colombia. Cenicafé. 63(2):79-89.

Sadeghian, K. S. y Arias, S. E. 2017, Lixiviación del potasio en suelos de la zona cafetera y su relación con la textura. Biotecnología en el Sector Agropecuario y Agroindustrial. 16(1):34-40.

Salazar-Ledesma, M.; Prado, B.; Zamora, O. and Siebe, C. 2018. Mobility of atrazine in soils of a wastewater irrigated maize field. Agric. Ecosys. Environ. 255:73-83.

Saxton, K. E. and Rawls, W. J. 2006. Soil water characteristic estimates by texture and organic matter for hydrologic solutions. Soil Sci. Society Am. J. 70(5):1569-1578.

Shakesby, R. A.; Doerr, S. H. and Walsh, R. P. D. 2000. The erosional impact of soil hydrophobicity: current problems and future research directions. J. Hydrol. 231-232:178-191.

Shukla, M. K.; Ellsworth, T. R.; Hudson, R. J. and Nielsen, D. R. 2003. Effect of water flux on solute velocity and dispersion. Soil Sci. Soc. Ame. J. 67(2):449-457.

SIAP. 2019. Cierre de la producción agrícola por estado. Puebla.

Villada-Canela, M. 2006. Estimación del riesgo de contaminar mantos freáticos por infiltración de herbicida atrazina en distritos de riego de México. 150 p.

Volverás-Mambuscay, B.; Amézquita, C. É. y Campo, Q. J. M. 2016. Recursos biofísicos. Corpoica. 17(3):361-377.

Wan Ting, L.; Hai-Zhen W.; Jian-Ming X. and Yan-Zheng, G. 2005. Sorption of dissolved organic matter and its effects on the atrazine sorption on soils. J. Environ Sci. 17(3):478-482.

Weber, J. 1995. Physicochemical and mobility studies with pesticides. Chapter 10, in Agrochemical environmental fate. 99-112 pp.