Intercomparison of reactive transport models applied to UO2 oxidative dissolution and uranium migration

Oxidative dissolution of uranium dioxide (UO₂) and the subsequent migration of uranium in a subsurface environment and an underground waste disposal have been simulated with reactive transport models. In these systems, hydrogeological and chemical processes are closely entangled and their interdependency has been analyzed in detail, notably with respect to redox reactions, kinetics of mineralogical evolution and hydrodynamic migration of species of interest.Different codes, where among CASTEM, CHEMTRAP and HYTEC, have been used as an intercomparison and verification exercise. Although the agreement between codes is satisfactory, it is shown that the discretization method of the transport equation (i.e. finite elements (FE) versus mixed-hybrid FE and finite differences) and the sequential coupling scheme may lead to systematic discrepancies.     

Intercomparison between TRIO-EF and IMPACT codes with reference to experimental strontium migration data

This work presents an intercomparison exercise between two geochemical migration codes, TRIO-EF (an object-oriented finite element code) and IMPACT (a chemical engineering code using mixing cells in series). The predictions of the two codes are compared with the reference experimental results obtained in a previous study of strontium transport in soil columns. This simulated geochemical system is well documented and includes ion exchange and dissolution-precipitation reactions. The solution transport is simulated by a one-dimensional advection-dispersion model. The predictions of TRIO-EF and IMPACT are both in good agreement with the experimental results. However, slight differences can be observed between the two codes, especially when concentration discontinuities are involved, such as precipitation fronts or changes in boundary conditions. These discrepancies between the two codes can mainly be attributed to the different discretisation approaches.     

The migration of uranium into peat-rich soils at Broubster, Caithness, Scotland, U.K.

Uranium is being actively transported from uraniferous sedimentary rocks into a peat bog at the Broubster natural analogue site in Caithness, Scotland. Massive calcareous sandstone within the Caithness Flags sequence is the main source of uranium which resides primarily within diagenetic apatite and dispersed U---Si---Ti phases. Supergene weathering processes have decalcified the sandstone and are effective in mobilising uranium by groundwater leaching, primarily along a fault zone. Uranium transport in solution by means of groundwater and surface flows is effectively terminated by retardation within 4-kyr-old peat deposits laid down on boulder clay. This process of secondary fixation has resulted in a young uranium anomaly in the peat which comprises in excess of 0.1 wt% U. The site has been investigated comprehensively to define the geometry of the anomaly together with the hydrogeology, hydrochemistry, petrology, mineralogy and the nature of the peat sink-term. The main physical and geochemical properties of the system, including the uranium decay series radionuclide distributions in water and solid samples, are documented in this paper. From these data, the processes governing the distribution of uranium have been quantified using a three-dimensional groundwater flow package and an equilibrium speciation model incorporating a recently developed electrostatic surface complexation model to account for cation-organic interactions. The results described form part of a coordinated project on natural radionuclide migration undertaken to improve confidence in predictive methods used for radiological assessment.     

Kinetic modeling of microbially-driven redox chemistry of radionuclides in subsurface environments: coupling transport, microbial metabolism and geochemistry

Microbial reactions play an important role in regulating pore water chemistry as well as secondary mineral distribution in many subsurface systems and, therefore, may directly impact radionuclide migration in those systems. This paper presents a general modeling approach to couple microbial metabolism, redox chemistry, and radionuclide transport in a subsurface environment. To account for the likely achievement of quasi-steady state biomass accumulations in subsurface environments, a modification to the traditional microbial growth kinetic equation is proposed. The conditions for using biogeochemical models with or without an explicit representation of biomass growth are clarified. Based on the general approach proposed in this paper, the couplings of uranium reactions with biogeochemical processes are incorporated into computer code BIORXNTRN Version 2.0. The code is then used to simulate a subsurface contaminant migration scenario, in which a water flow containing both uranium and a complexing organic ligand is recharged into an oxic carbonate aquifer. The model simulation shows that Mn and Fe oxyhydroxides may vary significantly along a flow path. The simulation also shows that uranium(VI) can be reduced and therefore immobilized in the anoxic zone created by microbial degradation.     

Database uncertainty as a limiting factor in reactive transport prognosis

The effect of uncertainties in thermodynamic databases on prediction performances of reactive transport modeling of uranium (VI) is investigated with a Monte Carlo approach using the transport code TReaC. TReaC couples the transport model to the speciation code PHREEQC by a particle tracking method. A speciation example is given to illustrate the effect of uncertainty in thermodynamic data on the predicted solution composition. The transport calculations consequently show the prediction uncertainty resulting from uncertainty in thermodynamic data. A conceptually simple scenario of elution of uranium from a sand column is used as an illustrating example. Two different cases are investigated: a carbonate-enriched drinking water and an acid mine water associated with uranium mine remediation problems. Due to the uncertainty in the relative amount of positively charged and neutral solution species, the uncertainty in the thermodynamic data also infers uncertainty in the retardation behavior. The carbonated water system shows the largest uncertainties in speciation calculation. Therefore, the model predictions of total uranium solubility have a broad range. The effect of data uncertainty in transport prediction is further illustrated by a prediction of the time when eluted uranium from the column exceeds a threshold value. All of these Monte Carlo transport calculations consume large amounts of computing time.     

Analytical solutions for multiple species reactive transport in multiple dimensions

Many numerical computer codes used to simulate multi-species reactive transport and biodegradation have been developed in recent years. Such numerical codes must be validated by comparison of the numerical solutions with an analytical solution. In this paper, a method for deriving analytical solutions of the partial differential equations describing multiple species multi-dimensional transport with first-order sequential reactions is presented. Although others have developed specific solutions of multi-species transport equations, here a more general analytical approach, capable of describing any number of reactive species in multiple dimensions is derived. A substitution method is used to transform the multi-species reactive transport problem to one that can be solved using previously published single-species solutions for various initial and boundary conditions. One- and three-dimensional examples are presented to illustrate the steps involved in extending single-species solutions to a four-species system with sequential first-order reactions.     

Sorption and bioreduction of hexavalent uranium at a military facility by the Chesapeake Bay

Directly adjacent to the Chesapeake Bay lies the Aberdeen Proving Ground, a U.S. Army facility where testing of armor-piercing ammunitions has resulted in the deposition of >70,000 kg of depleted uranium (DU) to local soils and sediments. Results of previous environmental monitoring suggested limited mobilization in the impact area and no transport of DU into the nation's largest estuary. To determine if physical and biological reactions constitute mechanisms involved in limiting contaminant transport, the sorption and biotransformation behavior of the radionuclide was studied using geochemical modeling and laboratory microcosms (500 ppb U(VI) initially). An immediate decline in dissolved U(VI) concentrations was observed under both sterile and non-sterile conditions due to rapid association of U(VI) with natural organic matter in the sediment. Reduction of U(VI) to U(IV) occurred only in non-sterile microcosms. In the non-sterile samples, intrinsic bioreduction of uranium involved bacteria of the order Clostridiales and was only moderately enhanced by the addition of acetate (41% vs. 56% in 121 days). Overall, this study demonstrates that the migration of depleted uranium from the APG site into the Chesapeake Bay may be limited by a combination of processes that include rapid sorption of U(VI) species to natural organic matter, followed by slow, intrinsic bioreduction to U(IV).     

Uranium transport around the reactor zone at Bangombé and Okélobondo (Oklo): examples of hydrogeological and geochemical model integration and data evaluation

The sites at Bangombé and Okélobondo (Oklo) in Gabon provide a unique opportunity to study the behaviour of products from natural nuclear reactions in the vicinity of reactor zones which were active around two billion years ago. The Commission of the European Communities initiated the Oklo Natural Analogue Programme. One of the principal aims was to study indications of present time migration of elements from the reactor zones under ambient conditions. The hydrogeological and hydrochemical data from the Oklo sites were modelled in order to better understand the geochemical behaviour of radionuclides in the natural system, by using independent models and by comparing the modelling outcome. Two modelling approaches were used: M3 code (hydrochemical mixing and mass balance model), developed by the Swedish Nuclear Fuel and Waste Management Company (SKB) and HYTEC (reactive transport model) developed by Ecole des Mines de Paris.Two different reactor zones were studied: Bangombé, a shallow site, the reactor being at 11 m depth, and OK84 at Okélobondo, situated at about 450 m depth, more comparable with a real repository location. This allowed the validation of modelling tools in two different sedimentary environments: one shallow, with a more homogeneous layering situated in an area of meteoric alteration, and the other offering the opportunity to study radionuclide migration from the reaction zone over a distance of 450 m through very heterogeneous sedimentary layers.The modeling results indicate that the chemical reactions retarding radionuclide transport are very different at the two sites. At Bangombé, the decomposition of organic material consumes oxygen and at Okélobondo the oxygen is consumed by inorganic reactions resulting, in both cases, in uranium retardation. Both modelling approaches (statistic with M3 code and deterministic with HYTEC code) could describe this situation.The goal of this exercise is to test codes which can help to describe and understand the processes taking place at the sites, validate the models with in situ data, and thus build confidence in the tools used for future site characterization. Ultimately, this allows identifying and selecting processes and parameters that can be used as input into repository performance assessment calculations and modelling exercises.     

Uranium transport around the reactor zone at Bangombé and Okélobondo (Oklo): examples of hydrogeological and geochemical model integration and data evaluation

The sites at Bangombé and Okélobondo (Oklo) in Gabon provide a unique opportunity to study the behaviour of products from natural nuclear reactions in the vicinity of reactor zones which were active around two billion years ago. The Commission of the European Communities initiated the Oklo Natural Analogue Programme. One of the principal aims was to study indications of present time migration of elements from the reactor zones under ambient conditions. The hydrogeological and hydrochemical data from the Oklo sites were modelled in order to better understand the geochemical behaviour of radionuclides in the natural system, by using independent models and by comparing the modelling outcome. Two modelling approaches were used: M3 code (hydrochemical mixing and mass balance model), developed by the Swedish Nuclear Fuel and Waste Management Company (SKB) and HYTEC (reactive transport model) developed by Ecole des Mines de Paris. Two different reactor zones were studied: Bangombé, a shallow site, the reactor being at 11 m depth, and OK84 at Okélobondo, situated at about 450 m depth, more comparable with a real repository location. This allowed the validation of modelling tools in two different sedimentary environments: one shallow, with a more homogeneous layering situated in an area of meteoric alteration, and the other offering the opportunity to study radionuclide migration from the reaction zone over a distance of 450 m through very heterogeneous sedimentary layers. The modeling results indicate that the chemical reactions retarding radionuclide transport are very different at the two sites. At Bangombé, the decomposition of organic material consumes oxygen and at Okélobondo the oxygen is consumed by inorganic reactions resulting, in both cases, in uranium retardation. Both modelling approaches (statistic with M3 code and deterministic with HYTEC code) could describe this situation. The goal of this exercise is to test codes which can help to describe and understand the processes taking place at the sites, validate the models with in situ data, and thus build confidence in the tools used for future site characterization. Ultimately, this allows identifying and selecting processes and parameters that can be used as input into repository performance assessment calculations and modelling exercises.     

Bacterial interactions with uranium: An environmental perspective

The presence of actinides in radioactive wastes is of major concern because of their potential for migration from the waste repositories and long-term contamination of the environment. Studies have been and are being made on inorganic processes affecting the migration of radionuclides from these repositories to the environment but it is becoming increasingly evident that microbial processes are of importance as well. Bacteria interact with uranium through different mechanisms including, biosorption at the cell surface, intracellular accumulation, precipitation, and redox transformations (oxidation/reduction). The present study is intended to give a brief overview of the key processes responsible for the interaction of actinides e.g. uranium with bacterial strains isolated from different extreme environments relevant to radioactive repositories. Fundamental understanding of the interaction of these bacteria with U will be useful for developing appropriate radioactive waste treatments, remediation and long-term management strategies as well as for predicting the microbial impacts on the performance of the radioactive waste repositories.     

Stochastic-convective transport with nonlinear reactions and mixing: finite streamtube ensemble formulation for multicomponent reaction systems with intra-streamtube dispersion

An effective streamtube ensemble method is developed to upscale convective-dispersive transport with multicomponent nonlinear reactions in steady nonuniform flow. The transport is cast in terms of a finite ensemble of independent discrete streamtubes that approximate convective transport along macroscopically averaged pathlines and dispersive transport longitudinally as microscopic mixing within streamtubes. The representation of fate and transport via a finite ensemble of effective linear streamtubes, allows the treatment of arbitrarily complex reaction systems involving both homogeneous and heterogeneous reactions, and longitudinal dispersive/diffusive mixing within streamtubes. This allows the use of reactive-transport codes designed to solve such problems in an Eulerian framework, as opposed to reliance on closed-form (convolutional or canonical) expressions for reactive transport in exclusively convective streamtubes. The approach requires both reactive-transport solutions for a representative ensemble of one-dimensional convective-dispersive-reactive streamtubes and the distribution of flux over the streamtube ensemble variants, and it does not allow for lateral mixing between streamtubes. Here, the only ensemble variant is travel time. The discussion details the way that the conventional Eulerian fate and transport model is converted first into an ensemble of transports along three-dimensional streamtubes of unknown geometry, and then to approximate one-dimensional streamtubes that are designed to honor the important global properties of the transport. Conditions under which such an 'equivalent' ensemble of one-dimensional streamtubes are described. The breakthrough curve of a nonreactive tracer in the ensemble is expressed as a combined Volterra-Fredholm integral equation, which serves as the basis for estimation of the distribution of flux over the variant of the ensemble, travel time. Transient convective speed and the effects of errors in flux distributions are described, and the method is applied to a demonstration problem involving nonlinear multicomponent reaction kinetics and strongly nonuniform flow.     

Reactive transport of uranium(VI) and phosphate in a goethite-coated sand column: an experimental study

The migration of uranium(VI) in subsurface environments is strongly influenced by its adsorption/desorption reactions at the solid/solution interface. Phosphate is often present in subsurface systems and was shown to significantly affect U(VI) adsorption in previous batch experiments. In this study, column experiments were conducted to investigate the effects of phosphate on U(VI) adsorption and transport under flow conditions. The adsorption of U(VI) and phosphate was very low on pure quartz sand with negligible effects on U(VI) and phosphate transport. However, U(VI) and phosphate transport was retarded in a column packed with goethite-coated sand. The presence of phosphate, either as a co-solute with U(VI) or pre-adsorbed, greatly increased U(VI) adsorption and retardation. U(VI) and phosphate adsorption in our column experiments were rate-limited, and the adsorption of U(VI) and phosphate was not reversible, with kinetic limitations more pronounced for desorption than for adsorption. This study demonstrated the importance of phosphate in controlling U(VI) mobility in subsurface environments and helped illustrate some phenomena potentially applicable to U(VI) adsorption and transport in natural systems, especially where U(VI) adsorption is rate-limited.     

Simulating bioremediation of uranium-contaminated aquifers; uncertainty assessment of model parameters

Bioremediation of trace metals and radionuclides in groundwater may require the manipulation of redox conditions via the injection of a carbon source. For example, after nitrate has been reduced, soluble U(VI) can be reduced simultaneously with other electron acceptors such as Fe(III) or sulfate to U(IV), which may precipitate as a solid (uraninite).To simulate the numerous biogeochemical processes that will occur during the bioremediation of trace-metal-contaminated aquifers, a time-dependent one-dimensional reactive transport model has been developed. The model consists of a set of coupled mass balance equations, accounting for advection, hydrodynamic dispersion, and a kinetic formulation of the biological or chemical transformations affecting an organic substrate, electron acceptors, corresponding reduced species, and trace metal contaminants of interest, uranium in this study. This set of equations is solved numerically, using a finite difference approximation. The redox conditions of the domain are characterized by estimating the pE, based on the concentration of the dominant terminal electron acceptor and its corresponding reduced species. This pE and the concentrations of relevant species are then used by a modified version of MINTEQA2, which calculates the speciation/sorption and precipitation/dissolution of the species of interest under equilibrium conditions. Kinetics of precipitation/dissolution processes are described as being proportional to the difference between the actual and calculated equilibrium concentration. A global uncertainty assessment, determined by Random Sampling High Dimensional Model Representation (RS-HDMR), was performed to attain a phenomenological understanding of the origins of output variability and to suggest input parameter refinements as well as to provide guidance for field experiments to improve the quality of the model predictions. By decomposing the model output variance into its different input contributions, RS-HDMR can identify the model inputs with the most influence on various model outputs, as well as their behavior pattern on the model output. Simulations are performed to illustrate the effect of biostimulation on the fate of uranium in a saturated aquifer, and to identify the key processes that need to be characterized with the highest accuracy prior to designing a uranium bioremediation scheme.     

Numerical validation of a Eulerian hydrochemical code using a 1D multisolute mass transport system involving heterogeneous kinetically controlled reactions

It is demonstrated that at steady state, the 1D thermo-kinetic hydrochemical Eulerian mass balance equations in pure advective mode are indeed identical to the governing mass balance equations of a single reaction path (or geochemical) code in open system mode. Thus, both calculated reaction paths should be theoretically identical whatever the chemical complexity of the water-rock system (i.e., multicomponent, multireaction zones kinetically and equilibrium-controlled). We propose to use this property to numerically test the thermo-kinetic hydrochemical Eulerian codes and we employ it to verify the algorithm of the 1D finite difference code KIRMAT. Compared to the other methods to perform such numerical tests (i.e., comparisons with analytical, semi-analytical solutions, between two Eulerian hydrochemical codes), the advantage of this new method is the absence of constraints on the chemical complexity of the modelled water-rock systems. Moreover, the same thermo-kinetic databases and geochemical functions can be easily and mechanically used in both calculations, when the numerical reference comes from the Eulerian code with no transport terms (u and D = 0) and modify to be consistent with the definition of the open system mode in geochemical modelling. The ability of KIRMAT to treat multicomponent pure advective transport, subjected to several kinetically equilibrium-controlled dissolution and precipitation reactions, and to track their boundaries has been successfully verified with the property of interest. The required numerical validation of the reference calculations is bypassed in developing the Eulerian code from an already checked single reaction path code. A forward time-upstream weighting scheme (a mixing cell scheme) is used in this study. An appropriate choice of grid spacing allows to calculate within the grid size uncertainty the correct mineral reaction zone boundaries, despite the presence of numerical dispersion. Its correction enables us to improve the convergence and to extend the numerical test to mixed advective-dispersive mass transport. However, the skewness factor involves numerical oscillations that prevent to compute different grid spacing. The use of a different chemically controlled time step constraint in both calculations induces some inconsistencies into the validation tests. This numerical validation method may be applied as well as to check a thermo-kinetic hydrochemical finite element based code, from a 1D heterogeneous systems, and 2D-3D systems provided that they are designed so as to be 1D equivalent. A one-step algorithm and the use of a numerical reference coming from the Eulerian code to be tested ensure the potential success (accuracy) of the numerical validation method.     

Variably saturated flow and multicomponent biogeochemical reactive transport modeling of a uranium bioremediation field experiment

Three-dimensional, coupled variably saturated flow and biogeochemical reactive transport modeling of a 2008 in situ uranium bioremediation field experiment is used to better understand the interplay of transport and biogeochemical reactions controlling uranium behavior under pulsed acetate amendment, seasonal water table variation, spatially variable physical (hydraulic conductivity, porosity) and geochemical (reactive surface area) material properties. While the simulation of the 2008 Big Rusty acetate biostimulation field experiment in Rifle, Colorado was generally consistent with behaviors identified in previous field experiments at the Rifle IFRC site, the additional process and property detail provided several new insights. A principal conclusion from this work is that uranium bioreduction is most effective when acetate, in excess of the sulfate-reducing bacteria demand, is available to the metal-reducing bacteria. The inclusion of an initially small population of slow growing sulfate-reducing bacteria identified in proteomic analyses led to an additional source of Fe(II) from the dissolution of Fe(III) minerals promoted by biogenic sulfide. The falling water table during the experiment significantly reduced the saturated thickness of the aquifer and resulted in reactants and products, as well as unmitigated uranium, in the newly unsaturated vadose zone. High permeability sandy gravel structures resulted in locally high flow rates in the vicinity of injection wells that increased acetate dilution. In downgradient locations, these structures created preferential flow paths for acetate delivery that enhanced local zones of TEAP reactivity and subsidiary reactions. Conversely, smaller transport rates associated with the lower permeability lithofacies (e.g., fine) and vadose zone were shown to limit acetate access and reaction. Once accessed by acetate, however, these same zones limited subsequent acetate dilution and provided longer residence times that resulted in higher concentrations of TEAP reaction products when terminal electron donors and acceptors were not limiting. Finally, facies-based porosity and reactive surface area variations were shown to affect aqueous uranium concentration distributions with localized effects of the fine lithofacies having the largest impact on U(VI) surface complexation. The ability to model the comprehensive biogeochemical reaction network, and spatially and temporally variable processes, properties, and conditions controlling uranium behavior during engineered bioremediation in the naturally complex Rifle IFRC subsurface system required a subsurface simulator that could use the large memory and computational performance of a massively parallel computer. In this case, the eSTOMP simulator, operating on 128 processor cores for 12h, was used to simulate the 110-day field experiment and 50 days of post-biostimulation behavior.     

Humic colloid-borne migration of uranium in sand columns

Column experiments were carried out to investigate the influence of humic colloids on subsurface uranium migration. The columns were packed with well-characterized aeolian quartz sand and equilibrated with groundwater rich in humic colloids (dissolved organic carbon (DOC): 30 mg dm(-3)). U migration was studied under an Ar/1% CO2 gas atmosphere as a function of the migration time, which was controlled by the flow velocity or the column length. In addition, the contact time of U with groundwater prior to introduction into a column was varied. U(VI) was found to be the dominant oxidation state in the spiked groundwater. The breakthrough curves indicate that U was transported as a humic colloid-borne species with a velocity up to 5% faster than the mean groundwater flow. The fraction of humic colloid-borne species increases with increasing prior contact time and also with decreasing migration time. The migration behavior was attributed to a kinetically controlled association/dissociation of U onto and from humic colloids and also a subsequent sorption of U onto the sediment surface. The column experiments provide an insight into humic colloid-mediated U migration in subsurface aquifers.     

Modelling of nonreactive tracer dipole tests in a shear zone at the Grimsel test site

The investigation of the migration of a high pH plume in a fractured shear zone is foreseen by a long-term experiment at the Grimsel rock laboratory. In order to characterise the initial conditions for the long-term experiment and to evaluate an optimal hydraulic in situ set-up, several dipole experiments with nonreacting tracers have been performed. The dipole experiments differ in geometry, pumping rates and orientation to the background water flow. Several single and double-porosity models have been applied to fit the results of these dipole tracer tests in order to extract values for some transport parameters and discriminate for certain transport processes. A two-dimensional porous medium approach was successfully used to fit tracer breakthrough curves measured for a dipole experiment. A model based on a one-dimensional dual porous medium approach was also successful, although the applied hydraulic dipole, with similar injection and extraction rates, suggests the existence of an extended two-dimensional flow field. For the two-dimensional porous medium approach, tracer breakthrough could only be fitted with a complex flow field geometry within the heterogeneous fractured shear zone. The heterogeneity was generated by heterogeneous porosity and hydraulic permeability distributions. Predictions for further dipole geometries and a sorbing tracer have been calculated by means of both models using the flow and transport parameters deduced from fits for a single dipole experiment. This allows for comparison with the measured breakthrough of sorbing tracers. The foreseen experiment with sorbing (radionuclide) tracers will help decide on the appropriate approach that should be used to describe such dipole experiments in this shear zone. Additionally, the migration and spreading of a solution with high pH has been calculated taking into account mineral dissolution and precipitation in a two-dimensional porous medium approach in order to estimate the amount and character of the mineral reactions induced by the interaction between the high pH solution and the rock.     

Mineralogical compositions of aquifer matrix as necessary initial conditions in reactive contaminant transport models

Mineralogical compositions and their spatial distributions are important initial conditions for reactive transport modeling. However, popular Kd-based "reactive" transport models only require contaminant concentrations in the pore fluids as initial conditions, and minerals implicitly represent infinite sources and sinks in these models. That situation results in a general neglect of mineralogical characterization in site investigations. This study uses a coupled multi-component reactive mass transport model to predict the natural attenuation of a ground water plume at a uranium mill tailings site in western USA. Numerous ground water geochemistry data are available at this site, but mineralogical data are sketchy. Even given the well-defined pore fluid chemistry, variations of secondary mineral species and mineral abundances in the aquifer resulted in significantly different modeling outcomes. Results show that the amount of calcite in the aquifer determines the distances of plume migration. The possible presence of jurbanite, an aluminum sulfate phase, can store acidity temporarily but cause more severe contamination on a later date. The surfaces of iron oxyhydroxides can store significant amounts of sulfate and protons and serve as a second source for prolonged contamination. These simulations under field conditions illustrate that mineralogical compositions are an essential requirement for accurate prediction of contaminant fate and transport.