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.     

Heterogeneity of chlorinated hydrocarbon sorption properties in a sandy aquifer

Hydraulic conductivity and sorption coefficients for chlorinated hydrocarbons (chloroform, carbon tetrachloride and tetrachloroethylene) were evaluated for 216 sediment samples collected across a 15 m transect and a 21 m depth interval in a contaminated aquifer near Schoolcraft, Michigan. Relationships between hydraulic conductivity, linear sorption partition coefficients, grain size classes, and spatial location were investigated using linear regression analysis and geostatistical techniques. Clear evidence of layering was found in sorption properties, hydraulic conductivity and grain sizes. Conductivity correlated well with grain size, as expected, but sorption varied inversely with grain size, contrary to some previous reports. No significant correlation was found between sorption properties and hydraulic conductivity. This is likely due to the unexpected presence of small amounts of highly sorptive coal-like solids, which dominate the sorption behavior but have little effect on conductivity. The results demonstrate that recent findings regarding the high sorption capacity of coal materials found in soils can exert a controlling influence on contaminant transport. Designers of in situ remediation systems should be cautioned that 1) it is not reasonable to assume that sorption capacity and hydraulic conductivity are related, 2) sorption capacity and hydraulic conductivity are critical measurements for contaminant site characterization and subsequent transport modeling, 3) estimating sorption capacity from organic carbon measurement may lead to greater errors than performing sorption isotherms, and 4) it is more important to characterize vertical heterogeneity rather than horizontal heterogeneity because both sorption and hydraulic conductivity are correlated across longer distances in the horizontal plane.     

Colloid-facilitated transport of radionuclides through fractured media

The sorption of radionuclides on natural colloids may significantly modify their transport behaviour through fractured media, since radionuclides bound to colloids may not be subject to the important retardation mechanisms of matrix diffusion and sorption onto pore surfaces. This paper reports on theoretical and experimental work aimed at assessing the relevance of colloid-facilitated transport to repository safety analyses, with specific reference to the Swiss case. Transport models are presented, developed in conjunction with field- and laboratory-based studies of deep groundwater in the crystalline basement of northern Switzerland, in which colloid size distributions, compositions and sorption properties have been measured. Various potential mechanisms giving rise to both reversible and irreversible sorption are discussed. In the first case, a simple approach is examined which is based on previously reported models of colloid transport and assumes reversible, linear sorption on colloids, for which experimental data have been obtained. It is shown that transport of radionuclides would not, in general, be significantly enhanced because of this process. A more recently developed and more complex model is then described incorporating irreversible sorption, in which case the transport of radionuclides tends to be strongly dependent on the extent of colloid-fracture wall interaction.     

A sorption kinetics model for arsenic adsorption to magnetite nanoparticles

Introduction  

Arsenic is a well known water contaminant that causes toxicological and carcinogenic effects. In this work magnetite nanoparticles were examined as possible arsenic sorbents. The objective of this work was to develop a sorption kinetics model, which could be used to predict the amount of arsenic adsorbed by magnetite nanoparticles in the presence of naturally occurring species using a first-order rate equation, modified to include adsorption, described by a Langmuir isotherm.     

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.     

Kinematic analysis of solute mass flows in rock fractures with spatially random parameters

Field data of physical properties in heterogeneous crystalline bedrock, like porosity and fracture aperture, is associated with uncertainty that can have a significant impact on the analysis of solute transport in rock fractures. Solutions to the central temporal moments of the residence time probability density function (PDF) are derived in a closed form for a solute Dirac pulse. The solutions are based on a model that takes into account advection along the fracture plane, diffusion into the rock matrix and sorption kinetics in the rock matrix. The most relevant rock properties including fracture aperture and several matrix properties as well as flow velocity are assumed to be spatially random along transport pathways. The mass transport is first solved in a general form along one-dimensional pathways, but the results can be extended to multi-dimensional flows simply by substituting the expected travel time for inert water parcels. Based on data obtained with rock samples taken at Asp? Hard Rock Laboratory in Sweden, the solutions indicate that the heterogeneity of the rock properties contributes to increasing significantly both the variance and the skewness of the residence time probability density function for a pulse travelling in a fracture. The Asp? data suggests that the bias introduced in the variance of the residence time PDF by neglecting the effect of heterogeneity of the rock properties on the radionuclide migration is very large for fractures thinner than a few tenths of a millimetre.     

Mixing-controlled biodegradation in a toluene plume--results from two-dimensional laboratory experiments

Various abiotic and biotic processes such as sorption, dilution, and degradation are known to affect the fate of organic contaminants, such as petroleum hydrocarbons in saturated porous media. Reactive transport modeling of such plumes indicates that the biodegradation of organic pollutants is, in many cases, controlled by mixing and therefore occurs locally at the plume's fringes, where electron donors and electron-acceptors mix. Herein, we aim to test whether this hypothesis can be verified by experimental results obtained from aerobic and anaerobic degradation experiments in two-dimensional sediment microcosms. Toluene was selected as a model compound for oxidizable contaminants. The two-dimensional microcosm was filled with quartz sand and operated under controlled flow conditions simulating a contaminant plume in otherwise uncontaminated groundwater. Aerobic degradation of toluene by Pseudomonas putida mt-2 reduced a continuous 8.7 mg L(-1) toluene concentration by 35% over a transport distance of 78 cm in 15.5 h. In comparison, under similar conditions Aromatoleum aromaticum strain EbN1 degraded 98% of the toluene infiltrated using nitrate (68.5+/-6.2 mg L(-1)) as electron acceptor. A major part of the biodegradation activity was located at the plume fringes and the slope of the electron-acceptor gradient was steeper during periods of active biodegradation. The distribution of toluene and the significant overlap of nitrate at the plume's fringe indicate that biokinetic and/or microscale transport processes may constitute additional limiting factors. Experimental data is corroborated with results from a reactive transport model using double Monod kinetics. The outcome of the study shows that in order to simulate degradation in contaminant plumes, detailed data sets are required to test the applicability of models. These will have to deal with the incorporation of existing parameters coding for substrate conversion kinetics and microbial growth.     

Effects of soil moisture and physical-chemical properties of organic pollutants on vapor-phase transport in the vadose zone

Vapor-phase transport of organic pollutants is one of the important pathways in the distribution and attenuation of volatile organic compounds in the vadose zone. In this study, the impact of vapor-phase partitioning and of the physical-chemical properties of organic pollutants on vapor-phase transport was assessed. An experimentally derived relationship to predict vapor sorption for a variety of soil types under varying soil moisture conditions was incorporated into the two-dimensional finite-element model, Vocwaste. The revised model was then used to simulate the transport of volatile organics. Vapor-phase partitioning in the model accounted for vapor uptake by sorption onto moist mineral surfaces as well as sorption at the liquid-solid interface and dissolution into soil water. Under dry conditions, vapor-phase sorption of volatile organic pollutants was shown to have a retarding effect on transport of organic vapors. However, for shallow, contaminated soils, volatilization was controlled by vapor diffusion, even under dry conditions where vapor-phase sorption was high. The influence of Henry's law constant and of the aqueous-phase (solid-liquid) partition coefficient for volatile organic pollutants was considered in the simulations. Volatilization of organic vapors was shown to be favored for contaminants with high Henry's law constants and low aqueous-phase partitioning coefficients. Because of the interdependence of these two physical-chemical properties, individual properties of the contaminant should not be considered in isolation in the evaluation of vapor transport.     

Competing TCE and cis-DCE degradation kinetics by zero-valent iron-experimental results and numerical simulation

The successful dechlorination of mixtures of chlorinated hydrocarbons with zero-valent metals requires information concerning the kinetics of simultaneous degradation of different contaminants. This includes intraspecies competitive effects (loading of the reactive iron surface by a single contaminant) as well as interspecies competition of several contaminants for the reactive sites available. In columns packed with zero-valent iron, the degradation behaviour of trichloroethylene (TCE), cis-dichloroethylene (DCE) and mixtures of both was measured in order to investigate interspecies competition. Although a decreasing rate of dechlorination is to be expected, when several degradable substances compete for the reactive sites on the iron surface, TCE degradation is nearly unaffected by the presence of cis-DCE. In contrast, cis-DCE degradation rates decrease significantly when TCE is added. A new modelling approach is developed in order to identify and quantify the observed competitive effects. The numerical model TBC (Transport, Biochemistry and Chemistry, Sch?fer et al., 1998a) is used to describe adsorption, desorption and dechlorination in a mechanistic way. Adsorption and degradation of a contaminant based on a limited number of reactive sites leads to a combined zero- and first-order degradation kinetics for high and low concentrations, respectively. The adsorption of several contaminants with different sorption parameters to a limited reactive surface causes interspecies competition. The reaction scheme and the parameters required are successfully transferred from Arnold and Roberts (2000b) to the model TBC. The degradation behaviour of the mixed contamination observed in the column experiments can be related to the adsorption properties of TCE and cis-DCE. By predicting the degradation of the single substances TCE and cis-DCE as well as mixtures of both, the calibrated model is used to investigate the effects of interspecies competition on the design of permeable reactive iron barriers. Even if TCE is present in only small concentrations (>3% of molar cis-DCE concentration) it is the contaminant limiting the residence time and the required thickness of the iron barrier.     

Sorption of pentachlorophenol on pine bark

The minimization of pentachlorophenol (PCP) transport in the environment driven by industrial wastewater discharges can be accomplished by sorption in natural, available and low cost by-products like pine bark. Taking into account that PCP is a chemical which behaviour is highly dominated by the surrounding features, this work intended to evaluate the sorption kinetics and equilibrium parameters according to the pH and temperature as well as the pine bark particle size. The PCP uptake by pine bark showed to be faster in the initial phase followed by a slower process, being 24 h the suitably time to reach the sorption equilibrium in the range of pH studied. The neutral PCP species showed to have higher binding capacity to pine bark than the anionic PCP, which was reflected in a decrease in the distribution coefficient (Kd) of the linear sorption isotherm with the increase of solution pH from 2 to 7. On the other hand, between 10 degrees C and 35 degrees C, the temperature does not seem to play a significant role in the PCP sorption by pine bark, while the sorbent size is a key parameter to enhance the overall process.     

Sorption of polycyclic aromatic hydrocarbons (PAHs) to carbonaceous materials in a river floodplain soil

We report on sorption isotherm of phenanthrene (Phe) for river floodplain soil associated with carbonaceous materials, with particular attention being devoted to the natural loading of Phe. Our sorption experiments with original soil samples, size, and density sub-fractions showed that the light fraction had the highest sorption capacity comparable to low rank coals. In addition, the light fraction contributed most for the sorption of Phe in total soil samples. K_oc values for all fractions were in the same range, thus indicating that coal and coal-derived particles in all samples are responsible for the enhanced sorption for Phe. Sorption was strongly nonlinear and the combined partitioning and pore-filling model gave a better fit than the Freundlich sorption model. In addition, the spiked PAHs did not show the same behavior as the naturally aged ones, therefore the accessibility of indigenous background organic contaminants was reduced when coal and coal-derived particles are associated with the soils.     

Prediction of some in situ tracer tests with sorbing tracers using independent data

Some recent converging tracer tests with sorbing tracers at the Asp? Hard Rock Laboratory in Sweden, the TRUE tests, have been predicted using only laboratory data and hydraulic data from borehole measurements. No model parameters were adjusted to obtain a better fit with the experiments. The independent data were fracture frequency and transmissivity data obtained in the field and laboratory data on sorption and matrix diffusion. Transmissivity measurements in five boreholes in the rock volume containing the region surrounding the injection and collection points show that there is a high frequency of water conducting fractures. Of 162 packed off sections with 0.5 m packer distances, 112 were found to have a transmissivity above the detection limit. The specific flow-wetted surface (FWS) of the rock mass could be estimated from these data. The transmissivities were found to be reasonably well described by a lognormal distribution. Laboratory data on diffusion and sorption properties together with the hydraulic data were used to "predict" the residence time distribution (RTD) of the sorbing tracers. The results were compared with the experimental breakthrough curves. In these experiments, the water residence time is very small compared to the residence time of the sorbing tracers due to their diffusion and sorption within the rock matrix. We thus could neglect the influence of the water residence time in our predictions. Therefore, no information on water residence times or on "dispersion" was needed. The dispersion of the sorbing tracers is caused by the different sorbing tracer residence times in different pathways. The sorbing tracer residence time is determined by the ratio of flowrate to the flow-wetted surface in the different pathways and not by the water residence time. Assuming a three-dimensional flow pattern and using the observed fracture frequency and flowrate distribution, breakthrough curves for three strongly sorbing tracers were predicted. Only the laboratory data, the transmissivity measurements and the pumping flowrate were used in the predictions. No information on the water residence time as obtained by the nonsorbing tracers was used. The predictions were surprisingly accurate.     

Influence of rate-limited sorption on the cleanup of layered soils by vapor extraction

The conditions under which rate-limited sorption is important for cleanup of layered soils by vapor extraction are investigated. The investigation includes two steps: (a) the cleanup time is estimated for a number of scenario cases by means of a numerical model and (b) the numerical results are approximated using analytical solutions derived for simplified models. In this way, equations are derived, which give insight into the influence of different parameters characterizing the properties of the soil, the geometry of the formation, the mass transfer mechanisms in it, and the distribution of the contaminant mass in the different phases (gas phase, water phase and solid phase). The numerical model used is based on the advection-dispersion differential equations for Darcian isothermal airflow, local equilibrium contaminant mass transfer between gas phase and soil water and first-order kinetics for mass transfer between soil water and solid phase. The numerical results are approximated combining an analytical solution to estimate cleanup time in layered formations for local equilibrium sorption, which has been presented in a previous work (J. Contam. Hydrol., 36 (1999) 105). with an analytical solution based on the well-mixed reservoir model under consideration of rate-limited sorption. The analytical approximation of the cleanup time is in reasonable agreement with the numerical results and allows its estimation with small computational effort.     

Organic compounds adsorption onto activated carbon: the effect of association between dissolved humic substances and pesticides

The quantitative determination of pesticide binding to dissolved humic substances is relevant to both water treatment operation using activated carbon adsorption process and the application of transport models that predict the environmental distribution patterns of a given hydrophobic contaminant. In this study and in a first set of experiments, the extent of binding between (i) three pesticides of environmental concern, aldicarb, lindane and pentachlorophenol, and (ii) dissolved commercial humic acid and soil extracted fulvic acid, was determined using dialysis experiments and water solubility enhancement tests. In a second set of experiments, the influence of dissolved humic substances or pesticide on the retention of the other co-adsorbate onto activated carbon was investigated in binary systems. It was found that association was negligible for aldicarb and that the pesticide sorption onto activated carbon was not affected by humic acid (8.5 mg liter(-1) DOC). The association constants K for lindane and pentachlorophenol were identical in the presence of fulvic acid (logK=4.1) but lower than that observed with humic acid. In the presence of humic acid, binding affinity for pentachlorophenol (logK=4.6) was higher than the one observed for lindane (logK=4.4), despite its much higher water solubility. This observation suggests that the aromatic character of the pentachlorophenol molecule contributes to association interactions with humic acid. From co-adsorption experiments onto activated carbon it was found that fulvic acid (7.7 mg litre(-1) DOC) slightly enhances sorption kinetics of pentachlorophenol. Lindane (1 mg litre(-1)) does not affect sorption kinetics for fulvic acid but markedly enhances both the sorption kinetics and adsorptive capacity for humic acid. Activated carbon retention of dissolved humic substances or pesticide appears to be enhanced by the association potential that exists between these co-adsorbates in some binary systems.     

Sorption kinetics and its effects on retention and leaching.

Sorption of pesticides to substrates used in biopurification systems is important as it controls the system's efficiency. Ideally, pesticide sorption should occur fast so that leaching of the pesticide in the biopurification system is minimized. Although modeling of pesticide transport commonly assumes equilibrium, this may not always be true in practice. Sorption kinetics have to be taken into account. This study investigated the batch sorption kinetics of linuron, isoproturon, metalaxyl, isoxaben and lenacil on substrates commonly used in a biopurification system, i.e. cow manure, straw, willow chopping, sandy loam soil, coconut chips, garden waste compost and peat mix. The first-order sorption kinetics model was fitted to the observed pesticide concentrations versus time resulting in an estimated kinetic rate constant alpha. Sorption appeared to be fast for the pesticides linuron and isoxaben, pesticides which were classified as immobile, while less mobile pesticides displayed an overall slower sorption. However, the substrate does not seem to be the main parameter influencing the sorption kinetics. Coconut chips, which is a substrate with a high organic matter content showed slow sorption for most of the pesticides. The effect of different estimated alpha values on the breakthrough of pesticides through a biopurification system was evaluated using the HYDRUS 1D model. Significant differences in leaching behavior were observed as a result of the obtained differences in sorption kinetics.     

Effect of clays, metal oxides, and organic matter on rhamnolipid biosurfactant sorption by soil

Rhamnolipids produced by Pseudomonas aeruginosa have been proposed as soil washing agents for enhanced removal of metal and organic contaminants from soil. A potential limitation for the application of rhamnolipids is sorption by soil matrix components. The objective of this study is to empirically determine the contribution of representative soil constituents (clays, metal oxides, and organic matter) to sorption of the rhamnolipid form most efficient at metal complexation (monorhamnolipid). Sorption studies show that monorhamnolipid (R1) sorption is concentration dependent. At low R1 concentrations that are relevant for enhancing organic contaminant biodegradation, R1 sorption followed the order: hematite (Fe(2)O(3))>kaolinite>MnO(2) approximately illite approximately Ca-montmorillonite>gibbsite (Al(OH)(3))>humic acid-coated silica. At high R1 concentrations, relevant for use in complexation/removal of metals or organics, R1 sorption followed the order: illite>humic acid-coated silica>Ca-montmorillonite>hematite>MnO(2)>gibbsite approximately kaolinite. These results allowed prediction of R1 sorption by a series of six soils. Finally, a comparison of R1 and R2 (dirhamnolipid) shows that the R1 form sorbs more strongly alone than when in a mixture of both the R1 and R2 forms. The information presented can be used to estimate, on an individual soil basis, the extent of rhamnolipid sorption. This is important for determining: (1) whether rhamnolipid addition is a feasible remediation option and (2) the amount of rhamnolipid required to efficiently remove the contaminant.     

Decontaminating soil organic pollutants with manufactured nanoparticles

Organic pollutants in soils might threaten the environmental and human health. Manufactured nanoparticles are capable to reduce this risk efficiently due to their relatively large capacity of sorption and degradation of organic pollutants. Stability, mobility, and reactivity of nanoparticles are prerequisites for their efficacy in soil remediation. On the basis of a brief introduction of these issues, this review provides a comprehensive summary of the application and effectiveness of various types of manufactured nanoparticles for removing organic pollutants from soil. The main categories of nanoparticles include iron (oxides), titanium dioxide, carbonaceous, palladium, and amphiphilic polymeric nanoparticles. Their advantages (e.g., unique properties and high sorption capacity) and disadvantages (e.g., high cost and low recovery) for soil remediation are discussed with respect to the characteristics of organic pollutants. The factors that influence the decontamination effects, such as properties, surfactants, solution chemistry, and soil organic matter, are addressed.     

Development of a continuously stirred flow cell for investigating sorption mass transfer

This paper introduces a new reversible-flow design for a continuously stirred reactor used to study sorption mass transfer in soil and solvent systems. The stirred reactor has potential advantages over conventional packed column or batch reactors because it isolates intraparticle sorption rate limitations from advective-dispersive transport, yet allows changes to flux through the reactor for analysis of sorption kinetics under dynamic conditions. Previously, stirred reactors have often failed due to clogging of sediment on the effluent frit. The reverse-flow backwashing design allows longer life and higher confidence in maintaining mixed conditions than previous designs. Mass transfer 'rate coefficients estimated from stirred and column experiments are compared; both techniques produced results consistent with a published correlation. The data also show that fitted sorption mass transfer coefficients can be strongly dependent on the choice of equilibrium partition coefficient (i.e. batch or first-moment derived values), and that the conventional two-site sorption kinetics model fails to accurately predict sorption mass transfer in the presence of changing solvent velocity through the reactor.