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

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

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

Mechanisms of inorganic carbon-14 attenuation in contaminated groundwater: Effect of solution pH on isotopic exchange and carbonate precipitation reactions

Radioactive 14C is a significant contaminant associated with nuclear fuels and wastes that is potentially highly mobile in the environment as dissolved inorganic carbonate species. This study investigated the mechanisms by which dissolved inorganic 14C is retained in surface and groundwater environments via precipitation and isotopic exchange reactions. Precipitation of calcite in the presence and absence of nucleation sites is considered along with isotopic exchange with both atmospheric CO2 and solid carbonates. Precipitation occurs at calcite supersaturation values of SICAL > 1.5 in the absence of nucleation sites and SICAL > 0–0.5 in the presence of nucleation sites, suggesting that precipitation of 14C-bearing carbonates is much more likely in subsurface environments where nucleation sites are abundant. The maximum 14C removal in solid isotopic exchange experiments occurred after approximately 2 weeks equilibration. In these experiments the amount of 14C removed from solution was proportional to the amount of calcite surface area present, and removal from solution was equivalent to rapid equalisation of the isotope ratio in an 8–10 Å active surface layer. Although the reactivity of natural carbonates may be lower than the calcite samples used in this study, these results suggest isotopic exchange with solids will be an important 14C retardation mechanism in subsurface environments containing only modest TIC concentrations. These results suggest that if inorganic 14C is released into sub-surface environments, both precipitation and solid phase isotopic exchange can result in non-conservative 14C-DIC transport and 14C contamination may persist in groundwater for decades following accidental releases. In contrast, in experiments open to atmosphere with pH values below 9.3, complete loss of dissolved inorganic 14C was very rapid and occurred with timescales of 10's of hours. 14C loss was due to a rapid exchange of dissolved 14C species with 12CO2 (g) and the kinetics of 14C removal increased as pH values were lowered (i.e. atmospheric isotopic exchange was first order with respect to the concentration of carbonic acid present). Thus these results suggest that release of inorganic 14C to surface waters with pH values <9.3 would result in rapid exchange with 12CO2 (g) and 14C would not persist in the aqueous environment, whereas 14C-DIC released to saturated subsurface environments may persist close to the release site for decades due to precipitation and solid phase exchange reactions preventing/retarding transport with the groundwater

Toxicity of Sediment-Associated Pesticides to Chironomus dilutus and Hyalella azteca

Two hundred sediment samples were collected and their toxicity evaluated to aquatic species in a previous study in the agriculturally dominated Central Valley of California, United States. Pyrethroid insecticides were the main contributors to the observed toxicity. However, mortality in approximately one third of the toxic samples could not be explained solely by the presence of pyrethroids in the matrices. Hundreds of pesticides are currently used in the Central Valley of California, but only a few dozen are analyzed in standard environmental monitoring. A significant amount of unexplained sediment toxicity may be due to pesticides that are in widespread use that but have not been routinely monitored in the environment, and even if some of them were, the concentrations harmful to aquatic organisms are unknown. In this study, toxicity thresholds for nine sediment-associated pesticides including abamectin, diazinon, dicofol, fenpropathrin, indoxacarb, methyl parathion, oxyfluorfen, propargite, and pyraclostrobin were established for two aquatic species, the midge Chironomus dilutus and the amphipod Hyalella azteca. For midges, the median lethal concentration (LC50) of the pesticides ranged from 0.18 to 964 μg/g organic carbon (OC), with abamectin being the most toxic and propargite being the least toxic pesticide. A sublethal growth endpoint using average individual ash-free dry mass was also measured for the midges. The no–observable effect concentration values for growth ranged from 0.10 to 633 μg/g OC for the nine pesticides. For the amphipods, fenpropathrin was the most toxic, with an LC50 of 1–2 μg/g OC. Abamectin, diazinon, and methyl parathion were all moderately toxic (LC50s 2.8–26 μg/g OC). Dicofol, indoxacarb, oxyfluorfen, propargite, and pyraclostrobin were all relatively nontoxic, with LC50s greater than the highest concentrations tested. The toxicity information collected in the present study will be helpful in decreasing the frequency of unexplained sediment toxicity in agricultural waterways

Mineralogical attenuation for metallic remediation in a passive system for mine water treatment

Passive systems with constructed wetlands have been consistently used to treat mine water from abandoned mines. Long-term and cost-effective remediation is a crucial expectation for these water treatment facilities. To achieve that, a complex chain of physical, chemical, biological, and mineralogical mechanisms for pollutants removal must be designed to simulate natural attenuation processes. This paper aims to present geochemical and mineralogical data obtained in a recently constructed passive system (from an abandoned mine, Jales, Northern Portugal). It shows the role of different solid materials in the retention of metals and arsenic, observed during the start-up period of the treatment plant. The mineralogical study focused on two types of materials: (1) the ochre-precipitates, formed as waste products from the neutralization process, and (2) the fine-grained minerals contained in the soil of the wetlands. The ochre-precipitates demonstrated to be poorly ordered iron-rich material, which gave rise to hematite upon artificial heating. The heating experiments also provided mineralogical evidence for the presence of an associated amorphous arsenic-rich compound. Chemical analysis on the freshly ochre-precipitates revealed high concentrations of arsenic (51,867 ppm) and metals, such as zinc (1,213 ppm) and manganese (821 ppm), indicating strong enrichment factors relative to the water from which they precipitate. Mineralogical data obtained in the soil of the wetlands indicate that chlorite, illite, chlorite–vermiculite and mica–vermiculite mixedlayers, vermiculite, kaolinite and goethite are concentrated in the fine-grained fractions (<20 and <2 μm). The chemical analyses show that high levels of arsenic (up to 3%) and metals are also retained in these fractions, which may be enhanced by the low degree of order of the clay minerals as suggested by an XRD study. The obtained results suggest that, although the treatment plant has been receiving water only since 2006, future performance will be strongly dependent on these identified mineralogical pollutant hosts.Fundação para a Ciência e a Tecnologia (FCT