Effect of Soil Solution Properties and Cu2+ Co-Existence on the Adsorption of Sulfadiazine onto Paddy Soil

Paddy soils are globally distributed and saturated with water long term, which is different from most terrestrial ecosystems. To better understand the environmental risks of antibiotics in paddy soils, this study chose sulfadiazine (SDZ) as a typical antibiotic. We investigated its adsorption behavior and the influence of soil solution properties, such as pH conditions, dissolved organic carbon (DOC), ionic concentrations (IC), and the co-existence of Cu2+. The results indicated that (1) changes in soil solution pH and IC lower the adsorption of SDZ in paddy soils. (2) Increase of DOC facilitated the adsorption of SDZ in paddy soils. (3) Cu2+ co-existence increased the adsorption of SDZ on organic components, but decreased the adsorption capacity of clay soil for SDZ. (4) Further FTIR and SEM analyses indicated that complexation may not be the only form of Cu2+ and SDZ co-adsorption in paddy soils. Based on the above results, it can be concluded that soil solution properties and co-existent cations determine the sorption behavior of SDZ in paddy soils.


Introduction
Overuse and the uncontrolled disposal of antibiotics have caused severe environmental problems, especially for the soil environment [1][2][3]. As one of the main crops in the world, rice is widely planted and feeds the majority of the world's population; while over 92% of rice production, as a primary staple food, is in Asia [4]. In China, the planting area of rice is near 30 million ha, which occupies over 20% of the total farmlands [5]. Therefore, paddy soils have a high chance of exposure to antibiotic pollutants, such as sulfonamides [1]. However, most studies focused on the adsorption and transport behaviors of antibiotics in upland farmland, while few investigations were about the behavior of antibiotics in paddy soils.
The co-existence of heavy metals and antibiotic pollution in the soil environment is receiving more and more attention. With their wide antibacterial spectrum, excellent curative effect, and low cost, sulfonamide antibiotics (SAs) are widely used in disease treatment and prevention for humans and animals [6]. Sulfadiazine (SDZ), one of the most commonly used SA chemicals, has been widely detected in various environmental mediums, especially in soils [1,7,8]. Previous studies proved that Cu ions could coordinate with SDZ [9][10][11], altering their molecular speciation and environmental behavior [11,12]. For example, it was reported that Cu 2+ co-addition could improve the adsorption of antibiotics to organic matters [13] and soils [14]. On the other hand, there was competitive adsorption between co-existent metal ions and the ionizable antibiotics, which inhibited the adsorption ability of antibiotics on soil components [15,16] and increase their transportation risk in the environment [17]. However, limited investigations have been conducted on the effect of Cu 2+ on SA sorption in paddy soils.
To obtain more knowledge about the environmental risk of SAs in paddy soils, we chose Sulfadiazine (SDZ) as a representative antibiotic. The batch sorption experiments were conducted in soil suspensions with different pH, ion concentrations (IC), and dissolved organic carbon (DOC), as well as co-existent Cu 2+ . The aim of this study was: (1) to explore the sorption kinetics and isotherm of SDZ in paddy soils under different soil solution conditions; (2) to characterize the sorption behavior and mechanisms of SDZ on various soil components; and (3) to investigate the effects of co-existent Cu 2+ on the sorption of SAs in paddy soils.

Paddy Soil Sampling and Characterization
Five soil samples were collected from the rice production area located in Huang gang, Hubei, of China. About 20 kg of topsoil (0-10 cm) was collected from sample sites occupying an area of approximately 500 m 2 and transported to the lab in a cabinet with an ice pack. After thorough mixing, soils were air-dried, ground, and sieved through a 2-mm sieve for the following soil property analysis and sorption batch experiments. The soil texture was silty loam with 2.55% clay, 89.67% silt, and 7.79% sand. Soil pH, organic matter (OM), cation exchange capacity (CEC), and Zeta potential are listed in Table 1.

Preparation of Soil Organic Particles and Soil Clay
In this study, the soil organic particle fraction was obtained with the wet sieving method suggested by Elliott [19]. Briefly, put 500 g of air-dried soil on a sieve of 2 mm, and separate aggregates by shaking the sieve up and down for 2 min with 50 repetitions, after being submerged in water for 5 min. Aggregates passed through a 0.25-mm sieve were collected and designated as a soil organic particle fraction. The physical and chemical properties of the prepared organic particles are listed in Table 1.
The soil clay fraction was extracted with 0.1 M Na 4 P 2 O 7 (pH = 7) and 0.1 M H 2 O 2 according to [20]. After being treated with Na 4 P 2 O 7 (Sigma, Shanghai, China) and 0.1 M H 2 O 2 (Sinopharm, Shanghai, China), the sample mixture was centrifuged and separated into three distinct layers. After removing the top layer using suction, the residual material was separated into two fractions by repeat centrifuging. After removing the top liquid layer, recover the clay mineral fraction of 1 M NaCl (Sinopharm, Shanghai, China), after repeat centrifuging. The physical and chemical properties of the prepared clay are listed in Table 1.

Batch Sorption Experiment
Sorption experiments were conducted in 100-mL plastic centrifuge tubes with Teflon lids, according to the OECD 106 method and Jiang [21]. Briefly, weigh 1.000 g dry paddy soil into 25 mL of 0.01 M CaCl 2 (Sinopharm, Shanghai, China) containing a specific concentration of SDZ. All sorption experiments were carried out with an oscillator ( Changzhou Guohua Electric Appliance Co., Ltd., Changzhou, China) at 150 rpm in darkness. The sample tubes containing an initial concentration of 8 µmol L −1 SDZ were shaken from 0 h to 72 h for the kinetic sorption experiment, and those containing a series of initial concentrations of 4 µmol L −1 , 8 µmol L −1 , 12 µmol L −1 , and 16 µmol L −1 were shaken for 72 h for the isotherm experiment. At the end of the sorption experiments, 2 mL supernatant was collected from each tube and centrifuged at 10000 r for 1 min. Then, the supernatant was filtered through a 0.45 µm Whatman filter (Whatman TM , Germany) for the subsequent quantification with HLPC. Each experiment was conducted in triplicates, and a blank treatment of SDZ solution without soil was used as a control, which was used to evaluate the loss of SDZ caused by sorption onto the tube walls and its degradation.
Environmental factors: The pH of soil solution was adjusted to 5.0, 7.0, and 9.0, respectively, with 1 M HCl (Sinopharm, Shanghai, China) and NaOH (Sinopharm, Shanghai, China) to analyze the influence of the medium pH. Fulvic acid (FA, ≥90%, CAS: 1415-93-6, Aladdin, Shanghai, China) was added at rates of 1 g L −1 , 3 g L −1 , and 5 g L −1 to the soil solution to analyze the effect of soil DOC on the adsorption of SDZ on soils. At the same time, the effect of solution ionic concentration was evaluated with Ca 2+ ions at the concentrations of 0.05 M and 0.1 M. The sorption characteristics of SDZ on various soil fractions was tested with the same procedure as the above soils, and the solid phase was replaced with the soil organic particles and clays prepared in 1.3.
Effect of co-existent Cu 2+ : For the sorption experiments of SDZ with Cu 2+ co-existence, the added concentrations of Cu 2+ were 200 and 500 mg L −1 , respectively. The others were the same as the soil sorption experiment.

Analysis Method
The pH, CEC, and OM of the soil samples were analyzed according to the description of Lu (1999) [22]. The zeta potentials of the crude soil, organic particle fraction, and soil clay fraction were measured with a zeta potential analyzer (JS94H, Shanghai Zhongchen Digital Technic Apparatus Co., Ltd., Shanghai, China) after being dispersed in a solution of 0.01 M NaCl (Sinopharm, Shanghai, China).
The concentration of SDZ was increased to 40 µmol L −1 in the soil solutions and they were then centrifuged after the adsorption reached equilibrium. The centrifuged samples were washed three times with ultra-water, dried, and powdered with a mortar and pestle for the subsequent scanning with FTIR spectrum and SEM visualization. The FTIR scanning was conducted with a FTIR Microscope-Spectrometer (Nicolet5700, Thermo, Shanghai, China) in the range 4000-400 cm −1 . The microstructure observation of soil samples was conducted with a SEM (Quanta250, FEI, Shanghai, China) at magnifications of ×20,000.

Data Calculation
The amount of adsorbed SDZ was calculated using the following equation: where Y is the absorbed amount of SDZ (µmol kg −1 ); V (L) is the volume of soil solution; C 0 and C e (µmol L −1 ) represent the concentration of SDZ at the beginning and end of the equilibrium sorption experiments in solution, respectively; and M (kg) is the dry soil weight added to the background solution.
The sorption kinetics of SDZ in paddy soils in various soil solutions were fitted with the following kinetic model equations: The sorption isotherm of SDZ in paddy soils in various soil solutions were fitted with the Linear model equation (4) and Freundlich model equation (5): For Equations (2)-(5), q t (µmol kg −1 ) is the amount of SDZ sorption to paddy soil at time t (h); q e (µmol kg −1 ) is the amount of SDZ sorption to paddy soil when sorption reaches the equilibrium; and k 1 /k 2 is the constant of the kinetics sorption velocity. C e (µmol L −1 ) is the concentration of SDZ in the soil supernatant when sorption reaches equilibrium; K d (L kg −1 ) is the coefficient of SDZ distribution between the liquid and solid phases in the equilibrium system; K F is the Freundlich sorption coefficient; and n is the nonlinearity factor.

Effect of Soil Solution Properties on the Sorption of SDZ
Batch sorption experiments were conducted to analyze the sorption characteristics of SDZ on paddy soils under different experimental systems. The results are shown in Figure 1. REVIEW 4 The amount of adsorbed SDZ was calculated using the following equation:

Int. J. Environ. Res. Public Health 2021, 18, x FOR PEER
where Y is the absorbed amount of SDZ (μmol kg -1 ); V (L) is the volume of soil tion; C0 and Ce (μmol L −1 ) represent the concentration of SDZ at the beginning and e the equilibrium sorption experiments in solution, respectively; and M (kg) is the dry weight added to the background solution.
The sorption kinetics of SDZ in paddy soils in various soil solutions were fitted the following kinetic model equations: The sorption isotherm of SDZ in paddy soils in various soil solutions were fitted the Linear model equation (4) and Freundlich model equation (5): , qt (μmol kg -1 ) is the amount of SDZ sorption to paddy s time t (h); qe (μmol kg -1 ) is the amount of SDZ sorption to paddy soil when sorption rea the equilibrium; and k1/k2 is the constant of the kinetics sorption velocity. Ce (μmol L the concentration of SDZ in the soil supernatant when sorption reaches equilibrium; kg -1 ) is the coefficient of SDZ distribution between the liquid and solid phases in the librium system; KF is the Freundlich sorption coefficient; and n is the nonlinearity fa

Effect of Soil Solution Properties on the Sorption of SDZ
Batch sorption experiments were conducted to analyze the sorption characterist SDZ on paddy soils under different experimental systems. The results are shown in F 1.  As plotted in Figure 1b, the adsorption capacity of SDZ in crude soils was decre with the pH changing to 5.0, 7.0, and 9.0, and qe was over 60 μmol kg -1 in crude so lower than 50 μmol kg -1 with the pH change ( Table 2). As for the ion strength, the ad tion capacity of SDZ on paddy soils was significantly reduced, to about 20 μmol kg -1 , the CaCl2 concentration increased to 0.05 M and 0.1 M CaCl2 (Figure 1d). The variati the fitting parameters of the adsorption isotherms (Table 3) indicated that the soil sol pH and ion strength could change the sorption affinity of SDZ in paddy soil. How there was a significant increase in SDZ adsorption capacity with increasing soil sol DOM ( Figure 1f and Table 2). In addition, the values of K1 and K2 were enhanced wi increasing of pH and DOM (Table 2), which indicated that the increase of pH and accelerated the reaction process and increased the velocity of SDZ adsorption on p soil (Figure 1a,1c and 1e). Whatever the changes in soil solution properties in the pr study, the adsorptions of SDZ had a good fit with the pseudo-first, pseudo-second-o and the Freundlich equation, as the R parameter of each selected fitting model was al over 0.99 at a level of p < 0.05 (list in Tables 2 and 3). This might suggest that there more than one dominant mechanism responsible for the adsorption of SDZ. Table 2. Fitting results of kinetics data to pseudo-first and pseudo-second-order equations fo adsorption in various soil solutions (mean values ± standard error).

Treatments
The  As plotted in Figure 1b, the adsorption capacity of SDZ in crude soils was decreased, with the pH changing to 5.0, 7.0, and 9.0, and q e was over 60 µmol kg −1 in crude soil but lower than 50 µmol kg −1 with the pH change ( Table 2). As for the ion strength, the adsorption capacity of SDZ on paddy soils was significantly reduced, to about 20 µmol kg −1 , when the CaCl 2 concentration increased to 0.05 M and 0.1 M CaCl 2 (Figure 1d). The variation in the fitting parameters of the adsorption isotherms (Table 3) indicated that the soil solution pH and ion strength could change the sorption affinity of SDZ in paddy soil. However, there was a significant increase in SDZ adsorption capacity with increasing soil solution DOM (Figure 1f and Table 2). In addition, the values of K 1 and K 2 were enhanced with the increasing of pH and DOM (Table 2), which indicated that the increase of pH and DOM accelerated the reaction process and increased the velocity of SDZ adsorption on paddy soil (Figure 1a,c,e). Whatever the changes in soil solution properties in the present study, the adsorptions of SDZ had a good fit with the pseudo-first, pseudo-second-order, and the Freundlich equation, as the R parameter of each selected fitting model was always over 0.99 at a level of p < 0.05 (list in Tables 2 and 3). This might suggest that there was more than one dominant mechanism responsible for the adsorption of SDZ.

Effect of Cu 2+ Co-Existing on the Sorption of SDZ
The adsorption kinetics and isotherms of SDZ were compared for the crude soil, organic particles, and clay with the presence or absence of Cu 2+ (Figure 2). Whether Cu 2+ was present or not, the adsorption of SDZ on the three soil fractions could reach equilibrium within about 48 h (Figure 2a,c,e). However, with the co-existence of Cu 2+ , the adsorption capacity of SDZ on the crude soil and the organic fraction was significantly increased, which did not happen with the clay (Figure 2f). The values of K F and n were also enhanced with Cu 2+ co-addition for the crude soils and organic particles. Morover, with the co-existence of Cu 2+ , the adsorption characteristics of SDZ on each soil fraction still fit well with the selected kinetic equations and isotherm equations, as most of the coefficients of R were >0.99 (Tables 4 and 5), although the adsorption of SDZ was decreased for clay soils (Table 5). Moreover, the subsequent linear model fitting also found that the adsorption affinity (k d , Table 5) of clay to SDZ was much lower than the crude soil and organic particles, and the presence of Cu further decreased this adsorption affinity. Nevertheless, this was increased for the crude soil and organic particles.

Effect of Cu 2+ Co-Existing on the Sorption of SDZ
The adsorption kinetics and isotherms of SDZ were compared for the crude soil, organic particles, and clay with the presence or absence of Cu 2+ (Figure 2). Whether Cu 2+ was present or not, the adsorption of SDZ on the three soil fractions could reach equilibrium within about 48 h (Figures 2a, 2c and 2e). However, with the co-existence of Cu 2+ , the adsorption capacity of SDZ on the crude soil and the organic fraction was significantly increased, which did not happen with the clay (Figure 2f). The values of KF and n were also enhanced with Cu 2+ co-addition for the crude soils and organic particles. Morover, with the co-existence of Cu 2+ , the adsorption characteristics of SDZ on each soil fraction still fit well with the selected kinetic equations and isotherm equations, as most of the coefficients of R were >0.99 (Tables 4 and 5), although the adsorption of SDZ was decreased for clay soils (Table 5). Moreover, the subsequent linear model fitting also found that the adsorption affinity (kd, Table 5) of clay to SDZ was much lower than the crude soil and organic particles, and the presence of Cu further decreased this adsorption affinity. Nevertheless, this was increased for the crude soil and organic particles.   13.76 ± 0.14 0.66 ± 0.14 0.9998 15.08 ± 0.13 0.22 ± 0.05 0.9999     Figure 3 shows the FTIR spectra of soil OC and clay particles with/without SDZ/Cu adsorption. Compared to the treatment without SDZ, attenuation was observed at wavenumbers from 1050 to 950 cm −1 (C-O-C) on the OC particles after SDZ adsorption. Compared to the single adsorption of SDZ, a stretching vibration in the range 1700-1600 cm −1 (C=C or -NH) in OC components became stronger with Cu co-addition. These results suggest that the -NH, C-O-C, or C=C groups play an important role in the SDZ adsorption to soil particles [23]. For the clays, a weak peak, ranging from 3000 cm −1 to 2900 cm −1 (-CH), occurred after SDZ adsorption, and this signal became stronger in clay particles with the co-presence of Cu. Moreover, an attenuation of the wavenumbers from 1700 to 1600 cm −1 and an obvious -OH stretching vibration (3678.0 cm −1 ) were observed for the Cu-SDZ co-adsorption, compared to the single adsorption of SDZ. This might suggest that hydrogen bonding plays an essential role in the co-adsorption of Cu and SDZ on soil clays. 1.00 ± 0.15 0.8920 8.57 ± 0.35 1.50 ± 0.14 0.9445 Figure 3 shows the FTIR spectra of soil OC and clay particles with/without SDZ/Cu adsorption. Compared to the treatment without SDZ, attenuation was observed at wavenumbers from 1050 to 950 cm -1 (C-O-C) on the OC particles after SDZ adsorption. Compared to the single adsorption of SDZ, a stretching vibration in the range 1700-1600 cm -1 (C=C or -NH) in OC components became stronger with Cu co-addition. These results suggest that the -NH, C-O-C, or C=C groups play an important role in the SDZ adsorption to soil particles [23]. For the clays, a weak peak, ranging from 3000 cm -1 to 2900 cm -1 (-CH), occurred after SDZ adsorption, and this signal became stronger in clay particles with the co-presence of Cu. Moreover, an attenuation of the wavenumbers from 1700 to 1600 cm -1 and an obvious -OH stretching vibration (3678.0 cm -1 ) were observed for the Cu-SDZ coadsorption, compared to the single adsorption of SDZ. This might suggest that hydrogen bonding plays an essential role in the co-adsorption of Cu and SDZ on soil clays. The surface morphologies of the soil OC and clay particles were scanned using SEM (Figure 4). The SEM images show that soil OC particles surface was rougher than the clay soils, with an irregular and bulky shape. However, there was no obvious change in the surface morphology of the OC and clay particles with the sorption of Cu and SDZ. According to the mapping results, the Cu was very hard to detect in OC and clay without Cu 2+ treatment, but it covered the surface of the samples treated with Cu 2+ , with the signal becoming stronger for the combined organic soil treated with SDZ. This indicated that Cu 2+ adsorbs easily on organic soil and clays, and SDZ co-addition increases the sorption The surface morphologies of the soil OC and clay particles were scanned using SEM (Figure 4). The SEM images show that soil OC particles surface was rougher than the clay soils, with an irregular and bulky shape. However, there was no obvious change in the surface morphology of the OC and clay particles with the sorption of Cu and SDZ. According to the mapping results, the Cu was very hard to detect in OC and clay without Cu 2+ treatment, but it covered the surface of the samples treated with Cu 2+ , with the signal becoming stronger for the combined organic soil treated with SDZ. This indicated that Cu 2+ adsorbs easily on organic soil and clays, and SDZ co-addition increases the sorption of Cu on organic soil. From the above results of the FTIR analysis and the batch adsorption experiments, it was indicated that co-adsorption by the complexation of Cu 2+ -SDZ is one method of Cu 2+ and SDZ adsorption to soil particles. While the co-adsorption of Cu 2+ and SDZ might mainly be conducted by their independent adsorption to different sites on the soil particles. becoming stronger for the combined organic soil treated with SDZ. This indicated that Cu 2+ adsorbs easily on organic soil and clays, and SDZ co-addition increases the sorption of Cu on organic soil. From the above results of the FTIR analysis and the batch adsorption experiments, it was indicated that co-adsorption by the complexation of Cu 2+ -SDZ is one method of Cu 2+ and SDZ adsorption to soil particles. While the co-adsorption of Cu 2+ and SDZ might mainly be conducted by their independent adsorption to different sites on the soil particles.

Effects of Soil Solution Properties on the Sorption of SDZ in Paddy Soils
In addition to the basic properties of antibiotics, soil physical and chemical properties also play an essential role in the sorption of antibiotics to soil constituents [24]. In this study, the adsorption kinetics and isotherms processes of SDZ were all well fitted (R > 0.99) by the first/second-order equations and the Freundlich equation, whatever the variation of the soil properties. This indicated that the variation of soil solution properties did not change the process and mechanism of SDZ adsorption on paddy soils, which might be ascribed to the hydrophobic distribution that contributed to the sorption of SDZ on paddy soil [21,25], in view of the chemical character of SDZ, which was amphoteric and weak acid polar, with little water-solubility. It was also proven that the adsorption process of SDZ in soils was driven by weak hydrophobic forces, as in neutral or anionic specie of natural soil. The present study also indicated further that the SDZ chemicals were easily adsorbed to the soil organic matter, as the SDZ adsorption increased with soil organic matter addition (Figure 1f), and there was more obvious signal noise for soil OC constituents than the clay after being treated with SDZ ( Figure 3). These results indirectly suggests that hydrophobic distribution has a beneficial effect on the adsorption of SDZ in paddy soils.
The pH determined whether SDZ existed as cations, zwitterions, or anions in the soils [26]. Hence, SDZ exhibited pH-dependent adsorption on soil constituents [27,28]. The variation of soil solution pH significantly decreased the adsorption capacity of SDZ on paddy soils, based on the adsorption coefficients of SDZ (Table 3); which is similar to the previous results, that the increase of pH from 4 to 8 significantly decreased adsorption coefficients of SAs from 30 to 1 [29]. This might be explained by the fact that the disturbance of

Effects of Soil Solution Properties on the Sorption of SDZ in Paddy Soils
In addition to the basic properties of antibiotics, soil physical and chemical properties also play an essential role in the sorption of antibiotics to soil constituents [24]. In this study, the adsorption kinetics and isotherms processes of SDZ were all well fitted (R > 0.99) by the first/second-order equations and the Freundlich equation, whatever the variation of the soil properties. This indicated that the variation of soil solution properties did not change the process and mechanism of SDZ adsorption on paddy soils, which might be ascribed to the hydrophobic distribution that contributed to the sorption of SDZ on paddy soil [21,25], in view of the chemical character of SDZ, which was amphoteric and weak acid polar, with little water-solubility. It was also proven that the adsorption process of SDZ in soils was driven by weak hydrophobic forces, as in neutral or anionic specie of natural soil. The present study also indicated further that the SDZ chemicals were easily adsorbed to the soil organic matter, as the SDZ adsorption increased with soil organic matter addition (Figure 1f), and there was more obvious signal noise for soil OC constituents than the clay after being treated with SDZ ( Figure 3). These results indirectly suggests that hydrophobic distribution has a beneficial effect on the adsorption of SDZ in paddy soils.
The pH determined whether SDZ existed as cations, zwitterions, or anions in the soils [26]. Hence, SDZ exhibited pH-dependent adsorption on soil constituents [27,28]. The variation of soil solution pH significantly decreased the adsorption capacity of SDZ on paddy soils, based on the adsorption coefficients of SDZ (Table 3); which is similar to the previous results, that the increase of pH from 4 to 8 significantly decreased adsorption coefficients of SAs from 30 to 1 [29]. This might be explained by the fact that the disturbance of the acid-base balance of crude soils would increase the electrostatic repulsion between anionic SDZand the negatively charged soil surface [25,30], and by the lower lipophilic interactions between the uncharged chemical molecules and soil particles [31]. As for the ion strength factor, a lower adsorption potential of SDZ was found with increasing Ca 2+ concentration (Figure 1), which was dissimilar to the previous reports that multivalence ions can improve antibiotic adsorption by covalent bonding action [32]. This inconsistent result might indicate that a high ion strength would weaken the function of the ion-bridge for SDZ adsorption on paddy soils, as it has a higher OM content (Table 1).

Effect of Cu 2+ Co-Existing on the Sorption of SDZ
With increased detection of heavy metals and antibiotics in soils, more attention has been paid to the process of antibiotics related to the co-existent metal ions, aside from the various soil properties [33][34][35]. It is widely thought that heavy metal cations may influence the mobility of antibiotics in soils, through ion exchange or complexation. In the present study, the adsorption isotherms of SDZ were well described by the Freundlich model (R = 0.9445-1.0000) in a binary system of Cu-SDZ. Moreover, both the adsorption capacity and affinity of SDZ were augmented, as Cu 2+ was increased, based on the changing trend of k F and n values (Table 5). These results indicated that the presence of co-existent Cu 2+ did not change the mechanism of SDZ adsorption to soil particles but could alter the adsorption capacity of SDZ to soil particles with different compositions.
The further FTIR and SEM analyses also detected the co-adsorption of Cu and SDZ on soil constituents (Figures 3 and 4). It was suspected that a Cu 2+ bridge between soil particles and SDZ [11,36] or the formation of SDZ-Cu complexes [37,38] occurred for the co-adsorption of Cu and SDZ on paddy soils. However, this might be mainly carried out by their independent adsorption to different sites of the soil particles, as the weak signal changes of the FTIR spectrum were found primarily in soil OC constituents after co-treating with Cu and SDZ ( Figure 3). It was reported that the presence of heavy metal ions could result in changes in soil properties [39][40][41], which might shape the sorption of antibiotics on soils. Therefore, the decrease of SDZ adsorption on clay in the present study might be attributed to the precipitation of soil mineral colloid induced by the presence of Cu 2+ [39], which decreased the size or blocked the pores of the soil clay particles [42].

Conclusions
The change of soil solution pH and ion concentration lowered the SDZ adsorption on paddy soil, but it was increased with the addition of organic matter. Whatever the changes of soil properties, the sorption kinetics of SDZ on paddy soils could be well described by the pseudo-first and second-order equations, and the Freundlich equation could fit the sorption isotherms well. Cu 2+ coexistence in the soil increased the adsorption of SDZ on crude soil and its organic components, but decreased the adsorption capacity of clay soil for SDZ. Based on a further analysis of FTIR and SEM, it could be concluded that the coadsorption of Cu 2+ and SDZ on soil constituents might be conducted by the complexation of Cu 2+ and SDZ, in addition to their independent adsorption to different active sites of the soil particles.