Hydrogeophysical data integration at larger scales: Application of Bayesian sequential simulation for the characterization of heterogeneous alluvial aquifers

Knowledge of the spatial distribution of hydraulic conductivity (K) within an aquifer is critical for reliable predictions of solute transport and the development of effective groundwater management and/or remediation strategies. While core analyses and hydraulic logging can provide highly detailed information, such information is inherently localized around boreholes that tend to be sparsely distributed throughout the aquifer volume. Conversely, larger-scale hydraulic experiments like pumping and tracer tests provide relatively low-resolution estimates of K in the investigated subsurface region. As a result, traditional hydrogeological measurement techniques contain a gap in terms of spatial resolution and coverage, and they are often alone inadequate for characterizing heterogeneous aquifers. Geophysical methods have the potential to bridge this gap. The recent increased interest in the application of geophysical methods to hydrogeological problems is clearly evidenced by the formation and rapid growth of the domain of hydrogeophysics over the past decade (e.g., Rubin and Hubbard, 2005).

Anatomy of contaminated aquifers of an industrial site: insights from the stable isotope compositions of waters and dissolved inorganic carbon

The hydrogen and oxygen isotopes of water and the carbon isotope composition of dissolved inorganic carbon (DIC) from different aquifers at an industrial site, highly contaminated by organic pollutants representing residues of the former gas production, have been used as natural tracers to characterize the hydrologic system. On the basis of their stable isotope compositions as well as the seasonal variations, different groups of waters (precipitation, surface waters, groundwaters and mineral waters) as well as seasonably variable processes of mixing between these waters can clearly be distinguished. In addition, reservoir effects and infiltration rates can be estimated. In the northern part of the site an influence of uprising mineral waters within the Quaternary aquifers, presumably along a fault zone, can be recognized. Marginal infiltration from the Neckar River in the cast and surface water infiltration adjacent to a steep hill on the western edge of the site with an infiltration rate of about one month can also be resolved through the seasonal variation. Quaternary aquifers closer to the centre of the site show no seasonal variations, except for one borehole close to a former mill channel and another borehole adjacent to a rain water channel. Distinct carbon isotope compositions and concentrations of DIC for these different groups of waters reflect variable influence of different components of the natural carbon cycle: dissolution of marine carbonates in the mineral waters, biogenic, soil-derived CO2 in ground- and surface waters, as well as additional influence of atmospheric CO2 for the surface waters. Many Quaternary aquifer waters have, however, distinctly lower delta(13)C(DIC) values and higher DIC concentrations compared to those expected for natural waters. Given the location of contaminated groundwaters at this site but also in the industrially well-developed valley outside of this site, the most likely source for the low C-13(DIC) values is a biodegradation of anthropogenic organic substances, in particular the tar oils at the site.

Permeability estimates from poro-elastic dispersion analyses of sonic logs

Modern sonic logging tools designed for shallow environmental and engineering applications allow for P-wave phase velocity measurements over a wide frequency band. Methodological considerations indicate that, for saturated unconsolidated sediments in the silt to sand range and source frequencies ranging from approximately 1 to 30 kHz, the observable poro-elastic P-wave velocity dispersion is sufficiently pronounced to allow for reliable first-order estimations of the underlying permeability structure. These predictions have been tested on and verified for a surficial alluvial aquifer. Our results indicate that, even without any further calibration, the thus obtained permeability estimates as well as their variabilities within the pertinent lithological units are remarkably close to those expected based on the corresponding granulometric characteristics.

Radiocarbon dating of modern groundwater: the role of the unsaturated zone

Biological and physical processes occurring in soils may lead to significant isotopic changes between the isotopic compositions of atmospheric CO2 and of soil CO2. Also, during water and gas transport from the soil surface to the water table, isotopic changes likely occur due to numerous physical processes such as gas production and diffusion, water advection, and gas-water-rock interactions. In most cases, these changes are not included in the correction models developed for groundwater dating, whereas they can significantly impact the calculation of the 14C age. We explore the role of these processes using: i) experimental data from two aquifer sites (Fontainebleau sands and Astian sands, France), ii) a distributed model to simulate the 14C activities of soil CO2, and iii) numerical simulations in order to highlight the role of the physical processes.¦The 13C content in soil CO2 showed seasonal variations and highlighted the competition between CO2 production and CO2 diffusion. Their respective contributions played a significant role in defining the isotopic composition of CO2 at the water table. On both study sites, variations of the 14C activity in soil CO2 reflect the competition between the fluxes of root derived-CO2 and organic matter derived-CO2. Since the nuclear weapon tests in the fifties and sixties, soil CO2 became significantly depleted in 14C compared to modern atmospheric CO2. Models that take into account this 14C depletion in soil CO2 for dating modern groundwater would lead to apparent younger 14C ages than models that only consider the 14C activity in atmospheric CO2. Moreover, since 2000-2005, the inverse effect is observed as soil CO2 is enriched in 14C compared to atmospheric CO2.¦Therefore, we conclude that the isotopic composition of CO2 at the water table have to be taken into account for the dating of modern groundwater. This requires a systematic sampling of soil CO2 and the measurement of its 13C and 14C contents. We used this information in a numerical simulation to calculate the evolution of isotopic composition of CO2 from the soil surface to the water table. This simulation integrated physical processes in the unsaturated zone (e.g. CO2 production and diffusion, water advection, etc.) and gas-water-rock interactions.

The role of basalt replenishment in the generation of basaltic andesites of the ongoing activity at Arenal volcano, Costa Rica: evidence from clinopyroxene and spinel

The bulk composition of magma erupted from Volcan Arenal has remained nearly constant (SiO2 = 53.6-54.9 wt%; MgO = 5.0-4.5 wt%) during almost 30 years of continuous activity (1969-1996). None the less, clinopyroxene (cpx) phenocrysts and their spinel inclusions record a much more complex open-system evolution in which steady-state production of the erupted basaltic andesitic magma is linked to episodic injections of basalt into Arenal's magma conduit/reservoir system. High-resolution major element zoning profiles (electron microprobe) on a large number of phenocrysts (>14,000 analyses), tied to back-scattered electron (BSE) images, have been used to assess the compositional characteristics of the magmatic end members as well as the timing and dynamics of magma replenishment events. No two cpx phenocrysts have exactly the same zoning profile. The vast majority of our analyses record the crystallization of cpx (Cr2O3 < 0.12 wt%; Mg# = 65-79; Al/Ti = 2-7) from a liquid comparable to or more evolved than erupted magma compositions. However, half of all cpx grains are cored by high-Cr cpx (Cr2O3 = 0.2-0.72 wt%) or contain similar basaltic compositions as abrupt growth bands in phenocrysts with and without high-Cr cores; phenocrysts with high-Cr cpx occur throughout the ongoing activity. In a few cases, high-Cr cpx occurs very near the outer margin of the grain without an apparent growth hiatus, particularly in 1968/69 and 1992/93. The main conclusions are: (1) all basaltic andesitic lavas erupted at Arenal during the ongoing activity that began in July, 1968, are the products of magma mixing, (2) clinopyroxenes record multiple replenishment events of basaltic magma in contrast to the near constancy of erupted bulk compositions, (3) some phenocrysts preserve records of multiple interactions with basaltic magmas requiring magmatic processes to operate on time-scales shorter than residence times of some phenocrysts, (4) multiple occurrences of clinopyroxene with high-Cr rims suggest that basalt replenishment events have occurred with sub-decadal frequency and may predate eruption by months or less. From this we infer that Arenal volcano is underlain by a continuously active, small-volume magmatic reservoir maintained in quasi-steady state by basalt recharge over several decades. The monotony of erupting Arenal magmas implies that fractionation, recharge, ascent, and eruption are well balanced in order for magmas to be essentially uniform while containing phenocrysts with vastly different growth histories at the time of eruption.

La gravimétrie expéditive appliquée à la recherche d'aquifères en zone aride. Cas de la nappe alluviale du Teloua (Agadez, Niger)

This paper illustrates the practicality and efficiency of gravimetry for aquifer prospecting in arid zones. Known for the long and tedious data-processing it requires, this method becomes expeditious when simplified as presented here. Its use is then fully justified in a survey of this kind. During the study of the Teloua alluvial aquifer (Agadez, Niger), several ancient channels were clearly and rapidly located. Comparison of the results obtained here with those from previous studies demonstrates anew that for comprehensive prospecting, several complementary geophysical methods should always be employed.

Constraints on the permeability structure of alluvial aquifers from P-wave sonic logs

Modern sonic logging tools designed for shallow environmental and engineering applications allow for P-wave phase velocity measurements over a wide frequency band. Methodological considerations indicate that, for saturated unconsolidated sediments in the silt to sand range and source frequencies ranging from approximately 1 to 30 kHz, the observable poro-elastic P-wave velocity dispersion is sufficiently pronounced to allow for reliable first-order estimations of the underlying permeability structure. These predictions have been tested on and verified for a surficial alluvial aquifer. Our results indicate that, even without any further calibration, the thus obtained permeability estimates as well as their variabilities within the pertinent lithological units are remarkably close to those expected based on the corresponding granulometric characteristics.

Water chemistry and isotope composition of the Acquarossa thermal system, Ticino, Switzerland

The thermal springs of Acquarossa and the nearby mineral springs of Soia have outlet temperatures of 12 degrees to 25 degrees C, TDS of 2290 to 3000 mg/kg and Ca-SO4 to Ca-SO4-HCO3 composition. Chemical geothermometers suggest reservoir temperatures close to 60 degrees C. P-CO2 values at depth are estimated to range from 0.3 to 2 bar. delta D and delta(18)O values indicate a meteoric origin and recharge elevations of 1600 +/- 150 m above sea level (a.s.l.) for these thermal and mineral waters. All these waters discharge from the overturned limb of the Simano nappe, probably dose to the contact between basement and underlying cover rocks. They therefore represent rain waters that descend slowly, heat at depth and locally rise relatively quickly to the surface, preserving part of their physical and chemical characteristics. (C) 1999 CNR. Published by Elsevier Science Ltd. All rights reserved.

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.

Geophysical Investigations for Hydrogeological Characterization at the Kappelen Test Site (Switzerland)

Vertical electric soundings, 2D resistivity imaging and several logging measurements were performed at Kappelen test site to identify the various geolelectric facies that allowed determining the tabular and horizontal structure of the aquifer. The surface-based geoelectric methods allowed for a reliable characterization of the overall structure and the geometry of the aquifer, while geophysical logging methods allowed for inferring detailed hydrogeophysical characteristics, such as the electrical resistivity, total porosity, global and matrix density and hydraulic conductivity. The synoptic interpretation and integration of this broad and diverse database allows for constraining the key hydrological characteristics and hence forms the basis for the detailed hydraulic modelling of flow and transport process.