Reactive transport modelling as a toolbox to compare remediation strategies for aquifers impacted by uranium in situ recovery.

More than 60% of worldwide uranium production is based on the In Situ Recovery mining technique. This exploitation method directly falls within the scope of the applications of reactive transport modelling to optimize uranium production and limit its associated environmental impact. We propose a modelling approach which is able to represent the natural evolution of an aquifer impacted by an ISR test performed using sulfuric acid. The model is calibrated on a 12 year-long data series obtained from 12 monitoring wells surrounding an ISR pilot cell. Through this process-based approach, we simulate the impact of several remediation strategies that can be considered in these contexts. In particular, we model the impact of Pump & Treat combined with reverse osmosis, as well as the circulation of non-impacted fluids through the reservoir with different operating strategies. Our approach allows to compare the effectiveness of these strategies. For this small-scale ISR pilot, monitored natural attenuation constitutes an interesting approach due to its faster pH recovery time with respect to Pump & Treat (5-10 years to pH ∼ 6), whose efficiency can be improved by the addition of exchangeable cations. Circulation of unimpacted fluids can reduce pH recovery times if performed for periods longer than the ISR exploitation and/or deployed with a delay. Combined with an economic evaluation of their deployment, this modelling approach can help the mining operator select and design optimal remediation strategies from an environmental and economical standpoint.

Influence of clay minerals on pH and major cation concentrations in acid-leached sands: Column experiments and reactive-transport modeling.

A series of laboratory experiments are conducted to simulate the acidification and subsequent recovery of a sand aquifer exploited by in situ recovery (ISR) mining. A sulfuric acid solution (pH 2) is first injected into a column packed with sand from the Zoovch Ovoo uranium roll front deposit (Mongolia). Solutions representative of local groundwater or enriched in cations (Na+, Mg2+) are then circulated through the column to simulate the inflow of aquifer water. pH and major ion concentrations (Na+, Cl-, SO42-, Ca2+, Mg2+, K+) measured at the column outlet reproduce the overall evolution of porewater chemistry observed in the field. The presence of minor quantities of swelling clay minerals (≈6 wt% smectite) is shown to exert an important influence on the behavior of inorganic cations, particularly H+, via ion-exchange reactions. Numerical models that consider ion-exchange on smectite as the sole solid-solution interaction are able to reproduce variations in pH and cation concentrations in the column experiments. This highlights the importance of clay minerals in controlling H+ mobility and demonstrates that sand from the studied aquifer can be described to a first order as an ion-exchanger. The present study confirms the key role of clay minerals in controlling water chemistry in acidic environments through ion-exchange processes. In a context of managing the long-term environmental footprint of industrial and mining activities (ISR, acid mine drainage…), this work will bring insights for modeling choices and identification of key parameters to help operators to define their production and/or remediation strategies.

Chemical and isotopic analyses of hydrocarbon-bearing fluid inclusions in olivine-rich rocks.

We examined the mineralogical, chemical and isotopic compositions of secondary fluid inclusions in olivine-rich rocks from two active serpentinization systems: the Von Damm hydrothermal field (Mid-Cayman Rise) and the Zambales ophiolite (Philippines). Peridotite, troctolite and gabbroic rocks in these systems contain abundant CH4-rich secondary inclusions in olivine, with less abundant inclusions in plagioclase and clinopyroxene. Olivine-hosted secondary inclusions are chiefly composed of CH4 and minor H2, in addition to secondary minerals including serpentine, brucite, magnetite and carbonates. Secondary inclusions in plagioclase are dominated by CH4 with variable amounts of H2 and H2O, while those in clinopyroxene contain only CH4. We determined hydrocarbon abundances and stable carbon isotope compositions by crushing whole rocks and analysing the released volatiles using isotope ratio monitoring-gas chromatography mass spectrometry. Bulk rock gas analyses yielded appreciable quantities of CH4 and C2H6 in samples from Cayman (4-313 nmol g-1 CH4 and 0.02-0.99 nmol g-1 C2H6), with lesser amounts in samples from Zambales (2-37 nmol g-1 CH4 and 0.004-0.082 nmol g-1 C2H6). Mafic and ultramafic rocks at Cayman exhibit δ13CCH4 values of -16.7‰ to -4.4‰ and δ13CC2H6 values of -20.3‰ to +0.7‰. Ultramafic rocks from Zambales exhibit δ13CCH4 values of -12.4‰ to -2.8‰ and δ13CC2H6 values of -1.2‰ to -0.9‰. Similarities in the carbon isotopic compositions of CH4 and C2H6 in plutonic rocks, Von Damm hydrothermal fluids, and Zambales gas seeps suggest that leaching of fluid inclusions may provide a significant contribution of abiotic hydrocarbons to deep-sea vent fluids and ophiolite-hosted gas seeps. Isotopic compositions of CH4 and C2H6 from a variety of hydrothermal fields hosted in olivine-rich rocks that are similar to those in Von Damm vent fluids further support the idea that a significant portion of abiotic hydrocarbons in ultramafic-influenced vent fluids is derived from fluid inclusions. This article is part of a discussion meeting issue 'Serpentinite in the Earth system'.

Abiotic methane synthesis and serpentinization in olivine-hosted fluid inclusions.

The conditions of methane (CH4) formation in olivine-hosted secondary fluid inclusions and their prevalence in peridotite and gabbroic rocks from a wide range of geological settings were assessed using confocal Raman spectroscopy, optical and scanning electron microscopy, electron microprobe analysis, and thermodynamic modeling. Detailed examination of 160 samples from ultraslow- to fast-spreading midocean ridges, subduction zones, and ophiolites revealed that hydrogen (H2) and CH4 formation linked to serpentinization within olivine-hosted secondary fluid inclusions is a widespread process. Fluid inclusion contents are dominated by serpentine, brucite, and magnetite, as well as CH4(g) and H2(g) in varying proportions, consistent with serpentinization under strongly reducing, closed-system conditions. Thermodynamic constraints indicate that aqueous fluids entering the upper mantle or lower oceanic crust are trapped in olivine as secondary fluid inclusions at temperatures higher than ∼400 °C. When temperatures decrease below ∼340 °C, serpentinization of olivine lining the walls of the fluid inclusions leads to a near-quantitative consumption of trapped liquid H2O. The generation of molecular H2 through precipitation of Fe(III)-rich daughter minerals results in conditions that are conducive to the reduction of inorganic carbon and the formation of CH4 Once formed, CH4(g) and H2(g) can be stored over geological timescales until extracted by dissolution or fracturing of the olivine host. Fluid inclusions represent a widespread and significant source of abiotic CH4 and H2 in submarine and subaerial vent systems on Earth, and possibly elsewhere in the solar system.