Long-term geochemical evolution of the near field repository: insights from reactive transport modelling and experimental evidences

The KBS-3 underground nuclear waste repository concept designed by the Swedish Nuclear Fuel and Waste Management Co. (SKB) includes a bentonite buffer barrier surrounding the copper canisters and the iron insert where spent nuclear fuel will be placed. Bentonite is also part of the backfill material used to seal the access and deposition tunnels of the repository. The bentonite barrier has three main safety functions: to ensure the physical stability of the canister, to retard the intrusion of groundwater to the canisters, and in case of canister failure, to retard the migration of radionuclides to the geosphere. Laboratory experiments (< 10 years long) have provided evidence of the control exerted by accessory minerals and clay surfaces on the pore water chemistry. The evolution of the pore water chemistry will be a primordial factor on the long-term stability of the bentonite barrier, which is a key issue in the safety assessments of the KBS-3 concept.In this work we aim to study the long-term geochemical evolution of bentonite and its pore water in the evolving geochemical environment due to climate change. In order to do this, reactive transport simulations are used to predict the interaction between groundwater and bentonite which is simulated following two different pathways: (1) groundwater flow through the backfill in the deposition tunnels, eventually reaching the top of the deposition hole, and (2) direct connection between groundwater and bentonite rings through fractures in the granite crosscutting the deposition hole. The influence of changes in climate has been tested using three different waters interacting with the bentonite: present-day groundwater, water derived from ice melting, and deep-seated brine. Two commercial bentonites have been considered as buffer material, MX-80 and Deponit CA-N, and one natural clay (Friedland type) for the backfill. They show differences in the composition of the exchangeable cations and in the accessory mineral content. Results from the simulations indicate that pore water chemistry is controlled by the equilibrium with the accessory minerals, especially carbonates. pH is buffered by precipitation/dissolution of calcite and dolomite, when present. The equilibrium of these minerals is deeply influenced by gypsum dissolution and cation exchange reactions in the smectite interlayer. If carbonate minerals are initially absent in bentonite, pH is then controlled by surface acidity reactions in the hydroxyl groups at the edge sites of the clay fraction, although its buffering capacity is not as strong as the equilibrium with carbonate minerals. The redox capacity of the bentonite pore water system is mainly controlled by Fe(II)-bearing minerals (pyrite and siderite). Changes in the groundwater composition lead to variations in the cation exchange occupancy, and dissolution–precipitation of carbonate minerals and gypsum. The most significant changes in the evolution of the system are predicted when ice-melting water, which is highly diluted and alkaline, enters into the system. In this case, the dissolution of carbonate minerals is enhanced, increasing pH in the bentonite pore water. Moreover, a rapid change in the population of exchange sites in the smectite is expected due to the replacement of Na for Ca.     

Improved interpretation of in-diffusion measurements with confined swelling clays

Diffusion is considered the principal transport mechanism of radio-nuclides and other low-molecular-weight pollutants in compacted clays used as barriers at various disposal and storage sites, for example, at projected deep repositories for radioactive waste. Porous filters are routinely used to confine swelling clays in diffusion studies of radio-tracers. The presence of the filter gives rise to considerable mass-transfer limitations at the clay boundary that result in erroneous diffusion parameters. We have solved the problem of in-diffusion with due account for this phenomenon by means of Fourier transforms. By using literature data on the in-diffusion of traces of radioactive cesium in an argillaceous rock (Opalinus clay) and a compacted bentonite (FEBEX bentonite), we have demonstrated that taking into account the mass-transfer limitations considerably improves the quality of the theoretical fit of the time evolution of radio-tracer concentration in the reservoir. Besides that, we have shown that ignoring the mass-transfer limitations leads to a noticeable underestimation of both the effective diffusion coefficient and the specific sorption capacity of the clay.     

Impact of drying and composting procedures on the concentrations of 4-nonylphenols, di-(2-ethylhexyl)phthalate and polychlorinated biphenyls in anaerobically digested sewage sludge

Enhanced treatments of sewage sludge produce a more manageable product for agricultural use by stabilizing the material, removing water, and reducing the possibility of pathogen transfer. We investigated the impact of pilot-scale composting and drying of sludge on physicochemical characteristics and on the concentrations of some organic contaminants. During the 143 day composting procedure, organic matter fell 22% and moisture by half. Concentrations of 4-nonylphenols (4-NPs) fell by 88% and di-(2-ethylhexyl) phthalate (DEHP) by 60%; losses continued throughout the procedure. Losses of total polychlorinated biphenyls (PCBs) were 11%, mostly from the lower molecular weight congeners, suggesting volatilization as the most likely loss mechanism. The drying process was much shorter, 40 days, yet organic matter content decreased by 27% and moisture by 85%. Losses of 4-NPs (39%) and DEHP (22%) were less than in composting and stopped when moisture content became constant. There were no losses of PCBs. Both treatments are simple, practical procedures that reduce the volume of waste and are applicable in situ on farms. Composting would be the method of choice for reducing organic contaminants but requires much longer times than drying.     

Inverse modeling of multicomponent reactive transport through single and dual porosity media

Compacted bentonite is foreseen as buffer material for high-level radioactive waste in deep geological repositories because it provides hydraulic isolation, chemical stability, and radionuclide sorption. A wide range of laboratory tests were performed within the framework of FEBEX (Full-scale Engineered Barrier EXperiment) project to characterize buffer properties and develop numerical models for FEBEX bentonite. Here we present inverse single and dual-continuum multicomponent reactive transport models of a long-term permeation test performed on a 2.5 cm long sample of FEBEX bentonite. Initial saline bentonite porewater was flushed with 5.5 pore volumes of fresh granitic water. Water flux and chemical composition of effluent waters were monitored during almost 4 years. The model accounts for solute advection and diffusion and geochemical reactions such as aqueous complexation, acid-base, cation exchange, protonation/deprotonation by surface complexation and dissolution/precipitation of calcite, chalcedony and gypsum. All of these processes are assumed at local equilibrium. Similar to previous studies of bentonite porewater chemistry on batch systems which attest the relevance of protonation/deprotonation on buffering pH, our results confirm that protonation/deprotonation is a key process in maintaining a stable pH under dynamic transport conditions. Breakthrough curves of reactive species are more sensitive to initial porewater concentration than to effective diffusion coefficient. Optimum estimates of initial porewater chemistry of saturated compacted FEBEX bentonite are obtained by solving the inverse problem of multicomponent reactive transport. While the single-continuum model reproduces the trends of measured data for most chemical species, it fails to match properly the long tails of most breakthrough curves. Such limitation is overcome by resorting to a dual-continuum reactive transport model.     

Kinetics of metribuzin release from bentonite-polymer composites in water

A series of bentonite polymer-composites (BPCs) loaded with metribuzin were studied for their controlled release in aqueous medium. The release of active ingredient from BPCs was significantly lower as compared to commercial metribuzin formulation. The results revealed that the cumulative metribuzin release was highest (81%) from the BPCs containing 8% clay (commercial bentonite) and 2% metribuzin which correspond to the lowest (14 days) half-life values i.e., time required for 50% release of active ingredient (t_1/2). The metribuzin release from the BPCs decreased with increased concentration of clays in polymer matrix and the release was further decreased with BPCs prepared with pure nano-bentonite. BPCs containing 12% clay and 2% metribuzin showed maximum t_1/2 values i.e., 25 and 51 days for commercial bentonite and pure nano-bentonite as clay sources, respectively. The differential behaviour in the metribuzin release rates from BPCs was ascribed due to variations in crosslinking of metribuzin in the composites. As metribuzin release was found to be slower in BPCs compared to commercial formulation, it could be used for control of weeds tailored to different crops.     

Effect of temperature on the containment properties of argillaceous rocks: The case study of Callovo-Oxfordian claystones

Heat generated by high level radioactive wastes could alter the performance of a clay repository. It was intended to investigate the effect of such a thermal period on the diffusive properties of Callovo-Oxfordian claystones. Thus, through-diffusion experiments with HTO, Cl-36, Na-22 and Cs-137 were performed before, during and after stages of heating at 80°C that lasted for up to one year. A special attention was paid to limit the occurrence of any chemical disturbance. Therefore (i) the temperature was raised to 80°C, then progressively brought back to 21°C, thanks to three intermediate temperature stages, and (ii) specific synthetic solutions were used for each temperature, chemistry of which being close to the equilibrium state, especially with respect to the carbonate and sulphate minerals. It was found that experiments carried out at 80°C showed a clear increase of the effective diffusion coefficient values for the four tracers with respect to those obtained at 21°C (by a factor of 3 for HTO and Cl-36, 5 for Na-22 and 2 for Cs-137). On the other hand, the porosity and rock capacity values did not exhibit any significant discrepancy between 21°C and 80°C, indicating no observable damage of both the pore conducing network and the sorption properties of clay minerals. The Stokes-Einstein relationship, based on the temperature dependency of the viscosity of bulk water, could be used to describe the temperature dependence of the diffusion of HTO and Cl-36 but failed to describe the diffusive evolution of the two sorbing cations, Na-22 and Cs-137. Furthermore, experiments performed after the thermal period led to diffusive properties well matching those obtained before heating. All these results suggest that at the lab scale the heating of rock samples would not alter the claystone containment properties.     

Non-aqueous-phase-liquid breakthrough during evaporative drying of clay barriers

In this study, an attempt has been made to model a real field scenario, whereby an initially almost saturated clay liner in a waste site is gradually drying, due to evaporation at its lower boundary. A detailed conceptual model that deals with the penetration and breakthrough of non-aqueous-phase-liquid (NAPL) in clay liners is introduced. Water content of clay samples was monitored during ambient evaporation through apertures at the base of sample holders. Clay drying rate served as the primary parameter for the NAPL breakthrough study. The interconnection between drying rates, structural damage formation (cracks and suction) and NAPL penetration is especially addressed. The processes taking place in the clay samples during drying appear to be associated with the capillary effects between the different fluid phases in the vicinity of either the NAPL-clay or the clay-air boundaries. A conceptual model of NAPL penetration and breakthrough of the clay layer has been considered, based on both indirect and direct observations of structural damages produced on either clay boundaries. A mutual interaction between these two boundaries is suggested and discussed. NAPL breakthrough is suggested to take place through cracks initiated on the upper soil surface.     

An integrated sorption-diffusion model for the calculation of consistent distribution and diffusion coefficients in compacted bentonite

A thermodynamic sorption model and a diffusion model based on electric double layer (EDL) theory are integrated to yield a surface chemical model that treats porewater chemistry, surface reactions, and the influence of charged pore walls on diffusing ions in a consistent fashion. The relative contribution of Stern and diffuse layer to the compensation of the permanent surface charge represents a key parameter; it is optimized for the diffusion of Cs in Kunipia-F bentonite, at a dry density of 400 kg/m3. The model is then directly used to predict apparent diffusivities (Da) of Cs, Sr, Cl-, I- and TcO4- and corresponding distribution coefficients (Kd) of Cs and Sr in different bentonites as a function of dry density, without any further adjustment of surface chemical and EDL parameters. Effective diffusivities (De) for Cs, HTO, and TcO4- are also calculated. All calculated values (Da, De, Kd) are fully consistent with each other. A comparison with published, measured data shows that the present model allows a good prediction and consistent explanation of (i) apparent and effective diffusivities for cations, anions, and neutral species in compacted bentonite, and of (ii) Kd values in batch and compacted systems.     

Bacterial deposition in unsaturated porous media as related to surface properties

Transport of Escherichia coli, Pseudomonas fluorescens and Bacillus subtilis in silica sand under water-unsaturated conditions was investigated using column experiments. It was hypothesised that bacterial deposition was due to their interactions within the pore environment, which was a function of their surface physicochemical properties as well as pore water chemistry. Surface thermodynamic properties of these bacterial strains were measured independently by means of contact angle measurements under different water saturation conditions using variable lawn moisture contents. Bacterial interactions with the liquid-gas interface and the porous media were calculated based on their surface properties and were related to their transport observations.     

Simultaneous sorption of lead and chlorobenzene by organobentonite

Clays or organoclays have been used as a barrier to prevent the transport of hazardous contaminants in landfills. However, clays are known to effectively sorb mostly inorganic contaminants, while organoclays are mainly used for organic contaminants. Since the organoclays are basically clay particles modified with cationic surfactants, there might exist an optimal coverage of cationic surfactant on the clay particles to sorb both inorganic and organic contaminants. In order to determine the optimal mass of cationic surfactants on the bentonites, sodium bentonites were treated with various ratios of hexadecyltrimethylammonium (HDTMA) to bentonites. Chlorobenzene and lead were selected as representative contaminants. When either chlorobenzene or lead exists as a single contaminant, chlorobenzene sorption increased with increasing HDTMA to bentonite ratios, and lead sorption decreased with increasing HDTMA to bentonite ratios. Sorption of chlorobenzene was a function of HDTMA coverage on the bentonites, while lead sorption was much more influenced by the initial lead concentration rather than the mass of HDTMA added to the bentonites.     

Climate change impacts on the leaching of a heavy metal contamination in a small lowland catchment

The Keersop catchment (43km(2)) in the south of The Netherlands has been contaminated by the emissions of four zinc ore smelters. The objective of this study was to assess the effects of future projected climate change on the hydrology and the leaching of heavy metals (i.e. Cd and Zn) in the catchment. The numerical, quasi-2D, unsaturated zone Soil Water Atmosphere Plant model was used with 100-year simulated daily time series of precipitation and potential evapotranspiration. The time series are representative of stationary climates for the periods 1961-1990 ("baseline") and 2071-2100 ("future"). The time series of future climate were obtained by downscaling the results of eight regional climate model (RCM) experiments, driven by the SRES A2 emissions scenario, using change factors for a series of climate statistics and applying them to stochastic weather generator models. The time series are characterized by increased precipitation in winter, less precipitation in summer, and higher air temperatures (between 2°C and 5°C) throughout the year. Future climate scenarios project higher evapotranspiration rates, more irrigation, less drainage, lower discharge rates and lower groundwater levels, due to increased evapotranspiration and a slowing down of the groundwater system. As a result, lower concentrations of Cd and Zn in surface water are projected. The reduced leaching of heavy metals, due to drying of the catchment, showed a positive impact on a limited aspect of surface water quality.     

A coupled THC model of the FEBEX in situ test with bentonite swelling and chemical and thermal osmosis

The performance assessment of a geological repository for radioactive waste requires quantifying the geochemical evolution of the bentonite engineered barrier. This barrier will be exposed to coupled thermal (T), hydrodynamic (H), mechanical (M) and chemical (C) processes. This paper presents a coupled THC model of the FEBEX (Full-scale Engineered Barrier EXperiment) in situ test which accounts for bentonite swelling and chemical and thermal osmosis. Model results attest the relevance of thermal osmosis and bentonite swelling for the geochemical evolution of the bentonite barrier while chemical osmosis is found to be almost irrelevant. The model has been tested with data collected after the dismantling of heater 1 of the in situ test. The model reproduces reasonably well the measured temperature, relative humidity, water content and inferred geochemical data. However, it fails to mimic the solute concentrations at the heater-bentonite and bentonite-granite interfaces because the model does not account for the volume change of bentonite, the CO(2)(g) degassing and the transport of vapor from the bentonite into the granite. The inferred HCO(3)(-) and pH data cannot be explained solely by solute transport, calcite dissolution and protonation/deprotonation by surface complexation, suggesting that such data may be affected also by other reactions.     

Experimental and modeling studies on sorption and diffusion of radium in bentonite

The sorption and desorption behavior of radium on bentonite and purified smectite was investigated as a function of pH, ionic strength and liquid to solid ratio by batch experiments. The distribution coefficients (Kd) were in the range of 10(2) to > 10(4) ml g-1 and depended on ionic strength and pH. Most of sorbed Ra was desorbed by 1 M KCl. The results for purified smectite indicated that Ra sorption is dominated by ion exchange at layer sites of smectite, and surface complexation at edge sites may increase Ra sorption at higher pH region. Reaction parameters between Ra and smectite were determined based on an interaction model between smectite and groundwater. The reaction parameters were then used to explain the results of bentonite by considering dissolution and precipitation of minerals and soluble impurities. The dependencies of experimental Kd values on pH, ionic strength and liquid to solid ratio were qualitatively explained by the model. The modeling result for bentonite indicated that sorption of Ra on bentonite is dominated by ion exchange with smectite. The observed pH dependency was caused by changes of Ca concentration arising from dissolution and precipitation of calcite. Diffusion behavior of Ra in bentonite was also investigated as a function of dry density and ionic strength. The apparent diffusion coefficients (Da) obtained in compacted bentonite were in the range of 1.1 x 10(-11) to 2.2 x 10(-12) m2 s-1 and decreased with increasing in dry density and ionic strength. The Kd values obtained by measured effective diffusion coefficient (De) and modeled De were consistent with those by the sorption model in a deviation within one order of magnitude.     

Scoping analyses of geochemical sealing of early cracks in a waste container and associated drip shield, Yucca Mountain, Nevada

Early after final emplacement of the nuclear waste containers at the proposed Yucca Mountain, Nevada, high-level-waste repository, tiny cracks (less than 200 microm wide, 1 to 2 cm deep, and a few centimeters long at most) could appear in the containers and in the drip shield protecting them. Modeling calculations were performed to understand how fast those cracks could be sealed. Under dripping conditions, they are expected to be bridged with water. If cracks are located in the drip shield, any further dripping on the waste containers located underneath will be limited. If cracks are located in a container, potentially harmful radionuclides could only travel by diffusion. In addition, water-bridged cracks will be sealed through at least two processes: precipitation of calcite with minor silica following evaporative concentration of the water residing in the cracks and continuous corrosion of the crack walls. The sealing rate is calculated as the intersection of the time of emergence of the cracks, the water dripping rate, and the decreasing evaporation rate. The evaporative driving force declines as short-lived radioactive elements, having given up much of the heat affecting the repository, are progressively depleted from the waste. Depending on the crack initiation time and environmental conditions, crack sealing varies from a few tens of years to a few thousand years. Because environmental conditions in the vicinity of the cracks and at the crack scale have not been produced, a parametric method scaling drift scale conditions is used.     

Controlled release and retarded leaching of pesticides by encapsulating in carboxymethyl chitosan /bentonite composite gel

A novel composite gel composed of carboxymethyl-chitosan (CM-chit) and bentonite (H-bent) was used as the carrier for encapsulating atrazine and imidacloprid to control their release in water and retard their leaching in soil. Strong interactions between CM-chit and H-bent in the composite were confirmed by FT-IR, and good dispersion of pesticides in the carrier was observed by SEM. According to the results of release experiments in water, the CM-chit/H-bent composite carrier showed double advantages of both encapsulation by the polymer and sorption by the bentonite. The time taken for 50 % of active ingredients to be released, t??, was prolonged to 572 h for atrazine and 24 h for imidacloprid, respectively. The difference between the two pesticides on release behavior was related to their hydrophobicity and water solubility. Leaching experiments through a soil layer showed that this novel carrier reduced the amount of pesticides available for leaching, and would be useful for diminishing the environmental pollution of pesticides.     

Competitive adsorption of alkylbenzene and water vapors on predominantly mineral surfaces

Competitive adsorption of ethylbenzene (EB) and water on bentonite and of p-xylene (pXYL) and water on kaolin and silica gel was studied using a technique that allowed the amount of adsorbed water and alkylbenzene to be measured independently. EB adsorption on bentonite was not affected by water at relative humidities (RH) near 0.23 but was reduced significantly at RH's near 0.50. pXYL adsorption on kaolin and silica gel decreased with increasing RH, especially above a RH of about 0.2. Increasing RH not only decreased the amount of alkylbenzene adsorption but also resulted in a change from Type-II adsorption isotherms to ones that were essentially linear. Linear isotherms for the adsorption of hydrophobic organic compounds on hydrated soil have generally been attributed to partitioning into organic carbon (OC). However, since the clays and oxide used here had very low to trace amounts of OC, it is suggested that processes involving only mineral surfaces can give rise to linear isotherms. Based on solubility considerations alone, partitioning of EB and pXYL into adsorbed water films was not considered to be an important adsorption mechanism in this study. The effect of cation hydration on the amount of water adsorbed from a mixture of water and pXYL vapors was evaluated by comparing adsorption on Li- and Na-saturated kaolin.     

Diffusion of radionuclides in concrete and concrete-bentonite systems

Various construction materials are under consideration for nuclear waste repositories. Two important materials are concrete and bentonite clay, which will act as mechanical barriers and prevent convective water flow. These barriers will also retard transport (diffusion controlled) of dissolved radionuclides by a combination of mechanical constraints and chemical interactions with the solid.An important issue is the possible change of the initial sodium bentonite into the calcium form due to interaction with calcium from the concrete. The initial leaching of concrete was studied using radioactive spiked concrete in contact with compacted bentonite.Measurement were made of the diffusion of Cs, Am and Pu into 5 different types of concrete in contact with pore water. The diffusivity measured for Cs agrees reasonably well with data found in the literature. No movement could be measured for Am and Pu (< 0.2 mm), even though the contact times were extremely long (2.5 and 5 yr, respectively). The diffusion of Na, Ca and Cs from concrete into bentonite was also measured.     

Sensitivity analysis of radionuclide migration in compacted bentonite: a mechanistic model approach

Mechanistic model calculations for the migration of Cs, Ra, Am and Pb in compacted bentonite have been carried out to evaluate sensitivities with respect to different parameter variations. A surface chemical speciation/electric double layer model is used to calculate: (i) porewater composition and radionuclide speciation in solution and at the bentonite surface, yielding the distribution of mobile and sorbed species and (ii) interaction of diffusing species with negatively charged pore walls to obtain diffusion parameters. The basic scenario considers the interaction of compacted bentonite with a fresh-type groundwater; variations include the presence of bentonite impurities and saline groundwater. It is shown that these scenarios result in significant variations of porewater composition that affect migration via three mechanisms that can partly compensate each other: (1) effects on sorption through radionuclide complexation in solution, and competition of major cations for surface sites; (2) changes in radionuclide solution speciation leading to different diffusing species under different conditions; (3) effects on diffusion through changes in the electric double layer properties of the clay pores as a function of ionic strength.     

Impact of clay mineral on air oxidation of PAH-contaminated soils

This work investigated the impact of a clay mineral (bentonite) on the air oxidation of the solvent extractable organic matters (EOMs) and the PAHs from contaminated soils. EOMs were isolated from two coking plant soils and mixed with silica sand or bentonite. These samples, as well as raw soils and bentonite/soil mixtures, were oxidized in air at 60 and 100 °C for 160 days. Mineralization was followed by measuring the CO₂ produced over the experiments. EOM, polycyclic aromatic compound (PAC), including PAH, contents were also determined. Oxidation led to a decrease in EOM contents and PAH concentrations, these diminutions were enhanced by the presence of bentonite. Transfer of carbon from EOM to insoluble organic matter pointed out a condensation phenomenon leading to a stabilization of the contamination. Higher mineralization rates, observed during the oxidation of the soil/bentonite mixtures, seem to indicate that this clay mineral had a positive influence on the transformation of PAC into CO₂.