Abstract: Clay mineral-rich sedimentary formations are currently under investigation to evaluate their potential use as host formations for installation of deep underground disposal facilities for radioactive waste (e.g. Boom Clay (BE), Opalinus Clay (CH), CallovoâOxfordian argillite (FR)). The ultimate safety of the cor- responding repository concepts depends largely on the capacity of the host formation to limit the flux towards the biosphere of radionuclides (RN) contained in the waste to acceptably low levels. Data for dif- fusion-driven transfer in these formations shows extreme differences in the measured or modelled behaviour for various radionuclides, e.g. between halogen RN (36Cl, 129I) and actinides (238,235U, 237Np, 232Th, etc.), which result from major differences between RN of the effects on transport of two phenom- ena: diffusion and sorption. This paper describes recent research aimed at improving understanding of these two phenomena, focusing on the results of studies carried out during the EC Funmig IP on clayrocks from the above three formations and from the Boda formation (HU).
Project results regarding phenomena governing water, cation and anion distribution and mobility in the pore volumes influenced by the negatively-charged surfaces of clay minerals show a convergence of the modelling results for behaviour at the molecular scale and descriptions based on electrical double layer models. Transport models exist which couple ion distribution relative to the clayâsolution interface and differentiated diffusive characteristics. These codes are able to reproduce the main trends in behav- iour observed experimentally, e.g. De(anion) < De(HTO) < De(cation) and De(anion) variations as a function of ionic strength and material density. These trends are also well-explained by models of transport through ideal porous matrices made up of a charged surface material. Experimental validation of these models is good as regards monovalent alkaline cations, in progress for divalent electrostatically-interact- ing cations (e.g. Sr2+) and still relatively poor for âstrongly sorbingâ, high Kd cations.
Funmig results have clarified understanding of how clayrock mineral composition, and the correspond- ing organisation of mineral grain assemblages and their associated porosity, can affect mobile solute (anions, HTO) diffusion at different scales (mm to geological formation). In particular, advances made in the capacity to map clayrock mineral grain-porosity organisation at high resolution provide additional elements for understanding diffusion anisotropy and for relating diffusion characteristics measured at different scales. On the other hand, the results of studies focusing on evaluating the potential effects of heterogeneity on mobile species diffusion at the formation scale tend to show that there is a minimal effect when compared to a homogeneous property model. Finally, the results of a natural tracer-based study carried out on the Opalinus Clay formation increase confidence in the use of diffusion parameters measured on laboratory scale samples for predicting diffusion over geological timeâspace scales.
Much effort was placed on improving understanding of coupled sorptionâdiffusion phenomena for sorbing cations in clayrocks. Results regarding sorption equilibrium in dispersed and compacted materi- als for weakly to moderately sorbing cations (Sr2+, Cs+, Co2+) tend to show that the same sorption model probably holds in both systems. It was not possible to demonstrate this for highly sorbing elements such as Eu(III) because of the extremely long times needed to reach equilibrium conditions, but there does not seem to be any clear reason why such elements should not have similar behaviour. Diffusion experiments carried out with Sr2+, Cs+ and Eu(III) on all of the clayrocks gave mixed results and tend to show that cou- pled diffusionâsorption migration is much more complex than expected, leading generally to greater mobility than that predicted by coupling a batch-determined Kd and Fick Ìs law based on the diffusion behaviour of HTO. If the Kd measured on equivalent dispersed systems holds as was shown to be the case for Sr, Cs (and probably Co) for Opalinus Clay, these results indicate that these cations have a De value higher than HTO (up to a factor of 10 for Cs+). Results are as yet very limited for very moderate to strongly sorbing species (e.g. Co(II), Eu(III), Cu(II)) because of their very slow transfer characteristics.
Abstract: Diffusion of cations and other contaminants through clays is of central interest, because clays and clay rocks are widely considered as barrier materials for waste disposal sites. An intriguing experimental observation has been made in this context: Often, the diffusive flux of cations at trace concentrations is much larger and the retardation smaller than expected based on their sorption coefficients. So-called surface diffusion of sorbed cations has been invoked to explain the observations but remains a controversial issue. Moreover, the corresponding surface diffusion coefficients are largely unknown. Here we show that, by an appiopriate scaling, published diffusion data covering a broad range of cations, clays, and chemical conditions can all be modeled satisfactorily by a surface diffusion model. The average mobility of sorbed cations seems to be primarily an intrinsic property of each cation that follows inversely its sorption affinity. With these surface mobilities, cation diffusion coefficients can now be estimated from those of water tracers. In pure clays at low salinities, surface diffusion can reduce the cation retardation by a factor of more than 1000.
Abstract: Materials with low permeability such as clays are major components in engineered or natural barriers in waste disposal sites. Diffusion of dissolved chemicals through such materials is a relatively slow process that can strongly limit the spreading of contaminants. The macroscopic diffusive transport of chemicals through barrier materials originates from Brownian motion of molecules and ions in the solution and their interaction with the surface of the solid medium. Up-scaling these molecular phenomena to the continuum scale, which is required for large-scale and long-time predictions, is particularly challenging because the considered scales differ by orders of magnitude. To address the up-scaling problem we developed a two-step simulation approach which enables us to derive macroscopic diffusion coefficients of water and ions for continuum equations from pore scale molecular diffusion coefficients. Our starting points for the up-scaling procedure were local pore diffusion coefficients derived in molecular dynamics simulations for specific local environments, such as the interlayer or edge regions of clay particles. We then assigned these local diffusion coefficients to different types of porosity of a model clay structure and obtained the structure-averaged diffusion coefficient of the sample by random walk simulations. We demonstrated the validity of the suggested approach by comparing the results of direct molecular simulations for a stack of pyrophyllite nanoparticles with the two-step up-scaling approach. We then investigated the effects of mineralogical heterogeneities and of anion exclusion on the larger-scale diffusion coefficients for a model clay structure which resembles a typical clay material. Our up-scaling concept is general and can be used for up-scaling molecular diffusion coefficients for porous materials with almost arbitrarily complex structures.
Abstract: Reactive transport modelling was used to simulate solute transport, thermodynamic reactions, ion exchange and biodegradation in the Porewater Chemistry (PC) experiment at the Mont Terri Rock Laboratory. Simulations show that the most important chemical processes controlling the fluid composition within the borehole and the surrounding formation during the experiment are ion exchange, biodegradation and dissolution/precipitation reactions involving pyrite and carbonate minerals. In contrast, thermodynamic mineral dissolution/precipitation reactions involving alumo-silicate minerals have little impact on the fluid composition on the time-scale of the experiment. With the accurate description of the initial chemical condition in the formation in combination with kinetic formulations describing the different stages of bacterial activities, it has been possible to reproduce the evolution of important system parameters, such as the pH, redox potential, total organic C. dissolved inorganic C and SO(4) concentration. Leaching of glycerol from the pH-electrode may be the primary source of organic material that initiated bacterial growth, which caused the chemical perturbation in the borehole. Results from these simulations are consistent with data from the over-coring and demonstrate that the Opalinus Clay has a high buffering capacity in terms of chemical perturbations caused by bacterial activity. This buffering capacity can be attributed to the carbonate system as well as to the reactivity of clay surfaces. (C) 2011 Elsevier Ltd. All rights reserved.
Abstract: Argillaceous formations generally act as aquitards because of their low hydraulic conductivities. This property, together with the large retention capacity of clays for cationic contaminants, has brought argillaceous formations into focus as potential host rocks for the geological disposal of radioactive and other waste. In several countries, programmes are under way to characterise the detailed transport properties of such formations at depth. In this context, the interpretation of profiles of natural tracers in pore waters across the formations can give valuable information about the large-scale and long-term transport behaviour of these formations. Here, tracer-profile data, obtained by various methods of pore-water extraction for nine sites in central Europe, are compiled. Data at each site comprise some or all of the conservative tracers: anions (Cl(-), Br(-)), water isotopes (delta(18)O, delta(2)H) and noble gases (mainly He). Based on a careful evaluation of the palaeo-hydrogeological evolution at each site, model scenarios are derived for initial and boundary pore-water compositions and an attempt is made to numerically reproduce the observed tracer distributions in a consistent way for all tracers and sites, using transport parameters derived from laboratory or in situ tests. The comprehensive results from this project have been reported in Mazurek et al. (2009). Here the results for three sites are presented in detail, but the conclusions are based on model interpretations of the entire data set. In essentially all cases, the shapes of the profiles can be explained by diffusion acting as the dominant transport process over periods of several thousands to several millions of years and at the length scales of the profiles. Transport by advection has a negligible influence on the observed profiles at most sites, as can be shown by estimating the maximum advection velocities that still give acceptable fits of the model with the data. The advantages and disadvantages of different conservative tracers are also assessed. The anion Cl(-) is well suited as a natural tracer in aquitards, because its concentration varies considerably in environmental waters. It can easily be measured, although the uncertainty regarding the fraction of the pore space that is accessible to anions in clays remains an issue. The stable water isotopes are also well suited, but they are more difficult to measure and their values generally exhibit a smaller relative range of variation. Chlorine isotopes (delta(37)Cl) and He are more difficult to interpret because initial and boundary conditions cannot easily be constrained by independent evidence. It is also shown that the existence of perturbing events such as the activation of aquifers due to uplift and erosion, leading to relatively sharp changes of boundary conditions, can be considered as a pre-requisite to obtain well-interpretable tracer signatures. On the other hand, gradual changes of boundary conditions are more difficult to parameterise and so may preclude a clear interpretation. (C) 2011 Elsevier Ltd. All rights reserved.
Abstract: The experimental verification of matrix diffusion in crystalline rocks largely relies on indirect methods performed in the laboratory. Such methods are prone to perturbations of the rock samples by collection and preparation and therefore the laboratory-derived transport properties and fluid composition might not represent in situ conditions. We investigated the effects induced by the drilling process and natural rock stress release by mass balance considerations and sensitivity analysis of analytical out-diffusion data obtained from originally saturated, large-sized drillcore material from two locations drilled using traced drilling fluid. For in situ stress-released drillcores of quartz-monzodiorite composition from the Aspo HRL, Sweden, tracer mass balance considerations and 1D and 2D diffusion modelling consistently indicated a contamination of <1% of the original pore water. This chemically disturbed zone extends to a maximum of 0.1 mm into the drillcore (61.8 mm x 180.1 mm) corresponding to about 0.66% of the total pore volume (0.77 vol.%). In contrast, the combined effects of stress release and the drilling process, which have influenced granodioritic drillcore material from 560 m below surface at Forsmark. Sweden, resulted in a maximum contamination of the derived porewater Cl(-) concentration of about 8%. The mechanically disturbed zone with modified diffusion properties covers the outermost similar to 6 mm of the drillcore (50 mm x 189 mm), whereas the chemically disturbed zone extends to a maximum of 0.3 mm based on mass balance considerations, and to 0.15 mm to 0.2 mm into the drillcore based on fitting the observed tracer data. This corresponds to a maximum of 2.4% of the total pore volume (0.62 vol.%) being affected by the drilling-fluid contamination. The proportion of rock volume affected initially by drilling fluid or subsequently with experiment water during the laboratory diffusion and re-saturation experiments depends on the size of the drillcore material and will become larger the smaller the sample used for the experiment. The results are further in support of matrix diffusion taking place in the undisturbed matrix of crystalline rocks at least in the cm range. (C) 2011 Elsevier B.V. All rights reserved.
Abstract: The Opalinus Clay formation in North Switzerland is a potential host rock for a deep underground radioactive waste repository. The distribution of U-238, U-234 and Th-230 was studied in rock samples of the Opalinus Clay from an exploratory borehole at Benken (Canton of Zurich) using MC-ICP-MS. The aim of U-234 was to assess the in situ, long-term migration behaviour in this rock. Very low hydraulic conductivities of the Opalinus Clay, reducing potential of the pore water and its chemical equilibrium with the host rock are expected to render both U-238 and Th-230 immobile. If U is heterogeneously distributed in the Opalinus Clay, gradients in the supply of U-234 from the rock matrix to the pore water by the decay of U-238 will be established. Diffusive redistribution separates U-234 from its immobile parent U-238 resulting in bulk rock U-234/U-238 activity disequilibria. These may provide a means of estimating the mobility of U-234 in the rock if the diffusion rate of U-234 is significant compared to its decay rate. Sampling was carried out on two scales. Drilling of cm-spaced samples from the drill-core was done to study mobility over short distances and elucidate possible small-scale lithological control. Homogenized 25-cm-long portions of a 2-m-long drill-core section were prepared to provide information on transport over a longer distance. Variations in U and/or Th content on the cm-scale between clays and carbonate-sandy layers are revealed by beta-scanning, which shows that the (dominant) clay is richer in both elements. Samples were digested using aqua regia followed by total HF dissolution, yielding two fractions. in all studied samples U was found to be concentrated in the HF digestion fraction. It has a high U/Th ratio and a study by SEM-EDS points to sub-mu m up to several mu m in size zircon grains as the main U-rich phase. This fraction consistently has U-234/U-238 activity ratios below unity. The minute zircon grains constitute the major reservoir of U in the rock and act as constant rate suppliers of U-234 into the rock matrix and the pore water. The aqua regia leach fraction was found to be enriched in Th, and complementary to the HF fraction, having U-234/U-238 activity ratios above unity. It is believed that these U activity ratios reflect the surplus of having U-234 delivered from the zircon grains. Some cm-spaced samples show bulk rock U-234/U-238 activity ratios that are markedly out of equilibrium. In most of them a striking negative correlation between the total U content and the bulk rock U-234/U-238 activity ratios is observed. This is interpreted to indicate net U-234 transfer from regions of higher supply of U-234 towards those of lower supply which is, in most cases, equivalent to transfer from clayey towards carbonate/sandy portions of the rock. In contrast, the 25 cm averaged samples all have uniform bulk rock U-234/U-238 activity ratios in equilibrium, indicating U immobility in the last 1-1.5 Ma on this spatial scale. It is concluded that the small-scale lithological variations which govern U spatial distribution in the Opalinus Clay are the major factor determining U-234 in situ supply rates, regulating its diffusive fluxes and controlling the observed bulk rock U-234/U-238 activity ratios. A simple box-model is presented to simulate the measured bulk rock U-234/U-238 activity ratios and to give an additional insight into the studied system. (C) 2008 Elsevier Ltd. All rights reserved.
Abstract: Diffusion of water and solutes through compacted clays or claystones is important when assessing the barrier function of engineered or geological barriers in waste disposal. The shape and the connectivity of the pore network as well as electrostatic interactions between the diffusant and the charged clay surfaces or cations compensating negative surface charges affect the resistance of the porous medium to diffusion. Comparing diffusion measurements performed at different spatial or time scales allows identification and extraction of the different factors. We quantified the electrostatic constraint q for five different, highly compacted clays (rho(b) = 1.85 +/- 0.05 g/cm(3)) using quasielastic neutron scattering (DENS) data. We then compared the DENS data with macroscopic diffusion data for the same clays and could derive the true geometric tortuosities G of the samples. Knowing the geometric and electrostatic factors for the different clays is essential when trying to predict diffusion coefficients for other conditions. We furthermore compared the activation energies E-a for diffusion at the two measurement scales. Because E-a is mostly influenced by the local, pore scale surroundings of the water, we expected the results to be similar at bath scales. This was indeed the case for the nonswelling clays kaolinite and illite, which had E-a values lower than that of bulk water, but not for montmorillonite, which had values lower than that in bulk water at the microscopic scale, but larger at the macroscopic scale. The differences could be connected to the strongly temperature dependent mobility of the cations in the clays, which may act as local barriers in the narrow pores at low temperatures.
Abstract: The effect of temperature and ionic strength on the diffusion of HTO parallel to the direction of compaction through 5 highly compacted clay minerals (bulk dry density, P-b.d = 1.90 +/- 0.05 Mg/m(3)), namely montmorillonite (Na- and Ca-form), illite (Na- and Ca-form), and kaolinite, was studied. The diffusion experiments were carried out at temperatures between 0 degrees C and 60 degrees C and at ionic strengths of 0.01 M and 1 M NaCl for the Na-form clays and kaolinite, and of 0.005 M and 0.5 M CaCl2 for the Ca-form. The ionic strength had an insignificant influence on the values of the effective diffusion coefficient (variation by less than 10%) for the clays under study at this degree of compaction. The effective diffusion coefficients followed the order Na-montmorillonite<Ca-montmorillonite<Caillite < Na-illite <= kaolinite. It is thought that the differences between Na- and Ca-montmorillonite originate from the larger size particles, and thus the lower tortuosity of the latter; whereas the differences between Na- and Ca-illite are related to the different degree of solvation of the Na and Ca cations. The activation energies were successfully calculated using the Arrhenius law. Swelling clays (Na- and Ca-montmorillonite) had slightly larger activation energy values (20 kJ/mol) compared to bulk water (17 kJ/mol): Ca-illite (16 kJ/mol), Na-illite (13 kJ/mol) and kaolinite (14.4 kJ/mol) lower values than that of bulk water. The low activation energies of the last three clays may be related to weaker H-bonds between water and the clay surfaces compared to those in bulk water. (c) 2008 Elsevier Ltd. All rights reserved.
Abstract: Argillaceous rocks are considered in Switzerland and in other countries as potential host rocks for the deep geologic disposal of radioactive waste. Opalinus Clay, a formation of Jurassic age, is at present the most favored candidate in Switzerland. It has a low hydraulic conductivity, no or only few active fractures, and a large retention capacity for sorbing solutes, which make this sediment well suited as an additional barrier for the spreading of contaminants. Transport through Opalinus Clay is typically dominated by diffusion. Accordingly, the diffusion of tracers through this material is intensively studied. In the underground research facility at Mont Terri in Switzerland, field experiments are performed that aim at investigating the diffusion and retention of sorbing tracers at intermediate scales and under relevant in-situ conditions. The field experiments use a borehole, from which tracers diffuse into the surrounding rock. The tracer cocktail in the borehole is continuously circulated, which allows to monitor the tracer decrease over time. When modeling these experiments, care has to be taken to correctly represent the inlet system. In this paper, it is shown how the numerical representation of the borehole and inlet system and the spatial discretization affect the simulation results of mobile and sorbing tracers. For mobile tracers, it is generally sufficient to use an effective diffusion coefficient for the circulated fluid about 30 times larger than that in the rock to mimic the continuous mixing. In contrast, for sorbing tracers even a 7000 times larger diffusion coefficient may, at early times, not lead to homogeneous borehole concentrations. This is because the equilibrium sorption on the rock quickly and drastically reduces the tracer concentrations at the interface. As a consequence, the simulated flux into the rock becomes too small, and the calculated average decrease of the borehole concentration is much too slow even for larger times. The slow decrease can be similar to that simulated for a much less sorbing solute, which of course would critically affect parameter estimation when fitting the model to observed data. Increasing the borehole diffusion coefficient to very large values to obtain a complete mixing can lead to a too fast simulated concentration decrease, if the spatial discretization is insufficient. Thus, a careful checking of the numerical results is required for the strongly sorbing tracers.
Abstract: The migration of radioactive and chemical contaminants in clay materials and argillaceous host rocks is characterised by diffusion and retention processes. Valuable information on such processes can be gained by combining diffusion studies at laboratory scale with field migration tests. In this work, the outcome of a multi-tracer in situ migration test performed in the Opalinus Clay formation in the Mont Terri underground rock laboratory (Switzerland) is presented. Thus, 1.16 x 10(5) Bq/L of HTO, 3.96 x 10(3) Bq/L of Sr-85, 6.29 x 10(2) Bq/L of Co-60, 2.01 x 10(-3) mol/L Cs, 9.10 x 10(-4) mol/L I and 1.04 x 10(-3) mol/L Br were injected into the borehole. The decrease of the radioisotope concentrations in the borehole was monitored using in situ gamma-spectrometry. The other tracers were analyzed with state-of-the-art laboratory procedures after sampling of small water aliquots from the reservoir. The diffusion experiment was carried out over a period of one year after which the interval section was overcored and analyzed. Based on the experimental data from the tracer evolution in the borehole and the tracer profiles in the rock, the diffusion of tracers was modelled with the numerical code CRUNCH. The results obtained for HTO (H-3), I- and Br- confirm previous lab and in situ diffusion data. Anionic fluxes into the formation were smaller compared to HTO because of anion exclusion effects. The migration of the cations Sr-85(2+), Cs+ and Co-60(2+) was found to be governed by both diffusion and sorption processes. For Sr-85(2+), the slightly higher diffusivity relative to HTO and the low sorption value are consistent with laboratory diffusion measurements on small-scale samples. In the case of Cs+, the numerically deduced high diffusivity and the Freundlich-type sorption behaviour is also supported by ongoing laboratory data. For Co, no laboratory diffusion data were yet available for comparison; however, the modelled data suggests that Co-60(2+) sorption was weaker than would be expected from available batch sorption data. Overall, the results demonstrate the feasibility of the experimental setup for obtaining high-quality diffusion data for conservative and sorbing tracers. (C) 2007 Elsevier Ltd. All rights reserved.
Abstract: The freezing behavior of water confined in compacted charged and uncharged clays (montmorillonite in Na-and Ca-forms, illite in Na-and Ca-forms, kaolinite and pyrophyllite) was investigated by neutron scattering. Firstly, the amount of frozen (immobile) water was measured as a function of temperature at the IN16 backscattering spectrometer, Institute Laue-Langevin (ILL). Water in uncharged, partly hydrophobic (kaolinite) and fully hydrophobic (pyrophyllite) clays exhibited a similar freezing and melting behavior to that of bulk water. In contrast, water in charged clays which are hydrophilic could be significantly supercooled. To observe the water dynamics in these clays, further experiments were performed using quasielastic neutron scattering. At temperatures of 250, 260 and 270 K the diffusive motion of water could still be observed, but with a strong reduction in the water mobility as compared with the values obtained above 273 K. The diffusion coefficients followed a non-Arrhenius temperature dependence well described by the Vogel-Fulcher-Tammann and the fractional power relations. The fits revealed that Na-and Ca-montmorillonite and Ca-illite have similar Vogel-Fulcher-Tammann temperatures (T-VFT, often referred to as the glass transition temperature) of similar to 120 K and similar temperatures at which the water undergoes the âstrong-fragileâ transition, T-s similar to 210 K. On the other hand, Na-illite had significantly larger values of T-VFT similar to 180 K and T-s similar to 240 K. Surprisingly, Ca-illite has a similar freezing behavior of water to that of montmorillonites, even though it has a rather different structure. We attribute this to the stronger hydration of Ca ions as compared with the Na ions occurring in the illite clays.
Abstract: The water diffusion in four different, highly compacted clays [montmorillonite in the Na- and Ca-forms, illite in the Na- and Ca-forms, kaolinite, and pyrophyllite (bulk dry density rho(b)=1.85 +/- 0.05 g/cm(3))] was studied at the atomic level by means of quasielastic neutron scattering. The experiments were performed on two time-of-flight spectrometers and at three different energy resolutions [FOCUS at SINQ, PSI (3.65 and 5.75 A), and TOFTOF at FRM II (10 A)] for reliable data analysis and at temperatures between 27 and 95 degrees C. Two different jump diffusion models were used to describe the translational motion. Both models describe the data equally well and give the following ranking of diffusion coefficients: Na-montmorillonite <= Ca-montmorillonite < Ca-illite < Na-illite < water <= pyrophyllite <= kaolinite. Uncharged clays had slightly larger diffusion coefficients than that of bulk water due to their hydrophobic surfaces. The time between jumps, tau(t), follows the sequence: Ca-montmorillonite >= Na-montmorillonite>Ca-illite>Na-illite >= kaolinite>pyrophyllite >= water, in both jump diffusion models. For clays with a permanent layer charge (montmorillonite and illite) a reduction in the water content by a factor of 2 resulted in a decrease in the self-diffusion coefficients and an increase in the time between jumps as compared to the full saturation. The uncharged clay kaolinite exhibited no change in the water mobility between the two hydration states. The rotational relaxation time of water was affected by the charged clay surfaces, especially in the case of montmorillonite; the uncharged clays presented a waterlike behavior. The activation energies for translational diffusion were calculated from the Arrhenius law, which adequately describes the systems in the studied temperature range. Na- and Ca-montmorillonite (similar to 11-12 kJ/mol), Na-illite (similar to 13 kJ/mol), kaolinite and pyrophyllite (similar to 14 kJ/mol), and Ca-illite (similar to 15 kJ/mol) all had lower activation energies than bulk water (similar to 17 kJ/mol in this study). This may originate from the reduced number and strength of the H-bonds between water and the clay surfaces, or ions, as compared to those in bulk water. Our comparative study suggests that the compensating cations in swelling clays have only a minor effect on the water diffusion rates at these high densities, whereas these cations influence the water motion in non-swelling clays.
Abstract: [ 1] In order to characterize the large-scale transport properties of the Opalinus Clay formation, the pore water isotope composition (delta O-18 and delta H-2) was determined on samples from the deep borehole Benken ( northeastern Switzerland) across Jurassic argillaceous rocks. The sequence of claystones and marls, delimited by two aquifers, is located at depth from about 400 to 700 m and exhibits very low hydraulic conductivities ( below 10(-13) m s(-1)). The isotope data of the pore water were obtained from core samples by diffusive vapor equilibration, vacuum distillation, and squeezing. Compared with the other methods, vacuum distillation led to too low values. To evaluate the large-scale transport properties of the formation, we performed a series of advective-dispersive model calculations and compared them with the experimental data. In accordance with the hydrogeological history, we varied initial and boundary conditions as well as model parameters. The main results can be summarized as follows: ( 1) Molecular diffusion to the underlying aquifer can explain the general features of the isotope profiles, ( 2) no signatures of advective flow could be detected, ( 3) the evolution time is of the order of 0.51 Ma ( relying on laboratory diffusion coefficients) with a possible range of about 0.22 Ma, which is geologically plausible, and ( 4) parameters measured on small scales ( centimeters or meters and months) are also plausible at the formation scale ( tens of meters and millions of years) for the sediments investigated.
Abstract: The translational diffusion of water in compacted clays at a high hydration level has been investigated by quasielastic neutron scattering at a time-of-flight spectrometer FOCUS (SINQ). Four compacted clays with systematic structural differences have been studied: Na-montmorillonite, Na-illite, kaolinite and pyrophyllite. The QENS experiments were performed using two different incident wavelengths in order to access a larger Q range and verify the data analysis. The translational diffusion coefficient for water in Na-montmorillonite and Na-illite are lower than those for bulk water, whereas the preliminary results for kaolinite and pyrophyllite show larger diffusion coefficient.
Abstract: The Opalinus Clay in Northern Switzerland has been identified as a potential host rock formation for the disposal of radioactive waste. Comprehensive understanding of gas transport processes through this low-permeability formation forms a key issue in the assessment of repository performance. Field investigations and laboratory experiments suggest an intrinsic permeability of the Opalinus Clay in the order of 10(-20) to 10(-21) m(2) and a moderate anisotropy ratio < 10. Porosity depends on clay content and burial depth; values of similar to 0.12 are reported for the region of interest. Porosimetry indicates that about 10-30% of voids can be classed as macropores, corresponding to an equivalent pore radius > 25 nm. The determined entry pressures are in the range of 0.4-10 MPa and exhibit a marked dependence on intrinsic permeability. Both in situ gas tests and gas permeameter tests on drillcores demonstrate that gas transport through the rock is accompanied by porewater displacement, suggesting that classical flow concepts of immiscible displacement in porous media can be applied when the gas entry pressure (i.e. capillary threshold pressure) is less than the minimum principal stress acting within the rock. Essentially, the pore space accessible to gas flow is restricted to the network of connected macropores, which implies a very low degree of desaturation of the rock during the gas imbibition process. At elevated gas pressures (i.e. when gas pressure approaches the level of total stress that acts on the rock body), evidence was seen for dilatancy controlled gas transport mechanisms. Further field experiments were aimed at creating extended tensile fractures with high fracture transmissivity (hydro- or gasfracs). The test results lead to the conclusion that gas fracturing can be largely ruled out as a risk for post-closure repository performance.
Abstract: The diffusion properties of the Opalinus Clay were studied in the underground research laboratory at Mont Terri (Canton Jura, Switzerland) and the results were compared with diffusion data measured in the laboratory on small-scale samples. The diffusion of HTO, Na-22(+), Cs+ and I- were investigated for a period of 10 months. The diffusion equipment used in the field experiment was designed in such a way that a solution of tracers was circulated through a sintered metal screen placed at the end of a borehole drilled in the formation. The concentration decrease caused by the diffusion of tracers into the rock could be followed with time and allowed first estimations of the effective diffusion coefficient. After 10 months, the diffusion zone was over-cored and the tracer profiles measured. From these profiles, effective diffusion coefficients and rock capacity factors Could be extracted by applying a two-dimensional transport model including diffusion and sorption. The simulations were done with the reactive transport code CRUNCH. In addition, results obtained from through-diffusion experiments oil small-sized samples with HTO, Cl-36(-) and Na-22(+) are presented and compared with the in situ data. In all cases. excellent agreement between the two data sets exists. Results for Cs+ indicated five times higher diffusion rates relative to HTO. Corresponding laboratory diffusion measurements are still lacking. However. our Cs+ data are in qualitative agreement wish through-diffusion data for Callovo-Oxfordian argillite rock samples. which also indicate significantly higher effective diffusivities for Cs+ relative to HTO.
Abstract: [1] Natural soil profiles may be interpreted as an arrangement of parts which are characterized by properties like hydraulic conductivity and water retention function. These parts form a complicated structure. Characterizing the soil structure is fundamental in subsurface hydrology because it has a crucial influence on flow and transport and defines the patterns of many ecological processes. We applied an image analysis method for recognition and classification of visual soil attributes in order to model flow and transport through a man-made soil profile. Modeled and measured saturation-dependent effective parameters were compared. We found that characterizing and describing conductivity patterns in soils with sharp conductivity contrasts is feasible. Differently, solving flow and transport on the basis of these conductivity maps is difficult and, in general, requires special care for representation of small-scale processes.
Abstract: Water flow and solute transport through soils are strongly influenced by the spatial arrangement of soil materials with different hydraulic and chemical properties. Knowing the specific or statistical arrangement of these materials is considered as a key toward improved predictions of solute transport. Our aim was to obtain two-dimensional material maps from photographs of exposed profiles. We developed a segmentation and classification procedure and applied it to the images of a very heterogeneous sand tank, which was used for a series of flow and transport experiments. The segmentation was based on thresholds of soil color, estimated from local median gray values, and of soil texture, estimated from local coefficients of variation of gray values. Important steps were the correction of inhomogeneous illumination and reflection, and the incorporation of prior knowledge in filters used to extract the image features and to smooth the results morphologically. We could check and confirm the success of our mapping by comparing the estimated with the designed sand distribution in the tank. The resulting material map was used later as input to model flow and transport through the sand tank. Similar segmentation procedures may be applied to any high-density raster data, including photographs or spectral scans of field profiles.
Abstract: The diffusion of radionuclides is an important safety aspect for nuclear waste disposal in argillaceous host rocks. A long-term diffusion experiment, termed DI-A, is being carried out at the Mont Terri Rock Laboratory in the Opalinus Clay formation. The aim of this experiment is the understanding of the migration and sorption behaviour of cationic and anionic species in consolidated clays. This study reports on the experimental layout and the first results obtained from the DI-A experiment, which include the investigation of HTO, Na-22(+), Cs+, and I- migration during a period of 1 year by analysing these tracers in the water circulating in the borehole. In addition, results obtained from through-diffusion experiments on small-sized samples with HTO, I-, and Cl-36(-) are presented. The decrease of tracer concentrations in the borehole is fastest for Cs+, followed by Na-22(+), HTO, and finally I-. The chemical composition of the artificial pore water in the borehole shows very little variation with time, thus indicating almost no chemical disturbance around the borehole. Through-diffusion experiments in the laboratory that were performed parallel to the bedding plane with two different methods yielded effective diffusion coefficients for HTO of 4-5 X 10(-11) m(2) s(-1) and significantly lower ones for anions Cl- and I- (0.7-1.6 X 10(-11) m(2) s(-1)). The results indicate the importance of anion exclusion effects arising from the negatively charged clay surfaces. Furthermore, they demonstrate the anisotropic diffusion properties of the clay formation with significantly increased diffusion rates parallel to bedding relative to the perpendicular direction. The tracer data of the in situ experiment were successfully described with 2D diffusion models using diffusion and sorption parameters obtained from the above mentioned and other laboratory studies. The modelling results indicate that HTO and I- diffused with no retardation. The retardation of Na+ and Cs+ could be described by empirical sorption expressions from previously derived batch sorption (Cs+) or diffusion (Na+) experiments. Overall, the obtained results demonstrate the feasibility of the technical concept to study the diffusion of nonsorbing and sorbing tracers in consolidated clays. (C) 2004 Elsevier B.V. All rights reserved.
Abstract: Based on a dye tracer experiment in a sand tank we addressed the problem of local dispersion of conservative tracers in the unsaturated zone. The sand bedding was designed to have a defined spatial heterogeneity with a strong anisotropy. We estimated the parameters that characterize the local dispersion and dilution from concentration maps of a high spatial and temporal resolution obtained by image analysis. The plume spreading and mixing behavior was quantified on the basis of the coefficient of variation of the concentration and of the dilution index. The heterogeneous structure modified the flow pattern depending on water saturation. The shape of the tracer plumes revealed the structural signature of the sand bedding at low saturation only. In this case pronounced preferential flow was observed. At higher flow rates the structure remained hidden by a spatially almost homogeneous behavior of the plumes. In this context, we mainly discuss the mechanism of re-distributing a finite mass of inert solutes over a large volume, due to macro- and micro-heterogeneities of the structure. (C) 2001 Elsevier Science Ltd. AU rights reserved.
Abstract: Field soils show rather different spreading behavior at different water saturations, frequently caused by layering of the soil material. We performed tracer experiments in a laboratory sand tank. Such experiments complement and help comprehension of field investigations. We estimated, by image analysis, the first two moments of small plumes traveling through a two-dimensional, heterogeneous medium with strongly anisotropic correlation structure. Three steady state regimes were analyzed. Two main conclusions were drawn. First, low saturation led to very large heterogeneity and to strong preferential flow. Thus the description of the flow paths and the prediction of the solute arrival times require, in this case, more accurate knowledge about the topological structure. Second, saturation-dependent macroscopic anisotropy is an essential element of transport in unsaturated media. For this reason, small structural soil features should be properly upscaled to give appropriate effective soil parameters to be input in transport models.
Abstract: Geological and pedological processes rarely form isotropic media as is usually assumed in transport studies. Anisotropy at the Darcy or field scale may be detected directly by measuring flow parameters or may become indirectly evident from movement and shape of solute plumes. Anisotropic behavior of a soil at one scale may, in many cases, be related to the presence of lower-scale directional structures. Miller similitude with different pore-scale geometries of the basic element is used to model macroscopic flow and transport behavior. Analytical expressions for the anisotropic conductivity tensor are derived based on the dynamic law that governs the flow problem at the pore scale. The effects of anisotropy on transport parameters are estimated by numerical modeling.
Abstract: In many field or laboratory situations, well-mixed reservoirs like, for instance, injection or detection wells and gas distribution or sampling chambers define boundaries of transport domains. Exchange of solutes or gases across such boundaries can occur through advective or diffusive processes. First we analyzed situations, where the inlet region consists of a well-mixed reservoir, in a systematic way by interpreting them in terms of injection type. Second, we discussed the mass balance errors that seem to appear in case of resident injections. Mixing cells (MC) can be coupled mathematically in different ways to a domain where advective-dispersive transport occurs: by assuming a continuous solute flux at the interface (flux injection, MC-FI), or by assuming a continuous resident concentration (resident injection). In the latter case, the flux leaving the mixing cell can be defined in two ways: either as the value when the interface is approached from the mixing-cell side (MC-RT -), or as the value when it is approached from the column side (MC-RT +). Solutions of these injection types with constant or-in one case-distance-dependent transport parameters were compared to each other as well as to a solution of a two-layer system, where the first layer was characterized by a large dispersion coefficient. These solutions differ mainly at small Peclet numbers. For most real situations, the model for resident injection MC-RI + is considered to be relevant. This type of injection was modeled with a constant or with an exponentially varying dispersion coefficient within the porous medium. A constant dispersion coefficient will be appropriate for gases because of the Eulerian nature of the usually dominating gaseous diffusion coefficient, whereas the asymptotically growing dispersion coefficient will be more appropriate for solutes due to the Lagrangian nature of mechanical dispersion, which evolves only with the fluid flow. Assuming a continuous resident concentration at the interface between a mixing cell and a column, as in case of the MC-RI + model, entails a flux discontinuity. This flux discontinuity arises inherently from the definition of a mixing cell: the mixing process is included in the balance equation, but does not appear in the description of the flux through the mixing cell. There, only convection appears because of the homogeneous concentration within the mixing cell. Thus, the solute flux through a mixing cell in close contact with a transport domain is generally underestimated. This leads to (apparent) mass balance errors, which are often reported for similar situations and erroneously used to judge the validity of such models. Finally, the mixing cell model MC-RI + defines a universal basis regarding the type of solute injection at a boundary. Depending on the mixing cell parameters, it represents, in its limits, flux as well as resident injections. (C) 1998 Elsevier Science B.V. All rights reserved.
Abstract: Safe disposal of toxic wastes in geologic formations requires minimal water and gas movement in the vicinity of storage areas, Ventilation of repository tunnels or caverns built in solid rock can desaturate the near field up to a distance of meters from the rock surface, even when the surrounding geological formation is saturated and under hydrostatic pressures. A tunnel segment at the Grimsel test site located in the Aare granite of the Bernese Alps (central Switzerland) has been subjected to a resaturation and, subsequently, to a controlled desaturation, Using thermocouple psychrometers (TP) and time domain reflectometry (TDR), the water potentials psi and water contents theta were measured within the unsaturated granodiorite matrix near the tunnel wall at depths between 0 and 160 cm. During the resaturation the water potentials in the first 30 cm from the rock surface changed within weeks from values of less than -1.5 MPa to near saturation. They returned to the negative initial values during desaturation, The dynamics of this saturation-desaturation regime could be monitored very sensitively using the thermocouple psychrometers, The TDR measurements indicated that water contents changed dose to the surface, but at deeper installation depths the observed changes were within the experimental noise. The field-measured data of the desaturation cycle were used to test the predictive capabilities of the hydraulic parameter functions that were derived from the water retention characteristics psi(theta) determined in the laboratory. A depth-invariant saturated hydraulic conductivity k(s) = 3.0 x 10(-11) m s(-1) was estimated from the psi(t) data at all measurement depths, using the one-dimensional, unsaturated water flow and transport model HYDRUS [Vogel er al., 1996], For individual measurement depths, the estimated k(s) varied between 9.8 x 10(-12) and 6.1 x 10(-11) The fitted k(s) values fell within the range of previously estimated k(s) for this location and led to a satisfactory description of the data, even though the model did not include transport of water vapor.
Abstract: A model is developed to describe transport and loss of methyl bromide (MeBr) in soil following application as a soil fumigant. The model is used to investigate the effect of soil and management factors on MeBr volatilization. Factors studied include depth of injection, soil water content, presence or absence of tarp, depth to downward barrier, and irrigation after injection. Of these factors, the most important was irrigation after injection followed by covering with the tarp, which increased the diffusive resistance of the soil and prevented early loss of MeBr. The model offers an explanation for the apparently contradictory observations of earlier field studies of MeBr volatilization from soils under different conditions. The model was also used to calculate the concentration-time index for various management alternatives, showing that the irrigation application did not make the surface soil more difficult to fumigate, except at very early times. Therefore, irrigation shows promise for reducing fumigant loss while at the same time permitting control of target organisms during fumigation.
Abstract: The occurrence of gaseous pollutants in soils has stimulated many experimental activities, including forced ventilation in the field as well as laboratory transport experiments with gases. The dispersion coefficient in advective-dispersive gas phase transport is often dominated by molecular diffusion, which leads to a large overall dispersivity gamma. Under such conditions it is important to distinguish between flux and resident modes of solute injection and detection. The influence of the inlet type oil the macroscopic injection mode was tested in two series of column experiments with gases at different mean flow velocities nu. First we compared infinite resident and flux injections, and second, semi-infinite resident and flux injections. It is shown that the macroscopically apparent injection condition depends on the geometry of the inlet section. A reduction of the cross-sectional area of the inlet relative to that of the column is very effective in excluding the diffusive solute input, thus allowing us to use the solutions for a flux Injection also at rather low mean flow velocities nu. If the whole cross section of a column is exposed to a large reservoir like that of ambient air, a semi-infinite resident injection is established, which can be distinguished from a flux injection even at relatively high velocities nu, depending on the mechanical dispersivity of the porous medium.
Abstract: A water desaturation zone develops around a tunnel in water-saturated rock when the evaporative water loss at the rock surface is larger than the water flow from the surrounding saturated region of restricted permeability. We describe the methods with which such water desaturation processes in rock materials can be quantified. The water retention characteristic theta(psi) of crystalline rock samples was determined with a pressure membrane apparatus. The negative water potential, identical to the capillary pressure, psi, below the tensiometric range (psi < -0.1 MPa) can be measured with thermocouple psychrometers (TP), and the volumetric water contents, theta, by means of time domain reflectometry (TDR). These standard methods were adapted for measuring the water status in a macroscopically unfissured granodiorite with a total porosity of approximately 0.01. The measured water retention curve of granodiorite samples from the Grimsel test site (central Switzerland) exhibits a shape which is typical for bimodal pore size distributions. The measured bimodality is probably an artifact of a large surface ratio of solid/voids. The thermocouples were installed without a metallic screen using the cavity drilled into the granodiorite as a measuring chamber. The water potentials observed in a cylindrical granodiorite monolith ranged between -0.1 and -3.0 MPa; those near the wall in a ventilated tunnel between -0.1 and -2.2 MPa. Two types of three-rod TDR Probes were used, one as a depth probe inserted into the rock, the other as a surface probe using three copper stripes attached to the surface for detecting water content changes in the rock-to-air boundary. The TDR signal was smoothed with a low-pass filter, and the signal length determined based on the first derivative of the trace. Despite the low porosity of crystalline rock these standard methods are applicable to describe the unsaturated zone in solid rock and may also be used in other consolidated materials such as concrete.
Abstract: Transport of volatile hydrocarbons in soils is largely controlled by interactions of vapours with the liquid and solid phase. Sorption on solids of gaseous or dissolved comPounds may be important. Since the contact time between a chemical and a specific sorption site can be rather short, kinetic or mass-transfer resistance effects may be relevant. An existing mathematical model describing advection and diffusion in the gas phase and diffusional transport from the gaseous phase into an intra-aggregate water phase is modified to include linear kinetic sorption on ps-solid and water-solid interfaces. The model accounts for kinetic mass transfer between all three phases in a soil. The solution of the Laplace-transformed equations is inverted numerically. We performed transient column experiments with 1,1,2-Trichloroethane, Trichloroethylene, and Tetrachloroethylene using air-dry solid and water-saturated porous glass beads. The breakthrough curves were calculated based on independently estimated parameters. The model calculations agree well with experimental data. The different transport behaviour of the three compounds in our system primarily depends on Henryâs constants.
