Inferring transport characteristics in a fractured rock aquifer by combining single-hole ground-penetrating radar reflection monitoring and tracer test data

Investigations of solute transport in fractured rock aquifers often rely on tracer test data acquired at a limited number of observation points. Such data do not, by themselves, allow detailed assessments of the spreading of the injected tracer plume. To better understand the transport behavior in a granitic aquifer, we combine tracer test data with single-hole ground-penetrating radar (GPR) reflection monitoring data. Five successful tracer tests were performed under various experimental conditions between two boreholes 6 m apart. For each experiment, saline tracer was injected into a previously identified packed-off transmissive fracture while repeatedly acquiring single-hole GPR reflection profiles together with electrical conductivity logs in the pumping borehole. By analyzing depth-migrated GPR difference images together with tracer breakthrough curves and associated simplified flow and transport modeling, we estimate (1) the number, the connectivity, and the geometry of fractures that contribute to tracer transport, (2) the velocity and the mass of tracer that was carried along each flow path, and (3) the effective transport parameters of the identified flow paths. We find a qualitative agreement when comparing the time evolution of GPR reflectivity strengths at strategic locations in the formation with those arising from simulated transport. The discrepancies are on the same order as those between observed and simulated breakthrough curves at the outflow locations. The rather subtle and repeatable GPR signals provide useful and complementary information to tracer test data acquired at the outflow locations and may help us to characterize transport phenomena in fractured rock aquifers.

Calibration of high-resolution geophysical data with tracer test measurements to improve hydrological predictions

Relationships between porosity and hydraulic conductivity tend to be strongly scale- and site-dependent and are thus very difficult to establish. As a result, hydraulic conductivity distributions inferred from geophysically derived porosity models must be calibrated using some measurement of aquifer response. This type of calibration is potentially very valuable as it may allow for transport predictions within the considered hydrological unit at locations where only geophysical measurements are available, thus reducing the number of well tests required and thereby the costs of management and remediation. Here, we explore this concept through a series of numerical experiments. Considering the case of porosity characterization in saturated heterogeneous aquifers using crosshole ground-penetrating radar and borehole porosity log data, we use tracer test measurements to calibrate a relationship between porosity and hydraulic conductivity that allows the best prediction of the observed hydrological behavior. To examine the validity and effectiveness of the obtained relationship, we examine its performance at alternate locations not used in the calibration procedure. Our results indicate that this methodology allows us to obtain remarkably reliable hydrological predictions throughout the considered hydrological unit based on the geophysical data only. This was also found to be the case when significant uncertainty was considered in the underlying relationship between porosity and hydraulic conductivity.

Stochastic inversion of tracer test and electrical geophysical data to estimate hydraulic conductivities.

Quantifying the spatial configuration of hydraulic conductivity (K) in heterogeneous geological environments is essential for accurate predictions of contaminant transport, but is difficult because of the inherent limitations in resolution and coverage associated with traditional hydrological measurements. To address this issue, we consider crosshole and surface-based electrical resistivity geophysical measurements, collected in time during a saline tracer experiment. We use a Bayesian Markov-chain-Monte-Carlo (McMC) methodology to jointly invert the dynamic resistivity data, together with borehole tracer concentration data, to generate multiple posterior realizations of K that are consistent with all available information. We do this within a coupled inversion framework, whereby the geophysical and hydrological forward models are linked through an uncertain relationship between electrical resistivity and concentration. To minimize computational expense, a facies-based subsurface parameterization is developed. The Bayesian-McMC methodology allows us to explore the potential benefits of including the geophysical data into the inverse problem by examining their effect on our ability to identify fast flowpaths in the subsurface, and their impact on hydrological prediction uncertainty. Using a complex, geostatistically generated, two-dimensional numerical example representative of a fluvial environment, we demonstrate that flow model calibration is improved and prediction error is decreased when the electrical resistivity data are included. The worth of the geophysical data is found to be greatest for long spatial correlation lengths of subsurface heterogeneity with respect to wellbore separation, where flow and transport are largely controlled by highly connected flowpaths.

Dispersivity Determination Through a Modeling Approach From a Tracer Test Based on Total Br Concentration in Soil Samples

Determination of reliable solute transport parameters is an essential aspect for the characterization of the mechanisms and processes involved in solute transport (e.g., pesticides, fertilizers, contaminants) through the unsaturated zone. A rapid inexpensive method to estimate the dispersivity parameter at the field scale is presented herein. It is based on the quantification by the X-ray fluorescence solid-state technique of total bromine in soil, along with an inverse numerical modeling approach. The results show that this methodology is a good alternative to the classic Br− determination in soil water by ion chromatography. A good agreement between the observed and simulated total soil Br is reported. The results highlight the potential applicability of both combined techniques to infer readily solute transport parameters under field conditions.

Solute transport in crystalline rocks at Äspö — II: Blind predictions, inverse modelling and lessons learnt from test STT1

Based on the results from detailed structural and petrological characterisation and on up-scaled laboratory values for sorption and diffusion, blind predictions were made for the STT1 dipole tracer test performed in the Swedish A¨ spo¨ Hard Rock Laboratory. The tracers used were nonsorbing, such as uranine and tritiated water, weakly sorbing 22Na+, 85Sr2 +, 47Ca2 +and more strongly sorbing 86Rb+, 133Ba2 +, 137Cs+. Our model consists of two parts: (1) a flow part based on a 2D-streamtube formalism accounting for the natural background flow field and with an underlying homogeneous and isotropic transmissivity field and (2) a transport part in terms of the dual porosity medium approach which is linked to the flow part by the flow porosity. The calibration of the model was done using the data from one single uranine breakthrough (PDT3). The study clearly showed that matrix diffusion into a highly porous material, fault gouge, had to be included in our model evidenced by the characteristic shape of the breakthrough curve and in line with geological observations. After the disclosure of the measurements, it turned out that, in spite of the simplicity of our model, the prediction for the nonsorbing and weakly sorbing tracers was fairly good. The blind prediction for the more strongly sorbing tracers was in general less accurate. The reason for the good predictions is deemed to be the result of the choice of a model structure strongly based on geological observation. The breakthrough curves were inversely modelled to determine in situ values for the transport parameters and to draw consequences on the model structure applied. For good fits, only one additional fracture family in contact with cataclasite had to be taken into account, but no new transport mechanisms had to be invoked. The in situ values for the effective diffusion coefficient for fault gouge are a factor of 2–15 larger than the laboratory data. For cataclasite, both data sets have values comparable to laboratory data. The extracted Kd values for the weakly sorbing tracers are larger than Swedish laboratory data by a factor of 25–60, but agree within a factor of 3–5 for the more strongly sorbing nuclides. The reason for the inconsistency concerning Kds is the use of fresh granite in the laboratory studies, whereas tracers in the field experiments interact only with fracture fault gouge and to a lesser extent with cataclasite both being mineralogically very different (e.g. clay-bearing) from the intact wall rock.

Calibration of tracer concentration measurement and results

A reliable assessment of relevant substance flows is very important for environmental risk assessments and efficiency analysis of measures to reduce or avoid emissions of micropollutants like drugs to water systems. Accordingly, a detailed preparation of monitoring campaigns should include an accuracy check for the sampling configuration to prove the reliability of the monitoring results and the subsequent data processing. The accuracy of substance flow analyses is expected to be particularly weak for substances having high short-term variations of concentrations in sewage. This is especially the case linked to the observation of substance flows close to source in waste water systems. The verification of a monitoring configuration in a hospital sewer in Luxembourg is in the centre of interest of the case study presented here. A tracer test in the sewer system under observation is an essential element of the suggested accuracy check and provides valuable information for an uncertainty analysis. The results illustrate the importance of accuracy checks as an essential element of the preparation of monitoring campaigns. Moreover the study shows that continuous flow proportional sampling enables a representative observation of short-term peak loads of the iodinated x-ray contrast media iobitridol close to source.

Integration of high-resolution geophysical data for the characterization of alluvial aquifers : new approaches and implications for hydrological predictions

Abstract Accurate characterization of the spatial distribution of hydrological properties in heterogeneous aquifers at a range of scales is a key prerequisite for reliable modeling of subsurface contaminant transport, and is essential for designing effective and cost-efficient groundwater management and remediation strategies. To this end, high-resolution geophysical methods have shown significant potential to bridge a critical gap in subsurface resolution and coverage between traditional hydrological measurement techniques such as borehole log/core analyses and tracer or pumping tests. An important and still largely unresolved issue, however, is how to best quantitatively integrate geophysical data into a characterization study in order to estimate the spatial distribution of one or more pertinent hydrological parameters, thus improving hydrological predictions. Recognizing the importance of this issue, the aim of the research presented in this thesis was to first develop a strategy for the assimilation of several types of hydrogeophysical data having varying degrees of resolution, subsurface coverage, and sensitivity to the hydrologic parameter of interest. In this regard a novel simulated annealing (SA)-based conditional simulation approach was developed and then tested in its ability to generate realizations of porosity given crosshole ground-penetrating radar (GPR) and neutron porosity log data. This was done successfully for both synthetic and field data sets. A subsequent issue that needed to be addressed involved assessing the potential benefits and implications of the resulting porosity realizations in terms of groundwater flow and contaminant transport. This was investigated synthetically assuming first that the relationship between porosity and hydraulic conductivity was well-defined. Then, the relationship was itself investigated in the context of a calibration procedure using hypothetical tracer test data. Essentially, the relationship best predicting the observed tracer test measurements was determined given the geophysically derived porosity structure. Both of these investigations showed that the SA-based approach, in general, allows much more reliable hydrological predictions than other more elementary techniques considered. Further, the developed calibration procedure was seen to be very effective, even at the scale of tomographic resolution, for predictions of transport. This also held true at locations within the aquifer where only geophysical data were available. This is significant because the acquisition of hydrological tracer test measurements is clearly more complicated and expensive than the acquisition of geophysical measurements. Although the above methodologies were tested using porosity logs and GPR data, the findings are expected to remain valid for a large number of pertinent combinations of geophysical and borehole log data of comparable resolution and sensitivity to the hydrological target parameter. Moreover, the obtained results allow us to have confidence for future developments in integration methodologies for geophysical and hydrological data to improve the 3-D estimation of hydrological properties.

Single-hole GPR reflection imaging of solute transport in a granitic aquifer

Identifying transport pathways in fractured rock is extremely challenging as flow is often organized in a few fractures that occupy a very small portion of the rock volume. We demonstrate that saline tracer experiments combined with single-hole ground penetrating radar (GPR) reflection imaging can be used to monitor saline tracer movement within mm-aperture fractures. A dipole tracer test was performed in a granitic aquifer by injecting a saline solution in a known fracture, while repeatedly acquiring single-hole GPR sections in the pumping borehole located 6 m away. The final depth-migrated difference sections make it possible to identify consistent temporal changes over a 30 m depth interval at locations corresponding to fractures previously imaged in GPR sections acquired under natural flow and tracer-free conditions. The experiment allows determining the dominant flow paths of the injected tracer and the velocity (0.4-0.7 m/min) of the tracer front. Citation: Dorn, C., N. Linde, T. Le Borgne, O. Bour, and L. Baron (2011), Single-hole GPR reflection imaging of solute transport in a granitic aquifer, Geophys. Res. Lett., 38, L08401, doi: 10.1029/2011GL047152.

Deposition and residues of azoxystrobin and imidacloprid on greenhouse lettuce with implications for human consumption

Lettuce greenhouse experiments were carried out from March to June 2011 in order to analyze how pesticides behave from the time of application until their intake via human consumption taking into account the primary distribution of pesticides, field dissipation, and post-harvest processing. In addition, experimental conditions were used to evaluate a new dynamic plant uptake model comparing its results with the experimentally derived residues. One application of imidacloprid and two of azoxystrobin were conducted. For evaluating primary pesticide distribution, two approaches based on leaf area index and vegetation cover were used and results were compared with those obtained from a tracer test. High influence of lettuce density, growth stage and type of sprayer was observed in primary distribution showing that low densities or early growth stages implied high losses of pesticides on soil. Washed and unwashed samples of lettuce were taken and analyzed from application to harvest to evaluate removal of pesticides by food processing. Results show that residues found on the Spanish preharvest interval days were in all cases below officially set maximum residue limits, although it was observed that time between application and harvest is as important for residues as application amounts. An overall reduction of 40–60% of pesticides residues was obtained from washing lettuce. Experimentally derived residues were compared with modeled residues and deviate from 1.2 to 1.4 for imidacloprid and azoxystrobin, respectively, presenting good model predictions. Resulting human intake fractions range from... for imidacloprid to ... for azoxystrobin.

Prokaryotic growth and presence of metabolic products and contamination tracers in rocks samples from ODP Leg 187 sites

The microbial population in samples of basalt drilled from the north of the Australian Antarctic Discordance (AAD) during Ocean Drilling Program Leg 187 were studied using deoxyribonucleic acid (DNA)-based methods and culturing techniques. The results showed the presence of a microbial population characteristic for the basalt environment. DNA sequence analysis revealed that microbes grouping within the Actinobacteria, green nonsulfur bacteria, the Cytophaga/Flavobacterium/Bacteroides (CFB) group, the Bacillus/Clostridium group, and the beta and gamma subclasses of the Proteobacteria were present in the basalt samples collected. The most dominant phylogenetic group, both in terms of the number of sequences retrieved and the intensities of the DNA bands obtained with the denaturing gradient gel electrophoresis analysis, was the gamma Proteobacteria. Enrichment cultures showed phylogenetic affiliation with the Actinobacteria, the CFB group, the Bacillus/Clostridium group, and the alpha, beta, gamma, and epsilon subclasses of the Proteobacteria. Comparison of native and enriched samples showed that few of the microbes found in native basalt samples grew in the enrichment cultures. Only seven clusters, two clusters within each of the CFB and Bacillus/Clostridium groups and five clusters within the gamma Proteobacteria, contained sequences from both native and enriched basalt samples with significant similarity. Results from cultivation experiments showed the presence of the physiological groups of iron reducers and methane producers. The presence of the iron/manganese-reducing bacterium Shewanella was confirmed with DNA analysis. The results indicate that iron reducers and lithotrophic methanogenic Archaea are indigenous to the ocean crust basalt and that the methanogenic Archaea may be important primary producers in this basaltic environment.

Tidal volume single breath washout of two tracer gases--a practical and promising lung function test

Background Small airway disease frequently occurs in chronic lung diseases and may cause ventilation inhomogeneity (VI), which can be assessed by washout tests of inert tracer gas. Using two tracer gases with unequal molar mass (MM) and diffusivity increases specificity for VI in different lung zones. Currently washout tests are underutilised due to the time and effort required for measurements. The aim of this study was to develop and validate a simple technique for a new tidal single breath washout test (SBW) of sulfur hexafluoride (SF6) and helium (He) using an ultrasonic flowmeter (USFM). Methods The tracer gas mixture contained 5% SF6 and 26.3% He, had similar total MM as air, and was applied for a single tidal breath in 13 healthy adults. The USFM measured MM, which was then plotted against expired volume. USFM and mass spectrometer signals were compared in six subjects performing three SBW. Repeatability and reproducibility of SBW, i.e., area under the MM curve (AUC), were determined in seven subjects performing three SBW 24 hours apart. Results USFM reliably measured MM during all SBW tests (n = 60). MM from USFM reflected SF6 and He washout patterns measured by mass spectrometer. USFM signals were highly associated with mass spectrometer signals, e.g., for MM, linear regression r-squared was 0.98. Intra-subject coefficient of variation of AUC was 6.8%, and coefficient of repeatability was 11.8%. Conclusion The USFM accurately measured relative changes in SF6 and He washout. SBW tests were repeatable and reproducible in healthy adults. We have developed a fast, reliable, and straightforward USFM based SBW method, which provides valid information on SF6 and He washout patterns during tidal breathing.

Using geophysical surveys to test tracer-based storage estimates in headwater catchments

Acknowledgements The authors are grateful to Stian Bradford, Chris Gabrielli, and Julie Timms for practical and logistical assistance. The provision of transport by Iain Malcolm and Ross Glover of Marine Scotland Science was greatly appreciated. We also thank the European Research Council ERC (project GA 335910 VEWA) for funding through the VeWa project and the Leverhulme Trust for funding through PLATO (RPG-2014-016).

Automatic Tracer-Dilution Method Used for Stage-Discharge Ratings and Streamflow Hydrographs on Small Iowa Streams, HR-285, 1990

An automatic system was designed to concurrently measure stage and discharge for the purpose of developing stage-discharge ratings and high flow hydrographs on small streams. Stage, or gage height, is recorded by an analog-to-digital recorder and discharge is determined by the constant-rate tracer-dilution method. The system measures flow above a base stage set by the user. To test the effectiveness of the system and its components, eight systems, with a variety of equipment, were installed at crest-stage gaging stations across Iowa. A fluorescent dye, rhodamine-WT, was used as the tracer. Tracer-dilution discharge measurements were made during 14 flow periods at six stations from 1986 through 1988 water years. Ratings were developed at three stations with the aid of these measurements. A loop rating was identified at one station during rapidly-changing flow conditions. Incomplete mixing and dye loss to sediment apparently were problems at some stations. Stage hydrographs were recorded for 38 flows at seven stations. Limited data on background fluorescence during high flows were also obtained.