Modelling radionuclide transport for time varying flow in a channel network

Water flowrates and flow directions may change over time in the subsurface for a number of reasons. In fractured rocks flow takes place in channels within fractures. Solutes are carried by the advective flow. In addition, solutes may diffuse in and out of stagnant waters in the rock matrix and other stagnant water regions. Sorbing species may sorb on fracture surfaces and on the micropore surfaces in the rock matrix. We present a method by which solute particles can be traced in flowing water undergoing changes in flowrate and direction in a complex channel network where the solutes can also interact with the rock by diffusion in the rock matrix. The novelty of this paper is handling of diffusion in the rock matrix under transient flow conditions. The diffusive processes are stochastic and it is not possible to follow a particle deterministically. The method therefore utilises the properties of a probability distribution function for a tracer moving in a fracture where matrix diffusion is active. The method is incorporated in a three dimensional channel network model. Particle tracking is used to trace out a multitude of flowpaths, each of which consists of a large number of channels within fractures. Along each channel the aperture and velocity as well as the matrix sorption properties can vary. An efficient method is presented whereby a particle can be followed along the variable property flowpath. For stationary flow conditions and a network of channels with advective flow and matrix diffusion, a simple analytical solution for the residence time distribution along each pathway can be used. Only two parameter groups need to be integrated along each path. For transient flow conditions, a time stepping procedure that incorporates a stochastic Monte-Carlo like method to follow the particles along the paths when flow conditions change is used. The method is fast and an example is used for illustrative purposes. It is exemplified by a case where land rises due to glacial rebound. It is shown that the effects of changing flowrates and directions can be considerable and that the diffusive migration in the matrix can have a dominating effect on the results.     

A stochastic multi-channel model for solute transport--analysis of tracer tests in fractured rock

Some of the basic assumptions of the advection-dispersion model (AD-model) are revisited. This model assumes a continuous mixing along the flowpath similar to Fickian diffusion. This implies that there is a constant dispersion length irrespective of observation distance. This is contrary to most field observations. The properties of an alternative model based on the assumption that individual water packages can retain their identity over long distances are investigated. The latter model is called the multi-channel model (MCh-model). Inherent in the latter model is that if the waters in the different pathways are collected and mixed, the "dispersion length" is proportional to distance. The conditions for when non-mixing between adjacent streams can be assumed are explored. The MCh- and AD-models are found to have very similar residence time distributions (RTD) for Peclet numbers larger than 3. A generalized relation between flowrate and residence time is developed, including the so-called cubic law and constant aperture assumptions. The two models extrapolate very differently when there is strong matrix interaction. The AD-model could severely underestimate the effluent concentration of a tracer pulse and overestimate the residence time. The conditions are explored for when in-filling particles in the fracture will not be equilibrated but will act as if there was seemingly a much larger flow wetted surface (FWS). It is found that for strongly sorbing tracers, relatively small particles can act in this way for systems and conditions that are typical of many tracer tests. The assumption that the tracer residence time found by cautiously injecting a small stream of traced water represents the residence time in the whole fracture is explored. It is found that the traced stream can potentially sample a much larger fraction of the fracture than the ratio between the traced flowrate and the total pumped flowrate. The MCh-model was used to simulate some recent tracer tests in what is assumed to be a single fracture at the Asp? Hard rock laboratory in Sweden. Non-sorbing tracers, HTO and Uranin were used to determine the mean residence time and its variance. Laboratory data on diffusion and sorption properties were used to "predict" the RTD of the sorbing tracers. At least 30 times larger FWS or 1000 times larger diffusion or sorption coefficients would be needed to explain the observed BTCs. Some possible reasons for such behavior are also explored.     

Prediction of some in situ tracer tests with sorbing tracers using independent data

Some recent converging tracer tests with sorbing tracers at the Asp? Hard Rock Laboratory in Sweden, the TRUE tests, have been predicted using only laboratory data and hydraulic data from borehole measurements. No model parameters were adjusted to obtain a better fit with the experiments. The independent data were fracture frequency and transmissivity data obtained in the field and laboratory data on sorption and matrix diffusion. Transmissivity measurements in five boreholes in the rock volume containing the region surrounding the injection and collection points show that there is a high frequency of water conducting fractures. Of 162 packed off sections with 0.5 m packer distances, 112 were found to have a transmissivity above the detection limit. The specific flow-wetted surface (FWS) of the rock mass could be estimated from these data. The transmissivities were found to be reasonably well described by a lognormal distribution. Laboratory data on diffusion and sorption properties together with the hydraulic data were used to "predict" the residence time distribution (RTD) of the sorbing tracers. The results were compared with the experimental breakthrough curves. In these experiments, the water residence time is very small compared to the residence time of the sorbing tracers due to their diffusion and sorption within the rock matrix. We thus could neglect the influence of the water residence time in our predictions. Therefore, no information on water residence times or on "dispersion" was needed. The dispersion of the sorbing tracers is caused by the different sorbing tracer residence times in different pathways. The sorbing tracer residence time is determined by the ratio of flowrate to the flow-wetted surface in the different pathways and not by the water residence time. Assuming a three-dimensional flow pattern and using the observed fracture frequency and flowrate distribution, breakthrough curves for three strongly sorbing tracers were predicted. Only the laboratory data, the transmissivity measurements and the pumping flowrate were used in the predictions. No information on the water residence time as obtained by the nonsorbing tracers was used. The predictions were surprisingly accurate.     

An interpretation of potential scale dependence of the effective matrix diffusion coefficient

Matrix diffusion is an important process for solute transport in fractured rock, and the matrix diffusion coefficient is a key parameter for describing this process. Previous studies have indicated that the effective matrix diffusion coefficient values, obtained from a large number of field tracer tests, are enhanced in comparison with local values and may increase with test scale. In this study, we have performed numerical experiments to investigate potential mechanisms behind possible scale-dependent behavior. The focus of the experiments is on solute transport in flow paths having geometries consistent with percolation theories and characterized by multiple local flow loops formed mainly by small-scale fractures. The water velocity distribution through a flow path was determined using discrete fracture network flow simulations, and solute transport was calculated using a previously derived impulse-response function and a particle-tracking scheme. Values for effective (or up-scaled) transport parameters were obtained by matching breakthrough curves from numerical experiments with an analytical solution for solute transport along a single fracture. Results indicate that a combination of local flow loops and the associated matrix diffusion process, together with scaling properties in flow path geometry, seems to be an important mechanism causing the observed scale dependence of the effective matrix diffusion coefficient (at a range of scales).     

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.     

Flow and transport in unsaturated fractured rock: effects of multiscale heterogeneity of hydrogeologic properties

The heterogeneity of hydrogeologic properties at different scales may have different effects on flow and transport processes in a subsurface system. A model for the unsaturated zone of Yucca Mountain, Nevada, is developed to represent complex heterogeneity at two different scales: (1) layer scale corresponding to geologic layering and (2) local scale. The layer-scale hydrogeologic properties are obtained using inverse modeling, based on the available measurements collected from the Yucca Mountain site. Calibration results show a significant lateral and vertical variability in matrix and fracture properties. Hydrogeologic property distributions in a two-dimensional, vertical cross-section of the site are generated by combining the average layer-scale matrix and fracture properties with local-scale perturbations generated using a stochastic simulation method. The unsaturated water flow and conservative (nonsorbing) tracer transport through the cross-section are simulated for different sets of matrix and fracture property fields. Comparison of simulation results indicates that the local-scale heterogeneity of matrix and fracture properties has a considerable effect on unsaturated flow processes, leading to fast flow paths in fractures and the matrix. These paths shorten the travel time of a conservative tracer from the source (repository) horizon in the unsaturated zone to the water table for small fractions of total released tracer mass. As a result, the local-scale heterogeneity also has a noticeable effect on global tracer transport processes, characterized by an average breakthrough curve at the water table, especially at the early arrival time of tracer mass. However, the effect is not significant at the later time after 20% tracer mass reaches the water table. The simulation results also verify that matrix diffusion plays an important role in overall solute transport processes in the unsaturated zone at Yucca Mountain.     

Field investigation into unsaturated flow and transport in a fault: model analyses

Results of a fault test performed in the unsaturated zone of Yucca Mountain, Nevada, were analyzed using a three-dimensional numerical model. The fault was explicitly represented as a discrete feature and the surrounding rock was treated as a dual-continuum (fracture-matrix) system. Model calibration against seepage and water-travel-velocity data suggests that lithophysal cavities connected to fractures can considerably enhance the effective fracture porosity and therefore retard water flow in fractures. Comparisons between simulation results and tracer concentration data also indicate that matrix diffusion is an important mechanism for solute transport in unsaturated fractured rock. We found that an increased fault-matrix and fracture-matrix interface areas were needed to match the observed tracer data, which is consistent with previous studies. The study results suggest that the current site-scale model for the unsaturated zone of Yucca Mountain may underestimate radionuclide transport time within the unsaturated zone, because an increased fracture-matrix interface area and the increased effective fracture porosity arising from lithophysal cavities are not considered in the current site-scale model.     

Variable-density groundwater flow and solute transport in heterogeneous porous media: approaches, resolutions and future challenges

In certain hydrogeological situations, fluid density variations occur because of changes in the solute or colloidal concentration, temperature, and pressure of the groundwater. These include seawater intrusion, high-level radioactive waste disposal, groundwater contamination, and geothermal energy production. When the density of the invading fluid is greater than that of the ambient one, density-driven free convection can lead to transport of heat and solutes over larger spatial scales and significantly shorter time scales than compared with diffusion alone. Beginning with the work of Lord Rayleigh in 1916, thermal and solute instabilities in homogeneous media have been studied in detail for almost a century. Recently, these theoretical and experimental studies have been applied in the study of groundwater phenomena, where the assumptions of homogeneity and isotropy rarely, if ever, apply. The critical role that heterogeneity plays in the onset as well as the growth and/or decay of convective motion is discussed by way of a review of pertinent literature and numerical simulations performed using a variable-density flow and solute transport numerical code. Different styles of heterogeneity are considered and range from continuously "trending" heterogeneity (sinusoidal and stochastic permeability distributions) to discretely fractured geologic media. Results indicate that both the onset of instabilities and their subsequent growth and decay are intimately related to the structure and variance of the permeability field. While disordered heterogeneity tends to dissipate convection through dispersive mixing, an ordered heterogeneity (e.g., sets of vertical fractures) allows instabilities to propagate at modest combinations of fracture aperture and separation distances. Despite a clearer understanding of the processes that control the onset and propagation of instabilities, resultant plume patterns and their migration rates and pathways do not appear amenable to prediction at present. The classical Rayleigh number used to predict the occurrence of instabilities fails, in most cases, when heterogeneous conditions prevail. The incorporation of key characteristics of the heterogeneous permeability field into relevant stability criteria and numerical models remains a challenge for future research.     

Modeling of non-reactive solute transport in fractured clayey till during variable flow rate and time

Fractures and biopores can act as preferential flow paths in clay aquitards and may rapidly transmit contaminants into underlying aquifers. Reliable numerical models for assessment of groundwater contamination from such aquitards are needed for planning, regulatory and remediation purposes. In this investigation, high resolution preferential water-saturated flow and bromide transport data were used to evaluate the suitability of equivalent porous medium (EPM), dual porosity (DP) and discrete fracture/matrix diffusion (DFMD) numerical modeling approaches for assessment of flow and non-reactive solute transport in clayey till. The experimental data were obtained from four large undisturbed soil columns (taken from 1.5 to 3.5 m depth) in which biopores and channels along fractures controlled 96-99% of water-saturated flow. Simulating the transport data with the EPM effective porosity model (FRACTRAN in EPM mode) was not successful because calibrated effective porosity for the same column had to be varied up to 1 order of magnitude in order to simulate solute breakthrough for the applied flow rates between 11 and 49 mm/day. Attempts to simulate the same data with the DP models CXTFIT and MODFLOW/MT3D were also unsuccessful because fitted values for dispersion, mobile zone porosity, and mass transfer coefficient between mobile and immobile zones varied several orders of magnitude for the different flow rates, and because dispersion values were furthermore not physically realistic. Only the DFMD modeling approach (FRACTRAN in DFMD mode) was capable to simulate the observed changes in solute transport behavior during alternating flow rate without changing values of calibrated fracture spacing and fracture aperture to represent the macropores.     

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.     

Long-term oxygen depletion from infiltrating groundwaters: model development and application to intra-glaciation and glaciation conditions

Processes that control the redox conditions in deep groundwaters have been studied. The understanding of such processes in a long-term perspective is important for the safety assessment of a deep geological repository for high-level nuclear waste. An oxidising environment at the depth of the repository would increase the solubility and mobility of many radionuclides, and increase the potential risk for radioactive contamination at the ground surface. Proposed repository concepts also include engineered barriers such as copper canisters, the corrosion of which increases considerably in an oxidising environment compared to prevailing reducing conditions. Swedish granitic rocks are typically relatively sparsely fractured and are best treated as a dual-porosity medium with fast flowing channels through fractures in the rock with a surrounding porous matrix, the pores of which are accessible from the fracture by diffusive transport. Highly simplified problems have been explored with the aim to gain understanding of the underlying transport processes, thermodynamics and chemical reaction kinetics. The degree of complexity is increased successively, and mechanisms and processes identified as of key importance are included in a model framework. For highly complex models, analytical expressions are not fully capable of describing the processes involved, and in such cases the solutions are obtained by numerical calculations. Deep in the rock the main source for reducing capacity is identified as reducing minerals. Such minerals are found inside the porous rock matrix and as infill particles or coatings in fractures in the rock. The model formulation also allows for different flow modes such as flow along discrete fractures in sparsely fractured rocks and along flowpaths in a fracture network. The scavenging of oxygen is exemplified for these cases as well as for more comprehensive applications, including glaciation considerations. Results show that chemical reaction kinetics control the scavenging of oxygen during a relatively short time with respect to the lifetime of the repository. For longer times the scavenging of oxygen is controlled by transport processes in the porous rock matrix. The penetration depth of oxygen along the flowpath depends largely on the hydraulic properties, which may vary significantly between different locations and situations. The results indicate that oxygen, in the absence of easily degradable organic matter, may reach long distances along a flow path during the life-time of the repository (hundreds to thousands of metres in a million years depending on e.g. hydraulic properties of the flow path and the availability of reducing capacity). However, large uncertainties regarding key input parameters exist leading to the conclusion that the results from the model must be treated with caution pending more accurate and validated data. Ongoing and planned experiments are expected to reduce these uncertainties, which are required in order to make more reliable predictions for a safety assessment of a nuclear waste repository.     

Fast method for simulation of radionuclide chain migration in dual porosity fracture rocks

In fractured rocks with a porous rock matrix such as granites, radionuclides will flow with the water in the fracture network. The nuclides will diffuse in and out the rock matrix where they can sorb and be considerably retarded compared to the water velocity. A water parcel entering the network will mix and split at the fracture intersections and parts of the original parcel will traverse a multitude of different fractures. The flowrates, velocities, sizes and apertures of the fractures can vary widely. Normally one must solve the transport equations for every fracture and use the effluent concentration as inlet condition to the next fracture and so on. It is shown that under some weakly simplified conditions it suffices to determine one single parameter group containing information on the flow wetted surface that a water parcel contacts along the entire path. It is also shown how this can be obtained. Then, solving the transport equations only once for time and location along the path gives the concentration and nuclide flux of every nuclide in the chain everywhere along a path. The same solution actually is valid for every path in the network. This dramatically reduces the computation effort. The same approach can be used for models based on streamtubes.     

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.     

Solute transport in crystalline rocks at Aspö--I: geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted Asp? Hard Rock Laboratory (SE Sweden). On a scale of decametres and larger, steeply dipping faults dominate and contain a variety of different fault rocks (mylonites, cataclasites, fault gouges). On a smaller scale, somewhat less regular fracture patterns were found. Conceptual models of the fault and fracture geometries and of the properties of rock types adjacent to fractures were derived and used as input for the modelling of in situ dipole tracer tests that were conducted in the framework of the Tracer Retention Understanding Experiment (TRUE-1) on a scale of metres. After the identification of all relevant transport and retardation processes, blind predictions of the breakthroughs of conservative to moderately sorbing tracers were calculated and then compared with the experimental data. This paper provides the geological basis and model calibration, while the predictive and inverse modelling work is the topic of the companion paper [J. Contam. Hydrol. 61 (2003) 175]. The TRUE-1 experimental volume is highly fractured and contains the same types of fault rocks and alterations as on the decametric scale. The experimental flow field was modelled on the basis of a 2D-streamtube formalism with an underlying homogeneous and isotropic transmissivity field. Tracer transport was modelled using the dual porosity medium approach, which is linked to the flow model by the flow porosity. Given the substantial pumping rates in the extraction borehole, the transport domain has a maximum width of a few centimetres only. It is concluded that both the uncertainty with regard to the length of individual fractures and the detailed geometry of the network along the flowpath between injection and extraction boreholes are not critical because flow is largely one-dimensional, whether through a single fracture or a network. Process identification and model calibration were based on a single uranine breakthrough (test PDT3), which clearly showed that matrix diffusion had to be included in the model even over the short experimental time scales, evidenced by a characteristic shape of the trailing edge of the breakthrough curve. Using the geological information and therefore considering limited matrix diffusion into a thin fault gouge horizon resulted in a good fit to the experiment. On the other hand, fresh granite was found not to interact noticeably with the tracers over the time scales of the experiments. While fracture-filling gouge materials are very efficient in retarding tracers over short periods of time (hours-days), their volume is very small and, with time progressing, retardation will be dominated by altered wall rock and, finally, by fresh granite. In such rocks, both porosity (and therefore the effective diffusion coefficient) and sorption K(d)s are more than one order of magnitude smaller compared to fault gouge, thus indicating that long-term retardation is expected to occur but to be less pronounced.     

Analysis of the migration of instantaneously injected cesium in artificial fractures of Lac du Bonnet granite, Manitoba, Canada

Breakthrough curves of ¹³⁷Cs and tritiated water injected instantaneously into artificial fractures in Lac du Bonnet granite were analyzed using the analytical solution for a single rock-fracture system and assuming the linear sorption isotherm of the solute. Parameters of nuclide diffusion and sorption in rock matrices, obtained by fitting, varied depending on the flow velocity in the fractures. According to theoretical calculations, different fracture flow velocities lead to different diffusion distances of nuclides in matrices at the same injection volume. As microscopic inhomogeneity is considered to exist in the rock matrix, the average diffusion-sorption characteristics of the matrix within the diffusion distance may have varied with the fracture flow velocity. Surface sorption was marked in fractures that had relatively high matrix sorption-diffusion capacities. The phenomenon was interpreted using the theoretical relationships developed between the surface sorption, matrix sorption and pore diffusion coefficient, and the porosity of matrices.The effect of the nonlinear sorption of solute was examined by numerically solving model equations that incorporate the nonlinear isotherm. This incorporation may contribute to the reduction of deviations between theoretical and experimental BTC's.     

Radionuclide transport and retardation in rock fracture and crushed rock column experiments

Transport and retardation of non-sorbing tritiated water and chloride and slightly sorbing sodium was studied in Syyry area SY-KR7 mica gneiss, in altered porous tonalite and in fresh tonalite. Experiments were performed using dynamic fracture and crushed rock column methods. Static batch method for sodium was introduced to compare retardation values from static and dynamic experiments. The ¹⁴C-PMMA method was used to study the pore structure of matrices. The pore aperture distribution was evaluated from Hg-porosimetry determinations and the surface areas were determined using the B.E.T. method. The flow characteristics and transport behavior of tracers were interpreted using a numerical compartment model for dispersion. The effect of matrix diffusion was calculated using an analytical solution to the advection-matrix diffusion problem in which surface retardation was taken into account. Radionuclide transport behavior in rock fractures was explained on the basis of rock structure.     

Characterization of fracture aperture field heterogeneity by electrical resistance measurement

We use electrical resistance measurements to characterize the aperture field in a rough fracture. This is done by performing displacement experiments using two miscible fluids of different electrical resistivity and monitoring the time variation of the overall fracture resistance. Two fractures have been used: their complementary rough walls are identical but have different relative shear displacements which create "channel" or "barrier" structures in the aperture field, respectively parallel or perpendicular to the mean flow velocity U(→). In the "channel" geometry, the resistance displays an initial linear variation followed by a tail part which reflects the velocity contrast between slow and fast flow channels. In the "barrier" geometry, a change in the slope between two linear zones suggests the existence of domains of different characteristic aperture along the fracture. These variations are well reproduced analytically and numerically using simple flow models. For each geometry, we present then a data inversion procedure that allows one to extract the key features of the heterogeneity from the resistance measurement.     

Applied tracer tests in fractured rock: Can we predict natural gradient solute transport more accurately than fracture and matrix parameters?

Applied tracer tests provide a means to estimate aquifer parameters in fractured rock. The traditional approach to analysing these tests has been using a single fracture model to find the parameter values that generate the best fit to the measured breakthrough curve. In many cases, the ultimate aim is to predict solute transport under the natural gradient. Usually, no confidence limits are placed on parameter values and the impact of parameter errors on predictions of solute transport is not discussed. The assumption inherent in this approach is that the parameters determined under forced conditions will enable prediction of solute transport under the natural gradient. This paper considers the parameter and prediction uncertainty that might arise from analysis of breakthrough curves obtained from forced gradient applied tracer tests. By adding noise to an exact solution for transport in a single fracture in a porous matrix we create multiple realisations of an initial breakthrough curve. A least squares fitting routine is used to obtain a fit to each realisation, yielding a range of parameter values rather than a single set of absolute values. The suite of parameters is then used to make predictions of solute transport under lower hydraulic gradients and the uncertainty of estimated parameters and subsequent predictions of solute transport is compared. The results of this study show that predictions of breakthrough curve characteristics (first inflection point time, peak arrival time and peak concentration) for groundwater flow speeds with orders of magnitude smaller than that at which a test is conducted can sometimes be determined even more accurately than the fracture and matrix parameters.     

Identification of key parameters controlling dissolved oxygen migration and attenuation in fractured crystalline rocks

In the crystalline rocks of the Canadian Shield, geochemical conditions are currently reducing at depths of 500-1000 m. However, during future glacial periods, altered hydrologic conditions could potentially result in enhanced recharge of glacial melt water containing a relatively high concentration of dissolved oxygen (O2). It is therefore of interest to investigate the physical and geochemical processes, including naturally-occurring redox reactions, that may control O2 ingress. In this study, the reactive transport code MIN3P is used in combination with 2k factorial analyses to identify the most important parameters controlling oxygen migration and attenuation in fractured crystalline rocks. Scenarios considered are based on simplified conceptual models that include a single vertical fracture, or a fracture zone, contained within a rock matrix that extends from the ground surface to a depth of 500 m. Consistent with field observations, Fe(II)-bearing minerals are present in the fractures (i.e. chlorite) and the rock matrix (biotite and small quantities of pyrite). For the parameter ranges investigated, results indicate that for the single fracture case, the most influential factors controlling dissolved O2 ingress are flow velocity in the fracture, fracture aperture, and the biotite reaction rate in the rock matrix. The most important parameters for the fracture zone simulations are flow velocity in the individual fractures, pO2 in the recharge water, biotite reaction rate, and to a lesser degree the abundance and reactivity of chlorite in the fracture zone, and the fracture zone width. These parameters should therefore receive increased consideration during site characterization, and in the formulation of site-specific models intended to predict O2 behavior in crystalline rocks.