Process-based reactive transport modeling of a permeable reactive barrier for the treatment of mine drainage

Reactive transport modeling of a permeable reactive barrier for the treatment of mine drainage was used to integrate a comprehensive data set including pore water chemistry and solid phase data from several sampling events over a >3-year time period. The simulations consider the reduction of sulfate by the organic carbon-based treatment material and the removal of sulfate and iron by precipitation of reduced mineral phases including iron monosulfides and siderite. Additional parameters constraining the model include dissolved H2S, alkalinity and pH, as well as a suite of solid phase S-fractions identified by extractions. Influences of spatial heterogeneity necessitated the use of a 2-dimensional modeling approach. Simulating observed seasonal fluctuations and long-term changes in barrier reactivity required the use of temperature dependent rate coefficients and a multimodal Monod-type rate expression accounting for the variable reactivity of different organic carbon fractions. Simulated dissolved concentrations of SO4, Fe, H2S, alkalinity and pH, as well as solid phase accumulations of reduced sulfur phases generally compare well to observed trends over 23 months. Spatial variations, seasonal fluctuations and the time-dependent decline in reactivity were also captured. The modeling results generally confirm, and further strengthen, the existing conceptual model for the site. Overall sulfate reduction and S-accumulation rates are constrained with confidence within a factor of 1.5.     

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.     

Reactive transport modeling of processes controlling the distribution and natural attenuation of phenolic compounds in a deep sandstone aquifer

Reactive solute transport modeling was utilized to evaluate the potential for natural attenuation of a contaminant plume containing phenolic compounds at a chemical producer in the West Midlands, UK. The reactive transport simulations consider microbially mediated biodegradation of the phenolic compounds (phenols, cresols, and xylenols) by multiple electron acceptors. Inorganic reactions including hydrolysis, aqueous complexation, dissolution of primary minerals, formation of secondary mineral phases, and ion exchange are considered. One-dimensional (1D) and three-dimensional (3D) simulations were conducted. Mass balance calculations indicate that biodegradation in the saturated zone has degraded approximately 1-5% of the organic contaminant plume over a time period of 47 years. Simulations indicate that denitrification is the most significant degradation process, accounting for approximately 50% of the organic contaminant removal, followed by sulfate reduction and fermentation reactions, each contributing 15-20%. Aerobic respiration accounts for less than 10% of the observed contaminant removal in the saturated zone. Although concentrations of Fe(III) and Mn(IV) mineral phases are high in the aquifer sediment, reductive dissolution is limited, producing only 5% of the observed mass loss. Mass balance calculations suggest that no more than 20-25% of the observed total inorganic carbon (TIC) was generated from biodegradation reactions in the saturated zone. Simulations indicate that aerobic biodegradation in the unsaturated zone, before the contaminant entered the aquifer, may have produced the majority of the TIC observed in the plume. Because long-term degradation is limited to processes within the saturated zone, use of observed TIC concentrations to predict the future natural attenuation may overestimate contaminant degradation by a factor of 4-5.     

A sensibility analysis of model selection in modeling the reactive transport of cesium in crushed granite

We performed a sensibility analysis of model selection in modeling the reactive transport of cesium in crushed granite through model calibration and validation. Based on some solid phase analysis data and kinetic batch experimental results, we hypothesized three two-site sorption models in the LEHGC reactive transport model to fit the breakthrough curves (BTCs) from the corresponding column experiments. The analysis of breakthrough curves shows that both the empirical two-site kinetic linear sorption model and the semi-mechanistic/semi-empirical two-site kinetic surface complexation model, regardless of their complexity, can match our experimental data fairly well under given test conditions. A numerical experiment to further compare the two models shows that they behave differently when the pore velocity is not of the same order of magnitude as our test velocities. This result indicates that further investigations to help determine a better model are needed. We suggest that a multistage column experiment, which tests over the whole range of practical flow velocities, should be conducted to help alleviate inadequate hypothesized models.     

Quantification of biotransformation of chlorinated hydrocarbons in a biostimulation study: added value via stable carbon isotope analysis

Stable carbon isotope analysis of chlorinated aliphatic compounds was performed at an in situ biostimulation pilot test area (PTA) at a site where 1,2-dichloroethane (1,2-DCA) and trichloroethene (TCE) were present in groundwater. Chlorinated products of TCE reductive dechlorination (cis-dichloroethene (cDCE) and vinyl chloride (VC)) were present at concentrations of 17.5 to 126.4 micromol/L. Ethene, a potential degradation product of both 1,2-DCA dihaloelimination and TCE reductive dechlorination was also present in the PTA. Emulsified soybean oil and lactate were added as electron donors to stimulate anaerobic dechlorination in the PTA. Stable carbon isotope analysis provided evidence that dechlorination was occurring in the PTA during biostimulation, and a means of monitoring changes in dechlorination efficiency over the 183 day monitoring period. Stable carbon isotope analysis was also used to determine if ethene production in the PTA was due to dechlorination of TCE, 1,2-DCA, or both. Fractionation factors (alpha) were determined in the laboratory during anaerobic biotransformation of 1,2-DCA via a dihaloelimination reaction in four separate enrichment cultures. These alpha values (as well as the previously published ranges of alpha for the dechlorination of TCE, cDCE and 1,2-DCA) were used, along with isotopic values measured during the pilot test, to derive quantitative estimates of biotransformation during the pilot test. Dechlorination was found to account for 10.7 to 35.9%, 21.9 to 74.9%, and 54.4 to 67.8% of 1,2-DCA, TCE and cDCE concentration loss respectively in the PTA. Stable carbon isotope analysis indicates that dechlorination of 1,2-DCA, TCE and cDCE were all significant processes during the pilot test, while ethene production during the pilot test was dominated by 1,2-DCA dihaloelimination. This study demonstrates how stable carbon isotope analysis can provide more conservative estimates of the extent of biotransformation than do conventional protocols. In addition, in a complex mixed plume such as this, compound specific isotope analysis is shown to be one of the few methods available for clarifying dominant biotransformation pathways where breakdown products are non-exclusive (i.e. ethene).     

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.     

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.     

Contaminant transport in groundwater in the presence of colloids and bacteria: Model development and verification

Colloids and bacteria (microorganisms) naturally exist in groundwater aquifers and can significantly impact contaminant migration rates. A conceptual model is first developed to account for the different physiochemical and biological processes, reaction kinetics, and different transport mechanisms of the combined system (contaminant–colloids–bacteria). All three constituents are assumed to be reactive with the reactions taking place between each constituent and the porous medium and also among the different constituents. A general linear kinetic reaction model is assumed for all reactive processes considered. The mathematical model is represented by fourteen coupled partial differential equations describing mass balance and reaction processes. Two of these equations describe colloid movement and reactions with the porous medium, four equations describe bacterial movement and reactions with colloids and the porous medium, and the remaining eight equations describe contaminant movement and its reactions with bacteria, colloids, and the porous medium. The mass balance equations are numerically solved for two-dimensional groundwater systems using a third-order, total variance-diminishing scheme (TVD) for the advection terms. Due to the complex coupling of the equations, they are solved iteratively each time step until a convergence criterion is met. The model is tested against experimental data and the results are favorable.     

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.     

Natural attenuation of chlorinated ethene compounds: model development and field-scale application at the Dover site

A multi-dimensional and multi-species reactive transport model was developed to aid in the analysis of natural attenuation design at chlorinated solvent sites. The model can simulate several simultaneously occurring attenuation processes including aerobic and anaerobic biological degradation processes. The developed model was applied to analyze field-scale transport and biodegradation processes occurring at the Area-6 site in Dover Air Force Base, Delaware. The model was calibrated to field data collected at this site. The calibrated model reproduced the general groundwater flow patterns, and also, it successfully recreated the observed distribution of tetrachloroethene (PCE), trichloroethene (TCE), dichloroethylene (DCE), vinyl chloride (VC) and chloride plumes. Field-scale decay rates of these contaminant plumes were also estimated. The decay rates are within the range of values that were previously estimated based on lab-scale microcosm and field-scale transect analyses. Model simulation results indicated that the anaerobic degradation rate of TCE, source loading rate, and groundwater transport rate are the important model parameters. Sensitivity analysis of the model indicated that the shape and extent of the predicted TCE plume is most sensitive to transmissivity values. The total mass of the predicted TCE plume is most sensitive to TCE anaerobic degradation rates. The numerical model developed in this study is a useful engineering tool for integrating field-scale natural attenuation data within a rational modeling framework. The model results can be used for quantifying the relative importance of various simultaneously occurring natural attenuation processes.     

Quantifying the physical and chemical mass transfer processes for the fate and transport of Co(II)EDTA in a partially-weathered limestone-shale saprolite

The objective of the research is to quantify the relative contributions of physical and chemical mass transfer to the movement of Co(II/III)EDTA (chelates of Cobalt and Ethylene Diamine Tetraacetic Acid or EDTA) through a limestone-shale saprolite soil. Saprolite is a collective term referring to partially-weathered bedrock. It exists extensively in the subsurface. Because the parent bedding structures are maintained during the weathering process, saprolite soils are characterized by intensive fractures and secondary deposits of minerals such as Al-, Fe- and Mn-oxides on the fracture surfaces. Movement of reactive species through the soils may be influenced by diffusion into the rock matrix, a physical mass transfer (PMT) process, and interfacial chemical reactions, a chemical mass transfer (CMT) process. The PMT and CMT processes are phenomenologically similar but mechanistically different. In this research, previous laboratory observations from a Br and Co(II)EDTA tracer injection into an undisturbed saprolite soil column were used. Mechanistic reactive transport models were formulated to quantify the PMT and CMT processes. The PMT process was independently characterized by using the non-reactive tracer Br. Model parameters thus obtained were subsequently used as constraints to quantify the CMT processes involving Co(II)EDTA and its oxidation product Co(III)EDTA. Our calculations indicated that the PMT rates of the less reactive Co(III)EDTA were comparable with their theoretical CMT rates. In contrast, for the more reactive species Co(II)EDTA, CMT rates are higher than PMT rates. Evaluations of alternative CMT process models further confirmed one of our hypotheses on the basis of previous experimental understandings. The hypothesis suggested that competition from Fe-oxide for Co(II)EDTA may account for the majority of the decrease of Co(III)EDTA effluent concentrations that resulted in the separation of total Co and Co(III)EDTA breakthrough curves. Because Co(III)EDTA is more mobile than Co(II)EDTA in the subsurface, the results of this research suggest independent quantifications of CoEDTA PMT and CMT processes if laboratory results are to be interpreted correctly and scaled up for field and predictive uses.     

Present state and future directions of modeling of geochemistry in hydrogeological systems

A first step towards understanding and controlling the fate and dissemination of radioactive waste is to create a concise and comprehensive theoretical framework for the rather non-linear processes involved--hence, the need for geochemical models. Two classes of geochemical models are commonly used, i.e., static and hydrodynamic models. In contrast to static models, hydrodynamic models combine geochemical reactions with hydrogeological processes such as ground-water flow, diffusion and dispersion. In this review, we examine the present state of geochemical models in terms of included processes, thermodynamic databases, missing phenomena, numerical behavior and performance. It is shown that over the past decade, significant progress has been made with respect to modeling of geochemistry in hydrodynamic systems: this is illustrated by describing several applications. Finally, we focus on the perspectives of geochemical modeling in the assessment of the safety of nuclear waste disposal.     

Dual-domain solute transfer and transport processes: evaluation in batch and transport experiments

Soil macropore networks establish a dual-domain transport scenario in which water and solutes are preferentially channeled through soil macropores while slowly diffusing into and out of the bulk soil matrix. The influence of macropore networks on intra-ped solute diffusion and preferential transport in a soil typical of subsurface-drained croplands in the Midwestern United States was studied in batch- and column-scale experiments. In the batch diffusion studies with soil aggregates, the estimated diffusion radius (length) of the soil aggregates corresponded to the half-spacing of the aggregate fissures, suggesting that the intra-ped fissures reduced the diffusion impedance and preferentially allowed solutes to diffuse into the soil matrix. In the column-scale solute transport experiments, the average diffusion radius (estimated from HYDRUS-2D simulations and a first-order diffusive transfer term) was nearly double that of the batch-scale study. This increase may be attributed to a loss of pore continuity and a compounding of the small diffusion impedance through macropores at the larger scale. The column-scale solute transport experiments also suggest that two preferential networks exist in the soil. At and near soil saturation, a primary network of large macropores (possibly root channels and earthworm burrows) dominate advective transport, causing a high degree of physical and sorption nonequilibrium and simultaneous breakthrough of a nonreactive (bromide) and a reactive (alachlor) solute. As the saturation level decreases, the primary network drains, while transport through smaller macropores (possibly intra-ped features) continues, resulting in a reduced degree of nonequilibrium and separation in the breakthrough curves of bromide and alachlor.     

Effect of cell physicochemical characteristics and motility on bacterial transport in groundwater

The influence of physicochemical characteristics and motility on bacterial transport in groundwater were examined in flow-through columns. Four strains of bacteria isolated from a crystalline rock groundwater system were investigated, with carboxylate-modified and amidine-modified latex microspheres and bromide as reference tracers. The bacterial isolates included a gram-positive rod (ML1), a gram-negative motile rod (ML2), a nonmotile mutant of ML2 (ML2m), and a gram-positive coccoid (ML3). Experiments were repeated at two flow velocities, in a glass column packed with glass beads, and in another packed with iron-oxyhydroxide coated glass beads. Bacteria breakthrough curves were interpreted using a transport equation that incorporates a sorption model from microscopic observation of bacterial deposition in flow-cell experiments. The model predicts that bacterial desorption rate will decrease exponentially with the amount of time the cell is attached to the solid surface. Desorption kinetics appeared to influence transport at the lower flow rate, but were not discernable at the higher flow rate. Iron-oxyhydroxide coatings had a lower-than-expected effect on bacterial breakthrough and no effect on the microsphere recovery in the column experiments. Cell wall type and shape also had minor effects on breakthrough. Motility tended to increase the adsorption rate, and decrease the desorption rate. The transport model predicts that at field scale, desorption rate kinetics may be important to the prediction of bacteria transport rates.     

Long-term geochemical evolution of the near field repository: insights from reactive transport modelling and experimental evidences

The KBS-3 underground nuclear waste repository concept designed by the Swedish Nuclear Fuel and Waste Management Co. (SKB) includes a bentonite buffer barrier surrounding the copper canisters and the iron insert where spent nuclear fuel will be placed. Bentonite is also part of the backfill material used to seal the access and deposition tunnels of the repository. The bentonite barrier has three main safety functions: to ensure the physical stability of the canister, to retard the intrusion of groundwater to the canisters, and in case of canister failure, to retard the migration of radionuclides to the geosphere. Laboratory experiments (< 10 years long) have provided evidence of the control exerted by accessory minerals and clay surfaces on the pore water chemistry. The evolution of the pore water chemistry will be a primordial factor on the long-term stability of the bentonite barrier, which is a key issue in the safety assessments of the KBS-3 concept.In this work we aim to study the long-term geochemical evolution of bentonite and its pore water in the evolving geochemical environment due to climate change. In order to do this, reactive transport simulations are used to predict the interaction between groundwater and bentonite which is simulated following two different pathways: (1) groundwater flow through the backfill in the deposition tunnels, eventually reaching the top of the deposition hole, and (2) direct connection between groundwater and bentonite rings through fractures in the granite crosscutting the deposition hole. The influence of changes in climate has been tested using three different waters interacting with the bentonite: present-day groundwater, water derived from ice melting, and deep-seated brine. Two commercial bentonites have been considered as buffer material, MX-80 and Deponit CA-N, and one natural clay (Friedland type) for the backfill. They show differences in the composition of the exchangeable cations and in the accessory mineral content. Results from the simulations indicate that pore water chemistry is controlled by the equilibrium with the accessory minerals, especially carbonates. pH is buffered by precipitation/dissolution of calcite and dolomite, when present. The equilibrium of these minerals is deeply influenced by gypsum dissolution and cation exchange reactions in the smectite interlayer. If carbonate minerals are initially absent in bentonite, pH is then controlled by surface acidity reactions in the hydroxyl groups at the edge sites of the clay fraction, although its buffering capacity is not as strong as the equilibrium with carbonate minerals. The redox capacity of the bentonite pore water system is mainly controlled by Fe(II)-bearing minerals (pyrite and siderite). Changes in the groundwater composition lead to variations in the cation exchange occupancy, and dissolution–precipitation of carbonate minerals and gypsum. The most significant changes in the evolution of the system are predicted when ice-melting water, which is highly diluted and alkaline, enters into the system. In this case, the dissolution of carbonate minerals is enhanced, increasing pH in the bentonite pore water. Moreover, a rapid change in the population of exchange sites in the smectite is expected due to the replacement of Na for Ca.     

Spatial and temporal features of density-dependent contaminant transport: Experimental investigation and numerical modeling

We investigate the spatial and temporal features of variable-density contaminant plumes migration in porous materials. Our analysis is supported by novel experimental results concerning concentration profiles inside a vertical column setup that has been conceived at CEA to this aim. The experimental method relies on X-ray spectrometry, which allows determining solute profiles as a function of time at several positions along the column. The salient outcomes of the measurements are elucidated, with focus on miscible fluids in homogeneous saturated media. The role of the injected solution molarity is evidenced. As molarity increases, the solutes plume transport progressively deviates from the usual Fickian behavior, and pollutants distribution becomes skewed in the direction dictated by gravity. By resorting to a finite elements approach, we numerically solve the nonlinear equations that rule the pollutants migration: a good agreement is found between the simulated profiles and the experimental data. At high molarity, a strong dependence on initial conditions is found. Finally, we qualitatively explore the (unstable) interfacial dynamics between the dense contaminant plume and the lighter resident fluid that saturates the column, and detail its evolution for finite-duration contaminant injections.     

Migration of uranium(IV)/(VI) in the presence of humic acids in quartz sand: a laboratory column study

The migration behavior of U(IV) and U(VI) in the presence of humic acid was studied in a quartz sand system. Laboratory column experiments were performed using humic acid, U(VI) in humic acid absence, U(IV) and U(VI) in humic acid presence, and for comparison a conservative tracer. In experiments using humic acid, both redox species of U migrate nearly as fast as the conservative tracer. Humic acid accelerates the U(VI) breakthrough compared to the humic acid-free system. There are strong indications for a similar effect on the U(IV) transport. At the same time, a part of U(IV) and U(VI) associated with the humic acid is immobilized in the quartz sand due to humic colloid filtration thus producing a delaying effect. Tailing at a low concentration level was observed upon tracer elution. The experimental breakthrough curves were described by reactive transport modeling using equations for equilibrium and kinetic reactions. The present study demonstrates that humic acids can play an important role in the migration of actinides. As natural organic matter is ubiquitous in aquifer systems, the humic colloid-borne transport of actinides is of high relevance in performance assessment.     

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

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

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.     

Influence of variable chemical conditions on EDTA-enhanced transport of metal ions in mildly acidic groundwater

Adsorption of Ni and Pb on aquifer sediments from Cape Cod, Massachusetts, USA increased with increasing pH and metal-ion concentration. Adsorption could be described quantitatively using a semi-mechanistic surface complexation model (SCM), in which adsorption is described using chemical reactions between metal ions and adsorption sites. Equilibrium reactive transport simulations incorporating the SCMs, formation of metal-ion-EDTA complexes, and either Fe(III)-oxyhydroxide solubility or Zn desorption from sediments identified important factors responsible for trends observed during transport experiments conducted with EDTA complexes of Ni, Zn, and Pb in the Cape Cod aquifer. Dissociation of Pb-EDTA by Fe(III) is more favorable than Ni-EDTA because of differences in Ni- and Pb-adsorption to the sediments. Dissociation of Ni-EDTA becomes more favorable with decreasing Ni-EDTA concentration and decreasing pH. In contrast to Ni, Pb-EDTA can be dissociated by Zn desorbed from the aquifer sediments. Variability in adsorbed Zn concentrations has a large impact on Pb-EDTA dissociation.