Microbially mediated redox processes in natural analogues for radioactive waste

Natural analogues allow scientists to investigate biogeochemical processes relevant to radioactive waste disposal that occur on time scales longer than those that may be studied by time-limited laboratory experiments. The Palmottu U-Th deposit in Finland and the Bangombé natural nuclear reactor in Gabon involve the study of natural uranium, and are both considered natural analogues for subsurface radioactive waste disposal. The microbial population naturally present in groundwater may affect the redox conditions, and hence, the radionuclide solubility and migration. Therefore, groundwater samples from the two sites were investigated for microbial populations. The total numbers of cells ranged from 10(4) to 10(6) cells ml(-1). Iron-reducing bacteria (IRB) were the largest culturable microbial population in the Palmottu groundwater and were present at up to 1.3 x 10(5) cells ml(-1). Sulfate-reducing bacteria (SRB) and acetogens could also be cultured from the Palmottu groundwater. The numbers of IRB and SRB were largest in groundwater with the lowest uranium concentrations. Removal of dissolved U(VI) from solution was concomitant with the growth of IRB enrichment cultures and the reduction of iron. The redox buffer in the Palmottu groundwater consists of iron and uranium species, both of which are affected by IRB. IRB and aerobic heterotrophs were cultured from the Bangombé groundwater, where redox potentials are buffered by iron and organic carbon species. Microbial populations similar to those found at Palmottu and Bangombé are found throughout the Fennoscandian Shield, a potential host rock for subsurface radioactive waste disposal. These results confirm that microorganisms can be expected to play a role in stabilizing radioactive waste disposed of in the subsurface by lowering redox potential and immobilizing radionuclides.     

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

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

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.     

Direct coupling of a genome-scale microbial in silico model and a groundwater reactive transport model

The activity of microorganisms often plays an important role in dynamic natural attenuation or engineered bioremediation of subsurface contaminants, such as chlorinated solvents, metals, and radionuclides. To evaluate and/or design bioremediated systems, quantitative reactive transport models are needed. State-of-the-art reactive transport models often ignore the microbial effects or simulate the microbial effects with static growth yield and constant reaction rate parameters over simulated conditions, while in reality microorganisms can dynamically modify their functionality (such as utilization of alternative respiratory pathways) in response to spatial and temporal variations in environmental conditions. Constraint-based genome-scale microbial in silico models, using genomic data and multiple-pathway reaction networks, have been shown to be able to simulate transient metabolism of some well studied microorganisms and identify growth rate, substrate uptake rates, and byproduct rates under different growth conditions. These rates can be identified and used to replace specific microbially-mediated reaction rates in a reactive transport model using local geochemical conditions as constraints. We previously demonstrated the potential utility of integrating a constraint-based microbial metabolism model with a reactive transport simulator as applied to bioremediation of uranium in groundwater. However, that work relied on an indirect coupling approach that was effective for initial demonstration but may not be extensible to more complex problems that are of significant interest (e.g., communities of microbial species and multiple constraining variables). Here, we extend that work by presenting and demonstrating a method of directly integrating a reactive transport model (FORTRAN code) with constraint-based in silico models solved with IBM ILOG CPLEX linear optimizer base system (C library). The models were integrated with BABEL, a language interoperability tool. The modeling system is designed in such a way that constraint-based models targeting different microorganisms or competing organism communities can be easily plugged into the system. Constraint-based modeling is very costly given the size of a genome-scale reaction network. To save computation time, a binary tree is traversed to examine the concentration and solution pool generated during the simulation in order to decide whether the constraint-based model should be called. We also show preliminary results from the integrated model including a comparison of the direct and indirect coupling approaches and evaluated the ability of the approach to simulate field experiment.     

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

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

Intercomparison 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.     

On the selection of primary variables in numerical formulation for modeling multiphase flow in porous media

Selecting the proper primary variables is a critical step in efficiently modeling the highly nonlinear problem of multiphase subsurface flow in a heterogeneous porous-fractured media. Current simulation and ground modeling techniques consist of (1) spatial discretization of mass and/or heat conservation equations using finite difference or finite element methods; (2) fully implicit time discretization; (3) solving the nonlinear, discrete algebraic equations using a Newton iterative scheme. Previous modeling efforts indicate that the choice of primary variables for a Newton iteration not only impacts computational performance of a numerical code, but may also determine the feasibility of a numerical modeling study in many field applications. This paper presents an analysis and general recommendations for selecting primary variables in simulating multiphase, subsurface flow for one-active phase (Richards' equation), two-phase (gas and liquid) and three-phase (gas, water and nonaqueous phase liquid or NAPL) conditions. In many cases, a dynamic variable switching or variable substitution scheme may have to be used in order to achieve optimal numerical performance and robustness. The selection of primary variables depends in general on the sensitivity of the system of equations to the variables selected at given phase and flow conditions. We will present a series of numerical tests and large-scale field simulation examples, including modeling one (active)-phase, two-phase and three-phase flow problems in multi-dimensional, porous-fractured subsurface systems.     

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

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

A serendipitous, long-term infiltration experiment: Water and tritium circulation beneath the CAMBRIC trench at the Nevada Test Site

Underground nuclear weapons testing at the Nevada Test Site introduced numerous radionuclides that may be used subsequently to characterize subsurface hydrologic transport processes in arid climates. In 1965, a unique, 16-year pumping experiment designed to examine radionuclide migration away from the CAMBRIC nuclear test, conducted in the saturated zone beneath Frenchman Flat, Nevada, USA, gave rise to an unintended second experiment involving radionuclide infiltration through the vadose zone, as induced by seepage of pumping effluents beneath an unlined discharge trench. The combined experiments have been reanalyzed using a detailed, three-dimensional numerical model of transient, variably saturated flow and mass transport in a heterogeneous subsurface, tailored specifically for large-scale and efficient calculations. Simulations have been used to estimate tritium travel and residence times in various parts of the system for comparison with observations in wells. Model predictions of mass transport were able to clearly demonstrate radionuclide recycling behavior between the trench and pumping well previously suggested by isotopic age dating information; match travel time estimates for radionuclides moving between the trench, the water table, and monitoring and pumping wells; and provide more realistic ways in which to interpret the pumping well elution curves. Collectively, the results illustrate the utility of integrating detailed numerical modeling with diverse observational data in developing more accurate interpretations of contaminant migration processes.     

Modeling in-situ uranium(VI) bioreduction by sulfate-reducing bacteria

We present a travel-time based reactive transport model to simulate an in-situ bioremediation experiment for demonstrating enhanced bioreduction of uranium(VI). The model considers aquatic equilibrium chemistry of uranium and other groundwater constituents, uranium sorption and precipitation, and the microbial reduction of nitrate, sulfate and U(VI). Kinetic sorption/desorption of U(VI) is characterized by mass transfer between stagnant micro-pores and mobile flow zones. The model describes the succession of terminal electron accepting processes and the growth and decay of sulfate-reducing bacteria, concurrent with the enzymatic reduction of aqueous U(VI) species. The effective U(VI) reduction rate and sorption site distributions are determined by fitting the model simulation to an in-situ experiment at Oak Ridge, TN. Results show that (1) the presence of nitrate inhibits U(VI) reduction at the site; (2) the fitted effective rate of in-situ U(VI) reduction is much smaller than the values reported for laboratory experiments; (3) U(VI) sorption/desorption, which affects U(VI) bioavailability at the site, is strongly controlled by kinetics; (4) both pH and bicarbonate concentration significantly influence the sorption/desorption of U(VI), which therefore cannot be characterized by empirical isotherms; and (5) calcium-uranyl-carbonate complexes significantly influence the model performance of U(VI) reduction.     

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

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

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.     

Modeling microbial-mediated reduction in batch reactors

The governing equations that depict microbially-mediated reduction of heavy metals in the subsurface include a system of coupled nonlinear partial differential equations (PDE's) that describe physical (transport), chemical (sorption), and microbial (reduction/oxidation) processes. The existence of nonlinear reaction terms makes numerical simulations more challenging; however, with the advent of time-splitting solution algorithms, nonlinear reaction terms can be isolated from the convective-dispersive components of the governing transport equations and then solved as a coupled system of nonlinear ordinary differential equations (ODE's). In this paper, four methods are evaluated for solving coupled systems of nonlinear ODE's that describe microbially-mediated reduction/oxidation processes. The evaluation involves a series of comparisons of transient simulations of electron donor oxidation, electron acceptor reduction, and microbial biomass accumulation. The methods evaluation is initiated with a comparison of simulation results obtained with the four methods to those generated with an analytical model. Next, laboratory observations, of nitrite consumption by Nitrobacter winogradski in batch reactors are used in a comparison of batch system simulations generated using each of the four methods and BIOKEMOD (biogeochemical kinetic/equilibrium reaction model). The evaluation finds one of the four methods, the quasi-steady-state approximation (QSSA), to be among the most accurate and easiest to implement. Final validation of the QSSA is performed simulating experimental results of microbially-mediated chromium reductions in batch cultures.     

非饱和土壤轻油污染多相流研究进展

轻油泄漏到地下后 ,在非饱和土壤中水、气、油三相相互作用、相互影响 ,共同构成多相流体系统。在简要分析轻油污染问题的产生、危害以及轻油泄漏在地表以下运移方式的基础上 ,介绍并总结了近些年来国内外土壤轻油污染有关的多相流研究方面的进展 ,包括国外的模型研究成果、国内外试验研究方法和研究成果 ,并指出目前研究中存在的问题和需要进一步探讨的方向     

Database uncertainty as a limiting factor in reactive transport prognosis

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

Humic colloid-borne migration of uranium in sand columns

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

A reaction-based paradigm to model reactive chemical transport in groundwater with general kinetic and equilibrium reactions

This paper presents a reaction-based water quality transport model in subsurface flow systems. Transport of chemical species with a variety of chemical and physical processes is mathematically described by M partial differential equations (PDEs). Decomposition via Gauss-Jordan column reduction of the reaction network transforms M species reactive transport equations into two sets of equations: a set of thermodynamic equilibrium equations representing N(E) equilibrium reactions and a set of reactive transport equations of M-N(E) kinetic-variables involving no equilibrium reactions (a kinetic-variable is a linear combination of species). The elimination of equilibrium reactions from reactive transport equations allows robust and efficient numerical integration. The model solves the PDEs of kinetic-variables rather than individual chemical species, which reduces the number of reactive transport equations and simplifies the reaction terms in the equations. A variety of numerical methods are investigated for solving the coupled transport and reaction equations. Simulation comparisons with exact solutions were performed to verify numerical accuracy and assess the effectiveness of various numerical strategies to deal with different application circumstances. Two validation examples involving simulations of uranium transport in soil columns are presented to evaluate the ability of the model to simulate reactive transport with complex reaction networks involving both kinetic and equilibrium reactions.     

Membrane Bioreactor Model for Removal of Organics from Wastewater

Abstract

The persistence of trace organics in wastewater effluent is a major environmental concern. Possible use of fixed microbial films in wastewater treatment processes is currently an active area of research that may be able to address many of these problems. In the waste effluent, the persistence of trace organics is attributed, in part, to the inability of microbial populations to extract energy from dilute environments at a rate fast enough to sustain themselves. To address this problem, a novel wastewater treatment scheme is considered. On the basis of previous hollow fiber biomass growth studies, we believe that an anaerobic biofilm supported by hollow fibers could achieve greater biomass density than a film grown on traditional impermeable supports. This in turn could lead to improved substrate removal efficiency in a reactor of a given volume. Using this concept, we developed a mathematical model to test the potential of hollow fiber membrane reactors for biodegradation of acetate solution. A computer simulation has shown that it would be possible to remove about 90% from feed solutions containing 0.1 mg-cm^-3 acetate with biomass density 25 mg-cm^-3 in the hollow fiber supported biofilm. More concentrated feeds could be effectively treated if sufficiently high biomass density could be attained. This process, therefore, shows promise in wastewater treatment. The advantages of hollow fiber membrane bioreactors are their high surface-to-volume ratio, separation of cells from flow, and high cell concentration. All of these are essential requirements for optimum utilization of biomass in wastewater treatment. The hollow fiber membrane bioreactor concept, therefore, may provide a new and unique approach to treating organics.     

Groundwater flow, multicomponent transport and biogeochemistry: development and application of a coupled process model

A research tool for modeling the reactive flow and transport of groundwater contaminants in multiple dimensions is presented. Arbitrarily complex coupled kinetic–equilibrium heterogeneous reaction networks, automatic code generation, transfer-function based solutions, parameter estimation, high-resolution methods for advection, and robust solvers for the mixed kinetic–equilibrium chemistry are some of the features of reactive flow and transport (RAFT) that make it a versatile research tool in the modeling of a wide variety of laboratory and field experiments. The treatment of reactions is quite general so that RAFT can be used to model biological, adsorption/desorption, complexation, and mineral dissolution/precipitation reactions among others. The integrated framework involving automated code generation and parameter estimation allows for the development, characterization, and evaluation of mechanistic process models. The model is described and used to solve a problem in competitive adsorption that illustrates some of these features. The model is also used to study the development of an in situ Fe(II)-zone by encouraging the growth of an iron-reducing bacterium with lactate as the electron donor. Such redox barriers are effective in sequestering groundwater contaminants such as chromate and TCE.     

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

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