Evidence of compound-dependent hydrodynamic and mechanical transverse dispersion by multitracer laboratory experiments

Evolution of plume geometry, dilution and reactive mixing in porous media under highly transient flow fields at the surface water-groundwater interface

Highly transient boundary conditions affect mixing of dissolved solutes in groundwater. An example of these transient boundary conditions occurs at the surface water-groundwater interface, where the water level in rivers can change rapidly due to the operation of hydropower plants, leading to a regime known as hydropeaking. Inspired by this phenomenon, this work studies at laboratory scale the effects of fluctuating surface water bodies on solute transport in aquifers. We performed flow-through experiments at two different flow velocities and under steady and transient flow conditions where a conservative tracer was injected in the system and its concentration measured with optical imaging methods. The experimental results were quantitatively interpreted with numerical simulations implementing a non-linear velocity-dependent dispersive transport model. We estimated plume dilution by computing the dilution index and its evolution as well as two key geometrical metrics of the transient plumes: the perimeter and the area. We further investigated reactive mixing and mixing enhancement considering mixing-controlled bimolecular reactions using different critical mixing ratios. In general, highly transient boundary conditions lead to a larger area, perimeter and plume dilution and the results show greater relative enhancement for the scenarios with low groundwater flow velocity. A linear relationship was observed between the evolution of the area and the dilution index of the plumes for the transient flow scenarios investigated. Considering reactive transport and mixing-limited reactions at the surface water-groundwater interface, we identified different dilution and reaction dominated regimes, characterized, respectively, by increasing and decreasing plume entropies at different mixing ratios of the reactants. Furthermore, reactive mixing was enhanced by transient flows leading to a faster degradation of contaminant plumes compared to corresponding steady flow conditions.

Integration of groundwater vulnerability with contaminants transport modeling in unsaturated zone, case study El-Sharqia, Egypt

Nowadays, irrigation uses large amount of marginal wastewater due to continuous decline in fresh water supply. As a consequence, using this wastewater for different purposes can cause some adverse environmental impacts. Anthropogenic activities such as septic tanks, sewage ponds, and polluted drains have large influence on deterioration of shallow groundwater aquifers. So, construction of many wastewater treatment plants in these areas is mandatory to control and mitigate this deterioration. Groundwater vulnerability assessment maps and contamination simulation in unsaturated zone can be beneficial in understanding contaminants pathways and groundwater quality evolution. This work is mainly focused on aquifer vulnerability assessment to pollution and the role of vadose zone in attenuation of contaminants transport through it prior to groundwater seepage. Therefore, about 56 drainage and groundwater samples were collected and analyzed for potentially toxic elements. The most vulnerable sector was determined using GOD method revealing that central parts of the study area are the most threatened zones with some scattered sporadic zone of sensitivity to pollution and this was verified through the zonation of Pb, Fe, and Mn spatial concentrations. The leakage of these elements through the unsaturated zone was further simulated using HYDRUS-1D model for the next 10-year period to determine the extent of the pollution plumes and maximum concentration of these elements that percolate to the groundwater directly. The concentration of Fe, Pb, and Mn at the end of the simulation reached low concentrations at the bottom layer of the unsaturated zone.

Surface complexation reactions in sandy porous media: Effects of incomplete mixing and mass-transfer limitations in flow-through systems

Although mixing and surface complexation reactions are key processes for solute transport in porous media, their coupling has not been extensively investigated. In this work, we study the impact of mass-transfer limitations on heterogeneous reactions taking place at the solid-solution interface of a natural sandy porous medium under advection-dominated flow-through conditions. A comprehensive set of 36 column experiments with different grain sizes (0.64, 1.3 and 2.3 mm), seepage velocities (1, 30 and 90 m/day), and hydrochemical conditions were performed. The injection of NaBr solutions of different concentrations (1-100 mM) led to the release of protons via deprotonation reactions of the quartz surface. pH and solute concentration breakthrough curves were measured at the outlet of the columns and the propagation of pH fronts in the column setups was tracked inside the porous medium with non-invasive optode sensors. The experimental results show that the deprotonation of the reactive surfaces, resulting from their interactions with the injected ionic species, strongly depends on the hydrodynamic conditions and differs among the tested porous media despite their apparent similar surface properties. Reactive transport modeling was used to quantitatively interpret the experimental results and to analyze the effects of mass-transfer limited physical processes on surface complexation reactions, propagation of pH fronts, transport of major ions and spatio-temporal evolution of surface composition. A dual domain mass transfer formulation (DDMT) combined with a surface complexation model (SCM) allowed capturing the effects of incomplete mixing on the surface reactions and to reproduce the experimental observations collected in the experiments with high flow velocities. The SCM was parameterized with a single set of surface complexation parameters, accounting for the similar surface properties of the porous media, and was capable of describing the surface complexation mechanisms and their impact on the hydrochemistry over the large range of tested ionic strengths.

Linking mixing and flow topology in porous media: An experimental proof

Transport processes in porous media are controlled by the characteristics of the flow field which are determined by the porous material properties and the boundary conditions of the system. This work provides experimental evidence of the relation between mixing and flow field topology in porous media at the continuum scale. The setup consists of a homogeneously packed quasi-two-dimensional flow-through chamber in which transient flow conditions, dynamically controlled by two external reservoirs, impact the transport of a dissolved tracer. The experiments were performed at two different flow velocities, corresponding to Péclet numbers of 191 and 565, respectively. The model-based interpretation of the experimental results shows that high values of the effective Okubo-Weiss parameter, driven by the changes of the boundary conditions, lead to high rates of increase of the Shannon entropy of the tracer distribution and, thus, to enhanced mixing. The comparison between a hydrodynamic dispersion model and an equivalent pore diffusion model demonstrates that despite the spatial and temporal variability in the hydrodynamic dispersion coefficients, the Shannon entropy remains almost unchanged because it is controlled by the Okubo-Weiss parameter. Overall, our work demonstrates that under highly transient boundary conditions, mixing dynamics in homogeneous porous media can also display complex patterns and is controlled by the flow topology.

Hyporheic transverse mixing zones and dispersivity: Laboratory and numerical experiments of hydraulic controls

Mixing of surface water and groundwater in shallow sediments is important to biogeochemical cycling and contaminant migration, and is often used to define the hyporheic zone. Yet knowledge of mixing processes in hyporheic zones is supported by surprisingly few rigorous lab or field observations, and differ from those in deeper groundwater by presence of enhanced head gradients, sediment heterogeneity, and temporal fluctuations. In a laboratory sediment (sand) tank we photographed a conservative dye to analyze transverse mixing zones between upwelling groundwater and bidirectional hyporheic exchange flows. We then conducted numerical modeling to investigate processes behind observed phenomena and estimate dispersivities. We found that transverse mixing zones were thin (i.e. mixing thickness measured in direction of steepest concentration gradient, δ, less than 5 cm), consistent with a small calibrated transverse dispersivity (~0.1 mm) and prior lab studies conducted at similar scales. In steady-state experiments and simulations, δ and estimated dispersion coefficients increased with the surface water head drop driving exchange flows. Given relatively constant deeper groundwater heads, increased Δh led to increased mixing zone length for both steady-state and transient conditions, indicating larger bedforms or weaker gaining conditions enhance subsurface mixing. However, Peclet number and flux-related dilution index simultaneously increased and decreased, respectively, indicating that enhancement of subsurface advection outpaced that of dispersion. In transient experiments and simulations, δ was greater than for steady-state, probably from temporary addition of longitudinal dispersion. During transient experiments, δ exhibited temporal noise, perhaps due to the mixing zone moving past varying patterns of sediment packing. Our results provide basic knowledge of mixing zone behavior in hyporheic zones with implications for hyporheic zone definitions, solute transport, mixing-dependent reaction, and water quality.

Magnitude of Diffusion- and Transverse Dispersion-Induced Isotope Fractionation of Organic Compounds in Aqueous Systems

Determining whether aqueous diffusion and dispersion lead to significant isotope fractionation is important for interpreting the isotope ratios of organic contaminants in groundwater. We performed diffusion experiments with modified Stokes diaphragm cells and transverse-dispersion experiments in quasi-two-dimensional flow-through sediment tank systems to explore isotope fractionation for benzene, toluene, ethylbenzene, 2,6-dichlorobenzamide, and metolachlor at natural isotopic abundance. We observed very small to negligible diffusion- and transverse-dispersion-induced isotope enrichment factors (ε < -0.4 ‰), with changes in carbon and nitrogen isotope values within ±0.5‰ and ±1‰, respectively. Isotope effects of diffusion did not show a clear correlation with isotopologue mass with calculated power-law exponents β close to zero (0.007 < β < 0.1). In comparison to ions, noble gases, and labeled compounds, three aspects stand out. (i) If a mass dependence is derived from collision theory, then isotopologue masses of polyatomic molecules would be affected by isotopes of multiple elements resulting in very small expected effects. (ii) However, collisions do not necessarily lead to translational movement but can excite molecular vibrations or rotations minimizing the mass dependence. (iii) Solute-solvent interactions like H-bonds can further minimize the effect of collisions. Modeling scenarios showed that an inadequate model choice, or erroneous choice of β, can greatly overestimate the isotope fractionation by diffusion and, consequently, transverse dispersion. In contrast, available data for chlorinated solvent and gasoline contaminants at natural isotopic abundance suggest that in field scenarios, a potential additional uncertainty from aqueous diffusion or dispersion would add to current instrumental uncertainties on carbon or nitrogen isotope values (±1‰) with an additional ±1‰ at most.

Impact of diffuse layer processes on contaminant forward and back diffusion in heterogeneous sandy-clayey domains

Low-permeability aquitards can significantly affect the transport, distribution, and persistence of contaminant plumes in subsurface systems. Although such low-permeability materials are often charged, the key role of charge-induced electrostatic processes during contaminant transport has not been extensively studied. This work presents a detailed investigation exploring the coupled effects of heterogeneous distribution of physical, chemical and electrostatic properties on reactive contaminant transport in field-scale groundwater systems including spatially distributed clay zones. We performed an extensive series of numerical experiments in three distinct heterogeneous sandy-clayey domains with different levels of complexity. The flow and reactive transport simulations were performed by explicitly resolving the complex velocity fields, the small-scale electrostatic processes, the compound-specific diffusive/dispersive fluxes and the chemical processes utilizing a multi-continua based reactive transport code (MMIT-Clay). In each particular domain, numerical experiments were performed focusing on both the forward and back diffusion through the sandy-clayey interfaces. The results illuminate the control of microscopic electrostatic mechanisms on macroscopic mass transfer. Coulombic interactions in the clay's diffuse layer can significantly accelerate or retard a particular species depending on its charge. Furthermore, the chemical heterogeneity plays a major role in mass storage and release during reactive transport. Neglecting such processes can lead to substantial over- or underestimation of the overall transport behavior, which underlines the need for integrated physical, chemical and electrostatic approaches to accurately describe mass transfer processes in systems including low-permeability inclusions.

Intramolecular non-covalent isotope effects at natural abundance associated with the migration of paracetamol in solid matrices during liquid chromatography

Position-specific isotope analysis by Nuclear Magnetic Resonance spectrometry was employed to study the ¹³C intramolecular isotopic fractionation associated with the migration of organic substrates through different stationary phases chromatography columns. Liquid chromatography is often used to isolate compounds prior to their isotope analysis and this purification step potentially alters the isotopic composition of target compounds introducing a bias in the later measured data. Moreover, results from liquid chromatography can yield the sorption parameters needed in reactive transport models that predict the transport and fate of organic contaminants to in the environment. The aim of this study was to use intramolecular isotope analysis to study both ¹³C and ¹⁵N isotope effects associated with the elution of paracetamol (acetaminophen) through different stationary phases and to compare them to effects observed previously for vanillin. Results showed very different intramolecular isotope fractionation profiles depending on the chemical structure of the stationary phase. The data also demonstrate that both the amplitude and the distribution of measured isotope effects depend on the nature of the non-covalent interactions involved in the migration process. Results provided by theoretical calculation performed during this study also confirmed the direct link between observed intramolecular isotope fractionation and the nature of involved intermolecular interactions. It is concluded that the nature of the stationary phase through which the substrate passes has a major impact on the intramolecular isotopic composition of organic compounds isolated by chromatography methods..

Plume deformation, mixing, and reaction kinetics: An analysis of interacting helical flows in three-dimensional porous media

Heterogeneity and macroscopic anisotropy of porous media play an important role for dilution and reaction enhancement of conservative and reactive plumes. In this study, we perform numerical simulations to investigate steady-state flow and transport in three-dimensional heterogeneous porous media. We consider two macroscopic anisotropic inclusions resulting in helical flows with twisting streamlines in a three-dimensional flow-through domain. The inclusions are obtained by alternating two layers of angled slices of coarse and fine porous media with different hydraulic conductivity. We investigate flow and transport scenarios considering different geometry and relative position of the two anisotropic inclusions yielding helical flow fields with different extent of interaction. We use metrics of stretching and folding to characterize the flow field and entropy-based metrics for the analysis of the conservative and reactive transport problems. The outcomes show that the two helices result in different patterns of twisting streamlines, which cause distinct deformation of the plumes. However, mixing and reaction enhancement could not be directly related to the extent of the flow field deformation: Configurations with strong deformation can result in only moderate mixing enhancement, whereas configurations with limited deformation of the flow field can lead to significant mixing of the solute plume. Finally, we explore the impact of different degradation rates on reactive transport and the role of reaction kinetics on the entropy balance for a reactant undergoing transport and mixing-controlled degradation in the twisting flow fields. The results show that strong mixing enhancement due to helical flow increases the importance of the reaction kinetics that becomes the rate-limiting process for solute reactive transport.

Assessment of 2,4-Dinitroanisole Transformation Using Compound-Specific Isotope Analysis after In Situ Chemical Reduction of Iron Oxides

Ferrous iron-bearing minerals are important reductants in the contaminated subsurface, but their availability for the reduction of anthropogenic pollutants is often limited by competition with other electron acceptors including microorganisms and poor accessibility to Fe(II) in complex hydrogeologic settings. The supply of external electron donors through in situ chemical reduction (ISCR) has been proposed as one remediation approach, but the quantification of pollutant transformation is complicated by the perturbations introduced to the subsurface by ISCR. Here, we evaluate the application of compound specific isotope analysis (CSIA) for monitoring the reduction of 2,4-dinitroanisole (DNAN), a component of insensitive munitions formulations, by mineral-bound Fe(II) generated through ISCR of subsurface material from two field sites. Electron balances from laboratory experiments in batch and column reactors showed that 3.6% to 11% of the total Fe in the sediments was available for the reduction of DNAN and its partially reduced intermediates after dithionite treatment. The extent of DNAN reduction was successfully quantified from its N isotope fractionation measured in the column effluent based on the derivation of a N isotope enrichment factor, ε_N, derived from a comprehensive series of isotope fractionation experiments with numerous Fe(II)-bearing minerals as well as dithionite-reduced subsurface materials. Our observations illustrate the utility of CSIA as a robust approach to evaluate the success of in situ remediation through abiotic contaminant reduction.

Effects of local transverse dispersion on macro-scale coefficients of oxygen-limited biodegradation in a stratified formation

The correct characterization of macro-scale contaminant transport and transformation rates is an important issue for modeling reactive transport in heterogeneous aquifers. While previous studies have investigated field-scale heterogeneity of transport and biochemical properties, the effects of local transverse dispersion on macro-scale transport and transformation rates have not been well understood. In this paper, the process of oxygen-limited biodegradation in a stratified aquifer is analysed by spectral perturbation approach, and longitudinal macrodispersivity, effective biodegradation rate, effective retardation factor and effective velocity are derived for the coupled transport equations of a system consisting of a contaminant and an oxidizing agent (oxygen). The effects of local transverse dispersion on these macro-scale coefficients are studied. It is shown that local transverse dispersion can smooth the heterogeneity in biodegradation and sorption processes and enlarge effective biodegradation rate and retardation factor. The local transverse dispersion can also limit the effects of heterogeneity in biodegradation process on longitudinal macrodispersivities and effective velocities for the contaminant and dissolved oxygen. But the effects of heterogeneity in sorption process on the contaminant macrodispersivity is likely to be magnified by local transverse dispersion.

Recent Advances in Experimental Studies of Steady-State Dilution and Reactive Mixing in Saturated Porous Media

Transverse dispersive mixing plays an important role in controlling natural attenuation of contaminant plumes and the performance of engineered remediation strategies. The extent of transverse mixing can be significantly affected by porous media heterogeneity and anisotropy. For instance, flow focusing in the high-permeability inclusions leads to an enhancement of dilution and reactive mixing in steady-state solute transport. Numerous modeling studies have been performed to understand the mechanism of conservative and reactive transport in homogeneous and complex heterogeneous porous media. However, experimental investigations are necessary to show an intuitive phenomenon and to validate the modeling results. This paper briefly reviews recent laboratory experimental studies on dilution and reactive mixing of steady-state transport in saturated homogeneous and heterogeneous porous media. In this context, setups and measuring techniques are described in pore-scale and Darcy-scale experiments. Parameters quantifying dilution and reactive mixing in the experiments are also introduced. Finally, we discuss the further experimental works necessary to deepen our understanding of dilution and reactive mixing in natural aquifers.

Sensing Coated Iron-Oxide Nanoparticles with Spectral Induced Polarization (SIP): Experiments in Natural Sand Packed Flow-Through Columns

The development of nanoparticle-based soil remediation techniques is hindered by the lack of accurate in situ nanoparticle (NP) monitoring and characterization methods. Spectral induced polarization (SIP), a noninvasive geophysical technique, offers a promising approach to detect and quantify NPs in porous media. However, its successful implementation as a monitoring tool requires an understanding of the polarization mechanisms, the governing NP-associated SIP responses and their dependence on the stabilizing coatings that are typically used for NPs deployed in environmental applications. Herein, we present SIP responses (0.1-10 000 Hz) measured during injection of a poloxamer-coated superparamagnetic iron-oxide nanoparticle (SPION) suspension in flow-through columns packed with natural sand from the Borden aquifer. An advective-dispersive transport model is fitted to outflow SPION concentration measurements to compute average concentrations over the SIP spatial response domain (within the columns). The average SPION concentrations are compared with the real and imaginary components of the complex conductivity. Excellent correspondence is found between the average SPION concentrations in the columns and the imaginary conductivity values, suggesting that NP-mediated polarization (that is, charge storage) increases proportionally with increasing SPION concentration. Our results support the possibility of SIP monitoring of spatial and temporal NP distributions, which can be immediately deployed in bench-scale studies with the prospect of future real-world field applications.

Simulation of Trinitrogen Migration and Transformation in the Unsaturated Zone at a Desert Contaminant Site (NW China) Using HYDRUS-2D

The protection of an unsaturated zone is essential for groundwater-quality security. Neglecting pollutant changes in the saturated zone can affect the accuracy of groundwater-quality assessments. Unlike water sampling, the nonreproducibility of soil sampling complicates the observation of contaminant changes at different times in the same location. The HYDRUS-2D model, coupled with the Richards equation and the convection-dispersion equation, was applied to simulate the migration and transformation of high ammonia concentrations in wastewater in an unsaturated zone. Long-term field observations were carried out for trinitrogen (NH4+, NO2−, and NO3−) from 2015 to 2018 at a wastewater discharge site located in a desert area in northwest China. Samples were collected twice a month. The model was calibrated and validated using statistics and observation data. Variations in trinitrogen concentrations were simulated using the model and fitted well with the measured values. Simulation results for trinitrogen migration and transformation demonstrated that there was no enrichment on the ground surface. Contaminants attenuated rapidly in the unsaturated zone after wastewater discharge stopped. NH4+ was oxidized to NO2− and NO3− under nitrification, except in the anoxic subclay lenses. Subclay lenses were not considered in previous research. These lenses had high enrichment with contaminants and prevented secondary nitrification, which might have led to extremely low NO3− concentrations. The removal rate of contaminants by the unsaturated zone in natural conditions is as high as 76%, and contaminants could be degraded to acceptable levels within 10 years (3650 days) without artificial interventions. This indicates that the unsaturated zone can delay migration and degrade contaminants, and should be taken into consideration in groundwater-quality assessments.

Benchmarking numerical codes for tracer transport with the aid of laboratory-scale experiments in 2D heterogeneous porous media

We present a combined experimental and numerical modeling study that addresses two principal questions: (i) is any particular Eulerian-based method used to solve the classical advection-dispersion equation (ADE) clearly superior (relative to the others), in terms of yielding solutions that reproduce BTCs of the kind that are typically sampled at the outlet of a laboratory cell? and (ii) in the presence of matches of comparable quality against such BTCs, do any of these methods render different (or similar) numerical BTCs at locations within the domain? To address these questions, we obtained measurements from carefully controlled laboratory experiments, and employ them as a reference against which numerical results are benchmarked and compared. The experiments measure solute transport breakthrough curves (BTCs) through a square domain containing various configurations of coarse, medium, and fine quartz sand. The approaches to solve the ADE involve Eulerian-Lagrangian and Eulerian (finite volume, finite elements, mixed and discontinuous finite elements) numerical methods. Model calibration is not examined; permeability and porosity of each sand were determined previously through separate, standard laboratory tests, while dispersivities are assigned values proportional to mean grain size. We find that the spatial discretization of the flow field is of critical importance, due to the non-uniformity of the domain. Although simulated BTCs at the system outlet are observed to be very similar for these various numerical methods, computed local (point-wise, inside the domain) BTCs can be very different. We find that none of the numerical methods is able to fully reproduce the measured BTCs. The impact of model parameter uncertainty on the calculated BTCs is characterized through a set of numerical Monte Carlo simulations; in cases where the impact is significant, assessment of simulation matches to the experimental data can be ambiguous.

Longitudinal dispersion coefficients for numerical modeling of groundwater solute transport in heterogeneous formations

Most recent research on hydrodynamic dispersion in porous media has focused on whole-domain dispersion while other research is largely on laboratory-scale dispersion. This work focuses on the contribution of a single block in a numerical model to dispersion. Variability of fluid velocity and concentration within a block is not resolved and the combined spreading effect is approximated using resolved quantities and macroscopic parameters. This applies whether the formation is modeled as homogeneous or discretized into homogeneous blocks but the emphasis here being on the latter. The process of dispersion is typically described through the Fickian model, i.e., the dispersive flux is proportional to the gradient of the resolved concentration, commonly with the Scheidegger parameterization, which is a particular way to compute the dispersion coefficients utilizing dispersivity coefficients. Although such parameterization is by far the most commonly used in solute transport applications, its validity has been questioned. Here, our goal is to investigate the effects of heterogeneity and mass transfer limitations on block-scale longitudinal dispersion and to evaluate under which conditions the Scheidegger parameterization is valid. We compute the relaxation time or memory of the system; changes in time with periods larger than the relaxation time are gradually leading to a condition of local equilibrium under which dispersion is Fickian. The method we use requires the solution of a steady-state advection-dispersion equation, and thus is computationally efficient, and applicable to any heterogeneous hydraulic conductivity K field without requiring statistical or structural assumptions. The method was validated by comparing with other approaches such as the moment analysis and the first order perturbation method. We investigate the impact of heterogeneity, both in degree and structure, on the longitudinal dispersion coefficient and then discuss the role of local dispersion and mass transfer limitations, i.e., the exchange of mass between the permeable matrix and the low permeability inclusions. We illustrate the physical meaning of the method and we show how the block longitudinal dispersivity approaches, under certain conditions, the Scheidegger limit at large Péclet numbers. Lastly, we discuss the potential and limitations of the method to accurately describe dispersion in solute transport applications in heterogeneous aquifers.

Rethinking aqueous phase diffusion related isotope fractionation: Contrasting theoretical effects with observations at the field scale

Aqueous phase diffusion-related isotope fractionation (DRIF) was investigated through modelling to determine under what subsurface conditions carbon isotope DRIF effects would be observable using typical sampling approaches. A dispersive enrichment factor was defined based on heavy and light isotopologue dispersion coefficients. For a given ratio of source concentration (C₀) to method detection limit (MDL), the maximum DRIF in a system increased linearly with transverse dispersive enrichment factor. Using this linear relationship, the critical dispersion enrichment factor for which DRIF would not be observable was quantified. Dispersive enrichment factors for various contaminants (benzene, toluene, chlorinated compounds) were estimated using field scale transverse dispersion coefficients upscaled from compound specific or non-compound specific local scale dispersivity. All predicted dispersive enrichment factors with non-compound specific dispersivity are smaller than critical values even for high C₀/MDL ratios (e.g. 25,000), indicating DRIF would generally not be observable in systems where soil dispersivity is non-compound specific. To date, this finding has not been clearly articulated in the DRIF literature. While the calculated dispersive enrichment factors for some compounds with compound specific transverse dispersivity exceeded the critical values at which DRIF might become significant, the zones in which DRIF could potentially be observable were limited to bands below 0.5m wide on lateral plume edges. In aquifer-aquitard systems, DRIF was theoretically detectable only in thin aquifers (e.g. 0.5m) bounded by thick (e.g. meters) aquitards. DRIF due to back diffusion from aquitards would not be observable regardless of aquifer thickness. Simulations addressing the mixing effect in wells demonstrated that DRIF effects would be difficult to identify in the field without a sampling strategy (including smaller than industry norm well screen lengths, and fine sampling scales) expressly targeted towards that goal. The results of this study help identify what the required characteristics of such a field strategy might be.

Impact of diffusion on transverse dispersion in two-dimensional ordered and random porous media

Solute dispersion in fluid flow results from the interaction between advection and diffusion. The relative contributions of these two mechanisms to mass transport are characterized by the reduced velocity ν, also referred to as the Péclet number. In the absence of diffusion (i.e., when the solute diffusion coefficient D_{m}=0 and ν→∞), divergence-free laminar flow of an incompressible fluid results in a zero-transverse dispersion coefficient (D_{T}=0), both in ordered and random two-dimensional porous media. We demonstrate by numerical simulations that a more realistic realization of the condition ν→∞ using D_{m}≠0 and letting the fluid flow velocity approach infinity leads to completely different results for ordered and random two-dimensional porous media. With increasing reduced velocity, D_{T} approaches an asymptotic value in ordered two-dimensional porous media but grows linearly in disordered (random) structures depending on the geometrical disorder of a structure: a higher degree of heterogeneity results in a stronger growth of D_{T} with ν. The obtained results reveal that disorder in the geometrical structure of a two-dimensional porous medium leads to a growth of D_{T} with ν even in a uniform pore-scale advection field; however, lateral diffusion is a prerequisite for this growth. By contrast, in ordered two-dimensional porous media the presence of lateral diffusion leads to a plateau for the transverse dispersion coefficient with increasing ν.

Experimental investigation of transverse mixing in porous media under helical flow conditions

Plume dilution and transverse mixing can be considerably enhanced by helical flow occurring in three-dimensional heterogeneous anisotropic porous media. In this study, we perform tracer experiments in a fully three-dimensional flow-through chamber to investigate the effects of helical flow on plume spiraling and deformation, as well as on its dilution. Porous media were packed in angled stripes of materials with different grain sizes to create blocks with macroscopically anisotropic hydraulic conductivity, which caused helical flows. Steady-state transport experiments were carried out by continuously injecting dye tracers at different inlet ports. High-resolution measurements of concentration and flow rates were performed at 49 outlet ports. These measurements allowed quantifying the spreading and dilution of the solute plumes at the outlet cross section. Direct evidence of plume spiraling and visual proof of helical flow was obtained by freezing and slicing the porous media at different cross sections and observing the dye-tracer distribution. We simulated flow and transport to interpret our experimental observations and investigate the effects of helical flow on mixing-controlled reactive transport. The simulation results were evaluated using metrics of reactive mixing such as the critical dilution index and the length of continuously injected steady-state plumes. The results show considerable reaction enhancement, quantified by the remarkable decrease of reactive plume lengths (up to four times) in helical flows compared to analogous scenarios in uniform flows.

Experimental Evidence of Helical Flow in Porous Media

Helical flow leads to deformation of solute plumes and enhances transverse mixing in porous media. We present experiments in which macroscopic helical flow is created by arranging different materials to obtain an anisotropic macroscopic permeability tensor with spatially variable orientation. The resulting helical flow entails twisting streamlines which cause a significant increase in lateral mass exchange and thus a large enhancement of plume dilution (up to 235%) compared to transport in homogenous media. The setup may be used to effectively mix solutes in parallel streams similarly to static mixers, but in porous media.

Impact of heterogeneity on oxygen transfer in a fluctuating capillary fringe

We performed quasi-two-dimensional flow through laboratory experiments to study the effect of a coarse-material inclusion, located in the proximity of the water table, on flow and oxygen transfer in the capillary fringe. The experiments investigate different phases of mass transfer from the unsaturated zone to anoxic groundwater under both steady-state and transient flow conditions, the latter obtained by fluctuating the water table. Monitoring of flow and transport in the different experimental phases was performed by visual inspection of the complex flow field using a dye tracer solution, measurement of oxygen profiles across the capillary fringe, and determination of oxygen fluxes in the effluent of the flow-through chamber. Our results show significant effects of the coarse-material inclusion on oxygen transfer during the different phases of the experiments. At steady state, the oxygen flux across the unsaturated/saturated interface was considerably enhanced due to flow focusing in the fully water-saturated coarse-material inclusion. During drainage, a zone of higher water saturation formed in the fine material overlying the coarse lens. The entrapped oxygen-rich aqueous phase contributed to the total amount of oxygen supplied to the system when the water table was raised back to its initial level. In case of imbibition, pronounced air entrapment occurred in the coarse lens, causing oxygen to partition between the aqueous and gaseous phases. The oxygen mass supplied to the anoxic groundwater following the imbibition event was found to be remarkably higher (approximately seven times) in the heterogeneous system compared with a similar experiment performed in a homogeneous porous medium.

Fringe-controlled biodegradation under dynamic conditions: quasi 2-D flow-through experiments and reactive-transport modeling

Biodegradation in contaminated aquifers has been shown to be most pronounced at the fringe of contaminant plumes, where mixing of contaminated water and ambient groundwater, containing dissolved electron acceptors, stimulates microbial activity. While physical mixing of contaminant and electron acceptor by transverse dispersion has been shown to be the major bottleneck for biodegradation in steady-state plumes, so far little is known on the effect of flow and transport dynamics (caused, e.g., by a seasonally fluctuating groundwater table) on biodegradation in these systems. Towards this end we performed experiments in quasi-two-dimensional flow-through microcosms on aerobic toluene degradation by Pseudomonas putida F1. Plume dynamics were simulated by vertical alteration of the toluene plume position and experimental results were analyzed by reactive-transport modeling. We found that, even after disappearance of the toluene plume for two weeks, the majority of microorganisms stayed attached to the sediment and regained their full biodegradation potential within two days after reappearance of the toluene plume. Our results underline that besides microbial growth, also maintenance and dormancy are important processes that affect biodegradation performance under transient environmental conditions and therefore deserve increased consideration in future reactive-transport modeling.

Evaluation of dispersivity coefficients by means of a laboratory image analysis

This paper describes the application of an innovative procedure that allows the estimation of longitudinal and transverse dispersivities in an experimental plume devised in a laboratory sandbox. The phenomenon of transport in porous media is studied using sodium fluorescein as tracer. The fluorescent excitation was achieved by using blue light and the concentration data were obtained through the processing of side wall images collected with a high resolution color digital camera. After a calibration process, the relationship between the luminosity of the emitted fluorescence and the fluorescein concentration was determined at each point of the sandbox. The relationships were used to describe the evolution of the transport process quantitatively throughout the entire domain. Some check tests were performed in order to verify the reliability of the experimental device. Numerical flow and transport models of the sandbox were developed and calibrated comparing computed and observed flow rates and breakthrough curves. The estimation of the dispersivity coefficients was carried out by analyzing the concentration field deduced from the images collected during the experiments; the dispersivity coefficients were evaluated in the domain zones where the tracer affected the porous medium under the hypothesis that the transport phenomenon is described by advection-dispersion equation (ADE) and by computing the differential components of the concentration by means of a numerical leap-frog scheme. The values determined agree with the ones referred in literature for similar media and with the coefficients obtained by calibrating the numerical model. Very interesting considerations have been made from the analysis of the performance of the methodology at different locations in the flow domain and phases of the plume evolution.

Experimental investigation of compound-specific dilution of solute plumes in saturated porous media: 2-D vs. 3-D flow-through systems

Dilution of solute plumes in groundwater strongly depends on transverse mixing. Thus, the correct parameterization of transverse dispersion is of critical importance for the quantitative description of solute transport. In this study we perform flow-through laboratory experiments to investigate the influence of transport dimensionality on transverse mixing. We present a high-resolution experimental setup to study solute dilution and transverse dispersion in three-dimensional porous media. We conduct multi-tracer experiments in the new 3-D setup and compare the results with the outcomes of analogous tracer experiments performed in a quasi 2-D system. We work under steady-state flow and transport conditions and consider a range of velocities relevant for groundwater flow (0.5-8 m/day). Transverse dispersion coefficients are determined from high-resolution concentration profiles at the outlet of the flow-through chambers (7×7 ports in the 3-D setup and 7 ports in the quasi 2-D system), considering conservative tracers with significantly different aqueous diffusion coefficients, namely fluorescein and dissolved oxygen. To quantify dilution in the 2-D and 3-D systems, we experimentally determine the flux-related dilution index using the flow rates and the concentrations measured at the inlet and outlet ports, and we propose semi-analytical expressions to predict its evolution with travel distance in uniform groundwater flow. The experimental results in the quasi 2-D and 3-D flow-through systems are consistent and show a compound-specific behavior of the transverse dispersion coefficient and its non-linear dependence on the seepage velocity in both setups. The degree of dilution and the compound-specific effects of transverse dispersion are considerably more pronounced in 3-D than in quasi 2-D transport systems.

Reaction front formation in contaminant plumes

The formation of successive fronts in contaminated groundwater plumes by subsoil bacterial action is a commonly accepted feature of their propagation, but it is not obviously clear from a mathematical standpoint quite how such fronts are formed or propagate. In this paper we show that these can be explained by combining classical reaction-diffusion theory involving just two reactants (oxidant and reductant), and a secondary reaction in which a reactant on one side of such a front is (re-)formed on the other side of the front via diffusion of its product across the front. We give approximate asymptotic solutions for the reactant profiles, and the propagation rate of the front.

Experimental sensitivity analysis of oxygen transfer in the capillary fringe

Oxygen transfer in the capillary fringe (CF) is of primary importance for a wide variety of biogeochemical processes occurring in shallow groundwater systems. In case of a fluctuating groundwater table two distinct mechanisms of oxygen transfer within the capillary zone can be identified: vertical predominantly diffusive mass flux of oxygen, and mass transfer between entrapped gas and groundwater. In this study, we perform a systematic experimental sensitivity analysis in order to assess the influence of different parameters on oxygen transfer from entrapped air within the CF to underlying anoxic groundwater. We carry out quasi two-dimensional flow-through experiments focusing on the transient phase following imbibition to investigate the influence of the horizontal flow velocity, the average grain diameter of the porous medium, as well as the magnitude and the speed of the water table rise. We present a numerical flow and transport model that quantitatively represents the main mechanisms governing oxygen transfer. Assuming local equilibrium between the aqueous and the gaseous phase, the partitioning process from entrapped air can be satisfactorily simulated. The different experiments are monitored by measuring vertical oxygen concentration profiles at high spatial resolution with a noninvasive optode technique as well as by determining oxygen fluxes at the outlet of the flow-through chamber. The results show that all parameters investigated have a significant effect and determine different amounts of oxygen transferred to the oxygen-depleted groundwater. Particularly relevant are the magnitude of the water table rise and the grain size of the porous medium.

On the importance of diffusion and compound-specific mixing for groundwater transport: an investigation from pore to field scale

Mixing processes significantly affect and limit contaminant transport and transformation rates in the subsurface. The correct quantification of mixing in groundwater systems must account for diffusion, local-scale dispersion and the flow variability in heterogeneous flow fields (e.g., flow-focusing in high-conductivity and de-focusing in low-conductivity zones). Recent results of multitracer laboratory experiments revealed the significant effect of compound-specific diffusive properties on the physical displacement of dissolved species across a representative range of groundwater flow velocities. The goal of this study is to investigate the role of diffusion and compound-specific mixing for solute transport across a range of scales including: (i) pore-scale (~10⁻² m), (ii) laboratory bench-scale (~10⁰ m) and (iii) field-scale (~10² m). We investigate both conservative and mixing-controlled reactive transport using pore-scale modeling, flow-through laboratory experiments and simulations, and field-scale numerical modeling of complex heterogeneous hydraulic conductivity fields with statistical properties similar to the ones reported for the extensively investigated Borden aquifer (Ontario, Canada) and Columbus aquifer (Mississippi, USA, also known as MADE site). We consider different steady-state and transient transport scenarios. For the conservative cases we use as a metric of mixing the exponential of the Shannon entropy to quantify solute dilution either in a given volume (dilution index) or in a given solute flux (flux-related dilution index). The decrease in the mass and the mass-flux of the contaminant plumes is evaluated to quantify reactive mixing. The results show that diffusive processes, occurring at the small-scale of a pore channel, strongly affect conservative and reactive solute transport at larger macroscopic scales. The outcomes of our study illustrate the need to consider and properly account for compound-specific diffusion and mixing limitations in order to accurately describe and predict conservative and reactive transport in porous media.

Estimating pathway-specific contributions to biodegradation in aquifers based on dual isotope analysis: theoretical analysis and reactive transport simulations

At field sites with varying redox conditions, different redox-specific microbial degradation pathways contribute to total contaminant degradation. The identification of pathway-specific contributions to total contaminant removal is of high practical relevance, yet difficult to achieve with current methods. Current stable-isotope-fractionation-based techniques focus on the identification of dominant biodegradation pathways under constant environmental conditions. We present an approach based on dual stable isotope data to estimate the individual contributions of two redox-specific pathways. We apply this approach to carbon and hydrogen isotope data obtained from reactive transport simulations of an organic contaminant plume in a two-dimensional aquifer cross section to test the applicability of the method. To take aspects typically encountered at field sites into account, additional simulations addressed the effects of transverse mixing, diffusion-induced stable-isotope fractionation, heterogeneities in the flow field, and mixing in sampling wells on isotope-based estimates for aerobic and anaerobic pathway contributions to total contaminant biodegradation. Results confirm the general applicability of the presented estimation method which is most accurate along the plume core and less accurate towards the fringe where flow paths receive contaminant mass and associated isotope signatures from the core by transverse dispersion. The presented method complements the stable-isotope-fractionation-based analysis toolbox. At field sites with varying redox conditions, it provides a means to identify the relative importance of individual, redox-specific degradation pathways.

Model complexity needed for quantitative analysis of high resolution isotope and concentration data from a toluene-pulse experiment

Separating microbial- and physical-induced effects on the isotope signals of contaminants has been identified as a challenge in interpreting compound-specific isotope data. In contrast to simple analytical tools, such as the Rayleigh equation, reactive-transport models can account for complex interactions of different fractionating processes. The question arises how complex such models must be to reproduce the data while the model parameters remain identifiable. In this study, we reanalyze the high-resolution data set of toluene concentration and toluene-specific δ(13)C from the toluene-pulse experiment performed by Qiu et al. (this issue). We apply five reactive-transport models, differing in their degree of complexity. We uniquely quantify degradation and sorption properties of the system for each model, estimate the contributions of biodegradation-induced, sorption-induced, and transverse-dispersion-induced isotope fractionation to the overall isotope signal, and investigate the error introduced in the interpretation of the data when individual processes are neglected. Our results show that highly resolved data of both concentration and isotope ratios are needed for unique process identification facilitating reliable model calibration. Combined analysis of these highly resolved data demands reactive transport models accounting for nonlinear degradation kinetics and isotope fractionation by both reactive and physical processes such as sorption and transverse dispersion.

The behavior of effective rate constants for bimolecular reactions in an asymptotic transport regime

Previous research has shown that rate constants measured in batch tests (κ) may over-predict the amount of product formation when used in continuum models, and that these rate constants are often much greater than effective ones (κ(eff)) determined from upscaling studies. However, there is evidence that mixing is more important than the rate constants when using upscaled models. We use a numerical two-dimensional pore-scale porous medium with an approach similar to an experimental column test, and focus on the scenario of the displacement and mixing of two solutions with irreversible bimolecular reactions. Break-through curves of multiple cross-sectional averaged concentrations are analyzed for conservative and reactive transport, as well as the segregation of reactant species along the cross-sections. We compute effective parameters for the continuum scale in order to better understand the impact of using intrinsic rate constants in upscaled models. For a range of Damköhler numbers (Da), we compute effective reaction rate parameters and a reaction effectiveness factor; the latter is described by an empirical formula that depends on the Damköhler number and captures the upscaled system behavior. Our pore-scale results also confirm the segregation concept advanced by Kapoor et al. (1997). We find that for Da>1, κ(eff)<<κ, and yet the relative difference in total mass transformation between the pore-scale simulation and what is predicted by the upscaled continuum model using κ is about 10%. The explanation for this paradox is that the early transition of the regime from rate-limited to mixing-limited results in a model that is relatively insensitive to the rate constant because mixing controls the availability of reactants. Thus, the reaction-rate parameter used in the model has limited influence on the rate of product computed.

Multispecies hydrodynamic dispersion under high concentration gradients

Recent research suggests that when high concentration gradients (HCG) are present, resulting sharp density differences can cause the dispersive flux relationship to deviate from its classical Fickian form. This paper presents stable, upward, miscible displacement experiments conducted in two different types of porous media for a wide range of concentration differences between resident and displacing fluids. The considered groundwater velocities ranged from advection-dominated transport to velocities where the contribution of molecular diffusion is important, with the corresponding Peclet numbers ranging from 0.2 to 320. In addition to single component displacing fluids, mixtures consisting of multiple solutes were considered. The results of this study provide further evidence that classical Fick's law over-estimates the dispersion coefficient under HCG conditions. The decrease in the apparent dispersion coefficient is shown to be a nonlinear function of both concentration difference and groundwater velocity. This observation is attributed to gravitational effects at the sub-continuum scale which are not directly accounted for in classical variable density advection/dispersion models. Mixture experiments showed that the dispersive behaviors of individual components in a groundwater contaminant mixture are coupled.

Numerical simulation of isotope fractionation in steady-state bioreactive transport controlled by transverse mixing

Compound-specific stable isotope analysis (CSIA) has increasingly been used as a tool to assess intrinsic biodegradation at contaminated field sites. Typically, the Rayleigh equation is used to estimate the extent of biodegradation from measured isotope ratios of the contaminant. However, if the rate-limiting step in overall degradation is not the microbial reaction itself, the Rayleigh equation may no more be applicable. In this study we simulate biodegradation of continuously emitted petroleum hydrocarbons in groundwater systems. These contaminants are effectively degraded at the plume fringe where transverse dispersion makes them mix with dissolved electron acceptors present in the ambient groundwater. We simulate reactive transport to study the coupled effects of transverse mixing and biodegradation on the spatial patterns of carbon isotope signatures and their interpretation based on depth-integrated sampling which represents the most common setup in the assessment of contaminated sites. We present scenarios mimicking a hydraulically uniform laboratory experiment and a field-scale application considering heterogeneous conductivity fields. We compare cases in which isotopologue-specific transverse dispersion is considered to those where this additional fractionation process is neglected. We show that these effects cause significant shifts in the isotopic signals and may lead to overestimation of biodegradation. Moreover, our results provide evidence that the rate-limiting effect of transverse mixing on the overall degradation spatially varies along the length of a steady-state contaminant plume. The control of biodegradation by transverse dispersion and the fractionating effect of dispersion modulate the fractionation caused by the microbial reaction alone. As a consequence, significantly nonlinear isotopic patterns are observed in a Rayleigh plot. Simulations in heterogeneous flow domains show that these effects persist at larger field scales and are sensitive to the degree of mixing enhancement, determined by the heterogeneity of the hydraulic conductivity fields, and to the groundwater flow velocity.

Stochastic evaluation of mixing-controlled steady-state plume lengths in two-dimensional heterogeneous domains

We study plumes originating from continuous sources that require a dissolved reaction partner for their degradation. The length of such plumes is typically controlled by transverse mixing. While analytical expressions have been derived for homogeneous flow fields, incomplete characterization of the hydraulic conductivity field causes uncertainty in predicting plume lengths in heterogeneous domains. In this context, we analyze the effects of three sources of uncertainty: (i) The uncertainty of the effective mixing rate along the plume fringes due to spatially varying flow focusing, (ii) the uncertainty of the volumetric discharge through (and thus total mass flux leaving) the source area, and (iii) different parameterizations of the Darcy-scale transverse dispersion coefficient. The first two are directly related to heterogeneity of hydraulic conductivity. In this paper, we derive semi-analytical expressions for the probability distribution of plume lengths at different levels of complexity. The results are compared to numerical Monte Carlo simulations. Uncertainties in mixing and in the source strength result in a statistical distribution of possible plume lengths. For unconditional random hydraulic conductivity fields, plume lengths may vary by more than one order of magnitude even for moderate degrees of heterogeneity. Our results show that the uncertainty of volumetric flux through the source is the most relevant contribution to the variance of the plume length. The choice of different parameterizations for the local dispersion coefficient leads to differences in the mean estimated plume length.

Acid/base front propagation in saturated porous media: 2D laboratory experiments and modeling

We perform laboratory scale reactive transport experiments involving acid-basic reactions between nitric acid and sodium hydroxide. A two-dimensional experimental setup is designed to provide continuous on-line measurements of physico-chemical parameters such as pH, redox potential (Eh) and electrical conductivity (EC) inside the system under saturated flow through conditions. The electrodes provide reliable values of pH and EC, while sharp fronts associated with redox potential dynamics could not be captured. Care should be taken to properly incorporate within a numerical model the mixing processes occurring inside the electrodes. The available observations are modeled through a numerical code based on the advection-dispersion equation. In this framework, EC is considered as a variable behaving as a conservative tracer and pH and Eh require solving the advection dispersion equation only once. The agreement between the computed and measured pH and EC is good even without recurring to parameters calibration on the basis of the experiments. Our findings suggest that the classical advection-dispersion equation can be used to interpret these kinds of experiments if mixing inside the electrodes is adequately considered.

Transverse hydrodynamic dispersion effects on isotope signals in groundwater chlorinated solvents' plumes

The effects of transverse hydrodynamic dispersion on altering transformation-induced compound-specific isotope analysis (CSIA) signals within groundwater pollution plumes have been assessed with reactive transport modeling accommodating diffusion-induced isotope fractionation (DIF) and implementing different parameterizations of local transverse dispersion. The model reproduced previously published field data showing a negative carbon isotope pattern (-2 ‰) at the fringes of a nondegrading PCE plume. We extended the study to reactive transport scenarios considering vinyl chloride as a model compound and assessing, through a detailed sensitivity analysis, the coupled effects of transverse hydrodynamic dispersion (with and without DIF) and aerobic fringe degradation on the evolution of carbon and chloride isotope ratios. Transformation-induced positive isotope signals were increasingly attenuated with distance from the source and higher degradation rate. The effect of DIF on the overall isotope signal attenuation was greatest near the source and for low values of groundwater flow velocity, transverse dispersion coefficient, molecular weight, rate constant, and isotope fractionation factor, α, of the degradation reaction. Models disregarding DIF underestimate the actual α. The approximately twice larger DIF effect for chlorine than for carbon together with the low α for oxidation resulted in strong chlorine CSIA depletions for VC at the plume fringe.

Evaluation of transverse dispersion effects in tank experiments by numerical modeling: parameter estimation, sensitivity analysis and revision of experimental design

Transverse dispersion represents an important mixing process for transport of contaminants in groundwater and constitutes an essential prerequisite for geochemical and biodegradation reactions. Within this context, this work describes the detailed numerical simulation of highly controlled laboratory experiments using uranine, bromide and oxygen depleted water as conservative tracers for the quantification of transverse mixing in porous media. Synthetic numerical experiments reproducing an existing laboratory experimental set-up of quasi two-dimensional flow through tank were performed to assess the applicability of an analytical solution of the 2D advection-dispersion equation for the estimation of transverse dispersivity as fitting parameter. The fitted dispersivities were compared to the "true" values introduced in the numerical simulations and the associated error could be precisely estimated. A sensitivity analysis was performed on the experimental set-up in order to evaluate the sensitivities of the measurements taken at the tank experiment on the individual hydraulic and transport parameters. From the results, an improved experimental set-up as well as a numerical evaluation procedure could be developed, which allow for a precise and reliable determination of dispersivities. The improved tank set-up was used for new laboratory experiments, performed at advective velocities of 4.9 m d(-1) and 10.5 m d(-1). Numerical evaluation of these experiments yielded a unique and reliable parameter set, which closely fits the measured tracer concentration data. For the porous medium with a grain size of 0.25-0.30 mm, the fitted longitudinal and transverse dispersivities were 3.49×10(-4) m and 1.48×10(-5) m, respectively. The procedures developed in this paper for the synthetic and rigorous design and evaluation of the experiments can be generalized and transferred to comparable applications.

A high-resolution non-invasive approach to quantify oxygen transport across the capillary fringe and within the underlying groundwater

Oxygen transport across the capillary fringe is relevant for many biogeochemical processes. We present a non-invasive technique, based on optode technology, to measure high-resolution concentration profiles of oxygen across the unsaturated/saturated interface. By conducting a series of quasi two-dimensional flow-through laboratory experiments, we show that vertical hydrodynamic dispersion in the water-saturated part of the capillary fringe is the process limiting the mass transfer of oxygen. A number of experimental conditions were tested in order to investigate the influence of grain size and horizontal flow velocity on transverse vertical dispersion in the capillary fringe. In the same setup, analogous experiments were simultaneously carried out in the fully water-saturated zone, therefore allowing a direct comparison with oxygen transfer across the capillary fringe. The outcomes of the experiments under various conditions show that oxygen transport in the two zones of interest (i.e., the unsaturated/saturated interface and the saturated zone) is characterized by very similar transverse dispersion coefficients. An influence of the capillary fringe morphology on oxygen transport has not been observed. These results may be explained by the narrow grain size distribution used in the experiments, leading to a steep decline in water saturation at the unsaturated/saturated interface and to the absence of trapped gas in this transition zone. We also modeled flow (applying the van Genuchten and the Brooks-Corey relationships) and two-dimensional transport across the capillary fringe, obtaining simulated profiles of equivalent aqueous oxygen concentration that were in good agreement with the observations.

Isotopic fractionation by transverse dispersion: flow-through microcosms and reactive transport modeling study

Flow-through experiments were carried out to investigate the role of transverse dispersion on the isotopic behavior of an organic compound during conservative and bioreactive transport in a homogeneous porous medium. Ethylbenzene was selected as model contaminant and a mixture of labeled (perdeuterated) and light isotopologues was continuously injected in a quasi two-dimensional flow-through system. We observed a significant fractionation of ethylbenzene isotopologues during conservative transport at steady state. This effect was particularly pronounced at the plume fringe and contrasted with the common assumption that physical processes only provide a negligible contribution to isotope fractionation. Under the experimental steady state conditions, transverse hydrodynamic dispersion was the only process that could have caused the observed fractionation. Therefore, the measured isotope ratios at the outlet ports were interpreted with different parameterizations of the transverse dispersion coefficient. A nonlinear compound-specific parameterization showed the best agreement with the experimental data. Successively, bioreactive experiments were performed in two subsequent stages: a first oxic phase, involving a single strain of ethylbenzene degraders and a second phase with aerobic and anaerobic (i.e., ethylbenzene oxidation coupled to nitrate reduction) degradation. Significant fractionation through biodegradation occurred exclusively due to the metabolic activity of the anaerobic degraders. We performed analytical and numerical reactive transport simulations of the different experimental phases which confirmed that both the effects of physical processes (diffusion and dispersion) and microbially mediated reactions have to be considered to match the observed isotopic fractionation behavior.