Transport of polar and non-polar volatile compounds in polystyrene foam and oriented strand board

Transport of hexanal and styrene in polystyrene foam (PSF) and oriented strand board (OSB) was characterized. A microbalance was used to measure sorption/desorption kinetics and equilibrium data. While styrene transport in PSF can be described by Fickian diffusion with a symmetrical and reversible sorption/desorption process, hexanal transport in both PSF and OSB exhibited significant hysteresis, with desorption being much slower than sorption. A porous media diffusion model that assumes instantaneous local equilibrium governed by a nonlinear Freundlich isotherm was found to explain the hysteresis in hexanal transport. A new nonlinear sorption and porous diffusion emissions model was, therefore, developed and partially validated using independent chamber data. The results were also compared to the more conventional linear Fickian-diffusion emissions model. While the linear emissions model predicts styrene emissions from PSF with reasonable accuracy, it substantially underestimates the rate of hexanal emissions from OSB. Although further research and more rigorous validation is needed, the new nonlinear emissions model holds promise for predicting emissions of polar VOCs such as hexanal from porous building materials.     

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.     

Prediction of mass-transfer coefficient for solute transport in porous media

Several previously reported laboratory studies related to transport of solutes through packed columns were utilized to develop predictive relationships for mass-transfer rate coefficient. The data were classified into two groups: those obtained under rate-limited mass transfer between mobile and immobile water regions (physical nonequilibrium conditions), and those derived from rate-limited mass transfer between instantaneous and slow sorption sites (sorption nonequilibrium conditions). The mass-transfer coefficient in all these studies was obtained by fitting breakthrough data to a transport model employing a first-order rate limitations with a "constant" mass-transfer coefficient, independent of flow conditions. This study demonstrated that the mass-transfer coefficient in these models is dependent on system parameters including pore-water velocity, length-scale, retardation coefficient, and particle or aggregate size. Predictive relationships were developed, through regression analysis, relating mass-transfer coefficient to residence time. The developed relationships adequately estimated previously reported field mass-transfer values. Successful simulations of field desorption data reported by Bahr [J. Contam. Hydrol. 4 (1989) 205] further demonstrate the potential applicability of the developed relationships.     

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.     

A method for estimating distributions of mass transfer rate coefficients with application to purging and batch experiments

Mass transfer between aquifer material and groundwater is often modeled as first-order rate-limited sorption or diffusive exchange between mobile zones and immobile zones with idealized geometries. Recent improvements in experimental techniques and advances in our understanding of pore-scale heterogeneity demonstrate that two (or even a few) rate coefficients are insufficient in many cases. Here, we investigate a piece-wise linear model for a continuous distribution of rate coefficients, that has several advantages over previously used ‘statistical' distribution models (with functional form from gamma or lognormal PDF's): (1) distributions of arbitrary, even bimodal, shapes can be represented; (2) linear estimation methods can be applied to determine the distribution from experimental data; (3) the uncertainty in the distribution can be determined for each of its sections; and (4) the relationship between the time scales of available data and those of estimatable mass transfer processes can be investigated. A statistical model refinement algorithm is presented that reduces the number of parameters (sections of the piece-wise linear model) to the admissible minimum. We show that purging experiments allow estimation of a wider zone of the rate distribution than do batch experiments, and hence will provide predictions that are accurate over a wider range of time scales. Finally, in an application to TCE gas-purging desorption data, the piece-wise linear rate-distribution model has a higher probability of being adequate than those using a gamma or lognormal distribution or a single rate coefficient.     

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.     

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 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.     

Adsorption behavior of reactive dye in aqueous solution on chemical cross-linked chitosan beads

A batch system was applied to study the adsorption of reactive dye (reactive red 189) from aqueous solutions by cross-linked chitosan beads. The ionic cross-linking reagent sodium tripolyphosphate was used to obtain more rigid chitosan beads. To stabilize chitosan in acid solutions, chemical cross-linking reagent epichlorohydrin (ECH), glutaraldehyde and ethylene glycol diglycidyl ether was used and ECH shows a higher adsorption capacity. The Langmuir and Freundlich adsorption models were applied to describe the equilibrium isotherms at different particle sizes and isotherm constants were determined. The Langmuir model agrees very well with experimental data and its calculated maximum monolayer adsorption capacity has very large value of 1802-1840 (g/kg) at pH 3.0, 30 degrees C. The kinetics of the adsorption with respect to the initial dye concentration, temperature, pH, ionic strength, and wet/dry beads were investigated. The pseudo-first-order, second-order kinetic models and intraparticle diffusion model were used to describe the kinetic data and the rate constants were evaluated. The dynamical data fit well with the second-order kinetic model, except for the dry beads fitting better with the first-order model. The adsorption capacity increases largely with decreasing solution pH or with increasing initial dye concentration. Thermodynamic parameters such as change in free energy (deltaG(0)), enthalpy (deltaH(0)), entropy (deltaS(0)) and activation energy were also determined. The adsorption mechanism is shown to be the electrostatic interactions between the dye and chitosan beads. The desorption data shows that the removal percent of dye RR 189 from the cross-linked chitosan beads is 63% in NaOH solutions at pH 10.0, 30 degrees C. The desorbed chitosan beads can be reused to adsorb the dye and to reach the same capacity as that before desorption.     

Modeling the effects of concentration history on the slow desorption of trichloroethene from a soil at 100% relative humidity

To determine the effects of concentration history on slow sorption processes, desorption kinetic profiles for trichloroethene (TCE) were measured for a soil at 100% relative humidity subject to different exposure concentrations and exposure times. Exposure concentrations ranged from 1% to 80% of the saturation vapor pressure (Ps) for TCE, and exposure times (i.e., time allowed for sorption before desorption begins) ranged from 1 to 96 days. A spherical diffusion model based on a gamma distribution of sorption rates and a gamma distribution of desorption rates was developed and applied to the data. At 80% P/Ps, the entire gamma distributions of sorption and desorption rates were available for TCE. In accordance with a micropore filling mechanism, the fraction of these distributions available for TCE sorption decreased with decreasing P/Ps. Experimental results are consistent with a micropore-filling mechanism, where the amount of slow desorbing mass decreased with decreasing exposure time, and the fraction of slow desorbing sites filled decreased with decreasing exposure concentration. Simulation results suggest that diffusion limits the rates that micropores fill, and that rates of sorption and desorption for soil contaminated at smaller values of P/Ps are, on average, less than those at larger values of P/Ps (i.e., slow desorption rates are a function of exposure concentration). Simulation results also suggest that the model adequately describes the effects of exposure concentration and exposure time on the rates of sorption and desorption, but not on the capacity of the slow sites for TCE. This work is important because contaminant concentrations in the subsurface vary in space and time, and the proposed model represents a new and mechanistically based approach to capture the effects of this heterogeneity on slow desorption.     

Kinetic and isotherm error optimization studies for adsorption of atrazine and imidacloprid on bark of Eucalyptus tereticornis L.

The aim of this study was to establish the bark of Eucalyptus tereticornis L. (EB) as a low cost bio-adsorbent for the removal of imidacloprid and atrazine from aqueous medium. The pseudo-first-order (PFO), pseudo-second-order (PSO), Elovich and intra-particle diffusion (IPD) models were used to describe the kinetic data and rate constants were evaluated. Adsorption data was analysed using ten 2-, 3- and 4-parameter models viz. Freundlich, Jovanovic, Langmuir, Temkin, Koble–Corrigan, Redlich–Peterson, Sips, Toth, Radke–Prausnitz, and Fritz-Schluender isotherms. Six error functions were used to compute the best fit single component isotherm parameters by nonlinear regression analysis. The results showed that the sorption of atrazine was better explained by PSO model, whereas the sorption of imidacloprid followed the PFO kinetic model. Isotherm model optimization analysis suggested that the Freundlich along with Koble–Corrigan, Toth and Fritz-Schluender were the best models to predict atrazine and imidacloprid adsorption onto EB. Error analysis suggested that minimization of chi-square (χ²) error function provided the best determination of optimum parameter sets for all the isotherms.     

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.     

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.     

Sorption and desorption dynamics of aroclor 1242 to natural sediment

Sediment sorption and desorption processes are important in determining the movement and fate of persistent organic compounds in aquatic systems. Batch experiments show that after an initial one week uptake period, continual release of Aroclor 1242 from sediment occurs over a six-month period. These observations suggest that a two-stage kinetic model, rather than the conventional equilibrium model, is more appropriate for representing sediment uptake and release processes. Additional batch studies were used to measure short- and long-term rate coefficients for these processes. Simulation studies, with multiple sediment and contaminant inputs, indicate that over a 16-day period a kinetic model better matches the experimental data than do three other equilibrium-based sorption/desorption models. Further long-term simulations demonstrate that the kinetic model, rather than equilibrium models, more adequately account for the persistence of organic contaminants in sediment.     

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.     

Competitive sorption of cis-DCE and TCE in silica gel as a model porous mineral solid

The competitive sorption of 1,2-cis-dichloroethene (cis-DCE) and trichloroethene (TCE) was investigated by means of column experiments using a model porous mineral solid represented by silica gel. The experimental isotherms were obtained by employing a chromatographic method. The competitive sorption isotherms were modelled with the extended Freundlich and extended Langmuir isotherms, using the parameters from single-solute experiments. The breakthrough curves were modelled with the advection-dispersion transport equation coupled with the lumped pore diffusion model. The best results were obtained when the extended Freundlich isotherm was employed. The competitive sorption was revealed with the presence of an overshoot in the breakthrough curve of cis-DCE and a decrease in the degree of sorption of cis-DCE (20%) and TCE (12%). A linear dependency of the overshoot with an increase in the concentration of cis-DCE at a fixed concentration of TCE was observed, between 16% and 20%, and at least at concentrations <6 mg L(-1) in the liquid phase. The displaced molecules of cis-DCE by TCE were accumulated through the column causing its overshoot; thus short columns may hinder its observation. Thermodynamic analysis shows an exothermic adsorption process of -34 to -41 kJ mol(-1), which is enhanced by sorption in micropores. The Gibbs free energy is positive for cis-DCE in the multi-component case, due to its displacement by TCE.     

Rigorous uncertainty assessment in contaminant transport inverse modelling: a case study of fluoride diffusion through clay liners

Inverse methods used in assessing landfill liner design have not yet taken advantage of current developments in inverse procedures. Here, a method for inverting contaminant transport models is presented including a general error model and procedures for differentially weighted multiple response regression. General error models are employed in cases where the residuals are heteroscedastic and correlated, and lead to valid inference on model parameter and predictive uncertainty. The Shuffled Complex Evolution algorithm is used to optimise model parameters. Model parameter uncertainty is assessed by exploring the posterior probability distribution with the Metropolis algorithm, a Markov chain Monte Carlo sampling method. The inverse method is applied to simultaneously determine the sorption and diffusion parameters from laboratory diffusion cell experiments. In these experiments, fluoride migration through kaolin clays was measured by sampling the source and collector cells over time. To uniquely determine the transport model parameters, it was necessary to simultaneously fit the observed data from two independent diffusion cell experiments with different initial concentrations. The jointly fitted transport model parameters compared well with those fitted to independent batch experiments.     

Cu sorption on Phragmites australis leaf and stem litter: a kinetic study

Decaying organic matter plays an important role in the cycling of metals in wetland ecosystems. Sorption kinetics of Cu(II) on Phragmites australis leaf and stem litter were studied. Fresh leaf and stem litter was sampled from a surface flow wetland at the end of the growing season. The effect of decomposition stage was studied with litter that had been decomposing for a period of 5 months. The Lagergren pseudo-first-order model, the pseudo-second-order model, the Elovich equation and two diffusion models based on spherical intra-particle diffusion were fitted to the experimental data. The sorption capacity was significantly affected by the decomposition of the litter. The sorption process was best described by the pseudo-second-order kinetics (R(2) approximately 0.99) but the rate constant was strongly dependent on the initial Cu concentration. The intra-particle diffusion model fitted the data only slightly less (R(2)>0.95) than the pseudo-second-order model. A theoretical comparison revealed that the good fit with the pseudo-second-order kinetics could be indicative of intra-particle diffusion. Sorption kinetics observed for the leaf and stem litter at different metal concentrations showed a fast initial sorption followed by a slow sorption phase.     

Sorption and desorption by ideal two-compartment systems: unusual behavior and data interpretation problems

This paper examines the current practices of fitting curves to sorption, desorption, and equilibrium data obtained from laboratory experiments. Systems of equations incorporating Freundlich isotherms and first-order kinetics for two different idealized sorbents, one "fast" and one "slow," were solved numerically to produce "data". Two-compartment curves were then fit to the data by nonlinear regression, and the parameters computed by the regression are compared with the original parameters used to produce the data. The results show that a sorbent with fast kinetics will not steadily accumulate sorbate until it reaches the equilibrium value but will overshoot equilibrium, accumulating an excess of sorbate. This overshoot will cause the sorption rates for both sorbents and the distribution between the fast and slow sorbents to be estimated incorrectly. The system may appear to be at equilibrium by external measures, but sorbate will slowly be redistributing from the fast to the slow sorbent. An isotherm constructed from data acquired during this process will have an incorrect coefficient and exponent. Consequently, the meaning of the results obtained by curve fitting may often be questionable and may say little about the phenomena occurring within the sorbate-sorbent-liquid system. Possible physical explanations for the effects observed are offered.     

Effect of drying on the desorption of diuron and terbuthylazine from natural soils

This work was initiated to study the effects of climate induced soil water status variations which can reach extreme values under natural conditions on the sorption process of hydrophobic organic compounds. Based on the classical slurry batch methodology an approach is developed that allows the fast and careful complete drying of soil suspensions (microwave technique). Classical adsorption experiments were followed by three desorption steps with and without drying cycles. Drying and re-wetting enhanced the sorption-desorption hysteresis and Freundlich adsorption coefficients increased from 5.9 to 16 and 5.2 to 21 over three drying cycles for diuron and terbuthylazine respectively. Assuming the validity of a dual stage adsorption process, model evaluation suggests that drying is as a shrinking-like process leading to conformational changes of the dominant sorbent (soil organic matter) which restrict the intra-micro-particle diffusion. Rewetting only leads to a partial recovery of the diffusional pore space.