Numerical simulation of a natural gradient tracer experiment for the natural attenuation study: flow and physical transport

Relationship between dual-domain parameters and practical characterization data

Dual-domain solute transport models produce significantly improved agreement to observations compared to single-domain (advection-dispersion) models when used in an a posteriori data fitting mode. However, the use of dual-domain models in a general predictive manner has been a difficult and persistent challenge, particularly at field-scale where characterization of permeability and flow is inherently limited. Numerical experiments were conducted in this study to better understand how single-rate mass transfer parameters vary with aquifer attributes and contaminant exposure. High-resolution reference simulations considered 30 different scenarios involving variations in permeability distribution, flow field, mass transfer timescale, and contaminant exposure time. Optimal dual-domain transport parameters were empirically determined by matching to breakthrough curves from the high-resolution simulations. Numerical results show that mobile porosity increases with lower permeability contrast/variance, smaller spatial correlation length, lower connectivity of high-permeability zones, and flow transverse to strata. A nonzero non-participating porosity improves empirical fitting, and becomes larger for flow aligned with strata, smaller diffusion coefficient, and larger spatial correlation length. The non-dimensional mass transfer coefficient or Damkohler number tends to be close to 1.0 and decrease with contaminant exposure time, in agreement with prior studies. The best empirical fit is generally achieved with a combination of macrodispersion and first-order mass transfer. Quantitative prediction of ensemble-average dual-domain parameters as a function of measurable aquifer attributes proved only marginally successful.

Lessons learned from 25 years of research at the MADE site

Field studies at well-instrumented research sites have provided extensive data sets and important insights essential for development and testing of transport theories and mathematical models. This paper provides an overview of over 25 years of research and lessons learned at one of such field research sites on the Columbus Air Force Base in Mississippi, commonly known as the Macrodispersion Experiment (MADE) site. Since the mid-1980s, field data from the MADE site have been used extensively by researchers around the world to explore complex contaminant transport phenomena in highly heterogeneous porous media. Results from field investigations and modeling analyses suggested that connected networks of small-scale preferential flow paths and relative flow barriers exert dominant control on solute transport processes. The classical advection-dispersion model was shown to inadequately represent plume-scale transport, while the dual-domain mass transfer model was found to reproduce the primary observed plume characteristics. The MADE site has served as a valuable natural observatory for contaminant transport studies where new observations have led to better understanding and improved models have sprung out analysis of new data.

Investigation of small-scale preferential flow with a forced-gradient tracer test

A new tracer experiment (referred to as MADE-5) was conducted at the well-known Macrodispersion Experiment (MADE) site to investigate the influence of small-scale mass-transfer and dispersion processes on well-to-well transport. The test was performed under dipole forced-gradient flow conditions and concentrations were monitored in an extraction well and in two multilevel sampler (MLS) wells located at 6, 1.5, and 3.75 m from the source, respectively. The shape of the breakthrough curve (BTC) measured at the extraction well is strongly asymmetric showing a rapidly arriving peak and an extensive late-time tail. The BTCs measured at seven different depths in the two MLSs are radically different from one another in terms of shape, arrival times, and magnitude of the concentration peaks. All of these characteristics indicate the presence of a complex network of preferential flow pathways controlling solute transport at the test site. Field-experimental data were also used to evaluate two transport models: a stochastic advection-dispersion model (ADM) based on conditional multivariate Gaussian realizations of the hydraulic conductivity field and a dual-domain single-rate (DDSR) mass-transfer model based on a deterministic reconstruction of the aquifer heterogeneity. Unlike the stochastic ADM realizations, the DDSR accurately predicted the magnitude of the concentration peak and its arrival time (within a 1.5% error). For the multilevel BTCs between the injection and extraction wells, neither model reproduced the observed values, indicating that a high-resolution characterization of the aquifer heterogeneity at the subdecimeter scale would be needed to fully capture 3D transport details.

Highlighting the implications of selenium dispersion from disposal of Kahota Industrial Triangle area, Islamabad, Pakistan using three-dimension solute transport model

This study highlights the implications of selenium (Se) dispersion in groundwater flow regimes of Kahota Industrial Triangle area located adjacent to the Soan River, Islamabad. Initially, a regional groundwater 3-D flow model has been developed, calibrated to the known observed heads of 24 water wells, verified, and confirmed that convergence has actually arrived to satisfy the steady state condition. Later, the transient simulation was carried out adding in the known recharge, storage factor, porosity, and observed drawdown matched with the simulated drawdown that appears to fall in close agreement with a difference of 0.25 m. As such the steady state groundwater model has facilitated to understand the mechanism of groundwater flow regimes in reference to the implications of selenium dispersion from disposal of Kahota Industrial Triangle area. Thirty-five water samples were collected mainly from the industrial water wells for the evaluation of heavy metals. Selenium being the major contributor of pollution has been short listed to monitor its dispersion using a solute transport model modular three-dimensional transport model (MT3D). Chemical parameters related to selenium characteristics including horizontal and vertical transverse dispersivity/longitudinal dispersivity, effective molecular diffusion coefficient and bulk density of the porous medium of aquifers have been used in MT3D contaminant transport model. MT3D is run for 30 years in steady state condition. As usual first run did not produce the exact field conditions. Therefore, the contaminant transport model is calibrated against the 32 values of observed selenium concentrations in boreholes by minor adjustments in the chemical parameter values. The final calibration has been achieved with residual value of 3.88x10(-5) Kg/m3. Seven hypothetical observation wells are used to monitor the selenium concentrations over a long-term period of time.

Multi tracer test for the implementation of enhanced in-situ bioremediation at a BTEX-contaminated megasite

At the Centre for Environmental Research Leipzig-Halle (UFZ) research site in Zeitz, Germany, benzene contaminates the lower of two aquifers with concentrations of up to 20 mg/l. Since the benzene plume has a minimum length of approximately 1 km, enhanced natural attenuation measures are being considered as a remediation strategy. This study describes the performance and evaluation of a multi-species reactive tracer test using the tracers fluorescein and bromide as conservative tracers and toluene as reactive tracer. Sampling was performed over a period of six months using a detailed network of multilevel sampling wells. Toluene was only slightly retarded in comparison to bromide, whereas fluorescein was retarded considerably stronger. Therefore, it was not possible to use fluorescein as an in situ tracer for the determination of groundwater velocities. The ionic nature of fluorescein is assumed to be the major reason for its retardation. The results show that the infiltration conditions were suitable to produce a wide spreading of the tracer front along the full thickness of the aquifer. Thus, a large aquifer volume can be treated in future enhanced bioremediation measures. The total quantity of infiltrated toluene (24 l) was degraded under sulfate-reducing conditions over a flow path of 50 m. Benzylsuccinate was identified as a metabolite of toluene degradation under sulfate-reducing conditions at this site. The modelling results show that toluene degradation was described more accurately using Monod kinetics than first-order kinetics. Since toluene was only slightly retarded in comparison to bromide, sorption and desorption processes were considered to be negligible.

Geophysical constraints on contaminant transport modeling in a heterogeneous fluvial aquifer

Approximately 3000 measurements of hydraulic conductivity in over 50 flowmeter boreholes were available at the Macro-Dispersion Experiment (MADE) site in Columbus, Mississippi, USA to quantify the heterogeneity in hydraulic conductivity at the site scale. This high-density measurement approach is perhaps infeasible for time and expense in typical groundwater remediation sites. A natural-gradient tracer experiment from the MADE site was simulated by a groundwater flow and solute transport model incorporating direct-current (DC) resistivity data collected over the observed plume location. Hydraulic conductivity from one borehole collected during the original site characterization was used to calibrate the electrical resistivity data to hydraulic conductivity using a previously derived log-log relationship. Application of this relationship, using site-specific empirical constants determined from the data, transforms the 3D electrical resistivity data into a 3D description of hydraulic conductivity that can be used in groundwater models. The validity of this approach was tested by using the geophysically derived hydraulic conductivity representation in numerical simulations of the natural-gradient tracer experiment. The agreement between the simulated and observed tracer plumes was quantified to gauge the effectiveness of geophysically derived and flowmeter based representations of the hydraulic conductivity field. This study demonstrates that a highly heterogeneous aquifer can be modeled with minimal hydrological data supplemented with geophysical data at least as well as previous models of the site using purely hydrologic data.

Delineating alluvial aquifer heterogeneity using resistivity and GPR data

Conceptual geological models based on geophysical data can elucidate aquifer architecture and heterogeneity at meter and smaller scales, which can lead to better predictions of preferential flow pathways. The macrodispersion experiment (MADE) site, with >2000 measurements of hydraulic conductivity obtained and three tracer tests conducted, serves as an ideal natural laboratory for examining relationships between subsurface flow characteristics and geophysical attributes in fluvial aquifers. The spatial variation of hydraulic conductivity measurements indicates a large degree of site heterogeneity. To evaluate the usefulness of geophysical methods for better delineating fluvial aquifer heterogeneities and distribution of preferential flow paths, a surface grid of two-dimensional ground penetrating radar (GPR) and direct current (DC) resistivity data were collected. A geological model was developed from these data that delineate four stratigraphic units with distinct electrical and radar properties including (from top to bottom) (1) a meandering fluvial system (MFS); (2) a braided fluvial system (BFS); (3) fine-grained sands; and (4) a clay-rich interval. A paleochannel, inferred by other authors to affect flow, was mapped in the MFS with both DC resistivity and GPR data. The channel is 2 to 4 m deep and, based on resistivity values, is predominantly filled with clay and silt. Comparing previously collected hydraulic conductivity measurements and tracer-plume migration patterns to the geological model indicates that flow primarily occurs in the BFS and that the channel mapped in the MFS has no influence on plume migration patterns.

Effects of rate-limited mass transfer on water sampling with partially penetrating wells

Nonequilibrium concentration type curves are numerically developed and sensitivity analyses are performed to examine the relationships between effluent concentrations in partially penetrating monitoring/extraction wells, the vertical plume shape, and the mass transfer characteristics of the aquifer. The governing two-dimensional, axisymmetric nonequilibrium solute transport equation is solved in three stages using an operator-splitting approach. In the first two stages, the advection and dispersion terms are solved with the Eulerian-Lagrangian method, based on the backward method of characteristics for advection and the standard implicit Galerkin finite element method for dispersion. In the third step, the first-order, immobile-mobile domain mass transfer term is computed analytically for both two-site and lognormally distributed, multirate models. Effluent concentration variations with time and contour plots of the pore water concentration distribution in the aquifer are compared for a wide range of field- and laboratory-measured mass transfer rates, various plume shapes, and relevant physical/chemical parameter values, including pumping rate, vertical anisotropy ratio, retardation factor, and porosity. The simulation results show that rate-limited mass transfer can have a significant impact on sample and aquifer pore water concentrations during three-dimensional transport to a partially penetrating well. An alternative dimensionless form of the nonequilibrium solute transport equation is derived to illustrate the key parameter groupings that quantify rate-limited sorption effects and show the relative importance of individual parameters. A hypothetical field application example demonstrates the fitting of dimensional type curves to discrete-interval sampling data in order to evaluate the mass transfer characteristics of an aquifer and shows how type curve superposition can be used to model complex plume shapes.

Analysis of solute transport in flow fields influenced by preferential flowpaths at the decimeter scale

Several recent studies at the Macrodispersion Experiment (MADE) site in Columbus, Mississippi, have indicated that the relative preferential flowpaths and flow barriers resulting from decimeter-scale aquifer heterogeneities appear to have a dominant effect on plume-scale solute transport. Numerical experiments are thus conducted in this study to explore the key characteristics of solute transport in two-dimensional flow fields influenced by decimeter-scale preferential flowpaths. A hypothetical but geologically plausible network of 10 cm wide channels of high hydraulic conductivity is used to represent the relative preferential flowpaths embedded in an otherwise homogeneous aquifer. When the hydraulic conductivity in the channels is 100 times greater than that in the remaining portion of the aquifer, the calculated concentration distributions under three source configurations all exhibit highly asymmetrical, non-Gaussian patterns. These patterns, with peak concentrations close to the source and extensive spreading downgradient, resemble that observed at the MADE site tracer tests. When the contrast between the channel and nonchannel hydraulic conductivities is reduced to 30:1 from 100:1, the calculated mass distribution curve starts to approach a Gaussian one with the peak concentration near the central portion of the plume. Additional analysis based on a field-scale model demonstrates that the existence of decimeter-scale preferential flowpaths can have potentially far-reaching implications for ground water remediation. Failure to account for them in numerical simulation could lead to overestimation of the effectiveness of the remedial measure under consideration.

Numerical simulation of a natural attenuation experiment with a petroleum hydrocarbon NAPL source

A three-dimensional solute transport model with biological reactions is presented for simulating the natural attenuation study (NATS) at the Columbus Air Force Base in eastern Mississippi. NATS consisted of the release of a petroleum-based nonaqueous phase liquid (NAPL) and subsequent monitoring of BTEX (benzene, toluene, ethylbenzene, p-xylene), naphthalene, decane, and bromide in a shallow, unconfined aquifer. Conceptual and mathematical models were developed for NAPL source release, sequential aerobic/anaerobic biodegradation, and sorption during NATS. A multiple species, solute transport code (SEAM3D) was used to simulate fully three-dimensional transport and aerobic, nitrate-reducing, ferrogenic, and methanogenic hydrocarbon biodegradation. Simulation results matched individual BTEX concentration distributions collected five- and nine-months following NAPL release. SEAM3D mass-balance calculations at t = nine months indicated that 49% of the hydrocarbon mass that dissolved into the aqueous phase was consumed by biodegradation, 13% of this mass was sorbed, and the remaining 38% was present in the aqueous phase. Mass calculations at t = nine months further indicated that aerobic biodegradation accounted for the majority of hydrocarbon biodegradation (46% of the biodegraded mass), followed by ferrogenesis (28%), nitrate-reduction (21%), and methanogenesis (5%). Model results were particularly sensitive to the NAPL release rate, the initial ferric iron (Fe[III]) concentration, hydrocarbon utilization rates, initial condition for the anaerobic microbial populations, and dispersivity.