Next Article in Journal
Antibacterial and Antibiofilm Activity of Ficus carica-Mediated Calcium Oxide (CaONPs) Phyto-Nanoparticles
Next Article in Special Issue
Impact of Exposure to Pyraclostrobin and to a Pyraclostrobin/Boscalid Mixture on the Mitochondrial Function of Human Hepatocytes
Previous Article in Journal
Biological Activity and Structural Diversity of Steroids Containing Aromatic Rings, Phosphate Groups, or Halogen Atoms
Previous Article in Special Issue
Using Glyphosate on Guarana Seedlings in the Amazon
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Soil Organic Matter Composition and pH as Factors Affecting Retention of Carbaryl, Carbofuran and Metolachlor in Soil

by
Irmina Ćwieląg-Piasecka
Institute of Soil Science, Plant Nutrition and Environmental Protection, Wroclaw University of Environmental and Life Sciences, Grunwaldzka 53 St., 50-357 Wrocław, Poland
Molecules 2023, 28(14), 5552; https://doi.org/10.3390/molecules28145552
Submission received: 14 June 2023 / Revised: 18 July 2023 / Accepted: 19 July 2023 / Published: 20 July 2023

Abstract

:
The majority of studies concerning the environmental behavior of hydrophobic pollutants in soil consider soil organic matter (SOM) content as a main factor influencing chemical retention, whereas the composition of SOM and its individual fraction share are often neglected. In the present paper, carbaryl, carbofuran and metolachlor retention by loamy sand and loam topsoil materials is compared and referred to humic acids (CHA) and the residual carbon (CR) content of SOM. Additionally, the sorption-desorption behavior of agrochemicals in soils was tested at a pH of three to seven. Calculated isothermal parameters point to favorable, spontaneous and physical pesticide sorption. Groundwater ubiquity score (GUS) indexes confirmed the low leaching ability of metolachlor on soils and moderate of carbofuran. The high affinity of carbaryl to CR may explain its pronounced sorption in loam soil and the lowest percolation potential. Carbofuran retention in soils was associated with montmorillonite (Mt) and CR fractions. Meanwhile, metolachlor uptake was related to humic acid and Mt content of the soils. Lower pH enhanced retention of the agrochemicals, except for carbaryl sorption in sandy loam soil. Results of this study highlight that SOM composition and mutual share of individual organic carbon fractions alongside pH may play a crucial role in predicting non-ionic pesticide behavior in soil.

1. Introduction

The effect of soil composition on pesticide sorption has been widely studied since these substances were introduced into the environment [1,2,3,4]. There are many pedological factors governing their retention. Among them, native organic matter content [4,5] as well as the presence of exogenous sources of OM (biochar, compost, peat) [6,7,8] or inorganic colloidal soil fractions (clay minerals–CM, oxides) [9,10,11] are believed to be decisive in many pesticides immobilization in soil. It has been postulated that in soils of low organic carbon content, clay minerals and hydr(oxides) may significantly contribute to the sorption of particular agrochemicals [12]. In addition to that, an important factor in determining the availability of potential sorption sites is soil pH, as it controls the form and thus reactivity of both the potential sorbate and various inorganic and organic soil colloids. Nevertheless, in studies of non-ionic organic compounds, the pH influence is often neglected, as these neutral substances are present in their molecular forms in the whole pH range. Meanwhile, pH modulates surface charges in soil, together with solubility of various soil colloids, and is, at the same time, a limiting factor in terms of the possible interaction types and sorption strength among the pollutant and soil constituents.
According to many authors, uncharged chemicals are primarily bound to soil organic matter [13,14,15] and hydrophobic bonding is postulated as the main route of their uptake [14,16,17]. Several mechanisms have been proposed for the interaction of non-ionic chemicals with humic substances (HS) in soil [18,19], such as van der Waals forces, hydrogen bonding, ligand exchange and charge transfer complexes. They may operate during the interaction between a given pesticide and soil particles depending on the content of particular SOM fractions and the nature of the organo-mineral associations [5].
Recently, more attention has been paid to the composition of SOM in terms of its role as an indicator of landscape changes [20], soil physical quality [21], cultivation practices on arable lands [16,22], as well as predictors of pesticide affinity to soil [23]. Despite these paramount findings and their importance, by 2020, Novotny et al. [23] counted only a few papers referring to the relationship between the chemical structure of SOM and the organic contaminant sorption capacity of the whole soil. According to them, particular SOM fractions (e.g., humic acids, aromatic and aliphatic structures) were found to be useful proxies in understanding the interaction between SOM and pesticides [24,25,26,27]. However, to obtain a complete overview of the individual SOM fractions on the whole soil capacity towards organic chemicals, a joint approach is necessary. It should take into account not only the agrochemical affinity to individual SOM components but also the content and mutual share of particular organic carbon fractions of the soil, such as total extractable (TEC), humic acids (CHA), or residual carbon (CR) fractions. Information on the behavior of organic contaminants in soils related strictly to organic carbon (OC) composition and its individual fraction content is still limited, although it might be decisive in elucidating the soil capacity towards agrochemicals.
Non-ionic pesticides, namely carbaryl, carbofuran and metolachlor, were chosen for the study as their properties are intermediate between polar and non-polar compounds. They vary in water solubility, the type of functional groups present in their structure and hydrophobicity. Carbaryl and carbofuran belong to the carbamate insecticides group. Both compounds are non-persistent, moderately mobile in soil and of low (carbaryl) to moderate (carbofuran) water solubility, which strongly influences their environmental fate, including leaching potential [10,12,28,29]. Despite their suspected endocrine activity and disapproval for application in certain EU countries, they are still used in many parts of the world and detected in soils of some central and eastern European regions, including Poland [30]. Metolachlor is a selective pre-emergency herbicide [31], which belongs to the aniline derivatives group. Its moderate adsorption in soil is correlated with relatively high water solubility of the chemical [32], while its high lipophilicity places it in the group of substances suspected of endocrine activity [33], similar to carbaryl and carbofuran.
Both carbamate and metolachlor retention in soil is attributed to its organic carbon content. Carbaryl sorption in soil is mostly connected to its SOM content, as well as illite (Il) and goethite (Ge), all of which have been shown to greatly increase pesticide retention when added to soils [10]. Studying model clay minerals and hydroxides, Chen et al. [9] discovered that carbaryl was efficiently adsorbed on montmorillonite (Mt), kaolinite (K) and goethite. Carbofuran sorption in soil was also related to organic matter and clay content [12], although Singh and Srivastava [29] claim that the mineral content of the soil is more important than that of OC when considering agrochemical retention. Indeed, in some other studies it was found that goethite and montmorillonite were the major soil constituents responsible for carbofuran retention in loamy sand soil of low organic carbon content [10]. Organic matter was also postulated as the predominant adsorbent for metolachlor, which was sufficiently mobile in soils low in organic carbon [34]. In the previous study it was found that the sorptive potential of soil towards the chemical was greatly enhanced by its amendment with montmorillonite, humic acids, goethite and illite [10]. Meanwhile, studies on the sorption of these two groups of non-ionic pesticides under varying soil pH conditions are scarce and often contradictory. Rama Krishna and Philip [35] reported that adsorption of neutral molecules increases with an increase in pH, which was corroborated by research on carbofuran in lateritic soils [12]. Metolachlor sorption was the most effective at a soil pH of five [36], whilst Mondal et al. [37] reported that carbaryl was efficiently retained in alluvial soils at a pH of six. In general, the majority of carbamates studies considering their adsorption capacity under various pHs of soils are contradictory and limited, as most of them are related to the chemical hydrolysis of the insecticides under neutral and alkaline pH conditions [38,39,40,41,42].
Considering their importance, this study was aimed at investigating the sorption of non-ionic pesticides with endocrine activity (carbaryl, carbofuran, metolachlor) on agricultural soils of different SOM fractional compositions and under various pH ranges. In order to comprehensively analyze the environmental behavior of a particular pesticide in soil, the study was performed on the whole soil material, not a single-component system. The evaluation was undertaken based on pesticide sorption-desorption studies performed on loamy sand and loam topsoil materials and an estimation of their release from the tested soils under a pH range of three to seven. Isothermal parameters obtained from Langmuir, Freundlich, Temkin and Dubinin–Radushkievich models were then related to individual SOM fraction contents and shares in the OC of the soils. The overall aim was to determine which SOM fraction represents the potential sink for the tested non-ionic agrochemicals in soils and to establish the pH conditions for their most effective retention. Such information may significantly contribute to the risk assessment of surface and groundwater contamination by the tested chemicals and fill the gaps in understanding of the attenuation phenomena of organic xenobiotic off-site migration in the environment.

2. Results and Discussion

2.1. Sorption Isotherms and Their Mathematical Approximations

Tested pesticide classes were subjected to sorption studies on sandy loam (L) and loam (C) topsoil materials and resulting isotherms (Figure 1) were approached with Langmuir, Freundlich Temkin and Dubinin–Radushkievich (D–R) models (Figure 2). Linear sorption parameters calculated for the obtained Langmuir-type [43] isotherms are summarized in Table 1. Applied models revealed high correlation coefficients (R2 > 0.950), except for the D–R model. The resulting L-shape of the isotherms is due to adsorption site limitations with the increasing pesticide concentration and presumably the energetic heterogeneity of the studied soils, as in the case of most solids [44]. In fact, the experimental isotherms approached a saturation upper limit but did not reach a complete plateau. This implies that at high adsorbate concentrations, weak intermolecular forces among the studied non-ionic pesticides might occur and participate in multilayer formation [45,46].
The maximum monolayer adsorption capacity Qmax (mg kg−1) of the Langmuir model, evaluated for the hydrophobic pesticides on soil L, correlated well with their Freundlich Kf constants and followed the increasing order of carbaryl (15.25) < carbofuran (27.93) < metolachlor (32.47). C soil Qmax parameters for carbaryl and carbofuran were doubled compared to L soil, whereas no significant change in case of metolachlor was noticed.
Analogous trends to those of Qmax were observed on both studied soils for the theoretical monolayer sorption capacity (Qm) calculated with the D–R approach. It points to a significant increase in the sorption sites active in insecticide retention. In addition, the hydrophobic character and low water solubility of the studied carbamates may favor their diffusion within organic matter structures. According to the literature, diffusion and the following physical entrapment within the soil–organic matter matrix is a common sorption mechanism, especially in the case of pesticides independent of any traditional chemical adsorption mechanisms [47].
Freundlich and Langmuir isotherm parameters, namely, the adsorption intensity (nf) and separation factor (RL), ranged from 0.248 to 0.611 (Table 1) and up to 0.80 (Figure 3), respectively. This indicates favorable (0.1 < nf < 1.0 [48]) and reversible (0 < RL < 1 [49]) adsorption of the solute and suggests that the pesticide molecules preferred high-energy sites followed by sorption at lower-energy sites [50,51].
The Temkin isotherm constant related to the heat of adsorption (bT), and the binding constant of the pesticide to a soil (A) exhibited decreasing and increasing trends, respectively, when comparing soils L and C (Table 1). According to the model assumptions, bT would more often decrease than increase with increasing coverage [52]. Hence, it can be concluded that the lower values of the constant obtained for loamy soil indicate a higher degree of surface coverage with the pesticides compared to sandy loam soil material. Furthermore, from the bT value, the nature of the process may be elucidated [53]. In the case of physical sorption, it should be lower than 4.184 kJ mol−1 (1 kcal mol−1), while for chemical sorption it is above 41,840 kJ mol−1 (10 kcal mol−1). Thus, the values of the bT parameter (Table 1) suggest that the sorption of carbamates and metolachlor in L and C soils is of a physical nature. This was also confirmed by the Ea values obtained within D–R isotherms (Table 1), which were lower than 8 kJ mol−1 and may be partly due to the van der Waals interaction between the pesticides and the hydrophobic soil components.
Freundlich Kf parameters, which revealed the highest variation when comparing both studied soils, were used for the calculation of the standard Gibb’s free energy change (ΔG°) of the tested pesticides (Table 2). Prior to ΔG° calculations, Kf parameters were converted to the dimensionless standard equilibrium constants using the appropriate equations, as suggested in the literature [54,55,56].
The negative ΔG° values obtained for all studied pesticides on soils L and C point to the spontaneous character of the sorption. The significant increase and the highest negative ΔG° value obtained for carbaryl sorption in C soil indicates that in the case of this insecticide, the adsorption capacity would be higher for clay than for sandy soil, as for carbofuran (Table 2). Meanwhile, ΔG° of much similar magnitude was obtained for metolachlor sorption in both studied soils. The values of this thermodynamic parameter calculated for the studied pesticides are comparable to those found in the literature, although a limited number of papers performing ΔG° calculations using dimensionless Kc can still be found [55,57,58].
In addition, Kf parameters were used to calculate organic carbon partition coefficient (Koc) and groundwater ubiquity score (GUS, [59]) indexes, providing information about pesticide mobility and leachability, respectively (Table 3). All calculated Koc values (except for carbaryl sorption in C soil) were below 500, which points to minimum or no adsorption of the chemicals into the soil and indicates high probability of their runoff or leaching [1,60]. A similar mobility ranking by FAO [61] based on log Koc values classifies all of the pesticides as of moderate mobility, achieving log Koc values in the range of two to three. Interestingly, despite the fact that the metolachlor sorption magnitude is higher in the loam than in the sandy loam soil under study, the Koc value for herbicides in C soil is lower than in L soil. This might suggest that organic matter is not the main soil constituent responsible for metolachlor immobilization in soil.
Calculated ground ubiquity score (GUS) indexes point to the following decreasing order of leaching threat in L soil: carbofuran > carbaryl > metolachlor and C soil: carbofuran > metolachlor > carbaryl. The lowest risk of soil C percolation (GUS value ≤ 1.8) was obtained for carbaryl and metolachlor and the chloroacetanilide pesticide in L soil. Carbofuran potential for movement toward groundwater was moderate in soil C (2.47), approaching high (> 2.80) in loamy sand soil.
The retention of carbaryl, carbofuran and metolachlor on the investigated loamy sand and loam soil materials was evaluated by comparing their sorption-desorption magnitudes. Studied pesticide trends in L and C soils are presented in Figure 4. They vary in the magnitude of sorption on loam and sandy loam soil studied for a given hydrophobic agrochemical. Nevertheless, their common feature is a decline in sorption magnitude with the increase in initial concentration, which confirms the resulting L-type shape of isotherms (Figure 1).
In general, carbaryl was adsorbed up to 44.51% and 71.07% in L and C soils, respectively (Figure 4a). Carbofuran was characterized by a similar pattern; however, its maximum sorption magnitude was 33.50% of pesticide initial dose introduced to soil L and 55.0% in the case of soil C (Figure 4b). Both carbamates revealed a great discrepancy in sorption magnitude among sandy loam and loam soils. On the contrary, metolachlor sorption in both soils was very similar, with the highest magnitude among the studied pesticides (84.06% and 71.96% for soils C and L, respectively, Figure 4c).
Meanwhile, desorption of carbamates and metolachlor achieved higher values in L rather than the C soil studied (Figure 5, expressed as a mass initially sorbed and desorbed). Carbaryl and carbofuran release from loamy sand soil was increasing along with the initial pesticide dose introduced to soil, achieving a maximum value of 62.6% and 87.0% of the pesticide dose initially sorbed, respectively.
Desorption in C soil was lower in the case of carbaryl (20.50–69.06% of the initial dose) and much higher for carbofuran in loam soil (from 49.76% of the lowest dose applied, up to approx. 74.70% in case of the highest initial concentration), making its sorption reversible [63] and enhancing a great threat of carbofuran percolation in soil. These results are consistent with the calculated GUS indexes (Table 3), confirming the low and slightly moderate potential of carbaryl to leach in soils C and L, respectively, and moderate potential for carbofuran in both studied soils. Metolachlor desorption in L soil was in a range of 31.20% to 45.89% and from 10.49% to 13.87% in C soil. This acetanilide pesticide release was the lowest among the group of the uncharged pesticides studied, as were the calculated GUS values (1.46 and 1.60 in soils L and C, respectively), suggesting its scarce leaching potential in both studied soils (Table 3). The strong adsorption and low desorption characteristics of metolachlor suggest that for the soils studied, the migration risk to groundwater should be lower than for the other uncharged pesticides examined.
In general, studied pesticides sorbed less to loamy sand than the loam soil (Figure 4), which contained more organic carbon and clay fractions.The highest negative standard Gibb’s free energy value obtained for carbaryl in soil C (Table 2) together with its lowest GUS index (1.12, Table 3) may indicate that the forces holding the insecticide in soil are relatively strong. Similarly, the relatively high and negative ΔG°, the lowest desorption magnitude (Figure 5e,f) and the highest log p value among the studied pesticides predestinate metolachlor to be efficiently retained in both soils.
Taking into account the parameters approximated from sorption isotherms as well as the results of desorption studies, it can be concluded that the molecular interactions between the studied pesticides and the active soil constituents are of a physical nature. In general, this type of sorption for uncharged pesticides in soils is primarily governed by van der Waals forces, hydrophobic interactions and hydrogen bonding among the suitable, complementary functional groups of both interacting moieties [15,18,41].

2.2. Influence of Fractional SOM Composition on Carbaryl, Carbofuran and Metolachlor Retention in Studied Soils

According to the previous study on adsorption of the investigated uncharged pesticides in L and C soils amended with various soil colloids, agrochemical retention was attributed to their distribution in soil organic matter and inorganic colloids [10]. Both soils revealed similar mineralogical compositions according to which montmorillonite and kaolinite were major clay minerals. However, substantial variations were revealed in the fractional compositions of organic matter of soils L and C (Table 4). According to them, loamy sand material is rich in humic acid fractions (CHAs), which constitutes 40.07% of total organic carbon (TOC), whereas in loam soil residual carbon fractions (CRs, also called stable organic matter or humin (HN) fraction) represent 59.71% with a dominant share in TOC. This discrepancy in the mutual share of particular SOM fractions and their chemodiversity will have a direct impact on the magnitude and stability of interactions with pesticides [16]. Humic substances (HS) may react with uncharged pesticides via different modes, of which hydrophobic retention is regarded a common adsorption mechanism [18,64]. This may occur on hydrophobic active sites of HS where the competition for water is minimized [64]. These regions include aliphatic side-chains (common in humin fraction) or lipid constituents, as well as lignin-derived moieties with a high carbon content and a small number of polar groups. Van der Waals interactions also play an important role in the adsorption of non-ionic pesticides on complementary sites of HA. Since these forces are additive, their contribution increases with the size of the interacting molecule and with its capacity to adapt to the HA surface. This SOM fraction presumably consists of micelles of a polymeric nature, where the basic structure consists of an aromatic ring (di- or trihydroxy-phenyl type) bridged by -O-, -NH-, -N=, -S- and other functional groups that contain both free -OH and -COOH groups and the double linkages of quinones [18]. These functionalities may be involved in adsorption of non-ionic polar pesticides directly via hydrogen bonding or indirectly by affecting the hydrophilic/hydrophobic balance of the SOM surface [64]. Therefore, it is well documented that natural SOM fractions (humic and fulvic acids, humin), as well as different biopolymers (e.g., cellulose, lignin, cutin), due to the structural differences, may exhibit various affinities to organic contaminants, which can be also evaluated by their various Koc coefficients [23].
When we compare the calculated Koc values for carbamates in the investigated L and C soils (Table 3), an increase in the soil C affinity to the pesticides can be elucidated. This is particularly evident in the case of carbaryl, where an over 5.7 times higher Koc value for loamy soil is obtained, whereas for carbofuran the increase is by a factor of 1.25. In addition, the maximum monolayer adsorption capacities Qmax of carbamates approximated for C soil were twice as much as for L soil, which correlates well to the doubled (compared to L soil) content of the residual carbon fraction of C soil (Table 4). Since the diverse share of OC fractions of the studied soils will evoke structural changes in SOM, this should be also reflected in their degree of sorption site heterogeneity [65]. According to a dual-mode sorption model by Xing and Pignatello [66], the more “glassy” or condensed nature of the organic matter, such as in the case of its humin (HN) or residual carbon (CR) fraction, the higher the nonlinearity of the sorption isotherm. Furthermore, a lower nf value for carbaryl sorption in C soil, compared to sandy soil, points to a higher degree of sorption site heterogeneity of loam soil components active in insecticide uptake, which in turn may be due to its high CR fraction content. Since the pH of both tested soil materials was similar, the enhancement of C soil’s adsorption observed for the two carbamates can be related to its residual carbon and partly to humic acid carbon contents. This may be due to the pronounced hydrophobic retention of the carbamate by the CR fraction, supported by hydrogen bonding via the carbaryl’s amide carbonyl group and carboxylic groups of the HA fraction [67].
The postulated particularly high affinity of carbaryl to the soil residual carbon fraction is in agreement with the studies of Murthy and Raghu (1991, [68]), who found that nearly 78.9% of carbaryl was accumulated in the humin fraction of sandy loam soil, followed by the fulvic acid fraction (19.2%), whereas the least amount of pesticide was measured in the humic acid fraction (1.9%). Worth mentioning is the pH (6.8) of the investigated soil, which was very similar to that of the L soil studied herein. To further corroborate this finding, according to the literature, among the OC fractions the residual stable lignin fraction (of aromatic backbone, containing 64–65% of carbon) followed by CFA and CHA accounted for greater adsorption of some pesticides, including carbaryl [27]. In another study [41] it was also proven that carbaryl sorption on humic acid as a single sorbent was scarce in comparison to a more aromatic biochar, where a hydrophobic effect together with van der Waals forces were postulated to be responsible for insecticide effective retention, similar to the interactions claimed herein for soil C. Hence, the high CR content of C soil may significantly contribute to the generally higher sorption magnitude of carbaryl compared to sandy loam material.
According to the literature [69], carbofuran predominantly incorporates into the recalcitrant humin fraction of SOM, constituting approximately 50% of the initial insecticide dose introduced into the soil. In a previous study [41], it was proven that retention of non-ionic carbamates on HA was much weaker (35.4% of carbofuran and 10.2% of carbaryl initial dose) than on more hydrophobic wheat-straw biochar. Therefore, hydrophobic attraction may also govern carbofuran uptake by the particular SOM fractions. In addition, studies performed on the same L and C soils investigated herein [10] have also shown that carbofuran retention was mainly controlled by goethite, expandable CMs, followed by HA (a slight but significant sorption increase) on sandy loam soil. This is in accordance with the study of Singh et al. [70], who found that carbofuran adsorption was strongly correlated with clay content followed by the organic matter content of soils. However, insecticide uptake on the loam soil material was attributed solely to montmorillonite, presumably due to the enhanced physical diffusion of the pesticide into Mt interlayers, governed by van der Waals forces [71]. It is commonly observed that the particulate organic (or organo-mineral) soil fractions form stable microaggregates where the retention of pesticides could depend more on steric factors and diffusion mechanisms than on the type and amounts of oxygen-containing functional groups [72]. Therefore, both organic (mainly CR) and inorganic (Mt) soil colloids contribute to carbofuran retention in the studied soils via hydrophobic bonding and van der Waals forces. However, care should be taken when assessing their mutual superiority and share in the insecticide uptake by soils. Nevertheless, the results obtained herein clearly demonstrate that carbofuran molecules were bound weakly to loam and sandy loam soils, confirming the physical nature of insecticide sorption and a great leaching potential.
In the present study, metolachlor uptake by soil C was slightly but significantly higher than that of soil L (Figure 4c). Its affinity for organic carbon (Koc) of soil C (Table 3) decreased nearly 1.5 times in comparison to sandy soil. Taking that into account, it can be concluded that the higher herbicide Kf calculated for C soil (Table 1) indicates that apart from particular SOM fractions, other soil constituents must also contribute to its retention in loamy soil. This is in agreement with the studies of Weber et al. [73], who related metolachlor uptake by soil to its organic matter, humic matter (HM) and clay contents. In general, they postulated that the herbicide sorption by soil was of a physical nature and was based on the following mechanisms: hydrophobic bonding to lipophilic SOM sites, charge–transfer interactions, van der Waals forces and H-bonds on polar surfaces of clay minerals. It is in line with other studies where it was proven that metolachlor exhibits a great affinity for the HA polar functional groups, which renders the formation of H-bonding via the carbonyl site and for hydrophobic bonding sites such as aromatic rings [41], which may even preclude H-bonding [18,74]. The recent results of the herbicide sorption on Chernozem soil as well as its HA and HN fractions [75] corroborate the higher affinity of the herbicide for humic acids. Kozak et al. [76] studied the adsorption of metolachlor by various SOM constituents, including humic substances (HA and FA), humin and nonoxidizable soil organic matter. Based on the estimated partition coefficients normalized to organic carbon, they proved that fulvic and humic acids exhibited the highest affinity for the herbicide among the soil OC constituents. On the contrary, studies of Ding et al. [65] revealed that the humin fraction was superior to humic acids in metolachlor uptake. However, the HN fractions they isolated from soils were, uncommonly, much more aromatic than the HA, and the qualitative differences in their structural characteristics were correlated with the long-term tillage management of the soil. This finding is in contradiction with the latest studies on the nature of humin, which demonstrate that its major components are predominantly aliphatic hydrocarbon functionalities [77]. Nevertheless, it must be kept in mind that the type of extraction or isolation method will strongly influence the chemodiversity of this stable organic matter fraction [16,75,78], as well as the results of the studies conducted on the single-sorbent systems. Among the inorganic soil constituents, particularly clay minerals, bentonites and montmorillonites were found to be effective sorbents of metolachlor [79]. According to the literature, the herbicide sorption takes place dominantly as molecular species on hydrophobic microsites of Mt surfaces [8,80]. Presumably, metolachlor is adsorbed at the uncharged and hydrophobic siloxane groups of the clays silicate surface and held there by accepting hydrogen bonds from water molecules around the interlayer cations. Hence, in the present study, it can be postulated that in the case of C soil, the overall sorption magnitude of herbicide can be attributed in part to the humic acid fraction content (decreased by a factor of two in relation to L soil) and increased clay minerals content (mainly Mt [10]).

2.3. pH-Dependent Sorption of Carbamates and Metolachlor on L and C Soils

Sorption of carbaryl by soil L was noted to increase along with the soil pH, appreciably up to pH seven (Figure 6a). According to some studies, carbaryl sorption is due, in part, to site-specific interactions between carbamate functional group and exchangeable cations [27,81], e.g., via Ca2+ bridging. Thus, the trend observed for sandy loam material may be explained by the increasing competition between hydrogen and Ca2+ ions for active sites on the adsorbent surface with pH decline. In such case, available adsorption domains for carbaryl on clay minerals (kaolinite and montmorillonite) and SOM fractions (Table 4) on soil L will decrease. Furthermore, it was found that an increase in soil pH decreases the amount of extractable carbaryl residues [68], increasing bound residues at the same time, which can be retained by the CR fraction of soil OC. The sorption magnitude found for carbaryl on soil C is in a narrow range of 34.19% to 49.55% and seemingly declines with a pH increase, but the differences are not statistically significant (Figure 6a). The higher sorption magnitude observed for carbaryl on C soil of a pH of three to four compared to L soil is presumably related to its higher organic carbon content, as well as SOM composition (Table 4). Over the range of experimental pH values, some part of surface functional groups (e.g., carbonyl, hydroxyl, carboxyl) should become more protonated with a pH decrease. This results in less negatively charged SOM fractions abundant in hydrophilic moieties at a lower pH [82], which together with the high share of CR (59.71% of TOC) may favor adsorption of the 1-naphthyl-N-methyl moiety of carbaryl to the active sorption sites via hydrophobic bonding [28].
Meanwhile, desorption of the agrochemical from the two soils studied was the lowest at pH three; however, magnitude of the phenomena varied with the soil (Table 5). In loam, material carbaryl desorbed up to 25.34% (at pH 7) of the dose initially adsorbed, whereas in loamy sand soil its release ranged from 47.23% to 59.08%. Two conclusions can be drawn from these observations. Firstly, low pH reduces carbaryl desorption from loamy sand and loam soils studied, enhancing its retention. Secondly, carbaryl is more efficiently retained by loamy material, suggesting that the residual carbon fraction may be essential in the stabilization and immobilization of agrochemical residues in the soil environment.
Carbofuran sorption exhibited an inverse trend with a pH increase for both studied soils (Figure 6b). Its sorption magnitude was significantly higher for loam soil material, achieving values of 27.36% and 32.69% of the carbofuran initial dose introduced for soils L and C, respectively. This discrepancy in the insecticide uptake by soils may be related to their different SOM quantities and compositions (Table 4), which significantly influence the number of specific potential sorption centers for the adsorbate. However, these results are in contradiction to the work of Rama Krishna and Philip [35], who found that with a pH increase from two to eight there was an increase in the sorption of carbofuran, methyl parathion and lindane in sandy, red, clayey and compost soils. A similar observation was found by Hsieh and Kao [12], who reported the enhanced adsorption of carbofuran in lateritic soils with a pH increase. Nevertheless, lateritic soils, rich in aluminum and iron minerals, are formed in wet, hot tropical areas, contrary to the soil materials of temperate climate and low inorganic colloids content investigated herein. Another study by Gupta et al. [82] revealed the inverse to pH adsorption trend for carbofuran, corroborating the results presented in this paper. Similarly, Salman et al. reported a favorable adsorption of carbofuran in acidic soils. Meanwhile, desorption of the insecticide was higher in sandy loam soil material (up to 98.1% of carbofuran initial dose at pH seven, Table 5), suggesting that the forces holding it there were weak, enhancing a serious threat of groundwater contamination. Carbofuran release in soil C was significantly lower (maximum of 33.05%), suggesting the high efficiency of soil colloids in its immobilization in soil.
Metolachlor sorption in soil L is also pH dependent, and similarly to carbofuran, it declines inversely with pH. The same behavior of pesticide can be observed in loamy soil (Figure 6c), which reveals an overall higher sorption capacity under the experimental conditions than the loamy sand soil material. This can be explained by the fact that a lower pH facilitates hydrogen bonds between the carbonyl oxygen of metolachlor and H atoms of carboxyl and hydroxyl groups of the organic moieties and charge–transfer (π-π) bonds between the aromatic nucleus of the herbicide and the aromatic rings on the organic matter surfaces [83]. Meanwhile, the desorption of agrochemicals (Table 5) increases with pH in sandy loam soil and is hindered by lower pH (three to five) conditions in case of loam soil material. Nevertheless, its magnitude is much higher in L soil (61.67%), whereas it seems to be limited by some C soil constituents (maximum desorption of 18.84% at pH seven), as they retain agrochemicals more efficiently than loamy sand.

3. Materials and Methods

3.1. Chemicals

Pesticides chosen for this study were carbaryl (1-naphthyl methylcarbamate), carbofuran (2,2-Dimethyl-2,3-dihydro-1-benzofuran-7-yl methylcarbamate) and metolachlor ((RS)-2-Chloro-N-(2-ethyl-6-methyl-phenyl)-N-(1-methoxypropan-2-yl)acetamide). All the active agrochemicals substances were purchased from Merck (Darmstadt, Germany) with a reported purity of >98%. In the case of metolachlor, its stock standard solution in acetonitrile (100 mg L−1) was used (HPLC grade). All reagents were of analytical grade. Working standard mixed solutions were prepared weakly by diluting each individual stock solution in 10 mM CaCl2 prepared in Milli-Q water and stored at 4 °C. The pesticides studied and some of their physical–chemical properties are presented in Table 6.

3.2. Arable Soils under Study

The representative 2 topsoil samples (0–30 cm) from farmlands in Ligota Piękna (near Wrocław, Poland) were randomly collected for the study (at 12 points within each location). They were selected due to the agricultural land use and texture, as sandy loam (L soil material) represents 13.3% and loam (C soil) represents 16.8% of the topsoil in south-west Poland [84]. They are the 2 most abundant groups among light- and heavy-textured soils of the region. According to the WRB classification [85], the soils were mostly Phaeozems and Umbrisols. Prior to the analyses, soil materials were air-dried, crushed gently and sieved through a sieve with a mesh size of 2 mm. The fraction of the soils <2 mm was collected and used in the study. Table 7 shows a summary of the physicochemical properties of the investigated soil materials and was derived from a previous study [10].
In the fraction <2 µm of both soils, montmorillonite was the major clay mineral, accompanied by kaolinite and small amounts of highly dispersed quartz. In addition, both soils contained low level of carbonates at 0.6% and 1.1% in L and C soil, respectively. Ca2+ was the main exchangeable cation in both tested soils.

3.2.1. Chemical Fractionation of the Humic Material

In order to correlate the soils affinity for the studied pesticides to a particular soil organic matter fraction, the extraction and chemical fractionation of the humic substances from investigated L and C soils was conducted. It was undertaken according to the method recommended by the International Humic Substances Society (IHSS), described in Swift [86], with some modification [87,88,89,90]. The first step in the procedure was sample decalcitation with H2SO4 [91] in order to obtain a low-molecular-weight fulvic fraction (CAC). In the next step, soil was exhaustively extracted with 0.1M NaOH until the supernatant was completely discolored. As a result, a fraction called total extractable carbon (TEC) was obtained. It consisted of humic acids (HA) and fulvic acids (FA). In the subsequent step, the alkaline extract was acidified to pH 2 to precipitate the humic acid fraction (CHA), which was then centrifuged from the CFA fraction remaining in the solution.
The carbon of humic acid (CHA) content was calculated from the difference between total extractable carbon (TEC) measured in alkali (NaOH) extract and carbon of fulvic acid (CFA) fraction after HA precipitation at a pH of less than 2. Carbon of nonhydrolyzable humic substances (also called residual carbon or humin fraction, CR) was calculated from the difference between total organic carbon (TOC) and sum of humic and fulvic acid fractions (CR = TOC − (CHA + CFA)).
Total organic carbon and concentrations of organic carbon in the subsequent fractions were measured after each extraction step on an Enviro TOC/TN analyzer (Elementar, Langenselbold, Germany).

3.2.2. Preparation of L and C Soils of Various pH Levels

To set a series of soil L and C at various pHs in the range of 3–7, the following procedure was implemented. An aliquot (200 g) of each soil (L or C) was placed in a beaker, soaked with 50 mL of deionized water and stirred for 12 h. Small volumes of 1 M HCl solution were added dropwise to a soil slurry and stirred to obtain the desired final pH values. The samples were equilibrated for another 12h followed by the pH measurement. The procedure was repeated until each soil series contained samples of pH adjusted to approximately: 3, 4, 5 and 7. Obtained soils were air-dried, ground and subjected to the final pH test. Each soil sample (n = 3) was equilibrated with 10 mM CaCl2 solution in a ratio of 1:2 (used in the batch studies) to ascertain that the experimental pH values would be maintained.
The pH range was set to cover the values typical for undisturbed soils, excluding values higher than 7. This was due to the fact that investigated carbamates undergo hydrolysis in alkaline conditions [38,42] and, furthermore, in pH values exceeding 7, humic acid fraction occurs in its dissolved form. Therefore, to obtain reliable results concerning carbamate concentration in a solution after sorption and to avoid the loss of potential sorbent from the solid phase, the upper limit for the pH range investigated was set to 7.

3.3. Batch Equilibrium Sorption Studies

OECD 106 and USEPA guidelines for the testing of chemicals were adopted for the sorption-desorption studies [63]. They were conducted separately for the sandy loam (L) and loam (C) soils with each of the 3 pesticides tested. Performing batch studies on whole soils allowed for the preservation of complex chemical compositions and mutual interactions among constituents of soil organic matter and inorganic colloidal fractions, as well as their natural chemodiversity. The preliminary investigation showed a soil to solution (10 mM CaCl2) ratio of 1:2 was the most appropriate for all samples [10]. Therefore, subsequent experiments were carried out using this ratio.
Triplicate soil aliquots (5 g of air-dried soil L or C) were equilibrated with 10 mL of each pesticide aqueous solution by gently rotating in an end-over shaker at 20 ± 2 °C for 24 h [10]. The concentration ranges for the batch equilibrium studies were 0–100 mg L−1 for carbofuran, 0–80 mg L−1 for carbaryl and 0–50 mg L−1 for metolachlor on both studied soils. After equilibration, the soil suspensions were centrifuged for 15 min at 10,000 rpm. The concentration of the tested pesticide in the supernatant (Ce; mg L−1) was measured, and the amount of pesticide sorbed on the soils (Qs; mg kg−1) was calculated from the concentration difference between the initial pesticide solution (C0) and Ce:
Qs = [(C0 − Ce)V]/m,
where V is the volume of solution (L) and m is the mass of soil (kg).
To evaluate the influence of pH on the selected carbamates and metolachlor adsorption in L and C soils, a single-concentration sorption-desorption experiment [63] was performed. The concentrations of the tested pesticides were as follows: carbaryl—10 mg L−1 and 20 mg L−1 (for L and C soils, respectively); carbofuran—15 mg L−1 and metolachlor 10 mg L−1 for both studied soils. Single point concentrations were chosen based on the batch studies, as they represent the pesticide initial concentration for which the sorption magnitude was maintained above 40% [63]. The soil to solution ratio, equilibration time, measuring conditions and methods used to calculate the sorption magnitude were identical to those in the batch experiment described above. Desorption studies were performed straight after taking up the supernatant after 24 h of equilibration. To the same vessel with the soil, an aliquot of background (10 mM CaCl2) solution was added and stirred for another 24 h. After centrifugation, the supernatants were filtered (0.45 μm) and subjected to quantitative analysis with LC/MS-MS (Thermo Scientific TSQ Quantum Access MAX) [41]. Briefly, 10 μL of the sample was injected into the analytical column (Thermo Scientific Hypersil GOLD column, 50 mm × 2.1 mm, with 1.9 μm particle size; flow rate of 0.50 mL min−1). Deionized water with 0.1% of formic acid constituted the mobile phase A, and 90% acetonitrile, 10% deionized water and 0.1% of formic acid was the mobile phase B. Gradient elution was applied as follows: 100% of eluent A from the first minute, which was replaced linearly up to 100% by eluent B in 20 min. This composition was maintained for a further 3 min before 100% of eluent A was applied, followed by a re-equilibration time of 7 min to give a total run time of 30 min. The mass spectrometer was run in single reaction monitoring (SRM) mode. Argon (99.995%) with a pressure of 1.87 × 10−3 mbar was applied in the collision cell. A capillary voltage of 4 kV was used in positive ionization mode. The ion transfer capillary temperature was set to 380 °C and the vaporizer of the ESI probe to 400 °C.
Desorption (D) was calculated using the following formula:
D (%) = (Qde/Qad) × 100%,
where Qde is the quantity of pesticide that desorbed (mg kg−1) and Qad is the quantity adsorbed after pesticide addition to the soil [92]. The desorption magnitude was expressed in relation to the pesticide amount initially sorbed during the first 24 h of contact with L or C soil [6].

3.4. Mathematical Modeling of Sorption Isotherms

From the batch studies performed, isotherms of an adsorption process were constructed using the Langmuir, Freundlich, Temkin and Dubinin–Radushkievich (D–R) models. They describe the distribution of adsorbate molecules between the liquid and the solid phase [93], showing a relationship between the amounts of pesticide adsorbed per unit weight of the soil (Qs) and the pesticide concentration in the solution at equilibrium (Ce). The adsorption curves and values of important variables of the linearized forms of isotherm equations for the studied pesticides in the soil samples were also determined.

3.4.1. Langmuir Model

This model determines the saturation of the available adsorption sites and assumes that the adsorbate can form a monolayer of molecules interacting with homogeneous adsorption sites of equivalent energies on the adsorbent surface and that no interactions between adsorbed species occur [1,94]. The retained amount increases until a limited value corresponding to the surface coating by a monolayer is reached.
The expression for the model is given by the following equation:
Qe = (Qmax KLCe/1 + KLCe),
Qe is the amount of adsorbate at equilibrium (mg kg−1); Ce is the equilibrium concentration of the adsorbate (mg L−1); Qmax is the maximum adsorption capacity (mg kg−1); KL is the Langmuir constant (the binding energy associated with adsorbent and adsorbate, L mg−1).
The linearized form of the Langmuir equation derived on the basis of Equation (3) is expressed as:
Ce/Qe = 1/(KLQmax) + Ce/Qmax,
Additionally, the Langmuir parameter (RL), also called a separation factor, was determined for the studied initial concentration range:
RL=1/(1+ KLC0),
Adsorption can be classified as irreversible if RL = 0; favorable if 0 < RL < 1; linear if RL = 1; or unfavorable if RL > 1, informing about the change in soil adsorption strength, the type of isotherm and adsorption state [95,96].

3.4.2. Freundlich Model

The Freundlich isotherm is an empirical model describing multilayer adsorption with non-uniform distribution of adsorption enthalpy and affinities onto the heterogeneous adsorbent surface without lateral interaction [92,97]. The energetically favored binding sites are theorized to be occupied first and the binding strength decreases sequentially with increased coverage of the sites. The model can be described by the following equation:
Qe = Kf Cenf,
where Qe is the amount of adsorbate at equilibrium (mg kg−1); Ce is the equilibrium concentration of the adsorbate (mg L−1) and Kf is the Freundlich affinity coefficient or Freundlich solid-water distribution coefficient (mg (1-n) Ln kg−1). Kf represents the pesticide affinity for soil, with high values indicating a stronger adsorption for the pesticide and also suggesting a lower mobility of the substance in the soil profile [98,99]. nf is the exponential coefficient (constant) associated with the energy distribution of the adsorption site, also being a measure of favorability of sorption (for the values of nf between 0 and 1) [100].
The pesticide concentration in the soil and soil solutions are fitted to the linearized Freundlich sorption isotherm equation:
ln Qe = ln Kf + nf ln Ce,
The Freundlich equation in a linear mode yields a linear plot in which the slope is nf and characterizes the isotherm nonlinearity.
Based on the constants obtained from the Freundlich approximation of the experimental data, a thermodynamic parameter determining the spontaneity of the adsorbate–adsorbent interactions, namely, standard Gibb’s free energy (ΔG°), can be calculated:
ΔG°= −RT ln Kc,
where T is the temperature (K), R is the gas constant (8.314 J mol−1 K−1) and Kc is the standard equilibrium constant (dimensionless), which can be obtained by Kf conversion using the proposed equation [54,55,56]:
K c = K f   ρ 1000   10 6 ρ 1 n f ,
where ρ is the density of pure water (1.0 g/mL). For other experimental equilibrium constants KL, Kf and KD, the strategies for converting them into Kc were adapted [101,102,103,104] and must not be omitted when elaborating ΔG° calculations.
A negative ΔG° value indicates that the process is thermodynamically favorable and can conduct spontaneously. The greater the absolute magnitude of the ΔG° value, the higher is the extent to which the adsorption reaction may take place [35].
Additionally, based on the obtained Kf parameters for each pesticide, the organic-carbon normalized Freundlich distribution coefficients Kfoc (in the text indicated as Koc—organic carbon-water partition coefficient) were calculated according to the equation below:
Kfoc = (Kf/OC) × 100%,
where OC is the organic carbon content of the soil (%).
The organic carbon partition coefficient (Koc) indicates a magnitude of pesticide affinity for the organic fraction of the soil [105]. Its high values suggest that the pesticide is firmly fixed in the soil organic matter fraction, so only a small amount of the compound could undergo runoff to surface waters or leaching to aquifers.
Pesticide leachability (or mobility) in the soil profile and the risk of groundwater contamination can also be calculated by determining the groundwater ubiquity score (GUS) index [59]:
GUS index = (log t1/2) × (4 − log(Koc)),
where t 1/2 is the half-life of pesticide and Koc is the distribution coefficient normalized to the organic carbon content of soil. GUS index values below 1.8 characterize non-leachers.

3.4.3. Temkin Model

The Temkin isotherm assumes a linear rather than logarithmic decrease in heat of adsorption of all molecules in the layer, while ignoring extremely low and very high concentrations of solute [92]. It also assumes uniform distribution of bounding energy up to some maximum. It is expressed by the following equation:
Qe= (RT/bT) ln Ce + (RT/ bT) ln A,
Qe = B ln Ce + B ln A,
where Qe is the amount of chemical adsorbed at equilibrium (mg kg−1) and Ce is the concentration of adsorbate in the solution at equilibrium (mg L−1). B is a constant defined by the expression B = RT/bT, bT is the Temkin constant related to the heat of adsorption (J mol−1), T is the absolute temperature (K), R is the gas constant (8.314 J mol−1 K−1) and A is the Temkin isotherm binding constant (L kg−1). From the plot of Qe vs. ln Ce (Equation (12)), B and A can be calculated from the slopes (B) and intercepts (B lnA), respectively.

3.4.4. Dubinin–Radushkievich Model

The equilibrium data can also be applied to the Dubinin–Radushkievich (D–R) model to determine whether the sorption is of physical or chemical type. It is an empirical adsorption model that is generally used to express adsorption mechanisms with Gaussian energy distribution onto heterogeneous surfaces [106,107]. The linear form of the D–R isotherm is presented in the form of the following equation:
lnQe = ln(Qm) − βε2,
where Qe is the amount of pesticide adsorbed into the soil (mg kg−1 ), Qm is the theoretical monolayer sorption capacity (mg kg−1 ), β is the constant (mol2 J −2), which is related to the average energy of sorption per mole of the sorbate as it is transferred to the surface of the solid from an infinite distance in the solution, and ε is the Polanyi potential, which is described by Equation (15), where T is the solution temperature (K) and R is the gas constant.
ε = RT ln(1 + 1/Ce),
The value of mean sorption energy, Ea (kJ mol−1), can be calculated from D–R parameter β, (obtained from the relationship of ln Qe and ε2, Equation (14)), according to the equation:
Ea= 1/(2 β)1/2
The value of the mean sorption energy provides information about chemical and physical sorption. The Ea value ranges from approximately 1 kJ mol−1 to 8 kJ mol−1 for physical sorption and from 8 kJ mol−1 to 16 kJ mol−1 for chemical sorption [53].

3.5. Statistical Analysis

Data are expressed as the mean values of 3 replicates ± standard deviation (SD). The values of pesticide sorption and desorption magnitudes on L and C soils, including pH influence (both sorption-desorption data), were statistically compared using Student’s t-test (p < 0.05), prior to which the Shapiro–Wilk normality test was performed. The data were compiled using Microsoft Excel 2019 and GraphPad Prism 9.5.1.733 software.

4. Conclusions

Carbaryl, carbofuran and metolachlor sorption were better retained by loam than sandy loam topsoil material. Langmuir, Freundlich and Temkin models fitted well to the obtained experimental L-type isotherms. Approximated parameters point to favorable and spontaneous multilayer adsorption with a higher degree of surface coverage with the pesticides on C soil compared to soil L. The Dubinin–Radushkievich approach confirmed the physical nature of the studied agrochemical adsorption. GUS indexes indicate low leaching potential of metolachlor and moderate potential of carbofuran on both soils. Carbaryl is a non-leacher in loam, while on sandy soil it reveals a moderate percolation ability.
Sandy loam material was more abundant in humic acid (CHA) fractions, whereas loam soil was abundant in residual carbon (CR). These differences in SOM composition were postulated to strongly influence the studied pesticide retention in soils. The residual carbon fraction of SOM is presumably superior in immobilization of carbaryl in the soil environment. Carbofuran retention in both soils was related to clay minerals (montmorillonite) as well as CR fraction content. The joint effect of humic acids and Mt was responsible for metolachlor uptake in both soils. Lower pH enhanced non-ionic pesticide retention in studied soils, except for carbaryl sorption in sandy soil. Therefore, SOM fractional composition and the mutual share of individual OC fractions together with pH are crucial for predicting the behavior of non-ionic pesticides in the soil environment.

Author Contributions

Conceptualization, methodology, software, validation, formal analysis, investigation, resources, data curation, writing—original draft preparation, writing—review and editing, visualization, project administration, funding acquisition: I.Ć.-P. All authors have read and agreed to the published version of the manuscript.

Funding

The research was funded by the National Science Center, Poland, grant number DEC-2012/05 /D/ST10/02223, and by Wroclaw University of Environmental and Life Sciences from the statutory funds for research. The APC is co-financed by Wroclaw University of Environmental and Life Sciences.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data will be made available upon request.

Acknowledgments

I would like to express special thanks to Bogna Kosyk and Anna Fleszar for their assistance in pesticide LC-MS/MS measurements.

Conflicts of Interest

The author declares no conflict of interest.

Sample Availability

Not applicable.

References

  1. Copaja, S.V.; Gatica, P. Effects of clay content in soil on pesticides sorption process. J. Chil. Chem. Soc. 2021, 66, 5086–5092. [Google Scholar] [CrossRef]
  2. Kah, M.; Brown, C.D. Adsorption of Ionisable Pesticides in Soils. In Reviews of Environmental Contamination and Toxicology: Continuation of Residue Reviews; Ware, G.W., Whitacre, D.M., Albert, L.A., de Voogt, P., Gerba, C.P., Hutzinger, O., Knaak, J.B., Mayer, F.L., Morgan, D.P., Park, D.L., et al., Eds.; Springer: New York, NY, USA, 2006; pp. 149–217. ISBN 978-0-387-32964-2. [Google Scholar]
  3. Rasool, S.; Rasool, T.; Gani, K.M. A review of interactions of pesticides within various interfaces of intrinsic and organic residue amended soil environment. Chem. Eng. J. Adv. 2022, 11, 100301. [Google Scholar] [CrossRef]
  4. Spark, K.M.; Swift, R.S. Effect of soil composition and dissolved organic matter on pesticide sorption. Sci. Total Environ. 2002, 298, 147–161. [Google Scholar] [CrossRef] [PubMed]
  5. Barchańska, H.; Czaplicka, M.; Kyzioł-Komosińska, J. Interaction of selected pesticides with mineral and organic soil components. Arch. Environ. Prot. 2020, 46, 80–91. [Google Scholar] [CrossRef]
  6. Ćwieląg-Piasecka, I.; Jamroz, E.; Medyńska-Juraszek, A.; Bednik, M.; Kosyk, B.; Polláková, N. Deashed Wheat-Straw Biochar as a Potential Superabsorbent for Pesticides. Materials 2023, 16, 2185. [Google Scholar] [CrossRef]
  7. Haberhauer, G.; Pfeiffer, L.; Gerzabek, M.H.; Kirchmann, H.; Aquino, A.J.A.; Tunega, D.; Lischka, H. Response of sorption processes of MCPA to the amount and origin of organic matter in a long-term field experiment. Eur. J. Soil Sci. 2001, 52, 279–286. [Google Scholar] [CrossRef] [Green Version]
  8. Cox, L.; Celis, R.; Hermosín, M.C.; Cornejo, J.; Zsolnay, A.; Zeller, K. Effect of Organic Amendments on Herbicide Sorption as Related to the Nature of the Dissolved Organic Matter. Environ. Sci. Technol. 2000, 34, 4600–4605. [Google Scholar] [CrossRef] [Green Version]
  9. Chen, H.; He, X.; Rong, X.; Chen, W.; Cai, P.; Liang, W.; Li, S.; Huang, Q. Adsorption and biodegradation of carbaryl on montmorillonite, kaolinite and goethite. Appl. Clay Sci. 2009, 46, 102–108. [Google Scholar] [CrossRef]
  10. Ćwieląg-Piasecka, I.; Debicka, M.; Medyńska-Juraszek, A. Effectiveness of Carbaryl, Carbofuran and Metolachlor Retention in Soils under the Influence of Different Colloid. Minerals 2021, 11, 924. [Google Scholar] [CrossRef]
  11. Ozbay, B.; Akyol, N.H.; Akyol, G.; Ozbay, I. Sorption and desorption behaviours of 2,4-D and glyphosate in calcareous soil from Antalya, Turkey. Water Environ. J. 2018, 32, 141–148. [Google Scholar] [CrossRef]
  12. Hsieh, T.-L.; Kao, M.-M. Adsorption of carbofuran on lateritic soils. J. Hazard. Mater. 1998, 58, 275–284. [Google Scholar] [CrossRef]
  13. Wu, D.; Ren, C.; Wu, C.; Li, Y.; Deng, X.; Li, Q. Mechanisms by which different polar fractions of dissolved organic matter affect sorption of the herbicide MCPA in ferralsol. J. Hazard. Mater. 2021, 416, 125774. [Google Scholar] [CrossRef] [PubMed]
  14. Sadegh-Zadeh, F.; Abd Wahid, S.; Jalili, B. Sorption, degradation and leaching of pesticides in soils amended with organic matter: A review. Adv. Environ. Technol. 2017, 3, 119–132. [Google Scholar] [CrossRef]
  15. Li, Y.-L.; He, W.; Liu, W.-X.; Kong, X.-Z.; Yang, B.; Yang, C.; Xu, F.-L. Influences of binding to dissolved organic matter on hydrophobic organic compounds in a multi-contaminant system: Coefficients, mechanisms and ecological risks. Environ. Pollut. 2015, 206, 461–468. [Google Scholar] [CrossRef]
  16. Ukalska-Jaruga, A.; Bejger, R.; Smreczak, B.; Podlasiński, M. Sorption of Organic Contaminants by Stable Organic Matter Fraction in Soil. Molecules 2023, 28, 429. [Google Scholar] [CrossRef]
  17. Vischetti, C.; Monaci, E.; Casucci, C.; De Bernardi, A.; Cardinali, A. Adsorption and Degradation of Three Pesticides in a Vineyard Soil and in an Organic Biomix. Environments 2020, 7, 113. [Google Scholar] [CrossRef]
  18. Senesi, N. Binding mechanisms of pesticides to soil humic substances. Behav. Pestic. Soil Environ. 1992, 123–124, 63–76. [Google Scholar] [CrossRef]
  19. Murphy, E.M.; Zachara, J.M.; Smith, S.C.; Phillips, J.L.; Wietsma, T.W. Interaction of Hydrophobic Organic Compounds with Mineral-Bound Humic Substances. Environ. Sci. Technol. 1994, 28, 1291–1299. [Google Scholar] [CrossRef] [Green Version]
  20. Becher, M.; Tołoczko, W.; Godlewska, A.; Pakuła, K.; Żukowski, E. Fractional Composition of Organic Matter and Properties of Humic Acids in the Soils of Drained Bogs of the Siedlce Heights in Eastern Poland. J. Ecol. Eng. 2022, 23, 208–222. [Google Scholar] [CrossRef]
  21. Pulido-Moncada, M.; Lozano, Z.; Delgado, M.; Dumon, M.; Van Ranst, E.; Lobo, D.; Gabriels, D.; Cornelis, W.M. Using soil organic matter fractions as indicators of soil physical quality. Soil Use Manag. 2018, 34, 187–196. [Google Scholar] [CrossRef] [Green Version]
  22. Song, G.; Simpson, A.J.; Hayes, M.H.B. Compositional changes in the humin fraction resulting from the long-term cultivation of an Irish grassland soil: Evidence from FTIR and multi-NMR spectroscopies. Sci. Total Environ. 2023, 880, 163280. [Google Scholar] [CrossRef] [PubMed]
  23. Novotny, E.H.; Turetta, A.P.D.; Resende, M.F.; Rebello, C.M. The quality of soil organic matter, accessed by 13C solid state nuclear magnetic resonance, is just as important as its content concerning pesticide sorption. Environ. Pollut. 2020, 266, 115298. [Google Scholar] [CrossRef] [PubMed]
  24. Ahangar, A.G.; Smernik, R.J.; Kookana, R.S.; Chittleborough, D.J. Clear effects of soil organic matter chemistry, as determined by NMR spectroscopy, on the sorption of diuron. Chemosphere 2008, 70, 1153–1160. [Google Scholar] [CrossRef] [PubMed]
  25. Chefetz, B.; Xing, B. Relative Role of Aliphatic and Aromatic Moieties as Sorption Domains for Organic Compounds: A Review. Environ. Sci. Technol. 2009, 43, 1680–1688. [Google Scholar] [CrossRef]
  26. Mitchell, P.J.; Simpson, M.J. High Affinity Sorption Domains in Soil Are Blocked by Polar Soil Organic Matter Components. Environ. Sci. Technol. 2013, 47, 412–419. [Google Scholar] [CrossRef] [PubMed]
  27. Raturi, S.; Islam, K.R.; Caroll, M.J.; Hill, R.L. Carbaryl, 2,4-D, and Triclopyr Adsorption in Thatch-Soil Ecosystems. J. Environ. Sci. Health Part B 2005, 40, 697–710. [Google Scholar] [CrossRef]
  28. Singh, R.P.; Singh, S.; Srivastava, G. Adsorption Thermodynamics of Carbaryl onto Four Texturally Different Indian Soils. Adsorpt. Sci. Technol. 2011, 29, 277–288. [Google Scholar] [CrossRef]
  29. Singh, R.P.; Srivastava, G. Adsorption and Movement of Carbofuran in Four Different Soils Varying in Physical and Chemical Properties. Adsorpt. Sci. Technol. 2009, 27, 193–203. [Google Scholar] [CrossRef]
  30. Ukalska-Jaruga, A.; Smreczak, B.; Siebielec, G. Assessment of Pesticide Residue Content in Polish Agricultural Soils. Molecules 2020, 25, 587. [Google Scholar] [CrossRef] [Green Version]
  31. Kouame, K.B.-J.; Savin, M.C.; Willett, C.D.; Bertucci, M.B.; Butts, T.R.; Grantz, E.; Roma-Burgos, N. S-metolachlor persistence in soil as influenced by within-season and inter-annual herbicide use. Environ. Adv. 2022, 9, 100318. [Google Scholar] [CrossRef]
  32. Singh, N.; Kloeppel, H.; Klein, W. Sorption behavior of metolachlor, isoproturon, and terbuthylazine in soils. J. Environ. Sci. Health Part B 2001, 36, 397–407. [Google Scholar] [CrossRef] [PubMed]
  33. Mathias, F.T.; Romano, R.M.; Sleiman, H.K.; de Oliveira, C.A.; Romano, M.A. Herbicide Metolachlor Causes Changes in Reproductive Endocrinology of Male Wistar Rats. ISRN Toxicol. 2012, 2012, 130846. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Obrigawitch, T.; Hons, F.M.; Abernathy, J.R.; Gipson, J.R. Adsorption, Desorption, and Mobility of Metolachlor in Soils. Weed Sci. 1981, 29, 332–336. [Google Scholar] [CrossRef]
  35. Rama Krishna, K.; Philip, L. Adsorption and desorption characteristics of lindane, carbofuran and methyl parathion on various Indian soils. J. Hazard. Mater. 2008, 160, 559–567. [Google Scholar] [CrossRef]
  36. Spongberg, A.L.; Ganliang, L. Adsorption of atrazine and metolachlor in three soils from Blue Creek wetlands, Waterville, Ohio. Sci. Soils 2000, 5, 1–9. [Google Scholar] [CrossRef]
  37. Mondal, N.K.; Chattoraj, S.; Sadhukhan, B.; Das, B. Evaluation of carbaryl sorption in alluvial soil. Songklanakarin J. Sci. Technol. 2013, 35, 727–738. [Google Scholar]
  38. Zhang, P.; Sun, H.; Yu, L.; Sun, T. Adsorption and catalytic hydrolysis of carbaryl and atrazine on pig manure-derived biochars: Impact of structural properties of biochars. J. Hazard. Mater. 2013, 244–245, 217–224. [Google Scholar] [CrossRef]
  39. Getzin, L.W. Persistence and Degradation of Carbofuran in Soil1. Environ. Entomol. 1973, 2, 461–468. [Google Scholar] [CrossRef]
  40. Trotter, D.M.; Kent, R.A.; Wong, M.P. Aquatic fate and effect of Carbofuran. Crit. Rev. Environ. Control 1991, 21, 137–176. [Google Scholar] [CrossRef]
  41. Ćwieląg-Piasecka, I.; Medyńska-Juraszek, A.; Jerzykiewicz, M.; Dębicka, M.; Bekier, J.; Jamroz, E.; Kawałko, D. Humic acid and biochar as specific sorbents of pesticides. J. Soils Sediments 2018, 18, 2692–2702. [Google Scholar] [CrossRef] [Green Version]
  42. Brahmaprakash, G.P.; Sethunathan, N. Metabolism of carbaryl and carbofuran in soil planted to rice. Agric. Ecosyst. Environ. 1985, 13, 33–42. [Google Scholar] [CrossRef]
  43. Farrell, J.; Reinhard, M. Desorption of halogenated organics from model solids, sediments, and soil under unsaturated conditions. 1. Isotherms. Environ. Sci. Technol. 1994, 28, 53–62. [Google Scholar] [CrossRef] [PubMed]
  44. Chiou, C. Fundamentals of the Adsorption Theory. In Partition and Adsorption of Organic Contaminants in Environmental Systems; Wiley Online Library: New Jersey, NJ, USA, 2003; pp. 39–52. ISBN 978-0-471-23325-1. [Google Scholar]
  45. Hammond, K.D.; Conner, W.C. Chapter One—Analysis of Catalyst Surface Structure by Physical Sorption. In Advances in Catalysis; Gates, B.C., Jentoft, F.C., Eds.; Academic Press: Cambridge, MA, USA, 2013; Volume 56, pp. 1–101. ISBN 0360-0564. [Google Scholar]
  46. Nakhli, A.; Khalfaoui, M.; Aguir, C.; Manel, B.; M’henni, M.; Ben Lamine, A. Statistical Physics Studies of Multilayer Adsorption on Solid Surface: Adsorption of Basic Blue 41 Dye onto Functionalized Posidonia Biomass. Sep. Sci. Technol. 2014, 49, 2525–2533. [Google Scholar] [CrossRef]
  47. Koskinen, W.C.; Harper, S.S. The Retention Process: Mechanisms. In Pesticides in the Soil Environment: Processes, Impacts and Modeling; SSSA Book Series; Soil Science Society of America: Madison WI, USA, 1990; pp. 51–77. ISBN 978-0-89118-861-2. [Google Scholar]
  48. Lyubchik, S.; Lyubchik, A.; Lygina, O.; Lyubchik, S.; Fonseca, I. Comparison of the Thermodynamic Parameters Estimation for the Adsorption Process of the Metals from Liquid Phase on Activated Carbons. In Thermodynamics; Moreno-Pirajan, J.C., Ed.; IntechOpen: Rijeka, Croatia, 2011; Chapter 4. [Google Scholar]
  49. Zanella, O.; Tessaro, I.C.; Féris, L.A. Study of CaCl2 as an agent that modifies the surface of activated carbon used in sorption/treatment cycles for nitrate removal. Braz. J. Chem. Eng. 2014, 31. [Google Scholar] [CrossRef] [Green Version]
  50. Foo, K.Y.; Hameed, B.H. Insights into the modeling of adsorption isotherm systems. Chem. Eng. J. 2010, 156, 2–10. [Google Scholar] [CrossRef]
  51. Meftaul, I.M.; Venkateswarlu, K.; Dharmarajan, R.; Annamalai, P.; Megharaj, M. Movement and Fate of 2,4-D in Urban Soils: A Potential Environmental Health Concern. ACS Omega 2020, 5, 13287–13295. [Google Scholar] [CrossRef] [PubMed]
  52. Franus, M.; Bandura, L.; Madej, J. Mono and Poly-Cationic Adsorption of Heavy Metals Using Natural Glauconite. Minerals 2019, 9, 470. [Google Scholar] [CrossRef] [Green Version]
  53. Ugrina, M.; Jurić, A.; Nuić, I.; Trgo, M. Modeling, Simulation, Optimization, and Experimental Verification of Mercury Removal onto Natural and Sulfur-Impregnated Zeolite Clinoptilolite—Assessment of Feasibility for Remediation of Mercury-Contaminated Soil. Processes 2023, 11, 606. [Google Scholar] [CrossRef]
  54. Olawale, S.A.; Okafor, C. Comparing the Different Methods Used in the Determination of Thermodynamic Parameters Using Adsorption of Pb(II) on to Chicken Feather as an Example. J. Mater. Sci. Res. Rev. 2020, 6, 1–11. [Google Scholar]
  55. Masini, J.C.; Abate, G. Guidelines to Study the Adsorption of Pesticides onto Clay Minerals Aiming at a Straightforward Evaluation of Their Removal Performance. Minerals 2021, 11, 1282. [Google Scholar] [CrossRef]
  56. Tran, H.N.; You, S.-J.; Hosseini-Bandegharaei, A.; Chao, H.-P. Mistakes and inconsistencies regarding adsorption of contaminants from aqueous solutions: A critical review. Water Res. 2017, 120, 88–116. [Google Scholar] [CrossRef] [PubMed]
  57. Marbawi, H.; Othman, A.R.; Halmi, M.I.E.; Gansau, J.A.; Sabullah, K.; Yasid, N.A. Re-evaluation of Thermodynamic Parameters using Dimensionless Langmuir Constants from Published Data on Glyphosate Adsorption onto Activated Carbon Loaded with Manganese and Iron. J. Environ. Bioremed. Toxicol. 2020, 3, 6–10. [Google Scholar] [CrossRef]
  58. Vaikosen, E.N.; Davidson, C.M.; Olu-Owolabi, B.I.; Gibson, L.T.; Agunbiade, F.O.; Kashimawo, A.J.; Adebowale, K.O. Kinetic and isotherm studies on the adsorption–desorption of technical-grade endosulfan in loamy soils under Theobroma cacao L cultivation, Southwestern Nigeria. Environ. Sci. Adv. 2023, 2, 257–277. [Google Scholar] [CrossRef]
  59. Gustafson, D.I. Groundwater ubiquity score: A simple method for assessing pesticide leachability. Environ. Toxicol. Chem. 1989, 8, 339–357. [Google Scholar] [CrossRef]
  60. Shan, Z.J.; Zhu, Z.L.; Hua, X.M. Investigation of effect of herbicide-suola on groundwater. Acta Sci. Circumstantine 1994, 14, 72–78. [Google Scholar]
  61. Assessing Soil Contamination A Reference Manual; FAO Pesticide Disposal Series 8; Editorial Group FAO Information Division: Rome, Italy, 2000.
  62. Augustijn-Beckers, P.W.M.; Hornsby, A.G.; Wauchope, R.D. The SCS/ARS/CES Pesticide Properties Database for Environmental Decision-Making. II. Additional Compounds. In Reviews of Environmental Contamination and Toxicology: Continuation of Residue Reviews; Ware, G.W., Ed.; Springer: New York, NY, USA, 1994; pp. 1–82. ISBN 978-1-4612-2662-8. [Google Scholar]
  63. Test, No. 106: Adsorption—Desorption Using a Batch Equilibrium Method; OECD Guidelines for the Testing of Chemicals, Section 1; OECD: Paris, France, 2000; ISBN 978-92-64-06960-2.
  64. Leenheer, J.A.; Ahlrichs, J.L. A Kinetic and Equilibrium Study of the Adsorption of Carbaryl and Parathion upon Soil Organic Matter Surfaces. Soil Sci. Soc. Am. J. 1971, 35, 700–705. [Google Scholar] [CrossRef]
  65. Ding, G. Soil Organic Matter and Metolachlor Sorption Characteristics as Affected by Soil Management. Ph. D. Thesis, University of Massachusetts Amherst Community, Amherst, MA, USA, 2002. [Google Scholar]
  66. Xing, B.; Pignatello, J.J. Dual-Mode Sorption of Low-Polarity Compounds in Glassy Poly(Vinyl Chloride) and Soil Organic Matter. Environ. Sci. Technol. 1997, 31, 792–799. [Google Scholar] [CrossRef]
  67. Zhang, Y.; Zhu, D.; Yu, H. Sorption of Aromatic Compounds to Clay Mineral and Model Humic Substance–Clay Complex: Effects of Solute Structure and Exchangeable Cation. J. Environ. Qual. 2008, 37, 817–823. [Google Scholar] [CrossRef]
  68. Murthy, N.B.K.; Raghu, K. Fate of14C-carbaryl in soils as a function of pH. Bull. Environ. Contam. Toxicol. 1991, 46, 374–379. [Google Scholar] [CrossRef]
  69. Benicha, M.; Mrabet, R.; Azmani, A. Characterization of carbofuran bound residues and the effect of ageing on their distribution and bioavailability in the soil of a sugar beet field in north-western Morocco. Eur. J. Environ. Sci. 2016, 6, 57–63. [Google Scholar] [CrossRef] [Green Version]
  70. Singh, R.P.; Kumari, K.; Singh, D. Influence of different factors on the adsorption of carbofuran (2,3-dihydro-2,2-dimethyl-7-benzofuranyl-N-methyl carbamate) on soils. Ecotoxicol. Environ. Saf. 1994, 29, 70–79. [Google Scholar] [CrossRef] [PubMed]
  71. Isak, N.; Xhaxhiu, K. A Review on the Adsorption of Diuron, Carbaryl, and Alachlor Using Natural and Activated Clays. Remediat. J. 2023, 1–15. [Google Scholar] [CrossRef]
  72. Dorado, J.; Almendros, G. Organo-Mineral Interactions Involved in Herbicide Sorption on Soil Amended with Peats of Different Maturity Degree. Agronomy 2021, 11, 869–886. [Google Scholar] [CrossRef]
  73. Weber, J.B.; McKinnon, E.J.; Swain, L.R. Sorption and Mobility of 14C-Labeled Imazaquin and Metolachlor in Four Soils As Influenced by Soil Properties. J. Agric. Food Chem. 2003, 51, 5752–5759. [Google Scholar] [CrossRef] [PubMed]
  74. Jayasundera, S.; Schmidt, W.F.; Hapeman, C.J.; Torrents, A. Influence of the Chemical Environment on Metolachlor Conformations. J. Agric. Food Chem. 1999, 47, 4435–4442. [Google Scholar] [CrossRef] [PubMed]
  75. Weber, J.; Jamroz, E.; Kocowicz, A.; Debicka, M.; Bekier, J.; Ćwieląg-Piasecka, I.; Ukalska-Jaruga, A.; Mielnik, L.; Bejger, R.; Jerzykiewicz, M. Optimized isolation method of humin fraction from mineral soil material. Environ. Geochem. Health 2022, 44, 1289–1298. [Google Scholar] [CrossRef]
  76. Kozak, J.; Weber, J.B.; Sheets, T.J. Adsorption of prometryn and metolachlor by selected soil organic matter fractions. Soil Sci. 1983, 136, 94–101. [Google Scholar] [CrossRef]
  77. Hayes, M.H.B.; Mylotte, R.; Swift, R.S. Chapter Two—Humin: Its Composition and Importance in Soil Organic Matter. In Advances in Agronomy; Sparks, D.L., Ed.; Academic Press: Cambridge, MA, USA, 2017; Volume 143, pp. 47–138. ISBN 0065-2113. [Google Scholar]
  78. Song, G.; Novotny, E.H.; Simpson, A.J.; Clapp, C.E.; Hayes, M.H.B. Sequential exhaustive extraction of a Mollisol soil, and characterizations of humic components, including humin, by solid and solution state NMR. Eur. J. Soil Sci. 2008, 59, 505–516. [Google Scholar] [CrossRef] [Green Version]
  79. Nennemann, A.; Mishael, Y.; Nir, S.; Rubin, B.; Polubesova, T.; Bergaya, F.; van Damme, H.; Lagaly, G. Clay-Based Formulations of Metolachlor with Reduced Leaching. Appl. Clay Sci. 2001, 18, 265–275. [Google Scholar] [CrossRef]
  80. Shattar, S.F.A.; Zakaria, N.A.; Foo, K.Y. Feasibility of montmorillonite-assisted adsorption process for the effective treatment of organo-pesticides. Desalination Water Treat. 2016, 57, 13645–13677. [Google Scholar] [CrossRef]
  81. Peng, X.; Wang, J.; Fan, B.; Luan, Z. Sorption of endrin to montmorillonite and kaolinite clays. J. Hazard. Mater. 2009, 168, 210–214. [Google Scholar] [CrossRef]
  82. Gupta, V.K.; Ali, I.; Suhas; Saini, V.K. Adsorption of 2,4-D and carbofuran pesticides using fertilizer and steel industry wastes. J. Colloid Interface Sci. 2006, 299, 556–563. [Google Scholar] [CrossRef] [PubMed]
  83. Weiping, L.; Pusino, A.; Gessa, C. Contribution of organic matter to metolachlor adsorption on some soils. J. Environ. Sci. 1995, 7, 121–125. [Google Scholar]
  84. Kabała, C.; Bekier, J.; Binczycki, T.; Bogacz, A.; Bojko, O.; Ćwieląg-Piasecka, I.; Debicka, M.; Cuske, M.; Gałka, B.; Gersztyn, L.; et al. Soils of Lower Silesia: Origins, Diversity and Protection; PTG, PTSH: Wrocław, Poland, 2015; ISBN 978-83-934096-4-8. [Google Scholar]
  85. IUSS Working Group WRB. World Reference Base for Soil Resources 2014, Update 2015. International Soil Classification System for Naming Soils and Creating Legends for Soil Maps; FAO: Rome, Italy, 2015. [Google Scholar]
  86. Swift, R.S. Organic matter characterization (chap 35). In Methods of Soil Analysis. Part 3. Chemical Methods; Sparks, D.L., Page, A.L., Helmke, P.A., Loeppert, R.H., Soltanpour, P.N., Tabatabai, M.A., Johnston, C.T., Sumner, M.E., Eds.; Soil Sci. Soc. Am. Book Series: 5; Soil Science Society of America Journal: Madison, WI, USA, 1996; pp. 1018–1020. [Google Scholar]
  87. Bekier, J.; Jamroz, E.; Walenczak-Bekier, K.; Uściła, M. Soil Organic Matter Composition in Urban Soils: A Study of Wrocław Agglomeration, SW Poland. Sustainability 2023, 15, 2277. [Google Scholar] [CrossRef]
  88. Jamroz, E.; Jerzykiewicz, M. Humic fractions as indicators of soil organic matter responses to clear-cutting in mountain and lowland conditions of southwestern Poland. Land Degrad. Dev. 2022, 33, 368–378. [Google Scholar] [CrossRef]
  89. Ukalska-Jaruga, A.; Smreczak, B.; Klimkowicz-Pawlas, A. Soil organic matter composition as a factor affecting the accumulation of polycyclic aromatic hydrocarbons. J. Soils Sediments 2019, 19, 1890–1900. [Google Scholar] [CrossRef] [Green Version]
  90. Mielnik, L.; Hewelke, E.; Weber, J.; Oktaba, L.; Jonczak, J.; Podlasiński, M. Changes in the soil hydrophobicity and structure of humic substances in sandy soil taken out of cultivation. Agric. Ecosyst. Environ. 2021, 319, 107554. [Google Scholar] [CrossRef]
  91. Senesi, N. Composted materials as organic fertilizers. Adv. Humic Subst. Res. 1989, 81–82, 521–542. [Google Scholar] [CrossRef]
  92. Inyinbor, A.A.; Adekola, F.A.; Olatunji, G.A. Kinetics, isotherms and thermodynamic modeling of liquid phase adsorption of Rhodamine B dye onto Raphia hookerie fruit epicarp. Water Resour. Ind. 2016, 15, 14–27. [Google Scholar] [CrossRef] [Green Version]
  93. Kyzas, G.Z.; Matis, K.A. Nanoadsorbents for pollutants removal: A review. J. Mol. Liq. 2015, 203, 159–168. [Google Scholar] [CrossRef]
  94. Piekarski, J.; Ignatowicz, K.; Dąbrowski, T. Application of an Adsorption Process on Selected Materials, Including Waste, as a Barrier to the Pesticide Penetration into the Environment. Materials 2022, 15, 4680. [Google Scholar] [CrossRef] [PubMed]
  95. Ogunmodede, O.T.; Ojo, A.A.; Adewole, E.; Adebayo, O.L. Adsorptive Removal of Anionic Dye from Aqueous Solutions by Algerian Kaolin: Characteristics, Isotherm, Kinetic and Thermodynamic Studies. Iranica J. Energy Environ. 2015, 6, 147–153. [Google Scholar] [CrossRef]
  96. Yu, H.; Liu, Y.; Shu, X.; Fang, H.; Sun, X.; Pan, Y.; Ma, L. Equilibrium, kinetic and thermodynamic studies on the adsorption of atrazine in soils of the water fluctuation zone in the Three-Gorges Reservoir. Environ. Sci. Eur. 2020, 32, 27. [Google Scholar] [CrossRef]
  97. Shikuku, V.O.; Kowenje, C.O.; Kengara, F.O. Errors in Parameters Estimation Using Linearized Adsorption Isotherms: Sulfadimethoxine Adsorption onto Kaolinite Clay. Chem. Sci. Int. J. 2018, 23, 1–6. [Google Scholar] [CrossRef] [Green Version]
  98. Kodešová, R.; Kočárek, M.; Kodeš, V.; Drábek, O.; Kozák, J.; Hejtmánková, K. Pesticide adsorption in relation to soil properties and soil type distribution in regional scale. J. Hazard. Mater. 2011, 186, 540–550. [Google Scholar] [CrossRef] [PubMed]
  99. Wang, X.Q.; Liu, J.; Zhang, N.; Yang, H. Adsorption, mobility, biotic and abiotic metabolism and degradation of pesticide exianliumi in three types of farmland. Chemosphere 2020, 254, 126741. [Google Scholar] [CrossRef] [PubMed]
  100. Chattoraj, S.; Mondal, N.K.; Das, B.; Roy, P.; Sadhukhan, B. Biosorption of carbaryl from aqueous solution onto Pistia stratiotes biomass. Appl. Water Sci. 2014, 4, 79–88. [Google Scholar] [CrossRef] [Green Version]
  101. Chen, T.; Da, T.; Ma, Y. Reasonable calculation of the thermodynamic parameters from adsorption equilibrium constant. J. Mol. Liq. 2021, 322, 114980. [Google Scholar] [CrossRef]
  102. Liu, Y. Is the Free Energy Change of Adsorption Correctly Calculated? J. Chem. Eng. Data 2009, 54, 1981–1985. [Google Scholar] [CrossRef]
  103. Tran, H.N. Improper Estimation of Thermodynamic Parameters in Adsorption Studies with Distribution Coefficient KD (qe/Ce) or Freundlich Constant (KF): Considerations from the Derivation of Dimensionless Thermodynamic Equilibrium Constant and Suggestions. Adsorpt. Sci. Technol. 2022, 2022, 5553212. [Google Scholar] [CrossRef]
  104. Zhou, X.; Zhou, X. The unit problem in the thermodynamic calculation of adsorption using the langmuir equation. Chem. Eng. Commun. 2014, 201, 1459–1467. [Google Scholar] [CrossRef]
  105. Wauchope, R.D.; Yeh, S.; Linders, J.B.H.J.; Kloskowski, R.; Tanaka, K.; Rubin, B.; Katayama, A.; Kördel, W.; Gerstl, Z.; Lane, M.; et al. Pesticide soil sorption parameters: Theory, measurement, uses, limitations and reliability. Pest Manag. Sci. 2002, 58, 419–445. [Google Scholar] [CrossRef] [PubMed]
  106. Travis, C.C.; Etnier, E.L. A Survey of Sorption Relationships for Reactive Solutes in Soil. J. Environ. Qual. 1981, 10, 8–17. [Google Scholar] [CrossRef]
  107. Ayawei, N.; Ebelegi, A.N.; Wankasi, D. Modelling and Interpretation of Adsorption Isotherms. J. Chem. 2017, 2017, 3039817. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Adsorption isotherms of carbaryl, carbofuran and metolachlor on soil L (a) and soil C (b). Error bars represent standard deviation of triplicate samples.
Figure 1. Adsorption isotherms of carbaryl, carbofuran and metolachlor on soil L (a) and soil C (b). Error bars represent standard deviation of triplicate samples.
Molecules 28 05552 g001
Figure 2. Comparison of the experimental data with the Langmuir, Freundlich, Temkin and Dubinin–Radushkevich isotherms for the adsorption of carbaryl (a,b), carbofuran (c,d) and metolachlor (e,f) in soils L and C.
Figure 2. Comparison of the experimental data with the Langmuir, Freundlich, Temkin and Dubinin–Radushkevich isotherms for the adsorption of carbaryl (a,b), carbofuran (c,d) and metolachlor (e,f) in soils L and C.
Molecules 28 05552 g002
Figure 3. Variation of the separation factor (RL) in relation to initial concentration range of the studied (a) carbaryl, (b) carbofuran, (c) metolachlor on soils L and C.
Figure 3. Variation of the separation factor (RL) in relation to initial concentration range of the studied (a) carbaryl, (b) carbofuran, (c) metolachlor on soils L and C.
Molecules 28 05552 g003
Figure 4. Sorption magnitude of carbaryl (a), carbofuran (b) and metolachlor (c) in soils L and C. Error bars represent standard deviation of triplicate samples. Different letters indicate significant differences (p < 0.05).
Figure 4. Sorption magnitude of carbaryl (a), carbofuran (b) and metolachlor (c) in soils L and C. Error bars represent standard deviation of triplicate samples. Different letters indicate significant differences (p < 0.05).
Molecules 28 05552 g004
Figure 5. Comparison of carbaryl (a,b), carbofuran (c,d) and metolachlor (e,f) mass sorbed and desorbed in relation to different initial pesticide concentrations added to soils L and C. Error bars represent standard deviation of triplicate samples.
Figure 5. Comparison of carbaryl (a,b), carbofuran (c,d) and metolachlor (e,f) mass sorbed and desorbed in relation to different initial pesticide concentrations added to soils L and C. Error bars represent standard deviation of triplicate samples.
Molecules 28 05552 g005
Figure 6. Sorption of non-ionic pesticides: (a) carbaryl, (b) carbofuran, (c) metolachlor, studied under a pH range of three to seven in L and C soils. Results are expressed as mean values ± standard deviation (n = 3). Different letters indicate significant differences (p < 0.05).
Figure 6. Sorption of non-ionic pesticides: (a) carbaryl, (b) carbofuran, (c) metolachlor, studied under a pH range of three to seven in L and C soils. Results are expressed as mean values ± standard deviation (n = 3). Different letters indicate significant differences (p < 0.05).
Molecules 28 05552 g006
Table 1. Linear sorption isotherm parameters of carbaryl, carbofuran and metolachlor, approximated from the results of batch experiments on L and C soil.
Table 1. Linear sorption isotherm parameters of carbaryl, carbofuran and metolachlor, approximated from the results of batch experiments on L and C soil.
SoilPesticideLangmuirFreundlichTemkinDubinin–Radushkevich
Qe = (Qmax KLCe/1 + KLCe)Qe = Kf CenfQe = (RT/bT)ln Ce + (RT/bT) lnAQe = Qm exp(−βɛ2)
QmaxKLR2KfnfR2bTAR2QmEaR2
LCarbaryl15.250.060.9881.270.6110.9940.7840.700.9787.820.5000.841
Carbofuran27.930.040.9922.660.4960.9590.3900.400.97318.430.2360.861
Metolachlor32.470.100.9643.290.6020.9890.3220.740.99420.720.5000.909
CCarbaryl32.570.200.99711.420.2480.9510.4675.370.97727.670.7070.937
Carbofuran61.350.040.9785.280.5150.9560.1940.460.97736.310.2890.883
Metolachlor33.780.110.9543.750.6090.9950.3130.910.98320.540.7070.880
Qe–amount of adsorbate at equilibrium (mg kg−1); Ce–equilibrium concentration of the adsorbate (mg L−1); Qmax–maximum adsorption capacity (mg kg−1); KL–Langmuir constant (L mg−1); Kf–Freundlich affinity coefficient (mg (1-n) Ln kg−1); nf–exponential coefficient associated with the isotherm nonlinearity; R–gas constant, 8.31 (J mol−1 K−1); T–the absolute temperature (K); bT–Temkin isotherm constant related to the heat of adsorption (kJ mol−1); A–equilibrium binding constant (L kg−1); Qm–theoretical monolayer sorption capacity (mg g−1); β–constant of the sorption energy (mol2 J −2); ε–Polanyi potential; Ea–mean sorption energy (kJ mol−1).
Table 2. Standard Gibb’s free energy change in tested pesticides on soils L and C.
Table 2. Standard Gibb’s free energy change in tested pesticides on soils L and C.
Soil/PesticideCarbarylΔG° [kJ·mol−1] CarbofuranMetolachlor
loamy sand−8.33−14.10−20.95
loam−35.02−23.97−19.90
Table 3. Leaching indexes calculated for the studied pesticides in soils L and C.
Table 3. Leaching indexes calculated for the studied pesticides in soils L and C.
PesticideL SoilC Soil
KfKoclog KocGUS *KfKoclog KocGUS *
carbaryl1.27134.152.131.8711.42761.512.881.12
carbofuran2.66280.002.452.645.28351.872.552.47
metolachlor3.29346.322.541.463.75250.002.401.60
* GUS-groundwater ubiquity score t ½ was taken from [62].
Table 4. Content of carbon fractions in the investigated loamy sand and loam soils, expressed as % of total organic carbon.
Table 4. Content of carbon fractions in the investigated loamy sand and loam soils, expressed as % of total organic carbon.
SoilTOCCACTECCHACFACR
g 100 g−1% of TOC
L0.956.8761.7640.0721.6931.37
C1.504.6435.6520.1715.4759.71
TOC-total organic carbon, CAC-carbon of low-molecular fulvic fraction, TEC-total extractable carbon consisting of humic (CHA) and fulvic acids (CFA), CR-carbon of humin fraction.
Table 5. Desorption of non-ionic pesticides studied under a pH range of three to seven in L and C soils. Results are expressed as mean values ± standard deviation (n = 3). Different letters indicate significant differences at p ≤ 0.05.
Table 5. Desorption of non-ionic pesticides studied under a pH range of three to seven in L and C soils. Results are expressed as mean values ± standard deviation (n = 3). Different letters indicate significant differences at p ≤ 0.05.
Desorption at:CarbarylCarbofuranMetolachlor
Soil LSoil CSoil LSoil CSoil LSoil C
pH 347.23 a ± 1.857.26 b ± 0.726.36 a ± 1.6110.23 b ± 0.2020.75 a ± 1.751.40 b ± 0.05
pH 451.98 c ± 1.0821.92 e ± 0.6755.86 c ± 0.2912.29 d ± 0.2449.77 c ± 2.761.62 b ± 0.57
pH 559.08 d ± 2.7923.56 e ± 1.6785.11 e ± 1.0524.07 f ± 1.1756.35 d ± 2.375.03 b ± 0.49
pH 756.03 d ± 0.3825.34 e ± 1.3598.01 g ± 1.9033.05 h ± 0.4961.67 d ± 6.7318.84 e ± 0.14
Table 6. Basic properties of pesticides under study 1.
Table 6. Basic properties of pesticides under study 1.
Property/pesticideCarbarylCarbofuranMetolachlor
Chemical structureMolecules 28 05552 i001Molecules 28 05552 i002Molecules 28 05552 i003
Water solubility (mg L−1) at 20 °C9.1322530
Octanol-water partition coefficient, Log P2.361.83.4
1 Data from PubChem open chemistry database and PPDB: Pesticide Properties Database, University of Hertfordshire.
Table 7. Basic physicochemical properties of tested arable L and C soils 1.
Table 7. Basic physicochemical properties of tested arable L and C soils 1.
SoilpHCtotCorgNCorg:NCECSandSiltClaySoil
Texture
(KCl)% cmolc·kg−1%
L7.21.30.950.17.126.476177Sandy Loam
C7.41.61.50.28.442.3473419Loam
1 Previously published by Ćwieląg-Piasecka et al. [10].
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ćwieląg-Piasecka, I. Soil Organic Matter Composition and pH as Factors Affecting Retention of Carbaryl, Carbofuran and Metolachlor in Soil. Molecules 2023, 28, 5552. https://doi.org/10.3390/molecules28145552

AMA Style

Ćwieląg-Piasecka I. Soil Organic Matter Composition and pH as Factors Affecting Retention of Carbaryl, Carbofuran and Metolachlor in Soil. Molecules. 2023; 28(14):5552. https://doi.org/10.3390/molecules28145552

Chicago/Turabian Style

Ćwieląg-Piasecka, Irmina. 2023. "Soil Organic Matter Composition and pH as Factors Affecting Retention of Carbaryl, Carbofuran and Metolachlor in Soil" Molecules 28, no. 14: 5552. https://doi.org/10.3390/molecules28145552

APA Style

Ćwieląg-Piasecka, I. (2023). Soil Organic Matter Composition and pH as Factors Affecting Retention of Carbaryl, Carbofuran and Metolachlor in Soil. Molecules, 28(14), 5552. https://doi.org/10.3390/molecules28145552

Article Metrics

Back to TopTop