Dynamics of supercooled water in highly compacted clays studied by neutron scattering

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.

Multi-scale micro-structure generation strategy for up-scaling transport in clays

Clays and claystones are used as backfill and barrier materials in the design of waste repositories, because they act as hydraulic barriers and retain contaminants. Transport through such barriers occurs mainly by molecular diffusion. There is thus an interest to relate the diffusion properties of clays to their structural properties. In previous work, we have developed a concept for up-scaling pore-scale molecular diffusion coefficients using a grid-based model for the sample pore structure. Here we present an operational algorithm which can generate such model pore structures of polymineral materials. The obtained pore maps match the rock’s mineralogical components and its macroscopic properties such as porosity, grain and pore size distributions. Representative ensembles of grains in 2D or 3D are created by a lattice Monte Carlo (MC) method, which minimizes the interfacial energy of grains starting from an initial grain distribution. Pores are generated at grain boundaries and/or within grains. The method is general and allows to generate anisotropic structures with grains of approximately predetermined shapes, or with mixtures of different grain types. A specific focus of this study was on the simulation of clay-like materials. The generated clay pore maps were then used to derive upscaled effective diffusion coefficients for non-sorbing tracers using a homogenization technique. The large number of generated maps allowed to check the relations between micro-structural features of clays and their effective transport parameters, as is required to explain and extrapolate experimental diffusion results. As examples, we present a set of 2D and 3D simulations and investigated the effects of nanopores within particles (interlayer pores) and micropores between particles. Archie’s simple power law is followed in systems with only micropores. When nanopores are present, additional parameters are required; the data reveal that effective diffusion coefficients could be described by a sum of two power functions, related to the micro- and nanoporosity. We further used the model to investigate the relationships between particle orientation and effective transport properties of the sample.

Effect of Nontronite Smectite Clay on the Chemical Evolution of Several Organic Molecules under Simulated Martian Surface Ultraviolet Radiation Conditions

Most of the phyllosilicates detected at the surface of Mars today are probably remnants of ancient environments that sustained long-term bodies of liquid water at the surface or subsurface and were possibly favorable for the emergence of life. Consequently, phyllosilicates have become the main mineral target in the search for organics on Mars. But are phyllosilicates efficient at preserving organic molecules under current environmental conditions at the surface of Mars? We monitored the qualitative and quantitative evolutions of glycine, urea, and adenine in interaction with the Fe3+-smectite clay nontronite, one of the most abundant phyllosilicates present at the surface of Mars, under simulated martian surface ultraviolet light (190-400 nm), mean temperature (218 +/- 2 K), and pressure (6 +/- 1 mbar) in a laboratory simulation setup. We tested organic-rich samples that were representative of the evaporation of a small, warm pond of liquid water containing a high concentration of organics. For each molecule, we observed how the nontronite influences its quantum efficiency of photodecomposition and the nature of its solid evolution products. The results reveal a pronounced photoprotective effect of nontronite on the evolution of glycine and adenine; their efficiencies of photodecomposition were reduced by a factor of 5 when mixed at a concentration of 2.6x10(-2) mol of molecules per gram of nontronite. Moreover, when the amount of nontronite in the sample of glycine was increased by a factor of 2, the gain of photoprotection was multiplied by a factor of 5. This indicates that the photoprotection provided by the nontronite is not a purely mechanical shielding effect but is also due to stabilizing interactions. No new evolution product was firmly identified, but the results obtained with urea suggest a particular reactivity in the presence of nontronite, leading to an increase of its dissociation rate. Key Words: Martian surface-Organic chemistry-Photochemistry-Astrochemistry-Nontronite-Phyllosilicates. Astrobiology 15, 221-237.

Field investigation of dried lakes in western United States as an analogue to desiccation fractures on Mars

Potential Desiccation Polygons (PDPs), tens to hundreds of meters in size, have been observed in numerous regions on Mars, particularly in ancient (>3Gyr old) terrains of inferred paleolacustrine/playa geologic setting, and in association with hydrous minerals such as smectites. Therefore, a better understanding of the conditions in which large desiccation polygons form could yield unique insight into the ancient climate on Mars. Many dried lakebeds/playas in western United States display large (>50m wide) desiccation polygons, which we consider to be analogues for PDPs on Mars. Therefore, we have carried out fieldwork in seven of these dried lakes in San Bernardino and the Death Valley National Park regions complemented with laboratory and spectral analysis of collected samples. Our study shows that the investigated lacustrine/playa sediments have (a) a soil matrix containing 40-75% clays and fine silt (by volume) where the clay minerals are dominated by illite/muscovite followed by smectite, (b) carbonaceous mineralogy with variable amounts of chloride and sulfate salts, and significantly, (c) roughly similar spectral signatures in the visible-near-infrared (VIS-NIR) range. We conclude that the development of large desiccation fractures is consistent with water table retreat. In addition, the comparison of the mineralogical to the spectral observations further suggests that remote sensing VIS-NIR spectroscopy has its limitations for detailed characterization of lacustrine/playa deposits. Finally, our results imply that the widespread distribution of PDPs on Mars indicates global or regional climatic transitions from wet conditions to more arid ones making them important candidate sites for future in situ missions.

Natural tracer profiles across argillaceous formations

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.

Analysis of polygonal cracking patterns in chloride-bearing terrains on Mars: Indicators of ancient playa settings

The ancient southern highlands on Mars (~3.5 Gyr old) contain > 600 regions that display spectral evidence in the infrared for the presence of chloride-bearing materials. Many of these locations were previously reported to display polygonal cracking patterns. We studied more than 80 of the chloride-bearing terrains using high-resolution (0.25-0.5 m/pixel) images, as well as near-infrared spectral data, to characterize the surface textures and the associated cracking patterns and mineralogies. Our study indicates that ~75% of the studied locations display polygonal cracks that resemble desiccation cracks, while some resemble salt expansion/thrust polygons. Furthermore, we detect, spectrally, the presence of smectites in association with ~30% of the studied fractured terrains. We note that smectites are a special class of swelling clay minerals that can induce formation of large desiccation cracks. As such, we suggest that the cracking patterns are indicative of the presence of smectite phyllosilicates even in the absence of spectral confirmation. Our results suggest that many chloride-bearing terrains have a lacustrine origin and a geologic setting similar to playas on Earth. Such locations would have contained ephemeral lakes that may have undergone repeated cycles of desiccation and recharging by a near-surface fluctuating water table in order to account for the salt-phyllosilicates associations. These results have notable implications for the ancient hydrology of Mars. We propose that the morphologies and sizes of the polygonal cracks can be used as paleoenvironmental, as well as lithological, indicators that could be helpful in planning future missions.

Geochemical synthesis for the Effingen Member in boreholes at Oftringen, Gösgen and Küttigen

The Effingen Member is a low-permeability rock unit of Oxfordian age (ca. 160 Ma) that occurs across northern Switzerland. It comprises sandy calcareous marls and (argillaceous) limestones. This report describes the hydrogeochemistry, mineralogy and supporting physical properties of the Effingen Member in three boreholes in the Jura-Südfuss area: Oftringen, Gösgen and Küttigen, where it is 220–240 m thick. The top of the Effingen Member is at 420, 66 and 32 m depths at the three sites. Core materials are available from Oftringen and Gösgen, whereas information from Küttigen is limited to cuttings, in-situ hydrogeological testing and geophysical logging. Hydrogeological boundaries of the Effingen Member vary between locations. Ground-water flows were identified during drilling at the top (Geissberg Member), but not at the base, of the Effingen Member at Oftringen, at the base (Hauptrogenstein Formation) of the Effingen Member at Gösgen, and in a limestone layer (Gerstenhübel unit) within the Effingen Member at Küttigen. The marls and limestones of the Effingen Member have carbonate contents of 46–91 wt.-% and clay-mineral contents of 5–37 wt.-%. Pyrite contents are up to 1.6 wt.-%, but no sulphate minerals were detected by routine analyses. Clay minerals are predominantly mixed-layer illite-smectite, illite and kaolinite, with sporadic traces of chlorite and smectite. Veins filled with calcite ± celestite occur through the Effingen Member at Oftringen but not at Gösgen or Küttigen. They formed at 50–70 ºC from externally derived fluids, probably of Miocene age. Water contents are 0.7–4.2 wt.-%, corresponding to a water-loss porosity range of 1.9–10.8 vol.-%. Specific surface areas, measured by the BET method, are 2–30 m2/g, correlating with clay-mineral contents. Water activity has been measured and yielded surprisingly low values down to 0.8. These cannot be explained by pore-water salinity alone and include other effects, such as changes in the fabric due to stress release or partial saturation. Observed variations in measurements are not fully understood. Cation exchange capacity (CEC) and exchangeable cation populations have been studied by the Ni-en method. CEC, derived from the consumption of the index cation Ni, is 9–99 meq/kgrock at a solid:liquid ratio of 1, correlating with the clay-mineral content. Cation concentrations in Ni-en extract solutions are in the order Na+≥Ca2+>Mg2+>K+>Sr2+. However, the analytical results from the Ni-en extractions have additional contributions from cations originating from pore water and from mineral dissolution reactions that occurred during extraction, and it was not possible to reliably quantify these contributions. Therefore, in-situ cation populations and selectivity coefficients could not be derived. A suite of methods have been used for characterising the chemical compositions of pore waters in the Effingen Member. Advective displacement was used on one sample from each Oftringen and Gösgen and is the only method that produces results that approach complete hydrochemical compositions. Aqueous extraction was used on core samples from these two boreholes and gives data only for Cl- and, in some cases, Br-. Out-diffusion was used on core samples from Oftringen and similarly gives data for Cl- and Br- only. For both aqueous extraction and out-diffusion, reaction of the experimental water with rock affected concentrations of cations, SO42 and alkalinity in experimental solutions. Another method, centrifugation, failed to extract pore water. Stable isotope ratios (δ18O and δ2H) of pore waters in core samples from Oftringen were analysed by the diffusive exchange method and helium contents of pore water in Oftringen samples were extracted for mass spectrometric analysis by quantitative outgassing of preserved core samples. Several lines of evidence indicate that drillcore samples might not have been fully saturated when opened and subsampled in the laboratory. These include comparisons of water-loss porosities with physical porosities, water-activity measurements, and high contents of dissolved gas as inferred from ground-water samples. There is no clear proof of partial saturation and it is unclear whether this might represent in-situ conditions or is due to exsolution of gas due to the pressure release since drilling. Partial saturation would have no impact on the recalculation of pore-water compositions from aqueous extraction experiments using water-loss porosity data. The largest uncertainty in the pore-water Cl- concentrations recalculated from aqueous extraction and out-diffusion experiments is the magnitude of the anion-accessible fraction of water-loss porosity. General experience of clay-mineral rich formations suggests that the anion-accessible porosity fraction is very often about 0.5 and generally in a range of 0.3 to 0.6 and tends to be inversely correlated with clay-mineral contents. Comparisons of the Cl- concentration in pore water obtained by advective displacement with that recalculated from aqueous extraction of an adjacent core sample suggests a fraction of 0.27 for an Oftringen sample, whereas the same procedure for a Gösgen sample suggests a value of 0.64. The former value for anion-accessible porosity fraction is presumed to be unrepresentative given the local mineralogical heterogeneity at that depth. Through-diffusion experiments with HTO and 36Cl- suggest that the anion-accessible porosity fraction in the Effingen Member at Oftringen and Gösgen is around 0.5. This value is proposed as a typical average for rocks of the Effingen Member, bearing in mind that it varies on a local scale in response to the heterogeneity of lithology and pore-space architecture. The substantial uncertainties associated with the approaches to estimating anion-accessible porosity propagate into the calculated values of in-situ pore-water Cl- concentrations. On the basis of aqueous extraction experiments, and using an anion-accessible porosity fraction of 0.5, Cl- concentrations in the Effingen Member at Oftringen reach a maximum of about 14 g/L in the centre. Cl- decreases upwards and downwards from that, forming a curved depth profile. Cl- contents in the Effingen Member at Gösgen increase with depth from about 3.5 g/L to about 14 g/L at the base of the cored profile (which corresponds to the centre of the formation). Out-diffusion experiments were carried out on four samples from Oftringen, distributed through the Effingen Member. Recalculated Cl- concentrations are similar to those from aqueous extraction for 3 out of the 4 samples, and somewhat lower for one sample. Concentrations of other components, i.e. Na+, K+, Ca2+, Mg2+, Sr2+, SO42- and HCO3- cannot be obtained from the aqueous extraction and out-diffusion experimental data because of mineral dissolution and cation exchange reactions during the experiments. Pore-water pH also is not constrained by those extraction experiments. The only experimental approach to obtain complete pore-water compositions for samples from Oftringen and Gösgen is advective displacement of pore water. The sample from Oftringen used for this experiment is from 445 m depth in the upper part of the Effingen Member and gave eluate with 16.5 g/L Cl- whereas aqueous extraction from a nearby sample indicated about 9 g/L Cl-. The sample from Gösgen used for advective displacement is from 123 m depth in the centre of the Effingen Member sequence and gave eluate with about 9 g/L Cl- whereas aqueous extraction gave 11.5 g/L Cl-. In both cases the pore waters have Na-(Ca)-Cl compositions and SO42- concentrations of about 1.1 g/L. The Gösgen sample has a Br/Cl ratio similar to that of sea water, whereas this ratio is lower for the Oftringen sample. Taking account of uncertainties in the applied experimental approaches, it is reasonable to place an upper limit of ca. 20 g/L on Cl- concentration for pore water in the Effingen Member in this area. There are major discrepancies between pore-water SO42- concentrations inferred from aqueous extraction or out-diffusion experiments and those obtained from advective displacement in both the Oftringen and Gösgen cases. A general conclusion is that all or at least part of the discrepancies are attributable to perturbation of the sulphur system and enhancement of SO42- by sulphate mineral dissolution and possibly minor pyrite oxidation during aqueous extraction and out-diffusion. Therefore, data for SO42- calculated from those pore-water sampling methods are considered not to be representative of in-situ conditions. A reference pore-water composition was defined for the Effingen Member in the Jura Südfuss area. It represents the probable upper limits of Cl- contents and corresponding anion and cation concentrations that are reasonably constrained by experimental data. Except for Cl- and possibly Na+ concentrations, this composition is poorly constrained especially with respect to SO42- and Ca2+ concentrations, and pH and alkalinity. Stable isotope compositions, δ18O and δ2H, of pore waters in the Effingen Member at Oftringen plot to the right of the meteoric water line, suggesting that 18O has been enriched by water-rock exchange, which indicates that the pore waters have a long residence time. A long residence time of pore water is supported by the level of dissolved 4He that has accumulated in pore water of the Effingen Member at Oftringen. This is comparable with, or slightly higher than, the amounts of 4He in the Opalinus Clay at Benken. Ground waters were sampled from flowing zones intersected by boreholes at the three locations. The general interpretation is that pore waters and ground-water solutes may have similar origins in Mesozoic and Cenozoic brackish-marine formations waters, but ground-water solutes have been diluted rather more than pore waters by ingress of Tertiary and Quaternary meteoric waters. The available hydrochemical data for pore waters from the Effingen Member at these three locations in the Jura-Südfuss area suggest that the geochemical system evolved slowly over geological periods of time, in which diffusion was an important mechanism of solute transport. The irregularity of Cl- and δ18O profiles and spatial variability of advective ground-water flows in the Malm-Dogger system suggests that palaeohydrogeological and hydrochemical responses to changing tectonic and surface environmental conditions were complex.

Translational diffusion of water in compacted clay systems

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.