Transferability of geoscientific information from various sources (study sites, underground rock laboratories, natural analogues) to support safety cases for radioactive waste repositories in argillaceous formations

In studies related to deep geological disposal of radioactive waste, it is current practice to transfer external information (e.g. from other sites, from underground rock laboratories or from natural analogues) to safety cases for specific projects. Transferable information most commonly includes parameters, investigation techniques, process understanding, conceptual models and high-level conclusions on system behaviour. Prior to transfer, the basis of transferability needs to be established. In argillaceous rocks, the most relevant common feature is the microstructure of the rocks, essentially determined by the properties of clay–minerals. Examples are shown from the Swiss and French programmes how transfer of information was handled and justified. These examples illustrate how transferability depends on the stage of development of a repository safety case and highlight the need for adequate system understanding at all sites involved to support the transfer.

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).

Experimental Characterization and Quantification of Cement-Bentonite Interaction using Core Infiltration Techniques Coupled with X-ray Tomography

Deep geological storage of radioactive waste foresees cementitious materials as reinforcement of tunnels and as backfill. Bentonite is proposed to enclose spent fuel canisters and as drift seals. Sand/bentonite (s/b) is foreseen as backfill material of access galleries or as drift seals. The emplacement of cementitious material next to clay material generates an enormous chemical gradient in pore-water composition that drives diffusive solute transport. Laboratory studies and reactive transport modeling predicted significant mineral alteration at and near interfaces, mainly resulting in a decrease of porosity in bentonite. The goal of this thesis was to characterize and quantify the cement/bentonite interactions both spatially and temporally in laboratory experiments. A newly developed mobile X-ray transparent core infiltration device was used to perform X-ray computed tomography (CT) scans without interruption of running experiments. CT scans allowed tracking the evolution of the reaction plume and changes in core volume/diameter/density during the experiments. In total 4 core infiltration experiments were carried out for this study with the compacted and saturated cores consisting of MX-80 bentonite and sand/MX-80 bentonite mixture (s/b; 65/35%). Two different high-pH cementitious pore-fluids were infiltrated: a young (early) ordinary Portland cement pore-fluid (APWOPC; K+–Na+–OH-; pH 13.4; ionic strength 0.28 mol/kg) and a young ‘low-pH’ ESDRED shotcrete pore-fluid (APWESDRED; Ca2+–Na+–K+–formate; pH 11.4; ionic strength 0.11 mol/kg). The experiments lasted between 1 and 2 years. In both bentonite experiments, the hydraulic conductivity was strongly reduced after switching to high-pH fluids, changing eventually from an advective to a diffusion-dominated transport regime. The reduction was mainly induced by mineral precipitation and possibly partly also by high ionic strength pore-fluids. Both bentonite cores showed a volume reduction and a resulting transient flow in which pore-water was squeezed out during high-pH infiltration. The outflow chemistry was characterized by a high ionic strength, while chloride in the initial pore water got replaced as main anionic charge carrier by sulfate, originating from gypsum dissolution. The chemistry of the high-pH fluids got strongly buffered by the bentonite, consuming hydroxide and in case of APWESDRED also formate. Hydroxide got consumed by mineral reactions (saponite and possibly talc and brucite precipitation), while formate being affected by bacterial degradation. Post-mortem analysis showed reaction zones near the inlet of the bentonite core, characterized by calcium and magnesium enrichment, consisting predominately of calcite and saponite, respectively. Silica got enriched in the outflow, indicating dissolution of silicate-minerals, identified as preferentially cristobalite. In s/b, infiltration of APWOPC reduced the hydraulic conductivity strongly, while APWESDRED infiltration had no effect. The reduction was mainly induced by mineral precipitation and probably partly also by high ionic strength pore-fluids. Not clear is why the observed mineral precipitates in the APWESDRED experiment had no effect on the fluid flow. Both s/b cores showed a volume expansion along with decreasing ionic strengths of the outflow, due to mineral reactions or in case of APWESDRED infiltration also mediated by microbiological activity, consuming hydroxide and formate, respectively. The chemistry of the high-pH fluids got strongly buffered by the s/b. In the case of APWESDRED infiltration, formate reached the outflow only for a short time, followed by enrichment in acetate, indicating most likely biological activity. This was in agreement to post-mortem analysis of the core, observing black spots on the inflow surface, while the sample had a rotten-egg smell indicative of some sulfate reduction. Post-mortem analysis showed further in both cores a Ca-enrichment in the first 10 mm of the core due to calcite precipitation. Mg-enrichment was only observed in the APWOPC experiment, originating from newly formed saponite. Silica got enriched in the outflow of both experiments, indicating dissolution of silicate-minerals, identified in the OPC experiment as cristobalite. The experiments attested an effective buffering capacity for bentonite and s/b, a progressing coupled hydraulic-chemical sealing process and also the preservation of the physical integrity of the interface region in this setup with a total pressure boundary condition on the core sample. No complete pore-clogging was observed but the hydraulic conductivity got rather strongly reduced in 3 experiments, explained by clogging of the intergranular porosity (macroporosity). Such a drop in hydraulic conductivity may impact the saturation time of the buffer in a nuclear waste repository, although the processes and geometry will be more complex in repository situation.

Interaction of corroding iron with bentonite in the ABM1 experiment at Äspö, Sweden: A microscopic approach

Bentonite and iron metals are common materials proposed for use in deep-seated geological repositories for radioactive waste. The inevitable corrosion of iron leads to interaction processes with the clay which may affect the sealing properties of the bentonite backfill. The objective of the present study was to improve our understanding of this process by studying the interface between iron and compacted bentonite in a geological repository-type setting. Samples of MX-80 bentonite samples which had been exposed to an iron source and elevated temperatures (up to 115ºC) for 2.5 y in an in situ experiment (termed ABM1) at the Äspö Hard Rock Laboratory, Sweden, were investigated by microscopic means, including scanning electron microscopy, μ-Raman spectroscopy, spatially resolved X-ray diffraction, and X-ray fluorescence. The corrosion process led to the formation of a ~100 mm thick corrosion layer containing siderite, magnetite, some goethite, and lepidocrocite mixed with the montmorillonitic clay. Most of the corroded Fe occurred within a 10 mm-thick clay layer adjacent to the corrosion layer. An average corrosion depth of the steel of 22–35 μm and an average Fe2+ diffusivity of 1–26×10–13 m2/s were estimated based on the properties of the Fe-enriched clay layer. In that layer, the corrosion-derived Fe occurred predominantly in the clay matrix. The nature of this Fe could not be identified. No indications of clay transformation or newly formed clay phases were found. A slight enrichment of Mg close to the Fe–clay contact was observed. The formation of anhydrite and gypsum, and the dissolution of some SiO2 resulting from the temperature gradient in the in situ test, were also identified. © 2014, Clay Minerals Society. All right reserved.

Experimental characterization of cement–bentonite interaction using core infiltration techniques and 4D computed tomography

Deep geological storage of radioactive waste foresees cementitious materials as reinforcement of tunnels and as backfill. Bentonite is proposed to enclose spent fuel drums, and as drift seals. The emplacement of cementitious material next to clay material generates an enormous chemical gradient in pore water composition that drives diffusive solute transport. Laboratory studies and reactive transport modeling predict significant mineral alteration at and near interfaces, mainly resulting in a decrease of porosity in bentonite. The goal of this project is to characterize and quantify the cement/bentonite skin effects spatially and temporally in laboratory experiments. A newly developed mobile X-ray transparent core infiltration device was used, which allows performing X-ray computed tomography (CT) periodically without interrupting a running experiment. A pre-saturated cylindrical MX-80 bentonite sample (1920 kg/m3 average wet density) is subjected to a confining pressure as a constant total pressure boundary condition. The infiltration of a hyperalkaline (pH 13.4), artificial OPC (ordinary Portland cement) pore water into the bentonite plug alters the mineral assemblage over time as an advancing reaction front. The related changes in X-ray attenuation values are related to changes in phase densities, porosity and local bulk density and are tracked over time periodically by non-destructive CT scans.

Paleozoic-aged brine and authigenic helium preserved in an Ordovician shale aquiclude

Consideration of the geosphere for isolation of nuclear waste has generated substantial interest in the origin, age, and movement of fl uids and gases in low-permeability rock formations. Here, we present profi les of isotopes, solutes, and helium in porewaters recovered from 860 m of Cambrian to Devonian strata on the eastern fl ank of the Michigan Basin. Of particular interest is a 240-m-thick, halite-mineralized, Ordovician shale and carbonate aquiclude, which hosts Br–-enriched, post-dolomitic brine (5.8 molal Cl) originating as evaporated Silurian seawater. Authigenic helium that has been accumulating in the aquiclude for more than 260 m.y. is found to be isolated from underlying allochthonous, 3He-enriched helium that originated from the rifted base of the Michigan Basin and the Canadian Shield. The Paleozoic age and immobility of the pore fl uids in this Ordovician aquiclude considerably strengthen the safety case for deep geological repositories, but also provide new insights into the origin of deep crustal brines and opportunities for research on other components of a preserved Paleozoic porewater system.

Dissolved CO2 and Alkane Gas in Clay Formations

Characterization of dissolved CO2 and alkane gas in clayrocks may help assessing the confinement properties of geological barriers considered as potential host rocks for a deep geological disposal as well as for caprocks of gas storages. A monitoring of alkanes with CO2, combined with carbon isotopes was performed on core samples coming from Underground Research Laboratories (Bure, Mont Terri, Tournemire) and the Schlattingen borehole in France and Switzerland. Composition of hydrocarbon gas and delta C-13 of methane strongly suggest a dominant thermogenic origin of methane which is mixed with a bacterial origin for the Toarcian shales, Pliensbachien and Callovian-Oxfordian clayrocks. Results also evidence the contrasted behavior of CO2, which is controlled by chemical equilibrium between pore water and carbonate mineralogy, compared to the alkanes which are present in the porosity as a stock of dissolved gases which can be depleted during degassing experiments. (C) 2015 The Authors. Published by Elsevier B.V.

Diffusion-driven transport in clayrock formations

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 corresponding 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 diffusion-driven transfer in these formations shows extreme differences in the measured or modelled behaviour for various radionuclides, e. g. between halogen RN (Cl-36, I-129) and actinides (U-238,U-235, Np-237, Th-232, etc.), which result from major differences between RN of the effects on transport of two phenomena: 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 behaviour observed experimentally, e.g. D-e(anion) < D-e(HTO) < D-e(cation) and D-e(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-interacting cations (e.g. Sr2+) and still relatively poor for 'strongly sorbing', high K-d cations. Funmig results have clarified understanding of how clayrock mineral composition, and the corresponding 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 materials 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 coupled diffusion-sorption migration is much more complex than expected, leading generally to greater mobility than that predicted by coupling a batch-determined K-d and Ficks law based on the diffusion behaviour of HTO. If the K-d 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 D-e 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. (C) 2011 Elsevier Ltd. All rights reserved.

Ediacaran 2,500-km-long synchronous deep continental subduction in the West Gondwana Orogen

The deeply eroded West Gondwana Orogen is a major continental collision zone that exposes numerous occurrences of deeply subducted rocks, such as eclogites. The position of these eclogites marks the suture zone between colliding cratons, and the age of metamorphism constrains the transition from subduction-dominated tectonics to continental collision and mountain building. Here we investigate the metamorphic conditions and age of high-pressure and ultrahigh-pressure eclogites from Mali, Togo and NE-Brazil and demonstrate that continental subduction occurred within 20 million years over at least a 2,500-km-long section of the orogen during the Ediacaran. We consider this to be the earliest evidence of large-scale deep-continental subduction and consequent appearance of Himalayan-scale mountains in the geological record. The rise and subsequent erosion of such mountains in the Late Ediacaran is perfectly timed to deliver sediments and nutrients that are thought to have been necessary for the subsequent evolution of sustainable life on Earth.

Global pulses of organic carbon burial in deep-sea sediments during glacial maxima

The burial of organic carbon in marine sediments removes carbon dioxide from the ocean–atmosphere pool, provides energy to the deep biosphere, and on geological timescales drives the oxygenation of the atmosphere. Here we quantify natural variations in the burial of organic carbon in deep-sea sediments over the last glacial cycle. Using a new data compilation of hundreds of sediment cores, we show that the accumulation rate of organic carbon in the deep sea was consistently higher (50%) during glacial maxima than during interglacials. The spatial pattern and temporal progression of the changes suggest that enhanced nutrient supply to parts of the surface ocean contributed to the glacial burial pulses, with likely additional contributions from more efficient transfer of organic matter to the deep sea and better preservation of organic matter due to reduced oxygen exposure. These results demonstrate a pronounced climate sensitivity for this global carbon cycle sink.

Geological and hydrological histories of the Argyre province, Mars

The geologic history of the multi-ringed Argyre impact basin and surroundings has been reconstructed on the basis of geologic mapping and relative-age dating of rock materials and structures. The impact formed a primary basin, rim materials, and a complex basement structural fabric including faults and valleys that are radial and concentric about the primary basin, as well as structurally-controlled local basins. Since its formation, the basin has been a regional catchment for volatiles and sedimentary materials as well as a dominant influence on the flow of surface ice, debris flows, and groundwater through and over its basement structures. The basin is interpreted to have been occupied by lakes, including a possible Mediterranean-sized sea that formed in the aftermath of the Argyre impact event The hypothesized lakes froze and diminished through time, though liquid water may have remained beneath the ice cover and sedimentation may have continued for some time. At its deepest, the main Argyre lake may have taken more than a hundred thousand years to freeze to the bottom even absent any heat source besides the Sun, but with impact-induced hydrothermal heat, geothermal heat flow due to long-lived radioactivities in early martian history, and concentration of solutes in sub-ice brine, liquid water may have persisted beneath thick ice for many millions of years. Existence of an ice-covered sea perhaps was long enough for life to originate and evolve with gradually colder and more hypersaline conditions. The Argyre rock materials, diverse in origin and emplacement mechanisms, have been modified by impact, magmatic, eolian, fluvial, lacustrine, glacial, periglacial, alluvial, colluvial, and tectonic processes. Post-impact adjustment of part of the impact-generated basement structural fabric such as concentric faults is apparent. Distinct basin-stratigraphic units are interpreted to be linked to large-scale geologic activity far from the basin, including growth of the Tharsis magmatic-tectonic complex and the growth into southern middle latitudes of south polar ice sheets. Along with the migration of surface and sub-surface volatiles towards the central part of the primaiy basin, the substantial difference in elevation with respect to the surrounding highlands and Tharsis and the Thaumasia highlands result in the trapping of atmospheric volatiles within the basin in the form of fog and regional or local precipitation, even today. In addition, the impact event caused long-term (millions of years) hydrothermal activity, as well as deep-seated basement structures that have tapped the internal heat of Mars, as conduits, for far greater time, possibly even today. This possibility is raised by the observation of putative open-system pingos and nearby gullies that occur in linear depressions with accompanying systems of faults and fractures. Long-term water and heat energy enrichment, complemented by the interaction of the nutrient-enriched primordial crustal and mantle materials favorable to life excavated to the surface and near-surface environs through the Argyre impact event, has not only resulted in distinct geomorphology, but also makes the Argyre basin a potential site of exceptional astrobiological significance. (C) 2015 Elsevier Inc. All rights reserved.

Scientific drilling projects in ancient lakes: Integrating geological and biological histories

Sedimentary sequences in ancient or long-lived lakes can reach several thousands of meters in thickness and often provide an unrivalled perspective of the lake's regional climatic, environmental, and biological history. Over the last few years, deep-drilling projects in ancient lakes became increasingly multi- and interdisciplinary, as, among others, seismological, sedimentological, biogeochemical, climatic, environmental, paleontological, and evolutionary information can be obtained from sediment cores. However, these multi- and interdisciplinary projects pose several challenges. The scientists involved typically approach problems from different scientific perspectives and backgrounds, and setting up the program requires clear communication and the alignment of interests. One of the most challenging tasks, besides the actual drilling operation, is to link diverse datasets with varying resolution, data quality, and age uncertainties to answer interdisciplinary questions synthetically and coherently. These problems are especially relevant when secondary data, i.e., datasets obtained independently of the drilling operation, are incorporated in analyses. Nonetheless, the inclusion of secondary information, such as isotopic data from fossils found in outcrops or genetic data from extant species, may help to achieve synthetic answers. Recent technological and methodological advances in paleolimnology are likely to increase the possibilities of integrating secondary information. Some of the new approaches have started to revolutionize scientific drilling in ancient lakes, but at the same time, they also add a new layer of complexity to the generation and analysis of sediment-core data. The enhanced opportunities presented by new scientific approaches to study the paleolimnological history of these lakes, therefore, come at the expense of higher logistic, communication, and analytical efforts. Here we review types of data that can be obtained in ancient lake drilling projects and the analytical approaches that can be applied to empirically and statistically link diverse datasets to create an integrative perspective on geological and biological data. In doing so, we highlight strengths and potential weaknesses of new methods and analyses, and provide recommendations for future interdisciplinary deep-drilling projects.