Performance assessment of NAPL remediation in heterogeneous alluvium

Over the last few years, more than 40 partitioning interwell tracer tests (PITTs) have been conducted at many different sites to measure nonaqueous phase liquid (NAPL) saturations in the subsurface. While the main goal of these PITTs was to estimate the NAPL volume in the subsurface, some were specifically conducted to assess the performance of remedial actions involving NAPL removal. In this paper, we present a quantitative approach to assess the performance of remedial actions to recover NAPL that can be used to assess any NAPL removal technology. It combines the use of PITTs (to estimate the NAPL volume in the swept pore volume between injection and extraction wells of a test area) with the use of several cores to determine the vertical NAPL distribution in the subsurface. We illustrate the effectiveness of such an approach by assessing the performance of a surfactant/foam flood conducted at Hill Air Force Base, UT, to remove a TCE-rich NAPL from alluvium with permeability contrasts as high as one order of magnitude. In addition, we compare the NAPL volumes determined by the PITTs with volumes estimated through geostatistical interpolation of aquifer sediment core data collected with a vertical frequency of 5-10 cm and a lateral borehole spacing of 0.15 m. We demonstrate the use of several innovations including the explicit estimation of not only the errors associated with NAPL volumes and saturations derived from PITTs but also the heterogeneity of the aquifer sediments based upon permeability estimates. Most importantly, we demonstrate the reliability of the     

Stochastic-convective transport with nonlinear reactions and mixing: finite streamtube ensemble formulation for multicomponent reaction systems with intra-streamtube dispersion

An effective streamtube ensemble method is developed to upscale convective-dispersive transport with multicomponent nonlinear reactions in steady nonuniform flow. The transport is cast in terms of a finite ensemble of independent discrete streamtubes that approximate convective transport along macroscopically averaged pathlines and dispersive transport longitudinally as microscopic mixing within streamtubes. The representation of fate and transport via a finite ensemble of effective linear streamtubes, allows the treatment of arbitrarily complex reaction systems involving both homogeneous and heterogeneous reactions, and longitudinal dispersive/diffusive mixing within streamtubes. This allows the use of reactive-transport codes designed to solve such problems in an Eulerian framework, as opposed to reliance on closed-form (convolutional or canonical) expressions for reactive transport in exclusively convective streamtubes. The approach requires both reactive-transport solutions for a representative ensemble of one-dimensional convective-dispersive-reactive streamtubes and the distribution of flux over the streamtube ensemble variants, and it does not allow for lateral mixing between streamtubes. Here, the only ensemble variant is travel time. The discussion details the way that the conventional Eulerian fate and transport model is converted first into an ensemble of transports along three-dimensional streamtubes of unknown geometry, and then to approximate one-dimensional streamtubes that are designed to honor the important global properties of the transport. Conditions under which such an 'equivalent' ensemble of one-dimensional streamtubes are described. The breakthrough curve of a nonreactive tracer in the ensemble is expressed as a combined Volterra-Fredholm integral equation, which serves as the basis for estimation of the distribution of flux over the variant of the ensemble, travel time. Transient convective speed and the effects of errors in flux distributions are described, and the method is applied to a demonstration problem involving nonlinear multicomponent reaction kinetics and strongly nonuniform flow.     

Upscaling heterogeneity in aquifer reactivity via exposure-time concept: forward model

Reactive properties of aquifer solid phase materials play an important role in solute fate and transport in the natural subsurface on time scales ranging from years in contaminant remediation to millennia in dynamics of aqueous geochemistry. Quantitative tools for dealing with the impact of natural heterogeneity in solid phase reactivity on solute fate and transport are limited. Here we describe the use of a structural variable to keep track of solute flux exposure to reactive surfaces. With this approach, we develop a non-reactive tracer model that is useful for determining the signature of multi-scale reactive solid heterogeneity in terms of solute flux distributions at the field scale, given realizations of three-dimensional reactive site density fields. First, a governing Eulerian equation for the non-reactive tracer model is determined by an upscaling technique in which it is found that the exposure time of solution to reactive surface areas evolves via both a macroscopic velocity and a macroscopic dispersion in the artificial dimension of exposure time. Second, we focus on the Lagrangian approach in the context of a streamtube ensemble and demonstrate the use of the distribution of solute flux over the exposure time dimension in modeling two-dimensional transport of a solute undergoing simplified linear reversible reactions, in hypothetical conditions following prior laboratory experiments. The distribution of solute flux over exposure time in a given case is a signature of the impact of heterogeneous aquifer reactivity coupled with a particular physical heterogeneity, boundary conditions, and hydraulic gradient. Rigorous application of this approach in a simulation sense is limited here to linear kinetically controlled reactions.     

Understanding the role of land use in urban stormwater quality management

Urbanisation significantly impacts water environments with increased runoff and the degradation of water quality. The management of quantity impacts are straight forward, but quality impacts are far more complex. Current approaches to safeguard water quality are largely ineffective and guided by entrenched misconceptions with a primary focus on 'end-of-pipe' solutions. The outcomes of a research study presented in the paper, which investigated relationships between water quality and six different land uses offer practical guidance in the planning of future urban developments. In terms of safeguarding water quality, high-density residential development which results in a relatively smaller footprint would be the preferred option. The research study outcomes bring into question a number of fundamental concepts and misconceptions routinely accepted in stormwater quality management. The research findings confirmed the need to move beyond customary structural measures and identified the key role that urban planning can play in safeguarding urban water environments.     

Upscaling heterogeneity in aquifer reactivity via the exposure-time concept: inverse model

A novel inverse technique is proposed to quantitatively characterize macroscopic variability in aquifer reactivity in a Lagrangian representation. Reactivity heterogeneity is expressed in terms of distributions of flux over cumulative time of exposure of the solution to reactive surface area, termed here 'cumulative reactivity'. In cases involving single aqueous species the combined effects of physical and reactivity heterogeneity on reactive solute transport can often be established and further investigated through joint distributions of flux over travel time and cumulative reactivity. The inverse technique requires the breakthrough curve of a passive tracer to determine the distribution of flux over travel time, and additional breakthrough curves of reactive tracers provide additional moments of the distribution of flux over cumulative reactivity given travel time. Thus breakthroughs of one passive and two reactive tracers can provide the mean and variance of the distribution of flux over cumulative reactivity. This Lagrangian characterization is achieved with knowledge of the types of reactive surfaces present, but not their spatial locations. The distributions can subsequently be applied via forward modeling using the same technique to predict breakthrough curves of other solutes undergoing first-order reactions in similar physically and chemically heterogeneous configurations.     

Approximation of a radial diffusion model with a multiple-rate model for hetero-disperse particle mixtures

An innovative method is proposed for approximation of the set of radial diffusion equations governing mass exchange between aqueous bulk phase and intra-particle phase for a hetero-disperse mixture of particles such as those occurring in suspension in surface water, in riverine/estuarine sediment beds, in soils and in aquifer materials. For this purpose the temporal variation of concentration at several uniformly distributed points within a normalized representative particle with spherical, cylindrical or planar shape is fitted with a 2-domain linear reversible mass exchange model. The approximation method is then superposed in order to generalize the model to a hetero-disperse mixture of particles. The method can reduce the computational effort needed in solving the intra-particle mass exchange of a hetero-disperse mixture of particles significantly and also the error due to the approximation is shown to be relatively small. The method is applied to describe desorption batch experiment of 1,2-dichlorobenzene from four different soils with known particle size distributions and it could produce good agreement with experimental data.     

Formulation of a soil-pesticide transport model based on a compartmental approach

A semianalytical soil-pesticide transport model is formulated based on a compartmental approach to determine spatial and temporal variations of pesticide residues across a soil profile. The compartmental model is implemented by drawing an analogy between a series of continuous-flow stirred tank reactors and a soil horizon that consists of multiple perfectly mixed compartments. The analogy is strengthened by exploiting a relation between the compartment series and the conventional convective-dispersive equation (CDE) for vertical transport in the soil. Consequently, the number of compartments in the model formulation is not free, but dictated as a function of transport parameters. The model formulation allows consideration of arbitrary boundary value specifications and also, for some cases, spatially varying initial concentration profiles. Sorption kinetics is represented via a two-site model that involves a linear sorption isotherm and a first-order irreversible sorption or a radial diffusive penetrating model. For these three cases, analysis of the compartmental model allows the resultant concentration profiles to be expressed in terms of the Poisson distribution. When a nonlinear kinetic sorption model is used to simulate the sorption processes, an analytical solution is not found and a numerical approach is required.     

Modeling colloid-facilitated transport of multi-species contaminants in unsaturated porous media

Colloid-facilitated transport has been recognized as a potentially important and overlooked contaminant transport process. In particular, it has been observed that conventional two phase sorption models are often unable to explain transport of highly sorbing compounds in the subsurface appropriately in the presence of colloids. In this study a one-dimensional model for colloid-facilitated transport of chemicals in unsaturated porous media is developed. The model has parts for simulating coupled flow, and colloid transport and dissolved and colloidal contaminant transport. Richards' equation is solved to model unsaturated flow, and the effect of colloid entrapment and release on porosity and hydraulic conductivity of the porous media is incorporated into the model. Both random sequential adsorption and Langmuir approaches have been implemented in the model in order to incorporate the effect of surface jamming. The concept of entrapment of colloids into the air-water interface is used for taking into account the effect of retardation caused due to existence of the air phase. A non-equilibrium sorption approach with options of linear and Langmuir sorption assumptions are implemented that can represent the competition and site saturation effects on sorption of multiple compounds both to the solid matrix and to the colloidal particles. Several demonstration calculations are performed and the conditions in which the non-equilibrium model can be approximated by an equilibrium model are also studied.     

Stochastic-convective transport with nonlinear reaction and mixing: application to intermediate-scale experiments in aerobic biodegradation in saturated porous media

Aerobic biodegradation of benzoate by Pseudomonas cepacia sp. in a saturated heterogeneous porous medium was simulated using the stochastic-convective reaction (SCR) approach. A laboratory flow cell was randomly packed with low permeability silt-size inclusions in a high permeability sand matrix. In the SCR upscaling approach, the characteristics of the flow field are determined by the breakthrough of a conservative tracer. Spatial information on the actual location of the heterogeneities is not used. The mass balance equations governing the nonlinear and multicomponent reactive transport are recast in terms of reactive transports in each of a finite number of discrete streamtubes. The streamtube ensemble members represent transport via a steady constant average velocity per streamtube and a conventional Fickian dispersion term, and their contributions to the observed breakthroughs are determined by flux-averaging the streamtube solute concentrations. The resulting simulations were compared to those from a high-resolution deterministic simulation of the reactive transport, and to alternative ensemble representations involving (i) effective Fickian travel time distribution function, (ii) purely convective streamtube transport, and (iii) streamtube ensemble subset simulations. The results of the SCR simulation compare favorably to that of a sophisticated high-resolution deterministic approach.