Adsorption of Heavy Metals on Soil Collected from Lixisol of Typical Karst Areas in the Presence of CaCO3 and Soil Clay and Their Competition Behavior

The content of heavy metals in the soil in Guizhou Province, which is a high-risk area for heavy metal exposure, is significantly higher than that in other areas in China. Therefore, the objective of this study was to evaluate the ability of CaCO3 and clay to accumulate heavy metals in topsoil sample collected from Lixisol using the method of indoor simulation. The results showed that the contents of Cu, Zn, Cd, Cr, Pb, Hg and As in the soil sample were 10.8 mg/kg, 125 mg/kg, 0.489 mg/kg, 23.5 mg/kg, 22.7 mg/kg, 58.3 mg/kg and 45.4 mg/kg, respectively. The soil pH values increased with the CaCO3 concentration in the soil, and the fluctuation of the soil pH values was weak after the CaCO3 concentrations reached 100 g/kg. The adsorption capacity of lime soil increased by approximately 10 mg/kg on average, and the desorption capacity decreased by approximately 300 mg/kg on average. The desorption of all heavy metals in this study did not change with increasing clay content. Pseudo-second-order kinetics were more suitable for describing the adsorption kinetics of heavy metals on the soil material, as evidenced by the higher R2 value. The Freundlich model can better describe the adsorption process of As on lime soil. The process of As, Cr, Cd and Hg adsorption on the soil sample was spontaneous and entropy-driven. Additionally, the process of Cu and Pb adsorption on the soil materials was spontaneous and enthalpy-driven. Generally, the adsorption and desorption of heavy metals in polluted soil increased and decreased, respectively, with increasing CaCO3 content. The effect of calcium carbonate on the accumulation of heavy metals in soil was greater than that of clay. In summary, CaCO3 and pH values in soil can be appropriately added in several areas polluted by heavy metals to enhance the crop yield and reduce the adsorption of heavy metals in soils.


Introduction
There are two main soil types, i.e., Lixisols and Acrisols, in Guizhou Province, of which the former are characterized by high content of calcium and magnesium carbonate (CaCO 3 , MgCO 3 ), and carbonate sample collected from Lixisol using kinetics, adsorption isotherms and thermodynamic study; and (3) investigate the effect of pH and electrolytes on adsorption capacity.

Materials
Lixisols, a typical soil in Guizhou Province, Southwest China, was selected in this study. All soil samples used in the experiment were collected from the topsoil with a stainless-steel shovel and packed back with cloth bags. After natural air drying at room temperature and removal of plants, gravel and other impurities in the soil, the soil was ground with a porcelain mortar, a container for crushing experimental materials in an experiment, and sifted to 100 mesh. After sifting, all the soil was collected in plastic bags for future use.

Analysis Method
The collected soil was divided into many parts, which was used to analyze the several indexes. For one of the soil samples, the pH of the soil was determined by potentiometry (soil-to-water ratio = 1:2.5) [22,23]. Soil organic matter was measured by the high-temperature external heat potassium dichromate oxidation volumetric method [24,25]. The calcium carbonate content in the soil was calculated using the volumetric titration method [26]. The soil particle size distribution was gauged by the hydrometer method [27]. The soil samples were digested by a microwave digester using HNO 3 -H 2 O 2 -HF (5:2:5), and acid was removed by a graphite furnace [28]. Hg and As were determined by an atomic fluorescence photometer (AFS-2100, Beijing Haiguang Instrument Co., Ltd., Beijing, China). Pb, Cu, Zn, Cd and Cr of soils were determined by inductively coupled plasma-optical emission spectroscopy (ICP-OES) (Prodigy XP, Leeman Labs Inc., Hudson, NY, USA).

Effect of Different Calcium Carbonate Contents on Adsorption and Desorption
Eleven soil samples of 0.5 g (accurate to 0.0001 g) were accurately weighed and added to 100 mL centrifuge tubes. Super pure CaCO 3 was added to each sample at 0, 0.25, 0.5, 1, 3, 5, 7, 10, 15, 20 and 25 g/kg. The standard heavy metal (Cd, Pb, Cr, Zn, As and Hg) mixed solutions (nitric acid medium), which were purchased from the National Center of Analysis and Testing for Non-ferrous Metals and Electronic Materials, of 50 mL with a concentration of 20 mg/L, were added to the above mixture. These compounds were oscillated on a constant-temperature oscillator (25 ± 1 • C). After 2 h, the cells were removed and cultured in a constant-temperature incubator (25 ± 1 • C) for 22 h and then centrifuged for 30 min (4000 r/min). The supernatant was taken and measured with the aid of a computer, and the formula used to calculate the adsorption capacity is as follows: where q e (mg/kg) is the capacity of soil adsorption/desorption in the solution at equilibrium; C i (mg/L) is the initial concentration of heavy metals in the supernatant; C e (mg/L) is the equilibrium concentration of heavy metals in the supernatant; v (mL) is the volume of the total mixed solution used; and m (g) is the weight of the added soil material. The supernatants of the above soil samples were removed, and 20 mL of 0.1 mol/L NaNO 3 was added. The mixtures were oscillated on a constant temperature oscillator (25 ± 1 • C) for 2 h. The solutions were cultured in a constant-temperature incubator (25 ± 1 • C) for 22 h after transfer to a volumetric flask and then centrifuged at 4000 r/min for 30 min. Finally, the concentration of heavy metals in the supernatant was determined using ICP-OES. The desorption capacity was the concentration difference of heavy metals in the supernatant before and after adding NaNO 3 . The formula for calculation is as follows: where S (mg/kg) is the desorption capacity of soil in the solution at equilibrium; C 1 and C 2 are the concentrations of heavy metals in the supernatant before and after adding NaNO 3 (mg/L), respectively; v (mL) is the volume of the total mixed solution used; and m (g) is the capacity of the added soils. Figure 1 shows the adsorption/desorption flow chart.

Particle Size Distribution of Soil
The particle size distribution of soil was measured by the hydrometer method. A soil sample of 50 g passing a 2 mm sieve was weighed into a triangular flask. The appropriate dispersant, 0.5 mol/L NaOH solution to acidic soil of 40 mL or sodium hexametaphosphate solution to alkaline soil of 60 mL, 250 mL pure water were successively added to the above triangular flask. After the mixture soaked overnight, the mixtures in flasks were boiled on the heating plate for 1 h (using the defoamer to prevent overflow). After cooling, the sieve was washed through 0.2 mm, and the suspension was washed into a 1 L measuring cylinder. After all the suspension had passed, the mixtures were sieved through 0.2 mm and the suspension was washed into a measuring cylinder of 1 L. After all the suspension had passed through, the sieve was washed with pure water. The particles left on the screen were washed into the aluminum box, dried and weighed. The part washed into the measuring cylinder was added with water to a constant volume of 1 L.
The solution in the measuring cylinder was stirred up and down repeatedly for one minute with an agitator. The time was counted at the end of mixing, which were respectively 30 s, 1 min, 2 min, 4 min, 8 min, 15 min, 30 min, 1 h, 2 h, 4 h, 8 h, 24 h, 48 h, and put it into the hydrometer to read, and measure and record the water temperature before reading. The particle size content was calculated according to the following equation, where ρ 1 is the hydrometer reading, g/L; ρ 2 is the correction value of hydrometer calibration curve, g/L; ρ 3 is the calibration temperature value of hydrometer, g/L; ρ 0 is the corrected value of dispersant for hydrometer reading, g/L; v is the volume of suspension, mL; m is the weight of soil sample, g. The diameter (d) of soil is calculated by the following Equation (4): where d is the diameter of soil sample, mm; η is the viscosity coefficient of water; L is the settlement depth of soil particles, g is the acceleration of gravity, cm/s; ρs is the soil particle density, g/cm 3 ; ρf is the density of water, g/cm 3 ; t is settlement time, s.

Adsorption Kinetic
Prepared solutions of 1 L contained 20 mg/L heavy metals and 0.01 mol/L NaNO 3 (as supporting electrolyte) in a beaker. A 25.00 g aliquot of soil material was added to each prepared solution of 1 L. At room temperature, the mixtures were stirred at a constant speed for 1 min, 2 min, 4 min, 6 (20 mL) was selected and centrifuged at 4000 r/min for 3 min. Pseudo-first-order kinetics, pseudo-second-order kinetics and particle internal diffusion were used to fit the data of the adsorption process of heavy metals on the soil sample [29,30], where q e (mg/g) and q t (mg/g) are the adsorption capacity at equilibrium and after time t (min), respectively; k is the kinetic adsorption rate constant. k int (mg/g·min −0.5 ) is the diffusion rate constant in particles.

Adsorption Isotherm
Six soil samples (0.5 g, accurate to 0.0001 g) were weighed into a 100 mL centrifuge tube, and then 20 mL of heavy metals (Cd, Pb, Cr, Zn, As, Hg) with concentrations of 0.00, 10.00, 20.00, 30.00, 40.00, and 50.00 mg/L were added to the centrifuge tube. NaNO 3 (0.01 mol/L) as the supporting electrolyte was dropped into these mixtures, which were oscillated at 25 • C for 24 h to reach adsorption equilibrium. After centrifugation and filtration, the supernatant was used to determine the concentration of heavy metals. The Langmuir and Freundlich models were frequently utilized to describe the isotherm adsorption effect at the soil/water interface [31,32]. Their equations are as follows: where q e (mg/g) is the adsorption capacity of adsorbate on the adsorbent surface at equilibrium; k (L/mg) is the Langmuir constant; Q m (mg/g) is the maximum adsorption capacity of adsorbate on the adsorbent surface; C e (mg/L) is the concentration of adsorbate in solution at equilibrium; k F (mg/g) is the Freundlich constant; and n is the index of anisotropy, which is usually greater than 1. A larger value indicates the significance of the nonlinear characteristic of the adsorption isotherm.

Estimation of Thermodynamic Parameters
The change in Gibbs free energy associated with a chemical reaction can be used to evaluate the spontaneity of the adsorption reaction [33]. Thermodynamic parameters can be obtained as follows [34]: where C s is the initial concentration of heavy metals on the soil material; C e is the concentration of heavy metals in solution at equilibrium; and R is the gas constant (8.31 J·mol −1 ·K −1 ).

Effect of Initial pH on Adsorption
Eight soil samples (0.5 g, accurate to 0.0001 g) were weighed into a centrifuge tube of 100 mL, A 20 mL aliquot of 20 mg/L heavy metal (0.01 mol/L NaNO 3 as supporting electrolyte) was added into the above centrifuge tube, and the pH values of the mixture were adjusted between 3.00 and 10.00 using 0.1 mol/L NaOH and HNO 3 . The pH values in this experiment were set as 8 gradients (3, 4, 5, 6, 7, 8, 9 and 10). After shaking, centrifuging and filtering, the supernatant was used to determine the concentration of heavy metals. Their adsorption capacity was calculated as in Section 2.3.

Effect of Temperature on Adsorption
The processing steps were the same as those described in Section 2.3, and the adsorption isothermal experiments were carried out at 50 • C, 45 • C, 40 • C, 35 • C and 30 • C. Their adsorption capacity was calculated as in Section 2.3.

Influence of Electrolyte on Adsorption
Five soil samples (0.5 g, accurate to 0.0001 g) were weighed into 100 mL centrifuge tubes. Different concentrations of NaNO 3 (10 −1 , 10 −2 , 10 −3 , 10 −4 and 10 −5 mol/L) were added to the above centrifuge tubes. The treatment process and adsorption capacity were calculated as in Section 2.3.

Physical and Chemical Indexes of Topsoil Sample Collected from Lixisol
The pH of topsoil in Guizhou Province was 7.64 in this study, which corresponds to weak alkaline soil ( Figure 2a). The contents of organic matter and calcium carbonate were 70.3 g/kg and 50.9 g/kg, respectively, which exhibited that the organic matter in soil of study area was rich according to Nutrient grading standard of the second national soil survey (China). The soil organic matter can inhibit or promote the adsorption of heavy metals in soil through a series of reactions, such as ion exchange adsorption, complexation and chelation between heavy metals and other functional groups, which affect the precipitation, migration, transformation and bioavailability of heavy metals [35]. The particle size distribution of the soil sample in our study was 37.78% sand, 26.79% powder and 35.44% clay, which pertained to loamy clay. The moisture content of calcareous soil was 3.58%.
pH values was weak after the CaCO3 concentrations reached 100 g/kg. In soil, the effect of CaCO3 on pH was due to the hydrolysis of carbonate (MCO3 + H2O ⇋ M 2+ + HCO -3 + OH -) [36]. Liu et al. [37] reported that CaCO3 concentrations slightly affected the soil pH when the content of calcium carbonate was 50-100 g/kg, which is consistent with our study. The different clay contents caused a difference in soil pH values, and the maximum and minimum pH values of the soil were 6.44 and 6.39, respectively.

Effect of Calcium Carbonate and Soil Clay on Adsorption and Desorption of Heavy Metals in Topsoil Sample Collected from Lixisol
Generally, the adsorption capacity of heavy metals in soil was higher than the desorption capacity, which implied that these heavy metals could be fixed by the soil sample, and few heavy metals were desorbed; thus, the migration capacity of fixed heavy metals will also be reduced. The soil material is extensively distributed in Guizhou Province, which is a reason for the higher heavy metal concentrations in the soil in Guizhou Province. As shown in Figure 3 and Table 1, the order of the maximum adsorption capacities for heavy metals is Cu > Pb > Zn > Cd > Hg > Cr > As. For Cu, the maximum adsorption capacity was 807 mg/kg, and the minimum adsorption capacity was 806 mg/kg. The maximum desorption capacity was 172 mg/kg, and the minimum desorption capacity was 18.1 mg/kg. The maximum and minimum adsorption capacities of Pb were 802 mg/kg and 799 mg/kg, respectively. The maximum and minimum desorption capacities of Pb were 3.45 mg/kg and −2.60 mg/kg, respectively. The maximum and minimum adsorption capacities of Zn were 802 mg/kg and 800 mg/kg, respectively. The maximum desorption capacity of Zn is 674 mg/kg and the minimum desorption capacity is 282 mg/kg. The maximum adsorption and desorption capacities of Cd were 795 mg/kg and 661 mg/kg, respectively. The minimum adsorption and desorption capacities of Cd were 790 mg/kg and 409 mg/kg, respectively. For Hg, the maximum and minimum adsorption capacities were 771 mg/kg and 758 mg/kg, respectively. The maximum desorption capacity was −0.436 mg/kg, and the minimum desorption capacity was −7.70 mg/kg. The maximum adsorption capacity of Cr was 767 mg/kg, and the minimum adsorption capacity was 760 mg/kg. The maximum desorption capacity of Cr was 510 mg/kg, and the minimum desorption capacity was 268 mg/kg. The maximum adsorption capacity of As was 763 mg/kg, and its minimum value was 666 mg/kg. The maximum and minimum desorption capacities of As were 8.90 mg/kg and −2.36 mg/kg. The desorption efficiency for Pb, Cd, Zn, Cu, Cr, Hg and As is 0.43%, 83.09%, 83.98%, 20.83%, 35.14%, not detected and 0.13%, respectively.   [36]. Liu et al. [37] reported that CaCO 3 concentrations slightly affected the soil pH when the content of calcium carbonate was 50-100 g/kg, which is consistent with our study. The different clay contents caused a difference in soil pH values, and the maximum and minimum pH values of the soil were 6.44 and 6.39, respectively.

Effect of Calcium Carbonate and Soil Clay on Adsorption and Desorption of Heavy Metals in Topsoil Sample Collected from Lixisol
Generally, the adsorption capacity of heavy metals in soil was higher than the desorption capacity, which implied that these heavy metals could be fixed by the soil sample, and few heavy metals were desorbed; thus, the migration capacity of fixed heavy metals will also be reduced. The soil material is extensively distributed in Guizhou Province, which is a reason for the higher heavy metal concentrations in the soil in Guizhou Province. As shown in Figure 3 and Table 1, the order of the maximum adsorption capacities for heavy metals is Cu > Pb > Zn > Cd > Hg > Cr > As. For Cu, the maximum adsorption capacity was 807 mg/kg, and the minimum adsorption capacity was 806 mg/kg. The maximum desorption capacity was 172 mg/kg, and the minimum desorption capacity was 18.1 mg/kg. The maximum and minimum adsorption capacities of Pb were 802 mg/kg and 799 mg/kg, respectively. The maximum and minimum desorption capacities of Pb were 3.45 mg/kg and −2.60 mg/kg, respectively. The maximum and minimum adsorption capacities of Zn were 802 mg/kg and 800 mg/kg, respectively. The maximum desorption capacity of Zn is 674 mg/kg and the minimum desorption capacity is 282 mg/kg. The maximum adsorption and desorption capacities of Cd were 795 mg/kg and 661 mg/kg, respectively. The minimum adsorption and desorption capacities of Cd were 790 mg/kg and 409 mg/kg, respectively. For Hg, the maximum and minimum adsorption capacities were 771 mg/kg and 758 mg/kg, respectively. The maximum desorption capacity was −0.436 mg/kg, and the minimum desorption capacity was −7.70 mg/kg. The maximum adsorption capacity of Cr was 767 mg/kg, and the minimum adsorption capacity was 760 mg/kg. The maximum desorption capacity of Cr was 510 mg/kg, and the minimum desorption capacity was 268 mg/kg. The maximum adsorption capacity of As was 763 mg/kg, and its minimum value was 666 mg/kg. The maximum and minimum desorption capacities of As were 8.90 mg/kg and −2.36 mg/kg. The desorption efficiency for Pb, Cd, Zn, Cu, Cr, Hg and As is 0.43%, 83.09%, 83.98%, 20.83%, 35.14%, not detected and 0.13%, respectively.   As shown in Table 2, the order of the maximum adsorption capacities for heavy metals is Pb > Cu > Hg > Zn > As > Cd > Cr. Firstly, the maximum and minimum adsorption capacities of Pb were 434 mg/kg and 399 mg/kg, respectively. The maximum desorption capacity of Pb was 328 mg/kg, and the minimum desorption capacity was 293 mg/kg. The maximum adsorption capacity of Cu was 431 mg/kg, and the minimum adsorption capacity was 378 mg/kg. The maximum desorption capacity was 412 mg/kg, and the minimum desorption capacity was 375 mg/kg. The maximum adsorption capacity of Hg was 415 mg/kg, and its minimum value was 389 mg/kg. The maximum and minimum desorption capacities were 98.7 mg/kg and 62.9 mg/kg, respectively. For the adsorption capacity of Zn, the maximum value was 372 mg/kg, and the minimum value was 288 mg/kg. The maximum desorption capacity of Zn was 391 mg/kg, and the minimum desorption capacity was 344 mg/kg. The maximum and minimum adsorption capacities of As were 368 mg/kg and 353 mg/kg, respectively. The maximum desorption capacity was 288 mg/kg, and the minimum desorption capacity was 272 mg/kg. The maximum adsorption capacity of Cd was 292 mg/kg, and the minimum adsorption capacity was 236 mg/kg. The maximum desorption capacity of Cd was 413 mg/kg, and the minimum desorption capacity was 378.1 mg/kg. For the adsorption capacity of Cr, its maximum and minimum values were 231 mg/kg and 145 mg/kg, respectively. The maximum desorption capacity was 448 mg/kg, and the minimum desorption capacity was 434 mg/kg. However, the adsorption concentration of heavy metals in soil increased with the contain of clay (from 382 g/kg to 454 g/kg). The phenomenon may be that the clay in the soil had a large surface energy, which made its adsorption capacity relatively strong. Therefore, the migration of pollutants was limited. Figure 4 shows that the effect of clay on heavy metals in lime soil was not as significant as that of CaCO3. Their adsorption capacities gradually increased with increasing clay content. However, the adsorption capacity of As was slightly decreased with increasing clay content. Generally, the desorption of all heavy metals in this study did not change with increasing clay content. The desorption efficiencies for Pb, Cd, Zn, Cu, Cr, Hg and As are 93.79%, 79.63%, 20.12%, 99.65%, 94.88%, 54.55% and 69.44%, respectively.  The adsorption capacity of Cu, Zn, Cd and Pb did not obviously change with the increase in CaCO 3 content in soil. However, the adsorption capacities of Hg and Cr slightly increased with increasing CaCO 3 content in the soil sample. A negative correlation between the adsorption capacities of As and the CaCO 3 content was observed because the adsorption capacity of As decreased significantly as the CaCO 3 content of the soil material increased. Zhong et al. [38] reported effects of CaCO 3 addition on uptake of heavy metals and arsenic in paddy fields, and the exchangeable arsenic content in the soils decreased significantly by adding CaCO 3 . The reason may be that Ca in soil can form insoluble precipitation with As, which was conducive to the fixation of As. The desorption capacities of Cu, Cd and Zn were decreased with increasing CaCO 3 content in the soil sample. However, with the increase in CaCO 3 content in the soil material, the desorption capacities of As and Pb did not significantly change. The desorption capacities of Cr were increased with increasing CaCO 3 content in soil sample. Therefore, with the increase in the CaCO 3 content in the soil sample, the adsorption and fixation of Cu, Cd and Zn were enhanced, and the adsorption and fixation of Cr and As were weakened, but Pb had no obvious effect. With the increase in CaCO 3 content, the adsorption capacity of soil material was increased by approximately 10 mg/kg on average, and the desorption capacity decreased by approximately 300 mg/kg on average. The pH value of soil increased after adding calcium carbonate, which promoted the precipitation reaction of exchangeable Cd 2+ , Cu 2+ , Pb 2+ and Zn 2+ with carbonate and hydroxide in soil. Meanwhile, the coprecipitation of Ca 2+ and metal ions is also conducive to the transformation of the exchange state to other forms. Additionally, the increase of soil pH can enhance the chelating capacity of clay and organic matter, promote the adsorption capacity of soil, and reduce the desorption of heavy metals, thus decreasing the solubility of heavy metals in soil [39]. In contrast, the adsorption capacity of the soil materials was continuously increased with the decrease in CaCO 3 content. This may be because with the decrease of pH value, the negative charges on the surface of clay minerals, hydrated oxides and organic matter in the soil decreased, and the adsorption capacity of heavy metals weakened. Meanwhile, the specific adsorption of heavy metals on the surface of oxides and the stability of soil organic matter metal complexes decreased with the decrease of pH value. With the decrease of pH, the concentrations of H + and Fe 2+ in soil increase, and the competitive adsorption with heavy metals increases, which is not conducive to the adsorption of heavy metals. The decrease of pH is not conducive to the formation of metal hydrate ions, and the affinity of hydrated ions on soil adsorption sites is significantly lower than that of heavy metal ions, which is not conducive to the adsorption of heavy metals on soil [40]. However, the effects on different elements are different, such as inhibiting the transport of nitrate nitrogen from the underground of wheat to the aboveground, promoting the uptake of chromium in the aerial parts of Chinese cabbage, and enhancing the enrichment of lead in lettuce stems and leaves.
As shown in Table 2, the order of the maximum adsorption capacities for heavy metals is Pb > Cu > Hg > Zn > As > Cd > Cr. Firstly, the maximum and minimum adsorption capacities of Pb were 434 mg/kg and 399 mg/kg, respectively. The maximum desorption capacity of Pb was 328 mg/kg, and the minimum desorption capacity was 293 mg/kg. The maximum adsorption capacity of Cu was 431 mg/kg, and the minimum adsorption capacity was 378 mg/kg. The maximum desorption capacity was 412 mg/kg, and the minimum desorption capacity was 375 mg/kg. The maximum adsorption capacity of Hg was 415 mg/kg, and its minimum value was 389 mg/kg. The maximum and minimum desorption capacities were 98.7 mg/kg and 62.9 mg/kg, respectively. For the adsorption capacity of Zn, the maximum value was 372 mg/kg, and the minimum value was 288 mg/kg. The maximum desorption capacity of Zn was 391 mg/kg, and the minimum desorption capacity was 344 mg/kg. The maximum and minimum adsorption capacities of As were 368 mg/kg and 353 mg/kg, respectively. The maximum desorption capacity was 288 mg/kg, and the minimum desorption capacity was 272 mg/kg. The maximum adsorption capacity of Cd was 292 mg/kg, and the minimum adsorption capacity was 236 mg/kg. The maximum desorption capacity of Cd was 413 mg/kg, and the minimum desorption capacity was 378.1 mg/kg. For the adsorption capacity of Cr, its maximum and minimum values were 231 mg/kg and 145 mg/kg, respectively. The maximum desorption capacity was 448 mg/kg, and the minimum desorption capacity was 434 mg/kg. However, the adsorption concentration of heavy metals in soil increased with the contain of clay (from 382 g/kg to 454 g/kg). The phenomenon may be that the clay in the soil had a large surface energy, which made its adsorption capacity relatively strong. Therefore, the migration of pollutants was limited.  Figure 4 shows that the effect of clay on heavy metals in lime soil was not as significant as that of CaCO 3 . Their adsorption capacities gradually increased with increasing clay content. However, the adsorption capacity of As was slightly decreased with increasing clay content. Generally, the desorption of all heavy metals in this study did not change with increasing clay content. The desorption efficiencies for Pb, Cd, Zn, Cu, Cr, Hg and As are 93.79%, 79.63%, 20.12%, 99.65%, 94.88%, 54.55% and 69.44%, respectively.  Note: AQ and DQ represent adsorption and desorption capacity, respectively.

Adsorption Kinetics of Heavy Metals Adsorbed on Topsoil Sample Collected from Lixisol
The rate of metal ions moving in solution and fixing on the surface of the soil determined the efficiency of the adsorption process. The residence time of a solute on the surface of the soil adsorbent was an important factor to judge whether the adsorption process was completed and to calculate the total amount of solute adsorption. The study of adsorption dynamic characteristics can be used to explore the possible adsorption mechanism from the perspective of the adsorption path.
The adsorption capacities of heavy metals on the soil sample dramatically increased with increasing contact time ( Figure 5). When the adsorption capacity reached a certain value, it weakly adsorption quality desorption quality a Adsorption capacity (mg/kg) Clay particle (g/kg) Clay particle (g/kg) adsorption quality desorption quality Clay particle (g/kg)

Adsorption Kinetics of Heavy Metals Adsorbed on Topsoil Sample Collected from Lixisol
The rate of metal ions moving in solution and fixing on the surface of the soil determined the efficiency of the adsorption process. The residence time of a solute on the surface of the soil adsorbent was an important factor to judge whether the adsorption process was completed and to calculate the total amount of solute adsorption. The study of adsorption dynamic characteristics can be used to explore the possible adsorption mechanism from the perspective of the adsorption path.
The adsorption capacities of heavy metals on the soil sample dramatically increased with increasing contact time ( Figure 5). When the adsorption capacity reached a certain value, it weakly increased or decreased with increasing contact time. The adsorption rate of the soil material can reach 50% in 5 min, 80% in 15 min, 95% in 90 min and 99% after 720 min. The adsorption process of heavy metals by the soil sample was mainly concentrated in the first 90 min, and the subsequent adsorption was relatively slow. The reason for this difference may be physical adsorption (fast adsorption stage) and chemical adsorption (slow adsorption stage). Generally, physical adsorption can achieve local equilibrium in a few milliseconds to a few seconds; however, chemical adsorption requires a certain activation energy to adsorb, resulting in a slow adsorption process.

Adsorption Isotherms of Heavy Metals Adsorbed on Topsoil Sample Collected from Lixisol
The shape and change rule of the adsorption isotherm can be used to understand the action intensity of the adsorbate and adsorbent, the state of the adsorbate on the interface and the structure of the adsorption layer. All heavy metals exhibit a larger slope at a lower equilibrium concentration (the heavy metal concentration which was no longer improved after the different initial concentrations of heavy metals adsorbed by the soil sample) ( Figure 6). The adsorption capacities of heavy metals were added with increasing initial heavy metal concentrations; however, the slope of the isotherm adsorption curve was gentle at higher equilibrium concentrations. The slope of As was the largest, followed by those of Hg, Cr and Pb, which were dramatically adsorbed and then slowly adsorbed. The slopes of Cu, Zn and Cd slowly increased with increasing initial concentration. The reason may be that the concentration of adsorbate in the solution was enhanced with increasing initial concentration, which was more favorable for the combination with adsorbate and adsorbent. However, the capacity of each heavy metal combined with the soil material was different since these heavy metals had different properties. Hence, the maximum adsorption capacity of each heavy metal has a higher difference, of which As was the most obvious with the increase in initial As concentration. In this study, As was selected for specific discussion due to its sensitivity to the initial concentration of heavy metals in soil sample (Figure 7). Sustainability 2020, 12, x FOR PEER REVIEW 12 of 18

Adsorption Isotherms of Heavy Metals Adsorbed on Topsoil Sample Collected from Lixisol
The shape and change rule of the adsorption isotherm can be used to understand the action intensity of the adsorbate and adsorbent, the state of the adsorbate on the interface and the structure of the adsorption layer. All heavy metals exhibit a larger slope at a lower equilibrium concentration (the heavy metal concentration which was no longer improved after the different initial concentrations of heavy metals adsorbed by the soil sample) ( Figure 6). The adsorption capacities of heavy metals were added with increasing initial heavy metal concentrations; however, the slope of the isotherm adsorption curve was gentle at higher equilibrium concentrations. The slope of As was the largest, followed by those of Hg, Cr and Pb, which were dramatically adsorbed and then slowly adsorbed. The slopes of Cu, Zn and Cd slowly increased with increasing initial concentration. The reason may be that the concentration of adsorbate in the solution was enhanced with increasing initial concentration, which was more favorable for the combination with adsorbate and adsorbent. However, the capacity of each heavy metal combined with the soil material was different since these heavy metals had different properties. Hence, the maximum adsorption capacity of each heavy metal has a higher difference, of which As was the most obvious with the increase in initial As concentration. In this study, As was selected for specific discussion due to its sensitivity to the initial concentration of heavy metals in soil sample (Figure 7). Four types of adsorption isotherm lines, i.e., "S", "L", "H" and "C", were used to discuss the adsorption isotherm in a study by Giles et al. [41]. The isotherm adsorption curve of As on lime soil had both "S" and "L" types. S-type adsorption was equivalent to a single multilayer reversible adsorption process occurring on the surface of a nonporous solid or macroporous solid, which was also consistent with the physical structure of soil [42]. The curve was steep first and then gentle with caused by monolayer adsorption or micropore filling of a small amount of silica gel in the clay. Figure 8 shows the Langmuir and Freundlich isotherm model fitting of As adsorption on lime soil. The isotherm adsorption parameters obtained by the fitting are shown in Table 4. The correlation coefficients (R 2 ) of the Langmuir and Freundlich models were 0.817 and 0.920, respectively, which implied that the Freundlich model can better describe the adsorption process of As on lime soil. In addition, the k and n values of the Freundlich model were 166 L/mg and 2.37.

Thermodynamic Study
The thermodynamics of the adsorption of heavy metals on lime soil can be used to study the exothermic, spontaneous or endothermic nature of heavy metal adsorption on the interface between the adsorbate and adsorbent [44]. The positive values of ΔS° suggested an increase in randomness at Four types of adsorption isotherm lines, i.e., "S", "L", "H" and "C", were used to discuss the adsorption isotherm in a study by Giles et al. [41]. The isotherm adsorption curve of As on lime soil had both "S" and "L" types. S-type adsorption was equivalent to a single multilayer reversible adsorption process occurring on the surface of a nonporous solid or macroporous solid, which was also consistent with the physical structure of soil [42]. The curve was steep first and then gentle with an increase in concentration at low concentrations, reaching the first inflection point to realize the single-layer adsorption saturation on the surface, which was the single-layer saturated adsorption capacity. The adsorbate diffused into the micropores of the adsorbent with increasing concentration. Because the surface area of the adsorbent was smaller than the internal surface area, the adsorption curve was concave upward. The adsorption capacity will not increase when the pores are completely filled and the adsorption sites are completely occupied, and the adsorption capacity will be controlled by the micropore volume, whose capillary diameter is slightly larger than the size of the adsorbate [43]. The adsorption process was mainly the result of volume filling (the adsorption capacity increased with increasing adsorbate concentration). When the binding sites or micropores on the adsorbent surface were completely filled, the adsorption reached saturation, and the adsorption capacity no longer increased with increasing adsorbate concentration. The L-type curve may be caused by monolayer adsorption or micropore filling of a small amount of silica gel in the clay. Figure 8 shows the Langmuir and Freundlich isotherm model fitting of As adsorption on lime soil. The isotherm adsorption parameters obtained by the fitting are shown in Table 4. The correlation coefficients (R 2 ) of the Langmuir and Freundlich models were 0.817 and 0.920, respectively, which implied that the Freundlich model can better describe the adsorption process of As on lime soil. In addition, the k and n values of the Freundlich model were 166 L/mg and 2.37.

Thermodynamic Study
The thermodynamics of the adsorption of heavy metals on lime soil can be used to study the exothermic, spontaneous or endothermic nature of heavy metal adsorption on the interface between the adsorbate and adsorbent [44]. The positive values of ΔS° suggested an increase in randomness at 10 20 30 40 50 200 Initial As concentration (mg/L)

Thermodynamic Study
The thermodynamics of the adsorption of heavy metals on lime soil can be used to study the exothermic, spontaneous or endothermic nature of heavy metal adsorption on the interface between the adsorbate and adsorbent [44]. The positive values of ∆S 0 suggested an increase in randomness at the solid/solution interface [45]. The positive value obtained for ∆H 0 indicates the endothermic nature of heavy metal adsorption, which demonstrates that this process consumes energy [46]. Therefore, the process was a spontaneous and entropy-driven process. The negative value of ∆G 0 showed the feasibility and spontaneity of heavy metal adsorption onto the soil sample [47]. Therefore, the process of As, Cr, Cd and Hg adsorption onto the soil material was spontaneous and entropy-driven (Table 5). Additionally, the process of Cu and Pb adsorbed by the soil sample was spontaneous and enthalpy-driven. The process of Zn adsorption onto the soil material may be doubly enthalpyand entropy-driven. The adsorption capacities were enhanced with increasing temperature, which may be due to the endothermic reaction. The adsorption capacity of heavy metals in soil was more easily adsorbed on the surface of soil particles with an increase in temperature, and there was no obvious desorption phenomenon due to the increase in temperature, which may be caused by stable chemical adsorption during the whole process. The adsorption capacities of heavy metals in the soil material fluctuated with increasing temperature; however, there was no significant increase or decrease. This may be because the adsorption has reached a balance (the binding sites of the surface and inner pores of the soil particles have been filled). Both adsorption and desorption were promoted with increasing temperature, which showed that the adsorption amount of heavy metals in lime soil did not obvious change with increasing temperature. The reason why the adsorption capacity of heavy metals in the soil material decreased with increasing temperature may be that adsorption was an exothermic reaction. The increase in temperature was an inhibiting factor for the total amount of heavy metals adsorbed in the soil sample, or it may be that after the adsorption reached a certain temperature, the adsorption process was changed from chemical adsorption to physical adsorption (the desorption rate of heavy metals from soil was greater than the adsorption rate).

Effect of pH and Electrolyte on Adsorption of Heavy Metals by Lime Soil
The adsorption capacities of Cr decreased with increasing pH, which may be caused by the unique physical properties of Cr. The adsorption capacities of Cu, Zn, Hg and Cd increased with increasing pH, which may be because the increase in pH effectively reduced the concentration of H + in the solution. However, soil colloids with a negative charge and heavy metal ions with a positive charge can more effectively increase the adsorption capacity of soil colloids for heavy metal ions; therefore, the adsorption capacities of heavy metal ions increased with increasing pH value. Surprisingly, the capacities of Pb and As were saturated at the beginning of adsorption (Figure 9a). Generally, an increase in pH can promote the adsorption of several heavy metals and inhibit the adsorption of Cr; however, its mechanism still needs to be further investigated. The maximum adsorption capacities of heavy metals in the soil sample were reached at pH = 8.00, which may be related to the physical and chemical properties of the soil material. It was found that the slow adsorption processes emerged at pH < 7.00, implying that the adsorption process of heavy metals in the soil sample was slow under acidic conditions. In contrast, the fast adsorption stage occurred at pH > 7.00, indicating that the adsorption process of heavy metals in the soil was fast under alkaline conditions. The influence of pH on the adsorption of heavy metals by soil particles was generally shown in the following four aspects: the influence on the solubility of heavy metals in solution [48], the influence on the adsorption characteristics of the natural colloid surface in solid particles [49], control of various adsorption reactions on the surface of solid particles, and the increase in pH, which promoted an increase in the adsorption point of the soil colloid since soil colloids generally have a negative charge [50]. Figure 9b exhibits the change in adsorption capacities followed by the different electrolyte concentrations (10 −5 , 10 −4 , 10 −3 , 10 −2 and 10 −1 mol/L). The adsorption capacities of heavy metals increased with increasing electrolyte concentrations, which may be due to the promotion of the ion activity of heavy metals in the electrolyte. Several heavy metals showed the phenomenon of decreasing first and increasing later with the increase in NaNO 3 concentration, which may be because the inert electrolyte was neutralized by the external negative charge of the soil, reducing the affinity towards heavy metals and the adsorption capacity of heavy metals. Generally, there are three approaches for a soil to adsorb heavy metals: the activity of free heavy metal ions will change because of the effect of ions or the pH by forming ions [51]; inert electrolytes will dissociate into anions and cations, and cations will compete with excessive heavy metal ions or heavy metal ions for adsorption [52]; inert electrolytes will affect the electrostatic point on the surface of the soil colloidal particles since the soil colloid is charged [53]. towards heavy metals and the adsorption capacity of heavy metals. Generally, there are three approaches for a soil to adsorb heavy metals: the activity of free heavy metal ions will change because of the effect of ions or the pH by forming ions [51]; inert electrolytes will dissociate into anions and cations, and cations will compete with excessive heavy metal ions or heavy metal ions for adsorption [52]; inert electrolytes will affect the electrostatic point on the surface of the soil colloidal particles since the soil colloid is charged [53].

Conclusions
Topsoil samples collected from Lixisol of Guizhou Province were selected in this study, and the adsorption behavior of heavy metals on the soil samples was explored by indoor simulation. The contents of Cu, Zn, Cd, Cr, Pb, Hg and As in the soil material were 10.8 mg/kg, 125 mg/kg, 0.489

Conclusions
Topsoil samples collected from Lixisol of Guizhou Province were selected in this study, and the adsorption behavior of heavy metals on the soil samples was explored by indoor simulation. The contents of Cu, Zn, Cd, Cr, Pb, Hg and As in the soil material were 10.8 mg/kg, 125 mg/kg, 0.489 mg/kg, 23.5 mg/kg, 22.7 mg/kg, 58.3 mg/kg and 45.4 mg/kg, respectively. The soil pH values increased with the increase of CaCO 3 concentration in the soil, and the fluctuation of the soil pH values was weak after the CaCO 3 concentration reached 100 g/kg. The adsorption of lime soil increased by approximately 10 mg/kg on average, and the desorption capacity decreased by approximately 300 mg/kg on average. The desorption of all heavy metals in this study did not change with increasing clay content.
Pseudo-second-order kinetics were more suitable to describe the adsorption kinetics of heavy metals on the soil sample, as evidenced by the higher R 2 value. Furthermore, the Freundlich model can better describe the adsorption process of As on lime soil. The process of As, Cr, Cd and Hg adsorption onto the soil material was spontaneous and entropy-driven. Additionally, the process of Cu and Pb adsorption onto the soil sample was spontaneous and enthalpy-driven. It was found that the slow adsorption processes occurred at pH < 7.00, implying that the adsorption process of heavy metals in the soil was slow under acidic conditions. In contrast, the fast adsorption stage occurred at pH > 7.00, indicating that the adsorption process of heavy metals in the soil was fast under alkaline conditions. Generally, the adsorption and desorption of heavy metals in polluted soil increased and decreased, respectively, with increasing CaCO 3 content. The effect of calcium carbonate on the accumulation of heavy metals in soil was greater than that of clay.

Conflicts of Interest:
The authors declare no competing interests.