Determining concentration fields of tracer plumes for layered porous media in flow-tank experiments

 In the laboratory, computer-assisted image analysis provides an accurate and efficient way to monitor tracer experiments. This paper describes the determination of detailed temporal concentration distributions of tracers in a flow-tank experiment by analyzing photographs of plumes of Rhodamine dye through the glass wall of the tank. The methodology developed for this purpose consists of four steps: (1) digitally scanning black and white negatives obtained from photographs of the flow–tank experiment; (2) calibrating and normalizing each digitized image to a standard optical-density scale by determining the relation between the optical density and pixel value for each image; (3) constructing standard curves relating the concentration in an optical density from five experimental runs with predetermined concentrations (2–97 mg/L); and (4) converting the optical density to concentration. The spatial distribution of concentration for two photographs was determined by applying these calibration and conversion procedures to all pixels of the digitized images. This approach provides an efficient way to study patterns of plume evolution and transport mechanisms. Received, March 1998 Revised, September 1998 Accepted, October 1998     

DNA tracers with information capacity and high detection sensitivity tested in groundwater studies

 In order to investigate the usefulness of unique synthetic DNA tracers in groundwater, a field experiment was conducted in Norway. DNA tracers and a sodium-chloride tracer were injected into an aquifer. The transport of DNA molecules was interpreted by comparing with the plume of chloride ions under forced-gradient steady-state flow conditions. Spatial concentration moments described the migration of conservative tracers. Mobility and migration of DNA in groundwater demonstrate that DNA tracers can be detected by using the polymerase chain reaction (PCR) and DNA sequence analysis. The results indicate that DNA tracers can be valuable tools as tracers in groundwater investigations. Received, June 1997 / Revised, July 1998, December 1998 / Accepted, January 1999     

Adaptive POD model reduction for solute transport in heterogeneous porous media

We study the applicability of a model order reduction technique to the solution of transport of passive scalars in homogeneous and heterogeneous porous media. Transport dynamics are modeled through the advection-dispersion equation (ADE) and we employ Proper Orthogonal Decomposition (POD) as a strategy to reduce the computational burden associated with the numerical solution of the ADE. Our application of POD relies on solving the governing ADE for selected times, termed snapshots. The latter are then employed to achieve the desired model order reduction. We introduce a new technique, termed Snapshot Splitting Technique (SST), which allows enriching the dimension of the POD subspace and damping the temporal increase of the modeling error. Coupling SST with a modeling strategy based on alternating over diverse time scales the solution of the full numerical transport model to its reduced counterpart allows extending the benefit of POD over a prolonged temporal window so that the salient features of the process can be captured at a reduced computational cost. The selection of the time scales across which the solution of the full and reduced model are alternated is linked to the Péclet number (P e), representing the interplay between advective and dispersive processes taking place in the system. Thus, the method is adaptive in space and time across the heterogenous structure of the domain through the combined use of POD and SST and by way of alternating the solution of the full and reduced models. We find that the width of the time scale within which the POD-based reduced model solution provides accurate results tends to increase with decreasing P e. This suggests that the effects of local-scale dispersive processes facilitate the POD method to capture the salient features of the system dynamics embedded in the selected snapshots. Since the dimension of the reduced model is much lower than that of the full numerical model, the methodology we propose enables one to accurately simulate transport at a markedly reduced computational cost.     

Modeling microbial processes in porous media

 The incorporation of microbial processes into reactive transport models has generally proceeded along two separate lines of investigation: (1) transport of bacteria as inert colloids in porous media, and (2) the biodegradation of dissolved contaminants by a stationary phase of bacteria. Research over the last decade has indicated that these processes are closely linked. This linkage may occur when a change in metabolic activity alters the attachment/detachment rates of bacteria to surfaces, either promoting or retarding bacterial transport in a groundwater-contaminant plume. Changes in metabolic activity, in turn, are controlled by the time of exposure of the microbes to electron acceptors/donor and other components affecting activity. Similarly, metabolic activity can affect the reversibility of attachment, depending on the residence time of active microbes. Thus, improvements in quantitative analysis of active subsurface biota necessitate direct linkages between substrate availability, metabolic activity, growth, and attachment/detachment rates. This linkage requires both a detailed understanding of the biological processes and robust quantitative representations of these processes that can be tested experimentally. This paper presents an overview of current approaches used to represent physicochemical and biological processes in porous media, along with new conceptual approaches that link metabolic activity with partitioning of the microorganism between the aqueous and solid phases. Received, January 1999 · Revised, June 1999, July 1999 · Accepted, October 1999     

The effect of distribution of iron-oxyhydroxide grain coatings on the transport of bacterial cells in porous media

 Among the demonstrated processes influencing the transport of bacteria through aquifers, the deposition of cells on mineral surfaces is one of the most important. For example, understanding the transport of introduced bacteria through aquifers is essential to designing some in situ bioremediation schemes. The impact of the presence and distribution of Fe(III)-oxyhydroxide-coated sand grains on bacterial transport through porous media was evaluated in column experiments in which bacteria (short rods; 1.2 μm length) were eluted through columns of quartz sand (0.5–0.6 mm in diameter) for several conditions of chemical heterogeneity of mineral substrate. Fe(III)-oxyhydroxide-coated sand was present as 10% of the mass, and it was arranged in three treatments: (1) homogeneously distributed, and present as a discrete layer (2) at the top and (3) at the bottom of 14-cm-long sand columns. A pulse input of 10⁸ cells ml^–1 was introduced in an artificial groundwater solution flowing at 14 cm h^–1 through the column, and eluted cells were counted. Peak breakthrough occurred at 1.0 pore volume. A large proportion of cells were retained; 14.7–15.8% of the cells were recovered after three pore volumes of solution had eluted through clean quartz sand, and only 2.1–4.0% were recovered from the Fe(III)-oxyhydroxide-coated sand mixtures. The three physical arrangements of the chemical heterogeneity resulted in essentially the same breakthrough of cells, indicating that the spatial distribution of iron coating does not affect the transport of bacteria. The results of the column transport experiments, which mimic hydrogeological conditions encountered in field problems, are consistent with our mechanistic understanding of bacterial sorption. Received: 10 April 1996 · Accepted: 17 February 1997     

Comparison of numerical simulation of solute transport with observed experimental data in a silt loam subsoil

 Solute transport experiments were conducted on loamy soils of north-eastern Iowa, USA, and the results were compared with a numerical solution of a classical advection-dispersion transport model developed in this study. Flow experiments in the laboratory on undisturbed soil columns showed a flow rate of water much higher than was estimated from the soil properties and grain-size analysis data, suggesting preferential flow regime in the soil. In contrast, the relative concentration peaks of Cl^– and Br^– in the effluent were only approximately 70% of those predicted by the classical advection-dispersion equation (ADE). In addition, the experimental breakthrough curves (BTCs) showed greater tailings of these ions than the model solution. These observations suggest a loss of solute mass during transport from the dynamic flowing regions to a stagnant, immobile water phase in the soil matrix. Experiments in small disturbed soil columns showed that movement of Cl^– and Br^– is in good agreement with predictions of the classical ADE when the tracers are applied as a continuous source. However, in the case of a pulse source, the BTCs of Cl^– and Br^– matched the model only in the ascending part of the curves. Such variation indicates greater retardation of these ions than that of simulation, probably caused by the decrease in soil permeability due to cation exchange reactions in the soil involving monovalent and divalent cation pairs such as K⁺–Ca²⁺ and K⁺–Mg²⁺. In addition, retardation occurred as a result of the continuous saturation of soil columns which seemed to have caused an expansion of clay minerals, thus resulting in decreased soil permeability. In both the continuous and the pulse-source experiments, K⁺ was not detected in the effluent samples, which seemed to have been lost in exchange reactions and adsorption. Received: 20 November 1998 · Accepted: 19 July 1999     

Relation between "terra rossa" from the Apulia aquifer of Italy and the radon content of groundwater: Experimental results and their applicability to radon occurrence in the aquifer

 The radon-222 (²²²Rn) activity in groundwater of the Apulian karstic aquifer in southern Italy is as great as 500 Becquerel per liter (Bq/L) locally. Normal radium-226 (²²⁶Ra) activity in the limestone and calcareous dolomites of the aquifer is not enough to explain such a high level. Laboratory investigations identified high ²²⁶Ra activity in the "terra rossa," the residuum occupying fissures and cavities in the bedrock, and also the relation between (1) ²²⁶Ra-bearing bedrock and "terra rossa" and (2) ²²²Rn in water. The "terra rossa" is the primary source of the radon in the groundwater. The experimental results show the need to characterize the "terra rossa" of Apulia on the basis of ²²⁶Ra activity and also to study the distribution and variations in ²²²Rn activity over time in the aquifer. Received, January 1996 · Revised, March 1997, September 1997, June 1998 · Accepted, July 1998     

Mobile sediment in an urbanizing karst aquifer: implications for contaminant transport

 Here we investigate geochemical characteristics of sediment in different compartments of a karst aquifer and demonstrate that mobile sediments in a karst aquifer can exhibit a wide range of properties affecting their contaminant transport potential. Sediment samples were collected from surface streams, sinkholes, caves, wells, and springs of a karst aquifer (the Barton Springs portion of the Edwards (Balcones Fault Zone) Aquifer, Central Texas) and their mineralogy, grain-size distribution, organic carbon content, and specific surface area analyzed. Statistical analysis of the sediments separated the sampling sites into three distinct groups: (1) streambeds, sinkholes, and small springs; (2) wells; and (3) caves. Sediments from the primary discharge spring were a mix of these three groups. High organic carbon content and high specific surface area gives some sediments an increased potential to transport contaminants; the volume of these sediments is likely to increase with continued urbanization of the watershed. Received: 13 April 1998 · Accepted: 6 October 1998     

Application of a solute transport model under variable velocity conditions in a conduit flow aquifer: Olalde karst system, Basque Country, Spain

A model based on numerical solutions, which allows for solving the dispersion equation under variable recharge and velocity conditions, is developed to simulate solute transport in conduit flow aquifers during flow recession periods. As an example, the evolution of a tracer in the little known karst conduit that links the sinking stream of Oma valley to the Olalde spring is investigated in the karstic region of Santa Eufemia-Ereñozar (Basque Country, Spain). The model, with different hypothetical structures, allows for obtaining series of tracer breakthrough curves, which are fitted to experimental data using an optimization algorithm. These results, although they can be used to simulate the tracer evolution between the two points considered, do not allow for determining the internal structure and spatial disposition of contributions in the aquifer.     

Long-range atmospheric transport of trace metals: the need for geoscience perspectives

 The scientific literature of the past 25 years indicates a lack of consensus over the relative significance of anthropogenic and natural sources of metals in rural and remote ecosystems, meaning areas located several hundred to several thousand kilometers from industrial point sources. Geoscience perspectives are needed to address the critical information gaps associated with this issue, particularly where information on the significance and variability of naturally occurring metals is required. Uncertainties in source apportionment need to be addressed by improving methods of fingerprinting sources of airborne particles, and establishing identities of individual particles on the basis of characteristic morphologies and chemical composition. Further research is required to obtain robust estimates of total metal emissions to the atmosphere and to the oceans in order to determine the relative contribution from anthropogenic sources. Such research must involve a systematic attempt to obtain spatially and temporally representative data on the natural release of metals at the global scale. Moreover, a clear understanding of all local natural and anthropogenic loadings to an ecosystem is required before firm conclusions can be drawn about the relative importance of long-range atmospheric transport. Received: 31 October 1996 · Accepted: 27 May 1997     

Characterization of solute transport parameters in leach ore: inverse modeling based on column experiments

Heap leaching is essentially a process in which metals are extracted from mine ores with lixiant. For a better understanding and modeling of this process, solute transport parameters are required to characterize the solute transport system of the leach heap. For porous media like leach ores, which contain substantial gravelly particles and have a broad range of particle size distributions, traditional small-scale laboratory experimental apparatus is not appropriate. In this paper, a 2.44 m long, 0.3 m inner diameter column was used for tracer test with boron as the tracer. Tracer tests were conducted for 2 bulk densities (1.92 and 1.62 g/cm³) and 2 irrigation rates (2 and 5 L/ (m²·h⁻¹)). Inverse modeling with two-region transport model using computer code CXTFIT was conducted based on the measured breakthrough curves to estimate the transport parameters. Fitting was focused on three parameters: dispersion coefficient D, partition coefficient β, and mass transfer coefficient ω. The results turned out to fall within reasonable ranges. Sensitivity analysis was conducted for the three parameters and showed that the order of sensitivity is β > ω > D. In addition, scaling of these parameters was discussed and applied to a real scale heap leach to predict the tracer breakthrough.     

Clarification of nonlinear retardation factors for colloid-enhanced transport in porous media

There have been a couple of contaminant retardation factors reported for the three-phase (aqueous, solid, and colloid) groundwater system. However, the retardation factor has often been presented by itself and not incorporated into the relevant transport equation, particularly when derived from the mass fraction approach. This may cause a misunderstanding of the retardation factor especially for the systems where multi-phases exist due to the presence of colloids and/or nonlinear sorption processes are involved. It is, therefore, necessary to clarify the form of the nonlinear retardation factor along with the relevant transport equation in the multi-phase system. Alternative forms of the retardation factor and relevant transport equation for specific conditions are presented in various combinations of the nonlinearity of involved sorption mechanisms. The retardation factors for specific conditions are compared with the ones available in the literature. The results indicate that more caution should be given in applying the retardation factor in order to explore contaminant transport in the multi-phase system where any nonlinear sorption is involved. Finally, presentation of the retardation factor along with the relevant transport equation in this study would help prevent possible misuse of the retardation factor in investigating contaminant transport in the multi-phase system.     

Elastic-wave velocities for porous media with power-law distribution of pore sizes

 A system consisting of an elastic matrix with randomly distributed pores of uniform shapes but with sizes governed by a power-law distribution is investigated in a theoretical way. The velocities of elastic waves propagating through such a system are calculated. The obtained results are that the velocities decrease with increasing total porosity and for a given porosity velocities moderately decrease with increasing absolute value of the exponent D in a power-law distribution of pore sizes. Received: 3 July 1995 / Accepted: 2 August 1996     

Variations in hydraulic conductivity with scale of measurement during aquifer tests in heterogeneous, porous carbonate rocks

 Previous studies have shown that hydraulic conductivity of an aquifer seems to increase as the portion of the aquifer tested increases. To date, such studies have all relied on different methods to determine hydraulic conductivity at each scale of interest, which raises the possibility that the observed increase in hydraulic conductivity is due to the measurement method, not to the scale. This study analyzes hydraulic conductivity with respect to scale during individual aquifer tests in porous, heterogeneous carbonate rocks in southeastern Wisconsin, USA. Results from this study indicate that hydraulic conductivity generally increases during an individual test as the volume of aquifer impacted increases, and the rate of this increase is the same as the rate of increase determined by using different measurement methods. Thus, scale dependence of hydraulic conductivity during single tests does not depend on the method of measurement. This conclusion is supported by 22 of 26 aquifer tests conducted in porous-flow-dominated carbonate units within the aquifer. Instead, scale dependency is probably caused by heterogeneities within the aquifer, a conclusion supported by digital simulation. All of the observed types of hydraulic-conductivity variations with scale during individual aquifer tests can be explained by a conceptual model of a simple heterogeneous aquifer composed of high-conductivity zones within a low-conductivity matrix. Received, January 1997 Revised, August 1997, November 1997 Accepted, November 1997     

Brine transport in porous media: On the use of Von Mises and similarity transformations

In this paper we use a Von Mises transformation to study brine transport in porous media. The model involves mass balance equations for fluid and salt, Darcy's law and an equation of state, relating the salt mass fraction to the fluid density. Application of the Von Mises transformation recasts the model equations into a single nonlinear diffusion equation. A further reduction is possible if the problem admits similarity. This yields a formulation in terms of a boundary value problem for an ordinary differential equation which can be treated by semi‐analytical means. Three specific similarity problems are considered in detail: (i) one‐dimensional, stable displacement of fresh water and brine in a porous column, (ii) flow of fresh water along the surface of a salt rock, (iii) mixing of parallel layers of brine and fresh water. This revised version was published online in July 2006 with corrections to the Cover Date.     

The use of connectivity clusters in stochastic modelling of highly heterogeneous porous media

Stochastic geostatistical techniques are essential tools for groundwater flow and transport modelling in highly heterogeneous media. Typically, these techniques require massive numbers of realizations to accurately simulate the high variability and account for the uncertainty. These massive numbers of realizations imposed several constraints on the stochastic techniques (e.g. increasing the computational effort, limiting the domain size, grid resolution, time step and convergence issues). Understanding the connectivity of the subsurface layers gives an opportunity to overcome these constraints. This research presents a sampling framework to reduce the number of the required Monte Carlo realizations utilizing the connectivity properties of the hydraulic conductivity distributions in a three-dimensional domain. Different geostatistical distributions were tested in this study including exponential distribution with the Turning Bands (TBM) algorithm and spherical distribution using Sequential Gaussian Simulation (SGSIM). It is found that the total connected fraction of the largest clusters and its tortuosity are highly correlated with the percentage of mass arrival and the first arrival quantiles at different control planes. Applying different sampling techniques together with several indicators suggested that a compact sample representing only 10% of the total number of realizations can be used to produce results that are close to the results of the full set of realizations. Also, the proposed sampling techniques specially utilizing the low conductivity clustering show very promising results in terms of matching the full range of realizations. Finally, the size of selected clusters relative to domain size significantly affects transport characteristics and the connectivity indicators.     

Numerical modeling of the flow and transport of radionuclides in heterogeneous porous media

This paper is concerned with numerical methods for the modeling of flow and transport of contaminant in porous media. The numerical methods feature the mixed finite element method over triangles as a solver to the Darcy flow equation and a conservative finite volume scheme for the concentration equation. The convective term is approximated with a Godunov scheme over the dual finite volume mesh, whereas the diffusion–dispersion term is discretized by piecewise linear conforming triangular finite elements. It is shown that the scheme satisfies a discrete maximum principle. Numerical examples demonstrate the effectiveness of the methodology for a coupled system that includes an elliptic equation and a diffusion–convection–reaction equation arising when modeling flow and transport in heterogeneous porous media. The proposed scheme is robust, conservative, efficient, and stable, as confirmed by numerical simulations.      

Detrital transport of metals by glaciers, an example from the Pinchi Mine, central British Columbia

 Abundant cinnabar (HgS) mineralization is associated with the Pinchi Fault in central British Columbia. Two formerly producing mercury mines have been developed on this fault: Pinchi and Bralorne Takla. The mercury content of till (a sediment type directly deposited by glaciers) in the area of this fault is primarily controlled by the occurrence of cinnabar mineralization in bedrock and the direction of ice flow. Cinnabar-bearing bedrock was eroded by glaciers, transported in the direction of ice flow, and deposited "down-ice" from its source. An example of such a dispersal train is documented for the Pinchi Mine area where mercury ore was transported over a distance of 12 km, as measured in the clay-sized fraction (< 0.002 mm) of till, and could have been transported over 24 km according to heavy mineral concentrates (specific gravity >3.3) of this same sediment. Antimony, chromium, and nickel dispersal trains were also detected in the region. These data indicate that natural glacial processes can result in the "mobilization" of metals in the surficial environment, a factor which has to be considered at mine sites in glaciated terrain, where mine reclamation and remediation measures are now required. Received: 31 October 1996 · Accepted: 27 May 1997