RETRASO, a code for modeling reactive transport in saturated and unsaturated porous media

The code RETRASO (REactive TRAnsport of SOlutes) simulates reactive transport of dissolved and gaseous species in non-isothermal saturated or unsaturated problems. Possible chemical reactions include aqueous complexation (including redox reactions), sorption, precipitation-dissolution of minerals and gas dissolution. Various models for sorption of solutes on solids are available, from experimental relationships (linear KD, Freundlich and Langmuir isotherms) to cation exchange and surface complexation models (constant capacitance, diffuse layer and triple layer models). Precipitation-dissolution and aqueous complexation can be modelled in equilibrium or according to kinetic laws. For the numerical solution of the reactive transport equations it uses the Direct Substitution Approach. The use of the code is demonstrated by three examples. The first example models various sorption processes in a smectite barrier. The second example models a complex chemical system in a two dimensional cross-section. The last example models pyrite weathering in an unsaturated medium

A new lattice Boltzmann model for interface reactions between immiscible fluids

In this paper, we describe a lattice Boltzmann model to simulate chemical reactions taking place at the interface between two immiscible fluids. The phase-field approach is used to identify the interface and its orientation, the concentration of reactant at the interface is then calculated iteratively to impose the correct reactive flux condition. The main advantages of the model is that interfaces are considered part of the bulk dynamics with the corrective reactive flux introduced as a source/sink term in the collision step, and, as a consequence, the model’s implementation and performance is independent of the interface geometry and orientation. Results obtained with the proposed model are compared to analytical solution for three different benchmark tests (stationary flat boundary, moving flat boundary and dissolving droplet). We find an excellent agreement between analytical and numerical solutions in all cases. Finally, we present a simulation coupling the Shan Chen multiphase model and the interface reactive model to simulate the dissolution of a collection of immiscible droplets with different sizes rising by buoyancy in a stagnant fluid

Use of groundwater lifetime expectancy for the performance assessment of a deep geologic waste repository: 1. Theory, illustrations, and implications

Long-term solutions for the disposal of toxic wastes usually involve isolation of the wastes in a deep subsurface geologic environment. In the case of spent nuclear fuel, if radionuclide leakage occurs from the engineered barrier, the geological medium represents the ultimate barrier that is relied upon to ensure safety. Consequently, an evaluation of radionuclide travel times from a repository to the biosphere is critically important in a performance assessment analysis. In this study, we develop a travel time framework based on the concept of groundwater lifetime expectancy as a safety indicator. Lifetime expectancy characterizes the time that radionuclides will spend in the subsurface after their release from the repository and prior to discharging into the biosphere. The probability density function of lifetime expectancy is computed throughout the host rock by solving the backward-in-time solute transport adjoint equation subject to a properly posed set of boundary conditions. It can then be used to define optimal repository locations. The risk associated with selected sites can be evaluated by simulating an appropriate contaminant release history. The utility of the method is illustrated by means of analytical and numerical examples, which focus on the effect of fracture networks on the uncertainty of evaluated lifetime expectancy.Comment: 11 pages, 8 figures; Water Resources Research, Vol. 44, 200

Assessing the joint impact of DNAPL source-zone behavior and degradation products on the probabilistic characterization of human health risk

The release of industrial contaminants into the subsurface has led to a rapid degradation of groundwater resources. Contamination caused by Dense Non-Aqueous Phase Liquids (DNAPLs) is particularly severe owing to their limited solubility, slow dissolution and in many cases high toxicity. A greater insight into how the DNAPL source zone behavior and the contaminant release towards the aquifer impact human health risk is crucial for an appropriate risk management. Risk analysis is further complicated by the uncertainty in aquifer properties and contaminant conditions. This study focuses on the impact of the DNAPL release mode on the human health risk propagation along the aquifer under uncertain conditions. Contaminant concentrations released from the source zone are described using a screening approach with a set of parameters representing several scenarios of DNAPL architecture. The uncertainty in the hydraulic properties is systematically accounted for by high-resolution Monte Carlo simulations. We simulate the release and the transport of the chlorinated solvent perchloroethylene and its carcinogenic degradation products in randomly heterogeneous porous media. The human health risk posed by the chemical mixture of these contaminants is characterized by the low-order statistics and the probability density function of common risk metrics. We show that the zone of high risk (hot spot) is independent of the DNAPL mass release mode, and that the risk amplitude is mostly controlled by heterogeneities and by the source zone architecture. The risk is lower and less uncertain when the source zone is formed mostly by ganglia than by pools. We also illustrate how the source zone efficiency (intensity of the water flux crossing the source zone) affects the risk posed by an exposure to the chemical mixture. Results display that high source zone efficiencies are counter-intuitively beneficial, decreasing the risk because of a reduction in the time available for the production of the highly toxic subspecies.Peer ReviewedPostprint (author's final draft

Debates—Stochastic subsurface hydrology from theory to practice: why stochastic modeling has not yet permeated into practitioners?

This is the peer reviewed version of the following article: [Sanchez-Vila, X., and D. Fernàndez-Garcia (2016), Debates—Stochastic subsurface hydrology from theory to practice: Why stochastic modeling has not yet permeated into practitioners?, Water Resour. Res., 52, 9246–9258, doi:10.1002/2016WR019302], which has been published in final form at http://onlinelibrary.wiley.com/doi/10.1002/2016WR019302/abstract. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-ArchivingWe address modern topics of stochastic hydrogeology from their potential relevance to real modeling efforts at the field scale. While the topics of stochastic hydrogeology and numerical modeling have become routine in hydrogeological studies, nondeterministic models have not yet permeated into practitioners. We point out a number of limitations of stochastic modeling when applied to real applications and comment on the reasons why stochastic models fail to become an attractive alternative for practitioners. We specifically separate issues corresponding to flow, conservative transport, and reactive transport. The different topics addressed are emphasis on process modeling, need for upscaling parameters and governing equations, relevance of properly accounting for detailed geological architecture in hydrogeological modeling, and specific challenges of reactive transport. We end up by concluding that the main responsible for nondeterministic models having not yet permeated in industry can be fully attributed to researchers in stochastic hydrogeology.Peer ReviewedPostprint (author's final draft

Influence of Oceanic Forcing on Fate of Nutrients in a Near-Shore-Aquifer

A reactive groundwater transport model has been developed to investigate the fate of nutrients (ammonium, nitrate, and phosphate) in a near-shore coastal aquifer subject to oceanic forcing (tides and waves) and their subsequent discharge to coastal waters. The model is developed by combining the variable-density groundwater flow model SEAWAT-2005 with the reactive multi-component transport model PHT3D v2.10. The influence of tides and waves are typically neglected in prior studies that have examined the transport and transformation of nutrients in coastal aquifers. Oceanic forcing however can induce a highly dynamic surficial salt-freshwater mixing and reaction zone in a near-shore aquifer and this may modify the transport pathways and concentrations of discharging nutrients. The reactions considered in the model include denitrification, nitrification, aerobic degradation of dissolved organic matter, iron oxidation, and phosphate adsorption. The reaction network implemented, including the kinetic rate expressions, has been verified previously by numerical simulations conducted for a near-shore aquifer not exposed to oceanic forcing. The simulations conducted reveal that oceanic forcing significantly modifies the discharge pathway of the groundwater-derived nutrients and the reactions occurring along this pathway. This alters the net production and consumption of nutrients in the near-shore aquifer and their subsequent loading rates to coastal waters. It is further shown that the fate of the nutrients is strongly controlled by the availability of chemical species including dissolved organic matter in seawater recirculating through the near-shore sediments. Moreover, for the conditions simulated, tides led to more intense salt-freshwater mixing in the near-shore aquifer and thus greater transformations of nutrients in the near-shore aquifer compared to regular wave forcing. This study significantly enhances conceptual understanding of the processes controlling the fate of nutrients in a near-shore aquifer and hence provides a valuable tool for improving prediction of nutrient loading rates to coastal waters