Reactive transport modeling of column experiments on the evolution of saline-alkaline waste solutions

Leakage of saline-alkaline tank waste solutions often creates a serious environmental contamination problem. To better understand the mechanisms controlling the fate of such waste solutions in the Hanford vadose zone, we simulated reactive transport in columns designed to represent local site conditions. The Pitzer ion interaction module was used, with principal geochemical processes considered in the simulation including quartz dissolution, precipitation of brucite, calcite, and portlandite, multi-component cation exchange, and aqueous complexation reactions. Good matches were observed between the simulated and measured column data at ambient temperature ( approximately 21 degrees C). Relatively good agreement was also obtained at high temperature ( approximately 70 degrees C). The decrease of pH at the plume front is examined through formation of secondary mineral phases and/or quartz dissolution. Substantial formation of secondary mineral phases resulting from multi-component cation exchange suggests that these phases are responsible for a decrease in pH within the plume front. In addition, a sensitivity analysis was conducted with respect to cation exchange capacity, selectivity coefficient, mineral assemblage, temperature, and ionic strength. This study could serve as a useful guide to subsequent experimental work, to thermodynamic models developed for the concentrated solutions at high ionic strength and to other types of waste plume studies.     

Reactive transport modelling of biogeochemical processes and carbon isotope geochemistry inside a landfill leachate plume

The biogeochemical processes governing leachate attenuation inside a landfill leachate plume (Banisveld, the Netherlands) were revealed and quantified using the 1D reactive transport model PHREEQC-2. Biodegradation of dissolved organic carbon (DOC) was simulated assuming first-order oxidation of two DOC fractions with different reactivity, and was coupled to reductive dissolution of iron oxide. The following secondary geochemical processes were required in the model to match observations: kinetic precipitation of calcite and siderite, cation exchange, proton buffering and degassing. Rate constants for DOC oxidation and carbonate mineral precipitation were determined, and other model parameters were optimized using the nonlinear optimization program PEST by means of matching hydrochemical observations closely (pH, DIC, DOC, Na, K, Ca, Mg, NH4, Fe(II), SO4, Cl, CH4, saturation index of calcite and siderite). The modelling demonstrated the relevance and impact of various secondary geochemical processes on leachate plume evolution. Concomitant precipitation of siderite masked the act of iron reduction. Cation exchange resulted in release of Fe(II) from the pristine anaerobic aquifer to the leachate. Degassing, triggered by elevated CO2 pressures caused by carbonate precipitation and proton buffering at the front of the plume, explained the observed downstream decrease in methane concentration. Simulation of the carbon isotope geochemistry independently supported the proposed reaction network.     

Large-scale modeling of reactive solute transport in fracture zones of granitic bedrocks

Final disposal of high-level radioactive waste in deep repositories located in fractured granite formations is being considered by several countries. The assessment of the safety of such repositories requires using numerical models of groundwater flow, solute transport and chemical processes. These models are being developed from data and knowledge gained from in situ experiments such as the Redox Zone Experiment carried out at the underground laboratory of Äspö in Sweden. This experiment aimed at evaluating the effects of the construction of the access tunnel on the hydrogeological and hydrochemical conditions of a fracture zone intersected by the tunnel. Most chemical species showed dilution trends except for bicarbonate and sulphate which unexpectedly increased with time. Molinero and Samper [Molinero, J. and Samper, J. Groundwater flow and solute transport in fracture zones: an improved model for a large-scale field experiment at Äspö (Sweden). J. Hydraul. Res., 42, Extra Issue, 157–172] presented a two-dimensional water flow and solute transport finite element model which reproduced measured drawdowns and dilution curves of conservative species. Here we extend their model by using a reactive transport which accounts for aqueous complexation, acid–base, redox processes, dissolution–precipitation of calcite, quartz, hematite and pyrite, and cation exchange between Na⁺ and Ca²⁺. The model provides field-scale estimates of cation exchange capacity of the fracture zone and redox potential of groundwater recharge. It serves also to identify the mineral phases controlling the solubility of iron. In addition, the model is useful to test the relevance of several geochemical processes. Model results rule out calcite dissolution as the process causing the increase in bicarbonate concentration and reject the following possible sources of sulphate: (1) pyrite dissolution, (2) leaching of alkaline sulphate-rich waters from a nearby rock landfill and (3) dissolution of iron monosulphides contained in Baltic seafloor sediments. Based on these results, microbially mediated processes are postulated as the most likely hypothesis to explain the measured increase of dissolved bicarbonates and sulphates after tunnel construction.     

Transport of Sr2+ and SrEDTA2- in partially-saturated and heterogeneous sediments

Strontium-90 has migrated deep into the unsaturated subsurface beneath leaking storage tanks in the Waste Management Areas (WMA) at the U.S. Department of Energy's (DOE) Hanford Reservation. Faster than expected transport of contaminants in the vadose zone is typically attributed to either physical hydrologic processes such as development of preferential flow pathways, or to geochemical processes such as the formation of stable, anionic complexes with organic chelates, e.g., ethylenediaminetetraacetic acid (EDTA). The goal of this paper is to determine whether hydrological processes in the Hanford sediments can influence the geochemistry of the system and hence control transport of Sr(2+) and SrEDTA(2-). The study used batch isotherms, saturated packed column experiments, and an unsaturated transport experiment in an undisturbed core. Isotherms and repacked column experiments suggested that the SrEDTA(2-) complex was unstable in the presence of Hanford sediments, resulting in dissociation and transport of Sr(2+) as a divalent cation. A decrease in sorption with increasing solid:solution ratio for Sr(2+) and SrEDTA(2-) suggested mineral dissolution resulted in competition for sorption sites and the formation of stable aqueous complexes. This was confirmed by detection of MgEDTA(2-), MnEDTA(2-), PbEDTA(2-), and unidentified Sr and Ca complexes. Displacement of Sr(2+) through a partially-saturated undisturbed core resulted in less retardation and more irreversible sorption than was observed in the saturated repacked columns, and model results suggested a significant reservoir (49%) of immobile water was present during transport through the heterogeneous layered sediments. The undisturbed core was subsequently disassembled along distinct bedding planes and subjected to sequential extractions. Strontium was unequally distributed between carbonates (49%), ion exchange sites (37%), and the oxide (14%) fraction. An inverse relationship between mass wetness and Sr suggested that sandy sediments of low water content constituted the immobile flow regime. Our results suggested that the sequestration of Sr(2+) in partially-saturated, heterogeneous sediments was most likely due to the formation of immobile water in drier regions having low hydraulic conductivities.     

Analysis of a deep well recharge experiment by calibrating a reactive transport model with field data

This paper describes the modeling of the hydrogeochemical effects of deep well recharge of oxic water into an anoxic pyrite-bearing aquifer. Kinetic expressions have been used for mineral dissolution-precipitation rates and organic matter oxidation. Hydrological and chemical parameters of the model were calibrated to field measurements. The results showed that oxidation of pyrite (FeS(2)) and, to a lesser extent, organic matter dominate the changes in quality of the recharged water during its passage through the aquifer. The recharge leads to the consumption of oxygen and nitrate and the formation of sulfate and ferrihydrite. Complexation reactions, cation exchange and precipitation and dissolution of calcite, siderite and rhodochrosite were also identified through the modeling. Despite problems of non-uniqueness of the calibrated parameters, the model was used successfully to depict the geochemical processes occurring in the aquifer. Non-uniqueness can be avoided by constraining the model as much as possible to measurements and/or data from literature, although they cannot be considered always as fixed values and should be considered as stochastic variables instead.     

Obtaining the porewater composition of a clay rock by modeling the in- and out-diffusion of anions and cations from an in-situ experiment

A borehole in the Callovo–Oxfordian clay rock in ANDRA's underground research facility was sampled during 1 year and chemically analyzed. Diffusion between porewater and the borehole solution resulted in concentration changes which were modeled with PHREEQC's multicomponent diffusion module. In the model, the clay rock's pore space is divided in free porewater (electrically neutral) and diffuse double layer water (devoid of anions). Diffusion is calculated separately for the two domains, and individually for all the solute species while a zero-charge flux is maintained. We explain how the finite difference formulas for radial diffusion can be translated into mixing factors for solutions. Operator splitting is used to calculate advective flow and chemical reactions such as ion exchange and calcite dissolution and precipitation. The ion exchange reaction is formulated in the form of surface complexation, which allows distributing charge over the fixed sites and the diffuse double layer. The charge distribution affects pH when calcite dissolves, and modeling of the experimental data shows that about 7% of the cation exchange capacity resides in the diffuse double layer. The model calculates the observed concentration changes very well and provides an estimate of the pristine porewater composition in the clay rock.     

Geochemical interactions between process-affected water from oil sands tailings ponds and North Alberta surficial sediments

In Northern Alberta, the placement of out-of-pit oil sands tailings ponds atop natural buried sand channels is becoming increasingly common. Preliminary modeling of such a site suggests that process-affected (PA) pond water will infiltrate through the underlying clay till aquitard, reaching the sand channel. However, the impact of seepage upon native sediments and groundwater resources is not known. The goal of this study is to investigate the role of adsorption and ion exchange reactions in the clay till and their effect on the attenuation or release of inorganic species. This was evaluated using batch sorption experiments (traditional and a recent modification using less disturbed sediment samples) and geochemical modeling with PHREEQC. The results show that clay till sediments have the capacity to mitigate the high concentrations of ingressing sodium (600 mg L(-1)), with linear sorption partitioning coefficients (K(d)) of 0.45 L kg(-1). Ion exchange theory was required to account for all other cation behaviour, precluding the calculation of such coefficients for other species. Qualitative evidence suggests that chloride will behave conservatively, with high concentrations remaining in solution (375 mg L(-1)). As a whole, system behaviour was found to be controlled by a combination of competitive ion exchange, dissolution and precipitation reactions. Observations, supported by PHREEQC simulations, suggest that the influx of PA water will induce the dissolution of pre-existing sulphate salts. Sodium present in the process-affected water will exchange with sediment-bound calcium and magnesium, increasing the divalent ions' pore fluid concentrations, and leading to the precipitation of a calcium-magnesium carbonate mineral phase. Thus, in similar tailings pond settings, particularly if the glacial till coverage is thin or altogether absent, it is reasonable to expect that high concentrations of sodium and chloride will remain in solution, while sulphate concentrations will exceed those of the ingressing plume (150 mg L(-1)).     

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.     

Manganese and trace-metal mobility under reducing conditions following in situ oxidation of TCE by KMnO4: a laboratory column experiment

The stability of Mn oxides, and the potential for mobilization of associated trace metals, were assessed by simulating the onset of microbially-mediated reducing conditions in a continuous-flow column experiment. The column had previously been used for an in situ chemical oxidation (ISCO) experiment in which trichloroethylene was reacted with permanganate in the presence of aqueous trace metals, which produced Mn oxyhydroxides (MnO(x)) that sequestered the trace metals and coated the column sand. The column influent solution represented the incursion of ambient groundwater containing dissolved organic carbon (DOC) into an ISCO treatment zone. The influx of DOC-containing groundwater initiated a series of cation-exchange, surface-complexation and reductive-dissolution reactions that controlled the release of aqueous metals from the system. Peak concentrations in the effluent occurred in the order Na, Mo, Cr, Zn, K, Mn, Fe, Pb, Mg, Ni, Cu and Ca. Manganese release from the column was controlled by a combination of cation exchange, reductive dissolution and precipitation of rhodochrosite. The trend in Fe concentrations was similar to that of Mn, and also resulted from a combination of reductive dissolution and cation exchange. Cation exchange and/or surface-complexation were the primary mechanisms controlling Cu, Ni, Mo and Pb release to solution, while Zn and Cr concentrations did not display coherent trends. Although metal release from the treatment zone was evident in the data, concentrations of trace metals remained below 0.05 mg L(-1) with the exception of Mo which reached concentrations on the order of 1 mg L(-1). The establishment of anaerobic conditions in ISCO-treated aquifers may result in a prolonged flux of aqueous Mn(II), but with the exception of MoO(4)(2-), it is unlikely that trace metals sequestered with MnO(x) during ISCO will be released to the groundwater in elevated concentrations.     

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.     

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.     

Effects of water content on reactive transport of 85Sr in Chernobyl sand columns

It is known that under unsaturated conditions, the transport of solutes can deviate from ideal advective-dispersive behaviour even for macroscopically homogeneous porous materials. Causes may include physical non-equilibrium, sorption kinetics, non-linear sorption, and the irregular distribution of sorption sites. We have performed laboratory experiments designed to identify the processes responsible for the non-ideality of radioactive Sr transport observed under unsaturated flow conditions in an Aeolian sandy deposit from the Chernobyl exclusion zone. Miscible displacement experiments were carried out at various water contents and corresponding flow rates in a laboratory model system. Results of our experiments have shown that breakthrough curves of a conservative tracer exhibit a higher degree of asymmetry when the water content decreases than at saturated water content and same Darcy velocity. It is possible that velocity variations caused by heterogeneities at the macroscopic scale are responsible for this situation. Another explanation is that molecular diffusion drives the solute mass transfer between mobile and immobile water regions, but the surface of contact between these water regions is small. At very low concentrations, representative of a radioactive Sr contamination of the pore water, sorption and physical disequilibrium dominate the radioactive Sr transport under unsaturated flow conditions. A sorption reaction is described by a cation exchange mechanism calibrated under fully saturated conditions. The sorption capacity, as well as the exchange coefficients are not affected by desaturation. The number of accessible exchange sites was calculated on the basis that the solid remained in contact with water and that the fraction of solid phase in contact with mobile water is numerically equal to the proportion of mobile water to total water content. That means that for this type of sandy soil, the nature of mineral phases is the same in advective and non-advective domains. So sorption reaction parameters can be estimated from more easily conducted saturated experiments, but hydrodynamic behaviour must be characterized by conservative tracer experiments under unsaturated flow conditions.     

Bioremediation of strontium (Sr) contaminated aquifer quartz sand based on carbonate precipitation induced by Sr resistant Halomonas sp

Contamination of aquifers or sediments by radioactive strontium (⁹⁰Sr) is a significant environmental problem. In the present study, microbially induced calcite precipitation (MICP) was evaluated for its potential to remediate strontium from aquifer quartz sand. A Sr resistant urease producing Halomonas sp. was characterized for its potential role in bioremediation. The bacterial strain removed 80% of Sr from soluble-exchangeable fraction of aquifer quartz sand. X-ray diffraction detected calcite, vaterite and aragonite along with calcite-strontianite (SrCO₃) solid solution in bioremediated sample with indications that Sr was incorporated into the calcite. Scanning electron micrography coupled with energy-dispersive X-ray further confirmed MICP process in remediation. The study showed that MICP sequesters soluble strontium as biominerals and could play an important role in strontium bioremediation from both ecological and greener point of view.     

Numerical simulation of reactive processes in an experiment with partially saturated bentonite

Bentonites are preferred materials for use as engineered barriers for high-level nuclear waste repositories. Simulation of geochemical processes in bentonite is therefore important for long-term safety assessment of those repositories. In this work, the porewater chemistry of a bentonite sample subject to simultaneous heating and hydration, as studied by Cuevas et al. [Cuevas, J., Villar, M., Fernández, A., Gómez, P., Martín, P., 1997. Porewaters extracted from compacted bentonite subjected to simultaneous heating and hydration. Applied Geochemistry 12, 473-481.], was assessed with a non-isothermal reactive transport model by coupling the geochemical software PHREEQC2 with the object-oriented FEM simulator GeoSys/RockFlow. Reactive transport modelling includes heat transport, two-phase flow, multicomponent transport and geochemical reactions in the liquid phase, i.e. ion exchange, mineral dissolution/precipitation and equilibrium reactions. Simulations show that the easily soluble minerals in bentonite determine the porewater chemistry. Temperature affects both two-phase flow and geochemical reactions. Porosity change due to dissolution/precipitation is low during the experiment. However, changes of the effective porosity caused by bentonite swelling can be very large. The simulated results agree well with the experimental data.     

Processes controlling metal transport and retention as metal-contaminated groundwaters efflux through estuarine sediments

Factors affecting the transport and retention of Cd, Cr, Cu, Ni, Pb and Zn in acidic groundwaters as they pass through estuarine sediments were investigated using column experiments. Acidic groundwaters caused the rapid dissolution of iron sulfide (AVS) and other iron and manganese phases from sediments that are important for metal binding and buffering. Metal breakthrough to overlying water occurred in the order of Ni>Zn>Cd>Cu>Cr/Pb. Metal transport increased as the sediment permeability increased, reflecting the low resistance to flow caused by larger sand-sized particles and the decreased abundance of metal adsorption sites on these materials. Metal mobility increased as the groundwater pH decreased, as flow rate or metal concentrations increased, and as the exposure duration increased. Groundwater Cr and Pb were promptly attenuated by the sediments, the mobility of Cu was low and decreased rapidly as sediment pH increased above 4.5, while Cd, Ni and Zn were the most easily transported to the surface sediments and released to the overlying waters. For groundwaters of pH 3, metal migration velocities through sandy sediments were generally 0.5-2% (Cr, Pb), 1-6% (Cu) and 4-13% (Cd, Ni, Zn) of the total groundwater velocity (9-700 m/yr). The oxidative precipitation of Fe(II) and Mn(II) in the groundwaters did not affect metal mobility through the sediments. The results indicated that the efflux of acidic and metal-contaminated groundwater through estuarine sediments would affect organisms resident in sandy sediments more greatly than organisms resident in fine-grained, silty, sediments.     

Colloid‐facilitated transport of metolachlor through intact soil columns

Abstract

This study evaluated the role of water dispersible colloids with diverse physicochemical and mineralogical characteristics in facilitating the transport of metolachlor through macropores of intact soil columns. The soil columns represented upper solum horizons of an Alfisol in the Bluegrass region of Kentucky. Three different colloid suspensions tagged with metolachlor [2‐chloro‐N‐(2‐ethyl‐6‐methylphenyl)‐N‐(2‐methoxy‐l‐methylethyl)acetamide] were introduced at a constant flux into undisturbed soil columns. The eluents were collected and analyzed periodically for colloid and metolachlor concentrations. Colloid recovery in the eluents ranged from 54 to 90 %. The presence of colloids enhanced the transport of metolachlor by 22 to 70 % depending on the colloid type and mobility. Colloids with higher pH, organic carbon, cation exchange capacity (CEC), total exchangeable bases (TEB), surface area (SA), and electrophoretic mobility (EM), showed better mobility, greater affinity for interaction with the herbicide and, thus, greater potential to co‐transport metolachlor. In contrast, increased level of kaolinite, Fe, and Al inhibited metolachlor adsorption and transport. In spite of the increased transportability of metolachlor by the presence of soil colloids, the colloid bound herbicide portion accounted for a very small part of the observed increase. This suggests that surface site exclusion mechanisms and preferential sorption induced by the presence of colloids are more important than ion exchange phenomena in promoting herbicide mobility in subsurface environments.     

A reaction-transport model for calcite precipitation and evaluation of infiltration fluxes in unsaturated fractured rock

The percolation flux in the unsaturated zone (UZ) is an important parameter addressed in site characterization and flow and transport modeling of the potential nuclear-waste repository at Yucca Mountain, NV, USA. The US Geological Survey (USGS) has documented hydrogenic calcite abundances in fractures and lithophysal cavities at Yucca Mountain to provide constraints on percolation fluxes in the UZ. The purpose of this study was to investigate the relationship between percolation flux and measured calcite abundances using reactive transport modeling. Our model considers the following essential factors affecting calcite precipitation: (1) infiltration, (2) the ambient geothermal gradient, (3) gaseous CO(2) diffusive transport and partitioning in liquid and gas phases, (4) fracture-matrix interaction for water flow and chemical constituents, and (5) water-rock interaction. Over a bounding range of 2-20 mm/year infiltration rate, the simulated calcite distributions capture the trend in calcite abundances measured in a deep borehole (WT-24) by the USGS. The calcite is found predominantly in fractures in the welded tuffs, which is also captured by the model simulations. Simulations showed that from about 2 to 6 mm/year, the amount of calcite precipitated in the welded Topopah Spring tuff is sensitive to the infiltration rate. This dependence decreases at higher infiltration rates owing to a modification of the geothermal gradient from the increased percolation flux. The model also confirms the conceptual model for higher percolation fluxes in the fractures compared to the matrix in the welded units, and the significant contribution of Ca from water-rock interaction. This study indicates that reactive transport modeling of calcite deposition can yield important constraints on the unsaturated zone infiltration-percolation flux and provide useful insight into processes such as fracture-matrix interaction as well as conditions and parameters controlling calcite deposition.     

Dissolution and sorption of hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) and 2,4,6-trinitrotoluene (TNT) residues from detonated mineral surfaces

Composition B (Comp B) is a commonly used military formulation composed of the toxic explosive compounds 2,4,6-trinitrotoluene (TNT), and hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX). Numerous studies of the temporal fate of explosive compounds in soils, surface water and laboratory batch reactors have been conducted. However, most of these investigations relied on the application of explosive compounds to the media via aqueous addition and thus these studies do not provide information on the real world loading of explosive residues during detonation events. To address this we investigated the dissolution and sorption of TNT and RDX from Comp B residues loaded to pure mineral phases through controlled detonation. Mineral phases included nontronite, vermiculite, biotite and Ottawa sand (quartz with minor calcite). High Performance Liquid Chromatography and Attenuated Total Reflectance Fourier Transform Infrared spectroscopy were used to investigate the dissolution and sorption of TNT and RDX residues loaded onto the mineral surfaces. Detonation resulted in heterogeneous loading of TNT and RDX onto the mineral surfaces. Explosive compound residues dissolved rapidly (within 9 h) in all samples but maximum concentrations for TNT and RDX were not consistent over time due to precipitation from solution, sorption onto mineral surfaces, and/or chemical reactions between explosive compounds and mineral surfaces. We provide a conceptual model of the physical and chemical processes governing the fate of explosive compound residues in soil minerals controlled by sorption-desorption processes.     

Thioarsenate formation upon dissolution of orpiment and arsenopyrite

Thioarsenates were previously determined as dominant species in geothermal and mineral waters with excess sulfide. Here, we used batch leaching experiments to determine their formation upon weathering or industrial leaching of the arsenic-sulfide minerals orpiment (As₂S₃) and arsenopyrite (FeAsS) under different pH and oxygen conditions. Under acidic conditions, as expected based on their known kinetic instability at low pH, no thioarsenates formed in either of the two mineral systems. Under neutral to alkaline conditions, orpiment dissolution yielded mono-, di- and trithioarsenate which accounted for up to 43-55% of total arsenic. Thioarsenate formation upon arsenopyrite dissolution was low at neutral (4%) but significant at alkaline pH, especially under suboxic to sulfidic conditions (20-43%, mainly as monothioarsenate). In contrast to orpiment, we postulate that recombination of arsenite and sulfide in solution is of minor importance for monothioarsenate formation during alkaline arsenopyrite dissolution. We propose instead that hydroxyl physisorption lead to formation of As-OH-S surface complexes by transposition of hydroxyl anions to arsenic or iron sites. Concurrently formed ironhydroxides could provide re-sorption sites for the freshly released monothioarsenate. However, sorption experiments with goethite showed slower sorption kinetics of monothioarsenate compared to arsenite, but comparable with arsenate. The discovery that thioarsenates are released by natural weathering and industrial leaching processes and that, once they are released, have a higher mobility than the commonly-investigated species arsenite and arsenate requires future studies to consider them when assessing arsenic release in sulfidic natural or mining-impacted environments.     

Thioarsenate formation upon dissolution of orpiment and arsenopyrite

Thioarsenates were previously determined as dominant species in geothermal and mineral waters with excess sulfide. Here, we used batch leaching experiments to determine their formation upon weathering or industrial leaching of the arsenic-sulfide minerals orpiment (As₂S₃) and arsenopyrite (FeAsS) under different pH and oxygen conditions. Under acidic conditions, as expected based on their known kinetic instability at low pH, no thioarsenates formed in either of the two mineral systems. Under neutral to alkaline conditions, orpiment dissolution yielded mono-, di- and trithioarsenate which accounted for up to 43–55% of total arsenic. Thioarsenate formation upon arsenopyrite dissolution was low at neutral (4%) but significant at alkaline pH, especially under suboxic to sulfidic conditions (20–43%, mainly as monothioarsenate). In contrast to orpiment, we postulate that recombination of arsenite and sulfide in solution is of minor importance for monothioarsenate formation during alkaline arsenopyrite dissolution. We propose instead that hydroxyl physisorption lead to formation of As-OH-S surface complexes by transposition of hydroxyl anions to arsenic or iron sites. Concurrently formed ironhydroxides could provide re-sorption sites for the freshly released monothioarsenate. However, sorption experiments with goethite showed slower sorption kinetics of monothioarsenate compared to arsenite, but comparable with arsenate. The discovery that thioarsenates are released by natural weathering and industrial leaching processes and that, once they are released, have a higher mobility than the commonly-investigated species arsenite and arsenate requires future studies to consider them when assessing arsenic release in sulfidic natural or mining-impacted environments.