Evaporation from a shallow water table: Diurnal dynamics of measured and simulated water and heat regime at the vicinity of a drying sand surface

Accurate estimates of water losses by evaporation from shallow water tables are important for hydrological, agricultural, and climatic purposes. An experiment was conducted in a weighing lysimeter to characterize the diurnal dynamics of evaporation under natural conditions. Sampling revealed a completely dry surface sand layer after 5 days of evaporation. Its thickness was <1 cm early in the morning, increasing to reach 4?5 cm in the evening. This evidence points out fundamental limitations of the approaches that assume hydraulic connectivity from the water table up to the surface, as well as those that suppose monotonic drying when unsteady conditions prevail. The computed vapor phase diffusion rates from the apparent drying front based on Fick's law failed to reproduce the measured cumulative evaporation during the sampling day. We propose that two processes rule natural evaporation resulting from daily fluctuations of climatic variables: (i) evaporation of water, stored during nighttime due to redistribution and vapor condensation, directly into the atmosphere from the soil surface during the early morning hours, that could be simulated using a mass transfer approach and (ii) subsurface evaporation limited by Fickian diffusion, afterward. For the conditions prevailing during the sampling day, the amount of water stored at the vicinity of the soil surface was 0.3 mm and was depleted before 11:00. Combining evaporation from the surface before 11:00 and subsurface evaporation limited by Fickian diffusion after that time, the agreement between the estimated and measured cumulative evaporation was significantly improved.

3D crosshole ERT for aquifer characterization and monitoring of infiltrating river water

The hydrogeological properties and responses of a productive aquifer in northeastern Switzerland are investigated. For this purpose, 3D crosshole electrical resistivity tomography (ERT) is used to define the main lithological structures within the aquifer (through static inversion) and to monitor the water infiltration from an adjacent river. During precipitation events and subsequent river flooding, the river water resistivity increases. As a consequence, the electrical characteristics of the infiltrating water can be used as a natural tracer to delineate preferential flow paths and flow velocities. The focus is primarily on the experiment installation, data collection strategy, and the structural characterization of the site and a brief overview of the ERT monitoring results. The monitoring system comprises 18 boreholes each equipped with 10 electrodes straddling the entire thickness of the gravel aquifer. A multi-channel resistivity system programmed to cycle through various four-point electrode configurations of the 180 electrodes in a rolling sequence allows for the measurement of approximately 15,500 apparent resistivity values every 7 h on a continuous basis. The 3D static ERT inversion of data acquired under stable hydrological conditions provides a base model for future time-lapse inversion studies and the means to investigate the resolving capability of our acquisition scheme. In particular, it enables definition of the main lithological structures within the aquifer. The final ERT static model delineates a relatively high-resistivity, low-porosity, intermediate-depth layer throughout the investigated aquifer volume that is consistent with results from well logging and seismic and radar tomography models. The next step will be to define and implement an appropriate time-lapse ERT inversion scheme using the river water as a natural tracer. The main challenge will be to separate the superposed time-varying effects of water table height, temperature, and salinity variations associated with the infiltrating water.

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.

Biogeochemical evolution of a marine shore tailings deposit during bioremediation

Summary : Mining activities produce enormous amounts of waste material known as tailings which are composed of fine to medium size particles. These tailings often contain sulfides, which oxidation can lead to acid and metal contamination of water; therefore they need to be remediated. In this work a tailings bioremediation approach was investigated by an interdisciplinary study including geochemistry, mineralogy and microbiology. The aim of the work was to study the effect of the implementation of wetland above oxidizing tailings on the hydrogeology and the biogeochemical element cycles, and to assess the system evolution over time. To reach these goals, biogeochemical processes occurring in a marine shore tailings deposit were investigated. The studied tailings deposit is located at the Bahìa de Ite, Pacific Ocean, southern Peru, where between 1940 and 1996 the tailings were discharged from the two porphyry copper mines Cuajone and Toquepala. After the end of deposition, a remediation approach was initiated in 1997 with a wetland implementation above the oxidizing tailings. Around 90% of the tailings deposits (total 16 km2) were thus remediated, except the central delta area and some areas close to the shoreline. The multi-stable isotope study showed that the tailings were saturated with fresh water in spite of the marine setting, due to the high hydraulic gradient resulting from the wetland implementation. Submarine groundwater discharge (SGD) was the major source of SO4 2-, C1-, Na+, Fe2+, and Mn2+ input into the tailings at the original shelf-seawater interface. The geochemical study (aquatic geochemistry and X-Ray diffraction (XRD) and sequential extractions from the solid fraction) showed that iron and sulfur oxidation were the main processes in the non-remediated tailings, which showed a top a low-pH oxidation zone with strong accumulation of efflorescent salts at the surface due to capillary upward transport of heavy metals (Fe, Cu, Zn, Mn, Cd, Co, and Ni) in the arid climate. The study showed also that the implementation of the wetland resulted in very low concentrations of heavy metals in solution (mainly under the detection limit) due to the near neutral pH and more reducing conditions (100-150 mV). The heavy metals, which were taken from solution, precipitated as hydroxides and sulfides or were bound to organic matter. The bacterial community composition analysis by Terminal Restriction Fragment Length Polymorphism (T-RFLP) and cloning and sequencing of 16S rRNA genes combined with a detailed statistical analysis revealed a high correlation between the bacterial distribution and the geochemical variables. Acidophilic autotrophic oxidizing bacteria were dominating the oxidizing tailings, whereas neutrophilic and heterotrophic reducing bacteria were driving the biogeochemical processes in the remediated tailings below the wetland. At the subsurface of the remediated tailings, an iron cycling was highlighted with oxidation and reduction processes due to micro-aerophilic niches provided by the plant rhizosphere in this overall reducing environment. The in situ bioremediation experiment showed that the main parameter to take into account for the effectiveness was the water table and chemistry which controls the system. The constructed remediation cells were more efficient and rapid in metal removal when saturation conditions were available. This study showed that the bioremediation by wetland implementation could be an effective and rapid treatment for some sulfidic mine tailings deposits. However, the water saturation of the tailings has to be managed on a long-term basis in order to guarantee stability. Résumé : L'activité minière produit d'énormes quantités de déchets géologiques connus sous le nom de « tailings » composées de particules de taille fine à moyenne. Ces déchets contiennent souvent des sulfures dont l'oxydation conduit à la formation d'effluents acides contaminés en métaux, d'où la nécessité d'effectuer une remédiation des sites de stockage concernés. Le but de ce travail est dans un premier temps d'étudier l'effet de la bio-remédiation d'un dépôt de tailings oxydés sur l'hydrogéologie du système et les cycles biogéochimiques des éléments et en second lieu, d'évaluer l'évolution du processus de remédiation dans le temps. Le site étudié dans ce travail est situé dans la Bahía de Ite, au sud du Pérou, au bord de l'Océan Pacifique. Les déchets miniers en question sont déposés dans un environnement marin. De 1940 à 1996, les déchets de deux mines de porphyre cuprifère - Cuajone et Toquepala - ont été acheminés sur le site via la rivière Locumba. En 1997, une première remédiation a été initiée avec la construction d'une zone humide sur les tailings. Depuis, environ 90% de la surface du dépôt (16 km2) a été traité, les parties restantes étant la zone centrale du delta du Locumba et certaines zones proches de la plage. Malgré la proximité de l'océan, les études isotopiques menées dans le cadre de ce travail ont montré que les tailings étaient saturés en eau douce. Cette saturation est due à la pression hydraulique résultant de la mise en place des zones humides. Un écoulement d'eau souterrain sous-marin a été à détecté à l'interface entre les résidus et l'ancien fond marin. En raison de la géologie locale, il constitue une source d'entrée de SO4 2-, Cl-, Na+, FeZ+, et Mn2+ dans le système. L'analyse de la géochimie aquatique, la Diffraction aux Rayons X (XRD) et l'extraction séquentielle ont montré que l'oxydation du fer et .des sulfures est le principal processus se produisant dans les déchets non remédiés. Ceci a entraîné le développement d'une zone d'oxydation à pH bas induisant une forte accumulation des sels efflorescents, conséquence de la migration capillaire des métaux lourds (Fe, Cu, Zn, Mn, Cd, Co et Ni) de la solution vers la surface dans ce climat aride. Cette étude a montré également que la construction de la zone humide a eu comme résultats une précipitation des métaux dans des phases minérales en raison du pH neutre et des conditions réductrices (100-150mV). Les métaux lourds ont précipité sous la forme d'hydroxydes et de sulfures ou sont adsorbés à la matière organique. L'analyse de la composition de la communauté bactérienne à l'aide la technique T-RFLP (Terminal Restriction Fragment Length Polymorphism) et par le clonage/séquençage des gènes de l'ARNr 16S a été combinée à une statistique détaillée. Cette dernière a révélé une forte corrélation entre la distribution de bactéries spécifiques et la géochimie : Les bactéries autotrophes acidophiles dominent dans les déchets oxydés non remédiés, tandis que des bactéries hétérotrophes neutrophiles ont mené les processus microbiens dans les déchets remédiés sous la zone humide. Sous la surface de la zone humide, nos analyses ont également mis en évidence un cycle du fer par des processus d'oxydoréduction rendus possibles par la présence de niches micro-aérées par la rhizosphère dans cet environnement réducteur. L'expérience de bio-remédiation in situ a montré que les paramètres clés qui contrôlent l'efficacité du traitement sont le niveau de la nappe aquifère et la chimie de l'eau. Les cellules de remédiation se sont montrées plus efficaces et plus rapides lorsque le système a pu être saturé en eau. Finalement, cette étude a montré que la bio-remédiation de déchets miniers par la construction de zones humides est un moyen de traitement efficace, rapide et peu coûteux. Cependant, la saturation en eau du système doit être gérée sur le long terme afin de garantir la stabilité de l'ensemble du système.

Self-potential investigations of a gravel bar in a restored river corridor

Self-potentials (SP) are sensitive to water fluxes and concentration gradients in both saturated and unsaturated geological media, but quantitative interpretations of SP field data may often be hindered by the superposition of different source contributions and time-varying electrode potentials. Self-potential mapping and close to two months of SP monitoring on a gravel bar were performed to investigate the origins of SP signals at a restored river section of the Thur River in northeastern Switzerland. The SP mapping and subsequent inversion of the data indicate that the SP sources are mainly located in the upper few meters in regions of soil cover rather than bare gravel. Wavelet analyses of the time-series indicate a strong, but non-linear influence of water table and water content variations, as well as rainfall intensity on the recorded SP signals. Modeling of the SP response with respect to an increase in the water table elevation and precipitation indicate that the distribution of soil properties in the vadose zone has a very strong influence. We conclude that the observed SP responses on the gravel bar are more complicated than previously proposed semi-empiric relationships between SP signals and hydraulic head or the thickness of the vadose zone. We suggest that future SP monitoring in restored river corridors should either focus on quantifying vadose zone processes by installing vertical profiles of closely spaced SP electrodes or by installing the electrodes within the river to avoid signals arising from vadose zone processes and time-varying electrochemical conditions in the vicinity of the electrodes.

The role of soil in vegetated gravelly river braid plains: more than just a passive response?

This paper reviews the role of alluvial soils in vegetated gravelly river braid plains. When considering decadal time scales of river evolution, we argue that it becomes vital to consider soil development as an emergent property of the developing ecosystem. Soil processes have been relatively overlooked in accounts of the interactions between braided river processes and vegetation, although soils have been observed on vegetated fluvial landforms. We hypothesise that soil development plays a major role in the transition (speed and pathway) from a fresh sediment deposit to a vegetated soil-covered landform. Disturbance (erosion and/or deposition), vertical sediment structure (process history), vegetation succession, biological activity and water table fluctuation are seen as the main controls on early alluvial soil evolution. Erosion and deposition processes may not only act as soil disturbing agents, but also as suppliers of ecosystem resources, because of their role in delivering and changing access (e.g. through avulsion) to fluxes of water, fine sediments and organic matter. In turn, the associated initial ecosystem may influence further fluvial landform development, such as through the trapping of fine-grained sediments (e.g. sand) by the engineering action of vegetation and the deposit stabilisation by the developing above and belowground biomass. This may create a strong feedback between geomorphological processes, vegetation succession and soil evolution which we summarise in a conceptual model. We illustrate this model by an example from the Allondon River (CH) and identify the research questions that follow.

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.

Inversion of multiple intersecting high-resolution crosshole GPR profiles for hydrological characterization at the Boise Hydrogeophysical Research Site

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. In particular, crosshole ground-penetrating radar (GPR) tomography has shown much promise in hydrology because of its ability to provide highly detailed images of subsurface radar wave velocity, which is strongly linked to soil water content. Here, we develop and demonstrate a procedure for inverting together multiple crosshole GPR data sets in order to characterize the spatial distribution of radar wave velocity below the water table at the Boise Hydrogeophysical Research Site (BHRS) near Boise, Idaho, USA. Specifically, we jointly invert 31 intersecting crosshole GPR profiles to obtain a highly resolved and consistent radar velocity model along the various profile directions. The model is found to be strongly correlated with complementary neutron porosity-log data and is further corroborated by larger-scale structural information at the BHRS. This work is an important prerequisite to using crosshole GPR data together with existing hydrological measurements for improved groundwater flow and contaminant transport modeling.

European climate variations over the past half-millennium reconstructed from groundwater

Temperature reconstructions for recent centuries are the basis of estimations of the natural variability in the climate system before and during the onset of anthropogenic perturbation. Here we present, for the first time, an independent and physically based reconstruction of mean annual temperature over the past half millennium obtained from groundwater in France. The reconstructed noble gas temperature (NGT) record suggests cooler than present climate conditions throughout the 16th-19th centuries. Periods of warming occur in the 17th-18th and 20th century, while cooling is reconstructed in the 19th century. A noticeable coincidence with other temperature records is demonstrated. Deuterium excess varies in parallel with the NGT, and indicates variation in the seasonality of the aquifer recharge; whereas high excess air in groundwater indicates periods with high oscillations of the water table.

Usages de l'eau à Gaya (Niger) : entre fortes potentialités et contraintes majeures

Résumé: Le département de Gaya, cadre de notre étude, est situé au sud-ouest de la république du Niger. Il dispose d'un important potentiel hydrique composé des eaux de surface (une centaine de mares permanentes, le fleuve Niger sur 106 km) et de sept aquifères superposés comprenant des nappes de subsurface (affleurantes par endroit) et des nappes artésiennes. L'étude sur les usages de l'eau à Gaya a été menée à travers plusieurs axes centrés sur l'estimation et la répartition spatiale des ressources en eau, le cadre juridique et institutionnel régulant leur mise en valeur, les différents secteurs d'utilisation de l'eau ainsi que les contraintes affectant cette utilisation. L'usage de la cartographie à travers les SIG dans le traitement et l'analyse des données, couplée à notre expérience d'une dizaine d'année de travaux dans la région, a permis de dresser des synthèses richement illustrées permettant de mieux comprendre tous les enjeux liés à la problématique des usages de l'eau dans cette partie du Niger. Contrairement à la vision que l'on a traditionnellement du Sahel où le manque d'eau constitue une des contraintes majeures au développement, ici des conditions locales particulières contredisent ce cliché et transposent le débat sur un autre plan. Il s'agit de la maîtrise de l'eau au niveau local à travers l'élaboration d'une politique appropriée qui tienne compte non seulement des spécificités locales de la ressource, mais aussi des différents types d'usages. La politique de l'eau au Niger, définie selon le Schéma directeur de mise en valeur et de gestion des ressources en eau, à travers la mise en place d'un important arsenal juridique et institutionnel, a eu le mérite de tracer un canevas sur la question, mais a montré ses limites au niveau pratique après dix ans d'essai. En effet au niveau de Gaya, ni l'Etat ni les partenaires au développement (bailleurs de fonds extérieurs) n'ont tenu compte des caractéristiques locales de la ressource ou du contexte socioéconomique particulier de la région. Ce qui a entraîné la réalisation d'infrastructures inadaptées aux réalités hydrogéologiques locales ainsi que des choix inappropriés au niveau de certains aménagements. En dépit de l'abondance de la ressource, son accès tant au niveau quantitatif que qualitatif, reste difficile pour une grande partie des acteurs ruraux. Les différents handicaps rencontrés dans la mise en valeur des ressources en eau résultent de cette incohérence de la politique nationale de l'eau, mais aussi de la difficulté de son application sur le terrain où persiste un pluralisme juridique caractérisé par la cohabitation de deux systèmes de régulation à savoir les droits coutumiers et la législation moderne. Ces différents éléments mis en évidence dans cette étude sur la zone de Gaya pourraient servir de base pour un meilleur aménagement des ressources en eau dans le cadre plus large d'une politique d'aménagement du territoire prenant en compte tous les facteurs tant physiques que socioéconomiques de la région. Abstract: The department of Gaya, in which this study was done, is located in the SW area of the Republic of Niger. It has an important hydrological potential composed of surface water (approximately 100 permanent ponds, 106 km of the Niger River) and 7 bodies of underground water sources including sub-surface and artresan wells. This study of the exploitation of wtaer in Gaya has been carried out employing several parameters based on: the estimation and spacial distribution of water ressources, the juridic and institutional rules governing their utilisation and the various constraints affecting this exploitation. The use of mapping when treating and analysing data, coupled with ten years personel field experience, resulted in a richly illustrated synthesis of this data. This, in turn, led to a better comprehension of all the factors related to problems of water utilisation in this particular region of Niger. Contrary to the generally accepted view that the lack of water ressources is a major contributing factor to the lack of development in the Sahel, in Gaya the local conditions contradict this statement. In this region, and at the local level, the proper use of water is based on the elaboration of an appropriate policy which takes into account not only the local specifics of water ressources but the various types of water utilsation as well. Local use of water and water ressources are dependant on established rules. Water policy in Niger is defined by the General Schema based on an important institutional and judicary arsenal of rules and regulations. However, after a ten-year trial period, this system was shown to have its limitations. In Gaya, neither the State nor the development agencies took into consideration local characteristics nor the socio-economic context of the region. This, in turn, resulted in putting in place infrastructures that were not adapted to local hydrogeological realities as well as inappropriate choices in land planning and development. In spite of the abundance of water ressources, access to them remains difficult for most of the rural population. The various difficulties encountered are the result of incoherent water policies on a national level as well as the lack of practical application in this area. This is due to a double judicary system where two regulatory systems co-exist:traditional laws and modern legislation. the different elements brought out by this study could serve as a basis for a better utilisation of water ressources on a larger scale in which land planning and development policies would take into consideration all the physcial as well as the socio-economical factors of this region.