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.     

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.     

A coupled THC model of the FEBEX in situ test with bentonite swelling and chemical and thermal osmosis

The performance assessment of a geological repository for radioactive waste requires quantifying the geochemical evolution of the bentonite engineered barrier. This barrier will be exposed to coupled thermal (T), hydrodynamic (H), mechanical (M) and chemical (C) processes. This paper presents a coupled THC model of the FEBEX (Full-scale Engineered Barrier EXperiment) in situ test which accounts for bentonite swelling and chemical and thermal osmosis. Model results attest the relevance of thermal osmosis and bentonite swelling for the geochemical evolution of the bentonite barrier while chemical osmosis is found to be almost irrelevant. The model has been tested with data collected after the dismantling of heater 1 of the in situ test. The model reproduces reasonably well the measured temperature, relative humidity, water content and inferred geochemical data. However, it fails to mimic the solute concentrations at the heater-bentonite and bentonite-granite interfaces because the model does not account for the volume change of bentonite, the CO(2)(g) degassing and the transport of vapor from the bentonite into the granite. The inferred HCO(3)(-) and pH data cannot be explained solely by solute transport, calcite dissolution and protonation/deprotonation by surface complexation, suggesting that such data may be affected also by other reactions.     

Multicomponent simulation of wastewater-derived nitrogen and carbon in shallow unconfined aquifers. II. Model application to a field site

A multicomponent reactive transport model as presented by MacQuarrie and Sudicky [MacQuarrie, K.T.B., Sudicky, E.A., this volume. Multicomponent simulation of wastewater-derived nitrogen and carbon in shallow unconfined aquifers: I. Model formulation and performance, J. Contam. Hydrol.] is applied to a well-studied wastewater plume in a sandy aquifer near Cambridge, Ontario. Domestic wastewater is released into the unsaturated zone via a drain field at a depth of about 0.8 m. The physical transport parameters for the model are obtained by simulating a non-reactive solute, while kinetic input data for the nitrogen and carbon reaction network are obtained from the literature. The model shows that the wastewater-loading rate has little influence on the moisture content in the unsaturated zone, thus oxygen diffusion in the air phase is an important transport mechanism. The model results are in general agreement with the field-determined moisture and oxygen profiles near the drain field. The simulation results show that oxidation of ammonium and dissolved organic carbon (DOC) goes to completion in the 1.5-m distance between the drain field and the water table, and that calcite dissolution limits the pH reduction to about 0.2 units. The model-predicted nitrate concentrations in the core of the plume are in the range of 20-25 mg N/l and are in good agreement with the field data. Overall, the results for the major reactive species from the model simulation agree well with the geochemical data obtained below the drain field and it is concluded that the major physical and biochemical processes have been correctly captured in the current model formulation.     

Sulfur and Carbon Cycling in a Flue Gas Desulfurization Sludge Disposal Site

Abstract

Products of a power plant flue gas desulfurization scrubber are discharged into a pond as sludge consisting of calcite (initial δ¹³C 3.2–3.8‰), gypsum (initial δ³⁴S 7.6–8.6‰), and aqueous solution. Reducing conditions exist below a boundary that appears to move vertically as a function of changes in pond water level. Under reducing conditions, bacteria partially reduce aqueous sulfate to low-δ³⁴S sulfide, consuming organic carbon and generating low-δ¹³C bicarbonate. Under oxidizing conditions, sulfide is converted to sulfate, leading to calcite dissolution, gypsum precipitation, and isotopic re-equilibration of remaining calcite with dissolved bicarbonate near the pond surface. The gypsum has δ³⁴S near 6‰, and calcite has δ¹³C as low as -1.7‰; the changes from initial values correspond to predictions based on isotopic balance and reaction stoichiometry. The pond largely contains the products of bacterial reduction. After the pond is abandoned, these products may adversely affect attempts to revegetate the site. Future bacterial reduction may be best controlled by dewatering and limiting the supply of organic matter in percolating surface water.     

Competitive sorption and diffusion of chromate and sulphate in a flow system with goethite in gel beads

Column experiments and model simulations were employed to evaluate the processes involved in multicomponent solute transport in a system with heterogeneous flow. Column experiments were performed with goethite embedded in polyacrylamide gel beads. The gel forms an immobile water region that can be accessed by diffusion. A two-region transport model with diffusion into spheres was combined with a surface complexation model to predict reactive transport in the goethite-gel bead system. Chromate and sulphate breakthrough curves were measured in a set of transport experiments, along with corresponding changes in the pH of the effluent. Sorption and transport of sulphate and chromate in separate columns were predicted from independently measured sorption parameters. The model overestimated the pH changes in the effluent, possibly because of proton buffering by the polyacrylamide gel. The effect of competitive sorption on transport was examined in experiments with both anions present. The model predicted the effect of competition very well in a system initially equilibrated with sulphate, followed by infiltration with chromate. However, when sulphate was infiltrated after equilibration with chromate, chromate desorption and sulphate adsorption were clearly overestimated by the transport model. The exchange between the more strongly bound chromate and the sulphate added subsequently may be too slow to cause a substantial chromate peak in the effluent. This suggests that the local equilibrium assumption was not applicable in this case.     

Parameters that control the cleanup of fractured permeable aquifers

This study develops a modeling approach for simulating and evaluating entrapped light nonaqueous-phase liquid (light NAPL-LNAPL) dissolution and transport of the solute in a fractured permeable aquifer (FPA). The term FPA refers to an aquifer made of porous blocks of high permeability that embed fractures. The fracture network is part of the domain characterized by high permeability and negligible storage. Previous studies show that sandstone aquifers often represent FPAs. The basic model developed in this study is a two-dimensional (2-D) model of permeable blocks that embed oblique equidistant fractures with constant aperture and orientation. According to this model, two major parameters govern NAPL dissolution and transport of the solute. These parameters are: 1) the dimensionless interphase mass transfer coefficient, K(f0), and 2) the mobility number, N(M0). These parameters represent measures of heterogeneity affecting flow, NAPL dissolution, and transport of the solute in the domain. The parameter K(f0) refers to the rate at which organic mass is transferred from the NAPL into the water phase. The parameter N(M0) represents the ratio of flow through the porous blocks to flow through the fracture network in regions free of entrapped NAPL. It also provides a measure of groundwater flow bypassing regions contaminated by entrapped NAPL. In regions contaminated by entrapped NAPL our simulations have often indicated very low permeability of the porous blocks, enabling a significant increase of the fracture flow at the expense of the permeable block flow. Two types of constitutive relationships also affect the rate of FPA cleanup: 1) the relationship between the saturation of the entrapped NAPL and the permeability of the porous blocks, and 2) the relationships representing effects of the entrapped NAPL saturation and the permeable block flow velocity on rates of interphase mass transfer. This study provides basic tools for evaluating the characteristics of pump-and-treat cleanup of FPAs by referring to sets of parameters and constitutive relationships typical of FPAs. The numerical simulations carried out in this study show that at high initial saturation of the entrapped NAPL, during initial stages of the FPA cleanup the contaminant concentration increases, but later it decreases. This phenomenon originates from significant groundwater bypassing the NAPL entrapped in the permeable blocks via the fracture network.     

Kinematic analysis of solute mass flows in rock fractures with spatially random parameters

Field data of physical properties in heterogeneous crystalline bedrock, like porosity and fracture aperture, is associated with uncertainty that can have a significant impact on the analysis of solute transport in rock fractures. Solutions to the central temporal moments of the residence time probability density function (PDF) are derived in a closed form for a solute Dirac pulse. The solutions are based on a model that takes into account advection along the fracture plane, diffusion into the rock matrix and sorption kinetics in the rock matrix. The most relevant rock properties including fracture aperture and several matrix properties as well as flow velocity are assumed to be spatially random along transport pathways. The mass transport is first solved in a general form along one-dimensional pathways, but the results can be extended to multi-dimensional flows simply by substituting the expected travel time for inert water parcels. Based on data obtained with rock samples taken at Asp? Hard Rock Laboratory in Sweden, the solutions indicate that the heterogeneity of the rock properties contributes to increasing significantly both the variance and the skewness of the residence time probability density function for a pulse travelling in a fracture. The Asp? data suggests that the bias introduced in the variance of the residence time PDF by neglecting the effect of heterogeneity of the rock properties on the radionuclide migration is very large for fractures thinner than a few tenths of a millimetre.     

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.     

Interpretation of injection-withdrawal tracer experiments conducted between two wells in a large single fracture

Tracer experiments conducted using a flow field established by injecting water into one borehole and withdrawing water from another are often used to establish connections and investigate dispersion in fractured rock. As a result of uncertainty in the uniqueness of existing models used for interpretation, this method has not been widely used to investigate more general transport processes including matrix diffusion or advective solute exchange between mobile and immobile zones of fluid. To explore the utility of the injection-withdrawal method as a general investigative tool and with the intent to resolve the transport processes in a discrete fracture, two tracer experiments were conducted using the injection-withdrawal configuration. The experiments were conducted in a fracture which has a large aperture (>500 microm) and horizontally pervades a dolostone formation. One experiment was conducted in the direction of the hydraulic gradient and the other in the direction opposite to the natural gradient. Two tracers having significantly different values of the free-water diffusion coefficient were used. To interpret the experiments, a hybrid numerical-analytical model was developed which accounts for the arcuate shape of the flow field, advection-dispersion in the fracture, diffusion into the matrix adjacent to the fracture, and the presence of natural flow in the fracture. The model was verified by comparison to a fully analytical solution and to a well-known finite-element model. Interpretation of the tracer experiments showed that when only one tracer, advection-dispersion, and matrix diffusion are considered, non-unique results were obtained. However, by using multiple tracers and by accounting for the presence of natural flow in the fracture, unique interpretations were obtained in which a single value of matrix porosity was estimated from the results of both experiments. The estimate of porosity agrees well with independent measurements of porosity obtained from core samples. This suggests that: (i) the injection-withdrawal method is a viable tool for the investigation of general transport processes provided all relevant experimental conditions are considered and multiple conservative tracers are used; and (ii) for the conditions of the experiments conducted in this study, the dominant mechanism for exchange of solute between the fracture and surrounding medium is matrix diffusion.     

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.     

Simulation and analysis of solute transport in 2D fracture/pipe networks: the SOLFRAC program

The Time Domain Random Walk (TDRW) method has been recently developed by Delay and Bodin [Delay, F. and Bodin, J., 2001. Time domain random walk method to simulate transport by advection-dispersion and matrix diffusion in fracture networks. Geophys. Res. Lett., 28(21): 4051-4054.] and Bodin et al. [Bodin, J., Porel, G. and Delay, F., 2003c. Simulation of solute transport in discrete fracture networks using the time domain random walk method. Earth Planet. Sci. Lett., 6566: 1-8.] for simulating solute transport in discrete fracture networks. It is assumed that the fracture network can reasonably be represented by a network of interconnected one-dimensional pipes (i.e. flow channels). Processes accounted for are: (1) advection and hydrodynamic dispersion in the channels, (2) matrix diffusion, (3) diffusion into stagnant zones within the fracture planes, (4) sorption reactions onto the fracture walls and in the matrix, (5) linear decay, and (6) mass sharing at fracture intersections. The TDRW method is handy and very efficient in terms of computation costs since it allows for the one-step calculation of the particle residence time in each bond of the network. This method has been programmed in C++, and efforts have been made to develop an efficient and user-friendly software, called SOLFRAC. This program is freely downloadable at the URL (labo.univ-poitiers.fr/hydrasa/intranet/telechargement.htm). It calculates solute transport into 2D pipe networks, while considering different types of injections and different concepts of local dispersion within each flow channel. Post-simulation analyses are also available, such as the mean velocity or the macroscopic dispersion at the scale of the entire network. The program may be used to evaluate how a given transport mechanism influences the macroscopic transport behaviour of fracture networks. It may also be used, as is the case, e.g., with analytical solutions, to interpret laboratory or field tracer test experiments performed in single fractures.     

Sulfur and carbon cycling in a flue gas desulfurization sludge disposal site

Products of a power plant flue gas desulfurization scrubber are discharged into a pond as sludge consisting of calcite (initial delta13C 3.2-3.8 per thousand), gypsum (initial delta34S 7.6-8.6 per thousand), and aqueous solution. Reducing conditions exist below a boundary that appears to move vertically as a function of changes in pond water level. Under reducing conditions, bacteria partially reduce aqueous sulfate to low-delta34S sulfide, consuming organic carbon and generating low-delta13C bicarbonate. Under oxidizing conditions, sulfide is converted to sulfate, leading to calcite dissolution, gypsum precipitation, and isotopic re-equilibration of remaining calcite with dissolved bicarbonate near the pond surface. The gypsum has delta34S near 6 per thousand, and calcite has delta13C as low as -1.7 per thousand; the changes from initial values correspond to predictions based on isotopic balance and reaction stoichiometry. The pond largely contains the products of bacterial reduction. After the pond is abandoned, these products may adversely affect attempts to revegetate the site. Future bacterial reduction may be best controlled by dewatering and limiting the supply of organic matter in percolating surface water.     

Evaluation of tritiated water diffusion through the Toarcian clayey formation of the Tournemire experimental site (France)

Through-diffusion experiments with tritiated water were performed on argillaceous samples from various zones of the Tournemire test site. It was intended to evaluate the homogeneity of the transport property of unfracturated samples and the influence of the orientation and the nature of the samples (presence of an opened fracture or a pre-existing tectonic fracture filled with calcite and pyrite). Homogeneous values of the tritiated water (HTO) effective diffusion coefficients were deduced from experiments carried out when diffusion occurred parallel to the stratigraphic bedding, with an apparent sensitivity to experimental conditions. Anisotropy was significant, De(HTO) perpendicular to the bedding being 1/3 lower than that parallel to the bedding. The observed fractures of the samples created by mechanical stress and partial dehydration during sawing and the presence of a pre-existing opened fracture did not affect the effective diffusion coefficients of tritiated water, which is probably due to the healing ability of the clayey medium during the re-saturation phases of the equilibrium steps performed prior to the diffusion experiments. On the contrary, a significant decrease of this transport parameter was induced by the occurrence of a pre-existing tectonic fracture, which was assigned to the dense structure of the filling phases.     

Utilization of fly ash to improve the quality of the acid mine drainage generated by oxidation of a sulphide-rich mining waste: column experiments

The production of Acid Mine Drainage (AMD) as a result of the oxidative dissolution of sulphides is one of the main pollution problems affecting natural watercourses in mining environments with sulphide-rich residues. In this work, the generation of AMD was prevented by means of the addition of fly ash to sulphide-rich residues in non-saturated column experiments. A column experiment filled with a pyrite-rich sludge with artificial irrigation leached acid drainages (pH approx. 2) containing high concentrations of sulphate, iron and other metals. However, non-saturated column experiments filled with pyritic-rich sludge and fly ash drained leachates characterized by alkaline pH (pH up to 10), low sulphate concentration, and lack of iron and other metals in solution. The pyrite oxidative dissolution at high pH, as a consequence of the leaching of fly ash, favours the metal precipitation inside the column (mainly iron), the coating of pyrite grains, and the attenuation of the oxidation process, resulting in a great improvement in the quality of the leachates.     

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.     

Flow and transport in fractured tuff at Yucca Mountain: numerical experiments on fast preferential flow mechanisms

Recent discovery of bomb-related ³⁶Cl at depth in fractured tuff in the unsaturated zone at the Yucca Mountain candidate high-level waste (HLW) repository site has called into question the usual modeling assumptions based on the equivalent continuum model (ECM). A dual continuum model (DCM) for simulating transient flow and transport at Yucca Mountain is developed. In order to ensure properly converged flow solutions, which are used in the transport simulation, a new flow solution convergence criteria is derived. An extensive series of simulation studies is presented which indicates that rapid movement of solute through the fractures will not occur unless there are intense episodic infiltration events. Movement of solute in the environs of the repository is enhanced if the properties of the tuff layer at the repository horizon are modified from current best-estimate values. Due to a large advective–dispersive coupling between the matrix and fractures, the matrix acts as a major buffer which inhibits rapid transport along the fractures. Consequently, fast movement of solutes through the fractures to the repository depth can only be explained if the matrix–fracture coupling term is significantly reduced from a value that would be calculated on the basis of data currently available.     

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.     

Inverse modeling of BTEX dissolution and biodegradation at the Bemidji, MN crude-oil spill site

The U.S. Geological Survey (USGS) solute transport and biodegradation code BIOMOC was used in conjunction with the USGS universal inverse modeling code UCODE to quantify field-scale hydrocarbon dissolution and biodegradation at the USGS Toxic Substances Hydrology Program crude-oil spill research site located near Bemidji, MN. This inverse modeling effort used the extensive historical data compiled at the Bemidji site from 1986 to 1997 and incorporated a multicomponent transport and biodegradation model. Inverse modeling was successful when coupled transport and degradation processes were incorporated into the model and a single dissolution rate coefficient was used for all BTEX components. Assuming a stationary oil body, we simulated benzene, toluene, ethylbenzene, m,p-xylene, and o-xylene (BTEX) concentrations in the oil and ground water, respectively, as well as dissolved oxygen. Dissolution from the oil phase and aerobic and anaerobic degradation processes were represented. The parameters estimated were the recharge rate, hydraulic conductivity, dissolution rate coefficient, individual first-order BTEX anaerobic degradation rates, and transverse dispersivity. Results were similar for simulations obtained using several alternative conceptual models of the hydrologic system and biodegradation processes. The dissolved BTEX concentration data were not sufficient to discriminate between these conceptual models. The calibrated simulations reproduced the general large-scale evolution of the plume, but did not reproduce the observed small-scale spatial and temporal variability in concentrations. The estimated anaerobic biodegradation rates for toluene and o-xylene were greater than the dissolution rate coefficient. However, the estimated anaerobic biodegradation rates for benzene, ethylbenzene, and m,p-xylene were less than the dissolution rate coefficient. The calibrated model was used to determine the BTEX mass balance in the oil body and groundwater plume. Dissolution from the oil body was greatest for compounds with large effective solubilities (benzene) and with large degradation rates (toluene and o-xylene). Anaerobic degradation removed 77% of the BTEX that dissolved into the water phase and aerobic degradation removed 17%. Although goodness-of-fit measures for the alternative conceptual models were not significantly different, predictions made with the models were quite variable.     

Migration behavior of naturally occurring radionuclides at the Nopal I uranium deposit, Chihuahua, Mexico

Oxidation of pyrite at the Nopal I uranium deposit, Peña Blanca district, Chihuahua, Mexico has resulted in the formation of Fe-oxides/hydroxides. Anomalous U concentrations (i.e. several hundred to several thousand ppm) measured in goethite, hematite, and amorphous Fe-oxyhydroxides in a major fracture that crosscuts the deposit and the absence of U minerals in the fracture suggest that U was retained during secondary mineral growth or sorbed on mineral surfaces. Mobilization and transport of U away from the deposit is suggested by decreasing U concentrations in fracture-infilling materials and in goethite and hematite with distance from the deposit. Greater than unity ²³⁴U/²³⁸U activity ratios measured in fracture-infilling materials indicate relatively recent ( < 1 Ma) U uptake from fluids that carried excess ²³⁴U. Systematic decreases in ²³⁴U/²³⁸U activity ratios of fracture materials with distance from the deposit suggest a multistage mobilization process, such as remobilization of U from ²³⁴U-enriched infill minerals or differential or diminished transport of U-bearing solutions containing excess ²³⁴U.