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.

Use of high-resolution geophysical data to characterize heterogeneous aquifers: influence of data integration method on hydrological predictions

The integration of geophysical data into the subsurface characterization problem has been shown in many cases to significantly improve hydrological knowledge by providing information at spatial scales and locations that is unattainable using conventional hydrological measurement techniques. The investigation of exactly how much benefit can be brought by geophysical data in terms of its effect on hydrological predictions, however, has received considerably less attention in the literature. Here, we examine the potential hydrological benefits brought by a recently introduced simulated annealing (SA) conditional stochastic simulation method designed for the assimilation of diverse hydrogeophysical data sets. We consider the specific case of integrating crosshole ground-penetrating radar (GPR) and borehole porosity log data to characterize the porosity distribution in saturated heterogeneous aquifers. In many cases, porosity is linked to hydraulic conductivity and thus to flow and transport behavior. To perform our evaluation, we first generate a number of synthetic porosity fields exhibiting varying degrees of spatial continuity and structural complexity. Next, we simulate the collection of crosshole GPR data between several boreholes in these fields, and the collection of porosity log data at the borehole locations. The inverted GPR data, together with the porosity logs, are then used to reconstruct the porosity field using the SA-based method, along with a number of other more elementary approaches. Assuming that the grid-cell-scale relationship between porosity and hydraulic conductivity is unique and known, the porosity realizations are then used in groundwater flow and contaminant transport simulations to assess the benefits and limitations of the different approaches.

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.

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.

Constraints on the permeability structure of alluvial aquifers from the poro-elastic inversion of multi-frequency P-wave velocity logs

We have explored the possibility of obtaining first-order permeability estimates for saturated alluvial sediments based on the poro-elastic interpretation of the P-wave velocity dispersion inferred from sonic logs. Modern sonic logging tools designed for environmental and engineering applications allow one for P-wave velocity measurements at multiple emitter frequencies over a bandwidth covering 5 to 10 octaves. Methodological considerations indicate that, for saturated unconsolidated sediments in the silt to sand range and typical emitter frequencies ranging from approximately 1 to 30 kHz, the observable velocity dispersion should be sufficiently pronounced to allow one for reliable first-order estimations of the permeability structure. The corresponding predictions have been tested on and verified for a borehole penetrating a typical surficial alluvial aquifer. In addition to multifrequency sonic logs, a comprehensive suite of nuclear and electrical logs, an S-wave log, a litholog, and a limited number laboratory measurements of the permeability from retrieved core material were also available. This complementary information was found to be essential for parameterizing the poro-elastic inversion procedure and for assessing the uncertainty and internal consistency of corresponding permeability estimates. Our results indicate that the thus obtained permeability estimates are largely consistent with those expected based on the corresponding granulometric characteristics, as well as with the available evidence form laboratory measurements. These findings are also consistent with evidence from ocean acoustics, which indicate that, over a frequency range of several orders-of-magnitude, the classical theory of poro-elasticity is generally capable of explaining the observed P-wave velocity dispersion in medium- to fine-grained seabed sediments

Simulated-annealing-based conditional simulation for the local-scale characterization of heterogeneous aquifers

Simulated-annealing-based conditional simulations provide a flexible means of quantitatively integrating diverse types of subsurface data. Although such techniques are being increasingly used in hydrocarbon reservoir characterization studies, their potential in environmental, engineering and hydrological investigations is still largely unexploited. Here, we introduce a novel simulated annealing (SA) algorithm geared towards the integration of high-resolution geophysical and hydrological data which, compared to more conventional approaches, provides significant advancements in the way that large-scale structural information in the geophysical data is accounted for. Model perturbations in the annealing procedure are made by drawing from a probability distribution for the target parameter conditioned to the geophysical data. This is the only place where geophysical information is utilized in our algorithm, which is in marked contrast to other approaches where model perturbations are made through the swapping of values in the simulation grid and agreement with soft data is enforced through a correlation coefficient constraint. Another major feature of our algorithm is the way in which available geostatistical information is utilized. Instead of constraining realizations to match a parametric target covariance model over a wide range of spatial lags, we constrain the realizations only at smaller lags where the available geophysical data cannot provide enough information. Thus we allow the larger-scale subsurface features resolved by the geophysical data to have much more due control on the output realizations. Further, since the only component of the SA objective function required in our approach is a covariance constraint at small lags, our method has improved convergence and computational efficiency over more traditional methods. Here, we present the results of applying our algorithm to the integration of porosity log and tomographic crosshole georadar data to generate stochastic realizations of the local-scale porosity structure. Our procedure is first tested on a synthetic data set, and then applied to data collected at the Boise Hydrogeophysical Research Site.

Ambient vertical flow in long-screen wells: a case study in the Fontainebleau Sands Aquifer (France)

A tritium (H-3) profile was constructed in a long-screened well (LSW) of the Fontainebleau Sands Aquifer (France), and the data were combined with temperature logs to gain insight into the potential effects of the ambient vertical flow (AVF) of water through the well on the natural aquifer stratification. AVF is commonly taken into account in wells located in fracture aquifers or intercepting two different aquifers with distinct hydraulic heads. However, due to the vertical hydraulic gradient of the flow lines intercepted by wells, AVF of groundwater is a common process within any type of aquifer. The detection of 3H in the deeper parts of the studied well ( approximate depth 50m), where H-3-free groundwater is expected, indicates that shallow young water is being transported downwards through the well itself. The temperature logs show a nearly zero gradient with depth, far below the mean geothermal gradient in sedimentary basins. The results show that the age distribution of groundwater samples might be biased in relation to the age distribution in the surroundings of the well. The use of environmental tracers to investigate aquifer properties, particularly in LSWs, is then limited by the effects of the AVF of water that naturally occurs through the well.

Groundwater dating and past climate reconstruction

In the past decade aquifers have increasingly become palaeoclimatic archives in their own right alongside ice cores, sediments and other proxy records. The main tool for this task has been the noble gas palaeo-thermometer in combination with quantitative groundwater dating using radionuclides. Noblegas radionuclides play a unique role as tracers in environmental studies due to their chemical inertness and low concentration making them ideal tracers. The same properties on the other hand make them difficult to measure on natural concentration levels. Therefore for decades low level counting (LLC) was the only method for detecting radioisotopes of argon and krypton at an atmospheric level. In recent times and with the increase of interest and potential applications the analytical efforts with novel detection methods have been intensified. In the talk noble gas groundwater dating techniques over times scales from decades to millions of years are discussed in relation to noble gas palaeo records at different locations in Europe and elsewhere.

L'exploitation des eaux souterraines en Syrie centrale : construction de diagnostics et politiques d'intervention

Syrian dry areas have been for several millennia a place of interaction between human populations and the environment. If environmental constraints and heterogeneity condition the human occupation and exploitation of resources, socio-political, economic and historical elements play a fundamental role. Since the late 1980s, Syrian dry areas are viewed as suffering a serious water crisis, due to groundwater overdraft. The Syrian administration and international development agencies believe that groundwater overexploitation is also leading to a decline of agricultural activities and to poverty increase. Action is thus required to address these problems.However, the overexploitation diagnosis needs to be reviewed. The overexploitation discourse appears in the context of Syria's opening to international organizations and to the market economy. It echoes the international discourse of "global water crisis". The diagnosis is based on national indicators recycling old Soviet data that has not been updated. In the post-Soviet era, the Syrian national water policy seems to abandon large surface water irrigation projects in favor of a strategy of water use rationalization and groundwater conservation in crisis regions, especially in the district of Salamieh.This groundwater conservation policy has a number of inconsistencies. It is justified for the administration and also probably for international donors, since it responds to an indisputable environmental emergency. However, efforts to conserve water are anecdotal or even counterproductive. The water conservation policy appears a posteriori as an extension of the national policy of food self-sufficiency. The dominant interpretation of overexploitation, and more generally of the water crisis, prevents any controversary approach of the status of resources and of the agricultural system in general and thus destroys any attempt to discuss alternatives with respect to groundwater management, allocation, and their inclusion in development programs.A revisited diagnosis of the situation needs to take into account spatial and temporal dimensions of the groundwater exploitation and to analyze the co-evolution of hydrogeological and agricultural systems. It should highlight the adjustments adopted to cope with environmental and economic variability, changes of water availability and regulatory measures enforcements. These elements play an important role for water availability and for the spatial, temporal, sectoral allocation of water resource. The groundwater exploitation in the last century has obviously had an impact on the environment, but the changes are not necessarily catastrophic.The current groundwater use in central Syria increases the uncertainty by reducing the ability of aquifers to buffer climatic changes. However, the climatic factor is not the only source of uncertainty. The high volatility of commodity prices, fuel, land and water, depending on the market but also on the will (and capacity) of the Syrian State to preserve social peace is a strong source of uncertainty. The research should consider the whole range of possibilities and propose alternatives that take into consideration the risks they imply for the water users, the political will to support or not the local access to water - thus involving a redefinition of the economic and social objectives - and finally the ability of international organizations to reconsider pre-established diagnoses.

Reorganization of a deeply incised drainage: role of deformation, sedimentation and groundwater flow

Deeply incised drainage networks are thought to be robust and not easily modified, and are commonly used as passive markers of horizontal strain. Yet, reorganizations (rearrangements) appear in the geologic record. We provide field evidence of the reorganization of a Miocene drainage network in response to strike-slip and vertical displacements in Guatemala. The drainage was deeply incised into a 50-km-wide orogen located along the North America-Caribbean plate boundary. It rearranged twice, first during the Late Miocene in response to transpressional uplift along the Polochic fault, and again in the Quaternary in response to transtensional uplift along secondary faults. The pattern of reorganization resembles that produced by the tectonic defeat of rivers that cross growing tectonic structures. Compilation of remote sensing data, field mapping, sediment provenance study, grain-size analysis and Ar(40)/Ar(39) dating from paleovalleys and their fill reveals that the classic mechanisms of river diversion, such as river avulsion over bedrock, or capture driven by surface runoff, are not sufficient to produce the observed diversions. The sites of diversion coincide spatially with limestone belts and reactivated fault zones, suggesting that solution-triggered or deformation-triggered permeability have helped breaching of interfluves. The diversions are also related temporally and spatially to the accumulation of sediment fills in the valleys, upstream of the rising structures. We infer that the breaching of the interfluves was achieved by headward erosion along tributaries fed by groundwater flow tracking from the valleys soon to be captured. Fault zones and limestone belts provided the pathways, and the aquifers occupying the valley fills provided the head pressure that enhanced groundwater circulation. The defeat of rivers crossing the rising structures results essentially from the tectonically enhanced activation of groundwater flow between catchments.

Stable isotope compositions of mammoth teeth from Niederweningen, Switzerland: Implications for the Late Pleistocene climate, environment, and diet

Oxygen and carbon isotope compositions of well-preserved mammoth teeth from the Middle Wurmian (40-70 ka) peat layer of Niederweningen, the most important mammoth site in Switzerland, were analysed to reconstruct Late Pleistocene palaeoclimatic and palaeoenvironmental conditions. Drinking water (delta(18)O values of approximately -12.3 +/- 0.9 parts per thousand were calculated front oxygen isotope compositions of mammoth tooth enamel apatite using a species-specific calibration for modern elephants. These delta(18)O(H2O) values reflect the mean oxygen isotope composition of the palaeo-precipitation and are similar to those directly measured for fate Pleistocene groundwater from aquifers in northern Switzerland and southern Germany. Using a present-day delta(18)O(H2)o-precipitation-air temperature relation for Switzerland, a mean annual air temperature (MAT) of around 4.3 +/- 2.1 degrees C can be calculated for the Middle Wurmian at this site. This MAT is in good agreement with palaeotemperature estimates on the basis of Middle Wurmian groundwater recharge temperatures and beetle assemblages. Hence, the climatic conditions in this region were around 4 degrees C cooler during the Middle Wurmian interstadial phase, around 45-50ka BP, than they are today. During this period the mammoths from Niederweningen lived in an open tundra-like, C(3) plant-dominated environment as indicated by enamel (delta(13)C values of -11.5 +/- 0.3 parts per thousand and pollen and macroplant fossils found in the embedding peat. The low variability of enamel delta(13)C and delta(18)O values from different mammoth teeth reflects similar environmental conditions and supports a relatively small time frame for the fossil assemblage. (C) 2006 Elsevier Ltd and INQUA. All rights reserved.

Open to closed system transition traced through the TDIC isotopic signature at the aquifer recharge stage, implications for groundwater 14C dating

14C dating models are limited when considering recent groundwater for which the carbon isotopic signature of the total dissolved inorganic carbon (TDIC) is mainly acquired in the unsaturated zone. Reducing the uncertainties of dating thus implies a better identification of the processes controlling the carbon isotopic composition of the TDIC during groundwater recharge. Geochemical interactions between gas, water and carbonates in the unsaturated zone were investigated for two aquifers (the carbonate-free Fontainebleau sands and carbonate-bearing Astian sands, France) in order to identify the respective roles of CO2 and carbonates on the carbon isotopic signatures of the TDIC; this analysis is usually approached using open or closed system terms. Under fully open system conditions, the seasonality of the 13C values in the soil CO2 can lead to important uncertainties regarding the so-called "initial 14C activity" used in 14C correction models. In a carbonate-bearing unsaturated zone such as in the Astian aquifer, we show that an approach based on fully open or closed system conditions is not appropriate. Although the chemical saturation between water and calcite occurs rapidly within the first metre of the unsaturated zone, the carbon isotopic contents (δ13C) of the CO2 and the TDIC evolve downward, impacted by the dissolution-precipitation of the carbonates. In this study, we propose a numerical approach to describe this evolution. The δ13C and the A 14C (radiocarbon activity) of the TDIC at the base of the carbonate-hearing unsaturated zone depends on (i) the δ13C and the A 14C of the TDIC in the soil determined by the soil CO2, (ii) the water's residence time in the unsaturated zone and (iii) the carbonate precipitation-dissolution fluxes. In this type of situation, the carbonate δ13C-A 14C evolutions indicate the presence of secondary calcite and permit the calculation of its accretion flux, equal to ~ 4.5 ± 0.5 x 10-9 mol grock-1 yr-1. More generally, for other sites under temperate climate and with similar properties to the Astian sands site, this approach allows for a reliable determination of the carbon isotopic composition at the base of the unsaturated zone as the indispensable "input function" data of the carbon cycle into the aquifer.

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.