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.     

Transport of Sr2+ and SrEDTA2- in partially-saturated and heterogeneous sediments

Strontium-90 has migrated deep into the unsaturated subsurface beneath leaking storage tanks in the Waste Management Areas (WMA) at the U.S. Department of Energy's (DOE) Hanford Reservation. Faster than expected transport of contaminants in the vadose zone is typically attributed to either physical hydrologic processes such as development of preferential flow pathways, or to geochemical processes such as the formation of stable, anionic complexes with organic chelates, e.g., ethylenediaminetetraacetic acid (EDTA). The goal of this paper is to determine whether hydrological processes in the Hanford sediments can influence the geochemistry of the system and hence control transport of Sr(2+) and SrEDTA(2-). The study used batch isotherms, saturated packed column experiments, and an unsaturated transport experiment in an undisturbed core. Isotherms and repacked column experiments suggested that the SrEDTA(2-) complex was unstable in the presence of Hanford sediments, resulting in dissociation and transport of Sr(2+) as a divalent cation. A decrease in sorption with increasing solid:solution ratio for Sr(2+) and SrEDTA(2-) suggested mineral dissolution resulted in competition for sorption sites and the formation of stable aqueous complexes. This was confirmed by detection of MgEDTA(2-), MnEDTA(2-), PbEDTA(2-), and unidentified Sr and Ca complexes. Displacement of Sr(2+) through a partially-saturated undisturbed core resulted in less retardation and more irreversible sorption than was observed in the saturated repacked columns, and model results suggested a significant reservoir (49%) of immobile water was present during transport through the heterogeneous layered sediments. The undisturbed core was subsequently disassembled along distinct bedding planes and subjected to sequential extractions. Strontium was unequally distributed between carbonates (49%), ion exchange sites (37%), and the oxide (14%) fraction. An inverse relationship between mass wetness and Sr suggested that sandy sediments of low water content constituted the immobile flow regime. Our results suggested that the sequestration of Sr(2+) in partially-saturated, heterogeneous sediments was most likely due to the formation of immobile water in drier regions having low hydraulic conductivities.     

Mobility of pharmaceuticals carbamazepine, diclofenac, ibuprofen, and propyphenazone in miscible-displacement experiments

Many pharmaceuticals pass the unsaturated zone before reaching an aquifer. Therefore, laboratory sand column transport experiments were conducted to study the transport behavior of carbamazepine, diclofenac, ibuprofen, and propyphenazone under unsaturated conditions. The test water was artificial sewage effluent to simulate the infiltration of reused wastewater. The test water was spiked with the pharmaceutically active compounds and the tracer LiCl. Afterwards it was passed through laboratory sand columns, one experiment for each pharmaceutical. The physical and chemical parameters were recorded and general ions measured. Pharmaceuticals were measured using solid phase extraction, derivatization, and detection with GC-MS. The column experiments indicate a significant elimination of ibuprofen (54%), propyphenazone (55%), and diclofenac (35%), whereas carbamazepine was not eliminated. Retardation factors varied between 1.84 for carbamazepine, 2.51 for propyphenazone, 3.00 for ibuprofen, and 4.80 for diclofenac. These results show that mobility and elimination of diclofenac, ibuprofen, and propyphenazone is about in the same range as for experiments under saturated conditions whereas carbamazepine had a significantly lower sorption and elimination under unsaturated conditions.     

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.     

Modelling the migration of contaminants through variably saturated dual-porosity, dual-permeability chalk

In the Hesbaye region in Belgium, tracer tests performed in variably saturated fissured chalk rocks presented very contrasting results in terms of transit times, according to artificially controlled water recharge conditions prevailing during the experiments. Under intense recharge conditions, tracers migrated across the partially or fully saturated fissure network, at high velocity in accordance with the high hydraulic conductivity and low effective porosity (fracture porosity). At the same time, a portion of the tracer was temporarily retarded in the almost immobile water located in the matrix. Under natural infiltration conditions, the fissure network remained inactive. Tracers migrated downward through the matrix, at low velocity in relation with the low hydraulic conductivity and the large porosity of the matrix. Based on these observations, Brouyère et al. (2004a) [Brouyère, S., Dassargues, A., Hallet, V., 2004a. Migration of contaminants through the unsaturated zone overlying the Hesbaye chalky aquifer in Belgium: a field investigation, J. Contam. Hydrol., 72 (1-4), 135-164, doi: 10.1016/j.conhyd.2003.10.009] proposed a conceptual model in order to explain the migration of solutes in variably saturated, dual-porosity, dual-permeability chalk. Here, mathematical and numerical modelling of tracer and contaminant migration in variably saturated fissured chalk is presented, considering the aforementioned conceptual model. A new mathematical formulation is proposed to represent the unsaturated properties of the fissured chalk in a more dynamic and appropriate way. At the same time, the rock water content is partitioned between mobile and immobile water phases, as a function of the water saturation of the chalk rock. The groundwater flow and contaminant transport in the variably saturated chalk is solved using the control volume finite element method. Modelling the field tracer experiments performed in the variably saturated chalk shows the adequacy and usefulness of the new conceptual, mathematical and numerical model.     

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.     

Characterization of flow and transport processes within the unsaturated zone of Yucca Mountain, Nevada, under current and future climates

This paper presents a large-scale modeling study characterizing fluid flow and tracer transport in the unsaturated zone of Yucca Mountain, Nevada, a potential repository site for storing high-level radioactive waste. The study has been conducted using a three-dimensional numerical model, which incorporates a wide variety of field data and takes into account the coupled processes of flow and transport in the highly heterogeneous, unsaturated fractured porous rock. The modeling approach is based on a dual-continuum formulation of coupled multiphase fluid and tracer transport through fractured porous rock. Various scenarios of current and future climate conditions and their effects on the unsaturated zone are evaluated to aid in the assessment of the proposed repository's system performance using different conceptual models. These models are calibrated against field-measured data. Model-predicted flow and transport processes under current and future climates are discussed.     

Flow and transport in unsaturated fractured rock: effects of multiscale heterogeneity of hydrogeologic properties

The heterogeneity of hydrogeologic properties at different scales may have different effects on flow and transport processes in a subsurface system. A model for the unsaturated zone of Yucca Mountain, Nevada, is developed to represent complex heterogeneity at two different scales: (1) layer scale corresponding to geologic layering and (2) local scale. The layer-scale hydrogeologic properties are obtained using inverse modeling, based on the available measurements collected from the Yucca Mountain site. Calibration results show a significant lateral and vertical variability in matrix and fracture properties. Hydrogeologic property distributions in a two-dimensional, vertical cross-section of the site are generated by combining the average layer-scale matrix and fracture properties with local-scale perturbations generated using a stochastic simulation method. The unsaturated water flow and conservative (nonsorbing) tracer transport through the cross-section are simulated for different sets of matrix and fracture property fields. Comparison of simulation results indicates that the local-scale heterogeneity of matrix and fracture properties has a considerable effect on unsaturated flow processes, leading to fast flow paths in fractures and the matrix. These paths shorten the travel time of a conservative tracer from the source (repository) horizon in the unsaturated zone to the water table for small fractions of total released tracer mass. As a result, the local-scale heterogeneity also has a noticeable effect on global tracer transport processes, characterized by an average breakthrough curve at the water table, especially at the early arrival time of tracer mass. However, the effect is not significant at the later time after 20% tracer mass reaches the water table. The simulation results also verify that matrix diffusion plays an important role in overall solute transport processes in the unsaturated zone at Yucca Mountain.     

Comparison of the solution behaviour of a pyrite-calcite mixture in batch and unsaturated sand column

A successful application of reaction transport algorithms to calculate the chemical evolution of natural systems requires accurate methods to compute the rates of mineral/fluid surface reactions. Regarding the transport of radio-nuclides in mining dumps the dissolution of minerals is of special importance. Using a kinetic rate law of the mineral dissolution verified for unsaturated conditions will allow a realistic modelling of the mineral weathering in the environment. Dissolution rates of minerals in an aqueous solution are determined by several characteristics. These are surface reaction rates, morphology of the mineral's surface and, in case it is the unsaturated zone, the degree of the water saturation. For this process, the quantity of the particle surfaces which are in contact with percolating water is most decisive. In order to study the differences of mineral dissolution under saturated and unsaturated conditions batch and column experiments were carried out with a pyrite-calcite mixture. The experimental results were verified by calculations. Comparing the dissolution in batch with those in the column experiment, which was performed with a water flow velocity of 0.64 cm/day and was analyzed in the region of a water saturation of 0.11, one can conclude that only a small portion of about 5% of the grain surface is chemically reactive in this unsaturated flow.     

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.     

Interpretation of injection-withdrawal tracer experiments conducted between two wells in a large single fracture

Tracer experiments conducted using a flow field established by injecting water into one borehole and withdrawing water from another are often used to establish connections and investigate dispersion in fractured rock. As a result of uncertainty in the uniqueness of existing models used for interpretation, this method has not been widely used to investigate more general transport processes including matrix diffusion or advective solute exchange between mobile and immobile zones of fluid. To explore the utility of the injection-withdrawal method as a general investigative tool and with the intent to resolve the transport processes in a discrete fracture, two tracer experiments were conducted using the injection-withdrawal configuration. The experiments were conducted in a fracture which has a large aperture (>500 microm) and horizontally pervades a dolostone formation. One experiment was conducted in the direction of the hydraulic gradient and the other in the direction opposite to the natural gradient. Two tracers having significantly different values of the free-water diffusion coefficient were used. To interpret the experiments, a hybrid numerical-analytical model was developed which accounts for the arcuate shape of the flow field, advection-dispersion in the fracture, diffusion into the matrix adjacent to the fracture, and the presence of natural flow in the fracture. The model was verified by comparison to a fully analytical solution and to a well-known finite-element model. Interpretation of the tracer experiments showed that when only one tracer, advection-dispersion, and matrix diffusion are considered, non-unique results were obtained. However, by using multiple tracers and by accounting for the presence of natural flow in the fracture, unique interpretations were obtained in which a single value of matrix porosity was estimated from the results of both experiments. The estimate of porosity agrees well with independent measurements of porosity obtained from core samples. This suggests that: (i) the injection-withdrawal method is a viable tool for the investigation of general transport processes provided all relevant experimental conditions are considered and multiple conservative tracers are used; and (ii) for the conditions of the experiments conducted in this study, the dominant mechanism for exchange of solute between the fracture and surrounding medium is matrix diffusion.     

In situ tracer tests to determine retention properties of a block scale fracture network in granitic rock at the Aspö Hard Rock Laboratory, Sweden

Experiments were conducted at the Asp? Hard Rock Laboratory in order to improve the understanding of radionuclide retention properties of fractured crystalline bedrock in the 10-100 m scale (TRUE Block Scale Project, jointly funded by ANDRA, ENRESA, Nirex, JNC, Posiva and SKB). A series of tracer experiments were performed using sorbing tracers in three different flow paths. The different flow paths had Euclidian lengths of 14, 17 and 33 m, respectively, and one to three water conducting structures. Four tests were performed using different cocktails made up of radioactive sorbing tracers (22,24Na+, 42K+, 47Ca2+, 85Sr2+, 83,86Rb+, 131,133Ba2+ and 134,137Cs+). For each tracer injection, the breakthrough of sorbing tracers was compared to the breakthrough of a conservative tracer, 82Br-, 131I-, HTO and 186ReO4-, respectively. In the two longer flow paths, no breakthrough of 83Rb+ and 137Cs+ was observed after 8 months of pumping. Selected tracer tests were subject to basic modelling in which a one-dimensional (1D) advection-dispersion model, including surface sorption, and an unlimited matrix diffusion were used for the interpretation of the results. The results of the modelling indicated that there is a slightly higher mass transfer into a highly porous material in the block-scale experiment compared with in situ experiments performed over shorter distances and significantly higher than what would have been expected from laboratory data obtained from studies of the interactions in nonaltered intact rock.     

A high-resolution non-invasive approach to quantify oxygen transport across the capillary fringe and within the underlying groundwater

Oxygen transport across the capillary fringe is relevant for many biogeochemical processes. We present a non-invasive technique, based on optode technology, to measure high-resolution concentration profiles of oxygen across the unsaturated/saturated interface. By conducting a series of quasi two-dimensional flow-through laboratory experiments, we show that vertical hydrodynamic dispersion in the water-saturated part of the capillary fringe is the process limiting the mass transfer of oxygen. A number of experimental conditions were tested in order to investigate the influence of grain size and horizontal flow velocity on transverse vertical dispersion in the capillary fringe. In the same setup, analogous experiments were simultaneously carried out in the fully water-saturated zone, therefore allowing a direct comparison with oxygen transfer across the capillary fringe. The outcomes of the experiments under various conditions show that oxygen transport in the two zones of interest (i.e., the unsaturated/saturated interface and the saturated zone) is characterized by very similar transverse dispersion coefficients. An influence of the capillary fringe morphology on oxygen transport has not been observed. These results may be explained by the narrow grain size distribution used in the experiments, leading to a steep decline in water saturation at the unsaturated/saturated interface and to the absence of trapped gas in this transition zone. We also modeled flow (applying the van Genuchten and the Brooks-Corey relationships) and two-dimensional transport across the capillary fringe, obtaining simulated profiles of equivalent aqueous oxygen concentration that were in good agreement with the observations.     

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.     

Modeling hydrocarbon biodegradation in tidal aquifers with water-saturation and heat inhibition effects.

A model is developed for hydrocarbon biodegradation, which includes saturated and unsaturated flow, multi-species transport, heat transport, and bacterial growth processes. Numerical accuracy of the model was tested against analytical solutions. The model was also verified against laboratory results for a saturated-flow problem and reasonable match was obtained. Expressions are proposed for inhibition due to water content and temperature fluctuations. Bioactivities under cyclic water content variation were studied under no-flow conditions. A quantitative approach was used to reconcile some of the apparent contradictory conclusions regarding the efficiency of biodegradation of soils under wetting and drying conditions. The efficiency depends on the nature of the oxygenation process. For cases involving the presence of dissolved oxygen and the absence of O2 vapor, subjecting the soil to constant water content close to its optimal value for degradation is most efficient. However, wetting and drying can enhance degradation if O2 is only provided through aeration or direct contact between air and the medium. Also presented are the results of a typical field application of the model and a discussion of the effects of tides, saturation inhibition, and heat inhibition. Other inhibition factors, such as pH or salinity, can be easily incorporated in the formulation. The quantitative approach developed here can be used in assessing bioremediation not only in tidal aquifers but also in areas where water-table or temperature effects are of significance. The approach can be useful in the design of remediation strategies under water-flow or no-flow conditions involving water content and temperature fluctuations.     

Migration of contaminants through the unsaturated zone overlying the Hesbaye chalky aquifer in Belgium: a field investigation

This paper presents the results of a detailed field investigation that was performed for studying groundwater recharge processes and solute downward migration mechanisms prevailing in the unsaturated zone overlying a chalk aquifer in Belgium. Various laboratory measurements were performed on core samples collected during the drilling of boreholes in the experimental site. In the field, experiments consisted of well logging, infiltration tests in the unsaturated zone, pumping tests in the saturated zone and tracer tests in both the saturated and unsaturated zones. Results show that gravitational flows govern groundwater recharge and solute migration mechanisms in the unsaturated zone. In the variably saturated chalk, the migration and retardation of solutes is strongly influenced by recharge conditions. Under intense injection conditions, solutes migrate at high speed along the partially saturated fissures, downward to the saturated zone. At the same time, they are temporarily retarded in the almost immobile water located in the chalk matrix. Under normal recharge conditions, fissures are inactive and solutes migrate slowly through the chalk matrix. Results also show that concentration dynamics in the saturated zone are related to fluctuations of groundwater levels in the aquifer. A conceptual model is proposed to explain the hydrodispersive behaviour of the variably saturated chalk. Finally, the vulnerability of the chalk to contamination issues occurring at the land surface is discussed.     

Evaluation of saline tracer performance during electrical conductivity groundwater monitoring

Saline solutions are the most commonly used hydrological tracers, because they can be easily and economically monitored by in situ instrumentation such as electrical conductivity (EC) loggers in wells or by geoelectrical measurements. Unfortunately, these low-cost techniques only provide information on the total concentration of ions in solution, i.e., they cannot resolve the ionic composition of the aqueous solution. This limitation can introduce a bias in the estimation of aquifer parameters where sorption phenomena between saline tracers and sediments become relevant. In general, only selected anions such as Cl(-) and Br(-) are recognised to be transported unretarded and they are referred to as conservative tracers or mobile anions. However, cations within the saline tracer may interact with the soil matrix through a range of processes such as ion exchange, surface complexation and via physical mass-transfer phenomena. Heterogeneous reactions with minerals or mineral surfaces may not be negligible where aquifers are composed of fine alluvial sediments. The focus of the present study was to examine and to quantify the bias between the aquifer parameters estimated during model-based interpretation of experimental data of EC measurements of saline tracer relative to the aquifer parameters found by specific measurements (i.e. via ionic chromatography, IC) of truly conservative species. To accomplish this, column displacement experiments with alluvial aquifer materials collected from the Po lowlands (Italy) were performed under water saturated conditions. The behaviour of six selected, commonly used saline tracers (i.e., LiCl, KCl, and NaCl; LiBr, KBr, and NaBr) was studied and the data analysed by inverse modelling. The results demonstrate that the use of EC as a tracer can lead to an erroneous parameterisation of the investigated porous media, if the reactions between solute and matrix are neglected. In general, errors were significant except for KCl and KBr, which is due to the weak interaction between dissolved K(+) and the sediment material. The study shows that laboratory scale pre-investigations can help with tracer selection and to optimise the concentration range targeted for in situ multilevel monitoring by unspecific geoelectrical instrumentation.     

Diffusion of tritiated water and 22Na+ through non-degraded hardened cement pastes

Diffusion experiments through hardened cement pastes (HCP) using tritiated water (HTO) and 22Na(+), considered to be conservative tracers, have been carried out in triplicates in a glove box under a controlled nitrogen atmosphere. Each experiment consisted of a through-diffusion test followed by an out-diffusion test. The experimental data were inversely modelled applying an automated Marquardt-Levenberg procedure. The analysis of the through-diffusion data allowed the extraction of values for the effective diffusion coefficients, D(e), and the rock capacity factor, alpha. Good agreement between measured and calculated tracer breakthrough curves was achieved using both a simple diffusion model without sorption and a diffusion/linear sorption model. The best-fit K(d)-values were found to be consistent with R(d)-values measured in previous batch-sorption experiments. The best-fit values from the through-diffusion tests were then used to predict the results of subsequent out-diffusion experiments. Good agreement between experimental data and predictions was achieved only for the case of linear sorption. Isotopic exchange can only partially account for both the amount of tracer taken up in the batch-sorption tests and the measured retardation in the diffusion experiments and, hence, additional mechanisms have to be invoked to explain the data.     

Strontium migration in a crystalline medium: effects of the presence of bentonite colloids

The effects of bentonite colloids on strontium migration in fractured crystalline medium were investigated. We analyzed first the transport behaviour of bentonite colloids alone at different flow rates; then we compared the transport behaviour of strontium as solute and of strontium previously adsorbed onto stable bentonite colloids at a water velocity of approximately 7.1·10(-6)m/s-224m/yr. Experiments with bentonite colloids alone showed that - at the lowest water flow rate used in our experiments (7.1·10(-6)m/s) - approximately 70% of the initially injected colloids were retained in the fracture. Nevertheless, the mobile colloidal fraction, moved through the fracture without retardation, at any flow rate. Bentonite colloids deposited over the fracture surface were identified during post-mortem analyses. The breakthrough curve of strontium as a solute, presented a retardation factor, R(f)~6, in agreement with its sorption onto the granite fracture surface. The breakthrough curve of strontium in the presence of bentonite colloids was much more complex, suggesting additional contributions of colloids to strontium transport. A very small fraction of strontium adsorbed on mobile colloids moved un-retarded (R(f)=1) and this fraction was much lower than the expected, considering the quantity of strontium initially adsorbed onto colloids (90%). This behaviour suggests the hypothesis of strontium sorption reversibility from colloids. On the other hand, bentonite colloids retained within the granite fracture played a major role, contributing to a slower strontium transport in comparison with strontium as a solute. This was shown by a clear peak in the breakthrough curve corresponding to a retardation factor of approximately 20.