Model coupling intraparticle diffusion/sorption, nonlinear sorption, and biodegradation processes

Diffusion, sorption and biodegradation are key processes impacting the efficiency of natural attenuation. While each process has been studied individually, limited information exists on the kinetic coupling of these processes. In this paper, a model is presented that couples nonlinear and nonequilibrium sorption (intraparticle diffusion) with biodegradation kinetics. Initially, these processes are studied independently (i.e., intraparticle diffusion, nonlinear sorption and biodegradation), with appropriate parameters determined from these independent studies. Then, the coupled processes are studied, with an initial data set used to determine biodegradation constants that were subsequently used to successfully predict the behavior of a second data set. The validated model is then used to conduct a sensitivity analysis, which reveals conditions where biodegradation becomes desorption rate-limited. If the chemical is not pre-equilibrated with the soil prior to the onset of biodegradation, then fast sorption will reduce aqueous concentrations and thus biodegradation rates. Another sensitivity analysis demonstrates the importance of including nonlinear sorption in a coupled diffusion/sorption and biodegradation model. While predictions based on linear sorption isotherms agree well with solution concentrations, for the conditions evaluated this approach overestimates the percentage of contaminant biodegraded by as much as 50%. This research demonstrates that nonlinear sorption should be coupled with diffusion/sorption and biodegradation models in order to accurately predict bioremediation and natural attenuation processes. To our knowledge this study is unique in studying nonlinear sorption coupled with intraparticle diffusion and biodegradation kinetics with natural media.     

Comparison of nonideal sorption formulations in modeling the transport of phthalate esters through packed soil columns

Sorption of dimethyl phthalate (DMP), diethyl phthalate (DEP) and dipropyl phthalate (DPP) to two soil materials that vary in organic matter content was investigated using miscible displacement experiments under saturated flow conditions. Generated breakthrough curves (BTCs) were inversely simulated using linear, equilibrium sorption (LE), nonlinear, equilibrium sorption (NL), linear, first-order nonequilibrium sorption (LFO), linear, radial diffusion (LRD), and nonlinear, first-order nonequilibrium sorption (NFO) models. The Akaike information criterion was utilized to determine the preferred model. The LE model could not adequately describe phthalate ester (PE) BTCs in higher organic matter soil or for more hydrophobic PEs. The LFO and LRD models adequately described the BTCs but a slight improvement in curve-fitting was gained in some cases when the NFO model was used. However, none of the models could properly describe the desorptive tail of DPP for the high organic matter soil. Transport of DPP through this soil was adequately predicted when degradation or sorption hysteresis was considered. Using the optimized parameter values along with values reported by others it was shown that the organic carbon distribution coefficient (K(oc)) of PEs correlates well with the octanol/water partition coefficient (K(ow)). Also, a strong relationship was found between the first-order sorption rate coefficient normalized to injection pulse size and compound residence time. A similar trend of timescale dependence was found for the rate parameter in the radial diffusion model. Results also revealed that the fraction of instantaneous sorption sites is dependent on K(ow) and appears to decrease with the increase in the sorption rate parameter.     

The effect of local-scale physical heterogeneity and nonlinear,rate-limited sorption/desorption on contaminant transport in porous media

Nonideal transport of contaminants in porous media has often been observed in laboratory characterization studies. It has long been recognized that multiple processes associated with both physical and chemical factors can contribute to this nonideal transport behavior. To fully understand system behavior, it is important to determine the relative contributions of these multiple factors when conducting contaminant transport and fate studies. In this study, the relative contribution of physical-heterogeneity-related processes versus those of nonlinear, rate-limited sorption/desorption to the observed nonideal transport of trichloroethene in an undisturbed aquifer core was determined through a series of miscible-displacement experiments. The results of experiments conducted using the undisturbed core, collected from a Superfund site in Tucson, AZ, were compared to those obtained from experiments conducted using the same aquifer material packed homogeneously. The results indicate that both physical and chemical factors, specifically preferential flow and associated rate-limited diffusive mass-transfer and rate-limited sorption/desorption, respectively, contributed to the nonideal behavior observed for trichloroethene transport in the undisturbed core. A successful prediction of trichloroethene transport in the undisturbed core was made employing a mathematical model incorporating multiple sources of nonideal transport, using independently determined model parameters to account for the multiple factors contributing to the nonideal transport behavior. The simulation results indicate that local-scale physical heterogeneity controlled the nonideal transport behavior of trichloroethene in the undisturbed core, and that nonlinear, rate-limited sorption/desorption were of secondary importance.     

Sorption-desorption behaviour of 2,4-dichlorophenol by marine sediments

Batch kinetic and isotherm experiments were conducted to determine the sorption-desorption behavior of 2,4-dichlorophenol from seawater solutions by marine sediments containing various amounts of organic carbon (from 1.02% to 12.72% dry weight). The results indicated linear type isotherms for sorption and desorption in all marine sediments studied. The observed difference in linear sorption coefficients between sorption and desorption was indicative of sorption hysteresis. The kinetic experiments showed that equilibrium was established in less than 20 h. The study is significant with respect to sediment remediation in contaminated harbors and coastal areas.     

Modelling of diffusion-limited retardation of contaminants in hydraulically and lithologically nonuniform media

A new reactive transport modelling approach and examples of its application are presented, dealing with the impact of sorption/desorption kinetics on the spreading of solutes, e.g. organic contaminants, in groundwater. Slow sorption/desorption is known from the literature to be strongly responsible for the retardation of organic contaminants. The modelling concept applied in this paper quantifies sorption/desorption kinetics by an intra-particle diffusion approach. According to this idea, solute uptake by or release from the aquifer material is modelled at small scale by a "slow" diffusion process where the diffusion coefficient is reduced as compared to the aqueous diffusion coefficient due to (i) the size and shape of intra-particle pores and (ii) retarded transport of solutes within intra-particle pores governed by a nonlinear sorption isotherm. This process-based concept has the advantage of requiring only measurable model parameters, thus avoiding fitting parameters like first-order rate coefficients.In addition, the approach presented here allows for modelling of slow sorption/desorption in lithologically nonuniform media. Therefore, it accounts for well-known experimental findings indicating that sorptive properties depend on (i) the grain size distribution of the aquifer material and (ii) the lithological composition (e.g. percentage of quartz, sandstone, limestone, etc.) of each grain size fraction. The small-scale physico-chemical model describing sorption/desorption is coupled to a large-scale model of groundwater flow and solute transport. Consequently, hydraulic heterogeneities may also be considered by the overall model. This coupling is regarded as an essential prerequisite for simulating field-scale scenarios which will be addressed by a forthcoming publication.This paper focuses on mathematical model formulation, implementation of the numerical code and lab-scale model applications highlighting the sorption and desorption behavior of an organic contaminant (Phenanthrene) with regard to three lithocomponents exhibiting different sorptive properties. In particular, it is shown that breakthrough curves (BTCs) for lithologically nonuniform media cannot be obtained via simple arithmetic averaging of breakthrough curves for lithologically uniform media. In addition, as no analytical solutions are available for model validation purposes, simulation results are compared to measurements from lab-scale column experiments. The model results indicate that the new code can be regarded as a valuable tool for predicting long-term contaminant uptake or release, which may last for several hundreds of years for some lithocomponents. In particular, breakthrough curves simulated by pure forward modelling reproduce experimental data much better than a calibrated standard first-order kinetics reactive transport model, thus indicating that the new approach is of high quality and may be advantageously used for supporting the design of remediation strategies at contaminated sites where some lithocomponents and/or grain size classes may provide a long-term pollutant source.     

Cesium migration in saturated silica sand and Hanford sediments as impacted by ionic strength

Large amounts of 137Cs have been accidentally released to the subsurface from the Hanford nuclear site in the state of Washington, USA. The cesium-containing liquids varied in ionic strengths, and often had high electrolyte contents, mainly in the form of NaNO3 and NaOH, reaching concentrations up to several moles per liter. In this study, we investigated the effect of ionic strengths on Cs migration through two types of porous media: silica sand and Hanford sediments. Cesium sorption and transport was studied in 1, 10, 100, and 1000 mM NaCl electrolyte solutions at pH 10. Sorption isotherms were constructed from batch equilibrium experiments and the batch-derived sorption parameters were compared with column breakthrough curves. Column transport experiments were analyzed with a two-site equilibrium-nonequilibrium model. Cesium sorption to the silica sand in batch experiments showed a linear sorption isotherm for all ionic strengths, which matched well with the results from the column experiments at 100 and 1000 mM ionic strength; however, the column experiments at 1 and 10 mM ionic strength indicated a nonlinear sorption behavior of Cs to the silica sand. Transport through silica sand occurred under one-site sorption and equilibrium conditions. Cesium sorption to Hanford sediments in both batch and column experiments was best described with a nonlinear Freundlich isotherm. The column experiments indicated that Cs transport in Hanford sediments occurred under two-site equilibrium and nonequilibrium sorption. The effect of ionic strength on Cs transport was much more pronounced in Hanford sediments than in silica sands. Effective retardation factors of Cs during transport through Hanford sediments were reduced by a factor of 10 when the ionic strength increased from 100 to 1000 mM; for silica sand, the effective retardation was reduced by a factor of 10 when ionic strength increased from 1 to 1000 mM. A two order of magnitude change in ionic strength was needed in the silica sand to observe the same change in Cs retardation as in Hanford sediments.     

Sorption and desorption kinetics of diuron, fluometuron, prometryn and pyrithiobac sodium in soils

The sorption and desorption characteristics of four herbicides (diuron, fluometuron, prometryn and pyrithiobac-sodium) in three different cotton growing soils of Australia was investigated. Kinetics and equilibrium sorption and desorption isotherms were determined using the batch equilibrium technique. Sorption was rapid (> 80% in 2 h) and sorption equilibrium was achieved within a short period of time (ca 4 h) for all herbicides. Sorption isotherms of the four herbicides were described by Freundlich equation with an r2 value > 0.98. The herbicide sorption as measured by the distribution coefficient (Kd) values ranged from 3.24 to 5.71 L/kg for diuron, 0.44 to 1.13 L/kg for fluometuron, 1.78 to 6.04 L/kg for prometryn and 0.22 to 0.59 L/kg for pyrithiobac-sodium. Sorption of herbicides was higher in the Moree soil than in Narrabri and Wee Waa soils. When the Kd values were normalised to organic carbon content of the soils (Koc), it suggested that the affinity of the herbicides to the organic carbon increased in the order: pyrithiobac-sodium < fluometuron < prometryn < or = diuron. The desorption isotherms were also adequately described by the Freundlich equation. For desorption, all herbicides exhibited hysteresis and the hysteresis was stronger for highly sorbed herbicides (diuron and prometryn) than the weakly sorbed herbicides (fluometuron and pyrithiobac-sodium). Hysteresis was also quantified as the percentage of sorbed herbicides which is not released during the desorption step (omega = [nad/nde - 1] x 100). Soil type and initial concentration had significant effect on omega. The effect of sorption and desorption properties of these four herbicides on the off-site transport to contaminate surface and groundwater are also discussed in this paper.     

Competitive sorption and diffusion of chromate and sulphate in a flow system with goethite in gel beads

Column experiments and model simulations were employed to evaluate the processes involved in multicomponent solute transport in a system with heterogeneous flow. Column experiments were performed with goethite embedded in polyacrylamide gel beads. The gel forms an immobile water region that can be accessed by diffusion. A two-region transport model with diffusion into spheres was combined with a surface complexation model to predict reactive transport in the goethite-gel bead system. Chromate and sulphate breakthrough curves were measured in a set of transport experiments, along with corresponding changes in the pH of the effluent. Sorption and transport of sulphate and chromate in separate columns were predicted from independently measured sorption parameters. The model overestimated the pH changes in the effluent, possibly because of proton buffering by the polyacrylamide gel. The effect of competitive sorption on transport was examined in experiments with both anions present. The model predicted the effect of competition very well in a system initially equilibrated with sulphate, followed by infiltration with chromate. However, when sulphate was infiltrated after equilibration with chromate, chromate desorption and sulphate adsorption were clearly overestimated by the transport model. The exchange between the more strongly bound chromate and the sulphate added subsequently may be too slow to cause a substantial chromate peak in the effluent. This suggests that the local equilibrium assumption was not applicable in this case.     

Sorption behavior of dibutyl phthalate and dioctyl phthalate by aged refuse

Sorption is a fundamental process controlling the transformation, fate, degradation, and biological activity of hydrophobic organic contaminants in the environment. We investigated the kinetics, isotherms, and potential mechanisms for the sorption of two phthalic acid esters (PAEs), dibutyl phthalate (DBP) and dioctyl phthalate (DOP), on aged refuse. A two-compartment first-order model performed better than a one-compartment first-order model in describing the kinetic sorption of PAEs, with a fast sorption process dominating. Both the Freundlich and Dubinin–Astakhov (DA) models fit the sorption isotherms of DBP and DOP, with the DA model being of a better fit over the range of apparent equilibrium concentrations. The values of the fitting parameters (n, b, E) of the PAEs suggest nonlinear sorption characteristics. Higher predicted partition coefficient values and saturated sorption capacity existed in refuse containing larger quantities of organic matter. The sorption capacity of DOP was significantly higher than that of DBP. PAE sorption was dependent on liquid phase pH. Desorption hysteresis occurred in PAE desorption experiments, especially for the long-chain DOP. PAEs may therefore be a potential environmental risk in landfill.     

The impact of contact time on pyrene sorptive behavior by a sandy-loam soil

Batch experiments with pyrene (PYR) were conducted to quantify the effect of contact time on its sorption and desorption behavior by a sandy-loam soil. Twenty-four and 48 h contact times were chosen for the nonequilibrium conditions and 240 h for the pseudoequilibrium study. All times was selected based on the kinetic results. The nonlinear, pseudoequilibrium sorption isotherm was fit to a two-stage Freundlich model: 3-7 mg/l for the first stage and 7-15 mg/l for the second stage. A substantial fraction of the sorbed PYR was not desorbed within the given desorption time. The reason of hysteresis was found to be a sorption enhancement due to soil hydration which provided more sorption sites. A desorption enhancement at 240-h desorption steps was attributed to the increased dissolved organic matter evolution. This study also found that both soil organic matter and clay materials had an equal role in PYR sorption enhancement and desorption resistance.     

Sorption and desorption kinetics of diuron,fluometuron, prometryn and pyrithiobac sodium in soils

Abstract

The sorption and desorption characteristics of four herbicides (diuron, fluometuron, prometryn and pyrithiobac‐sodium) in three different cotton growing soils of Australia was investigated. Kinetics and equilibrium sorption and desorption isotherms were determined using the batch equilibrium technique. Sorption was rapid (> 80% in 2 h) and sorption equilibrium was achieved within a short period of time (ca 4 h) for all herbicides. Sorption isotherms of the four herbicides were described by Freundlich equation with an r² value > 0.98. The herbicide sorption as measured by the distribution coefficient (K_d) values ranged from 3.24 to 5.71 L/kg for diuron, 0.44 to 1.13 L/kg for fluometuron, 1.78 to 6.04 L/kg for prometryn and 0.22 to 0.59 L/kg for pyrithiobac‐sodium. Sorption of herbicides was higher in the Moree soil than in Narrabri and Wee Waa soils. When the K_d values were normalised to organic carbon content of the soils (K_oC), it suggested that the affinity of the herbicides to the organic carbon increased in the order: pyrithiobac‐sodium < fluometuron < prometryn < diuron. The desorption isotherms were also adequately described by the Freundlich equation. For desorption, all herbicides exhibited hysteresis and the hysteresis was stronger for highly sorbed herbicides (diuron and prometryn) than the weakly sorbed herbicides (fluometuron and pyrithiobac‐sodium). Hysteresis was also quantified as the percentage of sorbed herbicides which is not released during the desorption step ω = [n_ad / n_de ‐1] x 100). Soil type and initial concentration had significant effect on ω. The effect of sorption and desorption properties of these four herbicides on the off‐site transport to contaminate surface and groundwater are also discussed in this paper.     

Capability and limitations of first-order and diffusion approaches to describe long-term sorption of chlortoluron in soil

This paper compares the capability of a first-order and a spherical diffusion model to describe and predict long-term sorption and desorption processes of chlortoluron in two soils. Chlortoluron sorption was investigated at different time scales utilizing one rate experiment (120 days) and two sorption/desorption experiments. Experimental periods for sorption and desorption were set to 1 day (five desorption steps) and 30 days (three desorption steps), respectively. Upon fitting, the two models satisfactorily described the whole set of data. The spherical diffusion model performed better than the first-order model. We then tested the predictive capability of the models by predicting 30-day sorption/desorption data using kinetic parameters fitted on 1-day sorption/desorption data only. While the spherical diffusion model was able to predict the 30-day data set, the first-order model failed completely. Fitting both models to subsets of the data corresponding to different experimental time scales revealed that the rate parameter as well as the Freundlich coefficient of the first-order model are strongly time-dependent--a property that is not shared by parameters of the spherical diffusion model. The apparent stability of the spherical diffusion model with regard to time dependency of its parameters indicates that sorptive uptake may be diffusion-controlled. This also explains the models greater predictive power across different time scales compared to the first-order model. Finally, we investigate the suitability of solute class specific log-linear relationships between the first-order rate parameter and the Freundlich coefficient presented by earlier researchers in the light of the time dependency observed for the parameters of the first-order model.     

Nonideal transport of reactive contaminants in heterogeneous porous media: 7. distributed-domain model incorporating immiscible-liquid dissolution and rate-limited sorption/desorption

The purpose of this work is to present a distributed-domain mathematical model incorporating the primary mass-transfer processes that mediate the transport of immiscible organic liquid constituents in water-saturated, locally heterogeneous porous media. Specifically, the impact of grain/pore-scale heterogeneity on immiscible-liquid dissolution and sorption/desorption is represented in the model by describing the system as comprising a continuous distribution of mass-transfer domains. With this conceptualization, the distributions of the initial dissolution rate coefficient and the sorption/desorption rate coefficient are represented as probability density functions. Several sets of numerical experiments are conducted to examine the effects of heterogeneous dissolution and sorption/desorption on contaminant transport and elution. Four scenarios with different combinations of uniform/heterogeneous rate-limited dissolution and uniform/heterogeneous rate-limited sorption/desorption are evaluated. The results show that both heterogeneous rate-limited sorption/desorption and heterogeneous rate-limited dissolution can significantly increase the time or pore volumes required to elute immiscible-liquid constituents from a contaminated porous medium. However, sorption/desorption has minimal influence on elution behavior until essentially all of the immiscible liquid has been removed. For typical immiscible-liquid constituents that have relatively low sorption, the asymptotic elution tailing produced by heterogeneous rate-limited sorption/desorption begins at effluent concentrations that are several orders of magnitude below the initial steady-state concentrations associated with dissolution of the immiscible liquid. Conversely, the enhanced elution tailing associated with heterogeneous rate-limited dissolution begins at concentrations that are approximately one-tenth of the initial steady-state concentrations. Hence, dissolution may generally control elution behavior of immiscible-liquid constituents in cases wherein grain/pore-scale heterogeneity significantly influences both dissolution and sorption/desorption.     

Transport of reactive colloids and contaminants in groundwater: effect of nonlinear kinetic interactions

Transport of reactive colloids in groundwater may enhance the transport of contaminants in groundwater. Often, the interpretation of results of transport experiments is not a simple task as both reactions of colloids with the solid matrix and reactions of contaminants with the solid matrix and mobile and immobile colloids may be time dependent and nonlinear. Further colloid transport properties may differ from solute transport properties. In this paper, a one-dimensional model for coupled and contaminant in a porous medium (COLTRAP) is presented together with simulation results. Calculated breakthrough curves (BTC's) during contamination and decontamination show systematically the effect of nonlinear and kinetic interactions on contaminant transport in the presence of reactive colloids, and the effect of colloid transport properties that differ from solute transport properties. It is shown that in case of linear kinetic reactions, the rate of exchange of mobile and immobile colloids have a large impact on the shape of BTC's even if the solid matrix is saturated with respect to colloids. BTC's during the contamination and decontamination phase have identical shapes in this case. Moreover, the slow reactions of contaminants and colloids may lead to unretarded breakthrough of contaminants. Independent of reaction rates, nonlinear reactions lead to BTC's that are steeper during contamination than in the linear case. A characteristic aspect of nonlinear sorption is that shapes of BTC's differ during the contamination and decontamination phase. It has been observed that shapes of some of the simulated adsorption and desorption curves are similar as shapes found in experiments reported in literature. This stresses the importance of incorporating both kinetics and nonlinearity in models for coupled colloid and contaminant transport and the capability of COLTRAP to interpret experimental results. Finally, to figure out whether nonlinear processes play a role, it is very important to consider both contamination and decontamination in transport experiments.     

Modeling colloid-facilitated transport of multi-species contaminants in unsaturated porous media

Colloid-facilitated transport has been recognized as a potentially important and overlooked contaminant transport process. In particular, it has been observed that conventional two phase sorption models are often unable to explain transport of highly sorbing compounds in the subsurface appropriately in the presence of colloids. In this study a one-dimensional model for colloid-facilitated transport of chemicals in unsaturated porous media is developed. The model has parts for simulating coupled flow, and colloid transport and dissolved and colloidal contaminant transport. Richards' equation is solved to model unsaturated flow, and the effect of colloid entrapment and release on porosity and hydraulic conductivity of the porous media is incorporated into the model. Both random sequential adsorption and Langmuir approaches have been implemented in the model in order to incorporate the effect of surface jamming. The concept of entrapment of colloids into the air-water interface is used for taking into account the effect of retardation caused due to existence of the air phase. A non-equilibrium sorption approach with options of linear and Langmuir sorption assumptions are implemented that can represent the competition and site saturation effects on sorption of multiple compounds both to the solid matrix and to the colloidal particles. Several demonstration calculations are performed and the conditions in which the non-equilibrium model can be approximated by an equilibrium model are also studied.     

BDE-47 sorption and desorption to soil matrix in single- and binary-solute systems

Three loamy-clay soil samples (LC1-3) with different properties were collected as the geosorbents to preliminarily investigate the sorption and desorption of 2,2′,4,4′-tetrabromodiphenyl ether (BDE-47) in single system and binary system with the presence of decabromodiphenyl ether (BDE-209), which can provide information in order to further understand the sorption mechanisms and evaluate the adsorption sites. A concentration of 10 μg L⁻¹ BDE-209 suppressed the sorption of BDE-47, and the trend became more and more significant with the increase of BDE-47 equilibrium concentration, however, BDE-47 caused no competitive effect on BDE-209 sorption, which was related with the better accessibility of more hydrophobic molecules to adsorption sites. In the binary system, nonlinearity of the BDE-47 sorption isotherms for the three samples changed in different ways, which originated from the varied soil properties. Desorption hysteresis was observed in all cases, which was estimated due to irreversible surface adsorption between sorbent and sorbate. BDE-209 made desorption of BDE-47 more hysteretic from soil samples, which was estimated to be ascribed to the accelerated sorbent state transition and new sites creation caused by BDE-209 sorption.     

Nonequilibrium sorption of phenols onto geosorbents: the impact of pH on intraparticle mass transfer

While the pH effect on sorption equilibrium of weak acids on natural sorbents was investigated in a number of studies, less is known about the pH dependence of sorption kinetics. This paper investigates the impact of pH on sorption kinetics during the transport of some selected phenols through a sandy aquifer material. Breakthrough curves measured in column experiments were analyzed using a mass transfer based nonequilibrium model designated as dispersed flow, film and particle diffusion model (DF-FPDM). In this model, the rate limiting intraparticle diffusion is characterized by the mass transfer coefficient, kSaV, which can be determined from breakthrough curves by curve fitting. The experimental results indicate that the kSaV is pH-dependent and inversely correlated with the pH-dependent distribution coefficient, K(d,app). Regression equations are presented that may be used to estimate approximate values of intraparticle mass transfer coefficients on the basis of experimentally determined or LFER predicted distribution coefficients.     

Application of the method of temporal moments to interpret solute transport with sorption and degradation

In this note, we applied the temporal moment solutions of [Das and Kluitenberg, 1996. Soil Sci. Am. J. 60, 1724] for one-dimensional advective-dispersive solute transport with linear equilibrium sorption and first-order degradation for time pulse sources to analyse soil column experimental data. Unlike most other moment solutions, these solutions consider the interplay of degradation and sorption. This permits estimation of a first-order degradation rate constant using the zeroth moment of column breakthrough data, as well as estimation of the retardation factor or sorption distribution coefficient of a degrading solute using the first moment. The method of temporal moment (MOM) formulae was applied to analyse breakthrough data from a laboratory column study of atrazine, hexazinone and rhodamine WT transport in volcanic pumice sand, as well as experimental data from the literature. Transport and degradation parameters obtained using the MOM were compared to parameters obtained by fitting breakthrough data from an advective-dispersive transport model with equilibrium sorption and first-order degradation, using the nonlinear least-square curve-fitting program CXTFIT. The results derived from using the literature data were also compared with estimates reported in the literature using different equilibrium models. The good agreement suggests that the MOM could provide an additional useful means of parameter estimation for transport involving equilibrium sorption and first-order degradation. We found that the MOM fitted breakthrough curves with tailing better than curve fitting. However, the MOM analysis requires complete breakthrough curves and relatively frequent data collection to ensure the accuracy of the moments obtained from the breakthrough data.     

Sorption of hydrophobic organic pollutants in saturated soil systems

Sorption equilibria and rates were characterized for a matrix of four aquifer sands and two slightly to moderately hydrophobic organic solutes (nitrobenzene and lindane), and the effects of sorption on the behavior of these solutes in saturated systems of the soils were determined. Experimental data were used to test and evaluate a variety of mathematical models for predicting contaminant fate and transport in groundwater systems.Observed equilibrium relationships between soil and solution phase solute concentrations were found to be described best by the nonlinear Freundlich isotherm model. It was further determined that the sorption process in the systems tested is rate controlled, requiring several days to approach equilibrium in completely mixed batch reactors. Subsequent modeling of solute transport in continuous flow soil column reactors was found to be most successful when rate-controlled models were used, the best results were obtained with a dual-resistance model incorporating the coupled mass transport steps of boundary-layer and intraparticle diffusion.     

Sorption-desorption of cadmium in aqueous palygorskite, sepiolite, and calcite suspensions: isotherm hysteresis

Sorption isotherms have been widely used to assess the heavy metal retention characteristics of soil particles. Desorption behavior of the retained metals, however, usually differ from that of sorption, leading to a lack of coincidence in the experimentally obtained sorption and desorption isotherms. In this study, we examine the nonsingularity of cadmium (Cd) sorption–desorption isotherms, to check the possible hysteresis and reversibility phenomena, in aqueous palygorskite, sepiolite and calcite systems. Sorption of Cd was carried out using a 24-h batch equilibration experiment with eight different Cd solution concentrations, equivalent to 20–100% of maximum sorption capacity of each mineral. Immediately after sorption, desorption took place using successive dilution method with five consecutive desorption steps. Both Cd sorption and desorption data were adequately described by Freundlich equation (0.81 < r² < 0.99). The sorption and desorption reactions, however, did not provide the same isotherms, indicating that hysteresis occurred in Cd sorption–desorption processes. The extent of hysteresis was quantified based on the differences obtained from sorption and desorption isotherms regarding the amount of Cd sorbed, the Freundlich exponent, and the Cd distribution coefficient. The results revealed that, sepiolite possessed the most hysteretic behavior among the minerals studied. Calcite showed much smaller hysteresis compared to the other two silicate clays at low Cd surface load, but its hysteresis indices significantly increased, and exceeded that of palygorskite, as the amount of Cd in the systems increased. The average amount of Cd released after five desorption steps, was 13.8%, 2.2% and 3.6% for the palygorskite, sepiolite and calcite, respectively, indicating that a large portion of Cd was irreversibly retained by the minerals.