Equilibrium, kinetic and thermodynamic studies on the adsorption of atrazine in soils of the water fluctuation zone in the Three-Gorges Reservoir

Environmental behavior of pesticide in soils is a key current research focus. Studying the adsorption characteristics of pesticides in soils as a parameter for evaluating the risk of groundwater pollution by pesticides is commonly applied in agriculture. To provide a theoretical basis for environment risk assessment and pollution remediation, the thermodynamics and kinetics of the equilibrium of atrazine adsorption in the Three-Gorges Reservoir area were assessed and analyzed via batch experiments. Results showed that the sorption of atrazine was an exothermic and spontaneous process at temperatures of 298–318 K. Atrazine was more easily adsorbed by soils at concentrations of 0–30 mg L−1, with low-temperature adsorption being more effective than high-temperature adsorption. The adsorption of atrazine to the two assessed soils was well fitted by the Freundlich and Langmuir models. The adsorption kinetics of atrazine in soils were consistent with the quasi-second-order kinetic model and intraparticle diffusion was found not to be the only control step. The monolayer adsorption occurred with non-uniform energy distribution on the soil surface, indicating that the adsorption of atrazine by the two kinds of soil was controlled by internal diffusion surface adsorption and liquid film diffusion, leading to the complexity of its adsorption kinetics. The values of standard free energy ΔrGm0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta rG_{m}^{0}$$\end{document} < 0, indicated that the adsorption of atrazine in soils was spontaneous and dominated by physical adsorption. Changes in standard enthalpy (ΔrHm0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {\text{r}}H_{m}^{0}$$\end{document}) indicated that the adsorption was exothermic. Atrazine exhibited a weak adsorption capacity in both soils, indicating it is highly mobile in the soil–water environment and can easily cause groundwater pollution. Therefore, much attention should be paid to the environmental behavior of pesticide soil moisture fluctuations, potentially leading to the broad transfer and spread of pollution.


Background
Pesticides are widely used in agriculture to increase crop yields and quality [1]. Atrazine is one of the most widely used pesticides worldwide, being a long-term and large-area herbicide suitable for the removal of broadleaved weed species [2][3][4][5]. Long-term atrazine residue accumulation in the soil accounts for 20%-70% of the applied dose during application and is one of the most commonly detected pesticides in soils and groundwater worldwide [6][7][8]. However, atrazine is one of the triazine herbicides which was banned by the European Union in 2004 (2004/248/EC) [9] and was added to the list of Open Access *Correspondence: lmma@tongji.edu.cn 1 College of Environmental Science and Engineering, Tongji University, Shanghai 200092, China Full list of author information is available at the end of the article priority substances (Directive 2013/39/EU) [10]. In recent years, atrazine and its decomposition intermediates have been regularly detected in groundwaters, lakes and bays in many countries [5,11,12]. Due to the persistence, high-water solubility and hazard associated with atrazine pollution, its large-scale and excessive use in agriculture presents a serious risk to both environmental and human health [12][13][14]. It has been reported that atrazine affects human immune, endocrine and reproductive systems; while in severe cases, it may lead to malformations and induce mutations, which are potentially carcinogenic to humans [4,5,15,16]. Therefore, atrazine's environmental behavior, particularly its transportation and fate, has received extensive research interest [7,16,17].
Atrazine undergoes a series of processes after entering the soil, such as soil adsorption-desorption, degradation, and crop absorption due to surface runoff and eluviation in surface and groundwater systems. Among these, adsorption-desorption is one of the key factors affecting atrazine's migration and transformation in the environment and, thus, directly affects the efficiency of atrazine pollution prevention and remediation strategies [18,19]. The adsorption behavior of atrazine in soil is also a key factor affecting the distribution, migration, bioavailability and toxicity of pollutants in soil, which has an important impact on the fate of pollutants in the environment [20][21][22][23]. Thus, the adsorption mechanism and influencing factors between pesticides and different soils have always been the focus of environmental studies [24]. There are many available studies on the adsorption, leaching and migration of atrazine in soils of different texture types, including paddy soil, black soil and loess soil, as well as the sorption and desorption processes [21,[25][26][27][28]. Studies have shown that the adsorption process of organic compounds is complex and affected by many factors, including pH, cation exchange capacity, surface area of the sorbent and ionic strength [29,30]. However, the adsorption kinetics and thermodynamic properties of atrazine in different soils remain unclear and further research is needed especially on the adsorption mechanisms of atrazine in soils of water fluctuation zones. The adsorption of atrazine in soils in different environments has been shown to vary significantly [30,31]. Therefore, it is important to study the adsorption behavior of atrazine in different soils, to analyze and evaluate the fate and risk of environmental migration and provide a reference for the treatment and control of atrazine pollution in soil [21,30].
The water fluctuation zone of a large reservoir is a combined land and water ecosystem, which varies within ecosystems in a cycle of alternating irrigation and drainage and is strongly disturbed by human activities [32,33]. The Three-Gorges Reservoir in China is a representative example with a vertical depth of 30 m and a total area of over 300 km 2 [34]. The water fluctuation zone of the Three-Gorges Reservoir is a traditional crop planting area. Agricultural production activities use a large amount of pesticides to increase production, which subsequently transfer to the soil. The absorbed pesticides in the soil are desorbed into the reservoir water body and serious environmental pollution caused by the release of pollutants from the soil becomes inevitable. Therefore, the purpose of this study was to investigate the adsorption mechanism and thermodynamics of atrazine in soils with a high intensity of anthropogenic activity, assessing the effect of parameters such as temperature and contact time. These results will help to better understand the dynamics and thermodynamics of atrazine adsorption processes, while providing methods for protecting drinking water resources from contamination and predicting the fate, transport and risk assessment of atrazine in the environment.

Soils
The assessed soil samples were collected from the bank of the Yangzi River from Wanzhou to Wushan in Chongqing, within the water fluctuation zone of the Three-Gorges Reservoir. The collected fresh soil samples were taken from the layer at a depth of 0-20 cm, naturally dried at room temperature, homogenized, with the removal of stones, plant roots and other non-soil debris, before being passed through a 2-mm mesh sieve for later use. The pH value was measured using a pH meter. Organic matter was treated by external heating with potassium dichromate. Soil particle size composition was determined via a Laser particle size analyzer (NKT6100, Shandong Nikeite Analytical Instrument Co., Ltd., Shandong, China). The cation exchange capacity (CEC) in soil was determined by cobalt extraction with hexaammonia trichloride by spectrophotometry (HJ 889-2017). The specific methods applied have previously been reported and all analyses were performed according to the previously reported methods without modification [35]. The measured properties of soils are outlined in Table 1.

Chemicals
The atrazine standard had a purity of > 99.9% (Dr. Ehrenstofer company; Germany), with the properties of atrazine outlined in Table 2. Methanol and acetonitrile were of chromatographic purity (China pharmaceutical group chemical reagent Co. Ltd., Shanghai, China). All other reagents were of analytical purity (ANPEL Scientific Instrument, Shanghai, Co., Ltd.). Ultrapure water was prepared using a Pulfil DZG-303A ultra-pure water system (Fusite instrument equipment co., Ltd., Shanghai).

Sorption experiments
Quantification of the sorption of atrazine in the two different soils was performed following the standard batch equilibrium method. Atrazine was weighed and dissolved in methanol, with a mole fraction of less than 0.001 used to prepare the pesticide solution. Studies have previously reported that methanol has no effect on adsorption [36][37][38][39]. Experiments were carried out in 50-mL polyethylene centrifuge tubes by adding 2.0 g of soil samples to 10-mL solutions containing different atrazine concentrations (0-30 mg L −1 ). During experiments, 0.01 mol calcium chloride was used as the background electrolyte of the soil solution, to minimize the change in ionic strength and contribute to the separation of solid and aqueous solutions after equilibrium. The centrifuge tube was shaken for 48 h at a constant speed of 200 r min −1 in a thermostatic vibrating screen, with the adsorption equilibrium reached at temperatures of 298, 308, and 318 K, respectively, and then centrifuged for 30 min at the same temperature. The same method was used for the blank experiment, without atrazine being added. Preliminary experiments have shown that the adsorption of pesticides to test tubes was negligible. The supernatant was collected and passed through a 0.45-um micron filter; then, the atrazine content was measured by HPLC. The equilibrium adsorption capacity q e (mg kg −1 ) of atrazine was calculated as follows in Eq. (1): where V represents the volume of solutions (L); C 0 is the initial concentration of atrazine; C e (mg L −1 ) is the equilibrium concentration of the adsorbate at time t (min); m is the mass of adsorbent (g); q e (mg kg −1 ) is the amount of atrazine adsorbed on soil.
All kinetic experiments were carried out at 25 °C under normal ambient pressures. A 2.0-g soil sample was placed in a 50-mL bottle containing 10-mL atrazine solution at a concentration of 30 mg L −1 and stirred using a constant temperature shaker. Samples were taken at set intervals and three repetitions were performed for each treatment. The supernatant was centrifuged at 8000 rpm for 5 min and then passed through a 0.45-µm micron filter for later use.

HPLC analysis
The concentrations of atrazine in soil samples were analyzed by HPLC (Agilent 1200), using an ACE Excle 5 C 18 analytical column (5 μm, 250 × 4.6 mm, Phenomenex, Guangzhou, China). HPLC analyses were performed using mobile phase elution with 40% (v/v) ultrapure water and 60% (v/v) acetonitrile and an ultraviolet detector to monitor the absorbance of atrazine at 220 nm.

Adsorption isotherms
To understand the interactions between atrazine and soil, as well as revealing the adsorption mechanism, adsorption isotherms were quantitatively established to analyze the process of atrazine transfer from the solid phase to the liquid phase [40]. The adsorption isotherm of atrazine in both soils was depicted at different temperatures, as shown in Fig. 1. Results revealed that the adsorption capacity of atrazine was significantly different in the two soils, which might be attributed to the difference in physical and chemical properties between the two soils [21,29,30]. The adsorption curves for atrazine in the two soils were both L-type, with absorption increasing sharply initially and then showing a gradually increasing trend at 288, 298 and 308 K (Fig. 1a, b). For each isotherm, the soil equilibrium adsorption capacity increased with increasing equilibrium concentration, due to the adsorption driving force increasing with the enhancement of equilibrium concentration [41]. The adsorption capacity of soil to atrazine decreased with the increase  in temperature within an initial atrazine concentration range from 0 to 30 mg L −1 , especially in the Wushan soil (Fig. 1b). In addition, the equilibrium adsorption capacity also decreased with increasing temperature, suggesting that the atrazine adsorption process in the two soils was an exothermic process [26]. Generally, the adsorption capacity of atrazine in both soils was not very large, demonstrating that atrazine can easily penetrate the soil vadose zone and pollute groundwater. The adsorption mechanism of atrazine in soil is likely to be related to soil properties. Monolayer sorption and heterogeneous energetic distribution of active sites on the surface of the soils were possible during the sorption process [41]. In the assessed range of atrazine concentrations, soil adsorption of atrazine exhibited a basic linear relationship, indicating that atrazine mainly adsorbs to soil organic matter via hydrophobic effects [42]. This type of distribution function is mainly related to the composition of soil organic matter and soil characteristics: (1) hydrogen bonds are formed between the neutral atrazine molecules and active centers on the surface of soil colloids; (2) ionic atrazine adsorbs on colloidal surfaces such as soil organic polyelectrolyte clay minerals; and (3) atrazine can form complexes with carboxyl groups and free hydroxyl groups in soil humic acid on the surface of soil colloidal micelles [26,43]. Therefore, adsorption models have been further used to analyze the mechanism of adsorption.
Linear, Friedrich or Langmuir adsorption models have commonly been used to describe analytical adsorption isotherms [44]. To clarify the overall adsorption process of atrazine in the two soils, Langmuir and Freundlich sorption models have been used. The Langmuir isotherm model [45] assumes that the energy adsorbed on the surface of the adsorbent is uniform, the adsorbent does not migrate in the surface plane and there is no interaction between the adsorbed molecules [46,47], the equation for which is described as follows in Eq. (2): where q e is the amount of atrazine adsorbed by soil (mg kg −1 ); q m is the maximum amount of atrazine adsorbed (mg kg −1 ); K L (L mg −1 ) is a constant related to the adsorption energy; and C e is the equilibrium concentration of atrazine (mg L −1 ). The Freundlich isotherm model is used to describe heterogeneous surface equations such as the heterogeneity of the adsorbent surface, the adsorption energy and the exponential distribution of the adsorption point. The adsorption data fit the Freundlich isotherm [48] and its expression is described in Eq. (3) as follows: where K F and n are the Freundlich constants that refer to relative capacity and adsorption intensity, respectively. The Freundlich and Langmuir models were used to fit the adsorption data of atrazine in soil, with the parameters depicted in Table 3. In all cases, the correlation coefficient (R 2 ) of the Langmuir and Freundlich isotherm models exceeded 0.99, indicating that both models fit well to the adsorption of atrazine in both soil samples. According to the Freundlich model adsorption constant (K F ), the adsorption capacity of the Wanzhou soil to atrazine was higher than Wushan soil and according to the Langmuir model, the sorption capacities ( q m ) of Wanzhou soil were higher than that of Wushan soil, which is consistent with the change in K F . These results also demonstrate that the adsorption mechanism of atrazine in soil is related to both physical and chemical properties of the soil [18,41].
The q m of the atrazine adsorbed by Wanzhou soil and Wushan soil decreased with increasing temperature, suggesting that the adsorption amount was negatively correlated with temperature. For the Freundlich model, n < 1 represents unfavorable absorption conditions and n > 1 represents favorable sorption conditions. In this study, the lowest n values for atrazine adsorbed by Wushan soil at the three temperatures of 298, 308, and 318 K were 1.49, 1.13, and 1.07, respectively. The n values were greater than 1 in all experiments, suggesting that sorption in both soil samples were favorable. The atrazine K F value for soil adsorption decreased with increasing temperature from 298 to 318 K. Moreover, the values of K F were less than 5 in all experiments, implying that atrazine has a weak adsorption capacity in these two soils, that atrazine is highly mobile in the soil-water environment and likely to cause groundwater pollution [27,49].
Based on further analysis of the Langmuir equation, the dimensionless parameter of the equilibrium or sorption intensity ( R L ) was calculated using Eq. (4) as follows [50,51]: where C 0 is the initial concentration of atrazine (mg L −1 ). The change in soil adsorption strength (R L ) indicates the type of isotherm and adsorption state, as per the following classification criteria: R L > 1 indicates unfavorable adsorption; R L = 1 indicates linear adsorption; 0 < R L < 1 indicates favorable adsorption; R L = 0 indicates irreversible adsorption [50,52].
In this study, the R L values for atrazine adsorption by Wanzhou soils were 0.37, 0.47 and 0.45; while for Wushan soils, they were 0.33, 0.31, and 0.17 at 298, 308, and 318 K, respectively. The R L values for atrazine indicate that it could be well adsorbed by the two assessed (4) soils. With increasing temperatures, the reaction rate decreased, indicating that the adsorption of atrazine in low temperature environments was easier than in high-temperature environments. Results imply that the two soils have a certain affinity for atrazine, but absorption was negatively correlated with temperature [49]. In addition, low-temperature conditions can promote the adsorption of atrazine in both soil samples, which is consistent with the results of the Freundlich model.

Adsorption kinetics
Adsorption kinetics are one of the main properties controlling the absorption rate and adsorption efficiency of solutes and are highly valuable when revealing the migration and transformation rate of atrazine in soil environments [27]. In this study, we used adsorption kinetics to assess the effects and dynamics of contact time on the adsorption capacity of atrazine in the two soil samples as shown in Fig. 2. As shown, the adsorption kinetics of atrazine in the two soils were similar, with both undergoing distinct stages of rapid adsorption, slow adsorption and equilibrium adsorption. Within 60 min, the adsorption rate of atrazine in the both soils was increased and  the adsorption capacity showed a sharp increase because of physical adsorption to mineral surfaces [49,53]. The faster initial sorption rate might be attributed to chemical bonds and hydrogen bonds between the atrazine and minerals [24]. With an increase in contact time, the adsorption rate decreased and the adsorption process gradually approached equilibrium within 48 h. The reason for this phenomenon may be due to the high initial atrazine concentration increasing the molecular collision frequency, allowing the dissolved state atrazine molecules to rapidly adsorb onto the surface of soil particles. In contrast, the molecular collision frequency decreased with low initial atrazine concentrations, resulting in a slower adsorption rate [26]. In addition, the rapid adsorption phase of atrazine to soil was determined by the large pores on the surface of particles; while, the slow adsorption phase was determined by micropores inside the soil particles. Slow adsorption processes may be related to the region where organic compounds gradually enter soil micropores or where soil organic matter is highly crosslinked, covering the adsorption sites on the internal surfaces of minerals [54]. This phenomenon is compatible with the reported results of adsorption reactions of organic pollutants on porous media particles, showing the reaction can generally be divided into a rapid reaction phase and a slow reaction phase [55]. The adsorption process of atrazine to soils is complicated, as a result of the interactions between organic matter and inorganic minerals in soil [24,56]. The rapid adsorption of organic pollutants can be attributed to their distribution in soil organic matter and on the surface of minerals. Physical adsorption in terms of the intermolecular interaction force is mainly manifested by the van der Waals force, dipole force, and hydrogen bonding force, with these effects usually completed in a relatively short time period [43]. To reach the adsorption site on the surface of soil particles, hydrophobic organic pollutants need to overcome the diffusion of water molecules on the surface of soil particles into soil micropores and diffuse into the pores of the soil particles through the matrix [43,57]. The rapid adsorption process of atrazine mainly occurs due to the presence of organic matter contained in the soil, as well as the distribution and physical adsorption of the mineral surface [43].
To describe the characteristics and mechanism of atrazine adsorption in soil, a pseudo-second-order kinetic model was employed, which can be described according to Eq. (5) as follows [58]: The integration of the Eq. (5) with the initial condition, q t = 0 at t = 0 leads to linear Eq. (6) as follows: where q t is the adsorption capacity of atrazine at time t (mg kg −1 ); k 2 is the pseudo-second-order reaction rate constant (kg (mg min) −1 ); q e is the amount of atrazine (mg kg −1 ) adsorbed by soils at equilibrium. The linearized pseudo-second-order kinetic model of atrazine adsorption to soil and the assessed parameters are presented in Fig. 3 and Table 4. According to the fitting correlation coefficient (R 2 ), the pseudo-second-order kinetic model exhibited good consistency with experimental data and could describe the dynamic processes of atrazine adsorption in both soils, with R 2 values almost equal to 1 for both soil samples. This indicates that the adsorption behavior of atrazine is obvious and that the rate during the adsorption process is controlled by chemical processes [58,59]. The pseudo-second-order kinetic model covers all processes of adsorption, such as external liquid film diffusion, internal particle diffusion and surface adsorption, among others [42], which can reflect the adsorption mechanism of atrazine in soil more accurately and comprehensively.  Generally, the adsorption kinetics of particles or molecules onto the surface of soils mainly includes two phases of transport and attachment separation [59][60][61]. The transport of particles is based on diffusion, which is caused by the attraction of the surface to the adsorbate [62]. However, the pseudo-second-order kinetic model does not identify the diffusion mechanism of atrazine soil adsorption.
Therefore, to better understand the factors affecting adsorption kinetics, the Weber-Morris intraparticle diffusion model [39] was applied to simulate experimental data, which was calculated according to Eq. (7) as follows: where k id (kg (mg min 0.5 ) −1 ) is the intraparticle diffusion rate constant; and c is a constant.
The linearized Weber-Morris intra-particle diffusion kinetic model and its parameters for atrazine sorption in soils are depicted in Table 3. According to previous studies, q t and t 0.5 in the interparticle diffusion model exhibit a linear relationship. Based on the intercept through the origin point, results indicated that the interparticle diffusion process of substances is not the only control step for adsorption rate, with other mechanisms potentially playing an important role [41,59,63,64]. From Fig. 4 and parameters in Table 3, it was found that the linear relationship between q t and t 0.5 was not reasonable and that plots do not pass through the point of origin. This indicates that the intragranular diffusion process of atrazine was not the main control process of adsorption in both soils. The adsorption process involves a certain degree of boundary layer control and the adsorption rate is also affected by extragranular diffusion processes, such as surface adsorption and liquid film diffusion [41,42,59].

Thermodynamic parameters of adsorption
Thermodynamic parameters help to describe the energy changes involved in the adsorption process [41]. To clarify the adsorption mechanism of atrazine in soil, the adsorption status of atrazine at different temperatures (288, 298 and 308 K) was analyzed and the influence of temperature on the equilibrium adsorption coefficient was assessed based on the standard Gibbs free energy ( rG 0 m ) formula (Eq. (8)). Furthermore, the relevant thermodynamic parameters of standard enthalpy ( rH 0 m ), and standard entropy ( rS 0 m ) were calculated using Eq. (9) as follows [59]: where K 0 is the adsorption equilibrium constant; R is the gas molar constant (8.314 J (mol K) −1 ); and T is absolute temperature. According to the slope and intercept of the plot ln K 0 versus 1/T, the values of rS 0 m and rH 0 m were obtained.
In general, the rG 0 m represents the spontaneity of the chemical reaction and is used to evaluate whether the adsorption reaction occurs spontaneously [27,65]. The thermodynamic parameters of atrazine adsorption by the two soils are shown in Table 5, with values of rG 0 m being < 0, indicating the sorption of atrazine by soils was spontaneous. The ranked order of the absolute values for atrazine absorption by soil at different temperatures was as follows: rG 0 m 298 K > rG 0 m 308 K > rG 0 m 318 K . The adsorption of atrazine in the soil increased with the increase in temperature, indicating adsorption at higher temperatures was easier than at low temperatures. It is generally believed that when the values of rG 0 m are between 0 and −20 kJ mol −1 , physical adsorption is occurring with van der Waals forces playing a dominant role, resulting in the adsorption effect being small and desorption occurring easily. When the values of rG 0 m are between − 80 and − 400 kJ mol −1 , chemical adsorption  is occurring with chemical bonding playing a dominant role, resulting in a larger adsorption energy and a higher tendency for irreversible adsorption [41,66]. In this study, the values of rG 0 m were between − 0.418 and − 0.849 kJ mol −1 for Wanzhou soil and between − 0.222 and − 0.794 kJ mol −1 for Wushan soil indicating that the adsorption of atrazine by both soil samples occurred by physical adsorption and confirming that the adsorption of the atrazine by soil is feasible and spontaneous.
The standard enthalpy changes ( rH 0 m ) for the adsorption of atrazine in WanZhou and WuShan soils were − 5.38 and − 9.27 kJ mol −1 , respectively. The value of rH 0 m was less than 0, indicating that the solid phase of atrazine adsorbed from water to the soils was via an exothermic process, releasing large amounts of heat [67]. This result showed that adsorption was mainly based on physical properties, which is consistent with the results of rG 0 m . Studies have previously reported that when the standard enthalpy ranges from 40 to 120 kJ mol −1 , adsorption occurs via a chemical process; otherwise, adsorption is via physical processes [68]. Therefore, the adsorption of the atrazine in this study appeared seems to occur via physical sorption. The standard entropy changes ( rS 0 m ) of the adsorption of atrazine in Wanzhou and Wushan were 16.9 and 28.7 J (mol K) −1 , respectively. The positive values of rS 0 m indicate that soils had an affinity for atrazine and also reflect changes in the structures of atrazine and soils. In addition, the positive rS 0 m value indicated that randomness of the solid-liquid interface increased during the adsorption process.

Conclusions
The effects of temperature and contact time on the adsorption of atrazine in the two soils (Wanzhou and Wushan) of the Three-Gorges Reservoir were investigated, showing that adsorption was favorable. Within the range of initial concentration variation, the adsorption of atrazine by soil decreased with increasing temperature (298-318 K). The atrazine adsorption data for both soil samples could be fitted by the Freundlich and Langmuir models, indicating that single-layer adsorption may occur, with heterogeneous active site energy distribution on the soil surface. The adsorption of atrazine can be described by the pseudo-second-order model and internal particle diffusion was found not to be the only ratecontrolling step. As the adsorption time progressed, the amount of adsorbed atrazine in the soil increased gradually, indicating that the high initial concentration solution can provide sufficient power for mass transfer resistance between the atrazine molecular reduction solution and the adsorbent. During the adsorption of atrazine in soils, rG 0 m < 0, rH 0 m < 0 and rS 0 m > 0, indicating that the adsorption process is a spontaneous exothermic process and that the degree of chaos in the whole adsorption process increased. Generally, the adsorption capacity was low, indicating that atrazine can easily penetrate the soil vadose zone to pollute groundwater. Therefore, much attention should be paid to the environmental behavior of pesticide exposed soil moisture fluctuations, leading to the potential broad transport of pollution.