Next Article in Journal
Mechanical Strength and Hydration Characteristics of Cement Mixture with Highly Concentrated Hydrogen Nanobubble Water
Next Article in Special Issue
Medical Plant Extract Purification from Cadmium(II) Using Modified Thermoplastic Starch and Ion Exchangers
Previous Article in Journal
Nanostructured Surface Finishing and Coatings: Functional Properties and Applications
Previous Article in Special Issue
Regular Polymeric Microspheres with Highly Developed Internal Structure and Remarkable Thermal Stability
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sorption Properties of Specific Polymeric Microspheres towards Desethyl-Terbuthylazine and 2-Hydroxy-Terbuthylazine: Batch and Column Studies

Faculty of Chemistry, Wrocław University of Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
*
Author to whom correspondence should be addressed.
Materials 2021, 14(11), 2734; https://doi.org/10.3390/ma14112734
Submission received: 28 April 2021 / Revised: 18 May 2021 / Accepted: 20 May 2021 / Published: 22 May 2021
(This article belongs to the Special Issue Structural and Thermal Properties of Polymeric Microspheres)

Abstract

:
This work investigates the sorption properties of poly(divinylbenzene) modified in the Diels–Alder reaction towards persistent and mobile metabolites of terbuthylazine. The batch experiments were carried out to examine the efficiency of desethyl-terbuthylazine and 2-hydroxy-terbuthylazine adsorption on the specific adsorbent and the impact of different factors on the adsorption process. Results fit well to a pseudo-second order kinetic model. It was confirmed that hydrogen bonds play an important role in the studied systems. Five times greater sorption of 2-hydroxy-terbuthylazine than desethyl-terbuthylazine was observed. The molecular structures of both metabolites exhibit complementarity to the arrangement of functional groups in the polymer but the differences in the physicochemical properties of the desethyl derivative make it a highly mobile compound with higher affinity to the aqueous phase. The equilibrium data in the batch study fit the Freundlich isotherm for 2-hydroxy-terbuthylazine, and for desethyl-terbuthylazine the Temkin and Dubinin–Radushkevich models were better. The adsorption capacities obtained under dynamic conditions were comparable with batch results. For column adsorption modeling the Bohart–Adams, Wolborska, Thomas and Yoon–Nelson models were used. The proposed microspheres can be reused easily with no significant decrease in adsorption capacity by using ethanol as eluent in the desorption.

Graphical Abstract

1. Introduction

Herbicides, as well as their degradation products, are toxic substances that can easily remain for many years in various environmental compartments due to their high chemical and biochemical stability in the environment. Agricultural soils are the first recipients of herbicides after their application. When plant protection products bind strongly to soil, which is mainly the case with hydrophobic and cationic compounds, their persistence there increases, while biodegradability becomes lower. On the other hand, the low adsorption of herbicides in soil greatly facilitates their migration to surface waters and groundwaters [1,2,3]. Therefore, herbicides are capable of creating residues that can be detected in the environment even after a long time [1].
Contaminants that are poorly degradable and easily moved to a wider area by spreading through the environment (groundwater, surface waters, bottom sediments) are defined as persistent and mobile compounds [4,5,6]. According to the first pan-European reconnaissance of the occurrence of polar organic persistent pollutants in European groundwater [7], terbuthylazine and desethyl-terbuthylazine were among the pollutants that exceeded the European Union (EU) groundwater standard of 0.1 μg L−1 [8]. A study was conducted for 59 selected compounds, including pharmaceuticals, pesticides and their transformation products, perfluorinated acids (PFAs), benzotriazoles and others. In total, 164 samples from 23 countries were analyzed, detecting the desethyl-terbuthylazine in 49% of samples (maximum concentration was 266 ng L−1) and the terbuthylazine in 34% of samples (maximum concentration was 716 ng L−1).
Terbuthylazine became one of the most widely used herbicides in crops protection, mainly due to the ban on the use of atrazine in Europe from 2004. It is known that this substance presents a potential risk due to its combination of persistence and toxicity to living organisms. There is a risk, that long-term exposure of aquatic organisms, mammals and soil invertebrates to high concentrations of this pesticide may cause negative effects such as DNA damage and carcinogenic effects [9,10]. Moreover, some studies show that the toxicity of the herbicide metabolites may be higher than their parent compound [11,12]. TBA mineralization and transformation rely on two processes—chlorohydrolysis, with associated depletion of the chlorine atom or N-dealkylation of the side chains of the s-triazine ring. The main degradation pathways of TBA in the soil are presented in [13]. Depending on the microbial activity, metabolites such as desethyl-terbuthylazine, hydroxy-terbuthylazine, desethylhydroxy-terbuthylazine, diamino-chlorotriazine), desbutylhydroxy-terbuthylazine, amino-dihydroxy-triazine may be formed. Desethyl-terbuthylazine and 2-hydroxy-terbuthylazine are the major products of degradation of terbuthylazine in soil [14] and the desethyl-terbuthylazine is the most detected terbuthylazine metabolite in surface waters [15,16]. Few studies show the negative effects of desethyl-terbuthylazine and 2-hydroxy-terbuthylazine on aquatic organisms [17,18,19]. However, as classification with R40 (carcinogen category 3) is proposed for terbuthylazine, the groundwater metabolites are considered toxicologically relevant [20].
The pesticide transformation products often show higher persistence and mobility than the active substance [21]. The persistence in the soil for terbuthylazine is medium to high—the period required for 50% dissipation (DT50) is 65–167 days (20 °C, soil moisture 13–36% w/w). DT50 for desethyl-terbuthylazine is 27–113 days (20 °C, soil moisture 11–29% w/w) and for 2-hydroxy-terbuthylazine is 207 to >1000 days (20 °C, soil moisture 11–29% w/w), so they exhibit moderate to high persistence and high to very high persistence, respectively [22]. Regarding mobility, 2-hydroxy-terbuthylazine and terbuthylazine are classified as medium mobile, while desethyl-terbuthylazine is highly mobile [22]. Therefore, terbuthylazine metabolites as persistent and mobile chemicals are frequently detected in the environmental compartments in Europe, in many cases at concentrations higher than the European drinking water limit according to Directive (EU) 2020/2184 [23] equal to 0.1 µg L−1, e.g., in the range of up to 9.24 μg L−1 in surface waters for 2-hydroxy-terbuthylazine and 2.90 μg L−1 and 3.50 μg L−1 for desethyl-terbuthylazine in surface water and groundwater, respectively [24,25,26,27]. Desethyl-terbuthylazine was also often found in tap water, for example in samples collected in France [28] and Croatia [29]. Regarding the above, Tasca and Fletcher signal that the range of emerging contaminants detected in drinking waters, and their seasonal and geographic variability, demands more accurate screening programs, and greater attention has to be focused on the fate of herbicide metabolites [30]. This is especially important in the case of toxic substances. However, it must be noted that also contamination with less toxic but persistent and mobile substances can become biologically relevant. Once released into the environment, these chemicals may spread easily to a wider area of the environment or adsorb on the bottom sediments of reservoirs or soil particles. On the one hand, adsorption is beneficial because it reduces the bioavailability of these compounds so that they are not readily absorbed by the organisms of a given ecosystem. On the other hand, adsorbed substances are accumulated and gradually released into water bodies. Therefore, developing the aquatic system and wastewater monitoring as the main source of persistent and mobile substances is one of the European Union’s objectives, which allows the risk for human health to be limited and the life cycle assessment of these substances in the environment to be determined.
However, a substantial analytical challenge exists related to the detection and quantification of persistent and mobile chemicals in the environmental samples. It is very important, especially for those substances that have not previously been covered in the routine monitoring, which also include desethyl-terbuthylazine and 2-hydroxy-terbuthylazine. As long as persistence and toxicity are not of great importance regarding chemical analysis, the mobility significantly affects analytical processes applied in the trace analysis of organic contaminants that are present in the water environment [31].
The analysis of trace impurities is preceded by the stage of preconcentration of targeted compound, most often with the use of the solid-phase extraction (SPE) method, using various types of adsorbent. It should also be emphasized that the persistent and mobile organic compounds are more of a concern for water as they are difficult to capture by sorption processes due to their high polarity [5]. Herbicide metabolites detection is also difficult due to their low concentrations in the environmental compartments (ng L−1) and much higher concentrations (mg L−1) of other pollutants occurring in environmental samples. In order to overcome all these difficulties, the specific sorptive materials for SPE should be used, the effectiveness of which will be increased by selective interactions with the functional groups of individual types/groups of persistent and mobile chemicals. Terbuthylazine metabolites belong to the group of triazines having both hydrogen acceptor and hydrogen donor, indicating a high potential for the formation of hydrogen bonds which, besides van der Waals interactions, play important roles in the sorption of this group of herbicides.
Nowadays, the most common method for the removal of contaminants from aquatic environments is adsorption on polymeric materials [32,33,34]. Our previous study revealed that poly(divinylbenzene) modified in the Diels–Alder reaction with maleic anhydride and subsequent hydrolysis can be successfully used for removal of terbuthylazine from water solutions [35,36], also from multi-component systems [37]. Carboxyl groups, implemented into the adsorbent surface, create centers of specific interactions such as hydrogen bonding between the adsorbate and the adsorbent, which resulted in intensification and selectivity of adsorption. The next study showed that the proposed microspheres was equally effective in real solutions, containing significant amounts of other substances typical for agricultural wastewater, which is important from the point of view of its practical application, for example, in systems of agricultural wastewater treatment [38].
Since desethyl-terbuthylazine and 2-hydroxy-terbuthylazine are derivatives of terbuthylazine and there is a visible structural similarity, using poly(divinylbenzene) modified with maleic anhydride for adsorption of these compounds therefore seems to be a suitable option. In this work, the batch and column experiments were carried out to examine the efficiency of desethyl-terbuthylazine and 2-hydroxy-terbuthylazine adsorption on specific polymeric microspheres and the impact of different factors such as pH and temperature on this process. Moreover, the aim of this study is to examine the kinetics and mechanism of terbuthylazine metabolites adsorption in the proposed adsorbate–adsorbent systems.

2. Materials and Methods

2.1. Materials

Analyzed herbicides: desethyl-terbuthylazine and 2-hydroxy-terbuthylazine (purities ≥ 99%) were obtained from Dr. Ehrenstorfer GmbH (Augsburg, Germany). Millipore deionized water was used for making the solutions during sorption experiments. Ethanol (96%) and sodium hydroxide were purchased from P.P.H. “STANLAB” Sp. J. (Lublin, Poland).
Polymeric microspheres were made in our laboratory by modification of the pendant vinyl groups of poly(divinylbenzene) by maleic acid anhydride through the Diels–Alder reaction. Polymer synthesis and modification were carried out according to the previously published procedures [39,40].

2.2. Sorption Studies

The physicochemical properties of analyzed terbuthylazine metabolites are presented in Table 1.

2.2.1. Batch Studies

Different amounts of wet-weighed polymer adsorbent were contacted with 20 mL of 10 mg L−1 ethanol/water solutions (1/9, v/v) of terbuthylazine metabolites at pH 3, 5 and 7 by continuous shaking at 20, 35 and 50 °C for 48 h. The concentrations of terbuthylazine metabolites were measured using ultraviolet–visible (UV/VIS) spectroscopy, Jasco V-630 apparatus (JASCO EUROPE s.r.l., Cremella, Italy) The wavelength was set at 215.0 for desethyl-terbuthylazine and 238.0 nm for 2-hydroxy-terbuthylazine. The value of terbuthylazine metabolite sorption capacity (qeq), expressed in milligrams of adsorbate per gram of dry adsorbent, was calculated from the mass balance and the initial concentration of metabolite in the solution. The distribution coefficients (K) were calculated as the ratio of the amount of metabolite adsorbed by 1 g of a dry polymer and the amount of metabolite remaining in 1 mL of solution at equilibrium after sorption. Sorption kinetics was examined at 20 °C by mixing adsorbent with herbicide metabolite solutions (C = 10 mg L−1, pH = 7), collecting samples at fixed times and analyzing the residual concentration of metabolite.

2.2.2. Fixed Bed Sorption Studies

Dynamic sorption was carried out in fixed bed columns (diameter 0.9 cm) filled with 0.3 g of tested sorbent. Sorption was carried out by passing downwards a terbuthylazine metabolites solutions (C = 10 mg L−1, pH = 7) with a flow rate 0.9 mL min−1. Samples were collected at regular time intervals (15 min) using a fraction collector and the outlet concentration was measured by UV/VIS spectroscopy. For desorption, a similar procedure was used using 96% ethanol. The sampling time was decreased to 5 min and the process was carried out until full elution of the metabolites. Five consecutive sorption/desorption cycles for desethyl-terbuthylazine and three sorption/desorption cycles for 2-hydroxy-terbuthylazine were used for determining the reusability of microspheres. All of the experiments were carried out at room temperature (20 °C).
In order to analyze experimental data, the terbuthylazine metabolites’ adsorption breakthrough curves were determined by plotting the effluent metabolite concentration (C) versus the effluent volume (V). The calculations of total and usable adsorptive capacity, mass transfer zone and mass exchange zone moving rate were made based on the methodology described in detail in reference [46].

3. Results and Discussion

3.1. Adsorbent Characterization

The proposed microspheres for terbuthylazine metabolites sorption were designed in such a way that they have sets of donor and acceptor atoms complementary to the structures of triazines molecules. For this purpose, synthesized by radical suspension polymerization poly(divinylbenzene) microspheres were modified in the Diels–Alder reaction (see Figure 1). Carboxyl groups, which are capable of creating specific interactions, were obtained after opening the ring of maleic anhydride with the use of base hydrolysis and cyclization with hydrochloric acid.
Selected adsorbent has a predominantly mesoporous structure with an average pore size 6.62 nm and a total pore volume 0.96. Its scanning electron microscope (SEM) micrographs (ZEISS, Zeiss EVO LS 15, Oberkochen, Germany) are shown in Figure 2. The studied polymeric microspheres are characterized by a well-developed surface area (578 m2 g−1) and good water regain (2.54 g g−1), which have a direct impact on the adsorption efficiency.
The carboxylic group content in the structure of the proposed adsorbent results from the effectiveness of modifying free vinyl groups of poly(divinylbenzene). Whereas the number of vinyl groups that can be modified depends on the polymerization conditions. Many attempts have been made to synthesize poly(divinylbenzene) microspheres in order to obtain as many sites as possible for reacting with maleic anhydride. Samples characterized by a lower amount of pendant vinyl groups, that could be modified, have a higher degree of crosslinking and at the same time lower carboxylic groups content. On the other hand, the degree of crosslinking is responsible for the appropriate structure and strength of the polymer microspheres. The SEM micrographs of the tested adsorbent (Figure 2) and the characteristics of its porous structure indicate that the selection of polymerization conditions was correct. The content of carboxyl groups in the tested polymer is 3.58 mmol g−1. It should be emphasized that the commercially available adsorbents have usually about 1 mmol g−1 of carboxylic groups (e.g., 0.74 mmol g−1 in Strata-X-CW 33u Polymeric Weak Cation from Phenomenex (Torrance, CA, USA) and so the content of carboxyl groups in the proposed microspheres is highly satisfactory. These groups can enlarge the capacity of sorption by the formation of specific interactions such as hydrogen bonds between the adsorbent and adsorbate. These interactions can be created between the hydroxyl hydrogen atom of the carboxylic group from modified polymer and the nitrogen atom from triazine, which contains free electron pair. Such interactions can be created also between the hydrogen atom of the triazine amino group and the carbonyl oxygen atom from the microspheres. The ability to the formation of these direct interactions indicates better adsorption and selectivity than that of traditional adsorbents.

3.2. Batch Sorption Experiments

Sorption experiments were conducted for two metabolites of terbuthylazine: desethyl-terbuthylazine and 2-hydroxy-terbuthylazine. Modified with maleic anhydride, poly(divinylbenzene) was selected for sorption of these triazines due to structural similarity of derivatives to terbuthylazine, the sorption of which on this material was very effective [35]. Molecular structures of examined herbicides show complementarity to the configuration of functional groups of modified microspheres. Removal of triazines may arise from the formation of not only non-specific interactions but also from creating specific ones.

3.2.1. Adsorption Kinetics

The time needed to reach the equilibrium in the metabolite adsorption process was determined in the batch adsorption experiment. The plots of the equilibrium sorption (qeq) versus time (t) are presented in Figure 3. It can be seen that sorption kinetics for the tested compounds are very different. The most intense adsorption shows 2-hydroxy-terbuthylazine for which the process efficiency is reached 56% after 5 min. For desethyl-terbuthylazine, a much slower process is observed. The sorption up to 56% is achieved after 10 h. Sorption equilibrium for both metabolites is reached after about 30 h. These kinetic results also confirm a much better affinity of the tested adsorbent for the hydroxyl metabolite than for terbuthylazine, the sorption kinetics of which was presented earlier [37].
The description of adsorption kinetics was carried out using two models: pseudo-first and pseudo-second order [48]. Equation (1) was used to determine the parameters of the pseudo-first order model:
l o g ( q e q q t ) = l o g ( q e q ) k 1 t / 2.303
where t is time (min), q e q is sorption capacity at equilibrium (mg g−1),   q t is sorption capacity in time (mg g−1), k 1 is pseudo-first order adsorption rate constant (min−1). The parameters of the pseudo-second order model were calculated using Equation (2):
t q t = 1 k 2 · q e q 2 + t / q e q
where t is time (min), q e q is sorption capacity in equilibrium (mg g−1),   q t is sorption capacity in time (mg g−1), k 2 is pseudo-second order adsorption rate constant (g mg−1 min−1). The results are presented in Table 2.
Experimental data better correspond to the pseudo-second order model of kinetics. It confirms that adsorption of herbicides is a result of the specific interactions between triazines and surface of the adsorbent. The value of qeq is similar to the one obtained during sorption experiments (see Table 3). The calculated correlation coefficients R2 for this model are higher. In the case of the pseudo-first order model, qeq results do not reflect experimental values. The pseudo-first order equation is greatly influenced by random errors resulting from the experiment. In the case of a pseudo-second order equation, the impact of these errors is inconsiderable.

3.2.2. Adsorption Isotherms

Based on the batch studies the dependency of metabolites sorption on equilibrium concentration was plotted and shown in Figure 4. Analysis of the isotherms enabled the determination of the maximum sorption (qeq) and retention rate (R) for both terbuthylazine metabolites. It also allows determining the distribution coefficients (K) of triazines between the adsorbent and the adsorbate solution. All data were calculated based on results obtained for 10 mg L−1 metabolites solutions at 20 °C. Results are shown in Table 3.
The best sorption is shown by 2-hydroxy-terbuthylazine (113 mg g−1). It is much higher than the sorption of terbuthylazine (64 mg g−1) [37], which is the parent compound. Adsorption of organic substances may depend on the physicochemical properties of the adsorbate such as molecular weight, size and shape of a molecule, kind of functional groups, and solubility. For the compounds in question, these properties are similar, however the presence of OH group in the structure of 2-hydroxy-terbuthylazine may generate additional hydrogen bonds and therefore increase sorption efficiency. Referring to desethyl-terbuthylazine sorption, it was observed [35] that with the reduction of the presence of methyl groups in the structure of triazine-based herbicide, the capacity of sorption decreases, and due to that sorption of terbuthylazine is higher than desethyl-terbuthylazine (21 mg g−1). Five times greater sorption of 2-hydroxy-terbuthylazine than desethyl-terbuthylazine may also be the result of a greater contribution of non-specific sorption resulting from much lower solubility of 2-hydroxy-terbuthylazine (7.19 mg L−1) in water than desethyl -terbuthylazine (327.1 mg L−1). The calculated retention rates for both metabolites are similar, but the distribution coefficient for 2-hydroxy-terbuthylazine is significantly higher.
Figure 4 also presents the sorption isotherms obtained at different temperatures (20, 35 and 50 °C). Adsorption is an exothermic process, which makes it less efficient as the temperature increases. Additionally, adsorption is affected not only by adsorbate-adsorbent interactions but also by interactions between the dissolved substance and the solvent, which are significantly influenced by temperature. As temperature increases, the solubility of the adsorbate increases, thus reducing its non-specific interactions with the adsorbent. It was observed for both terbuthylazine metabolites that with the increasing temperature of the process, the sorption efficiency decreased. This may also have been caused by the weakening of hydrogen bonds formed between groups in the structure of adsorbate and groups on the polymer surface. It can be seen that for 2-hydroxy-terbuthylazine the difference in adsorption at various temperatures is much bigger. It is likely that in the case of a metabolite having an additional group capable of forming hydrogen bonds (–OH), the specific interactions have a much greater influence on the effectiveness of adsorption process than non-specific interactions between adsorbent and adsorbate.
The role of hydrogen bonds in the sorption of terbuthylazine metabolites was also confirmed by a study in which the hydrogens of the adsorbent’s carboxyl groups were eliminated. For this purpose, 0.1 M sodium hydroxide was used, and after ion-exchange the polymer was washed using distilled water to neutral pH. The results of metabolites sorption on the microspheres having functional groups in the form of sodium salt and for comparison in the acidic form are presented in Figure 5.
It can be seen that reducing the possibility of hydrogen bonding between the hydroxyl hydrogen atom of the carboxyl group from modified polymer and the nitrogen atom containing the free electron pair from adsorbate lowers the sorption efficiency by about 50%. The same effect was observed for the sorption of terbuthylazine [37].
The pH of the adsorbate solution has a huge effect on the adsorption process due to the possibility of affecting the interactions between adsorbate and adsorbent. Figure 6 presents adsorption isotherms for terbuthylazine metabolites obtained at different pH.
In the case of desethyl-terbuthylazine, we do not observe any significant influence of pH on the sorption efficiency. Chlorinated triazines like desethyl-terbuthylazine are mainly weak bases and their logarithmic acid dissociation constant pKa is approximately around 2 [43], and therefore in the range of studied pH (3–7) desethyl-terbuthylazine occurs in the form of inert molecules. The carboxylic groups on the surface of the adsorbent are mostly not dissociated (the pKa value of the investigated polymer is around 6), so the possibilities for the formation of hydrogen interactions remain unchanged. In contrast, in the acidic pH range, 2-hydroxy-terbuthylazine displayed a markedly different behaviour than desethyl metabolite. The decrease in sorption efficiency at pH 5 and 3 is attributed to the high basicity of hydroxy triazines. With pKa values around 5 [43], these compounds are protonated at lower pH, resulting in cationic species that can interact more easily with water molecules, increasing the hydration shell and at the same time weakening the interaction with the adsorbent.

3.2.3. Sorption Models

Four models such as Freundlich, Langmuir, Temkin and Dubinin–Radushkevich were used to match experimental data. Table 4 presents parameters of adsorption modeling using these calculations. The Langmuir model assumes monolayer adsorption on the polymer surface, which contains many active sites with equal energies that can only adsorb one adsorbate molecule each. This adsorption is often called ideal because the surface-active centers in the polymer are evenly distributed. This model is expressed as Equation (3):
q e q = q m · K α · C e q / ( 1 + K α · C e q )
where: q e q is the amount of solute adsorbed per mass of adsorbent, C e q is the equilibrium concentration of solute in the solution. K α i q m are characteristic of Langmuir equation (shown in Table 4) and were calculated from linearized form of Equation (3). The high R2 coefficient for desethyl-terbuthylazine makes it possible to conclude that the Langmuir model is well suited to describe the adsorption process of this triazine. For 2-hydroxy-terbuthylazine this factor is slightly lower. However, the determined q e q values are far too high compared to the experimentally obtained values.
The Freundlich model is used to describe non-ideal, reversible adsorption and is represented by the Equation (4):
q = K F · C e q 1 / n
in which: q is sorption capacity at monolayer saturation (mg g−1), C e q is the equilibrium concentration of metabolite in aqueous solution, (mg L−1) K F is Freundlich’s constant and 1/n is the function of strength of the adsorption. K F determines the sorption capacity of the system and exponent 1/n indicates the variation of free enthalpy associated with adsorption from the solution. When 1/n is 1, the division between phases occurs independently of the herbicide concentration and the free enthalpy of the process is constant. A value of 1/n less than 1 suggests adsorption, in which the active centers have less and less free enthalpy. 1/n greater than 1 suggests cooperative adsorption, which assumes that increasing the number of adsorbate molecules on the surface of the adsorbent results in increased free enthalpy [49]. Characteristics were calculated using linearized form of Equation (4) and are presented in Table 4. The Freundlich model may also be suitable for describing the process of removing triazines from aqueous solutions due to the high correlation coefficients. For both studied metabolites, the 1/n factor is greater than 1, which suggests that cooperative sorption may occur in the system, in which the adsorbate absorbed earlier affects the adsorption of the remaining particles in the solution [49]. This phenomenon is very common in adsorption processes.
The Temkin isotherm model includes interactions between adsorbent and adsorbate. It assumes that there is a linear decrease in the heat of adsorption of molecules in the layer, which is a result of an increase in the surface coverage. This isotherm model usually corresponds to a continuous, infinite (unlimited by minimum or maximum) energy distribution of adsorption sites. It also assumes that the process of adsorption is characterized by a uniform distribution of binding energies, up to maximum value [50]. The linear form of this model can be expressed as Equation (5):
q e = B T l n ( A T ) + B T l n ( C e )
where: B T = R T / b T is a constant related to the heat of adsorption, and b T (J mol−1) is Temkin isotherm constant, q e (mg g−1) is an amount of adsorbate in adsorbent at equilibrium, A T (L g−1) is Temkin isotherm equilibrium binding constant, T is the temperature (293 K) and R is the universal gas constant (8.314 J mol−1 K−1). The calculated results are presented in Table 4. Correlation coefficient for desethyl-terbuthylazine is the highest in this adsorption model, which may mean that this model describes best experimental data. The correlation coefficient for 2-hydroxy-terbuthylazine is definitely lower, which may mean that this model is not suitable for this system. For both studied metabolites, the value of the b T is positive, which confirms exothermic character of adsorption process.
The Dubinin–Radushkevich isotherm model is mainly applied to the mechanism of adsorption with Gaussian energy distribution onto heterogeneous surfaces. In this model, adsorption follows a pore filling mechanism. It assumes a multilayer character of the surface and involves van der Waals forces. It is also based on the value of the energy and provides information about the adsorption, whether it is chemical or physical [51]. When the parameter of mean free energy E is below 8 kJ mol−1 the adsorption process is considered as physical, whereas when parameter E is between 8 kJ mol−1 and 16 kJ mol−1 the process of adsorption is chemical [52]. Dubinin–Radushkevich isotherms were calculated using Equations (6)–(8):
l n q e = l n q s ( K a d ε 2 )
ε = R T l n ( 1 + 1 C e )
E = 1 / 2 K a d
where: q e (mg g−1) is an amount of adsorbate in adsorbent at equilibrium, q s (mg g−1) is theoretical isotherm saturation capacity,   C e (mg L−1) is adsorbate concentration at equilibrium, E (J mol−1) is mean free energy per molecule of adsorbate, K a d (mol2 J−2) and ε are Dubinin–Radushkevich constants. The calculated results are presented in Table 4. Values of mean free energy obtained for metabolites are lower than 8 kJ mol−1, which suggests that the process of adsorption is physical in both systems. Value of q e corresponds very well with experimental data for desethyl-terbuthylazine and correlation coefficient is high. For the other metabolite, 2-hydroxy-terbuthylazine, correlation coefficient is lower, however, the value of the q e is similar to the value obtained experimentally. Figure 7 shows experimentally obtained isotherms (pH = 7, T = 20 °C) for both metabolites with predicted model curves for every theoretical model.
The SSE parameter is known as the sum of the square estimate of errors and was obtained using software Origin Pro [53]. The obtained results are included in Table 4. The lowest values were obtained for Freundlich and Langmuir models for both studied systems. However, when deciding which model is the most suitable, the correlation coefficient should also be taken into consideration.

3.3. Fixed-Bed Column Studies

3.3.1. Sorption Studies

Sorption experiments were carried out in dynamic mode for both metabolites: desethyl-terbuthylazine and 2-hydroxy-terbuthylazine. The characteristics of the fixed bed are described by the breakthrough curve, which is the plot of adsorbate concentration in the effluent (C) versus effluent volume (V). Breakthrough curves plotted from the results obtained are presented in Figure 8.
The analysis of the breakthrough curves clearly indicates that the tested specific adsorbent prefers sorption of 2-hydroxy-terbuthylazine, for which an extended breakthrough curve can be observed, indicating that a higher volume of metabolite solution could be treated. For desethyl-terbuthylazine the polymeric bed is rapidly saturated. Calculating the parameters of sorption on the basis of the breakthrough curve is an important issue because it provides the basic information for the design of the column adsorption system and can be used to determine the rational scale of the column adsorption for practical application. Breakthrough curves obtained for the tested metabolites can be described by mathematical equations presented in Table 5. Correlation coefficients are high and above 0.99 for both terbuthylazine derivatives. Characteristic parameters for dynamic sorption were calculated following the methodology described in the publication [46] and presented in Table 6.
Analysis of obtained data allows to determine not only the total (Pt) and useful (Pu) adsorption capacity, but also bed operation time until breakthrough (tb) and exhaustion (te). The capacities of adsorption obtained under column studies are comparable to those obtained in batch adsorption systems (see Table 3). There is visible difference in time necessary for breakthrough and exhaustion of the column bed for used herbicide metabolites. The longer the time needed to reach those points the better filling of the polymer bed is observed. Time needed until breakthrough and until exhaustion are longer for 2-hydroxy-terbuthylazine, which may mean that the used microspheres have higher affinity for this metabolite. In fixed-bed experiments, the adsorbate concentration in the mobile phase as well as in the solid phase varies with time and position in the polymer bed. The mass transfer zone in the column moves from the top and towards to the bottom of the bed. In effective adsorption systems, low heights (H0) of mass transfer fronts and low mass exchange moving rate (u) are observed. Small values of this parameters mean that the breakthrough curve is close to an ideal step with negligible mass-transfer resistance [46].
The best adsorption parameters in fixed bed column studies were obtained for 2-hydroxy-terbuthylazine, due to the lower mass exchange moving rate and smaller height of adsorption front. The highest mass exchange moving rate for desethyl-terbuthylazine results probably from the weak affinity of this metabolite to the investigated polymer mainly because of its high water affinity (high water solubility—see Table 1). Therefore, also, the breakthrough and exhaustion times have the smallest values and consequently, the saturation of the adsorbent occurs earlier.

3.3.2. Sorption Models

There are many fixed bed models, which can be used to describe the results obtained such as the Bohart–Adams model, Thomas model, Yoon–Nelson model or Wolborska model. Table 7 presents parameters of column adsorption modeling using these calculations.
The Bohart–Adams model can be expressed by Equation (9):
l n ( C t C 0 ) = k A B · C 0 · t k A B · N 0 · ( Z U 0 )
where: C 0 and C t (mg L−1) are the initial and breakthrough concentration, k A B (L mg−1 min−1) is the Bohart–Adams kinetic constant, N 0 (mg L−1) is the saturation concentration, Z (cm) is the bed depth of the column, U 0 is the ratio of the volumetric flow rate (cm3 min−1) to the cross-sectional area of the bed (cm2). This model describes the relationship Ct/C0 in a continuous system. It is usually used to express the initial part of the breakthrough curve. This model is assumed to follow a “step isotherm” in which the adsorption rate is proportional to the fractional adsorption capacity possible on the surface of the adsorbent. This concept is also consistent with the basic principle of surface reaction theory [54]. Values of k A B and N 0 were calculated from linearized plot of l n ( C t C 0 ) versus time (t). The value of the determined constant k A B is higher for desethyl-terbuthylazine, while the value of N 0 is higher for 2-hydroxy-terbuthylazine. Correlation coefficients R2 for both metabolites of terbuthylazine are low, which may mean that this model is not suitable for modeling of obtained data. It suggests that adsorption mechanism is complex and involves few rate limiting steps [55]. The Bohart–Adams model may be too simple to explain this adsorption mechanism.
The Thomas model describes the adsorption process indicating that external and internal diffusion is not the limiting step. Thomas model is consistent with Langmuir adsorption and desorption kinetics. It is commonly used to calculate column performance and to predict breakthrough curves. It assumes negligible axial dispersion. The model is based on second order kinetics, which is the main limitation of this model. It does not constrain sorption by chemical reaction and is controlled by mass transfer at the surface. The Thomas model can be represented by Equation (10):
l n ( C 0 C t 1 ) = k T h · q 0 · m / Q k T h · C 0 · t
where: C 0 and C t (mg L−1) are the initial and breakthrough concentration, k T h (mL min−1 mg−1) is Thomas constant, q 0 (mg g−1) is the adsorption capacity, t (min) is the total flow time, m (g) is the amount of adsorbent placed in the column, Q (cm3 min−1) is the flow rate. Values of k T h and q 0 were calculated from linearized plot of l n ( C 0 C t 1 ) versus time (t). Parameter q 0 does not correspond well with experimental data for both metabolites. The value of the Thomas kinetic constant k T h is greater for desethyl-terbuthylazine. Correlation coefficient for desethyl-terbuthylazine is high, however considering parameter q 0 this model is also not sufficient for obtained data. For 2-hydroxy-terbuthylazine correlation coefficient is definitely lower, which may mean this model is also not suitable.
The Yoon–Nelson model is applied to a range of effluent concentrations between the column breakthrough and saturation times. It assumes that the rate of decrease in the adsorption probability of each adsorbate molecule is proportional to the adsorption capacity of the adsorbate molecules and the breakthrough probability of the adsorbent column [56]. This model is considered simple, thus it does not require any detailed data about the adsorbent such as type of material or properties [57]. Using this model, it is possible to minimalize error from the use of the Thomas model, which may occur at lower or higher periods of time of the breakthrough. This model is most suitable for sorption from single component solutions [54,58]. The Yoon–Nelson model may be expressed by Equation (11):
l n ( C t C 0 C t ) = k Y n · t k Y n · τ
where:   k Y n (min−1) is the rate constant, τ (min) is the time needed for 50% adsorbate breakthrough. The linear plot of l n ( C t C 0 C t ) versus time (t) enables us to determine values of k Y n and τ . For desethyl-terbuthylazine the k Y n constant is higher than for 2-hydroxy-terbuthylazine. The time required to 50% breakthrough of the column ( τ ) was significantly longer for 2-hydroxy-terbuthylazine, for desethyl-terbuthylazine corresponds well with experimental data. The correlation coefficient is high only in the case of desethyl-terbuthylazine which may mean this model could be suitable for this metabolite.
The Wolborska model is mainly used to describe breakthrough curves at low concentrations of effluent [59]. Mass transfer may be represented by Equation (12):
l n ( C t C 0 ) = β · C 0 N 0 · t β · Z U 0
where: C 0 and C t (mg L−1) are the initial and breakthrough concentration, β (1 h−1) is the kinetic coefficient of the external mass transfer, N 0 (mg L−1) is the exchange capacity, Z (mm) is the height of the fixed bed, U 0 (mm h−1) is the superficial fluid velocity. Parameters β and N 0 were determined from linear plot: l n ( C t C 0 ) versus time (t). Determined parameters β and N 0 are higher for 2-hydroxy-terbuthylazine. Correlation coefficients are low for both metabolites, which may mean that kinetics systems were not dominated by external mass transfer [60]. Precise simulation of breakthrough curves may be critical for controlling the effectiveness of adsorbent. It may also be useful for determining the life of adsorption bed and regeneration time or even the adsorption capacities [61]. Figure 9 shows experimentally obtained curves (pH = 7, T = 20 °C) for both metabolites with predicted model curves for every theoretical fixed-bed column model.
The SSE parameters (defined as the sum of the square estimate of errors) were obtained using software Origin Pro [53] and are presented in Table 7. The values are high, which can suggest that analyzed models may not be suitable for such complex systems. It can be seen that for desethyl-terbuthylazine the Thomas and Yoon–Nelson models results may be suitable, whereas in the case of 2-hydroxy-terbuthylazine none of the examined models seem to be correct.
In summary of column adsorption modeling, the correlation coefficients are lowest for the Bohart–Adams and Wolborska models for both metabolites of terbuthylazine, whereas coefficients for the Thomas and Yoon–Nelson models are higher and comparable. Used models are considered simple and may not be effective in describing such complex systems. Accordingly, to obtain more detailed information, other models should be explored.

3.3.3. Desorption Studies

Desorption process was analyzed based on the mass balance of adsorbed herbicide and then desorbed from the column bed. Investigated herbicide metabolites show good solubility in ethanol, and due to that it was chosen as an eluent in the desorption process. Parameters such as desorption efficiency, amount of eluent and desorption time were determined and are shown in Table 8. Desorption efficiency is higher for desethyl-terbuthylazine and is associated with solubility of metabolites in ethanol. Desorption time was longer for 2-hydroxy-terbuthylazine (315 min) than for desethyl-terbuthylazine (185 min). It was noted that better solubility determines shorter time needed for elution of substances. During this study, it was also observed for both metabolites that after passing small volume of ethanol significant amounts of the herbicides were desorbed, but process was carried out until complete elution from the bed.
Figure 10 shows the desorption efficiency as a function of eluent volume passed through column bed. After passing 13.5 mL of eluent, about 50% of desethyl-terbuthylazine was desorbed, whereas for 2-hydroxy-terbuthylazine 27 mL were needed. 2-Hydroxy-terbuthylazine shows better affinity to used microspheres, which makes this metabolite more difficult to remove from the column bed than desethyl-terbuthylazine. However, we can assume that even the small amount of eluent is enough to regenerate the bed from both herbicides in industrial conditions, which makes the process economically advantageous.
Regeneration of polymer bed is important from environmental and also economical point of view. Due to the effective removal of both metabolites of terbuthylazine from the adsorbent it was decided to conduct sorption–desorption cycles. For desethyl-terbuthylazine five cycles were performed, whereas only three for 2-hydroxy-terbuthylazine due to the very long time needed for that process. Ethanol was used as eluent and after passing this solvent, the polymer bed was washed using distilled water. Analyzing the obtained results from the sorption–desorption cycles, it was noted that removal of both herbicides was still very effective. The results are shown in Figure 11.
In the case of desethyl-terbuthylazine, adsorption efficiencies in first three cycles are very high (up to 100% removal was achieved) and comparable. In the 4th and 5th cycles the adsorption efficiency is slightly reduced to about 90%. The efficiency of desorption is about 90% and is comparable in five cycles conducted. From sorption-desorption cycles of 2-hydroxy-terbuthylazine it can be noted that sorption efficiency is reduced by almost 20% in the second cycle. The third cycle showed adsorption efficiency reduced by 8% compared to the second cycle. Desorption efficiency is comparable in every cycle and is about 80%. Those results confirm the possibility of using poly(divinylbenzene) modified with maleic anhydride many times not only in the detection of herbicides, but also in water treatment.

4. Conclusions

The proposed polymeric microspheres are suitable for the sorption of selected terbuthylazine metabolites. We observed that the differences in the physicochemical properties of the compounds caused five times greater sorption of 2-hydroxy-terbuthylazine than desethyl-terbuthylazine. The kinetic data for both metabolites fit well to a pseudo-second order kinetic model, which suggests that the adsorption rate is dependent on the availability of adsorption sites more than on herbicide concentration. The equilibrium data in the batch study follow the Freundlich isotherm for 2-hydroxy-terbuthylazine, and for desethyl-terbuthylazine the Temkin and Dubinin–Radushkevich models are suitable. The Dubinin–Radushkevich isotherm model suggests that processes of adsorption are physical in both systems, while the Temkin model confirms their exothermic character. The adsorption of desethyl-terbuthylazine does not depend on the pH. For 2-hydroxy-terbuthylazine, which is even much better adsorbed than its parent compound, the intensity of adsorption obtained under acidic conditions was considerably smaller (comparable with sorption capacity for desethyl derivative); however, such results can still be considered satisfactory. The results of this study also show that hydrogen bonds have a significant impact on sorption efficiency. These specific interactions are also responsible for the selectivity of sorption. For this reason, the proposed microspheres should work effectively in environmental conditions. The adsorption capacities obtained under dynamic conditions were comparable with batch results. In column adsorption modeling using the Bohart–Adams, Wolborska, Thomas and Yoon–Nelson models, a good fit for the Yoon–Nelson model was observed for desethyl-terbuthylazine. For 2-hydroxy-terbuthylazine other models should be explored. Column studies have shown that adsorbed metabolites of terbuthylazine can be easily desorbed by ethanol, and due to that fact, the modified poly(divinylbenzene) can be used several times without a major change in adsorption capacity.
The satisfactory sorption efficiency both in static and dynamic conditions suggests that poly(divinylbenzene) modified with maleic anhydride may be successfully used in terbuthylazine metabolites detection and monitoring using solid-phase extraction, specifically as a filling of SPE columns. Probably, it could also be used as packing material for a chromatographic column dedicated to triazine herbicide measurements. Moreover, it can be also used successfully in removing these persistent and mobile compounds from aquatic environments.

Author Contributions

Conceptualization, S.R.; Methodology, S.R.; Formal analysis, S.R.; Investigation, W.B.; Writing—original draft preparation, S.R. and W.B.; Writing—review and editing, S.R.; Visualization, S.R. and W.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been supported by a statutory activity subsidy from the Polish Ministry of Science and Higher Education for the Faculty of Chemistry of Wrocław University of Science and Technology.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Hvězdová, M.; Kosubová, P.; Košíková, M.; Scherr, K.E.; Šimek, Z.; Brodský, L.; Šudoma, M.; Škulcová, L.; Sáňka, M.; Svobodová, M.; et al. Currently and recently used pesticides in Central European arable soils. Sci. Total Environ. 2018, 613–614, 361–370. [Google Scholar] [CrossRef]
  2. Khan, M.A.; Brown, C.D. Influence of commercial formulation on the sorption and leaching behaviour of propyzamide in soil. Sci. Total Environ. 2017, 578, 158–166. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Silva, T.S.; Souza, M.D.F.; Teófilo, T.M.D.S.; Dos Santos, M.S.; Porto, M.A.F.; Souza, C.M.M.; Dos Santos, J.B.; Silva, D.V. Use of neural networks to estimate the sorption and desorption coefficients of herbicides: A case study of diuron, hexazinone, and sulfometuron-methyl in Brazil. Chemosphere 2019, 236, 124333. [Google Scholar] [CrossRef] [PubMed]
  4. Arp, H.P.H.; Brown, T.N.; Berger, U.; Hale, S.E. Ranking REACH registered neutral, ionizable and ionic organic chemicals based on their aquatic persistency and mobility. Environ. Sci. Process. Impacts 2017, 19, 939–955. [Google Scholar] [CrossRef] [PubMed]
  5. Reemtsma, T.; Berger, U.; Arp, H.P.H.; Gallard, H.; Knepper, T.P.; Neumann, M.; Quintana, J.B.; De Voogt, P. Mind the Gap: Persistent and Mobile Organic Compounds—Water Contaminants That Slip Through. Environ. Sci. Technol. 2016, 50, 10308–10315. [Google Scholar] [CrossRef] [Green Version]
  6. Rüdel, H.; Körner, W.; Letzel, T.; Neumann, M.; Nödler, K.; Reemtsma, T. Persistent, mobile and toxic substances in the environment: A spotlight on current research and regulatory activities. Environ. Sci. Eur. 2020, 32, 1–11. [Google Scholar] [CrossRef]
  7. Loos, R.; Locoro, G.; Comero, S.; Contini, S.; Schwesig, D.; Werres, F.; Balsaa, P.; Gans, O.; Weiss, S.; Blaha, L.; et al. Pan-European survey on the occurrence of selected polar organic persistent pollutants in ground water. Water Res. 2010, 44, 4115–4126. [Google Scholar] [CrossRef] [PubMed]
  8. European Commission. Directive 2006/118/EC of the European Parliament and the Council of 12th of December 2006 on the protection of groundwater against pollution and deterioration. Off. J. Eur. Union 2006, L 372/19, 19–31. [Google Scholar]
  9. Mladinic, M.; Želježić, D.; Shaposhnikov, S.A.; Collins, A.R. The use of FISH-comet to detect c-Myc and TP 53 damage in extended-term lymphocyte cultures treated with terbuthylazine and carbofuran. Toxicol. Lett. 2012, 211, 62–69. [Google Scholar] [CrossRef]
  10. Lovaković, B.T.; Pizent, A.; Kašuba, V.; Kopjar, N.; Micek, V.; Mendaš, G.; Dvoršćak, M.; Mikolić, A.; Milić, M.; Žunec, S.; et al. Effects of sub-chronic exposure to terbuthylazine on DNA damage, oxidative stress and parent compound/metabolite levels in adult male rats. Food Chem. Toxicol. 2017, 108, 93–103. [Google Scholar] [CrossRef]
  11. Dos Santos, E.V.; Sáez, C.; Martínez-Huitle, C.A.; Cañizares, P.; Rodrigo, M.A. Combined soil washing and CDEO for the removal of atrazine from soils. J. Hazard. Mater. 2015, 300, 129–134. [Google Scholar] [CrossRef]
  12. Velisek, J.; Stara, A.; Zuskova, E.; Kouba, A. Effects of three triazine metabolites and their mixture at environmentally relevant concentrations on early life stages of marbled crayfish (Procambarus fallax f. virginalis). Chemosphere 2017, 175, 440–445. [Google Scholar] [CrossRef]
  13. Tasca, A.L.; Puccini, M.; Fletcher, A. Terbuthylazine and desethylterbuthylazine: Recent occurrence, mobility and removal techniques. Chemosphere 2018, 202, 94–104. [Google Scholar] [CrossRef] [Green Version]
  14. Gardo Gold Plant Protection Product Characterization. Available online: https://ctgb.blob.core.windows.net/documents/84656e937a8f90a5345e1eac8c26018c_20060801_13145_02.html (accessed on 4 March 2021).
  15. Čelić, M.; Jaén-Gil, A.; Briceño-Guevara, S.; Rodríguez-Mozaz, S.; Gros, M.; Petrović, M. Extended suspect screening to identify contaminants of emerging concern in riverine and coastal ecosystems and assessment of environmental risks. J. Hazard. Mater. 2021, 404, 124102. [Google Scholar] [CrossRef]
  16. Robles-Molina, J.; Gilbert-López, B.; García-Reyes, J.F.; Molina-Díaz, A. Monitoring of selected priority and emerging contaminants in the Guadalquivir River and other related surface waters in the province of Jaén, South East Spain. Sci. Total Environ. 2014, 479–480, 247–257. [Google Scholar] [CrossRef]
  17. Koutnik, D.; Stara, A.; Zusková, E.; Kouba, A.; Velisek, J. The chronic effects of terbuthylazine-2-hydroxy on early life stages of marbled crayfish (Procambarus fallax f. virginalis). Pestic. Biochem. Physiol. 2017, 136, 29–33. [Google Scholar] [CrossRef]
  18. Velisek, J.; Stara, A.; Koutnik, D.; Machova, J. Effect of Terbuthylazine-2-hydroxy at Environmental Concentrations on Early Life Stages of Common Carp (Cyprinus carpio L.). BioMed Res. Int. 2014, 2014, 621304. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Velisek, J.; Koutnik, D.; Zuskova, E.; Stara, A. Effects of the terbuthylazine metabolite terbuthylazine-desethyl on common carp embryos and larvae. Sci. Total Environ. 2016, 539, 214–220. [Google Scholar] [CrossRef]
  20. European Food Safety Authority. Conclusion on the peer review of the pesticide risk assessment of the active substance terbuthylazine. EFSA J. 2011, 9. [Google Scholar] [CrossRef] [Green Version]
  21. Kiefer, K.; Müller, A.; Singer, H.; Hollender, J. New relevant pesticide transformation products in groundwater detected using target and suspect screening for agricultural and urban micropollutants with LC-HRMS. Water Res. 2019, 165, 114972. [Google Scholar] [CrossRef] [Green Version]
  22. Abdourahime, H.; Anastassiadou, M.; Arena, M.; Auteri, D.; Barmaz, S.; Brancato, A.; Bura, L.; Cabrera, L.C.; Chaideftou, E.; Chiusolo, A.; et al. Updated peer review of the pesticide risk assessment for the active substance terbuthylazine in light of confirmatory data submitted. EFSA J. 2019, 17, e05817. [Google Scholar] [CrossRef]
  23. European Commission. Directive (EU) 2020/2184 of the European Parliament and of the Council of 16 December 2020 on the quality of water intended for human consumption (recast). Off. J. Eur. Union 2020, L 435/1, 1–62. [Google Scholar]
  24. Guzzella, L.; Pozzoni, F.; Giuliano, G. Herbicide contamination of surficial groundwater in Northern Italy. Environ. Pollut. 2006, 142, 344–353. [Google Scholar] [CrossRef] [PubMed]
  25. Hermosin, M.; Calderon, M.J.; Real, M.; Cornejo, J. Impact of herbicides used in olive groves on waters of the Guadalquivir river basin (southern Spain). Agric. Ecosyst. Environ. 2013, 164, 229–243. [Google Scholar] [CrossRef]
  26. Karlsson, A.S.; Lesch, M.; Weihermüller, L.; Thiele, B.; Disko, U.; Hofmann, D.; Vereecken, H.; Spielvogel, S. Pesticide contamination of the upper Elbe River and an adjacent floodplain area. J. Soils Sediments 2020, 20, 2067–2081. [Google Scholar] [CrossRef] [Green Version]
  27. Masiá, A.; Campo, J.; Navarro-Ortega, A.; Barceló, D.; Picó, Y. Pesticide monitoring in the basin of Llobregat River (Catalonia, Spain) and comparison with historical data. Sci. Total Environ. 2015, 503–504, 58–68. [Google Scholar] [CrossRef]
  28. Cotton, J.; Leroux, F.; Broudin, S.; Poirel, M.; Corman, B.; Junot, C.; Ducruix, C. Development and validation of a multiresidue method for the analysis of more than 500 pesticides and drugs in water based on on-line and liquid chromatography coupled to high resolution mass spectrometry. Water Res. 2016, 104, 20–27. [Google Scholar] [CrossRef]
  29. Fingler, S.; Mendaš, G.; Dvoršćak, M.; Stipičević, S.; Vasilić, Ž.; Drevenkar, V. Herbicide micropollutants in surface, ground and drinking waters within and near the area of Zagreb, Croatia. Environ. Sci. Pollut. Res. 2016, 24, 11017–11030. [Google Scholar] [CrossRef]
  30. Tasca, A.L.; Fletcher, A. State of the art of the environmental behaviour and removal techniques of the endocrine disruptor 3,4-dichloroaniline. J. Environ. Sci. Health Part A 2018, 53, 260–270. [Google Scholar] [CrossRef]
  31. Mesquita, D.P.; Quintelas, C.; Amaral, A.L.P.D.; Ferreira, E.C. Monitoring biological wastewater treatment processes: Recent advances in spectroscopy applications. Rev. Environ. Sci. Biotechnol. 2017, 16, 395–424. [Google Scholar] [CrossRef] [Green Version]
  32. Parlapiano, M.; Akyol, Ç.; Foglia, A.; Pisani, M.; Astolfi, P.; Eusebi, A.L.; Fatone, F. Selective removal of contaminants of emerging concern (CECs) from urban water cycle via Molecularly Imprinted Polymers (MIPs): Potential of upscaling and enabling reclaimed water reuse. J. Environ. Chem. Eng. 2021, 9, 105051. [Google Scholar] [CrossRef]
  33. Bagheri, A.R.; Aramesh, N.; Khan, A.A.; Gul, I.; Ghotekar, S.; Bilal, M. Molecularly imprinted polymers-based adsorption and photocatalytic approaches for mitigation of environmentally-hazardous pollutants—A review. J. Environ. Chem. Eng. 2021, 9, 104879. [Google Scholar] [CrossRef]
  34. El-Said, W.A.; El-Khouly, M.E.; Ali, M.H.; Rashad, R.T.; Elshehy, E.A.; Al-Bogami, A.S. Synthesis of mesoporous silica-polymer composite for the chloridazon pesticide removal from aqueous media. J. Environ. Chem. Eng. 2018, 6, 2214–2221. [Google Scholar] [CrossRef]
  35. Ronka, S. Removal of triazine-based herbicides on specific polymeric sorbent: Batch studies. Pure Appl. Chem. 2016, 88, 1167–1177. [Google Scholar] [CrossRef]
  36. Ronka, S. Removal of triazine-based herbicides on specific polymeric sorbent: Fixed bed column studies. Pure Appl. Chem. 2016, 88, 1179–1189. [Google Scholar] [CrossRef]
  37. Ronka, S.; Kujawska, M.; Juśkiewicz, H. Triazines removal by selective polymeric adsorbent. Pure Appl. Chem. 2014, 86, 1755–1769. [Google Scholar] [CrossRef]
  38. Ronka, S.; Kucharski, M. Application of novel polymeric, highly specific adsorbent for the removal of terbuthylazine from complex environmental samples. Int. J. Environ. Anal. Chem. 2020, 1–14. [Google Scholar] [CrossRef]
  39. Stranix, B.R.; Darling, G.D. Cycloaddition Functional Polymers from (Vinyl)polystyrene. U.S. Patent 6,534,611 B1, 18 March 2003. [Google Scholar]
  40. Kica, M.; Ronka, S. The Removal of Atrazine from Water using Specific Polymeric Adsorbent. Sep. Sci. Technol. 2014, 49, 1634–1642. [Google Scholar] [CrossRef]
  41. Schmitt, P.; Garrison, A.; Freitag, D.; Kettrup, A. Separation of s-triazine herbicides and their metabolites by capillary zone electrophoresis as a function of pH. J. Chromatogr. A 1996, 723, 169–177. [Google Scholar] [CrossRef]
  42. PPDB. Available online: https://sitem.herts.ac.uk/aeru/ppdb/en/Reports/1494.htm (accessed on 4 March 2021).
  43. EPA (United States Environmental Protection Agency). Available online: https://comptox.epa.gov/dashboard/dsstoxdb/results?search=DTXSID80184211#details (accessed on 4 March 2021).
  44. PPDB. Available online: https://sitem.herts.ac.uk/aeru/ppdb/en/Reports/1495.htm (accessed on 4 March 2021).
  45. LGC. Safety Data Sheet According to 1907/2006/EC. Article 31. Available online: https://assets.lgcstandards.com/sys-master%2Fpdfs%2Fh84%2Fh35%2F10141506240542%2FSDS_DRE-C17305000_ST-WB-MSDS-2922617-1-1-1.PDF (accessed on 4 March 2021).
  46. Kaleta, J. Removal of methylene blue from aqueous solution by adsorption. Arch. Environ. Prot. 2007, 33, 45–54. [Google Scholar]
  47. Stranix, B.R.; Darling, G.D. Functional Polymers from Vinylpolystyrene. Diels-Alder Reactions with Olefins. J. Org. Chem. 1997, 62, 9001–9004. [Google Scholar] [CrossRef]
  48. Simonin, J.-P. On the comparison of pseudo-first order and pseudo-second order rate laws in the modeling of adsorption kinetics. Chem. Eng. J. 2016, 300, 254–263. [Google Scholar] [CrossRef] [Green Version]
  49. Liu, S. Cooperative adsorption on solid surfaces. J. Colloid Interface Sci. 2015, 450, 224–238. [Google Scholar] [CrossRef] [PubMed]
  50. Kamga, F.T. Modeling adsorption mechanism of paraquat onto Ayous (Triplochiton scleroxylon) wood sawdust. Appl. Water Sci. 2018, 9, 1. [Google Scholar] [CrossRef] [Green Version]
  51. Itodo, A.U.; Itodo, H.U. Sorption energies estimation using Dubinin-Radushkevich and temkin adsorption isotherms. Life Sci. J. 2010, 7, 68–76. [Google Scholar]
  52. Chabani, M.; Amrane, A.; Bensmaili, A. Kinetic modelling of the adsorption of nitrates by ion exchange resin. Chem. Eng. J. 2006, 125, 111–117. [Google Scholar] [CrossRef]
  53. Software Origin Pro 2021, Version 9.8.0.200; OriginLab Corporation: Northampton, MA, USA, 2021.
  54. Singh, S.K.; Mehta, D.; Sehgal, D. Fixed Bed Column Study and Adsorption Modelling on the Adsorption of Malachite Green dye from wastewater using Acid Activated Sawdust. Int. J. Adv. Res. 2015, 3, 521–529. [Google Scholar]
  55. Majumdar, S.; Baishya, A.; Mahanta, D. Kinetic and Equilibrium Modeling of Anionic Dye Adsorption on Polyaniline Emeraldine Salt: Batch and Fixed Bed Column Studies. Fibers Polym. 2019, 20, 1226–1235. [Google Scholar] [CrossRef]
  56. Bharathi, K.S.; Ramesh, S.P.T. Fixed-bed column studies on biosorption of crystal violet from aqueous solution by Citrullus lanatus rind and Cyperus rotundus. Appl. Water Sci. 2013, 3, 673–687. [Google Scholar] [CrossRef] [Green Version]
  57. Homem, N.; Vieira, A.M.S.; Bergamasco, R.; Vieira, M.F. Low-cost biosorbent based on Moringa oleifera residues for herbicide atrazine removal in a fixed-bed column. Can. J. Chem. Eng. 2018, 96, 1468–1478. [Google Scholar] [CrossRef]
  58. Sharma, R.; Singh, B. Removal of Ni (II) ions from aqueous solutions using modified rice straw in a fixed bed column. Bioresour. Technol. 2013, 146, 519–524. [Google Scholar] [CrossRef] [PubMed]
  59. Djelloul, C.; Hamdaoui, O. Dynamic adsorption of methylene blue by melon peel in fixed-bed columns. Desalin. Water Treat. 2014, 56, 1–10. [Google Scholar] [CrossRef]
  60. Huang, M.-R.; Lu, H.-J.; Song, W.-D.; Li, X.-G. Dynamic Reversible Adsorption and Desorption of Lead Ions Through a Packed Column of Poly(m-phenylenediamine) Spheroids. Soft Mater. 2010, 8, 149–163. [Google Scholar] [CrossRef]
  61. Chittoo, B.S.; Sutherland, C. Column breakthrough studies for the removal and recovery of phosphate by lime-iron sludge: Modeling and optimization using artificial neural network and adaptive neuro-fuzzy inference system. Chin. J. Chem. Eng. 2020, 28, 1847–1859. [Google Scholar] [CrossRef]
Figure 1. Scheme of poly(divinylbenzene) modification using maleic anhydride [47].
Figure 1. Scheme of poly(divinylbenzene) modification using maleic anhydride [47].
Materials 14 02734 g001
Figure 2. Scanning electron microscope (SEM) images of modified poly(divinylbenzene) beads.
Figure 2. Scanning electron microscope (SEM) images of modified poly(divinylbenzene) beads.
Materials 14 02734 g002
Figure 3. Sorption kinetics for desethyl-terbuthylazine and 2-hydroxy-terbuthylazine on modified poly(divinylbenzene) microspheres.
Figure 3. Sorption kinetics for desethyl-terbuthylazine and 2-hydroxy-terbuthylazine on modified poly(divinylbenzene) microspheres.
Materials 14 02734 g003
Figure 4. Adsorption isotherms for desethyl-terbuthylazine and 2-hydroxy-terbuthylazine at different temperatures: 20 °C, 35 °C and 50 °C.
Figure 4. Adsorption isotherms for desethyl-terbuthylazine and 2-hydroxy-terbuthylazine at different temperatures: 20 °C, 35 °C and 50 °C.
Materials 14 02734 g004
Figure 5. Adsorption isotherms for desethyl-terbuthylazine and 2-hydroxy-terbuthylazine on microspheres having carboxyl groups in acidic (p(DVB)-COOH) and sodium (p(DVB)-COONa) forms.
Figure 5. Adsorption isotherms for desethyl-terbuthylazine and 2-hydroxy-terbuthylazine on microspheres having carboxyl groups in acidic (p(DVB)-COOH) and sodium (p(DVB)-COONa) forms.
Materials 14 02734 g005
Figure 6. Adsorption isotherms for desethyl-terbuthylazine and 2-hydroxy-terbuthylazineat different pH: 3, 5 and 7.
Figure 6. Adsorption isotherms for desethyl-terbuthylazine and 2-hydroxy-terbuthylazineat different pH: 3, 5 and 7.
Materials 14 02734 g006
Figure 7. Experimental isotherms with predicted batch sorption model curves for desethyl-terbuthylazine (a) and 2-hydroxy-terbuthylazine (b).
Figure 7. Experimental isotherms with predicted batch sorption model curves for desethyl-terbuthylazine (a) and 2-hydroxy-terbuthylazine (b).
Materials 14 02734 g007
Figure 8. Adsorption breakthrough curves obtained in studied systems.
Figure 8. Adsorption breakthrough curves obtained in studied systems.
Materials 14 02734 g008
Figure 9. Experimentally obtained breakthrough curves and fixed bed column modeling curves for desethyl-terbuthylazine (a) and 2-hydroxy-terbuthylazine (b).
Figure 9. Experimentally obtained breakthrough curves and fixed bed column modeling curves for desethyl-terbuthylazine (a) and 2-hydroxy-terbuthylazine (b).
Materials 14 02734 g009
Figure 10. Desorption of terbuthylazine metabolites using ethanol.
Figure 10. Desorption of terbuthylazine metabolites using ethanol.
Materials 14 02734 g010
Figure 11. Efficiency of sorption and desorption of studied terbuthylazine metabolites.
Figure 11. Efficiency of sorption and desorption of studied terbuthylazine metabolites.
Materials 14 02734 g011
Table 1. The physicochemical properties of analyzed terbuthylazine metabolites.
Table 1. The physicochemical properties of analyzed terbuthylazine metabolites.
HerbicideStructureMolecular WeightSolubility in Water
20 °C (mg L−1)
pKa * [41]logKow *
desethyl-terbuthylazine Materials 14 02734 i001201.66327.1 [42]~21.94 [43]
2-hydroxy-terbuthylazine Materials 14 02734 i002211.267.19 [44]~52.54 [45]
* pKa, acid dissociation constant; Kow, octanol/water partition coefficient.
Table 2. Parameters of pseudo-first order and pseudo-second order adsorption kinetic model.
Table 2. Parameters of pseudo-first order and pseudo-second order adsorption kinetic model.
Terbuthylazine MetabolitePseudo-First Kinetic ModelPseudo-Second Kinetic Model
R2qeq
(mg g−1)
k1
(min−1)
R2qeq
(mg g−1)
k2
(g mg−1 min−1)
desethyl-terbuthylazine0.96513.9 1.8 × 10 3 0.97516.4 1.2 × 10 4
2-hydroxy-terbuthylazine0.91042.4 1.2 × 10 3 0.991115 1.4 × 10 4
Table 3. Sorption capacity (qeq), retention rate (R) and distribution coefficient (K) for adsorbed terbuthylazine metabolites.
Table 3. Sorption capacity (qeq), retention rate (R) and distribution coefficient (K) for adsorbed terbuthylazine metabolites.
Terbuthylazine Metaboliteqeq (mg g−1)R (%)K (-)
desethyl-terbuthylazine21802800
2-hydroxy-terbuthylazine1138824,400
Table 4. Parameters of Langmuir, Freundlich, Temkin and Dubinin–Radushkevich isotherm models.
Table 4. Parameters of Langmuir, Freundlich, Temkin and Dubinin–Radushkevich isotherm models.
Isotherm ModelParameterMetabolite
desethyl-terbuthylazine2-hydroxy-terbuthylazine
Langmuirqm (mg g−1)263476
Kα (dm3 mg−1)0.010.05
R20.9690.916
SSE * 2 × 10 4 4 × 10 5
Freundlich1/n1.371.01
KF ((mg g−1)(L mg−1)1/n)1.9820.81
R20.9490.913
SSE * 0.0580.026
TemkinAT (L g−1)0.611.05
bT (J mol−1)20140
R20.9890.834
SSE * 4.031051.26
Dubinin-Radushkevichqs (mg g−1)21.694.9
Kad (mol2 J−2) 2 × 10 6 6 × 10 7
E (kJ mol−1)0.50.9
R20.9700.731
SSE * 0.120.23
* SSE, Sum of Square Estimate of Errors.
Table 5. Mathematical equations of the breakthrough curve for investigated terbuthylazine metabolites.
Table 5. Mathematical equations of the breakthrough curve for investigated terbuthylazine metabolites.
Terbuthylazine MetaboliteMathematical Equation of the Breakthrough CurveR2
desethyl-terbuthylaziney = 203.32x5 − 790.86x4 + 774.9x3 − 345.84x2 + 82.888x + 0.11310.996
2-hydroxy-terbuthylaziney = 0.0056x5 − 0.1269x4 + 1.0373x3 − 4.074x2 + 8.6167x − 1.181370.987
Table 6. Characteristic parameters of dynamic sorption obtained in the tested systems.
Table 6. Characteristic parameters of dynamic sorption obtained in the tested systems.
Terbuthylazine MetaboliteDynamic Adsorption
Capacity (mg g−1)
Bed OperationTime until
Breakthrough
tb (h)
Bed Operation Time until Exhausted
te (h)
Height of Adsorption
Front
H0 (cm)
Mass Exchange Moving Rate
u (cm min−1)
Pt (Total)Pu (Useful)
desethyl-terbuthylazine1010.255.753.13 9.5 × 10 3
2-hydroxy-terbuthylazine113398.5075.752.80 6.1 × 10 4
Table 7. Parameters of Adams–Bohart, Thomas, Yoon–Nelson and Wolborska models.
Table 7. Parameters of Adams–Bohart, Thomas, Yoon–Nelson and Wolborska models.
ModelParameterMetabolite
desethyl-terbuthylazine2-hydroxy-terbuthylazine
Adams–Bohart k A B (L mg−1 min−1) 2.60 × 10 4 4.00 × 10 5
N 0 (mg L−1)4254.545,294.4
R20.6350.385
SSE * 2.63146.21
Thomas k T h (mL mg−1 min−1) 8.40 × 10 4 6.00 × 10 5
q 0 (mg g−1) 1.13 × 10 4 2.73 × 10 5
R20.9150.606
SSE * 4.39185.40
Yoon–Nelson k Y n (min−1) 8.40 × 10 3 6.00 × 10 4
τ (min)103.62489.8
R20.9150.606
SSE * 4.93185.40
WolborskaΒ (1 h−1)0.01840.0332
N 0 (mg L−1)1.1815.24
R20.6350.385
SSE * 2.63146.45
* SSE, Sum of the Square Estimate of Errors.
Table 8. Characteristics of terbuthylazine metabolites desorption process.
Table 8. Characteristics of terbuthylazine metabolites desorption process.
Terbuthylazine MetaboliteDesorption Efficiency (%)Amount of Eluent (mL)Desorption Time (min)
desethyl-terbuthylazine95166185
2-hydroxy-terbuthylazine84284315
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ronka, S.; Bodylska, W. Sorption Properties of Specific Polymeric Microspheres towards Desethyl-Terbuthylazine and 2-Hydroxy-Terbuthylazine: Batch and Column Studies. Materials 2021, 14, 2734. https://doi.org/10.3390/ma14112734

AMA Style

Ronka S, Bodylska W. Sorption Properties of Specific Polymeric Microspheres towards Desethyl-Terbuthylazine and 2-Hydroxy-Terbuthylazine: Batch and Column Studies. Materials. 2021; 14(11):2734. https://doi.org/10.3390/ma14112734

Chicago/Turabian Style

Ronka, Sylwia, and Weronika Bodylska. 2021. "Sorption Properties of Specific Polymeric Microspheres towards Desethyl-Terbuthylazine and 2-Hydroxy-Terbuthylazine: Batch and Column Studies" Materials 14, no. 11: 2734. https://doi.org/10.3390/ma14112734

APA Style

Ronka, S., & Bodylska, W. (2021). Sorption Properties of Specific Polymeric Microspheres towards Desethyl-Terbuthylazine and 2-Hydroxy-Terbuthylazine: Batch and Column Studies. Materials, 14(11), 2734. https://doi.org/10.3390/ma14112734

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop