Multidimensional, multicomponent, subsurface reactive transport in nonuniform velocity fields: code verification using an advective reactive streamtube approach

High performance computing has made possible the development of high resolution, multidimensional, multicomponent reactive transport models that can be used to analyze complex geochemical environments. However, as increasingly complex processes are included in these models, the accuracy of the numerical formulation coupling the nonlinear processes becomes difficult to verify. Analytical solutions are not available for realistically complex problems and benchmark solutions are not generally available for specific problems. We present an advective reactive streamtube (ARS) transport technique that efficiently provides accurate solutions of nonlinear multicomponent reactive transport in nonuniform multidimensional velocity fields. These solutions can be compared with results from Eulerian-based advection-dispersion-reaction models to evaluate the accuracy of the numerical formulation used. The ARS technique includes mixed equilibrium and kinetic complexation and precipitation-dissolution reactions subject to the following assumptions: (1) transport is purely advective (i.e., no explicit diffusion or dispersion), and (2) chemistry is described by a canonical system of reactions that evolves with time and is unaffected by position in space. Results from the ARS technique are compared with results from the massively parallel, multicomponent reactive transport model MCTRACKER on a test problem involving irreversible oxidation of organic carbon and reaction of the oxidation products with two immobile mineral phases, gypsum and calcite, and fifteen aqueous complexes. Truncation error, operator splitting error, and the nonlinear transformation of these errors in the high-resolution reactive transport model are identified for this problem.     

Transport modeling of copper and cadmium with linear and nonlinear retardation factors

The predictive accuracy of using the one-dimensional advection–dispersion equation to evaluate the fate and transport of solute in a soil column is usually dependent on the proper determination of chemical retardation factors. Typically, the distribution coefficient (K_d) obtained by fitting the linear sorption isotherm has been extensively used to consider general geochemical reactions on solute transport in a low-concentration range. However, the linear distribution coefficient cannot be adequately utilized to describe the solute fate at a higher concentration level. This study employed the nonlinear equilibrium-controlled sorption parameters to determine the retardation factor used in column leaching experiments. Copper and cadmium transportation in a lateritic silty-clay soil column was examined. Through the explicit finite-difference calculations with a third-order total-variation-diminishing (TVD) numerical solution scheme, all results of the theoretical copper and cadmium breakthrough curves (BTCs) simulated by using the Freundlich nonlinear retardation factors revealed good agreement with the experimental observations.     

Modelling of physical and reactive processes during biodegradation of a hydrocarbon plume under transient groundwater flow conditions

Numerical experiments of non-reactive and reactive transport were carried out to quantify the influence of a seasonally varying, transient flow field on transport and natural attenuation at a hydrocarbon-contaminated field site. Different numerical schemes for solving advective transport were compared to assess their capability to model low transversal dispersivities in transient flow fields. For the field site, it is shown that vertical plume spreading is largely inhibited, particularly if sorption is taken into account. For the reactive simulations, a biodegradation reaction module for the geochemical transport model PHT3D was developed. Results of the reactive transport simulations show that under the site-specific conditions the temporal variations in groundwater flow do, to a modest extent, affect average biodegradation rates and average total (dissolved) contaminant mass in the aquifer. The model simulations demonstrate that the seasonal variability in groundwater flow only results in significantly enhanced biodegradation rates when a differential sorption of electron donor (toluene) and electron acceptor (sulfate) is assumed.     

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.     

Numerical modeling of humic colloid borne americium (III) migration in column experiments using the transport/speciation code K1D and the KICAM model

The humic colloid borne Am(III) transport was investigated in column experiments for Gorleben groundwater/sand systems. It was found that the interaction of Am with humic colloids is kinetically controlled, which strongly influences the migration behavior of Am(III). These kinetic effects have to be taken into account for transport/speciation modeling. The kinetically controlled availability model (KICAM) was developed to describe actinide sorption and transport in laboratory batch and column experiments. Application of the KICAM requires a chemical transport/speciation code, which simultaneously models both kinetically controlled processes and equilibrium reactions. Therefore, the code K1D was developed as a flexible research code that allows the inclusion of kinetic data in addition to transport features and chemical equilibrium. This paper presents the verification of K1D and its application to model column experiments investigating unimpeded humic colloid borne Am migration. Parmeters for reactive transport simulations were determined for a Gorleben groundwater system of high humic colloid concentration (GoHy 2227). A single set of parameters was used to model a series of column experiments. Model results correspond well to experimental data for the unretarded humic borne Am breakthrough.     

Numerical modeling of coupled variably saturated fluid flow and reactive transport with fast and slow chemical reactions

The couplings among chemical reaction rates, advective and diffusive transport in fractured media or soils, and changes in hydraulic properties due to precipitation and dissolution within fractures and in rock matrix are important for both nuclear waste disposal and remediation of contaminated sites. This paper describes the development and application of LEHGC2.0, a mechanistically based numerical model for simulation of coupled fluid flow and reactive chemical transport, including both fast and slow reactions in variably saturated media. Theoretical bases and numerical implementations are summarized, and two example problems are demonstrated. The first example deals with the effect of precipitation/dissolution on fluid flow and matrix diffusion in a two-dimensional fractured media. Because of the precipitation and decreased diffusion of solute from the fracture into the matrix, retardation in the fractured medium is not as large as the case wherein interactions between chemical reactions and transport are not considered. The second example focuses on a complicated but realistic advective-dispersive-reactive transport problem. This example exemplifies the need for innovative numerical algorithms to solve problems involving stiff geochemical reactions.     

Chemistry of surface water at a volcanic summit area, Norikura, central Japan: multivariate statistical approach.

Many hydrochemical studies on chemical formation of shallow ground water have been reported as results of water-rock interaction, and contamination of paleo-brine or human activities, whereas the preliminary formation of precipitation source in the recharged region has not been established yet. The purpose of this research work is to clarify the geochemical process of water formation from a water source unpolluted by seawater or human activity. Norikura volcano, located in western part of central Japan provided a suitable source for this research purpose, and hence chemical compositions of water samples from the summit and the mountainside area of Norikura volcano were determined. Most samples in the summit area showed very low electrical conductivity, and lower than 12 microS/cm. On the basis of the chemical compositions, principal component analysis (PCA) and factor analysis (FA), such as kinds of multivariate statistical techniques were used to extract geochemical factors affecting hydrochemical process. As a result, three factors were extracted. The first factor showed high loading on K+, Ca2+, SO2 and SiO2, and this factor was interpreted due to influence of the chemical interaction between acidic precipitated water and rocks. The second factor showed high loading on Na+ and Cl-, and it was assumed to be an influence of seawater salt. The third factor showed loading on NO3-, and it was interpreted to be caused by biochemical effect of vegetation. The proportionate contributions of these factors to the evolution of water chemical composition were found to be 45%, 20%, and 10% for factors 1, 2 and 3, respectively. The same exploration at the mountainside of Norikura volcano revealed that the chemical variances of the non-geothermal water samples were highly influenced by water-rock interactions. The silicate dissolution showed 45% contribution for all chemical variances, while the adsorption of Ca2+ and Mg2+ by precipitation or ion exchange showed 20% contribution. The seawater salt influence or biochemical effect was statistically negligible in this area. The clear differentiation of geochemical process on water formation was found between the summit area and the mountainside area.     

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.     

A numerical model for transient-hysteretic flow and solute transport in unsaturated porous media

A two-dimensional flow and transport model was developed for simulating transient water flow and nonreactive solute transport in heterogeneous, unsaturated porous media containing air and water. The model is composed of a unique combination of robust and accurate numerical algorithms for solving the Richards', Darcy flux, and advection-dispersion equations. The mixed form of Richards' equation is solved using a finite-element formulation and a modified Picard iteration scheme. Mass lumping is employed to improve solution convergence and stability behavior. The flow algorithm accounts for hysteresis in the pressure head-water content relationship. Darcy fluxes are approximated with a Galerkin and Petrov-Galerkin finite-element method developed for random heterogeneous porous media. The transport equation is solved using an Eulerian-Lagrangian method. A multi-step, fourth-order Runge-Kutta, reverse particle tracking technique and a quadratic-linear interpolation scheme are shown to be superior for determining the advective concentration. A Galerkin finite-element method is used for approximating the dispersive flux. The unsaturated flow and transport model was applied to a variety of rigorous problems and was found to produce accurate, mass-conserving solutions when compared to analytical solutions and published numerical results.     

Multi-component reactive transport modeling of natural attenuation of an acid groundwater plume at a uranium mill tailings site

Natural attenuation of an acidic plume in the aquifer underneath a uranium mill tailings pond in Wyoming, USA was simulated using the multi-component reactive transport code PHREEQC. A one-dimensional model was constructed for the site and the model included advective-dispersive transport, aqueous speciation of 11 components, and precipitation-dissolution of six minerals. Transport simulation was performed for a reclamation scenario in which the source of acidic seepage will be terminated after 5 years and the plume will then be flushed by uncontaminated upgradient groundwater. Simulations show that successive pH buffer reactions with calcite, Al(OH)3(a), and Fe(OH)3(a) create distinct geochemical zones and most reactions occur at the boundaries of geochemical zones. The complex interplay of physical transport processes and chemical reactions produce multiple concentration waves. For SO4(2-) transport, the concentration waves are related to advection-dispersion, and gypsum precipitation and dissolution. Wave speeds from numerical simulations compare well to an analytical solution for wave propagation.     

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.     

Contributions of advective and diffusive oxygen transport through multilayer composite caps over mine waste

The relative contributions of four mechanisms of oxygen transport in multilayer composite (MLC) caps placed over oxygen-consuming mine waste were evaluated using numerical and analytical methods. MLC caps are defined here as caps consisting of earthen and geosynthetic (polymeric) components where a composite barrier layer consisting of a geomembrane (1-2 mm thick polymeric sheet) overlying a clay layer is the primary barrier to transport. The transport mechanisms that were considered are gas-phase advective transport, gas-phase diffusive transport, liquid-phase advective transport via infiltrating precipitation and liquid-phase diffusive transport. A numerical model was developed to simulate gas-phase advective-diffusive transport of oxygen through a multilayer cap containing seven layers. This model was also used to simulate oxygen diffusion in the liquid phase. An approximate analytical method was used to compute the advective flux of oxygen in the liquid phase. The numerical model was verified for limiting cases using an analytical solution. Comparisons were also made between model predictions and field data for earthen caps reported by others. Results of the analysis show that the dominant mechanism for oxygen transport through MLC caps is gas-phase diffusion. For the cases that were considered, the gas-phase diffusive flux typically comprises at least 99% of the total oxygen flux. Thus, designers of MLC caps should focus on design elements and features that will limit diffusion of gas-phase oxygen.     

Operator-splitting errors in coupled reactive transport codes for transient variably saturated flow and contaminant transport in layered soil profiles

One possible way of integrating subsurface flow and transport processes with (bio)geochemical reactions is to couple by means of an operator-splitting approach two completely separate codes, one for variably-saturated flow and solute transport and one for equilibrium and kinetic biogeochemical reactions. This paper evaluates the accuracy of the operator-splitting approach for multicomponent systems for typical soil environmental problems involving transient atmospheric boundary conditions (precipitation, evapotranspiration) and layered soil profiles. The recently developed HP1 code was used to solve the coupled transport and chemical equations. For steady-state flow conditions, the accuracy was found to be mainly a function of the adopted spatial discretization and to a lesser extent of the temporal discretization. For transient flow situations, the accuracy depended in a complex manner on grid discretization, time stepping and the main flow conditions (infiltration versus evaporation). Whereas a finer grid size reduced the numerical errors during steady-state flow or the main infiltration periods, the errors sometimes slightly increased (generally less than 50%) when a finer grid size was used during periods with a high evapotranspiration demand (leading to high pressure head gradients near the soil surface). This indicates that operator-splitting errors are most significant during periods with high evaporative boundary conditions. The operator-splitting errors could be decreased by constraining the time step using the performance index (the product of the grid Peclet and Courant numbers) during infiltration, or the maximum time step during evapotranspiration. Several test problems were used to provide guidance for optimal spatial and temporal discretization.     

Determination of spatially-resolved porosity, tracer distributions and diffusion coefficients in porous media using MRI measurements and numerical simulations

This work is focused on measuring the concentration distribution of a conservative tracer in a homogeneous synthetic porous material and in heterogeneous natural sandstone using MRI techniques, and on the use of spatially resolved porosity data to define spatially variable diffusion coefficients in heterogeneous media. The measurements are made by employing SPRITE, a fast MRI method that yields quantitative, spatially-resolved tracer concentrations in porous media. Diffusion experiments involving the migration of H(2)O into D(2)O-saturated porous media are conducted. One-dimensional spatial distributions of H(2)O-tracer concentrations acquired from experiments with the homogeneous synthetic calcium silicate are fitted with the one-dimensional analytical solution of Fick's second law to confirm that the experimental method provides results that are consistent with expectations for Fickian diffusion in porous media. The MRI-measured concentration profiles match well with the solution for Fick's second law and provide a pore-water diffusion coefficient of 1.75×10(-9)m(2)s(-1). The experimental approach was then extended to evaluate diffusion in a heterogeneous natural sandstone in three dimensions. The relatively high hydraulic conductivity of the sandstone, and the contrast in fluid density between the H(2)O tracer and the D(2)O pore fluid, lead to solute transport by a combination of diffusion and density-driven advection. The MRI measurements of spatially distributed tracer concentration, combined with numerical simulations allow for the identification of the respective influences of advection and diffusion. The experimental data are interpreted with the aid of MIN3P-D - a multicomponent reactive transport code that includes the coupled processes of diffusion and density-driven advection. The model defines local diffusion coefficients as a function of spatially resolved porosity measurements. The D(e) values calculated for the heterogeneous sandstone and used to simulate diffusive and advective transport range from 5.4×10(-12) to 1.0×10(-10)m(2)s(-1). These methods have broad applicability to studies of contaminant migration in geological materials.     

Three-dimensional density-dependent flow and multicomponent reactive transport modeling of chlorinated solvent oxidation by potassium permanganate

A popular method for the treatment of aquifers contaminated with chlorinated solvents is chemical oxidation based on the injection of potassium permanganate (KMnO₄). Both the high density (1025 gL^− 1) and reactivity of the treatment solution influence the fate of permanganate (MnO₄) in the subsurface and affect the degree of contaminant treatment. The MIN3P multicomponent reactive transport code was enhanced to simulate permanganate-based remediation, to evaluate the pathways of MnO₄ utilization, and to assess the role of density contrasts for the delivery of the treatment solution. The modified code (MIN3P-D) provides a direct coupling between density-dependent fluid flow, solute transport, contaminant treatment, and geochemical reactions. The model is used to simulate a field trial of TCE oxidation in a sandy aquifer that is underlain by an aquitard. Three-dimensional simulations are conducted for a coupled reactive system comprised of ten aqueous components, two mineral phases, TCE (dissolved, adsorbed, and NAPL), reactive organic matter, and including ion exchange reactions. Model parameters are constrained by literature data and a detailed data set from the field site under investigation. The general spatial and transient evolution in observed concentrations of the oxidant, dissolved TCE, and reaction products are adequately reproduced by the simulations. The model elucidates the important role of density-induced flow and transport on the distribution of the treatment solution into NAPL containing regions located at the aquifer–aquitard interface. Model results further suggest that reactions that do not directly affect the stability of MnO₄ have a negligible effect on solution density and MnO₄ delivery.     

An analytical solution for vertical transport of volatile chemicals in the vadose zone

An analytical solution is presented for one-dimensional vertical transport of volatile chemicals through the vadose zone to groundwater. The solution accounts for the important transport mechanisms of the steady advection of water and gas, diffusion and dispersion in water and gas, as well as adsorption, and first-order degradation. By assuming a linear, equilibrium partitioning between water, gas and the adsorbed chemical phases, the dependent variable in the mathematical model becomes the total resident concentration. The general solution was derived for cases having a constant initial total concentration over a discrete depth interval and a zero initial total concentration elsewhere. A zero concentration gradient is assumed at the groundwater table. Examples are given to demonstrate the application of the new solution for calculating the case of a non-uniform initial source concentration, and estimating the transport of chemicals to the groundwater and the atmosphere. The solution was also used to verify a numerical code called VLEACH. We discovered an error in VLEACH, and found that the new solution agreed very well with the numerical results by corrected VLEACH. A simplified solution to predict the migration of volatile organic chemical due to the gas density effect has shown that a high source concentration may lead to significant downward advective gas-phase transport in a soil with a high air-permeability.     

Investigating the role of gas bubble formation and entrapment in contaminated aquifers: Reactive transport modelling

In many natural and contaminated aquifers, geochemical processes result in the production or consumption of dissolved gases. In cases where methanogenesis or denitrification occurs, the production of gases may result in the formation and growth of gas bubbles below the water table. Near the water table, entrapment of atmospheric gases during water table rise may provide a significant source of O(2) to waters otherwise depleted in O(2). Furthermore, the presence of bubbles will affect the hydraulic conductivity of an aquifer, resulting in changes to the groundwater flow regime. The interactions between physical transport, biogeochemical processes, and gas bubble formation, entrapment and release is complex and requires suitable analysis tools. The objective of the present work is the development of a numerical model capable of quantitatively assessing these processes. The multicomponent reactive transport code MIN3P has been enhanced to simulate bubble growth and contraction due to in-situ gas production or consumption, bubble entrapment due to water table rise and subsequent re-equilibration of the bubble with ambient groundwater, and permeability changes due to trapped gas phase saturation. The resulting formulation allows for the investigation of complex geochemical systems where microbially mediated redox reactions both produce and consume gases as well as affect solution chemistry, alkalinity, and pH. The enhanced model has been used to simulate processes in a petroleum hydrocarbon contaminated aquifer where methanogenesis is an important redox process. The simulations are constrained by data from a crude oil spill site near Bemidji, MN. Our results suggest that permeability reduction in the methanogenic zone due to in-situ formation of gas bubbles, and dissolution of entrapped atmospheric bubbles near the water table, both work to attenuate the dissolved gas plume emanating from the source zone. Furthermore, the simulations demonstrate that under the given conditions more than 50% of all produced CH(4) partitions to the gas phase or is aerobically oxidised near the water table, suggesting that these processes should be accounted for when assessing the rate and extent of methanogenic degradation of hydrocarbons.     

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.     

Modelling the closure-related geochemical evolution of groundwater at a former uranium mine

A newly developed reactive transport model was used to evaluate the potential effects of mine closure on the geochemical evolution in the aquifer downgradient from a mine site. The simulations were conducted for the K?nigstein uranium mine located in Saxony, Germany. During decades of operation, uranium at the former mine site had been extracted by in situ acid leaching of the ore underground, while the mine was maintained in a dewatered condition. One option for decommissioning is to allow the groundwater level to rise to its natural level, flooding the mine workings. As a result, pore water containing high concentrations of dissolved metals, radionuclides, and sulfate may be released. Additional contamination may arise due to the dissolution of minerals contained in the aquifer downgradient of the mine. On the other hand, dissolved metals may be attenuated by reactions within the aquifer. The geochemical processes and interactions involved are highly non-linear and their impact on the quality of the groundwater and surface water downstream of the mine is not always intuitive. The multicomponent reactive transport model MIN3P, which can describe mineral dissolution-precipitation reactions, aqueous complexation, and oxidation-reduction reactions, is shown to be a powerful tool for investigating these processes. The predictive capabilities of the model are, however, limited by the availability of key geochemical parameters such as the presence and quantities of primary and secondary mineral phases. Under these conditions, the model can provide valuable insight by means of sensitivity analyses.