Clay swelling, an important phenomenon in natural systems, can dramatically affect the properties of soils and sediments. Of particular interest in low-salinity, saturated systems are osmotic hydrates, forms of smectite in which the layer separation greatly exceeds the thickness of a single smectite layer due to the intercalation of water. In situ X-ray diffraction (XRD) studies have shown a strong link between ionic strength and average interlayer spacing in osmotic hydrates but also indicate the presence of structural disorder that has not been fully described. In the present study the structural state of expanded smectite in sodium chloride solutions was investigated by combining very low electron dose, highresolution cryogenic-transmission electron microscopy observations with XRD experiments. Wyoming smectite (SWy-2) was embedded in vitreous ice to evaluate clay structure in aqua. Lattice-fringe images showed that smectite equilibrated in aqueous, low-ionic-strength solutions, exists as individual smectite layers, osmotic hydrates composed of parallel layers, as well as disordered layer conformations. No evidence was found here for edge-to-sheet attractions, but significant variability in interlayer spacing was observed. Whether this variation could be explained by a dependence of the magnitude of long-range cohesive (van der Waals) forces on the number of layers in a smectite particle was investigated here. Calculations of the Hamaker constant for layer-layer interactions showed that van der Waals forces may span at least five layers plus the intervening water and confirmed that forces vary with layer number. Drying of the disordered osmotic hydrates induced re-aggregation of the smectite to form particles that exhibited coherent scattering domains. Clay disaggregation and restacking may be considered as an example of oriented attachment, with the unusual distinction that it may be cycled repeatedly by changing solution conditions.

Engineered barrier systems designed to isolate high-level radioactive waste utilize bentonite, a montmorillonite-rich material, to restrict contaminant transport due to its low hydraulic conductivity and high adsorption capacity. High temperatures (100–200 °C) near waste canisters resulting from radioactive decay may alter the clay's ability to adsorb contaminants. In this study, we examine U(VI) adsorption onto two different bentonite samples subjected to (1) 18 years of in situ heating during an underground experiment and (2) short-term (7-week), high temperature (300 °C) heating in the laboratory. Results show that U(VI) adsorption was lower for field-heated bentonite located closest to the heater, which experienced temperatures of approximately 95 °C, compared to a control sample, which experienced temperatures of approximately 20 °C, over a range of aqueous chemical conditions. On average, Kd values for U(VI) adsorption were 31% lower for 95 °C heated samples. By contrast, U(VI) adsorption onto intermediate-heated bentonite (50 °C) was indistinguishable from the adsorption onto the cold-zone (20 °C) sample. U(VI) adsorption onto lab-heated bentonite was also lower than onto control bentonite over the pH range 4.5–8.0, with an average decrease in Kd values of 50% after heating. Lower U(VI) adsorption to field-heated bentonite persisted after bentonite was purified to isolate the clay fraction. This allows us to rule out changes in pore-water chemistry or accessory mineral composition as causes of the lower adsorption. No evidence of montmorillonite illitization was observed in the heated samples. While some of the lower U(VI) adsorption in the lab-heated bentonite can be explained by changes in aqueous U(VI) speciation, we propose that lower U(VI) adsorption to field-heated bentonite may be primarily due to changes in the montmorillonite edge structure. The observed changes in U(VI) adsorption to bentonite after heating have implications for U(VI) diffusive transport through engineered barriers and must be considered when designing radioactive waste disposal repositories.

The reactive transport code CrunchClay was used to derive effective diffusion coefficients (D_e), clay porosities (ε), and adsorption distribution coefficients (K_D) from through-diffusion data while considering accurately the influence of unavoidable experimental biases on the estimation of these diffusion parameters. These effects include the presence of filters holding the solid sample in place, the variations in concentration gradients across the diffusion cell due to sampling events, the impact of tubing/dead volumes on the estimation of diffusive fluxes and sample porosity, and the effects of O-ring-filter setups on the delivery of solutions to the clay packing. Doing so, the direct modeling of the measurements of (radio)tracer concentrations in reservoirs is more accurate than that of data converted directly into diffusive fluxes. While the above-mentioned effects have already been described individually in the literature, a consistent modeling approach addressing all these issues at the same time has never been described nor made easily available to the community. A graphical user interface, CrunchEase, was created, which supports the user by automating the creation of input files, the running of simulations, and the extraction and comparison of data and simulation results. While a classical model considering an effective diffusion coefficient, a porosity and a solid/solution distribution coefficient (D_e-ε-K_D) may be implemented in any reactive transport code, the development of CrunchEase makes it easy to apply by experimentalists without a background in reactive transport modeling. CrunchEase makes it also possible to transition more easily from a D_e-ε-K_D modeling approach to a state-of-the-art process-based understanding modeling approach using the full capabilities of CrunchClay, which include surface complexation modeling and a multi-porosity description of the clay packing with charged diffuse layers.

The dissolution of CO2 in water leads to a pH decrease and a carbonate content increase in affected groundwater, which in turn can drive the mobilization of metals from sediments. The mechanisms of metal release postulated in various field and laboratory studies often differ. Drawing primarily on previously published results, we examine contrasting metal mobilization behaviors at two field tests and in one laboratory study, to investigate whether the same mechanisms could explain metal releases in these different experiments. Numerical modeling of the two field tests reveals that fast Ca-driven cation exchange (from calcite dissolution) can explain the release of most major and trace metal cations at both sites, and their parallel concentration trends. The dissolution of other minerals reacting more slowly (superimposed on cation exchange) also contributes to metal release over longer time frames, but can be masked by fast ambient groundwater velocities. Therefore, the magnitude and extent of mobilization depends not only on metal-mineral associations and sediment pH buffering characteristics, but also on groundwater flow rates, thus on the residence time of CO2-impacted groundwater relative to the rates of metal-release reactions. Sequential leaching laboratory tests modeled using the same metal-release concept as postulated from field experiments show that both field and laboratory data can be explained by the same processes. The reversibility of metal release upon CO2 degassing by de-pressurization is also explored using simple geochemical models, and shows that the sequestration of metals by resorption and re-precipitation upon CO2 exsolution is quite plausible and may warrant further attention.

One of the risks that CO2 geological sequestration imposes on the environment is the impact of potential CO2/brine leakage on shallow groundwater. The reliability of reactive transport models predicting the response of groundwater to CO2 leakage depends on a thorough understanding of the relevant chemical processes and key parameters affecting dissolved CO2 transport and reaction. Such understanding can be provided by targeted field tests integrated with reactive transport modeling. A controlled-release field experiment was conducted in Mississippi to study the CO2-induced geochemical changes in a shallow sandy aquifer at about 50 m depth. The field test involved a dipole system in which the groundwater was pumped from one well, saturated with CO2 at the pressure corresponding to the hydraulic pressure of the aquifer, and then re-injected into the same aquifer using a second well. Groundwater samples were collected for chemical analyses from four monitoring wells before, during and after the dissolved CO2 was injected. In this paper, we present reactive transport models used to interpret the observed changes in metal concentrations in these groundwater samples. A reasonable agreement between simulated and measured concentrations indicates that the chemical response in the aquifer can be interpreted using a conceptual model that encompasses two main features: (a) a fast-reacting but limited pool of reactive minerals that responds quickly to changes in pH and causes a pulse-like concentration change, and (b) a slow-reacting but essentially unlimited mineral pool that yields rising metal concentrations upon decreased groundwater velocities after pumping and injection stopped. During the injection, calcite dissolution and Ca-driven cation exchange reactions contribute to a sharp pulse in concentrations of Ca, Ba, Mg, Mn, K, Li, Na and Sr, whereas desorption reactions control a similar increase in Fe concentrations. After the injection and pumping stops and the groundwater flow rate decreases, the dissolution of relatively slow reacting minerals such as plagioclase drives the rising concentrations of alkali and alkaline earth metals observed at later stages of the test, whereas the dissolution of amorphous iron sulfide causes slowly increasing Fe concentrations.