Search USGS⌕ Search

SEARCH · Search USGS

Results for “Transport in Porous Media”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources
USGS1998

High-resolution Monte Carlo simulation of flow and conservative transport in heterogeneous porous media: 2. Transport results

In this, the second of two papers concerned with the use of numerical simulation to examine flow and transport parameters in heterogeneous porous media via Monte Carlo methods, results from the transport aspect of these simulations are reported on. Transport simulations contained herein assume a finite pulse input of conservative tracer, and the numerical technique endeavors to realistically simulate tracer spreading as the cloud moves through a heterogeneous medium. Medium heterogeneity is limited to the hydraulic conductivity field, and generation of this field assumes that the hydraulic-conductivity process is second-order stationary. Methods of estimating cloud moments, and the interpretation of these moments, are discussed. Techniques for estimation of large-time macrodispersivities from cloud second-moment data, and for the approximation of the standard errors associated with these macrodispersivities, are also presented. These moment and macrodispersivity estimation techniques were applied to tracer clouds resulting from transport scenarios generated by specific Monte Carlo simulations. Where feasible, moments and macrodispersivities resulting from the Monte Carlo simulations are compared with first- and second-order perturbation analyses. Some limited results concerning the possible ergodic nature of these simulations, and the presence of non-Gaussian behavior of the mean cloud, are reported on as well.

Water Resources Research↗
USGS2015

Anomalous solute transport in saturated porous media: Relating transport model parameters to electrical and nuclear magnetic resonance properties

The advection-dispersion equation (ADE) fails to describe commonly observed non-Fickian solute transport in saturated porous media, necessitating the use of other models such as the dual-domain mass-transfer (DDMT) model. DDMT model parameters are commonly calibrated via curve fitting, providing little insight into the relation between effective parameters and physical properties of the medium. There is a clear need for material characterization techniques that can provide insight into the geometry and connectedness of pore spaces related to transport model parameters. Here, we consider proton nuclear magnetic resonance (NMR), direct-current (DC) resistivity, and complex conductivity (CC) measurements for this purpose, and assess these methods using glass beads as a control and two different samples of the zeolite clinoptilolite, a material that demonstrates non-Fickian transport due to intragranular porosity. We estimate DDMT parameters via calibration of a transport model to column-scale solute tracer tests, and compare NMR, DC resistivity, CC results, which reveal that grain size alone does not control transport properties and measured geophysical parameters; rather, volume and arrangement of the pore space play important roles. NMR cannot provide estimates of more-mobile and less-mobile pore volumes in the absence of tracer tests because these estimates depend critically on the selection of a material-dependent and flow-dependent cutoff time. Increased electrical connectedness from DC resistivity measurements are associated with greater mobile pore space determined from transport model calibration. CC was hypothesized to be related to length scales of mass transfer, but the CC response is unrelated to DDMT.

Water Resources Research↗
USGS1987

Analytically-derived sensitivities in one-dimensional models of solute transport in porous media

Analytically-derived sensitivities are presented for parameters in one-dimensional models of solute transport in porous media. Sensitivities were derived by direct differentiation of closed form solutions for each of the odel, and by a time integral method for two of the models. Models are based on the advection-dispersion equation and include adsorption and first-order chemical decay. Boundary conditions considered are: a constant step input of solute, constant flux input of solute, and exponentially decaying input of solute at the upstream boundary. A zero flux is assumed at the downstream boundary. Initial conditions include a constant and spatially varying distribution of solute. One model simulates the mixing of solute in an observation well from individual layers in a multilayer aquifer system. Computer programs produce output files compatible with graphics software in which sensitivities are plotted as a function of either time or space. (USGS)

Open-File Report↗
USGS2010

Further comments on sensitivities, parameter estimation, and sampling design in one-dimensional analysis of solute transport in porous media

Sensitivities of solute concentration to parameters associated with first-order chemical decay, boundary conditions, initial conditions, and multilayer transport are examined in one-dimensional analytical models of transient solute transport in porous media. A sensitivity is a change in solute concentration resulting from a change in a model parameter. Sensitivity analysis is important because minimum information required in regression on chemical data for the estimation of model parameters by regression is expressed in terms of sensitivities. Nonlinear regression models of solute transport were tested on sets of noiseless observations from known models that exceeded the minimum sensitivity information requirements. Results demonstrate that the regression models consistently converged to the correct parameters when the initial sets of parameter values substantially deviated from the correct parameters. On the basis of the sensitivity analysis, several statements may be made about design of sampling for parameter estimation for the models examined: (1) estimation of parameters associated with solute transport in the individual layers of a multilayer system is possible even when solute concentrations in the individual layers are mixed in an observation well; (2) when estimating parameters in a decaying upstream boundary condition, observations are best made late in the passage of the front near a time chosen by adding the inverse of an hypothesized value of the source decay parameter to the estimated mean travel time at a given downstream location; (3) estimation of a first-order chemical decay parameter requires observations to be made late in the passage of the front, preferably near a location corresponding to a travel time of √ 2 times the half-life of the solute; and (4) estimation of a parameter relating to spatial variability in an initial condition requires observations to be made early in time relative to passage of the solute front.

Water Resources Research↗
USGS1993

Surface chemical effects on colloid stability and transport through natural porous media

Surface chemical effects on colloidal stability and transport through porous media were investigated using laboratory column techniques. Approximately 100 nm diameter, spherical, iron oxide particles were synthesized as the mobile colloidal phase. The column packing material was retrieved from a sand and gravel aquifer on Cape Cod, MA. Previous studies have indicated enhanced stability and transport of iron oxide particles due to specific adsorption of some inorganic anions on the iron oxide surface. This phenomenon was further evaluated with an anionic surfactant, sodium dodecyl sulfate. Surfactants constitute a significant mass of the contaminant loading at the Cape Cod site and their presence may contribute to colloidal transport as a significant transport mechanism at the site. Other studies at the site have previously demonstrated the occurrence of this transport mechanism for iron phosphate particles. Photon correlation spectroscopy, micro-electrophoretic mobility, and scanning electron microscopy were used to evaluate particle stability, mobility and size. Adsorption of negatively charged organic and inorganic species onto the surface of the iron oxide particles was shown to significantly enhance particle stability and transport through alterations of the electrokinetic properties of the particle surface. Particle breakthrough generally occurred simultaneously with tritiated water, a conservative tracer. The extent of particle breakthrough was primarily dependent upon colloidal stability and surface charge.

Massachusetts↗
USGS2013

Colloid transport in saturated porous media: Elimination of attachment efficiency in a new colloid transport model

A colloid transport model is introduced that is conceptually simple yet captures the essential features of colloid transport and retention in saturated porous media when colloid retention is dominated by the secondary minimum because an electrostatic barrier inhibits substantial deposition in the primary minimum. This model is based on conventional colloid filtration theory (CFT) but eliminates the empirical concept of attachment efficiency. The colloid deposition rate is computed directly from CFT by assuming all predicted interceptions of colloids by collectors result in at least temporary deposition in the secondary minimum. Also, a new paradigm for colloid re-entrainment based on colloid population heterogeneity is introduced. To accomplish this, the initial colloid population is divided into two fractions. One fraction, by virtue of physiochemical characteristics (e.g., size and charge), will always be re-entrained after capture in a secondary minimum. The remaining fraction of colloids, again as a result of physiochemical characteristics, will be retained “irreversibly” when captured by a secondary minimum. Assuming the dispersion coefficient can be estimated from tracer behavior, this model has only two fitting parameters: (1) the fraction of the initial colloid population that will be retained “irreversibly” upon interception by a secondary minimum, and (2) the rate at which reversibly retained colloids leave the secondary minimum. These two parameters were correlated to the depth of the Derjaguin-Landau-Verwey-Overbeek (DLVO) secondary energy minimum and pore-water velocity, two physical forces that influence colloid transport. Given this correlation, the model serves as a heuristic tool for exploring the influence of physical parameters such as surface potential and fluid velocity on colloid transport.

Water Resources Research↗
USGS2010

Gas transport in unsaturated porous media: The adequacy of Fick's law

The increasing use of natural unsaturated zones as repositories for landfills and disposal sites for hazardous wastes (chemical and radioactive) requires a greater understanding of transport processes in the unsaturated zone. For volatile constituents an important potential transport mechanism is gaseous diffusion. Diffusion, however, cannot be treated as an independent isolated transport mechanism. A complete understanding of multicomponent gas transport in porous media (unsaturated zones) requires a knowledge of Knudsen transport, the molecular and nonequimolar components of diffusive flux, and viscous (pressure driven) flux. The constitutive equations relating these flux components are available from the “dusty gas” model of Mason et al. (1967). This review presents a brief discussion of the underlying principles and interrelationships among each of the above flux mechanisms. Some aspects of these transport mechanisms are, to our knowledge, generally unrecognized in the Earth science literature. The principles underlying the transport mechanisms are illustrated with binary systems; the constitutive equations are then cast in forms thought to be most useful for the study of natural unsaturated zones. The viscous and diffusive fluxes are coupled in the constitutive equations through the Knudsen diffusivities; a knowledge of Knudsen diffusivities is necessary to calculate the viscous component of flux and pressure gradients. The Knudsen diffusivities can be calculated from measurements of the Klinkenberg effect. Two examples are presented showing that in natural systems, very small pressure gradients (1 Pa/m or less) can produce viscous fluxes greater than or equal to diffusive fluxes and that, conversely, pressure gradients of this magnitude can be generated by diffusive processes. The example calculations show that major concentration gradients can be developed for stagnant (zero flux, nonreactive) gases. A method is presented for approximating the viscous and diffusive flux components of gases in a multicomponent system from a knowledge of the concentration profiles of stagnant gases. In subsoil environments, argon and nitrogen are considered to be stagnant gases. Fick's laws are essentially, by definition, inadequate to deal with stagnant gases. In the examples presented, the error associated with estimating the total fluxes of nonstagnant gases by Fick's law, relative to stationary coordinates, ranges from a few percent to orders of magnitude.

Reviews of Geophysics↗
USGS1998

High-resolution Monte Carlo simulation of flow and conservative transport in heterogeneous porous media: 1. Methodology and flow results

In this, the first of two papers concerned with the use of numerical simulation to examine flow and transport parameters in heterogeneous porous media via Monte Carlo methods, various aspects of the modelling effort are examined. In particular, the need to save on core memory causes one to use only specific realizations that have certain initial characteristics; in effect, these transport simulations are conditioned by these characteristics. Also, the need to independently estimate length scales for the generated fields is discussed. The statistical uniformity of the flow field is investigated by plotting the variance of the seepage velocity for vector components in the x , y , and z directions. Finally, specific features of the velocity field itself are illuminated in this first paper. In particular, these data give one the opportunity to investigate the effective hydraulic conductivity in a flow field which is approximately statistically uniform; comparisons are made with first- and second-order perturbation analyses. The mean cloud velocity is examined to ascertain whether it is identical to the mean seepage velocity of the model. Finally, the variance in the cloud centroid velocity is examined for the effect of source size and differing strengths of local transverse dispersion.

Water Resources Research↗
USGS2007

Human enteric viruses in groundwater from a confined bedrock aquifer

Confined aquifers are overlain by low-permeability aquitards that are commonly assumed to protect underlying aquifers from microbial contaminants. However, empirical data on microbial contamination beneath aquitards is limited. This study determined the occurrence of human pathogenic viruses in well water from a deep sandstone aquifer confined by a regionally extensive shale aquitard. Three public water-supply wells were each sampled 10 times over 15 months. Samples were analyzed by reverse transcription−polymerase chain reaction (RT-PCR) for several virus groups and by cell culture for infectious enteroviruses. Seven of 30 samples were positive by RT-PCR for enteroviruses; one of these was positive for infectious echovirus 18. The virus-positive samples were collected from two wells cased through the aquitard, indicating the viruses were present in the confined aquifer. Samples from the same wells showed atmospheric tritium, indicating water recharged within the past few decades. Hydrogeologic conditions support rapid porous media transport of viruses through the upper sandstone aquifer to the top of the aquitard 61 m below ground surface. Natural fractures in the shale aquitard are one possible virus transport pathway through the aquitard; however, windows, cross-connecting well bores, or imperfect grout seals along well casings also may be involved. Deep confined aquifers can be more vulnerable to contamination by human viruses than commonly believed.

Wisconsin↗
USGS2018

VS2DRTI: Simulating heat and reactive solute transport in variably saturated porous media

Variably saturated groundwater flow, heat transport, and solute transport are important processes in environmental phenomena, such as the natural evolution of water chemistry of aquifers and streams, the storage of radioactive waste in a geologic repository, the contamination of water resources from acid-rock drainage, and the geologic sequestration of carbon dioxide. Up to now, our ability to simulate these processes simultaneously with fully coupled reactive transport models has been limited to complex and often difficult-to-use models. To address the need for a simple and easy-to-use model, the VS2DRTI software package has been developed for simulating water flow, heat transport, and reactive solute transport through variably saturated porous media. The underlying numerical model, VS2DRT, was created by coupling the flow and transport capabilities of the VS2DT and VS2DH models with the equilibrium and kinetic reaction capabilities of PhreeqcRM. Flow capabilities include two-dimensional, constant-density, variably saturated flow; transport capabilities include both heat and multicomponent solute transport; and the reaction capabilities are a complete implementation of geochemical reactions of PHREEQC. The graphical user interface includes a preprocessor for building simulations and a postprocessor for visual display of simulation results. To demonstrate the simulation of multiple processes, the model is applied to a hypothetical example of injection of heated waste water to an aquifer with temperature-dependent cation exchange. VS2DRTI is freely available public domain software.

Groundwater↗
USGS1997

Chemical factors influencing colloid-facilitated transport of contaminants in porous media

The effects of colloids on the transport of two strongly sorbing solutes a hydrophobic organic compound, phenanthrene, and a metal ion, Ni 2+ were studied in sand-packed laboratory columns under different pH and ionic strength conditions. Two types of column experiments were performed as follows: (i) sorption/mobilization experiments where the contaminant was first sorbed in the column under conditions where no colloids were released and mobilized under conditions where colloids were released as a result of ionic strength reduction in the influent; and (ii) transport experiments where the contaminant, dissolved or sorbed on colloids, was injected into columns packed with a strongly sorbing porous medium. In the first type of experiment, contaminant mobilization was significant only when all releasable colloids were flushed from the column. In all other cases, although high colloid particle concentrations were encountered, there was no marked effect on total contaminant concentrations. In the second type of experiment, colloid deposition efficiencies were shown to control the enhancement of transport. The deposition efficiency was a function of the pH (for a high organic content sand) and of the contaminant concentration (for a charged species such as Ni 2+ ).

Environmental Science & Technology↗
USGS2003

The effect of entrapped nonaqueous phase liquids on tracer transport in heterogeneous porous media: Laboratory experiments at the intermediate scale

This work considers the applicability of conservative tracers for detecting high-saturation nonaqueous-phase liquid (NAPL) entrapment in heterogeneous systems. For this purpose, a series of experiments and simulations was performed using a two-dimensional heterogeneous system (10??1.2 m), which represents an intermediate scale between laboratory and field scales. Tracer tests performed prior to injecting the NAPL provide the baseline response of the heterogeneous porous medium. Two NAPL spill experiments were performed and the entrapped-NAPL saturation distribution measured in detail using a gamma-ray attenuation system. Tracer tests following each of the NAPL spills produced breakthrough curves (BTCs) reflecting the impact of entrapped NAPL on conservative transport. To evaluate significance, the impact of NAPL entrapment on the conservative-tracer breakthrough curves was compared to simulated breakthrough curve variability for different realizations of the heterogeneous distribution. Analysis of the results reveals that the NAPL entrapment has a significant impact on the temporal moments of conservative-tracer breakthrough curves. ?? 2003 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology↗
USGS2004

Transport of Cryptosporidium oocysts in porous media: Role of straining and physicochemical filtration

The transport and filtration behavior of Cryptosporidium parvum oocysts in columns packed with quartz sand was systematically examined under repulsive electrostatic conditions. An increase in solution ionic strength resulted in greater oocyst deposition rates despite theoretical predictions of a significant electrostatic energy barrier to deposition. Relatively high deposition rates obtained with both oocysts and polystyrene latex particles of comparable size at low ionic strength (1 mM) suggest that a physical mechanism may play a key role in oocyst removal. Supporting experiments conducted with latex particles of varying sizes, under very low ionic strength conditions where physicochemical filtration is negligible, clearly indicated that physical straining is an important capture mechanism. The results of this study indicate that irregularity of sand grain shape (verified by SEM imaging) contributes considerably to the straining potential of the porous medium. Hence, both straining and physicochemical filtration are expected to control the removal of C. parvum oocysts in settings typical of riverbank filtration, soil infiltration, and slow sand filtration. Because classic colloid filtration theory does not account for removal by straining, these observations have important implications with respect to predictions of oocyst transport.

Environmental Science & Technology↗
USGS2002

Is it More Important to Characterize Heterogeneity or Differences in Hydraulic Conductivity Measurements?

As a first step toward understanding the role of sedimentary structures in flow and transport through porous media, this work deterministically examines how transport simulations compare to observed transport through simple, artificial structures in a laboratory experiment. Small-scale laboratory-measured values of hydraulic conductivity were used to simulate transport in an intermediate-scale (10-m long), two-dimensional, heterogeneous porous medium (s 2 lnK=1.26, mlnK = 4.18, where K is cm hr-1 ). Results were judged based on how well the simulated transport matched observed transport through the tank. Permeameter and column experiments produced laboratory measurements of hydraulic conductivity for each of the five sands used in the intermediate-scale experiments. Despite explicit numerical representation of the heterogeneity, predictions using the laboratory-measured values under-estimated the mean arrival time by as much as 35%. The significance of differences between simulated and observed mean arrival time was investigated by comparing variability of transport predictions using the different measurement methods to that produced by different realizations of the heterogeneous distribution. Results indicate that the variations in measured hydraulic conductivity were more important to transport than variations between realizations of the heterogeneous distribution of hydraulic conductivity.

Conference Paper↗
USGS1998

A note on conservative transport in anisotropic, heterogeneous porous media in the presence of small-amplitude transients

The late-time macrodispersion coefficients are obtained for the case of flow in the presence of a small-scale deterministic transient in a three-dimensional anisotropic, heterogeneous medium. The transient is assumed to affect only the velocity component transverse to the mean flow direction and to take the form of a periodic function. For the case of a highly stratified medium, these late-time macrodispersion coefficients behave largely as the standard coefficients used in the transport equation. Only in the event that the medium is isotropic is it probable that significant deviations from the standard coefficients would occur.

Water Resources Research↗
USGS2015

PHT3D-UZF: A reactive transport model for variably-saturated porous media

A modified version of the MODFLOW/MT3DMS-based reactive transport model PHT3D was developed to extend current reactive transport capabilities to the variably-saturated component of the subsurface system and incorporate diffusive reactive transport of gaseous species. Referred to as PHT3D-UZF, this code incorporates flux terms calculated by MODFLOW's unsaturated-zone flow (UZF1) package. A volume-averaged approach similar to the method used in UZF-MT3DMS was adopted. The PHREEQC-based computation of chemical processes within PHT3D-UZF in combination with the analytical solution method of UZF1 allows for comprehensive reactive transport investigations (i.e., biogeochemical transformations) that jointly involve saturated and unsaturated zone processes. Intended for regional-scale applications, UZF1 simulates downward-only flux within the unsaturated zone. The model was tested by comparing simulation results with those of existing numerical models. The comparison was performed for several benchmark problems that cover a range of important hydrological and reactive transport processes. A 2D simulation scenario was defined to illustrate the geochemical evolution following dewatering in a sandy acid sulfate soil environment. Other potential applications include the simulation of biogeochemical processes in variably-saturated systems that track the transport and fate of agricultural pollutants, nutrients, natural and xenobiotic organic compounds and micropollutants such as pharmaceuticals, as well as the evolution of isotope patterns.

Ground Water↗