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 37 records · Page 2Linked to original sources
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↗
USGS2005

Effect of ferric oxyhydroxide grain coatings on the transport of bacteriophage PRD1 and Cryptosporidium parvum oocysts in saturated porous media

To test the effect of geochemical heterogeneity on microorganism transport in saturated porous media, we measured the removal of two microorganisms, the bacteriophage PRD1 and oocysts of the protozoan parasite Cryptosporidium parvum , in flow-through columns of quartz sand coated by different amounts of a ferric oxyhydroxide. The experiments were conducted over ranges of ferric oxyhydroxide coating fraction of λ = 0−0.12 for PRD1 and from λ = 0−0.32 for the oocysts at pH 5.6−5.8 and 10 - 4 M ionic strength. To determine the effect of pH on the transport of the oocysts, experiments were also conducted over a pH range of 5.7−10.0 at a coating fraction of λ = 0.04. Collision (attachment) efficiencies increased as the fraction of ferric oxyhydroxide coated quartz sand increased, from α = 0.0071 to 0.13 over λ = 0−0.12 for PRD1 and from α = 0.059 to 0.75 over λ = 0−0.32 for the oocysts. Increasing the pH from 5.7 to 10.0 resulted in a decrease in the oocyst collision efficiency as the pH exceeded the expected point of zero charge of the ferric oxyhydroxide coatings. The collision efficiencies correlated very well with the fraction of quartz sand coated by the ferric oxyhydroxide for PRD1 but not as well for the oocysts.

Environmental Science & Technology↗
USGS2010

Transport of reacting solutes in porous media: Relation between mathematical nature of problem formulation and chemical nature of reactions

Examples involving six broad reaction classes show that the nature of transport-affecting chemistry may have a profound effect on the mathematical character of solute transport problem formulation. Substantive mathematical diversity among such formulations is brought about principally by reaction properties that determine whether (1) the reaction can be regarded as being controlled by local chemical equilibria or whether it must be considered as being controlled by kinetics, (2) the reaction is homogeneous or heterogeneous, (3) the reaction is a surface reaction (adsorption, ion exchange) or one of the reactions of classical chemistry (e.g., precipitation, dissolution, oxidation, reduction, complex formation). These properties, as well as the choice of means to describe them, stipulate, for instance, (1) the type of chemical entities for which a formulation's basic, mass-balance equations should be written; (2) the nature of mathematical transformations needed to change the problem's basic equations into operational ones. These and other influences determine such mathematical features of problem formulations as the nature of the operational transport-equation system (e.g., whether it involves algebraic, partial-differential, or integro-partial-differential simultaneous equations), the type of nonlinearities of such a system, and the character of the boundaries (e.g., whether they are stationary or moving). Exploration of the reasons for the dependence of transport mathematics on transport chemistry suggests that many results of this dependence stem from the basic properties of the reactions' chemical-relation (i.e., equilibrium or rate) equations.

Water Resources Research↗
USGS1979

A review of numerical simulation of hydrothermal systems

Many advances in simulating single and two-phase fluid flow and heat transport in porous media have recently been made in conjunction with geothermal energy research. These numerical models reproduce system thermal and pressure behaviour and can be used for heat-transport problems other than those associated with geothermal energy development, such as high-level radioactive waste disposal and heat-storage projects. Although these models are general, additional research is necessary before they can be applied to certain site-specific problems that are concerned with additional processes, such as mass transport and flow in fractured media.

Hydrological Sciences Bulletin↗
USGS2007

Role of chemotaxis in the transport of bacteria through saturated porous media

Populations of chemotactic bacteria are able to sense and respond to chemical gradients in their surroundings and direct their migration toward increasing concentrations of chemicals that they perceive to be beneficial to their survival. It has been suggested that this phenomenon may facilitate bioremediation processes by bringing bacteria into closer proximity to the chemical contaminants that they degrade. To determine the significance of chemotaxis in these processes it is necessary to quantify the magnitude of the response and compare it to other groundwater processes that affect the fate and transport of bacteria. We present a systematic approach toward quantifying the chemotactic response of bacteria in laboratory scale experiments by starting with simple, well-defined systems and gradually increasing their complexity. Swimming properties of individual cells were assessed from trajectories recorded by a tracking microscope. These properties were used to calculate motility and chemotaxis coefficients of bacterial populations in bulk aqueous media which were compared to experimental results of diffusion studies. Then effective values of motility and chemotaxis coefficients in single pores, pore networks and packed columns were analyzed. These were used to estimate the magnitude of the chemotactic response in porous media and to compare with dispersion coefficients reported in the field. This represents a compilation of many studies over a number of years. While there are certainly limitations with this approach for ultimately quantifying motility and chemotaxis in granular aquifer media, it does provide insight into what order of magnitude responses are possible and which characteristics of the bacteria and media are expected to be important.

Advances in Water Resources↗
USGS2010

Behavior of sensitivities in the one-dimensional advection-dispersion equation: Implications for parameter estimation and sampling design

The spatial and temporal variability of sensitivities has a significant impact on parameter estimation and sampling design for studies of solute transport in porous media. Physical insight into the behavior of sensitivities is offered through an analysis of analytically derived sensitivities for the one-dimensional form of the advection-dispersion equation. When parameters are estimated in regression models of one-dimensional transport, the spatial and temporal variability in sensitivities influences variance and covariance of parameter estimates. Several principles account for the observed influence of sensitivities on parameter uncertainty. (1) Information about a physical parameter may be most accurately gained at points in space and time with a high sensitivity to the parameter. (2) As the distance of observation points from the upstream boundary increases, maximum sensitivity to velocity during passage of the solute front increases and the consequent estimate of velocity tends to have lower variance. (3) The frequency of sampling must be “in phase” with the S shape of the dispersion sensitivity curve to yield the most information on dispersion. (4) The sensitivity to the dispersion coefficient is usually at least an order of magnitude less than the sensitivity to velocity. (5) The assumed probability distribution of random error in observations of solute concentration determines the form of the sensitivities. (6) If variance in random error in observations is large, trends in sensitivities of observation points may be obscured by noise and thus have limited value in predicting variance in parameter estimates among designs. (7) Designs that minimize the variance of one parameter may not necessarily minimize the variance of other parameters. (8) The time and space interval over which an observation point is sensitive to a given parameter depends on the actual values of the parameters in the underlying physical system.

Water Resources Research↗
USGS2008

Hydrocarbon lakes on Titan: Distribution and interaction with a porous regolith

Synthetic Aperture Radar (SAR) images of Titan's north polar region reveal quasi‐circular to complex features which are interpreted to be liquid hydrocarbon lakes. We investigate methane transport in Titan's hydrologic cycle using the global distribution of lake features. As of May 2007, the SAR data set covers ∼22% of the surface and indicates multiple lake morphologies which are correlated across the polar region. Lakes are limited to latitudes above 55°N and vary from <10 to more than 100,000 km 2 . The size and location of lakes provide constraints on parameters associated with subsurface transport. Using porous media properties inferred from Huygens probe observations, timescales for flow into and out of observed lakes are shown to be in the tens of years, similar to seasonal cycles. Derived timescales are compared to the time between collocated SAR observations in order to consider the role of subsurface transport in Titan's hydrologic cycle.

Geophysical Research Letters↗
USGS2001

Predictive modeling of flow and transport in a two‐dimensional intermediate‐scale, heterogeneous porous medium

As a first step toward understanding the role of sedimentary structures in flow and transport through porous media, this work deterministically examines how small‐scale laboratory‐measured values of hydraulic conductivity relate to in situ values of simple, artificial structures in an intermediate‐scale (10 m long), two‐dimensional, heterogeneous, laboratory experiment. Results were judged based on how well simulations using measured values of hydraulic conductivities matched measured hydraulic heads, net flow, and transport through the tank. Discrepancies were investigated using sensitivity analysis and nonlinear regression estimates of the in situ hydraulic conductivity that produce the best fit to measured hydraulic heads and net flow. Permeameter and column experiments produced laboratory measurements of hydraulic conductivity for each of the sands used in the intermediate‐scale experiments. Despite explicit numerical representation of the heterogeneity the laboratory‐measured values underestimated net flow by 12–14% and were distinctly smaller than the regression‐estimated values. The significance of differences in measured hydraulic conductivity values 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.

Water Resources Research↗
USGS2015

Simulation of the regional groundwater-flow system of the Menominee Indian Reservation, Wisconsin

A regional, two-dimensional, steady-state groundwater-flow model was developed to simulate the groundwater-flow system and groundwater/surface-water interactions within the Menominee Indian Reservation. The model was developed by the U.S. Geological Survey (USGS), in cooperation with the Menominee Indian Tribe of Wisconsin, to contribute to the fundamental understanding of the region&rsquo;s hydrogeology. The objectives of the regional model were to improve understanding of the groundwater-flow system, including groundwater/surface-water interactions, and to develop a tool suitable for evaluating the effects of potential regional water-management programs. The computer code GFLOW was used because of the ease with which the model can simulate groundwater/surface-water interactions, provide a framework for simulating regional groundwater-flow systems, and be refined in a stepwise fashion to incorporate new data and simulate groundwater-flow patterns at multiple scales. Simulations made with the regional model reproduce groundwater levels and stream base flows representative of recent conditions (1970&ndash;2013) and illustrate groundwater-flow patterns with maps of (1) the simulated water table and groundwater-flow directions, (2) probabilistic areas contributing recharge to high-capacity pumped wells, and (3) estimation of the extent of infiltrated wastewater from treatment lagoons. The groundwater-flow model described in this report simulates the major hydrogeologic features of the modeled area, including surficial unconsolidated aquifers, groundwater/surface-water interactions, and groundwater withdrawals from existing high-capacity production wells. Areas contributing recharge to pumped high-capacity wells on the Menominee Indian Reservation were delineated by tracking simulated water particles from the water table to wells in combination with Monte Carlo techniques, and maps of the probability of capture for each well nest were produced. Groundwater-agebased areas contributing recharge to wells were simulated by using the calibrated set of parameters and porosity values adjusted to account for bias in simulated saturated thickness. Simulations were performed for current (2013) pumping rates. The simulations show a range in sensitivity of the simulated areas contributing recharge to wells given the parameters evaluated through the Monte Carlo analysis. The areas contributing recharge to supply wells for the villages of Zoar and Neopit are long and narrow, with a sharp gradation from high to low probability of capture. The areas contributing recharge to supply wells for Middle Village and the village of Keshena exhibit a sharp gradation from high to low probability over a relatively small area between the well and a local groundwater mound. The highest probability areas contributing recharge to the supply wells for the Villages of Onekewat and Redwing are in the immediate vicinity of the wells. These wells also have an extensive area with low probability for capturing water that is likely due to a locally low hydraulic gradient and the large degree of uncertainty associated with the lakebed resistance parameters that control interaction between groundwater and local lakes. Additional field investigations and associated local model refinements would facilitate further reductions in uncertainty associated with simulated areas contributing recharge to the wells. The likely extent of the Neopit wastewater plume was simulated by using the groundwater-flow model and Monte Carlo techniques to evaluate the sensitivity of predictive simulations to a range of model parameter values. Wastewater infiltrated from the currently operating lagoons flows predominantly south toward Tourtillotte Creek. Some of the infiltrated wastewater is simulated as having a low probability of flowing beneath Tourtillotte Creek to the nearby West Branch Wolf River. Results for the probable extent of the wastewater plume are considered to be qualitative because the method only considers advective flow and does not account for processes affecting contaminant transport in porous media. Therefore, results for the probable extent of the wastewater plume are sensitive to the number of particles used to represent flow from the lagoon and the resolution of a synthetic grid used for the analysis. Nonetheless, it is expected that the qualitative results may be of use for identifying potential downgradient areas of concern that can then be evaluated using the quantitative &ldquo;area contributing recharge to wells&rdquo; method or traditional contaminant-transport simulations.

Wisconsin↗
USGS2010

Multiobjective sampling design for parameter estimation and model discrimination in groundwater solute transport

Sampling design for site characterization studies of solute transport in porous media is formulated as a multiobjective problem. Optimal design of a sampling network is a sequential process in which the next phase of sampling is designed on the basis of all available physical knowledge of the system. Three objectives are considered: model discrimination, parameter estimation, and cost minimization. For the first two objectives, physically based measures of the value of information obtained from a set of observations are specified. In model discrimination, value of information of an observation point is measured in terms of the difference in solute concentration predicted by hypothesized models of transport. Points of greatest difference in predictions can contribute the most information to the discriminatory power of a sampling design. Sensitivity of solute concentration to a change in a parameter contributes information on the relative variance of a parameter estimate. Inclusion of points in a sampling design with high sensitivities to parameters tends to reduce variance in parameter estimates. Cost minimization accounts for both the capital cost of well installation and the operating costs of collection and analysis of field samples. Sensitivities, discrimination information, and well installation and sampling costs are used to form coefficients in the multiobjective problem in which the decision variables are binary (zero/one), each corresponding to the selection of an observation point in time and space. The solution to the multiobjective problem is a noninferior set of designs. To gain insight into effective design strategies, a one-dimensional solute transport problem is hypothesized. Then, an approximation of the noninferior set is found by enumerating 120 designs and evaluating objective functions for each of the designs. Trade-offs between pairs of objectives are demonstrated among the models. The value of an objective function for a given design is shown to correspond to the ability of a design to actually meet an objective.

Water Resources Research↗
USGS1996

Effects of thermal vapor diffusion on seasonal dynamics of water in the unsaturated zone

The response of water in the unsaturated zone to seasonal changes of temperature ( T ) is determined analytically using the theory of nonisothermal water transport in porous media, and the solutions are tested against field observations of moisture potential and bomb fallout isotopic ( 36 Cl and 3 H) concentrations. Seasonally varying land surface temperatures and the resulting subsurface temperature gradients induce thermal vapor diffusion. The annual mean vertical temperature gradient is close to zero; however, the annual mean thermal vapor flux is downward, because the temperature‐dependent vapor diffusion coefficient is larger, on average, during downward diffusion (occurring at high T ) than during upward diffusion (low T ). The annual mean thermal vapor flux is shown to decay exponentially with depth; the depth (about 1 m) at which it decays to e −1 of its surface value is one half of the corresponding decay depth for the amplitude of seasonal temperature changes. This depth‐dependent annual mean flux is effectively a source of water, which must be balanced by a flux divergence associated with other transport processes. In a relatively humid environment the liquid fluxes greatly exceed the thermal vapor fluxes, so such a balance is readily achieved without measurable effect on the dynamics of water in the unsaturated zone. However, if the mean vertical water flux through the unsaturated zone is very small (<1 mm y −1 ), as it may be at many locations in a desert landscape, the thermal vapor flux must be balanced mostly by a matric‐potential‐induced upward flux of water. This return flux may include both vapor and liquid components. Below any near‐surface zone of weather‐related fluctuations of matric potential, maintenance of this upward flux requires an increase with depth in the annual mean matric potential; this theoretical prediction is supported by long‐term field measurements in the Chihuahuan Desert. The analysis also makes predictions, confirmed by the field observations, regarding the seasonal variations of matric potential at a given depth. The conceptual model of unsaturated zone water transport developed here implies the possibility of near‐surface trapping of any aqueous constituent introduced at the surface.

Water Resources Research↗
USGS1993

Fate and transport of bacteria injected into aquifers

Advances in our understanding of the fate and transport of bacteria introduced into aquifers, including the potential use of genetically engineered bacteria for biorestoration, are highlighted by new findings in the following areas: modeling of bacterial attachment during transport through porous media, the long-term survival of a chlorobenzoate-degrading bacterium injected into a contaminated sandy aquifer, and molecular techniques that may be used in tracking genetically engineered bacteria in groundwater environments.

Current Opinion in Biotechnology↗
USGS1986

Numerical simulation of advective-dispersive multisolute transport with sorption, ion exchange and equilibrium chemistry

A model was developed that can simulate the effect of certain chemical and sorption reactions simultaneously among solutes involved in advective-dispersive transport through porous media. The model is based on a methodology that utilizes physical-chemical relationships in the development of the basic solute mass-balance equations; however, the form of these equations allows their solution to be obtained by methods that do not depend on the chemical processes. The chemical environment is governed by the condition of local chemical equilibrium, and may be defined either by the linear sorption of a single species and two soluble complexation reactions which also involve that species, or binary ion exchange and one complexation reaction involving a common ion. Partial differential equations that describe solute mass balance entirely in the liquid phase are developed for each tenad (a chemical entity whose total mass is independent of the reaction process) in terms of their total dissolved concentration. These equations are solved numerically in two dimensions through the modification of an existing groundwater flow/transport computer code. (Author 's abstract)

Water-Resources Investigations Report↗
USGS2010

Use of a reactive gas transport model to determine rates of hydrocarbon biodegradation in unsaturated porous media

A mathematical model is presented that simulates the transport and reaction of any number of gaseous phase constituents (e.g. CO 2 , O 2 , N 2 , and hydrocarbons) in unsaturated porous media. The model was developed as part of a method to determine rates of hydrocarbon biodegradation associated with natural cleansing at petroleum product spill sites. The one-dimensional model can be applied to analyze data from column experiments or from field sites where gas transport in the unsaturated zone is approximately vertical. A coupled, non-Fickian constitutive relation between fluxes and concentration gradients, together with the capability of incorporating heterogeneity with respect to model parameters, results in model applicability over a wide range of experimental and field conditions. When applied in a calibration mode, the model allows for the determination of constituent production/consumption rates as a function of the spatial coordinate. Alternatively, the model can be applied in a predictive mode to obtain the distribution of constituent concentrations and fluxes on the basis of assumed values of model parameters and a biodegradation hypothesis. Data requirements for the model are illustrated by analyzing data from a column experiment designed to determine the aerobic degradation rate of toluene in sediments collected from a gasoline spill site in Galloway Township, New Jersey.

Water Resources Research↗