Remediation of a marine shore tailings deposit and the importance of water-rock interaction on element cycling in the coastal aquifer

We present the study of the geochemical processes associated with the first successful remediation of a marine shore tailings deposit in a coastal desert environment (Bahia de Ite, in the Atacama Desert of Peru). The remediation approach implemented a wetland on top of the oxidized tailings. The site is characterized by a high hydrauliz gradient produced by agricultural irrigation on upstream gravel terraces that pushed river water (similar to 500 mg/L SO(4)) toward the sea and through the tailings deposit. The geochemical and isotopic (delta(2)H(water) and delta(18)O(water), delta(34)S(sulfate) , delta(18)O(sulfate)) approach applied here revealed that evaporite horizons (anhydrite and halite) in the gravel terraces are the source of increased concentrations of SO(4), Cl, and Na up to similar to 1500 mg/L in the springs at the base of the gravel terraces. Deeper groundwater interacting with underlying marine sequences increased the concentrations of SO(4), Cl, and Na up to 6000 mg/L and increased the alkalinity up to 923 mg/L CaCO(3) eq. in the coastal aquifer. These waters infiltrated into the tailings deposit at the shelf-tailings interface. Nonremediated tailings had a low-pH oxidation zone (pH 1-4) with significant accumulations of efflorescent salts (10-20 cm thick) at the surface because of upward capillary transport of metal cations in the arid climate. Remediated tailings were characterized by neutral pH and reducing conditions (pH similar to 7, Eh similar to 100 mV). As a result, most bivalent metals such as Cu, Zn, and Ni had very low concentrations (around 0.01 mg/L or below detection limit) because of reduction and sorption processes. In contrast, these reducing conditions increased the mobility of iron from two sources in this system: (1) The originally Fe(III)-rich oxidation zone, where Fe(II) was reduced during the remediation process and formed an Fe(II) plume, and (2) reductive dissolution of Fe(III) oxides present in the original shelf lithology formed an Fe-Mn plume at 10-m depth. These two Fe-rich plumes were pushed toward the shoreline where more oxidizing and higher pH conditions triggered the precipitation of Fe(HI)hydroxide coatings on silicates. These coatings acted as a filter for the arsenic, which naturally infiltrated with the river water (similar to 500 mu g/L As natural background) into the tailings deposit.

Error models in hydrogeology applications

Notre consommation en eau souterraine, en particulier comme eau potable ou pour l'irrigation, a considérablement augmenté au cours des années. De nombreux problèmes font alors leur apparition, allant de la prospection de nouvelles ressources à la remédiation des aquifères pollués. Indépendamment du problème hydrogéologique considéré, le principal défi reste la caractérisation des propriétés du sous-sol. Une approche stochastique est alors nécessaire afin de représenter cette incertitude en considérant de multiples scénarios géologiques et en générant un grand nombre de réalisations géostatistiques. Nous rencontrons alors la principale limitation de ces approches qui est le coût de calcul dû à la simulation des processus d'écoulements complexes pour chacune de ces réalisations. Dans la première partie de la thèse, ce problème est investigué dans le contexte de propagation de l'incertitude, oú un ensemble de réalisations est identifié comme représentant les propriétés du sous-sol. Afin de propager cette incertitude à la quantité d'intérêt tout en limitant le coût de calcul, les méthodes actuelles font appel à des modèles d'écoulement approximés. Cela permet l'identification d'un sous-ensemble de réalisations représentant la variabilité de l'ensemble initial. Le modèle complexe d'écoulement est alors évalué uniquement pour ce sousensemble, et, sur la base de ces réponses complexes, l'inférence est faite. Notre objectif est d'améliorer la performance de cette approche en utilisant toute l'information à disposition. Pour cela, le sous-ensemble de réponses approximées et exactes est utilisé afin de construire un modèle d'erreur, qui sert ensuite à corriger le reste des réponses approximées et prédire la réponse du modèle complexe. Cette méthode permet de maximiser l'utilisation de l'information à disposition sans augmentation perceptible du temps de calcul. La propagation de l'incertitude est alors plus précise et plus robuste. La stratégie explorée dans le premier chapitre consiste à apprendre d'un sous-ensemble de réalisations la relation entre les modèles d'écoulement approximé et complexe. Dans la seconde partie de la thèse, cette méthodologie est formalisée mathématiquement en introduisant un modèle de régression entre les réponses fonctionnelles. Comme ce problème est mal posé, il est nécessaire d'en réduire la dimensionnalité. Dans cette optique, l'innovation du travail présenté provient de l'utilisation de l'analyse en composantes principales fonctionnelles (ACPF), qui non seulement effectue la réduction de dimensionnalités tout en maximisant l'information retenue, mais permet aussi de diagnostiquer la qualité du modèle d'erreur dans cet espace fonctionnel. La méthodologie proposée est appliquée à un problème de pollution par une phase liquide nonaqueuse et les résultats obtenus montrent que le modèle d'erreur permet une forte réduction du temps de calcul tout en estimant correctement l'incertitude. De plus, pour chaque réponse approximée, une prédiction de la réponse complexe est fournie par le modèle d'erreur. Le concept de modèle d'erreur fonctionnel est donc pertinent pour la propagation de l'incertitude, mais aussi pour les problèmes d'inférence bayésienne. Les méthodes de Monte Carlo par chaîne de Markov (MCMC) sont les algorithmes les plus communément utilisés afin de générer des réalisations géostatistiques en accord avec les observations. Cependant, ces méthodes souffrent d'un taux d'acceptation très bas pour les problèmes de grande dimensionnalité, résultant en un grand nombre de simulations d'écoulement gaspillées. Une approche en deux temps, le "MCMC en deux étapes", a été introduite afin d'éviter les simulations du modèle complexe inutiles par une évaluation préliminaire de la réalisation. Dans la troisième partie de la thèse, le modèle d'écoulement approximé couplé à un modèle d'erreur sert d'évaluation préliminaire pour le "MCMC en deux étapes". Nous démontrons une augmentation du taux d'acceptation par un facteur de 1.5 à 3 en comparaison avec une implémentation classique de MCMC. Une question reste sans réponse : comment choisir la taille de l'ensemble d'entrainement et comment identifier les réalisations permettant d'optimiser la construction du modèle d'erreur. Cela requiert une stratégie itérative afin que, à chaque nouvelle simulation d'écoulement, le modèle d'erreur soit amélioré en incorporant les nouvelles informations. Ceci est développé dans la quatrième partie de la thèse, oú cette méthodologie est appliquée à un problème d'intrusion saline dans un aquifère côtier. -- Our consumption of groundwater, in particular as drinking water and for irrigation, has considerably increased over the years and groundwater is becoming an increasingly scarce and endangered resource. Nofadays, we are facing many problems ranging from water prospection to sustainable management and remediation of polluted aquifers. Independently of the hydrogeological problem, the main challenge remains dealing with the incomplete knofledge of the underground properties. Stochastic approaches have been developed to represent this uncertainty by considering multiple geological scenarios and generating a large number of realizations. The main limitation of this approach is the computational cost associated with performing complex of simulations in each realization. In the first part of the thesis, we explore this issue in the context of uncertainty propagation, where an ensemble of geostatistical realizations is identified as representative of the subsurface uncertainty. To propagate this lack of knofledge to the quantity of interest (e.g., the concentration of pollutant in extracted water), it is necessary to evaluate the of response of each realization. Due to computational constraints, state-of-the-art methods make use of approximate of simulation, to identify a subset of realizations that represents the variability of the ensemble. The complex and computationally heavy of model is then run for this subset based on which inference is made. Our objective is to increase the performance of this approach by using all of the available information and not solely the subset of exact responses. Two error models are proposed to correct the approximate responses follofing a machine learning approach. For the subset identified by a classical approach (here the distance kernel method) both the approximate and the exact responses are knofn. This information is used to construct an error model and correct the ensemble of approximate responses to predict the "expected" responses of the exact model. The proposed methodology makes use of all the available information without perceptible additional computational costs and leads to an increase in accuracy and robustness of the uncertainty propagation. The strategy explored in the first chapter consists in learning from a subset of realizations the relationship between proxy and exact curves. In the second part of this thesis, the strategy is formalized in a rigorous mathematical framework by defining a regression model between functions. As this problem is ill-posed, it is necessary to reduce its dimensionality. The novelty of the work comes from the use of functional principal component analysis (FPCA), which not only performs the dimensionality reduction while maximizing the retained information, but also allofs a diagnostic of the quality of the error model in the functional space. The proposed methodology is applied to a pollution problem by a non-aqueous phase-liquid. The error model allofs a strong reduction of the computational cost while providing a good estimate of the uncertainty. The individual correction of the proxy response by the error model leads to an excellent prediction of the exact response, opening the door to many applications. The concept of functional error model is useful not only in the context of uncertainty propagation, but also, and maybe even more so, to perform Bayesian inference. Monte Carlo Markov Chain (MCMC) algorithms are the most common choice to ensure that the generated realizations are sampled in accordance with the observations. Hofever, this approach suffers from lof acceptance rate in high dimensional problems, resulting in a large number of wasted of simulations. This led to the introduction of two-stage MCMC, where the computational cost is decreased by avoiding unnecessary simulation of the exact of thanks to a preliminary evaluation of the proposal. In the third part of the thesis, a proxy is coupled to an error model to provide an approximate response for the two-stage MCMC set-up. We demonstrate an increase in acceptance rate by a factor three with respect to one-stage MCMC results. An open question remains: hof do we choose the size of the learning set and identify the realizations to optimize the construction of the error model. This requires devising an iterative strategy to construct the error model, such that, as new of simulations are performed, the error model is iteratively improved by incorporating the new information. This is discussed in the fourth part of the thesis, in which we apply this methodology to a problem of saline intrusion in a coastal aquifer.

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.

Enterobiasis and strongyloidiasis and associated co-infections and morbidity markers in infants, preschool- and school-aged children from rural coastal Tanzania: a cross-sectional study.

Background. There is a paucity of data pertaining to the epidemiology and public health impact of Enterobius vermicularis and Strongyloides stercoralis infections. We aimed to determine the extent of enterobiasis, strongyloidiasis, and other helminth infections and their association with asymptomatic Plasmodium parasitaemia, anaemia, nutritional status, and blood cell counts in infants, preschool-aged (PSAC), and school-aged children (SAC) from rural coastal Tanzania.MethodsA total of 1,033 children were included in a cross-sectional study implemented in the Bagamoyo district in 2011/2012. Faecal samples were examined for intestinal helminth infections using a broad set of quality controlled methods. Finger-prick blood samples were subjected to filariasis and Plasmodium parasitaemia testing and full blood cell count examination. Weight, length/height, and/or mid-upper arm circumference were measured and the nutritional status determined in accordance with age.Results E. vermicularis infections were found in 4.2% of infants, 16.7%, of PSAC, and 26.3% of SAC. S. stercoralis infections were detected in 5.8%, 7.5%, and 7.1% of infants, PSAC, and SAC, respectively. Multivariable regression analyses revealed higher odds of enterobiasis in children of all age-groups with a reported anthelminthic treatment history over the past six months (odds ratio (OR): 2.15; 95% confidence interval (CI): 1.22 - 3.79) and in SAC with a higher temperature (OR: 2.21; CI: 1.13 - 4.33). Strongyloidiasis was associated with eosinophilia (OR: 2.04; CI: 1.20-3.48) and with Trichuris trichiura infections (OR: 4.13; CI: 1.04-16.52) in children of all age-groups, and with asymptomatic Plasmodium parasitaemia (OR: 13.03; CI: 1.34 - 127.23) in infants. None of the investigated helminthiases impacted significantly on the nutritional status and anaemia, but moderate asymptomatic Plasmodium parasitaemia was a strong predictor for anaemia in children aged older than two years (OR: 2.69; 95% CI: 1.23 ¿ 5.86).Conclusions E. vermicularis and S. stercoralis infections were moderately prevalent in children from rural coastal Tanzania. Our data can contribute to inform yet missing global burden of disease and prevalence estimates for strongyloidiasis and enterobiasis. The association between S stercoralis and asymptomatic Plasmodium parasitaemia found here warrants further comprehensive investigations.

Fracture imaging within a granitic rock aquifer using multiple-offset single-hole and cross-hole GPR reflection data

The sparsely spaced highly permeable fractures of the granitic rock aquifer at Stang-er-Brune (Brittany, France) form a well-connected fracture network of high permeability but unknown geometry. Previous work based on optical and acoustic logging together with single-hole and cross-hole flowmeter data acquired in 3 neighbouring boreholes (70-100 m deep) has identified the most important permeable fractures crossing the boreholes and their hydraulic connections. To constrain possible flow paths by estimating the geometries of known and previously unknown fractures, we have acquired, processed and interpreted multifold, single- and cross-hole GPR data using 100 and 250 MHz antennas. The GPR data processing scheme consisting of timezero corrections, scaling, bandpass filtering and F-X deconvolution, eigenvector filtering, muting, pre-stack Kirchhoff depth migration and stacking was used to differentiate fluid-filled fracture reflections from source generated noise. The final stacked and pre-stack depth-migrated GPR sections provide high-resolution images of individual fractures (dipping 30-90°) in the surroundings (2-20 m for the 100 MHz antennas; 2-12 m for the 250 MHz antennas) of each borehole in a 2D plane projection that are of superior quality to those obtained from single-offset sections. Most fractures previously identified from hydraulic testing can be correlated to reflections in the single-hole data. Several previously unknown major near vertical fractures have also been identified away from the boreholes.

Reconstruction of metal pollution and recent sedimentation processes in Havana Bay (Cuba): A tool for coastal ecosystem management.

Since 1998 the highly polluted Havana Bay ecosystem has been the subject of a mitigation program. In order to determine whether pollution-reduction strategies were effective, we have evaluated the historical trends of pollution recorded in sediments of the Bay. A sediment core was dated radiometrically using natural and artificial fallout radionuclides. An irregularity in the (210)Pb record was caused by an episode of accelerated sedimentation. This episode was dated to occur in 1982, a year coincident with the heaviest rains reported in Havana over the XX century. Peaks of mass accumulation rates (MAR) were associated with hurricanes and intensive rains. In the past 60 years, these maxima are related to strong El Niño periods, which are known to increase rainfall in the north Caribbean region. We observed a steady increase of pollution (mainly Pb, Zn, Sn, and Hg) since the beginning of the century to the mid 90s, with enrichment factors as high as 6. MAR and pollution decreased rapidly after the mid 90s, although some trace metal levels remain high. This reduction was due to the integrated coastal zone management program introduced in the late 90s, which dismissed catchment erosion and pollution.

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.

A filtering method to correct time-lapse 3D ERT data and improve imaging of natural aquifer dynamics

We have developed a processing methodology that allows crosshole ERT (electrical resistivity tomography) monitoring data to be used to derive temporal fluctuations of groundwater electrical resistivity and thereby characterize the dynamics of groundwater in a gravel aquifer as it is infiltrated by river water. Temporal variations of the raw ERT apparent-resistivity data were mainly sensitive to the resistivity (salinity), temperature and height of the groundwater, with the relative contributions of these effects depending on the time and the electrode configuration. To resolve the changes in groundwater resistivity, we first expressed fluctuations of temperature-detrended apparent-resistivity data as linear superpositions of (i) time series of riverwater-resistivity variations convolved with suitable filter functions and (ii) linear and quadratic representations of river-water-height variations multiplied by appropriate sensitivity factors; river-water height was determined to be a reliable proxy for groundwater height. Individual filter functions and sensitivity factors were obtained for each electrode configuration via deconvolution using a one month calibration period and then the predicted contributions related to changes in water height were removed prior to inversion of the temperature-detrended apparent-resistivity data. Applications of the filter functions and sensitivity factors accurately predicted the apparent-resistivity variations (the correlation coefficient was 0.98). Furthermore, the filtered ERT monitoring data and resultant time-lapse resistivity models correlated closely with independently measured groundwater electrical resistivity monitoring data and only weakly with the groundwater-height fluctuations. The inversion results based on the filtered ERT data also showed significantly less inversion artefacts than the raw data inversions. We observed resistivity increases of up to 10% and the arrival time peaks in the time-lapse resistivity models matched those in the groundwater resistivity monitoring data.

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.

Zonation for 3D aquifer characterization based on joint inversions of multimethod crosshole geophysical data

Predictive groundwater modeling requires accurate information about aquifer characteristics. Geophysical imaging is a powerful tool for delineating aquifer properties at an appropriate scale and resolution, but it suffers from problems of ambiguity. One way to overcome such limitations is to adopt a simultaneous multitechnique inversion strategy. We have developed a methodology for aquifer characterization based on structural joint inversion of multiple geophysical data sets followed by clustering to form zones and subsequent inversion for zonal parameters. Joint inversions based on cross-gradient structural constraints require less restrictive assumptions than, say, applying predefined petro-physical relationships and generally yield superior results. This approach has, for the first time, been applied to three geophysical data types in three dimensions. A classification scheme using maximum likelihood estimation is used to determine the parameters of a Gaussian mixture model that defines zonal geometries from joint-inversion tomograms. The resulting zones are used to estimate representative geophysical parameters of each zone, which are then used for field-scale petrophysical analysis. A synthetic study demonstrated how joint inversion of seismic and radar traveltimes and electrical resistance tomography (ERT) data greatly reduces misclassification of zones (down from 21.3% to 3.7%) and improves the accuracy of retrieved zonal parameters (from 1.8% to 0.3%) compared to individual inversions. We applied our scheme to a data set collected in northeastern Switzerland to delineate lithologic subunits within a gravel aquifer. The inversion models resolve three principal subhorizontal units along with some important 3D heterogeneity. Petro-physical analysis of the zonal parameters indicated approximately 30% variation in porosity within the gravel aquifer and an increasing fraction of finer sediments with depth.

Full-waveform inversion of cross-hole ground-penetrating radar data to characterize a gravel aquifer close to the Thur River, Switzerland

Cross-hole radar tomography is a useful tool for mapping shallow subsurface electrical properties viz. dielectric permittivity and electrical conductivity. Common practice is to invert cross-hole radar data with ray-based tomographic algorithms using first arrival traveltimes and first cycle amplitudes. However, the resolution of conventional standard ray-based inversion schemes for cross-hole ground-penetrating radar (GPR) is limited because only a fraction of the information contained in the radar data is used. The resolution can be improved significantly by using a full-waveform inversion that considers the entire waveform, or significant parts thereof. A recently developed 2D time-domain vectorial full-waveform crosshole radar inversion code has been modified in the present study by allowing optimized acquisition setups that reduce the acquisition time and computational costs significantly. This is achieved by minimizing the number of transmitter points and maximizing the number of receiver positions. The improved algorithm was employed to invert cross-hole GPR data acquired within a gravel aquifer (4-10 m depth) in the Thur valley, Switzerland. The simulated traces of the final model obtained by the full-waveform inversion fit the observed traces very well in the lower part of the section and reasonably well in the upper part of the section. Compared to the ray-based inversion, the results from the full-waveform inversion show significantly higher resolution images. At either side, 2.5 m distance away from the cross-hole plane, borehole logs were acquired. There is a good correspondence between the conductivity tomograms and the natural gamma logs at the boundary of the gravel layer and the underlying lacustrine clay deposits. Using existing petrophysical models, the inversion results and neutron-neutron logs are converted to porosity. Without any additional calibration, the values obtained for the converted neutron-neutron logs and permittivity results are very close and similar vertical variations can be observed. The full-waveform inversion provides in both cases additional information about the subsurface. Due to the presence of the water table and associated refracted/reflected waves, the upper traces are not well fitted and the upper 2 m in the permittivity and conductivity tomograms are not reliably reconstructed because the unsaturated zone is not incorporated into the inversion domain.

Constraining 3-D electrical resistance tomography with GPR reflection data forimproved aquifer characterization

Surface-based ground penetrating radar (GPR) and electrical resistance tomography (ERT) are common tools for aquifer characterization, because both methods provide data that are sensitive to hydrogeologically relevant quantities. To retrieve bulk subsurface properties at high resolution, we suggest incorporating structural information derived from GPR reflection data when inverting surface ERT data. This reduces resolution limitations, which might hinder quantitative interpretations. Surface-based GPR reflection and ERT data have been recorded on an exposed gravel bar within a restored section of a previously channelized river in northeastern Switzerland to characterize an underlying gravel aquifer. The GPR reflection data acquired over an area of 240×40 m map the aquifer's thickness and two internal sub-horizontal regions with different depositional patterns. The interface between these two regions and the boundary of the aquifer with then underlying clay are incorporated in an unstructured ERT mesh. Subsequent inversions are performed without applying smoothness constraints across these boundaries. Inversion models obtained by using these structural constraints contain subtle resistivity variations within the aquifer that are hardly visible in standard inversion models as a result of strong vertical smearing in the latter. In the upper aquifer region, with high GPR coherency and horizontal layering, the resistivity is moderately high (N300 Ωm). We suggest that this region consists of sediments that were rearranged during more than a century of channelized flow. In the lower low coherency region, the GPR image reveals fluvial features (e.g., foresets) and generally more heterogeneous deposits. In this region, the resistivity is lower (~200 Ωm), which we attribute to increased amounts of fines in some of the well-sorted fluvial deposits. We also find elongated conductive anomalies that correspond to the location of river embankments that were removed in 2002.