Next Article in Journal
Macrocyclic Ionic Liquids with Amino Acid Residues: Synthesis and Influence of Thiacalix[4]arene Conformation on Thermal Stability
Next Article in Special Issue
Development and Validation of the LC–MS/MS Method for Determination of 130 Natural and Synthetic Cannabinoids in Cannabis Oil
Previous Article in Journal
Optical Anisotropy of Porphyrin Nanocrystals Modified by the Electrochemical Dissolution
Previous Article in Special Issue
A Comparative Study between Screen-Printed and Solid-Contact Electrodes for the Stability-Indicating Determination of Bromazepam
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Organic Matter on the Sorption of Cefdinir, Memantine and Praziquantel on Different Soil and Sediment Samples

by
Dragana Mutavdžić Pavlović
1,*,
Kristina Tolić Čop
1,
Helena Prskalo
1 and
Mislav Runje
2
1
Department of Analytical Chemistry, Faculty of Chemical Engineering and Technology, University of Zagreb, 10000 Zagreb, Croatia
2
Pliva Croatia TAPI R&D, 10000 Zagreb, Croatia
*
Author to whom correspondence should be addressed.
Molecules 2022, 27(22), 8008; https://doi.org/10.3390/molecules27228008
Submission received: 24 October 2022 / Revised: 10 November 2022 / Accepted: 15 November 2022 / Published: 18 November 2022

Abstract

:
Pharmaceuticals are known for their great effects and applications in the treatment and suppression of various diseases in human and veterinary medicine. The development and modernization of science and technologies have led to a constant increase in the production and consumption of various classes of pharmaceuticals, so they pose a threat to the environment, which can be subjected to the sorption process on the solid phase. The efficiency of sorption is determined by various parameters, of which the physicochemical properties of the compound and the sorbent are very important. One of these parameters that determine pharmaceutical mobility in soil or sediment is the soil–water partition coefficient normalized to organic carbon (Koc), whose determination was the purpose of this study. The influence of organic matter, suspended in an aqueous solution of pharmaceutical (more precisely: cefdinir, memantine, and praziquantel), was studied for five different types of soil and sediment samples from Croatia. The linear, Freundlich, and Dubinin–Raduskevich sorption isotherms were used to determine specific constants such as the partition coefficient Kd, which directly describes the strength of sorbate and sorbent binding. The linear model proved to be the best with the highest correlation coefficients, R2 > 0.99. For all three pharmaceuticals, a positive correlation between sorption affinity described by Kd and Koc and the amount of organic matter was demonstrated.

Graphical Abstract

1. Introduction

In the last 30 years, pharmaceuticals have become of great interest to many scientists. In fact, they are suspected of adversely affecting human and animal health because, due to a significant increase in production and consumption, they now ultimately enter the environment in large quantities, directly or indirectly, through various pathways. The physico-chemical properties (good water solubility and low degradability) of some of these substances allow them to pass through natural filters and water treatment plants. They are present in the environment either in the unchanged form of the parent compound or in the form of the corresponding active metabolites. In this way, they affect drinking water supplies, groundwater quality and the life of terrestrial and aquatic organisms. Since analytical techniques have not been able to achieve lower sensitivity (at the nano/micro scale) for many years, there is still a lack of information on the fate of many contaminants such as pharmaceuticals in the environment [1,2,3].
Environmental contamination by pharmaceuticals depends on several parameters, including the type and amount of a particular pharmaceutical. The transfer of pharmaceuticals from the treated area to surface water is primarily influenced by their water solubility, degradability, photostability, acid dissociation constant (pKa), partition coefficient (Kd), partition coefficient of organic carbon in water (Koc), and partition coefficient of octanol in water (Kow). Weather conditions and soil/sediment composition have a great influence on environmental behavior [1,4,5].
Sorption is an important process because it can significantly influence and control the fate and effects of organic compounds in the environment [6]. To describe the sorption affinity of pharmaceuticals to a solid, the partition coefficient Kd (L/kg) is commonly used. It is determined experimentally from the slope of the linear isotherm, representing the ratio of solid (qe) and liquid (Ce) phases according to the following equation [7]:
q e = K d × C e
In addition to the linear isotherm, the amount of pharmaceutical sorbed onto the solid samples can also be determined by Freundlich (Equation (2)) and Dubinin–Raduschkevich (Equation (3)) isotherms according to the following equations [5,8]:
l o g   q e = l o g K F × 1 n
q e = q m e x p β ε 2
where qe and Ce represent the sorbed amount of the pharmaceutical (µg/g) and the final concentration at equilibrium (µg/mL); KF is the Freundlich sorption coefficient ((μg/g)/(μg/mL)1/n); and n is the dimensionless Freundlich exponent describing the deviation from linearity between the sorbed pharmaceutical and the residual concentration at equilibrium; qm is the maximum sorption capacity (µg/g); β is a constant related to the average sorption energy (mol2/kJ2), and ε is the Polany potential obtained from equation (Equation (4)):
ε = R T l n 1 + 1 C e
The Freundlich isotherm is used to test whether the sorption process occurs on the heterogeneous surface of the sorbent at active sites with different levels of adsorption energy [9]. Since the parameter n indicates the intensity of sorption, its values can be used to infer the sorption. For example, when n < 1, the isotherm takes a concave shape because the free energy of binding between the sorbed molecules and the sorbent is lower. In the case of n = 1, the isotherm is linear because the free energy is constant at all concentrations of the tested molecules, and in the case of n > 1, the isotherm is convex because the free energy is higher as more molecules are sorbed near the adsorbent [10]. By determining the parameters of the Dubinin–Raduschkevich isotherm, information about the reaction mechanism and the nature of sorption can be obtained. If the value of the free energy of sorption E < is 8 kJ/mol, it can be assumed that physisorption predominates, in contrast to chemical sorption where the values of the free energy of sorption are between 8–16 kJ/mol [11,12].
Since the Kd value depends on the organic carbon (or organic matter) content, it is better to use the partition coefficient Koc, normalized to organic carbon, to model the distribution in the environment under consideration and obtain information about the potential risk of pharmaceuticals [4,13,14]. The relationship between Kd and Koc coefficients was defined by the following equation [15]:
K o c = K d % O C × 100 ;
% O M % O C = 1.724
Organic matter in soils and sediments is a complex, heterogeneous material. In addition to the familiar humic fractions such as humic acid (HA), fulvic acid, and humin, it also consists of a considerable amount of particulate organic matter. This organic matter is usually not bound to other soil/sediment components such as minerals, but the interaction of the particulate organic matter that accumulates on the soil/sediment surface with contaminants is important [16]. The situation is somewhat different for sediments because they are in contact with water most of the time. The water in contact with the sediment washes out many substances, including the organic matter, especially from the part that is constantly in the water. In sediments exposed to climatic influences due to water withdrawal during certain dry periods, a higher content of organic matter is observed. This content is still lower than that in soils from the same geographical area, but its importance is the same. Humic substances in sediments and organic particles in the water phase play an important role in the flocculation and accumulation of hydrophobic compounds, thus affecting the fate and dispersion of contaminants in the environment [17,18]. In addition, suspended organic matter can allow the formation of new active sorption sites, while higher concentrations of dissolved humic acid can cause desorption of pharmaceuticals from soil and facilitate their transfer to water [19].
However, few studies have investigated the interaction between organic matter (OM) and pharmaceutical agents such as diclofenac [20], sulfamethazine [21], and tetracyclines [16,22] and polycyclic aromatic hydrocarbons (PAH)/chlorophenol [17]. Additionally, Guo and coworkers [23] studied the sorption of tylosin (TYL) and sulfamethazine (SMT), two ionizable antibiotics commonly used in agriculture as veterinary therapeutics and growth promoters, to solid matter HA. The high sorption of these compounds to soils and sediments is often associated with a high HA content. OM is usually represented in the literature by humic acids, as these are commercially available. The literature emphasizes that the sorbing efficiency is related to the hydrophobicity of the sorbate, and that the interaction of HA with hydrophobic organic matter is mainly due to π–π interactions, hydrophobic effects and hydrogen bonding [24].
In this paper, three pharmaceutical target compounds from different therapeutic classes were selected. Praziquantel is a low-toxicity anthelmintic pharmaceutical known for the treatment of Schistosoma worms infections and other novel tropical diseases [25,26]. Memantine, an N-methyl-D-aspartate (NMDA) receptor antagonist is used to alleviate symptoms associated with Alzheimer’s disease [27]. Cefdinir is a cephalosporin antibiotic used to treat urinary tract infections, osteomyelitis, and meningitis. Its presence in wastewater has been confirmed at high concentrations (125–175 mg/L), making it and other cephalosporins a group of emerging contaminants [28]. Although environmental data were found only for cefdinir to the authors’ knowledge, all three pharmaceuticals tested are widely distributed and are expected to be present in the environment. Ribeiro and coworkers [29] have drawn attention to this problem with respect to cephalosporins, but also to other widely used antibiotics [30], particularly with respect to their behavior in different soil types [31]. The fact that cephalosporins can be sorbed well to soil and sediment particles is shown by the sorption potential of the aqueous cephalexin solution to interact with activated charcoal and grapevine, where recovery rates of more than 80% have been achieved with this method of removal from an aqueous medium.
The aim of this paper was to experimentally determine how the influence of organic matter affects the sorption affinity of all three pharmaceuticals by using natural soil and river sediment samples. These samples possessed different properties depending on the sampling location from which they were taken in different Croatian areas. Linear, Dubinin–Radushkevich and Freundlich isotherms were used to describe the sorption tendency. The mechanism of sorption was defined by kinetic models: Lagergren’s pseudo-first-order pseudo-model, second-order pseudo-model and intraparticle diffusion model. The thermodynamic parameter standard Gibbs energy (ΔG°) was calculated from the estimated Kd values at 25 °C. The results of the present study may contribute to a better understanding of the mobility and fate of the studied pharmaceuticals in ecosystems.

2. Materials and Methods

2.1. Materials

The following high-purity (>99%) pharmaceuticals were used in this study: cefdinir (CEF), memantine (MEM), and praziquantel (PRAZ). CEF and MEM were purchased from Sigma-Aldrich (Steinheim, Germany), whereas PRAZ was purchased from Genera d.d. (Kalinovica, Croatia). The formulas of the studied pharmaceuticals and their physico-chemical properties are listed in Supplementary Materials Table S1. Stock solutions of each pharmaceutical were prepared separately by dissolving a certain amount of the pharmaceutical into an appropriate volume of 0.01 M CaCl2 to obtain working standard solutions (2.0; 1.0; 0.5; 0.3; 0.2 and 0.1 mg/L) by serial dilution. The working solutions of the studied pharmaceuticals 2.0 mg/L also contain a maximum of 1% methanol for better dilution of the studied pharmaceuticals in aqueous solutions. All solutions were stored in a way that protected them from light at 4 °C. Acetonitrile HPLC grade was supplied by Baker (Deventer, The Netherlands) and formic acid was purchased from Merck (Darmstadt, Germany).
Humic acid, manufactured by Sigma-Aldrich (Steinheim, Germany), was obtained as a dry formulation powder. A total of 1 mg/L humic acid solution was prepared by dissolving 0.5 mg humic acid powder in 500 mL of 0.01 M CaCl2. The solution was placed in an ultrasonic bath to ensure complete dissolution of humic acid. The pH of this humic acid solution was checked to establish that is was 7.00, and adjusted accordingly if there were minor deviations.

2.2. Sediment and Soil Samples

Experiments were performed with five samples of river sediment and five natural soil samples. The sediment and soil samples were collected in six different Croatian regions: Sisak-Moslavina County (sediment 1—Glina River, sediment 2—Lonja River, sediment 3—Petrinjčica River), Zadar County (soil 1—Bruvno and soil 2—Gračac), Brod-Posavina County (soil 3—Dolina and soil 4—Ljupina), the Požega-Slavonia County (sediment 4—Pakra River), Osijek-Baranja County (soil 5—Josipovac) and the Primorsko-Goranska County (sediment 5—Studena River). All samples were collected in an area far from human activities (no farms or factories), so the samples should not contain pharmaceuticals, especially those tested. Major or minor differences in their physicochemical properties give a clear picture of which of the physicochemical parameters is crucial for the sorption of a particular tested pharmaceutical.
Samples of both sediments were collected during the summer by hand-held device, and with a trowel below the aqueous layer when the water level was lower, which greatly simplified the overall sampling procedure. All collected solid samples were dried, ground, sieved through a 2 mm sieve, and characterized in the previously described manner [13].
Table S2 shows the physicochemical properties of the sorbent samples used.

2.3. Batch Sorption Experiments and Data Analysis

All batch sorption experiments of the tested pharmaceuticals on river sediment and natural soils were performed according to the OECD 106 procedure [32]. The procedure is performed in triplicate by shaking with a laboratory shaker (Innova 4080 Incubator Shaker, NewBrunswick Scientific, Edison, NJ, USA), which allows continuous contact between the sediment or soil samples and the solutions of the tested pharmaceuticals. To avoid photolytic degradation, shaking is performed in the dark, and to avoid microbiological activity, all samples are sterilized beforehand.
To determine the contact time required to reach sorption equilibrium, preliminary experiments were performed by shaking soil/sediment samples with solutions of the tested pharmaceuticals at different time intervals (10, 20, 30, 40 and 50 min and 1, 2, 4, 6, 12, 18, and 24 h) at 25 °C. Derived from the preliminary experiments, a time period of 24 h was chosen for all further experiments. These results also form the basis for gaining knowledge about sorption kinetics.
The same sorption procedure is performed for each of the tested pharmaceuticals. The procedure consists of adding 10 mL of the solution of one of the tested pharmaceuticals at a known concentration (0.1–2.0 mg/L) in 50 mL of laboratory glassware to 1 g of air-dried sediment or soil samples. Based on the preliminary tests, the prepared suspension was shaken in a shaker (at 200 rpm) for 24 h at controlled temperature conditions (adjusted at 25 °C), filtered through 0.45 µm syringe filters, and added to HPLC vials. To investigate the effect of organic matter on the sorption of studied pharmaceuticals, experiments were performed with four concentration levels of humic acid (0; 1; 10 and 100 mg/L). Therefore, 9 mL of the prepared humic acid solution were added to all the laboratory glassware prepared for the series of experiments with a humic acid concentration of 1 mg/L of solution, as previously mentioned (the total volume of the solution in the laboratory glassware was 10 mL).
When investigating the effects of higher humic acid concentrations (10 mg/L and 100 mg L) on the sorption of the tested pharmaceuticals, the procedure is somewhat different. Since humic acids are sparingly soluble, it is not possible to prepare solutions containing 10 and 100 mg/L. Therefore, the required amount of humic acids is added to 10 mL of the pharmaceutical solution (which is already in contact with the soil and sediment) to obtain the desired concentration of humic acids. For these experiments, 0.10 mg (for 10 mg/L) and 1.00 mg (for 100 mg/L) of humic acid were dissolved in 10 mL of one of the six concentrations of each pharmaceutical. All experiments were performed in 0.01 M CaCl2 solution at the initial pH of pharmaceutical under study (in the case of CEF this is pH 6.0, in the case of PRAZ pH 6.5, and in the case of MEM pH 7.0).
The residual concentration of the pharmaceutical (CEF, MEM and PRAZ) in the remaining liquid phase after sorption was analyzed by the UHPLC-MS (Agilent 6490 coupled with Agilent Infinity UHPLC system Triple Quadrupole Mass Spectrometer, Santa Clara, CA, USA) with electrospray ionization according to the method described previously [33].

3. Results and Discussion

3.1. Effect of Contact Time and Initial Concentration of Pharmaceuticals

To determine the sorption coefficient, it is first necessary to determine the time within which the maximum sorption of the tested compounds can be achieved for the soil and sediment samples tested. Since the initial concentration of the tested compounds plays an important role, three different concentrations (0.1, 0.5 and 2.0 mg/L) were prepared for each tested pharmaceutical. Each of these solutions were contacted with the tested soil and sediment samples and shaken for different time intervals (10, 20, 30, 40, 50 min and 1, 2, 4, 6, 18 and 24 h). Some of the results can be seen in Supplementary information, Figure S1. Sorption kinetics of CEF, PRAZ, and MEM can be divided into two phases: “fast” and “slow” sorption [34]. During the first hours (at the lowest concentration) or 6 h after shaking (at the higher concentrations), intense sorption of the pharmaceuticals takes place on the tested samples. After the mentioned times, the sorbed amounts do not change in all cases until the final equilibrium is reached due to the gradual saturation of the sorption active sites [34]. Although in all cases the maximum sorption is reached after 18 h of shaking, to simplify the experiments all further experiments were performed with 24 h of shaking. The lowest percentage of sorption for CEF and PRAZ solutions is obtained at the lowest concentration (0.1 mg/L), while in the case of MEM the lowest sorption is obtained at the highest concentration (2.0 mg/L). For CEF sorption, the percentage of sorption in both samples shown changes only slightly with the change in concentration, while for PRAZ sorption the sorption effect is significant or at least more pronounced compared to CEF sorption.

3.2. Sorption Isotherms

In the determination of sorption coefficients, concentrations of the studied pharmaceuticals in the range of 0.1 mg/L to 2 mg/L were used, corresponding to a range of 1 to 20 mg/kg soil/sediment. The range of concentrations studied is consistent with the average concentrations found in studies of a similar nature (about 0.2–20 mg/kg) [21,35,36], although the ranges in environmental samples are much larger (for sediment 0.02–285 mg/kg, for soil 0.034–530 mg/kg) [37,38]. CEF, MEM, and PRAZ are relatively unexplored in environmental samples, and it is difficult to say whether the concentration chosen for this study corresponds to their actual concentration in the environment. Nevertheless, the lowest concentrations were chosen to allow for analytical detection [39].
During the performance of all experiments, the stability of the tested pharmaceuticals in solution and their possible sorption on the walls of the test vessels were investigated. For this purpose, the control samples CEF, MEM and PRAZ were used with the highest concentrations (2 mg/L) in 0.01 M CaCl2 and analyzed for LC-MS/MS. The obtained results showed that all pharmaceutical products are stable in aqueous solutions for the required time and no decrease in their concentration was observed during the performance of all experiments, which means that they were not sorbed on the walls of the test vessels. Additionally, no interferences of the ten different sediment and soil sample matrices were observed in the chromatograms under the specified experimental conditions, which gives the method sufficient selectivity. The characteristic mass spectra of the blank sediment samples in 0.01 M CaCl2 and 2 mg/L CEF, MEM and PRAZ in 0.01 M CaCl2 are shown in Figure 1.
It can be seen from Figure 1 that sediment sample 1 (Glina) did not contain any of the investigated pharmaceuticals (CEF, MEM and PRAZ). The absence of the characteristic mass spectra of the mentioned pharmaceuticals (specific precursor ion/production transitions for each of the investigated pharmaceuticals) is sufficient evidence for this. Identical results were obtained for other sediment and soil samples in support of the above results, i.e., the fact that the soil and sediment samples used do not contain pharmaceuticals, in previously publications [13,40]
The linear sorption isotherms of all three pharmaceuticals tested are shown in Supplementary Information, Figure S2. Table 1, Table 2 and Table 3 show the sorption coefficients and all other parameters derived from three different models of sorption isotherms (linear, Freundlich, and Dubinin–Radushkevich) for each soil and sediment tested.
From the obtained regression coefficients R2, it can be concluded that only the linear isotherm describes the sorption process with a value of R2 > 0.989 in the case of all tested pharmaceuticals, while the regression coefficients for the Freundlich isotherm generally range between 0.738–0.993. It is also evident from the presented tables that the Dubinin–Radushkevich sorption model does not describe the situation very well for all three tested pharmaceuticals (R2 ranges from 0.482–0.999), but this was not achieved with any other applied model for all pharmaceuticals tested simultaneously. The Dubinin–Radushkevich model, based on the value of free sorption energy E (kJ/mol), was chosen only because, unlike all other models, it can inform the reader about the sorption mechanism operating between the studied pharmaceuticals and the soil/sediment samples [9]. The Dubinin–Radushkevich isotherm shows the worst agreement with the experimental data, and its regression coefficients are generally in the range of 0.671–0.999, although there are data whose R2 is below the indicated range in the case of MEM. On the basis of the calculated E values shown in Table 1, Table 2 and Table 3, sorption between all tested pharmaceuticals and sorbents was characterized as a physical process based on weak van der Waals forces.
Based on the values obtained for MEM (Kd ranging from 0.75–38.51 mL/g, KF ranging from 0.91–66.07 ((μg/g)(mL/μg)1/n), and the qm values referring to the amount of maximum adsorbed substance), it is clear that MEM binds less to the tested sediments and soils compared to other tested pharmaceuticals. Comparing the Kd values in general for all three tested pharmaceuticals, it is observed that the highest Kd values were obtained for CEF (Kd = 10.69–972.73 mL/g), although the highest Kd values were still obtained between soil 1 and PRAZ. PRAZ showed better mobility in soil and sediment samples with lower Kd values compared to CEF (Kd = 0.27–69.70 mL/g) and bound less to the sediment and soil samples tested than the other pharmaceuticals compared to CEF, although it bound better than MEM.
From the values presented in Table 1, Table 2 and Table 3, it can be seen that for CEF in the experiments with the addition of 100 mg/L HA, generally all values of n are less than 1 (0.16–0.64). This indicates lower sorption for the tested soils and sediments compared to other pharmaceutical products tested, and also indicates higher heterogeneity of the sorbate surface. The n-values obtained for MEM are higher than 1 for seven samples, indicating higher sorption intensity. PRAZ showed constant sorption affinity over the entire concentration range, with n-values close to 1 for almost all soil and sediment samples used. Exceptions were observed in the cases of sediment 3 with 0 mg/L HA and sediment 4 with 100 mg/L HA (n = 2.584), sediment 4 (n = 1.736), and sediment 3 (n = 2.070), supporting the fact that more and more PRAZ molecules are sorbed near the soil/sediment and that the previously sorbed molecules enhance further sorption by modifying the mentioned surfaces. For CEF, n > 1 is only in the case for soil 1 (n = 1.5198), soil 2 (n = 1.8560), sediment 4 (n = 1.370) and sediment 3 (n = 1.1350), while MEM reaches all values for the soil/sediment samples studied (n = 1.5868–3.2020). In general, howeverm the Kd values indicate a low binding affinity of MEM.

3.3. Influence of Organic Matter and Composition of Soil/Sediment Samples

Soil organic matter is a complex, heterogeneous material [16] called humus. Humus consists of humic and fulvic acids and humin. The organic component of soil or sediment is one of the most influential parameters affecting the sorption of pharmaceuticals, which is confirmed by numerous works from the literature [21,24,41,42,43]. The effect of soil organic components on pharmaceutical mobility or sorption is twofold. This notion is supported by available research [44], where increased sorption was observed at a lower concentration of OM (1 mg/L), in contrast to a higher concentration of humic acid (10 mg/L), which resulted in an opposite effect. On the other hand, studies by Ling et al. [45] show that the binding of OM to sediment particles increases the organic carbon content and changes the surface morphology, which directly leads to an increase in adsorption sites on the sediment surface, i.e., increases the adsorption of the tested component (oxytetracycline). This is because organic matter can form complexes with relatively polar compounds or absorb hydrophobic organic contaminants, apparently increasing water solubility and possibly decreasing sorption. In addition, the presence of organic matter can promote adsorption of pharmaceuticals, and the higher the concentration of OM is, the more pronounced this effect becomes [41].
Since the pharmaceuticals used in this paper are relatively unexplored in terms of their behavior in solid environmental samples (soil, sediment), the question arises as to how they behave in the presence of different amounts of organic matter. To determine whether and how OM affects the sorption of CEF, MEM, and PRAZ to the soil and sediment samples studied, experiments were conducted with the addition of humic acids. Just as the partition coefficient (Kd) increases with the amount of organic matter, the normalized distribution constant of organic carbon content (Koc) also increases with an increasing amount of organic carbon. The results are shown in Figure 2.
Sediment 5 and soil 2 in the case of MEM have the lowest partition coefficients at 0 and 100 mg/L humic acids, indicating that these samples have the lowest affinity for sorption. On the other hand, sediment 2 and soil 5 have the highest affinity among the studied soils and sediments, as their partition coefficients are highest at 100 mg/L humic acids, being 13.10 and 38.51 mL/g, respectively. The Koc values in the case of MEM sorption range from 60.81 to 1557.9 mL/g.
Considering the other tested pharmaceuticals, we find that Kd reaches the highest values for CEF sorption rather than for MEM and PRAZ. The value is highest for sediment 2 (720.4 mg/L) and soil 5 (972.7 mg/L) and in the case of MEM. From the physico-chemical properties of soil 5, there is a very high humus content of 11.59%, the highest of all sediments and soils studied, as well as the highest conductivity (144.30 µS/cm) and a TDS value of 92.40 mg/L. The physico-chemical properties reported for soil 5 suggest that the high value of the above factors has a great influence on the sorption capacity of soil 5. Sediment 4 and soil 2 have the lowest values of partition coefficients of 100.7 and 92.25 mg/L, respectively. In addition to soil 5 and sediment 2, soils 3 and 4 also have large sorption capacities. Since large amounts of partition coefficients were recorded, CEF resulted in large Koc values, ranging from 1213.0 to 85,655.4 mL/g.
An increase in humic acid concentration results in greater sorption to the solid phase and with respect to PRAZ. This is reflected in a series of increases in the partition coefficients for all five soils studied, as well as for the sediment samples. The ratio of the equilibrium concentration of sorbed substance and residue in solution is highest for soil 1 with 116.42 mL/g and for sediment 1 with 31.45 mL/g. On the other hand, the sorption of PRAZ under organic influence is weakest in soil 4 because the distribution between solid and liquid phases is 9.609 mL/g at a HA concentration of 100 mg/L. In the calculation of Koc, values between 15.3 and 5606.37 mL/g were obtained depending on the organic matter content.
It is very important to point out that increasing the concentration of HA also affects other physicochemical properties of the soil/sediment. For example, there is a non-negligible effect on the pH of soil and sediment samples, which can easily shift towards more acidic values with an increase in the concentration of HA [24]. It is important to note that this change is not the same for all sediment and soil samples studied, as some soils and sediments are resistant to pH changes due to their high buffering capacity, which is due to the presence of various organic acids and other soil/sediment components that can form neutral salts.
The study of the influence of organic matter on sorption is certainly of greater importance. Humic acids should not be considered as the only organic component that can influence sorption and as such can be found in the soil/sediment. It is convenient to observe because it is a natural component of soil and sediment and is very readily available for the experimental part since it can be extracted and then easily tested. Humic acids here can stand for any organic component/substance that occurs in soil/sediment and thus affects the behavior of pharmaceuticals in the environment. For these three tested pharmaceuticals (CEF, MEM and PRAZ), it is observed that as the concentration of organic matter (in the form of humic acids) increases, their sorption in soil and sediment increases, which significantly affects their mobility in the environment. In this respect, CEF is the most immobile in the environment, as soils/sediments rich in organic matter are less hazardous to water than to aquatic organisms. At the same time, soils/sediments rich in organic matter are not able to bind PRAZ, especially MEM, which poses a potential hazard in water.
Since sorption depends not only on the physicochemical properties of the test substances, but also on the physicochemical properties of the studied soils and sediments, a correlation was performed for all studied pharmaceuticals to determine how changes in humic acid concentration affect the physicochemical parameters of the soils/sediments.
Based on the Pearson coefficient [7] from Table 4, it can be seen that the different physicochemical parameters do not have the same influence on the sorption of all three tested pharmaceuticals. In the case of CEF, there is a weak correlation between Kd and the texture of the soil/sediment samples (coarse sand, silt, and fine sand content), while other physicochemical parameters indicate that there is no significant linear correlation, as the r values generally range from 0 to ±0.25. With increasing humic acid concentration, the association increases in the case of silt and CEC content of Zn, Cu and Mn, although all this is not significant. However, in practice, it is known that the significance of the correlation coefficient is evaluated at the significance level. So, based on p < 0.05, we conclude that Pearson’s linear correlation coefficient is statistically significant in several cases, i.e., we can say that there is no statistically significant correlation between Kd CEF values and coarse sand, OM, CaCO3, CEC content and Cu content, as p < 0.05. For other parameters, we also cannot claim this with statistical certainty, since p is slightly higher than 0.05, and so a larger number of measurements should be performed. However, in the case of memantine, the situation is somewhat different. Based on the Pearson coefficients presented in Table 4, a weak correlation is observed between Kd for MEM and coarse sand, fine sand, CEC and CaCO3, this correlation being negative in the case of coarse sand and CaCO3 content. A moderate coupling is observed for pH, while the coupling is excellent in the case of EC or total dissolved salts (TDS). An increase in humic acid concentration leads to a decrease in the Pearson coefficient in almost all cases (the only exceptions are clay, OM and the percentage of Mn, Cu and Fe). In the case of EC and TDS, where the initial value of the Pearson coefficient (γ (HA) = 0 mg/L) indicates an excellent association with Kd for MEM, an increase in humic acid concentration (γ (HA) = 100 mg/L) also resulted in a decrease in the association to a nonsignificant level. The evaluation of the significance level of the calculated correlation coefficients showed that there is no correlation between the studied variables (r < ±0.25) for clay and the proportions of Zn, Cu, Fe and Mn, as the correlation coefficient is significant (p < 0.03). On the basis of p < 0.05, it can be claimed that the correlation between Kd and OM increases significantly (on the basis of the value of Pearson’s coefficient, the transition from non-significant to weak correlation was observed), while the situation is reversed for coarse sand and CE. In PRAZ, a strong correlation (r = 0.50–0.75) was observed between Kd and the content of clay, which changed slightly with the increasing concentration of humic acids (p < 0.05). The content of Fe and Mn in soils/sediments significantly affected the Kd value of PRAZ (p < 0.05). For other physicochemical parameters, no correlation was found between Kd values with or without an increase in humic acid concentration.
From the correlation analysis performed for all three tested pharmaceuticals (CEF, MEM and PRAZ), it can be concluded that pH does not change significantly with increasing humic acid concentration. This is supported by the fact that all tested soils and sediments have a high buffering capacity and resist pH changes due to the formation of neutral salts between the metal ions and the humic acids or organic acids present, which is confirmed by the above data.
The experimentally studied dependence of the organic component of the soil/sediment on the sorption coefficient Kd of the tested pharmaceuticals and the performed correlation analysis confirm again that sorption is a complex process that depends not only on the physicochemical properties of the tested pollutants, but also on the physicochemical properties of the soil/sediment.

3.4. Sorption Kinetics

The kinetic study of the pharmaceutical sorption process was performed at three concentration levels to cover the entire range used in other experiments. The sorption results were interpreted using two kinetic models—Lagergren’s (pseudo-first order) and Ho´s (pseudo-second order) (Equations (7) and (8)) [46,47].
ln q e q t = l n q e k 1 t
t q t = 1 k 2 q e 2 + t q e
where qt and qe represent the amount of sorbed analyte at time t and at equilibrium (μg/g), k1 is the constant rate pseudo-first-order value (min−1), and k2 is the constant rate of kinetic model pseudo-second-order value (g/μg min). The obtained experimental qe values were closer to the values calculated with the pseudo-second-order kinetics, which means that mentioned model is more suitable for describing the sorption process (which is also confirmed by R2 to 1).
Table 5 clearly shows that the pseudo-second-order kinetic model describes the sorption kinetics much better than the pseudo-first-order model (R2 in the range of 0.9996–1.0000) for all the pharmaceuticals tested. Using this model, the maximum concentrations (qe) of the sorbed tested components that are in equilibrium were calculated for the environmental samples tested. These qe values range from 16.668–19.484 μg/g for all tested pharmaceuticals (0.796–19.484 μg/g for CEF; 0.989–16.668 μg/g for MEM; 0.685–18.649 μg/g for PRAZ).

3.5. Sorption Thermodynamics

The obtained Kd values are inserted into Equation (9) to describe the thermodynamics of the process by calculating the Gibbs free energy [48].
Δ G ° = R T l n K d
As shown in Table 6, cefdinir, memantine, and praziquantel were spontaneously sorbed to the soil and sediment samples studied by physisorption due to weak Van der Waals forces, which is confirmed with negative Δ G ° values ranging from −17.046 to 0.05 kJ/mol [43,49]. These statements are also consistent with the energy calculated from Dubinin–Radushkevich isotherms (see Section 3.2).

4. Conclusions

For the first time, to the author’s knowledge, the influence of organic matter has been studied as an important component of solid environmental samples such as soils and sediments through the sorption of three different pharmaceuticals. R2 values close to one for a linear isotherm showed good agreement between the experimental data and the model, whose considerations can also be applied to the sorption system studied. The addition of humic acids in the sorption process between pharmaceutical and sorbent promoted the sorption affinity of CEF, MEM and PRAZ with soils and sediments. Depending on the sorbent used and its physicochemical properties, a greater or lesser change in the partition coefficient was observed. In general, a higher content of organic matter in the solution promoted binding between the tested pharmaceuticals and all soil and sediment samples. A high sorption affinity can be observed for CEF, with the highest Kd values ranging from 31.84 to 972.7 mL/g, while PRAZ and MEM generally show a high sorption behavior of CEF, and as PRAZ and MEM are expected to be more exposed in the environment, which was confirmed by low Kd, KF and qm values. MEM clearly stands out as the most mobile component among the tested pharmaceuticals, tending to the aquatic environment with a Kd value below 38.5 mL/g. Thus, its fate in the environment does not depend on the sorption process, but rather on other biotic processes such as photolysis and hydrolysis. The good agreement between the experimental results and Freundlich isotherm model describes sorption for all tested pharmaceuticals as a complex process occurring in both one and multiple layers of the sorbents used. The pseudo-second-order kinetic model showed the best agreement with the obtained results with an R2 greater than 0.99. The thermodynamics were spontaneous, with CEF, MEM and PRAZ being sorbed by physisorption to the soil and sediment samples. The energy values determined using the Dubinin–Raduschkevich isotherm <8 kJ/mol confirmed the physical nature of the sorption.
Finally, this research confirmed the complexity of sorption, which is influenced by various process parameters. By focusing on the characterization of the sorbent and the pharmaceuticals, researchers can promote or reduce the mobility of contaminants in the environment. This and all similar research contribute to a better understanding of the various theoretical aspects of contaminant fate and behavior in different environments.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules27228008/s1, Figure S1: Kinetics for pharmaceuticals on the investigated soil and sediment samples: (a) CEF- soil 1 (b) CEF- sediment 1 (c) MEM- soil 5 (d) MEM- sediment 1 (e) PRAZ- soil 5 and (f) PRAZ- sediment 2. T = 25 °C; Figure S2: Linear sorption isotherms of CEF, MEM and PRAZ in five sediment samples and in five soil samples at 25 °C; Table S1: Selected pharmaceuticals, their structures and physico-chemical properties [33]; Table S2: Physico-chemical characterization of sediment and soil samples.

Author Contributions

Conceptualization, D.M.P., methodology, D.M.P.; software, D.M.P. validation, D.M.P. and K.T.Č.; formal analysis, D.M.P., K.T.Č., H.P., M.R.; investigation D.M.P., K.T.Č., H.P.; resources, D.M.P.; data curation, D.M.P. and K.T.Č.; writing—original draft preparation, D.M.P.; writing—review and editing, D.M.P. and K.T.Č.; visualization, D.M.P., K.T.Č., H.P.; supervision, D.M.P.; project administration, D.M.P.; funding acquisition, D.M.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Croatian Science Foundation under the project Fate of pharmaceuticals in the environment and during advanced wastewater treatment (PharmaFate) (IP-09-2014-2353).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

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

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Not available.

References

  1. Díaz-Cruz, S.; Barceló, D. Occurrence and analysis of selected pharmaceuticals and metabolites as contaminants present in waste waters, sludge and sediments. Handb. Environ. Chem. 2004, 5, 227–260. [Google Scholar]
  2. De Santiago-Martín, A.; Meffe, R.; Teijón, G.; Martínez Hernández, V.; López-Heras, I.; Alonso, C.; Arenas Romasanta, M.; de Bustamante, I. Pharmaceuticals and trace metals in the surface water used for crop irrigation: Risk to health or natural attenuation? Sci. Total Environ. 2020, 705, 135825. [Google Scholar] [CrossRef]
  3. Patel, M.; Kumar, R.; Kishor, K.; Mlsna, T.; Pittman, C.U., Jr.; Mohan, D. Pharmaceuticals of emerging concern in aquatic systems: Chemistry, occurrence, effects, and removal methods. Chem. Rev. 2019, 119, 3510–3673. [Google Scholar] [CrossRef] [Green Version]
  4. Díaz-Cruz, M.S.; López de Alda, M.J.; Barceló, D. Environmental behavior and analysis of veterinary and human drugs in soil, sediments and sludge. Trends Anal. Chem. 2003, 22, 340–351. [Google Scholar] [CrossRef]
  5. Filep, T.; Szabó, L.; Csaba Kondor, A.; Jakab, G.; Szalai, Z. Evaluation of the effect of the intrinsic chemical properties of pharmaceutically active compounds (PhACs) on sorption behaviour in soils and goethite. Ecotoxicol. Environ. Saf. 2021, 215, 112120. [Google Scholar] [CrossRef]
  6. Radović, T.T.; Grujić, S.D.; Kovačević, S.R.; Laušević, M.D.; Dimkić, M.A. Sorption of selected pharmaceuticals and pesticides on different river sediments. Environ. Sci. Pollut. Res. 2016, 23, 25232–25244. [Google Scholar] [CrossRef]
  7. Tolić, K.; Mutavdžić Pavlović, D.; Židanić, D.; Runje, M. Nitrofurantoin in sediments and soils: Sorption, isotherms and kinetics. Sci. Total Environ. 2019, 681, 9–17. [Google Scholar] [CrossRef]
  8. Li, J.; Hu, J.; Sheng, G.; Zhao, G.; Huang, Q. Effect of pH, ionic strength, foreign ions and temperature on the adsorption of Cu(II) from aqueous solution to GMZ bentonite. Colloids Surf. A Physicochem. Eng. Asp. 2009, 349, 195–201. [Google Scholar] [CrossRef]
  9. Lunelli Bachmann, S.A.; Calvete, T.; Amaral Féris, L. Caffeine removal from aqueous media by adsorption: An overview of adsorbents evolution and the kinetic, equilibrium and thermodynamic studies. Sci. Total Environ. 2021, 767, 144229. [Google Scholar] [CrossRef]
  10. Dada, A.O.; Olalekan, A.P.; Olatunya, A.M.; Dada, O. Langmuir, Freundlich, Temkin and Dubinin–Radushkevich: Isotherms Studies of Equilibrium Sorption of Zn2+ Unto Phosphoric Acid Modified Rice Husk. J. Appl. Chem. 2012, 3, 38–45. [Google Scholar]
  11. Nnaji, C.C.; Agim, A.E.; Mama, C.N.; Emenike, P.G.C.; Ogarekpe, N.M. Equilibrium and thermodynamic investigation of biosorption of nickel from water by activated carbon made from palm kernel chaf. Sci. Rep. 2021, 11, 7808. [Google Scholar] [CrossRef]
  12. Hasani, N.; Selimi, T.; Mele, A.; Thaçi, V.; Halili, J.; Berisha, A.; Sadiku, M. Theoretical, Equilibrium, Kinetics and Thermodynamic Investigations of Methylene Blue Adsorption onto Lignite Coal. Molecules 2022, 27, 1856. [Google Scholar] [CrossRef]
  13. Mutavdžić Pavlović, D.; Ćurković, L.; Blažek, D.; Župan, J. The sorption of sulfamethazine on soil samples: Isotherms and error analysis. Sci. Total Environ. 2014, 497–498, 543–552. [Google Scholar] [CrossRef] [Green Version]
  14. Northcott, G.L.; Jones, K.C. Experimental approaches and analytical techniques for determining organic compound bound residues in soil and sediment. Environ. Pollut. 2000, 108, 19–43. [Google Scholar] [CrossRef]
  15. Doretto, K.M.; Rath, S. Sorption of sulfadiazine on Brazilian soils. Chemosphere 2013, 90, 2027–2034. [Google Scholar] [CrossRef] [Green Version]
  16. Liu, X.; Zhang, H.; Luo, Y.; Zhu, R.; Wang, H.; Huang, B. Sorption of oxytetracycline in particulate organic matter in soils and sediments: Roles of pH, ionic strength and temperature. Sci. Total Environ. 2020, 714, 136628. [Google Scholar] [CrossRef]
  17. De Paolis, F.; Kukkonen, J. Binding of organic pollutants to humic acid and fulvic acids: Influence of pH and the structure of humic material. Chemosphere 1997, 34, 1693–1704. [Google Scholar] [CrossRef]
  18. Zhang, Y.L.; Lin, S.S.; Dai, C.M.; Shi, L.; Zhou, X.F. Sorption–desorption and transport of trimethoprim and sulfonamide antibiotics in agricultural soil: Effect of soil type, dissolved organic matter, and pH. Environ. Sci. Pollut. Res. 2014, 21, 5827–5835. [Google Scholar] [CrossRef]
  19. Wegst-Uhrich, S.R.; Navarro, D.A.G.; Zimmerman, L.; Aga, D.S. Assessing antibiotic sorption in soil: A literature review and new case studies on sulphonamides and macrolides. Chem. Cent. J. 2014, 8, 5. [Google Scholar] [CrossRef] [Green Version]
  20. Margon, A.; Cantone, P.; Pastrello, A. Interaction between Diclofenac and Soil Humic Acids. Soil Sediment Contam. 2009, 18, 489–496. [Google Scholar] [CrossRef]
  21. Lertpaitoonpan, W.; Ong, S.K.; Moorman, T.B. Effect of organic carbon and pH on soil sorption of sulfamethazine. Chemosphere 2009, 76, 558–564. [Google Scholar] [CrossRef]
  22. Chen, Y.; Hu, C.; Deng, D.; Li, Y.; Luo, L. Factors affecting sorption behaviors of tetracycline to soils: Importance of soil organic carbon, pH and Cd contamination. Ecotoxicol. Environ. Saf. 2020, 197, 110572. [Google Scholar] [CrossRef]
  23. Guo, X.; Tu, B.; Ge, J.; Yang, C.; Song, X.; Dang, Z. Sorption of tylosin and sulfamethazine on solid humic acid. J. Environ. Sci. 2016, 43, 208–215. [Google Scholar] [CrossRef]
  24. Chianese, S.; Fenti, A.; Iovino, P.; Musmarra, D.; Salvestrini, S. Sorption of Organic Pollutants by Humic Acids: A Review. Molecules 2020, 25, 918. [Google Scholar] [CrossRef] [Green Version]
  25. Havlíkova, L.; Satínský, D.; Solich, P. Aspects of decontamination of ivermectin and praziquantel from environmental waters using advanced oxidation technology. Chemosphere 2016, 144, 21–28. [Google Scholar] [CrossRef]
  26. Šagud, I.; Zanolla, D.; Perissutti, B.; Passerini, N.; Škorić, I. Identification of degradation products of praziquantel during the mechanochemical activation. J. Pharm. Biomed. 2018, 159, 291–295. [Google Scholar] [CrossRef]
  27. Dong, J.; Zhou, M.; Wu, X.; Du, M.; Wang, X. Memantine combined with environmental enrichment improves spatial memory and alleviates Alzheimer’s disease-like pathology in senescence-accelerated prone-8 (SAMP8) mice. J. Biomed. Res. 2012, 26, 439–447. [Google Scholar] [CrossRef] [Green Version]
  28. Das, N.; Madhavan, J.; Selvi, A.; Das, D. An overview of cephalosporin antibiotics as emerging contaminants: A serious environmental concern. 3 Biotech 2019, 9, 231–244. [Google Scholar] [CrossRef]
  29. Ribeiro, A.R.; Sures, B.; Schmidt, T.C. Cephalosporin antibiotics in the aquatic environment: A critical review of the occurrence, fate, ecotoxicity and removal technologies. Environ. Pollut. 2018, 241, 1153–1166. [Google Scholar] [CrossRef]
  30. Pan, M.; Chu, L.M. Fate of antibiotics in soil and their uptake by edible crops. Sci. Total Environ. 2017, 599-600, 500–512. [Google Scholar] [CrossRef]
  31. Tetzner, N.F.; Guedes Maniero, M.; Rodrigues-Silva, C.; Rath, S. On-line solid phase extraction-ultra high performance liquid chromatography-tandem mass spectrometry as a powerful technique for the determination of sulfonamide residues in soils. J. Chromatogr. A 2016, 1452, 89–97. [Google Scholar] [CrossRef] [PubMed]
  32. OECD. Adsorption-desorption using a batch equilibrium method. In OECD Guideline for the Testing of Chemicals 106; Organization for Economic Cooperation and Development: Paris, France, 2000. [Google Scholar]
  33. Mutavdžić Pavlović, D.; Tolić Čop, K.; Barbir, V.; Gotovuša, M.; Lukač, I.; Lozančić, A.; Runje, M. Sorption of cefdinir, memantine, praziquantel and trimethoprim in sediment and soil samples. Environ. Sci. Pollut. Res. 2022, 29, 66841–66857. [Google Scholar] [CrossRef] [PubMed]
  34. Yue, L.; Ge, C.; Feng, D.; Yu, H.; Deng, H.; Fu, B. Adsorption–desorption behavior of atrazine on agricultural soils in China. J. Environ. Sci. 2017, 57, 180–189. [Google Scholar] [CrossRef] [PubMed]
  35. Kim, Y.K.; Lim, S.J.; Han, M.H.; Cho, J.Y. Sorption characteristics of oxytetracycline, amoxicillin, and sulfathiazole in two different soil types. Geoderma 2012, 185–186, 97–101. [Google Scholar] [CrossRef]
  36. Mrozik, W.; Stefańska, J. Adsorption and biodegradation of antidiabetic pharmaceuticals in soils. Chemosphere 2014, 95, 281–288. [Google Scholar] [CrossRef] [PubMed]
  37. Thiele-Bruhn, S. Pharmaceutical antibiotic compounds in soils—A review. J. Plant Nutr. Soil Sci. 2003, 166, 145–167. [Google Scholar] [CrossRef]
  38. Babić, S.; Mutavdžić Pavlović, D. Analysis of PhACs in solid environmental samples (soil, sediment, sludge). In Analysis, Removal, Effects and Risk of Pharmaceuticals in the Water Cycle: Occurrence and Transformation in the Environment; Elsevier: Amsterdam, The Netherlands, 2013; pp. 129–169. [Google Scholar]
  39. Pereira Leal, R.M.; Ferracciú Alleoni, L.R.; Tornisielo, V.L.; Borges Regitano, J. Sorption of fluoroquinolones and sulfonamides in 13 Brazilian soils. Chemosphere 2013, 92, 979–985. [Google Scholar] [CrossRef] [Green Version]
  40. Mutavdžić Pavlović, D.; Ćurković, L.; Grčić, I.; Šimić, I.; Župan, J. Isotherm, kinetic, and thermodynamic study of ciprofloxacin sorption on sediments. Environ. Sci. Pollut. Res. 2016, 24, 10091–10106. [Google Scholar] [CrossRef]
  41. Wang, Z.; Jiang, Q.; Wang, R.; Yuan, X.; Yang, S.; Wang, W.; Zhao, Y. Effects of Dissolved Organic Matter on Sorption of Oxytetracycline to Sediments. Geofluids 2018, 2018, 1254529. [Google Scholar] [CrossRef] [Green Version]
  42. Yang, F.; Wang, M.; Wang, Z. Sorption behavior of 17 phthalic acid esters on three soils: Effects of pH and dissolved organic matter, sorption coefficient measurement and QSPR study. Chemosphere 2013, 93, 82–89. [Google Scholar] [CrossRef]
  43. Xu, Y.; Yu, X.; Xu, B.; Peng, D.; Guo, X. Sorption of pharmaceuticals and personal care products on soil and soil components: Influencing factors and mechanisms. Sci. Total Environ. 2021, 753, 141891. [Google Scholar] [CrossRef] [PubMed]
  44. Kulshrestha, P.; Giese, R.F., Jr.; Aga, D.S. Investigating the molecular interactions of oxytetracycline in clay and organic matter: Insights on factors afecting its mobility in soil.Environ. Sci. Technol. 2004, 38, 4097–4105. [Google Scholar] [CrossRef] [PubMed]
  45. Ling, W.; Xu, J.; Gao, Y.; Wang, H. Infuence of dissolved organic matter (DOM) on environmental behaviors of organic pollutants in soils. Chin. J. Appl. Ecol. 2004, 15, 326–330. [Google Scholar]
  46. Namasivayam, C.; Kavitha, D. Removal of Congo Red from water by adsorption onto activated carbon prepared from coir pith, an agricultural solid waste. Dyes Pigments 2002, 54, 47–58. [Google Scholar] [CrossRef]
  47. Pholosi, A.; Naidoo, E.B.; Ofomaja, A.E. Intraparticle diffusion of Cr(VI) through biomass and magnetite coated biomass: A comparative kinetic and diffusion study. S. Afr. J. Chem. Eng. 2020, 32, 39–55. [Google Scholar] [CrossRef]
  48. Maszkowska, J.; Wagil, M.; Mioduszewska, K.; Kumirska, J.; Stepnowski, P.; Białk-Bielińska, A. Thermodynamic studies for adsorption of ionizable pharmaceuticals onto soil. Chemosphere 2014, 111, 568–574. [Google Scholar] [CrossRef]
  49. Xiang, L.; Xiao, T.; Mo, C.H.; Zhao, H.-M.; Li, Y.-W.; Li, H.; Cai, Q.-Y.; Zhou, D.-M.; Wong, M.-H. Sorption kinetics, isotherms, and mechanism of aniline aerofloat to agricultural soils with various physicochemical properties. Ecotoxicol. Environ. Saf. 2018, 154, 84–91. [Google Scholar] [CrossRef]
Figure 1. Characteristic MRM mass spectra of (A) standard solution of CEF: 396.0→226.9; MEM: 180.2→163.1; PRAZ: 313.0→203.0 and (B) investigated sediment matrix (Glina).
Figure 1. Characteristic MRM mass spectra of (A) standard solution of CEF: 396.0→226.9; MEM: 180.2→163.1; PRAZ: 313.0→203.0 and (B) investigated sediment matrix (Glina).
Molecules 27 08008 g001
Figure 2. Soil and sediment sorption coefficients for tested pharmaceuticals.
Figure 2. Soil and sediment sorption coefficients for tested pharmaceuticals.
Molecules 27 08008 g002
Table 1. The linear, Freundlich, and Dubinin–Radushkevich isotherm parameters for CEF.
Table 1. The linear, Freundlich, and Dubinin–Radushkevich isotherm parameters for CEF.
IsothermHA0 mg/L1 mg/L10 mg/L100 mg/L
Sediment/
Soil
12345123451234512345
12345123451234512345
LinearKd143.2393.7113.047.30117.3143.3399.7116.658.55242.4151.9413.0197.776.40341.1233.1720.4351.0100.8438.5
31.8410.69144.5297.6350.352.6824.39147.2304.2355.484.8257.53284.4499.9573.4176.792.26444.9813.1972.7
R20.9960.9980.9950.9960.9910.9980.9980.9960.9960.9940.9930.9930.9980.9990.9320.9980.9920.9920.9990.995
0.9920.9900.9980.9910.9980.9900.9910.9930.9920.9950.9930.9940.9930.9960.9900.9920.9920.9890.9920.990
FreundlichKF202.33342.092.2634.36169.0182.4261.854.5834.67105.482.04261.885.9031.33181.66854.94.44 × 1074.92 × 104353.23.54 × 105
23.619.89264.22741.62944.428.4527.54208.9335.7234.441.2132.58107.6295.8386.41798.9223.362.61 × 1052.35 × 1084.62 × 109
n0.870.571.141.370.900.911.151.581.591.481.271.201.401.751.350.390.200.330.570.27
1.521.860.790.550.571.711.030.900.971.111.591.541.441.181.090.470.640.280.180.16
R20.9380.9450.9530.9760.9410.9750.9940.9920.9840.9870.9410.9870.9760.9480.9750.9320.8220.8620.9560.894
0.9520.9790.9870.9230.9310.9620.9780.9870.9930.9860.9650.9620.9530.9820.9850.8740.8900.8080.8350.839
D-Rqm40.38168.1726.0014.7240.0636.1437.8516.9914.5225.7317.9940.3319.6211.4635.02294.43.07 × 103827.163.192.27 × 103
12.916.6245.31171.2165.111.0915.5440.1943.1033.0314.4112.0623.8541.3443.84151.347.932.13 × 1036.71 × 1043.24 × 105
β0.0300.0400.0220.0200.0300.0270.0150.0120.0150.0120.0150.0140.0130.0110.0130.0790.1400.0880.0610.11
0.0200.0170.0330.0460.0420.0120.0400.0270.0200.0160.0130.0140.0120.0140.0150.0680.0540.1000.150.17
E4.103.554.815.014.064.285.776.455.876.485.795.936.276.846.302.521.892.382.872.13
5.015.443.883.303.456.513.544.275.045.636.255.956.485.985.792.723.042.241.831.71
R20.9470.9620.9510.9340.9700.9750.9780.9370.9390.9780.8490.9840.9110.8280.9790.9590.8390.8820.9890.911
0.9280.8350.9980.9550.9530.8330.9900.9990.9760.9480.8940.8480.9210.9740.9650.9010.9040.8300.8450.849
Kd—mL/g; KF—(μg/g)(mL/μg)1/n; qm—μg/g; β—mol2/kJ2; E—kJ/mol. Blue numbers correspond to sediment samples; black number correspond to soil samples (see the list of used sediment and soil samples in Section 2.2).
Table 2. The linear, Freundlich, and Dubinin–Radushkevich isotherm parameters for MEM.
Table 2. The linear, Freundlich, and Dubinin–Radushkevich isotherm parameters for MEM.
IsothermHA0 mg/L1 mg/L10 mg/L100 mg/L
Sediment/
Soil
12345123451234512345
12345123451234512345
LinearKd5.112.492.971.430.986.574.446.303.843.158.7112.428.114.233.209.1113.1012.057.248.54
5.760.755.642.6212.669.991.108.814.2024.2912.422.5811.025.9531.7112.913.1512.006.1338.51
R20.9930.9910.9910.9920.9930.9940.9920.9960.9920.9900.9930.9940.9900.9900.9950.9930.9960.9980.9900.998
0.9930.9940.9930.9910.9910.9970.9930.9970.9910.9960.9930.9920.9910.9910.9900.9940.9940.9930.9900.991
FreundlichKF4.322.683.131.641.036.504.216.343.342.677.109.797.524.513.4212.686.9312.595.568.07
4.850.915.012.827.758.341.477.414.0321.239.803.299.686.3713.099.862.999.755.0966.07
n1.811.992.051.911.591.761.541.481.531.092.822.122.221.971.670.442.400.782.861.22
1.951.861.751.973.211.431.881.841.681.342.121.862.162.042.761.511.191.461.320.595
R20.8650.8880.8890.9180.8930.9700.9490.9820.8810.8760.9150.9490.9500.9810.9300.9170.7460.9890.8420.991
0.8660.8900.8930.8820.8670.9490.9520.9380.9060.9940.9490.9760.9860.9770.8570.9500.9760.9500.9040.988
D-Rqm3.252.342.681.521.014.833.454.912.702.274.886.115.144.042.8512.584.359.283.966.07
3.611.153.812.415.145.561.425.023.189.906.112.996.454.766.776.442.826.153.7425.16
β0.0200.0230.0210.0290.0370.0210.0300.0300.0290.0500.00770.0120.0120.0210.0270.170.00860.07330.00770.039
0.0180.0310.0220.0220.00550.0270.0310.0170.0250.0200.0120.0270.0130.0170.00550.0230.0510.0240.0330.078
E4.984.654.864.183.674.854.084.104.123.168.066.486.404.884.281.737.622.618.063.57
5.354.004.764.749.534.344.045.394.504.986.484.276.255.469.534.693.134.553.892.53
R20.5840.6370.6480.6480.6020.7820.7340.8220.6040.6180.6850.7500.7320.8150.6860.9700.4820.9470.5750.869
0.6040.5950.6590.6180.6420.7530.7170.7140.6460.8820.7500.8330.8240.7900.6460.7670.8110.7650.6540.990
Kd—mL/g; KF—(μg/g)(mL/μg)1/n; qm—μg/g; β—mol2/kJ2; E—kJ/mol. Blue numbers correspond to sediment samples; black number correspond to soil samples (see the list of used sediment and soil samples in Section 2.2).
Table 3. The linear, Freundlich, and Dubinin–Radushkevich isotherm parameters for PRAZ.
Table 3. The linear, Freundlich, and Dubinin–Radushkevich isotherm parameters for PRAZ.
IsothermHA0 mg/L1 mg/L10 mg/L100 mg/L
Sediment/
Soil
12345123451234512345
12345123451234512345
LinearKd11.732.293.681.018.0515.236.0612.323.699.3918.237.5213.214.479.7231.4613.4515.7611.0911.03
35.849.198.410.2727.7569.709.6214.816.8034.5870.299.7920.616.9035.41116.429.9621.439.6152.12
R20.9930.9980.9940.9910.9970.9930.9920.9920.9930.9940.9900.9920.9950.9950.9930.9970.9950.9920.9940.997
0.9920.9960.9860.9920.9970.9900.9930.9920.9950.9890.9940.9940.9990.9910.9980.9960.9960.9950.9910.996
FreundlichKF17.742.224.091.308.2416.566.8913.093.768.7117.467.8714.064.9510.4236.8114.9317.5011.6411.07
38.6411.1714.060.4939.45115.0811.3817.227.0343.65101.1612.0223.287.6440.00201.849.6823.9910.0263.53
n0.520.602.071.741.001.101.241.171.520.990.970.941.001.010.600.930.960.991.000.97
0.920.530.522.580.690.800.910.991.250.920.950.960.990.980.980.780.980.971.290.92
R20.8510.9250.9260.9430.9810.9780.9920.9810.9580.9640.9720.9540.9800.9760.9880.9820.9810.9700.9840.990
0.9470.9320.7380.9350.9350.9200.9070.9560.9880.9380.9460.9570.9880.9750.9820.9470.9920.9780.9540.967
D-Rqm16.923.113.251.276.9811.486.019.703.276.4611.056.9810.105.018.6118.0010.6511.708.578.29
17.9510.641.060.5020.4035.995.1812.515.8121.6231.4210.5813.916.8318.4041.877.4314.438.1024.10
β0.140.140.0180.0350.0580.0420.0440.0410.0340.0530.0480.0650.0510.0660.110.0430.0530.0480.0530.055
0.0430.140.150.0250.0720.0430.0050.0520.0420.0440.0390.120.0550.0610.0410.0370.0550.0470.0380.038
E1.901.917.473.792.943.433.363.483.863.083.232.783.142.752.153.393.063.233.093.02
3.411.922.564.482.643.4110.433.103.473.373.602.053.022.873.493.663.013.283.623.63
R20.9370.9720.6710.7180.9380.9790.9470.9700.7720.8500.9360.9260.9590.9620.9540.9900.9670.9590.9330.937
0.9410.9630.8810.6750.9530.9580.9230.9800.9040.9760.9710.9730.9810.9530.9840.9600.9240.9870.9350.986
Kd—mL/g; KF—(μg/g)(mL/μg)1/n; qm—μg/g; β—mol2/kJ2; E—kJ/mol. Blue numbers correspond to sediment samples; black number correspond to soil samples (see the list of used sediment and soil samples in Section 2.2).
Table 4. Change in Pearson’s correlation coefficients between soil/sediment physicochemical properties and Kd value of studied pharmaceuticals with change in OM concentration (from 0 to 100 mg/L humic acid).
Table 4. Change in Pearson’s correlation coefficients between soil/sediment physicochemical properties and Kd value of studied pharmaceuticals with change in OM concentration (from 0 to 100 mg/L humic acid).
ParameterCEFMEMPRAZ
Coarse sand −0.490 to −0.459−0.408 to −0.3090.209 to 0.109
Clay −0.278 to −0.233−0.127 to −0.1620.665 to 0.640
Silt 0.393 to 0.4540.223 to −0.003−0.371 to −0.288
Fine sand 0.490 to 0.4580.408 to 0.310−0.209 to −0.109
pH 0.125 to −0.0230.616 to 0.3770.228 to 0.382
EC 0.198 to 0.1000.753 to 0.370−0.274 to 0.133
TDS * 0.198 to 0.1000.753 to 0.370−0.274 to 0.133
CEC −0.184 to −0.2130.413 to −0.138−0.102 to 0.014
OM −0.257 to −0.258−0.035 to −0.279−0.015 to −0.111
CaCO3−0.182 to −0.189−0.404 to −0.2270.404 to −0.007
Zn −0.077 to −0.117−0.023 to −0.031−0.179 to −0.232
Cu −0.071 to −0.085−0.005 to −0.063−0.327 to −0.357
Fe 0.210 to 0.1840.093 to 0.175−0.556 to −0.430
Mn −0.007 to −0.0620.232 to 0.029−0.570 to −0.420
TDS * = Total Dissolved Solids.
Table 5. Sorption kinetic parameters of pharmaceuticals for soil 5.
Table 5. Sorption kinetic parameters of pharmaceuticals for soil 5.
Initial Concentration, mg/Lqe,exp, μg/gPseudo-First-OrderPseudo-Second-Order
qe,calc, μg/gk1, min−1R2qe,calc, μg/gk2, g/μg minR2
CEF2.019.4842.19181.15·10−30.806519.53130.0111.0000
0.54.7950.84841.15·10−30.90294.81930.0240.9999
0.10.7960.23491.2·10−40.36723.90780.0391.0000
MEM2.016.6685.74384.606·10−40.794016.75040.006400.9999
0.54.7720.83851.152·10−30.73054.78470.029061.0000
0.10.9890.06331.382·10−20.73460.98960.382271.0000
PRAZ2.018.6495.5071.152·10−30.906418.8320.033111.0000
0.54.0111.6494.606·10−40.82564.0360.021570.9997
0.10.6850.3771.612·10−40.41060.6860.343420.9996
Table 6. Standard Gibbs free energy calculated for different concentration of HA.
Table 6. Standard Gibbs free energy calculated for different concentration of HA.
ΔG°, kJ/mol
0 mg/L HA1 mg/L HA10 mg/L HA100 mg/L HA
CEF
sediment 1−12.30−12.30−12.45−13.51
sediment 2−14.80−14.84−14.92−16.30
sediment 3−11.71−11.79−13.10−14.52
sediment 4−9.55−10.08−10.74−11.43
sediment 5−11.81−13.60−14.59−15.07
soil 1−8.57−9.82−11.00−12.82
soil 2−5.87−7.91−10.04−11.21
soil 3−12.32−12.37−14.00−15.11
soil 4−14.11−14.17−15.40−16.60
soil 5−14.52−14.55−15.74−17.05
MEM
sediment 1−4.04−4.66−5.36−5.47
sediment 2−2.26−3.69−6.24−6.37
sediment 3−2.70−4.56−5.19−6.16
sediment 4−0.89−3.33−3.57−4.90
sediment 50.05−2.84−2.88−5.31
soil 1−4.34−5.70−6.24−6.34
soil 20.75−0.24−2.35−2.84
soil 3−4.29−5.39−5.95−6.16
soil 4−2.39−3.56−4.42−4.49
soil 5−6.29−7.90−8.56−9.05
PRAZ
sediment 1−6.10−6.75−7.19−8.54
sediment 2−2.06−4.46−5.00−6.44
sediment 3−3.23−6.22−6.39−6.83
sediment 4−0.02−3.23−3.71−5.96
sediment 5−5.17−5.55−5.63−5.95
soil 1−8.87−10.52−10.54−11.79
soil 2−5.50−5.61−5.65−5.69
soil 3−5.28−6.68−7.50−7.59
soil 43.24−4.75−4.79−5.61
soil 5−8.23−8.78−8.84−9.80
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mutavdžić Pavlović, D.; Tolić Čop, K.; Prskalo, H.; Runje, M. Influence of Organic Matter on the Sorption of Cefdinir, Memantine and Praziquantel on Different Soil and Sediment Samples. Molecules 2022, 27, 8008. https://doi.org/10.3390/molecules27228008

AMA Style

Mutavdžić Pavlović D, Tolić Čop K, Prskalo H, Runje M. Influence of Organic Matter on the Sorption of Cefdinir, Memantine and Praziquantel on Different Soil and Sediment Samples. Molecules. 2022; 27(22):8008. https://doi.org/10.3390/molecules27228008

Chicago/Turabian Style

Mutavdžić Pavlović, Dragana, Kristina Tolić Čop, Helena Prskalo, and Mislav Runje. 2022. "Influence of Organic Matter on the Sorption of Cefdinir, Memantine and Praziquantel on Different Soil and Sediment Samples" Molecules 27, no. 22: 8008. https://doi.org/10.3390/molecules27228008

APA Style

Mutavdžić Pavlović, D., Tolić Čop, K., Prskalo, H., & Runje, M. (2022). Influence of Organic Matter on the Sorption of Cefdinir, Memantine and Praziquantel on Different Soil and Sediment Samples. Molecules, 27(22), 8008. https://doi.org/10.3390/molecules27228008

Article Metrics

Back to TopTop