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.     

Reactive transport of 85Sr in a chernobyl sand column: static and dynamic experiments and modeling

The effects of nonlinear sorption and competition with major cations present in the soil solution on radioactive strontium transport in an eolian sand were examined. Three laboratory techniques were used to identify and quantify the chemical and hydrodynamic processes involved in strontium transport: batch experiments, stirred flow-through reactor experiments and saturated laboratory columns. The major goal was to compare the results obtained under static and dynamic conditions and to describe in a deterministic manner the predominant processes involved in radioactive strontium transport in such systems. Experiments under dynamic conditions, namely flow-through reactor and column experiments, were in very good agreement even though the solid/liquid ratio was very different. The experimental data obtained from the flow-through reactor study pointed to a nonlinear, instantaneous and reversible sorption process. Miscible displacement experiments were conducted to demonstrate the competition between stable and radioactive strontium and to quantify its effect on the 85Sr retardation factor. The results were modeled using the PHREEQC computer code. A suitable cation-exchange model was used to describe the solute/soil reaction. The model successfully described the results of the entire set of miscible displacement experiments using the same set of parameter values for the reaction calculations. The column study revealed that the stable Sr aqueous concentration was the most sensitive variable of the model, and that the initial state of the sand/solution system had also to be controlled to explain and describe the measured retardation factor of radioactive strontium. From these observations, propositions can be made to explain the discrepancies observed between some data obtained from static (batches) and dynamic (reactor and column) experiments. Desorbed antecedent species (stable Sr) are removed from the column or reactor in the flow system but continue to compete for sorption sites in the batch system. Batch experiments are simple and fast, and provide a very useful means of multiplying data. However, interpretation becomes difficult when different species compete for sorption sites in the soil/solution system. A combination of batches, flow-through reactor and column experiments, coupled with hydrogeochemical modeling, would seem to offer a very powerful tool for identifying and quantifying the predominant processes on a cubic decimeter scale (dm3) and for providing a range of radioactive strontium retardation factor as a function of the geochemistry of the soil/solution system.     

Adsorption and transport of arsenate in carbonate-rich soils: coupled effects of nonlinear and rate-limited sorption

The transport and fate of arsenate in carbonate-rich soil under alkaline conditions was investigated with multiple approaches combining batch, sequential extraction and column experiments as well as transport modeling studies. Batch experiments indicated that sorption isotherm was nonlinear over a wide range of concentration (0.1-200 mg L(-1)) examined. As(V) adsorption to the calcareous soil was initially fast but then continued at a slower rate, indicating the potential effect of rate-limited sorption on transport. Column experiments illustrated that transport of As(V) was significantly retarded compared to a non-reactive tracer. The degree of retardation decreased with increasing As(V) concentration. As(V) breakthrough curves exhibited nonideal transport behavior due to the coupled effects of nonlinear and rate-limited sorption on arsenate transport, which is consistent with the results of modeling studies. The contribution of nonlinear sorption to the arsenate retardation was negligible at low concentration but increased with increasing As(V) concentration. Sequential extraction results showed that nonspecifically sorbed (easily exchangeable, outer sphere complexes) fraction of arsenate is dominant with respect to the inner-sphere surface bound complexes of arsenate in the carbonate soil fraction, indicating high bioavailability and transport for arsenate in the carbonate-rich soils of which Fe and Al oxyhydroxide fractions are limited.     

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.     

Fate of seven pesticides in an aerobic aquifer studied in column experiments

The fate of selected pesticides (bentazone, isoproturon, DNOC, MCPP, dichlorprop and 2,4-D) and a metabolite (2,6-dichlorobenzamide (BAM)) was investigated under aerobic conditions in column experiments using aquifer material and low concentrations of pesticides (approximately 25 microg/l). A solute transport model accounting for kinetic sorption and degradation was used to estimate sorption and degradation parameters. Isoproturon and DNOC were significantly retarded by sorption, whereas the retardation of the phenoxy acids (MCPP, 2,4-D and dichlorprop), BAM and bentazone was very low. After lag periods of 16-33 days for the phenoxy acids and 80 days for DNOC, these pesticides were degraded quickly with 0.-order rate constants of 1.3-2.6 microg/l/day. None of the most probable degradation products were detected.     

Correlations of nonlinear sorption of organic solutes with soil/sediment physicochemical properties

The estimation of solute sorptive behaviors is essential when direct sorption data are unavailable and will provide a convenient way to assess the fate and the biological activity of organic solutes in soil/sediment environments. In this study, the sorption of 2,4-dichlorophenol (2,4-DCP) on 19 soil/sediment samples and the sorption of 13 organic solutes on one sediment were investigated. All sorption isotherms are nonlinear and can be described satisfactorily by a simple dual-mode model (DMM): q(e)=KpCe+Q0 . bCe/(1+bCe), where Kp (mlg(-1)) is the partition coefficient; Ce (microgml(-1)) is the equilibrium concentration; Q0 (microgg(-1)) is the maximum adsorption capacity; Q0 . b (mlg(-1)) is the Langmuir-type isotherm slope in the low concentration (Henry's law) range and b (mlmicrog(-1)) is a constant related to the affinity of the surface for the solute. Based on these nonlinear sorption isotherms and similar other nonlinear isotherms, it is observed that, for both polar 2,4-DCP and nonpolar phenanthrene, Kp, Q0 and Q0 . b are linearly correlated with soil/sediment organic carbon content (f(oc) in the range of 0.118-53.7%). The results indicate that the nonlinear sorption of organic solutes results primarily from interactions with soil/sediment organic matter. The K*oc K*oc=Kp/f(oc)), Qoc (Qoc=Q0/f(oc)), Loc (Loc=Q0 . b/f(oc)) and b for a given organic solute with different soils/sediments are largely invariant. Furthermore, logK*oc, logb and logLoc for various organic solutes are correlated significantly with the solute logKow or logSw (logKow in the range of 0.9 to 5.13 and logSw in the range of -6.176 to -0.070). A fundamental empirical equation was then established to calculate approximately the nonlinear sorption from soil/sediment f(oc) and solute Sw for a given solute equilibrium concentration.     

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.     

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.     

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.     

Assessment of sorbent/water ratio effect on adsorption using dimensional analysis and batch experiments

Many factors affect adsorption phenomena in solid-liquid systems. One of the most important factors is the sorbent/water (S/W) ratio in the system. However, the effect of varying S/W ratios on the adsorption is still unclear. In this study, batch experiments were examined to observe the adsorption of four contaminants (copper, cadmium, Butachlor, and Deltamethrin) in six soils with texture ranging from silty clay to loamy sand and with different S/W ratios. Dimensional analysis was used to assess the relationship between adsorption phenomena and S/W ratio. We have assumed that the total amount of sorbate sorbed in soil is a function of the equilibrium concentration, the volume of sorbate solution, and the sorbent amount in the system. A power function (Freundlich-like) model was obtained from the dimensional analysis. It can describe precisely the adsorption phenomena of different sorbents and sorbates in the moisture regime of paddy soils. Therefore, proper adsorption parameters can be obtained by this power function model regardless of the solids effect, which can then be utilized to describe the fate of solute in soil using solute transport models.     

Determination of bromacil transport as a function of water and carbon content in soils

This study was conducted to determine the significance of bromacil transport as a function of water and carbon content in soils and to explore the implications of neglecting sorption when making assessments of travel time of bromacil through the vadose zone. Equilibrium batch sorption tests were performed for loamy sand and sandy soil added with four different levels of powdered activated carbon (PAC) content (0, 0.01, 0.05, and 0.1%). Column experiments were also conducted at various water and carbon contents under steady-state flow conditions. The first set of column experiments was conducted in loamy sand containing 1.5% organic carbon under three different water contents (0.23, 0.32, and 0.41) to measure breakthrough curves (BTCs) of bromide and bromacil injected as a square pulse. In the second set of column experiments, BTCs of bromide and bromacil injected as a front were measured in saturated sandy columns at the four different PAC levels given above. Column breakthrough data were analyzed with both equilibrium and nonequilibrium (two-site) convection-dispersion equation (CDE) models to determine transport and sorption parameters under various water and carbon contents. Analysis with batch data indicated that neglect of the partition-related term in the calculation of solute velocity may lead to erroneous estimation of travel time of bromacil, i.e. an overestimation of the solute velocity by a factor of R. The column experiments showed that arrival time of the bromacil peak was larger than that of the bromide peak in soils, indicating that transport of bromacil was retarded relative to bromide in the observed conditions. Extent of bromacil retardation (R) increased with decreasing water content and increasing PAC content, supporting the importance of retardation in the estimation of travel time of bromacil even at small amounts of organic carbon for soils with lower water content.     

Sorption isotherms of brominated diphenyl ethers on natural soils with different organic carbon fractions

Sorption isotherms of BDE-28 and BDE-47 on natural soils with different contents of soil organic matter (SOM) were investigated. Due to low water solubility of BDEs and resulted narrow ranges of aqueous equilibrium concentration, the linear distribution model showed similar and good fitting efficiency to the linear portion of nonlinear Freundlich curve. For the same sample, the linear and nonlinear model fitting sorption coefficients were close. At the statistically significant level of 0.05 or 0.1, significant relationships of total organic carbon fraction (fOC) with the fitting sorption coefficients can be observed. As for BDE-28, the relationships of fOC and SOM fractions with the single point partition coefficients at different aqueous concentrations of BDEs were significant; while for BDE-47, the relationships became less significant or insignificant, especially at higher aqueous concentrations. The findings in this study may facilitate more understanding on transport and fate of studied BDEs in soil systems.     

A controlled field experiment on groundwater contamination by a multicomponent DNAPL: dissolved-plume retardation

A natural gradient emplaced-source (ES) controlled field experiment was conducted at the Borden aquifer research site, Ontario, to study the transport of dissolved plumes emanating from residual dense nonaqueous-phase liquid (DNAPL) source zones. The specific objective of the work presented here is to determine the effects of solute and co-solute concentrations on sorption and retardation of dissolved chlorinated solvent-contaminant plumes. The ES field experiment comprised a controlled emplacement of a residual multicomponent DNAPL below the groundwater table and intensive monitoring of dissolved-phase plumes of trichloromethane (TCM), trichloroethylene (TCE), and perchloroethylene (PCE) plumes continuously generated in the aquifer down gradient from gradual source dissolution. Estimates of plume retardation (and dispersion) were obtained from 3-D numerical simulations that incorporated transient source input and flow regimes monitored during the test. PCE, the most retarded solute, surprisingly exhibited a retardation factor approximately 3 times lower than observed in a previous Borden tracer test by Mackay et al. [Water Resour. Res. 22 (1986) 2017] conducted approximately 150 m away. Also, an absence of temporal trend in PCE retardation contrasted with the previous Borden test. Supporting laboratory studies on ES site core indicated that sorption was nonlinear and competitive, i.e. reduced sorption of PCE was observed in the presence of TCE. Consideration of the effects of relatively high co-solute (TCE) concentration (competitive sorption) in addition to PCE concentration effects (nonlinear sorption) was necessary to yield laboratory-based PCE retardation estimates consistent with the field plume values. Concentration- and co-solute-based sorption and retardation analysis was also applied to the previous low-concentration pulse injection test of Mackay et al. [Water Resour. Res. 22 (1986) 2017] and was able to successfully predict the temporal field retardation trends observed in that test. While it is acknowledged that other "nonideal transport" effects may contribute, our analysis predicts differences in the PCE retardation magnitude and trend between the two experiments that are consistent with field observations based on the marked solute concentration differences that resulted from contrasting source conditions. Solute and co-solute concentration effects have heretofore received little attention, but may have wide significance in aquifers contaminated by point-source pollutants because many plumes contain mixed solutes over wide concentration ranges in strata that are likely subject to nonlinear sorption.     

Retention behavior of hydrophobic organic chemicals as a function of temperature in soil leaching column chromatography

To study the transport mechanism of hydrophobic organic chemicals (HOCs) and the energy change in soil/solvent system, a soil leaching column chromatographic (SLCC) experiment at an environmental temperature range of 20-40 degrees C was carried out, which utilized a reference soil (SP 14696) packed column and a methanol-water (1:4 by volume ratio) eluent. The transport process quickens with the increase of column temperature. The ratio of retention factors at 30 and 40 degrees C (k'30/k'40) ranged from 1.08 to 1.36. The lower enthalpy change of the solute transfer in SLCC (from eluent to soil) than in conventional reversed-phase liquid chromatography (e.g., from eluent to C18) is consistent with the hypothesis that HOCs were dominantly and physically partitioned between solvent and soil. The results were also verified by the linear solvation energy relationships analysis. The chief factor controlling the retention was found to be the solute solvophobic partition, and the second important factor was the solute hydrogen-bond basicity, while the least important factors were the solute polarizability-dipolarity and hydrogen-bond acidity. With the increase of temperature, the contributions of the solute solvophobic partition and hydrogen-bond basicity gradually decrease, and the latter decreases faster than the former.     

Sorption of acidic organic solute onto kaolinitic soils from methanol-water mixtures

The fate of the acidic organic solute from the soil-water-solvent system is not well-understood. In this study, the effect of the acidic functional group of organic solute in the sorption from cosolvent system was evaluated. The sorption of naphthalene (NAP) and 1-naphthoic acid (1-NAPA) by three kaolinitic soils and two model sorbents (kaolinite and humic acid) were measured as functions of the methanol volume fractions (f (c) ≤ 0.4) and ionic compositions (CaCl(2) and KCl). The solubility of 1-NAPA was also measured in various ionic compositions. The sorption data were interpreted using the cosolvency-induced sorption model. The K (m) values (= the linear sorption coefficient) of NAP with kaolinitic soil for both ionic compositions was log linearly decreased with f (c). However, the K (m) values of 1-NAPA with both ionic compositions remained relatively constant over the f (c) range. For the model sorbent, the K (m) values of 1-NAPA with kaolinite for the KCl system and with humic acid for both ionic compositions decreased with f (c), while the sorption of 1-NAPA with kaolinite for the CaCl(2) system was increased with f (c). From the solubility data of 1-NAPA with f (c), no significant difference was observed with the different ionic compositions, indicating an insignificant change in the aqueous activity of the liquid phase. In conclusion, the enhanced 1-NAPA sorption, greater than that predicted from the cosolvency-induced model, was due to an untraceable interaction between the carboxylate and hydrophilic soil domain in the methanol-water system. Therefore, in order to accurately predict the environmental fate of acidic pesticides and organic solutes, an effort to quantitatively incorporate the enhanced hydrophilic sorption into the current cosolvency-induced sorption model is required.     

Sorption of imidazolium-based ionic liquids to aquatic sediments

Ionic liquids (ILs) have received much attention as "green" alternatives to traditional solvents because they do not evaporate, eliminating concerns over fugitive emissions. However, if ionic liquids are used in industrial applications, they may enter aquatic systems via effluent, and their fate and transport may be influenced by sorption to sediments. In this study, we conducted batch mixing experiments with four alkylmethylimidizolium-based ILs and four types of aquatic sediments to asses the capacity for natural aquatic sediments to remove these chemicals from the water column. The concentration isotherms were non linear with point estimates of the distribution coefficient (K(d)) decreasing with increasing concentration. Apparent distribution coefficients ranged from 7.9 to 95.7l kg(-1) at an initial concentration of 0.5mM and were positively related to sediment organic matter (SOM) content. These K(d) values indicate that the ILs did not sorb strongly to the tested sediments. Increased alkyl chain length did not lead to increased sorption suggesting that hydrophobic interactions were not the most important sorption mechanism. We conclude that aquatic sediments have a limited capacity to sorb alkylmethylimidazolium ILs and that the transport of these contaminants in aquatic systems will not be strongly attenuated by sediments.     

Column studies on transport of deicing additive benzotriazole in a sandy aquifer and a zerovalent iron barrier

Benzotriazole (BTA), a chemical with wide industrial applications, is a typical additive in deicer/anti-icer used at airport. To achieve a better understanding of the transport behaviour and environmental fate of BTA, laboratory column studies have been performed on subsoil samples from Oslo Airport, Gardermoen. To explore possibilities for aquifer remediation, BTA behaviour was also studied in a column of granular zerovalent iron (Fe(0)). The subsoil column study demonstrates a very limited retardation of BTA. Consecutive loadings of BTA of the subsoil column showed no change of the break-through curve (BTC) and complete desorption was observed. The sorption behaviour of BTA to metallic iron (Fe(0)) was rather complex. Considerable retardation was observed in the Fe(0) column and repeated BTA loading resulted in an earlier break-through. Between 20% and 50% of the input concentration was retained permanently in the iron (Fe(0)) column. The BTA sorption to metallic iron was found to be enhanced by chloride which lowered the break-through concentration (i.e the C/C(0) plateau). The fraction of BTA remaining in the iron column was found to vary with the flow rate, indicating a time dependant multilayer sorption mechanism. The steady increase in the amount of adsorbed BTA to the iron column during loading corresponds to a rather strong bonding of 4-15 BTA layers to the iron surface. A very slow desorption of BTA was observed; even after flushing with 753 pore volumes of BTA free water, 7.5% of the BTA remained in the column. A geochemical model was developed based on PHREEQC-2 to simulate the sorption and transport of BTA in the tested materials. The BTA sorption was modelled with Freundlich sorption isotherms, as earlier determined in batch experiments. A slight adjustment of the Freundlich parameters was required to fit the observed column break-through. However, our model was not able to simulate the long-term retainment of BTA in the granular iron columns. The simulations confirm the high mobility of BTA in groundwater aquifers and suggest that zerovalent iron could be used to retain a BTA plume, although oxidation of the sorbent might reduce the long-term performance of such a remediation scheme and slow desorption has been observed.     

Development of a simple, accurate SPME-based method for assay of VOCs in column breakthrough experiments

Solid-phase microextraction (SPME) with gas chromatography is to be used for assay of effluent liquid samples from soil column experiments associated with VOC fate/transport studies. One goal of the fate/transport studies is to develop accurate, highly reproducible column breakthrough curves for 1,2-cis-dichloroethylene (cis-DCE) and trichloroethylene (TCE) to better understand interactions with selected natural solid phases. For SPME, the influences of the sample equilibration time, extraction temperature and the ratio of volume of sample bottle to that of the liquid sample (V(T)/V(w)) are the critical factors that could influence accuracy and precision of the measured results. Equilibrium between the gas phase and liquid phase was attained after 200 min of equilibration time. The temperature must be carefully controlled due to variation of both the Henry's constant (K(h)) and the fibre/gas phase distribution coefficient (K(fg)). K(h) decreases with decreasing temperature while K(fg) increases. Low V(T)/V(w) yields better sensitivity but results in analyte losses and negative bias of the resultant assay. High V(T)/V(w) ratio yields reduced sensitivity but analyte losses were found to be minimal, leading to better accuracy and reproducibility. A fast SPME method was achieved, 5 min for SPME extraction and 3.10 min for GC analysis. A linear calibration function in the gas phase was developed to analyse the breakthrough curve data, linear between a range of 0.9-236 microgl(-1), and a detection limit lower than 5 microgl(-1).     

Application of continuous time random walk theory to nonequilibrium transport in soil

Continuous time random walk (CTRW) formulations have been demonstrated to provide a general and effective approach that quantifies the behavior of solute transport in heterogeneous media in field, laboratory, and numerical experiments. In this paper we first apply the CTRW approach to describe the sorbing solute transport in soils under chemical (or) and physical nonequilibrium conditions by curve-fitting. Results show that the theoretical solutions are in a good agreement with the experimental measurements. In case that CTRW parameters cannot be determined directly or easily, an alternative method is then proposed for estimating such parameters independently of the breakthrough curve data to be simulated. We conduct numerical experiments with artificial data sets generated by the HYDRUS-1D model for a wide range of pore water velocities (υ) and retardation factors (R) to investigate the relationship between CTRW parameters for a sorbing solute and these two quantities (υ, R) that can be directly measured in independent experiments. A series of best-fitting regression equations are then developed from the artificial data sets, which can be easily used as an estimation or prediction model to assess the transport of sorbing solutes under steady flow conditions through soil. Several literature data sets of pesticides are used to validate these relationships. The results show reasonable performance in most cases, thus indicating that our method could provide an alternative way to effectively predict sorbing solute transport in soils. While the regression relationships presented are obtained under certain flow and sorption conditions, the methodology of our study is general and may be extended to predict solute transport in soils under different flow and sorption conditions.