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.     

Field tracer-transport tests in unsaturated fractured tuff.

This paper presents the results of a field investigation in the unsaturated, fractured welded tuff within the Exploratory Studies Facility (ESF) at Yucca Mountain, NV. This investigation included a series of tests during which tracer-laced water was released into a high-permeability zone within a horizontal injection borehole. The tracer concentration was monitored in the seepage collected in an excavated slot about 1.6 m below the borehole. Results showed significant variability in the hydrologic response of fractures and the matrix. Analyses of the breakthrough curves suggest that flow and transport pathways are dynamic, rather than fixed, and related to liquid-release rates. Under high release rates, fractures acted as the predominant flow pathways, with limited fracture-matrix interaction. Under low release rates, fracture flow was comparatively less dominant, with a noticeable contribution from matrix flow. Observations of tracer concentrations rebounding in seepage water, following an interruption of flow, provided evidence of mass exchange between the fast-flowing fractures and slow- or non-flowing regions. The tests also showed the applicability of fluorinated benzoate tracers in situations where multiple tracers of similar physical properties are warranted.     

Tracer test at El Berrocal site

Tracer tests provide highly valuable information about the transport properties of saturated rocks which is essential to the characterization of a potential radioactive waste repository site. In the frame of El Berrocal project, a set of tracer tests was performed in a complex geometry of inclined boreholes, combined with highly fractured transmissive zones. The aims of the tracer test programme were to gain experience, knowledge and insight into field transport experiments. To achieve this a detailed programme of new instrumentation design, site characterization and laboratory tasks was developed. For field monitoring a new electronic system was developed. The system is able to measure up to 256 parameters per borehole, with surface equipment to control pumping rates and physical and chemical parameters at both injection and extraction boreholes. The experiments progressed from single borehole dilution tests under both natural flow and forced gradient conditions to convergent flow tracer tests. Dilution tests helped to discriminate the most suitable borehole sections at which to inject tracers. The tracers were selected by the results of the laboratory programme. Uranine (fluorescein), DTPA-gadolinium (diethylenetriaminopentacetic acid, gadolinium (III)), and deuterium were injected simultaneously in one borehole section and recovered at another borehole 20 m away, pumping at a flowrate of 0.1 1 min⁻¹. First results showed a thickness porosity of 1.2 × 10⁻³ m and a longitudinal dispersivity of 17.0 m using uranine data acquired over a period of 4 d, at which point the recovery concentration had reached a maximum. However, gadolinium and deuterium appeared to travel faster, arriving at peak values after only 2 d of injection.     

Using constant head step tests to determine hydraulic apertures in fractured rock

The initial step in the analysis of contaminant transport in fractured rock requires the consideration of groundwater velocity. Practical methods for estimating the average linear groundwater velocity (vˉ) in fractured rock require determination of hydraulic apertures which are commonly calculated by applying the cubic law using transmissivity (T) values and the number of hydraulically active fractures in the test interval. High-resolution, constant-head step injection testing of cored boreholes in a 100 m thick fractured dolostone aquifer was conducted using inflatable packers to isolate specific test intervals from the rest of the borehole. The steps in each test interval were gradually increased from very low to much higher injection rates. At smaller injection rates, the flow rate vs. applied pressure graph projects through the origin and indicates Darcian flow; non Darcian flow is evident at higher injection rates. Non-Darcian flow results in significantly lower calculated T values, which translates to smaller hydraulic aperture values. Further error in the calculated hydraulic aperture stems from uncertainty in the number of hydraulically active fractures in each test interval. This estimate can be inferred from borehole image and core logs, however, all of the fractures identified are not necessarily hydraulically active. This study proposes a method based on Reynolds number calculations aimed at improving confidence in the selection of the number of active fractures in each test interval.     

Fluorescent dye imaging of the volume sampled by single well forced-gradient tracer tests evaluated in a laboratory-scale aquifer physical model

This study presents a new method to visualise forced-gradient tracer tests in 2-D using a laboratory-scale aquifer physical model. Experiments were designed to investigate the volume of aquifer sampled in vertical dipole flow tracer tests (DFTT) and push-pull tests (PPT), using a miniature monitoring well and straddle packer arrangement equipped with solute injection and recovery chambers. These tests have previously been used to estimate bulk aquifer hydraulic and transport properties for the evaluation of natural attenuation and other remediation approaches. Experiments were performed in a silica glass bead-filled box, using a fluorescent tracer (fluorescein) to deduce conservative solute transport paths. Digital images of fluorescein transport were captured under ultraviolet light and processed to analyse tracer plume geometry and obtain point-concentration breakthrough histories. Inorganic anion mixtures were also used to obtain conventional tracer breakthrough histories. Concentration data from the conservative tracer breakthrough curves was compared with the digital images and a well characterised numerical model. The results show that the peak tracer breakthrough response in dipole flow tracer tests samples a zone of aquifer close to the well screen, while the sampling volume of push-pull tests is limited by the length of the straddle packers used. The effective sampling volume of these single well forced-gradient tests in isotropic conditions can be estimated with simple equations. The experimental approach offers the opportunity to evaluate under controlled conditions the theoretical basis, design and performance of DFTTs and PPTs in porous media in relation to measured flow and transport properties.     

Solute transport in crystalline rocks at Äspö — I: Geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted Äspö 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_ds 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.     

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.     

Anthropogenic contaminants as tracers in an urbanizing karst aquifer

Karst aquifers are uniquely vulnerable to contamination. In the Barton Springs segment of the karstic Edwards aquifer (Texas, U.S.A.), urban contaminants such as pesticides and volatile organic compounds frequently are detected in spring base flow. To determine whether contaminant concentrations change in response to storms, and if they therefore might act as tracers of focused recharge, samples were collected from Barton Springs at closely spaced intervals following three storms. Two herbicides (atrazine and simazine), two insecticides (carbaryl and diazinon), and a solvent (tetrachloroethene) described breakthrough curves over a 1-week period following one or more storms. The breakthrough curves were decomposed into two to five log-normal subcurves, which were interpreted as representing pulses of contaminants moving through the aquifer. Each subcurve could be used in the same way as an artificial tracer to determine travel time to and recovery at the spring. The contaminants have several advantages over artificial tracers: they represent the actual compounds of interest, they are injected essentially simultaneously at several points, and they are injected under those conditions when transport is of the most interest, i.e., following storms. The response of storm discharge, specific conductance, and contaminant loading at the spring depended on initial aquifer flow conditions, which varied from very low (spring discharge of 0.48 m3/s) to high (spring discharge of 2.7 m3/s): concentrations and recovery were the highest when initial aquifer flow conditions were low. This behavior provides information about aquifer structure and the influence of aquifer flow condition on transport properties.     

Prediction of contaminant transport in fractured carbonate aquifer types: a case study of the Permian Magnesian Limestone Group (NE England,UK)

Viruses and bacteria which are characterized by finite lives in the subsurface are rapidly transported via fractures and cavities in fractured and karst aquifers. Here, we demonstrate how the coupling of a robust outcrop characterization and hydrogeophysical borehole testing is essential for prediction of contaminant velocities and hence wellhead protection areas. To show this, we use the dolostones of the Permian Magnesian Limestone aquifer in NE England, where we incorporated such information in a groundwater flow and particle tracking model. Within this aquifer, flow in relatively narrow (mechanical aperture of ~?10^?1–1 mm) fractures is coupled with that in pipe cavities (~?0.20-m diameter) following normal faults. Karstic cavities and narrow fractures are hydraulically very different. Thus, the solutional features are represented within the model by a pipe network (which accounts for turbulence) embedded within an equivalent porous medium representing Darcian flowing fractures. Incorporation of fault conduits in a groundwater model shows that they strongly influence particle tracking results. Despite this, away from faulted areas, the effective flow porosity of the equivalent porous medium remains a crucial parameter. Here, we recommend as most appropriate a relatively low value of effective porosity (of 2.8?×?10^?4) based on borehole hydrogeophysical testing. This contrasts with earlier studies using particle tracking analyses on analogous carbonate aquifers, which used much higher values of effective porosity, typically ~?10² times higher than our value, resulting in highly non-conservative estimates of aquifer vulnerability. Low values of effective flow porosities yield modelled flow velocities ranging from ~?100 up to ~?500 m/day in un-faulted areas. However, the high fracturing density and presence of karstic cavities yield modelled flow velocities up to ~?9000 m/day in fault zones. The combination of such flow velocities along particle traces results in 400-day particle traces up to 8-km length, implying the need for large well protection areas and high aquifer vulnerability to slowly degrading contaminants.

     

Movement of leachate from beneath turkey litter sited over chalk in southern England

Farm waste stores are widespread in the UK, with many overlying the principal aquifer, the Chalk. The stores pose a threat to groundwater quality through the infiltration of high concentrations of nitrogen species and organic carbon together with pathogenic microbes. Two cored boreholes have been drilled into the unsaturated chalk to depths of 15 and 20 metres, respectively, through a site which has been used to store turkey litter for in excess of 20 years. Porewaters were extracted from the cores and analysed for a range of chemical elements. In addition, chalk core material was also taken for microbial examination. Both boreholes showed very high concentrations of nitrate-N (3000 mg/L), ammonia (5000 mg/L), organic carbon (3000 mg/L), and potassium (10,000 mg/L) in the top 5 metres of the profile. Below this depth concentrations declined dramatically. Highest concentrations were found in the borehole constructed in the middle of the site. The borehole constructed at the edge of the store showed much lower concentrations but did show a peak of nitrate around 10 metres below ground level. The apparent lack of movement beneath the centre of the store suggests the turkey litter is relatively impermeable and most leaching occurs where the covering of litter is thin or when the litter is annually cleared. If the leachate continues to migrate at this apparent rate, it will take more than 100 years to reach the water table.     

Diffusive partitioning tracer test for the quantification of nonaqueous phase liquid (NAPL) in the vadose zone: Performance evaluation for heterogeneous NAPL distribution

A partitioning tracer test based on gas-phase diffusion in the vadose zone yields estimates of the residual nonaqueous phase liquid (NAPL) saturation. The present paper investigates this technique further by studying diffusive tracer breakthrough curves in the vadose zone for a heterogeneous NAPL distribution. Tracer experiments were performed in a lysimeter with a horizontal layer of artificial kerosene embedded in unsaturated sand. Tracer disappearance curves at the injection point and tracer breakthrough curves at some distance from the injection point were measured inside and outside of the NAPL layer. A numerical code was used to generate independent model predictions based on the physicochemical sand, NAPL, and tracer properties. The measured and modeled tracer breakthrough curves were in good agreement confirming the validity of important modeling assumptions such as negligible sorption of chlorofluorocarbon (CFC) tracers to the uncontaminated sand and their fast reversible partitioning between the soil air and the NAPL phase. Subsequently, the model was used to investigate different configurations of NAPL contamination. The experimental and model results show that the tracer disappearance curves of a single-well diffusive partitioning tracer test (DPTT) are dominated by the near-field presence of NAPL around the tip of the soil gas probe. In contrast, breakthrough curves of inter-well tracer tests reflect the NAPL saturation in between the probes, although there is no unique interpretation of the tracer signals if the NAPL distribution is heterogeneous. Numerical modeling is useful for the planning of a DPTT application. Simulations suggest that several cubic meters of soil can be investigated with a single inter-well partitioning tracer test of 24-hour duration by placing the injection point in the center of the investigated soil volume and probes at up to 1 m distance for the monitoring of gaseous tracers.     

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.     

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.     

Nutrient dynamics as indicators of karst processes: comparison of the Chalk aquifer (Normandy, France) and the Edwards aquifer (Texas, U.S.A.)

Karst aquifers display a range of geologic and geomorphic characteristics in a wide range of climatic and land-use settings; identification of transport dynamics representative of karst aquifers in general could help advance our understanding of these complex systems. To this end, nutrient, turbidity, and major ion dynamics in response to storms were compared at multiple sites in two karst aquifers with contrasting characteristics and settings: the Chalk aquifer (Eure Department, Normandy, France) and the Barton Springs segment of the Edwards Aquifer (Texas, U.S.A.). The Chalk aquifer is typified by high matrix porosity, thick surficial deposits (up to 30 m thick), and agricultural land use; the Barton Springs segment is typified by low matrix porosity, outcropping limestone, and urban land use. Following one to three storms, from 5 to 16 samples from springs and wells were analyzed for major ions, and specific conductance and turbidity were monitored continuously. Comparison of the chemographs indicated some generalized responses, including an increase in turbidity and potassium concentrations and a decrease in major ion and nitrate concentrations with infiltrating storm runoff. Factor analysis of major ions and turbidity revealed strikingly similar behavior of the chemical variables for the two aquifers: The first two factors, explaining more than 75% of the variability, illustrate that dynamics of most major ions (including nitrate) are opposed to those of turbidity and of potassium. The results demonstrate that potassium and nitrate are effective tracers of infiltrating storm runoff and resident ground water, respectively, and the similar results for these two highly contrasting aquifers suggest that the dynamics identified might be applicable to karst systems in general.     

Nitrate-enhanced bioremediation of BTEX-contaminated groundwater: parameter estimation from natural-gradient tracer experiments

Two natural-gradient pulse tracer tests were conducted in a petroleum-contaminated aquifer to evaluate the potential for benzene, toluene, ethylbenzene, and xylenes (BTEX) biodegradation under enhanced nitrate-reducing conditions. Addition of nitrate resulted in loss of toluene, ethylbenzene, and m,p-xylenes (TEX) after an initial lag period of approximately 9 days. Losses of benzene were not observed over the 60-day monitoring period. Tracer breakthrough curves (BTCs) were analyzed to derive transport and biodegradation parameters, including advective velocities, retardation factors, dispersion coefficients, biodegradation rate constants, and nitrate utilization ratios. Using the parameters derived from the BTC analysis, numerical simulations of one of the tracer experiments were conducted using BIONAPL/3D [Molson, J., BIONAPL/3D User Guide, A 3D Coupled Flow and Multi-Component Reactive transport model. University of Waterloo, Waterloo, Ontario, Canada]. Simulations using the BTC-derived transport and biodegradation parameters successfully reproduced benzene, TEX, and nitrate concentrations measured during the tracer experiment. Comparisons of observed and simulated nitrate concentrations indicate that the mass ratio of nitrate-N utilized to TEX degraded increased over time during the experiment, reaching values many times that expected based on stoichiometry of TEX oxidation coupled to nitrate reduction. Excess nitrate loss is likely due to oxidation of other organics in addition to TEX.     

Experimental determination of sorption in fractured flow systems

Fracture "skins" are alteration zones on fracture surfaces created by a variety of biological, chemical, and physical processes. Skins increase surface area, where sorption occurs, compared to the unaltered rock matrix. This study examines the sorption of organic solutes on altered fracture surfaces in an experimental fracture-flow apparatus. Fracture skins containing abundant metal oxides, clays, and organic material from the Breathitt Formation (Kentucky, USA) were collected in a manner such that skin surface integrity was maintained. The samples were reassembled in the lab in a flow-through apparatus that simulated approximately 2.7 m of a linear fracture "conduit." A dual-tracer injection scheme was utilized with the sorbing or reactive tracer compared to a non-reactive tracer (chloride) injected simultaneously. Sorption was assessed from the ratio of the first temporal moments of the breakthrough curves and from the loss of reactive tracer mass and evaluated as a function of flow velocity and solute type. The breakthrough curves suggest dual-flow regimes in the fracture with both sorbing and non-sorbing flow fields. Significant sorption occurs for the reactive components, and sorption increased with decreasing flow rate and decreasing compound solubility. Based on moment analysis, however, there was little retardation of the center of solute mass. These data suggest that non-equilibrium sorption processes dominate and that slow desorption and boundary layer diffusion cause extensive tailing in the breakthrough curves.     

A method for determining volatile organic solvents in chalk pore waters (Southern and Eastern England) and its relevance to the evaluation of groundwater contamination

Conventional methods of sampling groundwater in Chalk aquifers sample the mobile water in the fissures only and not the non-mobile water held in the microporous matrix. However, the microporous matrix can play an important role in modifying groundwater quality by attenuating contaminant migration through diffusion exchange and, in the long term, by acting as a reservoir for contaminant retention.A method for extracting and analysing chalk pore waters for chlorinated solvents has been developed. Chalk pore water-depth profiles at the site of a spillage of tetrachloroethene demonstrated that the immiscible phase of the solvent had migrated to depths of 50 m, producing solvent concentrations in the matrix pore waters of up to 40,000 μg L⁻¹.In addition, the pore water-depth profile indicated that the immiscible phase had ponded upon a less permeable horizon within the Chalk aquifer at shallow depths.     

Combined use of field and laboratory testing to predict preferred flow paths in an heterogeneous aquifer

Elevated nitrate concentrations within a municipal water supply aquifer led to pilot testing of a field-scale, in situ denitrification technology based on carbon substrate injections. In advance of the pilot test, detailed characterization of the site was undertaken. The aquifer consisted of complex, discontinuous and interstratified silt, sand and gravel units, similar to other well studied aquifers of glaciofluvial origin, 15-40 m deep. Laboratory and field tests, including a conservative tracer test, a pumping test, a borehole flowmeter test, grain-size analysis of drill cuttings and core material, and permeameter testing performed on core samples, were performed on the most productive depth range (27-40 m), and the results were compared. The velocity profiles derived from the tracer tests served as the basis for comparison with other methods. The spatial variation in K, based on grain-size analysis, using the Hazen method, were poorly correlated with the breakthrough data. Trends in relative hydraulic conductivity (K/K(avg)) from permeameter testing compared somewhat better. However, the trends in transient drawdown with depth, measured in multilevel sampling points, corresponded particularly well with those of solute mass flux. Estimates of absolute K, based on standard pumping test analysis of the multilevel drawdown data, were inversely correlated with the tracer test data. The inverse nature of the correlation was attributed to assumptions in the transient drawdown packages that were inconsistent with the variable diffusivities encountered at the scale of the measurements. Collectively, the data showed that despite a relatively low variability in K within the aquifer under study (within a factor of 3), water and solute mass fluxes were concentrated in discrete intervals that could be targeted for later bioremediation.