Characterization of subsoil heterogeneity, estimation of grain size distribution and hydraulic conductivity at the Krauthausen test site using Cone Penetration Test

A Cone Penetration Test (CPT) survey with a high spatial resolution was performed in order to investigate the stratigraphy as well as the spatial variability of various soil properties of the Krauthausen test site. Analyses of the CPT measurements showed the subsurface to be dominated by a planar layered structure. Variogram analysis of the various CPT parameters disclosed that within each layer the soil properties have an anisotropic spatial correlation structure. A correlation analysis of the measured CPT data and co-located grain size distributions from soil samples was performed. Since the correlation coefficients were greater equal to 0.7, a reliable empirical relationship between the data sets could be developed. Based on this empirical relationship grain size distributions were estimated at CPT locations. The statistical processing of estimated and measured grain size distributions with respect to their spatial correlation structure disclosed good agreement between the data sets. The estimated grain size distributions from CPT data were used to estimate the hydraulic conductivity in the aquifer. The results provide detailed information of the spatial heterogeneity of the hydraulic conductivity at Krauthausen test site. The validation of these results, using a prior investigation of hydraulic conductivity statistics, suggests the CPT a fast and inexpensive tool for the estimation of three dimensional hydraulic conductivity fields with sufficient accuracy.     

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.     

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.     

Effects of source zone heterogeneity on surfactant-enhanced NAPL dissolution and resulting remediation end-points

The effectiveness of removal of nonaqueous phase liquids (NAPLs) from the entrapment source zone of the subsurface has been limited by soil heterogeneity and the inability to locate all entrapped sources. The goal of this study was to demonstrate the uncertainty of degree of source removal associated with aquifer heterogeneity. In this demonstration, source zone NAPL removal using surfactant-enhanced dissolution was considered. Model components that simulate the processes of natural dissolution in aqueous phase and surfactant-enhanced dissolution were incorporated into an existing code of contaminant transport. The dissolution modules of the simulator used previously developed Gilland-Sherwood type phenomenological models of NAPL dissolution to estimate mass transfer coefficients that are upscaleable to multidimensional flow conditions found at field sites. The model was used to simulate the mass removal from 10 NAPL entrapment zone configurations based on previously conducted two-dimensional tank experiments. These entrapment zones represent the NAPL distribution in spatially correlated random fields of aquifer hydraulic conductivity. The numerical simulations representing two-dimensional conditions show that effectiveness of mass removal depends on the aquifer heterogeneity that controls the NAPL entrapment and delivery of the surfactant to the locations of entrapped NAPLs. Flow bypassing resulting from heterogeneity and the reduction of relative permeability due to NAPL entrapment reduces the delivery efficiency of the surfactant, thus prolonging the remediation time to achieve desired end-point NAPL saturations and downstream dissolved concentrations. In some extreme cases, the injected surfactant completely bypassed the NAPL source zones. It was also found that mass depletion rates for different NAPL source configurations vary significantly. The study shows that heterogeneity result in uncertainties in the mass removal and achievable end-points that are directly related to dissolved contaminant plume development downstream of the NAPL entrapment zone.     

Impact of long-term wastewater irrigation on sorption and transport of atrazine in Mexican agricultural soils

In the Mezquital Valley, Mexico, crops have been irrigated with untreated municipal wastewater for more than a century. Atrazine has been applied to maize and alfalfa grown in the area for weed control for 15 years. Our objectives were to analyse (i) how wastewater irrigation affects the filtering of atrazine, and (ii) if the length of irrigation has a significant impact. We compared atrazine sorption to Phaeozems that have been irrigated with raw wastewater for 35 (P35) and 85 (P85) years with sorption to a non-irrigated (P0) Phaeozem soil under rainfed agriculture. The use of bromide as an inert water tracer in column experiments and the subsequent analysis of the tracers’ breakthrough curves allowed the calibration of the hydrodynamic parameters of a two-site non equilibrium convection-dispersion model. The quality of the irrigation water significantly altered the soils’ hydrodynamic properties (hydraulic conductivity, dispersivity and the size of pores that are hydraulically active). The impacts on soil chemical properties (total organic carbon content and pH) were not significant, while the sodium adsorption ratio was significantly increased. Sorption and desorption isotherms, determined in batch and column experiments, showed enhanced atrazine sorption and reduced and slower desorption in wastewater-irrigated soils. These effects increased with the length of irrigation. The intensified sorption-desorption hysteresis in wastewater-irrigated soils indicated that the soil organic matter developed in these soils had fewer high-energy, easily accessible sorption sites available, leading to lower and slower atrazine desorption rates. This study leads to the conclusion that wastewater irrigation decreases atrazine mobility in the Mezquital valley Phaeozems by decreasing the hydraulic conductivity and increasing the soil's sorption capacity.     

Impact of long-term wastewater irrigation on sorption and transport of atrazine in Mexican agricultural soils

In the Mezquital Valley, Mexico, crops have been irrigated with untreated municipal wastewater for more than a century. Atrazine has been applied to maize and alfalfa grown in the area for weed control for 15 years. Our objectives were to analyse (i) how wastewater irrigation affects the filtering of atrazine, and (ii) if the length of irrigation has a significant impact. We compared atrazine sorption to Phaeozems that have been irrigated with raw wastewater for 35 (P35) and 85 (P85) years with sorption to a non-irrigated (P0) Phaeozem soil under rainfed agriculture. The use of bromide as an inert water tracer in column experiments and the subsequent analysis of the tracers' breakthrough curves allowed the calibration of the hydrodynamic parameters of a two-site non equilibrium convection-dispersion model. The quality of the irrigation water significantly altered the soils' hydrodynamic properties (hydraulic conductivity, dispersivity and the size of pores that are hydraulically active). The impacts on soil chemical properties (total organic carbon content and pH) were not significant, while the sodium adsorption ratio was significantly increased. Sorption and desorption isotherms, determined in batch and column experiments, showed enhanced atrazine sorption and reduced and slower desorption in wastewater-irrigated soils. These effects increased with the length of irrigation. The intensified sorption-desorption hysteresis in wastewater-irrigated soils indicated that the soil organic matter developed in these soils had fewer high-energy, easily accessible sorption sites available, leading to lower and slower atrazine desorption rates. This study leads to the conclusion that wastewater irrigation decreases atrazine mobility in the Mezquital valley Phaeozems by decreasing the hydraulic conductivity and increasing the soil's sorption capacity.     

Description of hydrogeologic heterogeneity and evaluation of radionuclide transport at an underground nuclear test

Realistic models of contaminant transport in groundwater demand detailed characterization of the spatial distribution of subsurface hydraulic properties, while at the same time programmatic constraints may limit collection of pertinent hydraulic data. Fortunately, alternate forms of data can be used to improve characterization of spatial variability. We utilize a methodology that augments sparse hydraulic information (hard data) with more widely available hydrogeologic information to generate equiprobable maps of hydrogeologic properties that incorporate patterns of connected permeable zones. Geophysical and lithologic logs are used to identify hydrogeologic categories and to condition stochastic simulations using Sequential Indicator Simulation (SIS). The resulting maps are populated with hydraulic conductivity values using field data and Sequential Gaussian Simulation (SGS). Maps of subsurface hydrogeologic heterogeneity are generated for the purpose of examining groundwater flow and transport processes at the Faultless underground nuclear test, Central Nevada Test Area (CNTA), through large-scale, three-dimensional numerical modeling. The maps provide the basis for simulation of groundwater flow, while transport of radionuclides from the nuclear cavity is modeled using particle tracking methods. Sensitivity analyses focus on model parameters that are most likely to reduce the long travel times observed in the base case. The methods employed in this study have improved our understanding of the spatial distribution of preferential flowpaths at this site and provided the critical foundation on which to build models of groundwater flow and transport. The results emphasize that the impacts of uncertainty in hydraulic and chemical parameters are dependent on the radioactive decay of specific species, with rapid decay magnifying the effects of parameters that change travel time.     

Estimating the relevance of engineered carbonaceous nanoparticle facilitated transport of hydrophobic organic contaminants in porous media

Naturally occurring nanoparticles (NP) enhance the transport of hydrophobic organic contaminants (HOCs) in porous media. In addition, the debate on the environmental impact of engineered nanoparticles (ENP) has become increasingly important. HOC bind strongly to carbonaceous ENP. Thus, carbonaceous ENP may also act as carriers for contaminant transport and might be important when compared to existing transport processes. ENP bound transport is strongly linked to the sorption behavior, and other carbonaceous ENP-specific properties. In our analysis the HOC-ENP sorption mechanism, as well as ENP size and ENP residence time, was of major importance. Our results show that depending on ENP size, sorption kinetics and residence time in the system, the ENP bound transport can be estimated either as (1) negligible, (2) enhancing contaminant transport, or (3) should be assessed by reactive transport modeling. One major challenge to this field is the current lack of data for HOC-ENP desorption kinetics.     

Parameter and observation importance in modelling virus transport in saturated porous media-investigations in a homogenous system

This paper evaluates the importance of seven types of parameters to virus transport: hydraulic conductivity, porosity, dispersivity, sorption rate and distribution coefficient (representing physical-chemical filtration), and in-solution and adsorbed inactivation (representing virus inactivation). The first three parameters relate to subsurface transport in general while the last four, the sorption rate, distribution coefficient, and in-solution and adsorbed inactivation rates, represent the interaction of viruses with the porous medium and their ability to persist. The importance of four types of observations to estimate the virus-transport parameters are evaluated: hydraulic heads, flow, temporal moments of conservative-transport concentrations, and virus concentrations. The evaluations are conducted using one- and two-dimensional homogeneous simulations, designed from published field experiments, and recently developed sensitivity-analysis methods. Sensitivity to the transport-simulation time-step size is used to evaluate the importance of numerical solution difficulties. Results suggest that hydraulic conductivity, porosity, and sorption are most important to virus-transport predictions. Most observation types provide substantial information about hydraulic conductivity and porosity; only virus-concentration observations provide information about sorption and inactivation. The observations are not sufficient to estimate these important parameters uniquely. Even with all observation types, there is extreme parameter correlation between porosity and hydraulic conductivity and between the sorption rate and in-solution inactivation. Parameter estimation was accomplished by fixing values of porosity and in-solution inactivation.     

Contaminant degradation in physically and chemically heterogeneous aquifers

This paper examines the importance of the correlation between hydraulic conductivity (K) and degradation rate constant (k) during the transport of reactive contaminants in heterogeneous aquifers. We simulated reactive transport in an ensemble of two-dimensional heterogeneous aquifers. Two sets of transport simulations were conducted: one in which a perfect positive correlation was assumed between ln(K) and ln(k), and one in which a perfect negative correlation was assumed. We found that the sign of the correlation has important consequences for the contaminant transport. Qualitatively, a negative correlation leads to significantly more pronounced "fingering" of the contaminant plume than does a positive correlation, with potentially important consequences for downgradient receptors. Quantitatively, the expected behavior (as quantified by the contaminant mass remaining in the aquifer) is statistically different between the positive and negative cases: on average, more contaminant mass persists when K and k are negatively correlated. Also, the negative correlation leads to more variability between realizations of the ensemble, whereas a positive correlation induces relatively little variability between realizations. We discuss the implications of these findings for the management of contaminated aquifers.     

Heterogeneity of sorption and transport-related properties in a sand–gravel aquifer at Columbus, Mississippi

Properties related to sorption and transport of organic compounds have been determined on 126 sections of 17 cores taken in an aquifer at Columbus Air Force Base in Columbus, MS. Each core section was homogenized prior to analysis. Organic carbon content (OC), specific surface area (SA), distribution coefficient (K_d) for naphthalene, and particle size distribution were measured on each section. Hydraulic conductivity (K_h) for each section was calculated from the particle size distributions. K_h values obtained were comparable with those from earlier borehole flowmeter and pulse tracer tests. Frequency distributions for all properties were lognormal. The arithmetic means and standard deviations for all samples are: OC=0.028% (+0.031, −0.015), SA=4.02 m²/g (+3.95, −1.99), K_d=0.198 l/kg (+0.195, −0.098), K_h=0.00033 m/s (+0.00051, −0.00020). These standard deviations are asymmetrical about the mean because statistics were calculated using log-transformed data, and antilogarithms then taken to obtain the results in the units of property measurement. Variabilities, expressed as coefficients of variation, were similar for all properties. Correlations between the properties were investigated. A good correlation between naphthalene K_d and OC (r=0.78) was found, and other correlations were weak, thus indicating that organic carbon content may control sorption of nonpolar organic solutes in this low carbon aquifer. Autocorrelation (variogram) analysis indicated that, for all properties, correlation lengths were less than the distance between sample points, which were separated by about 20 m horizontally and 1 m vertically. Separate statistical analysis of two widely separated groups of wells showed the groups similar in all properties, except organic carbon. Large-scale inhomogeneity was not detected, although earlier tracer tests produced irregular plumes indicating inhomogeneity in observed solute transport. Implications of the results to site characterization, in situations where aquifers are heterogeneous on short length scales, are discussed.     

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

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

Evaluation of groundwater flow patterns around a dual-screened groundwater circulation well

Dual-screened groundwater circulation wells (GCWs) can be used to remove contaminant mass and to mix reagents in situ. GCWs are so named because they force water in a circular pattern between injection and extraction screens. The radial extent, flux and direction of the effective flow of this circulation cell are difficult to measure or predict. The objective of this study is to develop a robust protocol for assessing GCW performance. To accomplish this, groundwater flow patterns surrounding a GCW are assessed using a suite of tools and data, including: hydraulic head, in situ flow velocity, measured hydraulic conductivity data from core samples, chemical tracer tests, contaminant distribution data, and numerical flow and transport models. The hydraulic head data show patterns that are consistent with pumping on a dual-screened well, however, many of the observed changes are smaller than expected. In situ thermal perturbation flow sensors successfully measured horizontal flow, but vertical flow could not be determined with sufficient accuracy to be useful in mapping flow patterns. Two types of chemical tracer tests were utilized at the site and showed that much of the flow occurs within a few meters of the GCW. Flow patterns were also assessed based on changes in contaminant (trichloroethylene, TCE) concentrations over time. The TCE data clearly showed treated water moving away from the GCW at shallow and intermediate depths, but the circulation of that water back to the well, except very close to the well, was less clear. Detailed vertical and horizontal hydraulic conductivities were measured on 0.3 m-long sections from a continuous core from the GCW installation borehole. The measured vertical and horizontal hydraulic conductivity data were used to construct numerical flow and transport models, the results of which were compared to the head, velocity and concentration data. Taken together, the field data and modeling present a fairly consistent picture of flow and transport around the GCW. However, the time and expense associated with conducting all of those tests would be prohibitive for most sites. As a consequence, a sequential protocol for GCW characterization is presented here in which the number of tools used can be adjusted to meet the needs of individual sites. While not perfect, we believe that this approach represents the most efficient means for evaluating GCW performance.     

Tracer diffusion coefficients in sedimentary rocks: correlation to porosity and hydraulic conductivity

Matrix diffusion is an important transport process in geologic materials of low hydraulic conductivity. For predicting the fate and transport of contaminants, a detailed understanding of the diffusion processes in natural porous media is essential. In this study, diffusive tracer transport (iodide) was investigated in a variety of geologically different limestone and sandstone rocks. Porosity, structural and mineralogical composition, hydraulic conductivity, and other rock properties were determined. The effective diffusion coefficients were measured using the time-lag method. The results of the diffusion experiments indicate that there is a close relationship between total porosity and the effective diffusion coefficient of a rock (analogous to Archie's Law). Consequently, the tortousity factor can be expressed as a function of total porosity. The relationship fits best for thicker samples (> 1.0 cm) with high porosities (> 20%), because of the reduced influence of heterogeneity in larger samples. In general, these correlations appear to be a simple way to determine tortuosity and the effective diffusion coefficient from easy to determine rock porosity values.     

In situ modification of an aquifer material by a cationic surfactant to enhance retardation of organic contaminants

Sorption of hexadecyltrimethylammonium chloride (HDTMA), a cationic surfactant, on aquifer material from Columbus AFB, Mississippi, U.S.A., was examined. Transport studies using flow-through columns and a box model aquifer showed that an almost stationary zone of HDTMA-modified aquifer material could be produced in situ without a significant decrease in hydraulic conductivity.Perchloroethylene (PCE) and naphthalene sorption isotherms on the HDTMA-modified aquifer material were linear, and sorption coefficients were increased by over two orders of magnitude relative to the unmodified material. The retardation of PCE by insitu emplaced HDTMA zones within a column was examined. Agreement between batch- and column-derived sorption coefficients and breakthrough curve symmetry indicates that local equilibrium was attained. Significant retardation of a naphthalene plume by an in situ emplaced surfactant zone was demonstrated in the box model aquifer system.The experimental results indicate that it is feasible to create in situ a sorbent zone within an aquifer using cationic surfactants. In most situations, the sorbent zone concept needs to be coupled with contaminant degradation processes for sorbent emplacement to be a practical tool in the remediation of groundwater contamination sites. Sorbent zones may be of benefit in the engineering of suitable environments for microbial or abiotic degradation reactions and by providing time slow reactions to occur.     

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.     

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.     

Fate and mobility of pharmaceuticals in solid matrices

The sorption and mobility of six pharmaceuticals were investigated in two soil types with different organic carbon and clay content, and in bacterial biomass (aerobic and anaerobic). The pharmaceuticals examined were carbamazepine, propranolol, diclofenac sodium, clofibric acid, sulfamethoxazole and ofloxacin. The sorption experiments were performed according to the OECD test Guideline 106. The distribution coefficients determined by this batch equilibrium method varied with the pharmaceutical tested and the solid matrix type. Ofloxacin was particularly strongly adsorbed (except of the case of using anaerobic biomass for the solid matrix) while clofibric acid was found to be weakly adsorbed. The fate of pharmaceuticals in soil was also assessed using lysimeters. Important parameters that were studied were: the pharmaceutical loading rate and the hydraulic loading rate for adsorption and the rate and duration of a "rain" event for desorption. Major differences in the mobility of the six pharmaceuticals were observed and correlated with the adsorption/desorption properties of the compounds.     

Upscaling heterogeneity in aquifer reactivity via exposure-time concept: forward model

Reactive properties of aquifer solid phase materials play an important role in solute fate and transport in the natural subsurface on time scales ranging from years in contaminant remediation to millennia in dynamics of aqueous geochemistry. Quantitative tools for dealing with the impact of natural heterogeneity in solid phase reactivity on solute fate and transport are limited. Here we describe the use of a structural variable to keep track of solute flux exposure to reactive surfaces. With this approach, we develop a non-reactive tracer model that is useful for determining the signature of multi-scale reactive solid heterogeneity in terms of solute flux distributions at the field scale, given realizations of three-dimensional reactive site density fields. First, a governing Eulerian equation for the non-reactive tracer model is determined by an upscaling technique in which it is found that the exposure time of solution to reactive surface areas evolves via both a macroscopic velocity and a macroscopic dispersion in the artificial dimension of exposure time. Second, we focus on the Lagrangian approach in the context of a streamtube ensemble and demonstrate the use of the distribution of solute flux over the exposure time dimension in modeling two-dimensional transport of a solute undergoing simplified linear reversible reactions, in hypothetical conditions following prior laboratory experiments. The distribution of solute flux over exposure time in a given case is a signature of the impact of heterogeneous aquifer reactivity coupled with a particular physical heterogeneity, boundary conditions, and hydraulic gradient. Rigorous application of this approach in a simulation sense is limited here to linear kinetically controlled reactions.     

Mercury sorption to sediments: dependence on grain size, dissolved organic carbon, and suspended bacteria

A combination of laboratory scale derived correlations and measurements of grain size distribution, DOC (dissolved organic carbon) concentration, and density of suspended bacteria promises to be useful in estimating Hg(II) sorption in heterogeneous streambeds and groundwater environments. This was found by shaking intact sediment and fractions thereof (<63-2000mum) with solutions of HgCl(2) (1.0-10.0ngml(-1)). The intact sediment was also shaken with the Hg(II) solutions separately in presence of DOC (6.5-90.2mugml(-1)) or brought in contact with suspensions of a strain of groundwater bacteria (2x10(4)-2x10(6)cellsml(-1)). Hg(II) sorption was rather weak and positively correlated with the grain size, and the sorption coefficient (K(d)) varied between about 300 and 600mlg(-1). By using the relative surface areas of the fractions, K(d) for the intact sediment was back calculated with 2% deviation. K(d) was negatively correlated with the concentration of DOC and positively correlated with the number of bacteria. A multiple regression showed that K(d) was significantly more influenced by the number of bacteria than by the grain size. The findings imply that common DOC concentrations in groundwater and streambeds, 5-20mugml(-1), will halve the K(d) obtained from standard sorption assays of Hg(II), and that K(d) will almost double when the cell numbers are doubled at densities that are common in aquifers. The findings suggest that simultaneous measurements of surface areas of sediment particles, DOC concentrations, and bacterial numbers are useful to predict spatial variation of Hg(II) sorption in aquifers and sandy sediments.