3D-electrical resistivity tomography monitoring of salt transport in homogeneous and layered soil samples

Monitoring transport of dissolved substances in soil deposits is particularly relevant where safety is concerned, as in the case of geo-environmental barriers. Geophysical methods are very appealing, since they cover a wide domain, localising possible preferential flow paths and providing reliable links between geophysical quantities and hydrological variables. This paper describes a 3D laboratory application of Electrical Resistivity Tomography (ERT) used to monitor solute transport processes. Dissolution and transport tests on both homogeneous and heterogeneous samples were conducted in an instrumented oedometer cell. ERT was used to create maps of electrical conductivity of the monitored domain at different time intervals and to estimate concentration variations within the interstitial fluid. Comparisons with finite element simulations of the transport processes were performed to check the consistency of the results. Tests confirmed that the technique can monitor salt transport, infer the hydro-chemical behaviour of heterogeneous geomaterials and evaluate the performances of clay barrier

Inverse chemical modeling and radiocarbon dating of palaeogroundwaters: The Tertiairy Ledo-paniselian aquifer in Flanders, Belgium.

Groundwater samples from the Ledo-Paniselian aquifer have been interpreted for chemical reaction patterns

Verification and intercomparison of reactive transport codes to describe root-uptake

Several mathematical models have been developed to simulate processes and interactions in the plant rhizosphere. Most of these models are based on a rather simplified description of the soil chemistry and interactions of plant roots in the rhizosphere. In particular the feedback loops between exudation, water and solute uptake are mostly not considered, although their importance in the bioavailability of mineral elements for plants has been demonstrated. The aim of this work was to evaluate three existing coupled speciation-transport tools to model rhizosphere processes. In the field of hydrogeochemistry, such␣computational tools have been developed to␣describe acid-base and redox reactions, complexation and ion exchange, adsorption and precipitation of chemical species in soils and aquifers using thermodynamic and kinetic relationships. We implemented and tested a simple rhizosphere model with three geochemical computational tools (ORCHESTRA, MIN3P, and PHREEQC). The first step was an accuracy analysis of the different solution strategies by comparing the numerical results to the analytical solution of solute uptake (K or Ca) by a single cylindrical root. All models are able to reproduce the concentration profiles as well as the uptake flux. The relative error of the simulated concentration profile decreases with increasing distance from the root. The uptake flux was simulated for all codes with less than 5% error for K and less than 0.4% for Ca. The strength of the codes presented in this paper is that they can also be used to investigate more complex and coupled biogeochemical processes in rhizosphere models. This is shown exemplarily with simulations involving both exudation and uptake and the simultaneous uptake of solute and wate

Rapport final du projet européen CatClay sur les processus de migration des cations dans les roches argileuses indurées

International audienceIn the framework of the feasibility studies on the radioactive waste disposal in deep argillaceous formations, it isnow well established that the transport properties of solutes in clay rocks, i.e. parameter values for Fick’s law, are mainlygoverned by the negatively charged clay mineral surface. While a good understanding of the diffusive behaviour of non-reactiveanionic and neutral species is now achieved, much effort has to be placed on improving understanding of coupledsorption/diffusion phenomena for sorbing cations. Indeed, several cations known to form highly stable surface complexes withsites on mineral surfaces migrate more deeply into clay rock than expected. Therefore, the overall objective of the EC CatClayproject is to address this issue, using a ‘bottom-up’ approach, in which simpler, analogous systems (here a compacted clay,‘pure’ illite) are experimentally studied and modelled, and then the transferability of these results to more complex materials, i.e.the clay rocks under consideration in France, Switzerland and Belgium for hosting radioactive waste disposal facilities, isverified. The cations of interest were chosen for covering a representative range of cations families: from a moderately sorbingcation, the strontium, to three strongly sorbing cations, Co(II), Zn(II) and Eu(III). For the 4 years of this project, much effort wasdevoted to developing and applying specific experimental methods needed for acquiring the high precision, reliable data neededto test the alternative hypotheses represented by different conceptual-numerical models. The enhanced diffusion of the sorbingcations of interest was confirmed both in the simpler analogous illite system for Sr2+, Co(II) and Zn(II), but also in the naturalclay rocks, except for Eu(III). First modelling approach including diffusion in the diffuse double layer (DDL) promisinglysucceeded in reproducing the experimental data under the various conditions both in illite and clay rocks, even though someassumptions made have to be verified. In parallel, actual 3D geometrical pore size distributions of compacted illite, and in lessextent, clay rock samples, were successfully determined by combining TEM and FIB-nt analyses on materials maintained in awater-like saturation state by means of an extensive impregnation step. Based on this spatial distribution of pores, first numericaldiffusion experiments were carried at the pore scale through virtual illite, enabling a better understanding of how transferpathways are organized in the porous media. Finally, the EC CatClay project allowed a better understanding of the migration ofstrongly sorbing tracers through low permeability ‘clay rock’ formations, increasing confidence in our capacity to demonstratethat the models used to predict radionuclide migration through these rocks are scientifically sound

A global method for coupling transport with chemistry in heterogeneous porous media

Modeling reactive transport in porous media, using a local chemical equilibrium assumption, leads to a system of advection-diffusion PDE's coupled with algebraic equations. When solving this coupled system, the algebraic equations have to be solved at each grid point for each chemical species and at each time step. This leads to a coupled non-linear system. In this paper a global solution approach that enables to keep the software codes for transport and chemistry distinct is proposed. The method applies the Newton-Krylov framework to the formulation for reactive transport used in operator splitting. The method is formulated in terms of total mobile and total fixed concentrations and uses the chemical solver as a black box, as it only requires that on be able to solve chemical equilibrium problems (and compute derivatives), without having to know the solution method. An additional advantage of the Newton-Krylov method is that the Jacobian is only needed as an operator in a Jacobian matrix times vector product. The proposed method is tested on the MoMaS reactive transport benchmark.Comment: Computational Geosciences (2009) http://www.springerlink.com/content/933p55085742m203/?p=db14bb8c399b49979ba8389a3cae1b0f&pi=1

Dissolved noble gases and stable isotopes as tracers of preferential fluid flow along faults in the Lower Rhine Embayment, Germany

Groundwater in shallow unconsolidated sedimentary aquifers close to the Bornheim fault in the Lower Rhine Embayment (LRE), Germany, has relatively low δ2H and δ18O values in comparison to regional modern groundwater recharge, and 4He concentrations up to 1.7 × 10−4 cm3 (STP) g–1 ± 2.2 % which is approximately four orders of magnitude higher than expected due to solubility equilibrium with the atmosphere. Groundwater age dating based on estimated in situ production and terrigenic flux of helium provides a groundwater residence time of ∼107 years. Although fluid exchange between the deep basal aquifer system and the upper aquifer layers is generally impeded by confining clay layers and lignite, this study’s geochemical data suggest, for the first time, that deep circulating fluids penetrate shallow aquifers in the locality of fault zones, implying  that sub-vertical fluid flow occurs along faults in the LRE. However, large hydraulic-head gradients observed across many faults suggest that they act as barriers to lateral groundwater flow. Therefore, the geochemical data reported here also substantiate a conduit-barrier model of fault-zone hydrogeology in unconsolidated sedimentary deposits, as well as corroborating the concept that faults in unconsolidated aquifer systems can act as loci for hydraulic connectivity between deep and shallow aquifers. The implications of fluid flow along faults in sedimentary basins worldwide are far reaching and of particular concern for carbon capture and storage (CCS) programmes, impacts of deep shale gas recovery for shallow groundwater aquifers, and nuclear waste storage sites where fault zones could act as potential leakage pathways for hazardous fluids