Analyses of locally measured bromide breakthrough curves from a natural gradient tracer experiment at Krauthausen

Solute travel time distributions were derived from breakthrough curves (BTCs) of bromide concentrations, which were measured during a large-scale tracer experiment in a quaternary fluviatile aquifer at Krauthausen. Travel time distributions to a specific point in the aquifer were derived from locally measured BTCs, using averaged absolute concentrations ?abs(x1,t), normalized concentrations ?norm(x1,t), and velocity-weighted normalized concentrations ?vw(x1,t). The travel time distributions were characterized in terms of equivalent convective-dispersive transport parameters: the equivalent solute velocity and equivalent dispersivity. Parameters were derived from BTCs using moment analyses and least-squares fits of the 1-D convection-dispersion equation (CDE). Both local and averaged BTCs showed pronounced tailing which was not well described by the 1-D CDE and which indicates the presence of macroscopic regions with low velocities in the aquifer. Therefore, dispersivities derived from CDE fits were significantly smaller than those derived from time moments. The BTCs of ?abs(x1,t) were dominated by only a few local BTCs with high concentrations and were less representative for the travel time distribution than BTCs of averaged normalized concentrations. Dispersivities derived from ?norm(x1,t) and ?vw(x1,t) were very similar. Finally, estimates of dispersivities and vertical correlation length of lnK, gamma 3, from BTCs were in agreement with a first-order estimate of the dispersivity and gamma 3 based on grain size data and flow meter measurements.     

Estimating first-order reaction rate coefficient for transport with nonequilibrium linear mass transfer in heterogeneous media

A travel-time based approach is developed for estimating first-order reaction rate coefficients for transport with nonequilibrium linear mass transfer in heterogeneous media. Tracer transport in the mobile domain is characterized by a travel-time distribution, and mass transfer rates are described by a convolution product of concentrations in the mobile domain and a memory function rather than predefining the mass transfer model. A constant first-order reaction is assumed to occur only in the mobile domain. Analytical solutions in Laplace domain can be derived for both conservative and reactive breakthrough curves (BTCs). Temporal-moment analyses are presented by using the first and second moments of conservative and reactive BTCs and the mass consumption of the reactant for an inverse Gaussian travel-time distribution. In terms of moment matching, there is no need for one to specify the mass transfer model. With the same capacity ratio and the mean retention time, all mass transfer models will lead to the same moment-derived reaction rate coefficients. In addition, the consideration of mass transfer generally yields larger estimations of the reaction rate coefficient than models ignoring mass transfer. Furthermore, the capacity ratio and the mean retention time have opposite influences on the estimation of the reaction rate coefficient: the first-order reaction rate coefficient is positively linearly proportional to the capacity ratio, but negatively linearly proportional to the mean retention time.     

Transport of particulate material and dissolved tracer in a highly permeable porous medium: comparison of the transfer parameters

We are experimentally studying, by means of short-pulse injection, the transport and deposition kinetics of suspended particles (silts of the order of 10 microm) in a highly permeable medium consisting of a column of gravel. In our experiments, the breakthrough curves (BTCs) are well described by analytical solutions of a convection/dispersion model with first-order deposition kinetics. All the transport parameters calculated by the model for both particles and dissolved tracer depend on the flow rate. We demonstrate the existence of a critical flow rate, determined experimentally, beyond which the transfer time for the particles is longer than that for the tracer. This phenomenon is unusual in comparison with the results available in the literature. The increase in transfer time of particles in comparison to tracer leads us to assume a purely mechanical phenomenon, that is, collision between particles and grains of the medium with instantaneous reset in motion when the flow rate is sufficient to avoid settling. Thanks to the polydispersivity of the injected suspension and the control of grain size at the outlet, it can also be determined that the coarser particles are recovered before the finer particles, as expected when one considers the size-exclusion effect.     

Multi tracer test for the implementation of enhanced in-situ bioremediation at a BTEX-contaminated megasite

At the Centre for Environmental Research Leipzig-Halle (UFZ) research site in Zeitz, Germany, benzene contaminates the lower of two aquifers with concentrations of up to 20 mg/l. Since the benzene plume has a minimum length of approximately 1 km, enhanced natural attenuation measures are being considered as a remediation strategy. This study describes the performance and evaluation of a multi-species reactive tracer test using the tracers fluorescein and bromide as conservative tracers and toluene as reactive tracer. Sampling was performed over a period of six months using a detailed network of multilevel sampling wells. Toluene was only slightly retarded in comparison to bromide, whereas fluorescein was retarded considerably stronger. Therefore, it was not possible to use fluorescein as an in situ tracer for the determination of groundwater velocities. The ionic nature of fluorescein is assumed to be the major reason for its retardation. The results show that the infiltration conditions were suitable to produce a wide spreading of the tracer front along the full thickness of the aquifer. Thus, a large aquifer volume can be treated in future enhanced bioremediation measures. The total quantity of infiltrated toluene (24 l) was degraded under sulfate-reducing conditions over a flow path of 50 m. Benzylsuccinate was identified as a metabolite of toluene degradation under sulfate-reducing conditions at this site. The modelling results show that toluene degradation was described more accurately using Monod kinetics than first-order kinetics. Since toluene was only slightly retarded in comparison to bromide, sorption and desorption processes were considered to be negligible.     

Controlled release, blind test of DNAPL remediation by ethanol flushing

A dense nonaqueous phase liquid (DNAPL) source zone was established within a sheet-pile isolated cell through a controlled release of perchloroethylene (PCE) to evaluate DNAPL remediation by in-situ cosolvent flushing. Ethanol was used as the cosolvent, and the main remedial mechanism was enhanced dissolution based on the phase behavior of the water-ethanol-PCE system. Based on the knowledge of the actual PCE volume introduced into the cell, it was estimated that 83 L of PCE were present at the start of the test. Over a 40-day period, 64% of the PCE was removed by flushing the cell with an alcohol solution of approximately 70% ethanol and 30% water. High removal efficiencies at the end of the test indicated that more PCE could have been removed had it been possible to continue the demonstration. The ethanol solution extracted from the cell was recycled during the test using activated carbon and air stripping treatment. Both of these treatment processes were successful in removing PCE for recycling purposes, with minimal impact on the ethanol content in the treated fluids. Results from pre- and post-flushing partitioning tracer tests overestimated the treatment performance. However, both of these tracer tests missed significant amounts of the PCE present, likely due to inaccessibility of the PCE. The tracer results suggest that some PCE was inaccessible to the ethanol solution which led to the inefficient PCE removal rates observed. The flux-averaged aqueous PCE concentrations measured in the post-flushing tracer test were reduced by a factor of 3 to 4 in the extraction wells that showed the highest PCE removal compared to those concentrations in the pre-flushing tracer test.     

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.     

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

A novel inverse technique is proposed to quantitatively characterize macroscopic variability in aquifer reactivity in a Lagrangian representation. Reactivity heterogeneity is expressed in terms of distributions of flux over cumulative time of exposure of the solution to reactive surface area, termed here 'cumulative reactivity'. In cases involving single aqueous species the combined effects of physical and reactivity heterogeneity on reactive solute transport can often be established and further investigated through joint distributions of flux over travel time and cumulative reactivity. The inverse technique requires the breakthrough curve of a passive tracer to determine the distribution of flux over travel time, and additional breakthrough curves of reactive tracers provide additional moments of the distribution of flux over cumulative reactivity given travel time. Thus breakthroughs of one passive and two reactive tracers can provide the mean and variance of the distribution of flux over cumulative reactivity. This Lagrangian characterization is achieved with knowledge of the types of reactive surfaces present, but not their spatial locations. The distributions can subsequently be applied via forward modeling using the same technique to predict breakthrough curves of other solutes undergoing first-order reactions in similar physically and chemically heterogeneous configurations.     

Evaluating equilibrium and non-equilibrium transport of bromide and isoproturon in disturbed and undisturbed soil columns

In this study, displacement experiments of isoproturon were conducted in disturbed and undisturbed columns of a silty clay loam soil under similar rainfall intensities. Solute transport occurred under saturated conditions in the undisturbed soil and under unsaturated conditions in the sieved soil because of a greater bulk density of the compacted undisturbed soil compared to the sieved soil. The objective of this work was to determine transport characteristics of isoproturon relative to bromide tracer. Triplicate column experiments were performed with sieved (structure partially destroyed to simulate conventional tillage) and undisturbed (structure preserved) soils. Bromide experimental breakthrough curves were analyzed using convective-dispersive and dual-permeability (DP) models (HYDRUS-1D). Isoproturon breakthrough curves (BTCs) were analyzed using the DP model that considered either chemical equilibrium or non-equilibrium transport. The DP model described the bromide elution curves of the sieved soil columns well, whereas it overestimated the tailing of the bromide BTCs of the undisturbed soil columns. A higher degree of physical non-equilibrium was found in the undisturbed soil, where 56% of total water was contained in the slow-flow matrix, compared to 26% in the sieved soil. Isoproturon BTCs were best described in both sieved and undisturbed soil columns using the DP model combined with the chemical non-equilibrium. Higher degradation rates were obtained in the transport experiments than in batch studies, for both soils. This was likely caused by hysteresis in sorption of isoproturon. However, it cannot be ruled out that higher degradation rates were due, at least in part, to the adopted first-order model. Results showed that for similar rainfall intensity, physical and chemical non-equilibrium were greater in the saturated undisturbed soil than in the unsaturated sieved soil. Results also suggested faster transport of isoproturon in the undisturbed soil due to higher preferential flow and lower fraction of equilibrium sorption sites.     

Modeling in-situ uranium(VI) bioreduction by sulfate-reducing bacteria

We present a travel-time based reactive transport model to simulate an in-situ bioremediation experiment for demonstrating enhanced bioreduction of uranium(VI). The model considers aquatic equilibrium chemistry of uranium and other groundwater constituents, uranium sorption and precipitation, and the microbial reduction of nitrate, sulfate and U(VI). Kinetic sorption/desorption of U(VI) is characterized by mass transfer between stagnant micro-pores and mobile flow zones. The model describes the succession of terminal electron accepting processes and the growth and decay of sulfate-reducing bacteria, concurrent with the enzymatic reduction of aqueous U(VI) species. The effective U(VI) reduction rate and sorption site distributions are determined by fitting the model simulation to an in-situ experiment at Oak Ridge, TN. Results show that (1) the presence of nitrate inhibits U(VI) reduction at the site; (2) the fitted effective rate of in-situ U(VI) reduction is much smaller than the values reported for laboratory experiments; (3) U(VI) sorption/desorption, which affects U(VI) bioavailability at the site, is strongly controlled by kinetics; (4) both pH and bicarbonate concentration significantly influence the sorption/desorption of U(VI), which therefore cannot be characterized by empirical isotherms; and (5) calcium-uranyl-carbonate complexes significantly influence the model performance of U(VI) reduction.     

Soil attenuation of As(III,V) and Se(IV,VI) seepage potential at ash disposal facilities

Leachate from ash landfills is frequently enriched with As and Se but their off-site movement is not well understood. The attenuation potential of As and Se by soils surrounding selected landfills during leachate seepage was investigated in laboratory column studies using simulated ash leachate. As(III, V) and Se(IV, VI) concentrations as well as pH, flow rate, and a tracer were monitored in influent and effluent for up to 800 pore volumes followed by sequential desorption, extraction, and digestion of column segments. Column breakthrough curves (BTCs) were compared to predictions based on previously measured sorption isotherms. Early As(V) breakthrough and retarded As(III) breakthrough relative to predicted BTCs are indicative of oxidative transformation during seepage. For Se(VI), which exhibits linear sorption and the lowest sorption propensity, measured BTCs were predicted fairly well by equilibrium sorption isotherms, except for the early arrival of Se(IV) in one site soil, which in part, may be due to higher column pH values compared to batch isotherms. Most of the As and Se retained by soils during leaching was found to be strongly sorbed (60–90%) or irreversibly bound (10–40%) with <5% readily desorbable. Redox potential favoring transformation to the more sorptive valence states of As(V) and Se(IV) will invoke additional attenuation beyond equilibrium sorption-based predictions. With the exception of Se(IV) on one site soil, results indicate that attenuation by down-gradient soils of As and Se in ash landfill seepage will often be no less than what is predicted by equilibrium sorption capacity with further attenuation expected due to favorable redox transformation processes, thus mitigating contaminant plumes and associated risks.     

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.     

Using aliphatic alcohols as gaseous tracers in determination of water contents and air-water interfacial areas in unsaturated sands

A new type of gaseous tracer utilizing nontoxic aliphatic alcohols for the determination of water content and air-water interfacial area is tested on unsaturated sands of low water content. Alcohol vapors are generated at room temperature and passed through the experimental sand column. Breakthrough curves (BTCs) of these vapors are obtained by monitoring their effluent concentrations using GC-FID. The retardation factor with respect to each vapor transport process is obtained by optimizing BTCs data using the CXTFIT program in the reverse problem mode. The water content and the interfacial area are subsequently calculated from their retardation factors by both equilibrium and nonequilibrium transport models. Experimental results indicate that the pentanol tracer is feasible in the determination of water content at conditions when the degree of water saturation is low. In the determination of air-water interfacial area, decanol is selected due to its interfacial adsorption characteristics. By comparing to interfacial areas from theoretical predictions as well as other conventional tarcer methods, the ones determined from the decanol tracer tests are found to be close to the true interfacial areas when the water content is low.     

Numerical simulations of radon as an in situ partitioning tracer for quantifying NAPL contamination using push-pull tests

Presented here is a reanalysis of results previously presented by [Davis, B.M., Istok, J.D., Semprini, L., 2002. Push-pull partitioning tracer tests using radon-222 to quantify non-aqueous phase liquid contamination. J. Contam. Hydrol. 58, 129-146] of push-pull tests using radon as a naturally occurring partitioning tracer for evaluating NAPL contamination. In a push-pull test where radon-free water and bromide are injected, the presence of NAPL is manifested in greater dispersion of the radon breakthrough curve (BTC) relative to the bromide BTC during the extraction phase as a result of radon partitioning into the NAPL. Laboratory push-pull tests in a dense or DNAPL-contaminated physical aquifer model (PAM) indicated that the previously used modeling approach resulted in an overestimation of the DNAPL (trichloroethene) saturation (S(n)). The numerical simulations presented here investigated the influence of (1) initial radon concentrations, which vary as a function of S(n), and (2) heterogeneity in S(n) distribution within the radius of influence of the push-pull test. The simulations showed that these factors influence radon BTCs and resulting estimates of S(n). A revised method of interpreting radon BTCs is presented here, which takes into account initial radon concentrations and uses non-normalized radon BTCs. This revised method produces greater radon BTC sensitivity at small values of S(n) and was used to re-analyze the results from the PAM push-pull tests reported by Davis et al. The re-analysis resulted in a more accurate estimate of S(n) (1.8%) compared with the previously estimated value (7.4%). The revised method was then applied to results from a push-pull test conducted in a light or LNAPL-contaminated aquifer at a field site, resulting in a more accurate estimate of S(n) (4.1%) compared with a previously estimated value (13.6%). The revised method improves upon the efficacy of the radon push-pull test to estimate NAPL saturations. A limitation of the revised method is that 'background' radon concentrations from a non-contaminated well in the NAPL-contaminated aquifer are needed to accurately estimate NAPL saturation. The method has potential as a means of monitoring the progress of NAPL remediation.     

Water and solute movement in soil as influenced by macropore characteristics: 2. Macropore tortuosity

The paper describes the results of a laboratory study on the effects of macropore tortuosity on breakthrough curves BTCs and solute distribution in a Forman loam (fine loamy-mixed Udic Haploborolls) soil. BTC were obtained using 2-D columns (slab) containing artificial macropores of five different tortuosity levels. The BTCs were run under a constant hydraulic head of 0.08 m over an initially air dry soil. The input solutions contained 1190 mg l⁻¹ of potassium bromide, 10 mg l⁻¹ of Rhodamine WT, and 100 mg l⁻¹ of FD&C Blue #1. A soil column without macropores served as a control. The displacement of a non-adsorbed tracer was not affected by the tortuosity level. An increase in macropore tortuosity progressively increased the breakthrough time, increased the apparent retardation coefficient (R′), decreased the depth to the center of mass of a given adsorbed tracer, and increased the anisotropy in tracer distribution profile. The relative importance of macropore tortuosity increased with an increase in the adsorption coefficient of the tracer. Compared to macropore continuity, the macropore tortuosity had greater impact on solute distribution profile than in its leaching.     

Characterizing pesticide sorption and degradation in microscale biopurification systems using column displacement experiments

Biopurification systems treating pesticide contaminated water are very efficient, however they operate as a black box. Processes inside the system are not yet characterized. To optimize the performance, knowledge of degradation and retention processes needs to be generated. Therefore, displacement experiments were carried out for four pesticides (isoproturon, bentazone, metalaxyl, linuron) in columns containing different organic mixtures. Bromide, isoproturon and bentazone breakthrough curves (BTCs) were well described using the convection-dispersion equation (CDE) and a first-order degradation kinetic approach. Metalaxyl and linuron BTCs were well described using the CDE model expanded with Monod-type kinetics. Freundlich sorption, first-order degradation and Monod kinetics coefficients were fitted to the BTCs. Fitted values of the distribution coefficient K_f,column were much lower than those determined from batch experiments. Based on mobility, pesticides were ranked as: bentazone > metalaxyl - isoproturon > linuron. Based on degradability, pesticides were ranked as: linuron > metalaxyl - isoproturon > bentazone.     

Darcian preferential water flow and solute transport through bimodal porous systems: Experiments and modelling

Soils often exhibit a variety of small-scale heterogeneities such as inter-aggregate pores and voids which partition flow into separate regions. In this paper a methodological approach is discussed for characterizing the hydrological behaviour of a heterogeneous clayey–sandy soil in the presence of structural inter-aggregate pores. For the clay soil examined, it was demonstrated that, coupling the transfer function approach for analyzing BTCs and water retention data obtained with different methods from laboratory studies captures the bimodal geometry of the porous system along with the related existence of fast and slow flow paths. To be effectively and reliably applied this approach requires that the predominant effects of the soil hydrological behaviour near saturation be supported by accurate experimental data of both breakthrough curves (BTCs) and hydraulic functions for high water content values. This would allow the separation of flow phases and hence accurate identification of the processes and related parameters.     

Effects of degree of water saturation on dispersivity and immobile water in sandy soil columns

Three natural nonaggregated soil samples, with similar grain-size distributions, have been used to determine the dispersive behavior of porous media under steady, saturated and unsaturated flow conditions. Tritium was used as a tracer and was found to have no sorption on the solid matrix. Generated breakthrough curves (BTCs) for the unsaturated experiments were symmetrical with no evidence of tailing. The unsaturated experiments for two of the soils were adequately described by considering all the water in the pore volume as mobile. However, about 10% of the pore water, independent of the degree of saturation, was found to be immobile in the case of the third soil during unsaturated flow. For this soil, there was no mass transfer between the two water regions, indicating that the immobile water is essentially isolated from the flowing water fraction. For all three soils, dispersivity under unsaturated conditions was found to be higher, independent of the degree of water saturation, than the value determined for the saturated experiments. This is inconsistent with what would be expected from the simple bundle-of-capillary-tubes model and does not agree well with a more sophisticated conceptualization of the porous medium. The data, however, clearly indicate a wider range in pore-water velocities when these soils are desaturated.     

Microbial reduction of sulfate injected to gas condensate plumes in cold groundwater

Despite a rapid expansion over the past decade in the reliance on intrinsic bioremediation to remediate petroleum hydrocarbon plumes in groundwater, significant research gaps remain. Although it has been demonstrated that bacterial sulfate reduction can be a key electron accepting process in many petroleum plumes, little is known about the rate of this reduction process in plumes derived from crude oil and gas condensates at cold-climate sites (mean temperature <10 degrees C), and in complex hydrogeological settings such as silt/clay aquitards. In this field study, sulfate was injected into groundwater contaminated by gas condensate plumes at two petroleum sites in Alberta, Canada to enhance in-situ bioremediation. In both cases the groundwater near the water table had low temperature (6-9 degrees C). Monitoring data had provided strong evidence that bacterial sulfate reduction was a key terminal electron accepting process (TEAP) in the natural attenuation of dissolved hydrocarbons at these sites. At each site, water with approximately 2000 mg/L sulfate and a bromide tracer was injected into a low-sulfate zone within a condensate-contaminant plume. Monitoring data collected over several months yielded conservative estimates for sulfate reduction rates based on zero-order kinetics (4-6 mg/L per day) or first-order kinetics (0.003 and 0.01 day(-1)). These results favor the applicability of in-situ bioremediation techniques in this region, under natural conditions or with enhancement via sulfate injection.     

Use of tracer tests to investigate changes in flow and transport properties due to bioclogging of porous media

Tracer tests were conducted in three laboratory columns to study changes in the hydraulic properties of a porous medium due to bioclogging. About 30 breakthrough curves (BTCs) for each column were obtained. The BTCs were analyzed using analytical equilibrium and dual-porosity models, and estimates of the hydrodynamic dispersion and mass transfer coefficients were obtained by curve fitting. The change in transport properties developed in three stages: an initial phase (I) with no significant changes in transport properties, phase II with growth of biomass near the inlet of the columns causing changes in dispersivity, and phase III with added growth of micro-colonies deeper in the columns causing mass transfer of solutes from the water phase to the biophase. Tracer transport changed from being uniform to more non-uniform with increase in mass transfer of the tracer between the mobile phase and the immobile biomass. An increase in the bulk dispersivity value of up to one order of magnitude was observed. Numerical simulations suggest that local dispersivity values may be as much as 40 times higher in the more severe clogged areas inside the column. The bulk hydraulic conductivities of the columns decreased by up to three orders of magnitude. The hydraulic conductivity and dispersivity parameters were almost recovered after disinfection of the columns. Different models relating the changes of the hydraulic conductivity to the changes in the mobile porosity due to bioclogging were reviewed, and the micro-colony relation of Thullner et al. [Thullner, M., Zeyer, J., Kinzelbach, W., 2002. Influence of microbial growth on hydraulic properties of pore networks, Transport in Porous Media, 49, 99-122.] was found to best describe the relation between the bulk hydraulic parameters.