Laboratory migration experiments with radionuclides and natural colloids in a granite fracture

Natural colloids in groundwater could facilitate radionuclide transport, provided the colloids are mobile, are present in sufficient concentrations and can adsorb radionuclides. This paper describes the results of a laboratory migration study carried out with combinations of radionuclides and natural colloids within a fracture in a large granite block to experimentally determine the impact of colloids on radionuclide transport. The 85Sr used in this study is an example of a moderately sorbing radionuclide, while the 241Am is typical of a strongly sorbed radionuclide with very low solubility. The natural colloids used in this study were isolated from granite groundwater from Atomic Energy of Canada (AECL) Underground Research Laboratory (URL), and consisted of mostly 1-10 nm organic colloids, along with lesser amounts of 10-450 nm colloids (organics and aluminosilicates). The measured coefficients for radionuclide sorption onto these colloids were between 3 x 10(2) and 1 x 10(3) ml/g for 85Sr, and between 7 x 10(4) and 7 x 10(5) mg/l for 241Am. The 85Sr sorption on the natural colloids appeared to be reversible. Migration experiments in the granite block were carried out by establishing a flow field between two boreholes (out of a total of nine) intersecting a main horizontal fracture. These experiments showed that dissolved 85Sr behaved as a moderately sorbing tracer, while dissolved 241Am was completely adsorbed by the fracture surfaces and showed no evidence of transport. However, when natural colloids were injected together with dissolved 241Am, a small amount of 241Am transport was observed, demonstrating the ability of natural colloids to facilitate the transport of radionuclides with low solubility. Natural colloids had only a minor effect on the transport of 85Sr. In a separate experiment to test the effect of higher colloid concentrations on 85Sr migration, synthetic colloids were produced from Avonlea bentonite. The introduction of a relatively high concentration of bentonite colloids actually reduced 85Sr transport because, compared to natural colloids, the bentonite colloids were less mobile and they sorbed 85Sr more strongly.     

Anomalous transport of colloids and solutes in a shear zone

Transport experiments with colloids and radionuclides in a shear zone were conducted during the Colloid and Radionuclide Retardation experiment (CRR) at Nagra's Grimsel Test Site. Breakthrough curves of bentonite colloids and uranine, a non-sorbing solute, were measured in an asymmetric dipole flow field. The colloid breakthrough is earlier than that of uranine. Both breakthrough curves show anomalously long late time tails and the slope of the late time tails for the colloids is slightly higher. Anomalous late time tails are commonly associated with matrix diffusion processes; the diffusive interaction of solutes transported in open channels with the adjacent porous rock matrix or zones of stagnant water. The breakthrough curves for different colloid size classes are very similar and show no signs of fractionation due to their (size-dependent) diffusivity. It is proposed that tailing of the colloids is mainly caused by the structure of the flow field and that for the colloid transport, matrix diffusion is of minor importance. This has consequences for the interpretation of the uranine breakthrough. Comparisons of experimental results with numerical studies and with the evaluation of the colloid breakthrough with continuous time random theory imply that the tailing in the conservative solute breakthrough in this shear zone is not only caused by matrix diffusion. Part of the tailing can be attributed to advective transport in fracture networks and advection in low velocity regions. Models based on the advection-dispersion equation and matrix diffusion do not properly describe the temporal and spatial evolution of colloid and solute transport in such systems with a consistent set of parameters.     

Experimental determination of sorption in fractured flow systems

Fracture "skins" are alteration zones on fracture surfaces created by a variety of biological, chemical, and physical processes. Skins increase surface area, where sorption occurs, compared to the unaltered rock matrix. This study examines the sorption of organic solutes on altered fracture surfaces in an experimental fracture-flow apparatus. Fracture skins containing abundant metal oxides, clays, and organic material from the Breathitt Formation (Kentucky, USA) were collected in a manner such that skin surface integrity was maintained. The samples were reassembled in the lab in a flow-through apparatus that simulated approximately 2.7 m of a linear fracture "conduit." A dual-tracer injection scheme was utilized with the sorbing or reactive tracer compared to a non-reactive tracer (chloride) injected simultaneously. Sorption was assessed from the ratio of the first temporal moments of the breakthrough curves and from the loss of reactive tracer mass and evaluated as a function of flow velocity and solute type. The breakthrough curves suggest dual-flow regimes in the fracture with both sorbing and non-sorbing flow fields. Significant sorption occurs for the reactive components, and sorption increased with decreasing flow rate and decreasing compound solubility. Based on moment analysis, however, there was little retardation of the center of solute mass. These data suggest that non-equilibrium sorption processes dominate and that slow desorption and boundary layer diffusion cause extensive tailing in the breakthrough curves.     

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.     

Colloid-facilitated transport of radionuclides through fractured media

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

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.     

Applied tracer tests in fractured rock: Can we predict natural gradient solute transport more accurately than fracture and matrix parameters?

Applied tracer tests provide a means to estimate aquifer parameters in fractured rock. The traditional approach to analysing these tests has been using a single fracture model to find the parameter values that generate the best fit to the measured breakthrough curve. In many cases, the ultimate aim is to predict solute transport under the natural gradient. Usually, no confidence limits are placed on parameter values and the impact of parameter errors on predictions of solute transport is not discussed. The assumption inherent in this approach is that the parameters determined under forced conditions will enable prediction of solute transport under the natural gradient. This paper considers the parameter and prediction uncertainty that might arise from analysis of breakthrough curves obtained from forced gradient applied tracer tests. By adding noise to an exact solution for transport in a single fracture in a porous matrix we create multiple realisations of an initial breakthrough curve. A least squares fitting routine is used to obtain a fit to each realisation, yielding a range of parameter values rather than a single set of absolute values. The suite of parameters is then used to make predictions of solute transport under lower hydraulic gradients and the uncertainty of estimated parameters and subsequent predictions of solute transport is compared. The results of this study show that predictions of breakthrough curve characteristics (first inflection point time, peak arrival time and peak concentration) for groundwater flow speeds with orders of magnitude smaller than that at which a test is conducted can sometimes be determined even more accurately than the fracture and matrix parameters.     

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.     

Colloid and suspended particle migration experiments in a granite fracture

To determine the mobility of colloids (0.001–0.45 μm) and suspended particles (> 0.45 μm) in granite fractures, laboratory particle-migration and conservative tracer studies have been carried out in a natural fracture within a large granite block, with overall dimensions of 83×90×60 cm. Flow fields within this horizontal fracture were controlled through a set of 9 boreholes drilled orthogonally to the fracture. Laboratory experiments were performed using a range of average water velocities which contained values low enough to closely approximate the natural flow velocities of < 2 m yr⁻¹ in plutonic rocks of the Canadian Shield. The particles used had diameters between 0.02 and 22 μm, and included latex spheres, glass spheres and colloidal silica. Migration experiments were carried out with a filtered groundwater, ionic strength of 0.01 mol kg⁻¹, obtained from a granite fracture within the Whiteshell Research Area of Manitoba. Flushing experiments showed that suspended particles as large as 40 μm could be mobilized from the fracture surface. The mobility of suspended particles was significantly less than that of colloids. However, within the size range of colloids used in these studies (0.022–0.090 μm), colloid size did not affect colloid migration. Although, in general, colloids eluted ahead of the conservative tracer, colloid mobility was significantly reduced when the average groundwater velocity dropped below between 32 and 240 m yr⁻¹. Colloid transport was found to be very sensitive to flow path and flow direction.     

Transport of radionuclides in natural fractures: some aspects of laboratory migration experiments

Results are reported for a series of migration experiments performed in a hydraulically characterized, single, natural fracture in a block of granite with overall dimensions of 1 × 1 × 0.7 m (all approximate), using the conservative, poorly sorbing and strongly sorbing radionuclides ³H₂0, ¹³¹I, ²²Na, ⁸⁵Sr, ¹³⁷Cs, ⁶⁰Co, ^154,155Eu, ²³⁷Np, and ²³⁸Pu. The volumetric flow velocity of the transport solution was 3 ml h⁻¹, giving a residence time in the fracture of approximately 50 h. Elution profiles were obtained for ³H₂0, ¹³¹I, ²²Na, ⁸⁵Sr and ¹³⁷Cs but no evidence of the other radionuclides was observed in the eluent. Results from supporting static sorption measurements on crushed geological materials and granite coupons showed in general higher sorption on alteration minerals than on granite. Sorption was lowest for ²²Na and ⁸⁵Sr.The migration of ¹³¹I, ²²Na, ⁸⁵Sr through the fracture in real time was followed using end-window Geiger-Müller probes located in unused boreholes. Additional information, obtained by alpha and gamma scanning of the fracture surfaces after separating the block along the fracture, confirmed that transport had occurred along the flow path predicted from the hydraulic characterization of the fracture and that, over a 5.5 month period, the bulk of the injected ¹³⁷Cs had migrated only 75 cm along the flow path. The ⁶⁰Co, the rare earths and the actinides had not moved beyond the location of the injection borehole, suggesting that fracture infilling minerals played a major role in retarding radionuclide transport. Additional confirmation of the role of secondary minerals in radionuclide retardation was obtained using selective sequential extraction on the fracture surfaces. These observations support the inclusion of sorption data for fracture infilling minerals in the sorption database developed for the geosphere model for the Canadian Nuclear Fuel Waste Management Program.     

Transport modeling of copper and cadmium with linear and nonlinear retardation factors

The predictive accuracy of using the one-dimensional advection–dispersion equation to evaluate the fate and transport of solute in a soil column is usually dependent on the proper determination of chemical retardation factors. Typically, the distribution coefficient (K_d) obtained by fitting the linear sorption isotherm has been extensively used to consider general geochemical reactions on solute transport in a low-concentration range. However, the linear distribution coefficient cannot be adequately utilized to describe the solute fate at a higher concentration level. This study employed the nonlinear equilibrium-controlled sorption parameters to determine the retardation factor used in column leaching experiments. Copper and cadmium transportation in a lateritic silty-clay soil column was examined. Through the explicit finite-difference calculations with a third-order total-variation-diminishing (TVD) numerical solution scheme, all results of the theoretical copper and cadmium breakthrough curves (BTCs) simulated by using the Freundlich nonlinear retardation factors revealed good agreement with the experimental observations.     

Generation and stability of bentonite colloids at the bentonite/granite interface of a deep geological radioactive waste repository

The possible mechanisms of colloid generation at the near field/far field interface of a radioactive repository have been investigated by means of novel column experiments simulating the granite/bentonite boundary, both in dynamic and in quasi-static water flow conditions. It has been shown that solid particles and colloids can be detached from the bulk and mobilised by the water flow. The higher the flow rate, the higher the concentration of particles found in the water, according to an erosion process. However, the gel formation and the intrinsic tactoid structure of the clay play an important role in the submicron particle generation even in the compacted clay and in a confined system. In fact, once a bentonite gel is formed, in the regions where the clay is contacted with water, clay colloids can be formed even in quasi-static flow conditions. The potential relevance of these colloids in radionuclide transport has been studied by evaluating their stability in different chemical environments. The coagulation kinetics of natural bentonite colloids was experimentally studied as a function of the ionic strength and pH, by means of time-resolved light scattering techniques. It has been shown that these colloids are very stable in low saline (approximately 1 x 10(-3) M) and alkaline (pH > or = 8) waters.     

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.     

Modeling depth-variant and domain-specific sorption and biodegradation in dual-permeability media

A dual-permeability model (S_1D_DUAL) was developed to simulate the transport of land-applied pesticides in macroporous media. In this model, one flow domain was represented by the bulk matrix and the other by the preferential flow domain (PFD) where water and chemicals move at faster rates. The model assumed the validity of Darcian flow and the advective-dispersive solute transport in each of the two domains with inter-domain transfer of water and solutes due to pressure and concentration gradients. It was conceptualized that sorption and biodegradation rates vary with soil depth as well as in each of the two flow domains. In addition to equilibrium sorption, kinetic sorption was simulated in the PFD. Simulations were conducted to evaluate the combined effects of preferential flow, depth- and domain-variant sorption, and degradation on leaching of two pesticides: one with strong sorption potential (trifluralin) and the other with weak sorption potential (atrazine). Simulation results for a test case showed that water flux in the PFD was three times more than in the matrix for selected storm events. When equilibrium sorption was considered, the simulated profile of trifluralin in each domain was similar; however, the atrazine profile was deeper in the PFD than in the bulk matrix under episodic storm events. With an assumption of negligible sorption in the PFD, both the atrazine and the trifluralin profiles moved twice deeper into the PFD. The simulated concentrations of the chemicals were several orders higher in the PFD than in the matrix, even at deeper depths. The volume fraction of the macropores and the sorption and biodegradation properties of the chemicals could also affect the amount of pesticides leaving the root zone. For an intense storm event, slow sorption reaction rates in the PFD produced higher breakthrough concentrations of atrazine at the bottom of the simulated soil profile, thus posing the risk for breakthrough of chemicals from the root zone.     

Study of the contaminant transport into granite microfractures using nuclear ion beam techniques

Hydrated bentonite is a very plastic material and it is expected to enter in the rock microfractures at the granite/bentonite boundary of a deep geological high-level waste repository. This process is enhanced by the high swelling pressure of the clay. Since bentonite has a very good sorption capability for many radionuclides, the displacement of the clay might lead to a "clay-mediated" contaminant transport into the rock. The aim of this work is to study the contaminant transport into granite microfractures using nuclear ion beam techniques, and to determine to what extent the clay can favour it. To do so, bentonite previously doped with uranium, cesium and europium was put in contact with the surface of granite sheets. Granite sheets contacted with non-doped bentonite and with radionuclide solutions were also prepared as references. This allowed analysing the differences in the diffusion behaviour of the three systems: clay, radionuclides and clay plus radionuclides. A combination of Rutherford backscattering spectrometry (RBS) and other nuclear ion-beam techniques such as particle-induced X-ray emission (PIXE) and microPIXE was used to study the depth and lateral distribution of clay and contaminants inside granite. It was also tried to evaluate not only the diffusion depth and diffusion coefficients but also the different areas of the granite where the diffusants have a preferential access.     

Solute transport in crystalline rocks at Aspö--I: geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted Asp? Hard Rock Laboratory (SE Sweden). On a scale of decametres and larger, steeply dipping faults dominate and contain a variety of different fault rocks (mylonites, cataclasites, fault gouges). On a smaller scale, somewhat less regular fracture patterns were found. Conceptual models of the fault and fracture geometries and of the properties of rock types adjacent to fractures were derived and used as input for the modelling of in situ dipole tracer tests that were conducted in the framework of the Tracer Retention Understanding Experiment (TRUE-1) on a scale of metres. After the identification of all relevant transport and retardation processes, blind predictions of the breakthroughs of conservative to moderately sorbing tracers were calculated and then compared with the experimental data. This paper provides the geological basis and model calibration, while the predictive and inverse modelling work is the topic of the companion paper [J. Contam. Hydrol. 61 (2003) 175]. The TRUE-1 experimental volume is highly fractured and contains the same types of fault rocks and alterations as on the decametric scale. The experimental flow field was modelled on the basis of a 2D-streamtube formalism with an underlying homogeneous and isotropic transmissivity field. Tracer transport was modelled using the dual porosity medium approach, which is linked to the flow model by the flow porosity. Given the substantial pumping rates in the extraction borehole, the transport domain has a maximum width of a few centimetres only. It is concluded that both the uncertainty with regard to the length of individual fractures and the detailed geometry of the network along the flowpath between injection and extraction boreholes are not critical because flow is largely one-dimensional, whether through a single fracture or a network. Process identification and model calibration were based on a single uranine breakthrough (test PDT3), which clearly showed that matrix diffusion had to be included in the model even over the short experimental time scales, evidenced by a characteristic shape of the trailing edge of the breakthrough curve. Using the geological information and therefore considering limited matrix diffusion into a thin fault gouge horizon resulted in a good fit to the experiment. On the other hand, fresh granite was found not to interact noticeably with the tracers over the time scales of the experiments. While fracture-filling gouge materials are very efficient in retarding tracers over short periods of time (hours-days), their volume is very small and, with time progressing, retardation will be dominated by altered wall rock and, finally, by fresh granite. In such rocks, both porosity (and therefore the effective diffusion coefficient) and sorption K(d)s are more than one order of magnitude smaller compared to fault gouge, thus indicating that long-term retardation is expected to occur but to be less pronounced.     

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.     

Chromate transport through columns packed with surfactant-modified zeolite/zero valent iron pellets

Chromate transport through columns packed with zeolite/zero valent iron (Z/ZVI) pellets, either untreated or treated with the cationic surfactant hexadecyltrimethylammonium (HDTMA), was studied at different flow rates. In the presence of sorbed HDTMA, the chromate retardation factor increased by a factor of five and the pseudo first-order rate constant for chromate reduction increased by 1.5-5 times. The increase in rate constant from the column studies was comparable to a six-fold increase in the rate constant determined in a batch study. At a fast flow rate, the apparent delay in chromate breakthrough from the HDTMA modified Z/ZVI columns was primarily caused by the increase in chromate reduction rate constant. In contrast, at a slower flow rate, the retardation in chromate transport from the HDTMA modified Z/ZVI columns mainly originated from chromate sorption onto the HDTMA modified Z/ZVI pellets. Due to dual porosity, the presence of immobile water was responsible for the earlier breakthrough of chromate in columns packed with zeolite and Z/ZVI pellets. The results from this study further confirm the role of HDTMA in enhancing sorption and reduction efficiency of contaminants in groundwater remediation.     

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.     

Effects of water content on reactive transport of 85Sr in Chernobyl sand columns

It is known that under unsaturated conditions, the transport of solutes can deviate from ideal advective-dispersive behaviour even for macroscopically homogeneous porous materials. Causes may include physical non-equilibrium, sorption kinetics, non-linear sorption, and the irregular distribution of sorption sites. We have performed laboratory experiments designed to identify the processes responsible for the non-ideality of radioactive Sr transport observed under unsaturated flow conditions in an Aeolian sandy deposit from the Chernobyl exclusion zone. Miscible displacement experiments were carried out at various water contents and corresponding flow rates in a laboratory model system. Results of our experiments have shown that breakthrough curves of a conservative tracer exhibit a higher degree of asymmetry when the water content decreases than at saturated water content and same Darcy velocity. It is possible that velocity variations caused by heterogeneities at the macroscopic scale are responsible for this situation. Another explanation is that molecular diffusion drives the solute mass transfer between mobile and immobile water regions, but the surface of contact between these water regions is small. At very low concentrations, representative of a radioactive Sr contamination of the pore water, sorption and physical disequilibrium dominate the radioactive Sr transport under unsaturated flow conditions. A sorption reaction is described by a cation exchange mechanism calibrated under fully saturated conditions. The sorption capacity, as well as the exchange coefficients are not affected by desaturation. The number of accessible exchange sites was calculated on the basis that the solid remained in contact with water and that the fraction of solid phase in contact with mobile water is numerically equal to the proportion of mobile water to total water content. That means that for this type of sandy soil, the nature of mineral phases is the same in advective and non-advective domains. So sorption reaction parameters can be estimated from more easily conducted saturated experiments, but hydrodynamic behaviour must be characterized by conservative tracer experiments under unsaturated flow conditions.