Transport of volatile compounds in porous media in the presence of a trapped gas phase

The presence of an immobile gaseous phase in an otherwise-saturated porous medium affects the transport of volatile compounds. The linear theory of partitioning tracers suggests that a volatile tracer introduced into such a system should be retarded with a constant retardation factor. Using high concentrations, however, the saturation of the gaseous phase will change as an effect of the tracer test itself. Competitive gas transfer among all volatile compounds and the change of saturation may lead to tracer concentrations that are temporarily higher than those injected. We analyze the system in the framework of the coherence theory by Helfferich [Soc. Pet. Eng. J. 21 (1) (1981) 51]. The governing equations are formulated as functions of total concentration, i.e., the mass of solutes in all phases per pore volume. Neglecting dispersion and mass-transfer kinetics, we derive the characteristic form of the resulting system of hyperbolic equations. In a system with N volatile compounds, a variation of the concentrations splits up into N waves, each traveling with its own characteristic velocity. If the presence of a gaseous phase is sustained, one wave will be a standing one. We perform numerical model calculations for tracers with various Henry's law coefficients and show that the results agree with the semi-analytical solution obtained by coherence theory.     

Chromium transport, oxidation, and adsorption in manganese-coated sand

We examine how the processes of advection, dispersion, oxidation-reduction, and adsorption combine to affect the transport of chromium through columns packed with pyrolusite (beta-MnO2)-coated sand. We find that beta-MnO2 effectively oxidizes Cr(III) to Cr(VI) and that the extent of oxidation is sensitive to changes in pH, pore water velocity, and influent concentrations of Cr(III). Cr(III) oxidation rates, although initially high, decline well before the supply of beta-MnO2 is depleted, suggesting that a reaction product inhibits the conversion of Cr(III) to Cr(VI). Rate-limited reactions govern the weak adsorption of each chromium species, with Cr(III) adsorption varying directly with pH and Cr(VI) adsorption varying inversely with pH. The breakthrough data on chromium transport can be matched closely by calculations of a simple model that accounts for (1) advective-dispersive transport of Cr(III), Cr(VI), and dissolved oxygen, (2) first-order kinetics adsorption of the reduced and oxidized chromium species, and (3) nonlinear rate-limited oxidation of Cr(III) to Cr(VI). Our work supplements the limited database on the transport of redox-sensitive metals in porous media and provides a means for quantifying the coupled processes that contribute to this transport.     

Biodegradation during contaminant transport in porous media: 4. Impact of microbial lag and bacterial cell growth

Miscible-displacement experiments were conducted to examine the impact of microbial lag and bacterial cell growth on the transport of salicylate, a model hydrocarbon compound. The impacts of these processes were examined separately, as well as jointly, to determine their relative effects on biodegradation dynamics. For each experiment, a column was packed with porous medium that was first inoculated with bacteria that contained the NAH plasmid encoding genes for the degradation of naphthalene and salicylate, and then subjected to a step input of salicylate solution. The transport behavior of salicylate was non-steady for all cases examined, and was clearly influenced by a delay (lag) in the onset of biodegradation. This microbial lag, which was consistent with the results of batch experiments, is attributed to the induction and synthesis of the enzymes required for biodegradation of salicylate. The effect of microbial lag on salicylate transport was eliminated by exposing the column to two successive pulses of salicylate, thereby allowing the cells to acclimate to the carbon source during the first pulse. Elimination of microbial lag effects allowed the impact of bacterial growth on salicylate transport to be quantified, which was accomplished by determining a cell mass balance. Conversely, the impact of microbial lag was further investigated by performing a similar double-pulse experiment under no-growth conditions. Significant cell elution was observed and quantified for all conditions/systems. The results of these experiments allowed us to differentiate the effects associated with microbial lag and growth, two coupled processes whose impacts on the biodegradation and transport of contaminants can be difficult to distinguish.     

Transport modeling of copper and cadmium with linear and nonlinear retardation factors

The predictive accuracy of using the one-dimensional advection–dispersion equation to evaluate the fate and transport of solute in a soil column is usually dependent on the proper determination of chemical retardation factors. Typically, the distribution coefficient (K_d) obtained by fitting the linear sorption isotherm has been extensively used to consider general geochemical reactions on solute transport in a low-concentration range. However, the linear distribution coefficient cannot be adequately utilized to describe the solute fate at a higher concentration level. This study employed the nonlinear equilibrium-controlled sorption parameters to determine the retardation factor used in column leaching experiments. Copper and cadmium transportation in a lateritic silty-clay soil column was examined. Through the explicit finite-difference calculations with a third-order total-variation-diminishing (TVD) numerical solution scheme, all results of the theoretical copper and cadmium breakthrough curves (BTCs) simulated by using the Freundlich nonlinear retardation factors revealed good agreement with the experimental observations.     

Evaluating non-equilibrium solute transport in small soil columns

Displacement studies on leaching of bromide and two pesticides (atrazine and isoproturon) were conducted under unsaturated steady state flow conditions in 24 small undisturbed soil columns (5.7 cm in diameter and 10 cm long) each collected from two sites differing in soil structure and organic carbon content in North Germany. There were large and irregular variabilities in the characteristics of both soils, as well as in the shapes of breakthrough curves (BTCs) of different columns, including some with early breakthrough and increased tailing, qualitatively indicating the presence of preferential flow. It was estimated that one preferential flow column (PFC) at site A, and four at site B, contributed, respectively to 11% and 58% of the accumulated leached fraction and to more than 80% of the maximum observed standard deviation (SD) in the field-scale concentration and mass flux of pesticides at two sites. The bromide BTCs of two sites were analyzed with the equilibrium convection-dispersion equation (CDE) and a non-equilibrium two-region/mobile-immobile model. Transport parameters of these models for individual BTCs were determined using a curve fitting program, CXTFIT, and by the time moment method. For the CDE based equilibrium model, the mean values of retardation factor, R, considered separately for all columns, PFCs or non-preferential flow columns (NPFCs) were comparable for the two methods; significant differences were observed in the values of dispersion coefficients of two sites using the two estimation methods. It was inferred from the estimated parameters of non-equilibrium model that 5-12% of water at site A, and 12% at site B, was immobile during displacement in NPFCs. The corresponding values for PFCs of two sites were much larger, ranging from 25% to 51% by CXTFIT and from 24% to 72% by the moment method, suggesting the role of certain mechanisms other than immobile water in higher degrees of non-equilibrium in these columns. Peclet numbers in PFCs of both sites were consistently smaller than five, indicating the inadequacy of the non-equilibrium model to incorporate the effect of all forms of non-equilibrium in PFCs. Overall, the BTCs of individual NPFCs, PFCs and of field average concentration at the two sites were better reproduced with parameters obtained from CXTFIT than by the moment method. The moment method failed to capture the peak concentrations in PFCs, but tended to describe the desorption and tail branches of BTCs better than the curve fitting approach.     

Mineralogical compositions of aquifer matrix as necessary initial conditions in reactive contaminant transport models

Mineralogical compositions and their spatial distributions are important initial conditions for reactive transport modeling. However, popular Kd-based "reactive" transport models only require contaminant concentrations in the pore fluids as initial conditions, and minerals implicitly represent infinite sources and sinks in these models. That situation results in a general neglect of mineralogical characterization in site investigations. This study uses a coupled multi-component reactive mass transport model to predict the natural attenuation of a ground water plume at a uranium mill tailings site in western USA. Numerous ground water geochemistry data are available at this site, but mineralogical data are sketchy. Even given the well-defined pore fluid chemistry, variations of secondary mineral species and mineral abundances in the aquifer resulted in significantly different modeling outcomes. Results show that the amount of calcite in the aquifer determines the distances of plume migration. The possible presence of jurbanite, an aluminum sulfate phase, can store acidity temporarily but cause more severe contamination on a later date. The surfaces of iron oxyhydroxides can store significant amounts of sulfate and protons and serve as a second source for prolonged contamination. These simulations under field conditions illustrate that mineralogical compositions are an essential requirement for accurate prediction of contaminant fate and transport.     

On the hydro-dispersive equivalence between multi-layered mineral barriers

In the context of municipal solid waste and hazardous waste disposal, the notion of "equivalence" between different barrier designs appears in regulatory documents from several industrialized countries. While in the past, equivalence has been thought of mainly in terms of contaminant travel times, in recent years it has been defined more in terms of the magnitude of a disposal site's potential impact on groundwater resources. This paper presents some original analytical solutions to the problem of contaminant migration through a multi-layered mineral barrier. The solutions account for the two major mechanisms of subsurface contaminant migration, namely, advection and diffusion-dispersion. An example application using the proposed solutions and a numerical model illustrates how one multi-layered mineral barrier can be considered superior to another from a strictly hydro-dispersive viewpoint. The influence of partial saturation of the mineral barrier is investigated using a numerical solution to the Richards equation for unsaturated flow. It is emphasized that conclusions relative to the superiority of one multi-layered barrier, with respect to another, should not only consider hydro-dispersive aspects, but also other processes such as the mechanical and chemical evolutions of the different barrier components. Although such phenomena are poorly addressed by existing models, failure to take them into account, at least in a qualitative fashion, may lead to unconservative conclusions with respect to barrier equivalence.     

Migration of colloids in discretely fractured porous media: effect of colloidal matrix diffusion

Numerical simulations of colloid transport in discretely fractured porous media were performed to investigate the importance of matrix diffusion of colloids as well as the filtration and remobilization of colloidal particles in both the fractures and porous matrix. To achieve this objective a finite element numerical code entitled COLDIFF was developed. The processes that COLDIFF takes into account include advective-dispersive transport of colloids, filtration and remobilization of colloidal particles in both fractures and porous matrix, and diffusive interactions of colloids between the fractures and porous matrix. Three sets of simulations were conducted to examine the importance of parameters and processes controlling colloid migration. First, a sensitivity analysis was performed using a porous block containing a single fracture to determine the relative importance of various phenomenological coefficients on colloid transport. The primary result of the analysis showed that the porosity of the matrix and the process of colloid filtration in fractures play important roles in controlling colloid migration. Second, simulations were performed to replicate and examine the results of a laboratory column study using a fractured shale saprolite. Results of this analysis showed that the filtration of colloidal particles in the porous matrix can greatly affect the tailing of colloid concentrations after the colloid source was removed. Finally, field-scale simulations were performed to examine the effect of matrix porosity, fracture filtration and fracture remobilization on long-term colloid concentration and migration distance. The field scale simulations indicated that matrix diffusion and fracture filtration can significantly reduce colloid migration distance. Results of all three analyses indicated that in environments where porosity is relatively high and colloidal particles are small enough to diffuse out of fractures, the characteristics of the porous matrix that affect colloid transport become more important than those of the fracture network. Because the properties of the fracture network tend to have greater uncertainty due to difficulties in their measurement relative to those of the porous matrix, prediction uncertainties associated with colloid transport in discretely fractured porous media may be reduced.     

Reactive transport modeling of processes controlling the distribution and natural attenuation of phenolic compounds in a deep sandstone aquifer

Reactive solute transport modeling was utilized to evaluate the potential for natural attenuation of a contaminant plume containing phenolic compounds at a chemical producer in the West Midlands, UK. The reactive transport simulations consider microbially mediated biodegradation of the phenolic compounds (phenols, cresols, and xylenols) by multiple electron acceptors. Inorganic reactions including hydrolysis, aqueous complexation, dissolution of primary minerals, formation of secondary mineral phases, and ion exchange are considered. One-dimensional (1D) and three-dimensional (3D) simulations were conducted. Mass balance calculations indicate that biodegradation in the saturated zone has degraded approximately 1-5% of the organic contaminant plume over a time period of 47 years. Simulations indicate that denitrification is the most significant degradation process, accounting for approximately 50% of the organic contaminant removal, followed by sulfate reduction and fermentation reactions, each contributing 15-20%. Aerobic respiration accounts for less than 10% of the observed contaminant removal in the saturated zone. Although concentrations of Fe(III) and Mn(IV) mineral phases are high in the aquifer sediment, reductive dissolution is limited, producing only 5% of the observed mass loss. Mass balance calculations suggest that no more than 20-25% of the observed total inorganic carbon (TIC) was generated from biodegradation reactions in the saturated zone. Simulations indicate that aerobic biodegradation in the unsaturated zone, before the contaminant entered the aquifer, may have produced the majority of the TIC observed in the plume. Because long-term degradation is limited to processes within the saturated zone, use of observed TIC concentrations to predict the future natural attenuation may overestimate contaminant degradation by a factor of 4-5.     

Processes affecting the movement of organochlorine pesticides (OCPs) between soil and air in an industrial site in Turkey

Soil and atmospheric concentrations, dry deposition and soil-air gas exchange of organochlorine pesticides (OCPs) were investigated at an industrial site in Aliaga, Izmir, Turkey. Current-use pesticides, endosulfan and chlorpyrifos, had the highest atmospheric levels in summer and winter. Summertime total (gas + particle) OCP concentrations in air were higher, probably due to increased volatilization at higher temperatures and seasonal local/regional applications of current-use pesticides. Particle deposition fluxes were generally higher in summer than in winter. Overall average dry particle deposition velocity for all the OCPs was 4.9 ± 4.1 cm s⁻¹ (average ± SD). ΣDDXs (sum of p,p′-DDT, p,p′-DDD, and p,p′-DDE) were the most abundant OCPs in Aliaga soils (n = 48), probably due to their heavy historical use and persistence. Calculated fugacity ratios and average net gas fluxes across the soil-air interface indicated volatilization for α-CHL, γ-CHL, heptachlorepoxide, cis-nonachlor, trans-nonachlor, and p,p′-DDT in summer, and for α-CHL, γ-CHL, trans-nonachlor, endosulfan sulfate, and p,p′-DDT in winter. For the remaining OCPs, soil acted as a sink during both seasons. Comparison of the determined fluxes showed that dry particle, gas-phase, and wet deposition are significant OCP input mechanisms to the soil in the study area.