Water–montmorillonite systems: Neutron scattering and tracer throughdiffusion studies

Designs for deep geological respositories of nuclear waste include bentonite as a hydraulic and chemisorption buffer material to protect the biosphere from leakage of radionuclides. Bentonite is chosen because it is a cheap, naturally occurring material with the required properties. It consists essentially of montmorillonite, a swelling clay mineral. Upon contact with groundwater such clays can seal the repository by incorporating water in the interlayers of their crystalline structure. The intercalated water exhibits significantly different properties to bulk water in the surrounding interparticle pores, such as lower diffusion coefficients (González Sánchez et. al. 2008). This doctoral thesis presents water distribution and diffusion behavior on various time and space scales in montmorillonite. Experimental results are presented for Na- and Cs-montmorillonite samples with a range of bulk dry densities (0.8 to 1.7 g/cm3). The experimental methods employed were neutron scattering (backscattering, diffraction, time-of-flight), adsorption measurements (water, nitrogen) and tracer-through diffusion. For the tracer experiments the samples were fully saturated via the liquid phase under volume-constrained conditions. In contrast, for the neutron scattering experiments, the samples were hydrated via the vapor phase and subsequently compacted, leaving a significant fraction of interparticle pores unfilled with water. Owing to these differences in saturation, the water contents of the samples for neutron scattering were characterized by gravimetry whereas those for the tracer experiments were obtained from the bulk dry density. The amount of surface water in interlayer pores could be successfully discriminated from the amount of bulk-like water in interparticle pores in Na- and Csmontmorillonite using neutron spectroscopy. For the first time in the literature, the distribution of water between these two pore environments was deciphered as a function of gravimetric water content. The amount was compared to a geometrical estimation of the amount of interlayer and interparticle water determined by neutron diffraction and adsorption measurements. The relative abundances of the 1 to 4 molecular water layers in the interlayer were determined from the area ratios of the (001)-diffraction peaks. Depending on the characterization method, different fractions of surface water and interlayer water were obtained. Only surface and interlayer water exists in amontmorillonite with water contents up to 0.18 g/g according to spectroscopic measurements and up to 0.32 g/g according to geometrical estimations, respectively. At higher water contents, bulk-like and interparticle water also exists. The amounts increase monotonically, but not linearly, from zero to 0.33 g/g for bulk-like water and to 0.43 g/g for interparticle water. It was found that water most likely redistributes between the surface and interlayer sites during the spectroscopic measurements and therefore the reported fraction is relevant only below about -10 ºC (Anderson, 1967). The redistribution effect can explain the discrepancy in fractions between the methods. In a novel approach the fractions of water in different pore environments were treated as a fixed parameter to derive local diffusion coefficients for water from quasielastic neutron scattering data, in particular for samples with high water contents. Local diffusion coefficients were obtained for the 1 to 4 molecular water layers in the interlayer of 0.5·10–9, 0.9·10–9, 1.5·10–9 and 1.4·10–9 m²/s, respectively, taking account of the different water fractions (molecular water layer, bulk-like water). The diffusive transport of 22Na and HTO through Na-montmorillonite was measured on the laboratory experimental scale (i.e. cm, days) by tracer through-diffusion experiments. We confirmed that diffusion of HTO is independent of the ionic strength of the external solution in contact with the clay sample but dependent on the bulk dry density. In contrast, the diffusion of 22Na was found to depend on both the ionic strength of the pore solution and on the bulk dry density. The ratio of the pore and surface diffusion could be experimentally determined for 22Na from the dependence of the diffusion coefficient on the ionic strength. Activation energies were derived from the temperaturedependent diffusion coefficients via the Arrhenius relation. In samples with high bulk dry density the activation energies are slightly higher than those of bulk water whereas in low density samples they are lower. The activation energies as a function of ionic strengths of the pore solutions are similar for 22Na and HTO. The facts that (i) the slope of the logarithmic effective diffusion coefficients as a function of the logarithmic ionic strength is less than unity for low bulk dry densities and (ii) two water populations can be observed for high gravimetric water contents (low bulk dry densities) support the interlayer and interparticle porosity model proposed by Glaus et al. (2007), Bourg et al. (2006, 2007) and Gimmi and Kosakowski (2011).

The mobility of Nb in rutile-saturated NaCl- and NaF-bearing aqueous fluids from 1-6.5 GPa and 300-800 °C

Rutile (TiO2) is an important host phase for high field strength elements (HFSE) such as Nb in metamorphic and subduction zone environments. The observed depletion of Nb in arc rocks is often explained by the hypothesis that rutile sequesters HFSE in the subducted slab and overlying sediment, and is chemically inert with respect to aqueous fluids evolved during prograde metamorphism in the forearc to subarc environment. However, field observations of exhumed terranes, and experimental studies, indicate that HFSE may be soluble in complex aqueous fluids at high pressure (i.e., >0.5 GPa) and moderate to high temperature (i.e., >300 degrees C). In this study, we investigated experimentally the mobility of Nb in NaCl- and NaF-bearing aqueous fluids in equilibrium with Nb-bearing rutile at pressure-temperature conditions applicable to fluid evolution in arc environments. Niobium concentrations in aqueous fluid at rutile saturation were measured directly by using a hydrothermal diamond-anvil cell (HDAC) and synchrotron X-ray fluorescence (SXRF) at 2.1 to 6.5 GPa and 300-500 degrees C, and indirectly by performing mass loss experiments in a piston-cylinder (PC) apparatus at similar to 1 GPa and 700-800 degrees C. The concentration of Nb in a 10 wt% NaCl aqueous fluid increases from 6 to 11 mu g/g as temperature increases from 300 to 500 degrees C, over a pressure range from 2.1 to 2.8 GPa, consistent with a positive temperature dependence. The concentration of Nb in a 20 wt% NaCl aqueous fluid varies from 55 to 150 mu g/g at 300 to 500 degrees C, over a pressure range from 1.8 to 6.4 GPa; however, there is no discernible temperature or pressure dependence. The Nb concentration in a 4 wt% NaF-bearing aqueous fluid increases from 180 to 910 mu g/g as temperature increases from 300 to 500 degrees C over the pressure range 2.1 to 6.5 GPa. The data for the F-bearing fluid indicate that the Nb content of the fluid exhibits a dependence on temperature between 300 and 500 degrees C at >= 2 GPa, but there is no observed dependence on pressure. Together, the data demonstrate that the hydrothermal mobility of Nb is strongly controlled by the composition of the fluid, consistent with published data for Ti. At all experimental conditions, however, the concentration of Nb in the fluid is always lower than coexisting rutile, consistent with a role for rutile in moderating the Nb budget of arc rocks.

Viscous-Flow Approach to In Situ Infiltration and In Vitro Saturated Hydraulic Conductivity Determination

Infiltration is dominantly gravity driven, and a viscous-flow approach was developed. Laminar film flow equilibrates gravity with the viscous force and a constant flow velocity evolves during a period lasting 3/2 times the duration of a constant input rate, qS. Film thickness F and the specific contact area L of the film per unit soil volume are the key parameters. Sprinkler irrigation produced in situ time series of volumetric water contents, θ(z,t), as determined with TDR probes. The wetting front velocity v and the time series of the mobile water content, w(z,t) were deduced from θ(z,t). In vitro steady flow in a core of saturated soil provided volume flux density, q(z,t), and flow velocity, v, as determined from a heat front velocity. The F and L parameters of the in situ and the in vitro experiments were compared. The macropore-flow restriction states that, for a particular permeable medium, the specific contact area L must be independent from qS i.e., dL/dqS = 0. If true, then the relationship of qS ∝ v3/2 could scale a wide range of input rates 0 ≤ qS ≤ saturated hydraulic conductivity, Ksat, into a permeable medium, and kinematic-wave theory would become a versatile tool to deal with non-equilibrium flow. The viscous-flow approach is based on hydromechanical principles similar to Darcy’s law, but currently it is not suited to deduce flow properties from specified individual spatial structures of permeable media.

REE and istotopic composition of ODP Site 127-797 clays

X-ray diffraction analyses have been carried out on 128 samples of Miocene to Quaternary sediments from ODP Sites 794, 795 and 797. Some clay fractions of samples from Site 797 have also been studied for rare earth elements and by Nd isotopic analyses. These three sites display similar lithological and clay assemblages (with dominant chlorite, illite and smectite) showing that the sedimentation was homogeneous throughout the whole Japan Sea Basin. Three mineralogical zones are recognized. The first zone (Lower Miocene sandy clay of Sites 794 and 797) is mainly composed of chlorite resulting from hydrothermal transformation of arc-derived smectite, due to sill injections during the initial oceanic spreading stage. The second zone (Lower Miocene to Lower Pliocene siliceous claystone and diatomaceous silty clay) is dominated by arc-derived smectite; the abundance of this mineral decreases upwards while illite and chlorite increase. This trend reflects a change of detrital source, from an eastern arc-derived source (epsilon -Nd**t>-3.3); variable LREE enrichment) to a western continental crust source (epsilon-Nd**t<-9.4; shale-like REE patterns); climatic modifications in the current dynamics are proposed as a cause for this change. The third zone (Upper Pliocene to Recent silty clay with minor diatom oozes) is characterized at Site 797 by increasing amounts of illite and chlorite. This reflects a more and more important western supply which is assumed to be related to tectonic rejuvenations of the Asian margin or climatic modifications affecting the alteration conditions or the current dynamics. At Sites 794 and 795, the more or less sharp supply of chlorite seems to be driven by the incipient subduction zone on the eastern margin of the Japan Sea.

Geochemistry and mineralogy of interstitial waters and clays in sediments of ODP Leg 180 sites

Oxygen and strontium isotopes and Rb and Ba were determined in interstitial water (IW) collected from Sites 1109, 1115, and 1118 drilled on the Woodlark Rise during Ocean Drilling Program Leg 180. The trace element and mineralogical composition of the clay fraction of sediments isolated from the squeeze cakes corresponding to IW samples from Site 1109 was also determined.