Search

Contaminants like nitroaromatic compounds can be degraded in the subsurface to similar or even more toxic metabolites. Degradation or transformation rates are dependent on physical, chemical and biological properties which can be different in sedimentological layers or other heterogeneous structures of aquifers. Sediments with low hydraulic conductivity can even consist of immobile water. These regions are only accessible by diffusion. Most modelling approaches accounting for immobile water regions focused on the mathematical description of the transport and decay of the parent compound. The objective of this study was to develop an analytical model to quantify the transport and formation of a metabolite in dual-porosity media describing the exchange between mobile and immobile water regions based on the metabolite's diffusion coefficient. Column experiments with a well-defined immobile water region were performed under anoxic conditions at three different water flow velocities. The model compound 4-Cl-nitrobenzene was reduced to 4-Cl-aniline (4-Cl-An) by surface-bound Fe (II) species within the immobile water region. Transport and formation of the metabolite were quantified with a modified solution of the single fissure dispersion model assuming additionally for the region with immobile water first-order metabolite production, irreversible sorption and an instantaneous equilibrium sorption. The number of unknown fitting parameters was reduced to two (sorption rate and retardation factor) by stepwise parameter estimation using tracer and parent compound data. Experimental results of the metabolite for each water flow velocity were successfully described with a first-order production term (λ prod = 1.51 ± 0.08 h-1), retardation factor (R im = 2.94 ± 0.45) and first-order irreversible sorption rate (K im = 0.39 ± 0.16 h-1) within the immobile water region. Model results supported that 4-Cl-An was formed within the immobile water region. 4-Cl-An sorbed instantaneously onto the clay matrix while a fraction was irreversibly sorbed. Experimental results and the provided analytical solution help to improve the understanding about reactive transport and the formation of metabolites in dual-porosity media.

Variable density flow in porous media has been studied primarily using numerical models because it is a semi-chaotic and transient process. Most of these studies have been 2D, owing to the computational restrictions on 3D simulations, and the ability to observe variable density flow in 2D experimentation. However, it is recognised that variable density flow is a three-dimensional process. A 3D system may cause weaker variable density flow than a 2D system due to stronger dispersion, but may also result in bigger fingers and hence stronger variable density flow because of more space for fingers to coalesce. This study aimed to determine the representativeness of 2D modelling to simulate 3D variable density flow. 3D homogeneous sand column experiments were conducted at three different water flow velocities with three different bromide tracer solutions mixed with methanol resulting in different density ratios. Both 2D axisymmetric and 3D numerical simulations were performed to reproduce experimental data. Experimental results showed that the magnitude of variable density flow increases with decreasing flow rates and decreasing density ratios. The shapes of the observed breakthrough curves differed significantly from those produced by 2D axisymmetric and 3D simulations. Compared to 2D simulations, the onset of instabilities was delayed but the growth was more pronounced in 3D simulations. Despite this difference, both 2D axisymmetric and 3D models successfully simulated mass recovery with high efficiency (between 77% and 99%). This study indicates that 2D simulations are sufficient to understand integrated features of variable density flow in homogeneous sand column experiments.

Nitroaromatic compounds (NACs) are reduced by structural or surface bound Fe (II) species under anaerobic conditions in the subsurface. This reaction preferentially occurs on clay minerals which are mainly present in areas with low hydraulic conductivity containing nearly immobile water. Diffusion is the dominating transport process in these zones. Due to the complexity in such heterogeneous systems, the mathematical prediction of reactive solute transport taking into account diffusive mass exchange into immobile water regions still remains challenging. Therefore, the influence of immobile water regions on the fate of 4-Cl-Nitrobenzene (4-Cl-Nb) was quantified in dual-porosity column experiments at three different mean transit times under saturated anaerobic conditions in the presence of soluble Fe (II). A multi-tracer approach and a Single Fissure Dispersion Model (SFDM) were used to estimate input parameter to further model the transport of 4-Cl-Nb. Reactive solute transport of 4-Cl-Nb was quantified considering instantaneous sorption on to the clay matrix and a reduction within the immobile water region following first-order kinetics. The experimental results indicated that sorption onto the clay matrix enhanced the mass exchange of 4-Cl-Nb into immobile water region compared to nonreactive solutes. At the same time the abiotic reduction of 4-Cl-Nb limited the process of back diffusion to mobile water regions. Fitted retardation factors (Rim = 4.62 ± 0.68) and decay rates (k = 1.51 ± 0.08 h− 1) were independent on tested flow velocities. Findings of this study can advance the understanding on the fate of NACs in the subsurface which is essential for prediction of reactive solute transport at field scale.

Phosphorous (P) leaching (e.g., from effluents, fertilizers) and transport in highly permeable subsurface media can be an important pathway that contributes to eutrophication of receiving surface waters as groundwater recharges the base-flow of surface waters. Here we investigated attenuation and transport of orthophosphate-P in gravel aquifer and vadose zone media in the presence and absence of model colloids (Escherichia coli, kaolinite, goethite). Experiments were conducted using repacked aquifer media in a large column (2m long, 0.19m in diameter) and intact cores (0.4m long, 0.24m in diameter) of vadose zone media under typical field flow rates. In the absence of the model colloids, P was readily traveled through the aquifer media with little attenuation (up to 100% recovery) and retardation, and P adsorption was highly reversible. Conversely, addition of the model colloids generally resulted in reduced P concentration and mass recovery (down to 28% recovery), and increased retardation and adsorption irreversibility in both aquifer and vadose zone media. The degree of colloid-assisted P attenuation was most significant in the presence of fine material and Fe-containing colloids at low flow rate but was least significant in the presence of coarse gravels and E. coli at high flow rate. Based on the experimental results, setback distances of 49-53m were estimated to allow a reduction of P concentrations in groundwater to acceptable levels in the receiving water. These estimates were consistent with field observations in the same aquifer media. Colloid-assisted P attenuation can be utilized to develop mitigation strategies to better manage effluent applications in gravelly soils. To efficiently retain P within soil matrix and reduce P leaching to groundwater, it is recommended to select soils that are rich in iron oxides, to periodically disturb soil preferential flow paths by tillage, and to apply a low irrigation rate.

Diffusive mass exchange into immobile water regions within heterogeneous porous aquifers influences the fate of solutes. The percentage of immobile water is often unidentified in natural aquifers though. Hence, the mathematical prediction of solute transport in such heterogeneous aquifers remains challenging. The objective of this study was to find a simple analytical model approach that allows quantifying properties of mobile and immobile water regions and the portion of immobile water in a porous system. Therefore, the Single Fissure Dispersion Model (SFDM), which takes into account diffusive mass exchange between mobile and immobile water zones, was applied to model transport in well-defined saturated dual-porosity column experiments. Direct and indirect model validation was performed by running experiments at different flow velocities and using conservative tracer with different molecular diffusion coefficients. In another column setup, immobile water regions were randomly distributed to test the model applicability and to determine the portion of immobile water. In all setups, the tracer concentration curves showed differences in normalized maximum peak concentration, tailing and mass recovery according to their diffusion coefficients. These findings were more pronounced at lower flow rates (larger flow times) indicating the dependency of diffusive mass exchange into immobile water regions on tracers' molecular diffusion coefficients. The SFDM simulated all data with high model efficiency. Successful model validation supported the physical meaning of fitted model parameters. This study showed that the SFDM, developed for fissured aquifers, is applicable in porous media and can be used to determine porosity and volume of regions with immobile water. Copyright © 2015 John Wiley & Sons, Ltd.

Fate and transport of solutes in heterogeneous porous media is largely affected by diffusive mass exchange between mobile and immobile water zones. Since it is difficult to directly measure and determine the effect in the aquifers , multi-tracer experiments in combination with mathematical modeling are used to obtain quantitative information about unknown system parameters such as the effective mobile and immobile porosity, and the diffusive mass exchange between mobile and immobile water zones. The Single Fissure Dispersion Model (SFDM) describing nonreactive transport of solutes in saturated dual-porosity media, has been employed as a modeling approach to explain dual-porosity experiments in the field and laboratory (column experiments). SFDM optimization with conventional methods of minimization was immensely difficult due to its complex analytical form. Thus, previous studies used a trial and error procedure to fit it to the experimental observations. In this study, a rigorous optimization technique based on the newly developed scatter search method is presented that automatically minimizes the SFDM to find the optimal values of the hydrogeologically related parameters. The new program (OptSFDM) is accompanied with an easy-to-use graphical user interface (GUI) that is flexible and fully integrated. The program usability is showcased by a few, previously presented experimental case studies, and compared against the currently available, trial-and-error based, command-line executable, SFDM code.