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 73 records · Page 4Linked to original sources
USGS1996

Experimental investigation and modeling of uranium (VI) transport under variable chemical conditions

The transport of adsorbing and complexing metal ions in porous media was investigated with a series of batch and column experiments and with reactive solute transport modeling. Pulses of solutions containing U(VI) were pumped through columns filled with quartz grains, and the breakthrough of U(VI) was studied as a function of variable solution composition ( p H, total U(VI) concentration, total fluoride concentration, and p H-buffering capacity). Decreasing p H and the formation of nonadsorbing aqueous complexes with fluoride increased U(VI) mobility. A transport simulation with surface complexation model (SCM) parameters estimated from batch experiments was able to predict U(VI) retardation in the column experiments within 30%. SCM parameters were also estimated directly from transport data, using the results of three column experiments collected at different p H and U(VI) pulse concentrations. SCM formulations of varying complexity (multiple surface types and reaction stoichiometries) were tested to examine the trade-off between model simplicity and goodness of fit to breakthrough. A two-site model (weak- and strong-binding sites) with three surface complexation reactions fit these transport data well. With this reaction set the model was able to predict (1) the effects of fluoride complexation on U(VI) retardation at two different p H values and (2) the effects of temporal variability of p H on U(VI) transport caused by low p H buffering. The results illustrate the utility of the SCM approach in modeling the transport of adsorbing inorganic solutes under variable chemical conditions.

Water Resources Research↗
USGS2002

Transport of suspended solids from a karstic to an alluvial aquifer: The role of the karst/alluvium interface

This study focuses on the coupled transport of dissolved constituents and particulates, from their infiltration on a karst plateau to their discharge from a karst spring and their arrival at a well in an alluvial plain. Particulate markers were identified and the transport of solids was characterised in situ in porous and karstic media, based on particle size analyses, SEM, and traces. Transport from the sinkhole to the spring appeared to be dominated by flow through karst: particulate transport was apparently conservative between the two sites, and there was little difference in the overall character of the particle size distribution of the particulates infiltrating the sinkhole and of those discharging from the spring. Qualitatively, the mineralogy of the infiltrating and discharging material was similar, although at the spring an autochthonous contribution from the aquifer was noted (chalk particles eroded from the parent rock by weathering). In contrast, transport between the spring and the well appears to be affected by the overlying alluvium: particles in the water from the well, showed evidence of considerable size-sorting. Additionally, SEM images of the well samples showed the presence of particles originating from the overlying alluvial system; these particles were not found in samples from the sinkhole or the spring. The differences between the particulates discharging from the spring and the well indicate that the water pumped from the alluvial plain is coming from the karst aquifer via the very transmissive, complex geologic interface between the underlying chalk formation and the gravel at the base of the overlying alluvial system. ?? 2002 Elsevier Science B.V. All rights reserved.

Journal of Hydrology↗
USGS2008

Comment on “Macrodispersion in sand-shale sequences” by A. J. Desbarats

Desbarats [1990] used a particle-tracking scheme to investigate the physics of three-dimensional solute transport in aquifers composed of two porous media of different hydraulic conductivities. The spatially heterogeneous fluid velocity was assumed to be the only mechanism of solute movement; local or pore scale dispersion and molecular diffusion were assumed to be negligible. The particle-tracking scheme used by Desbarats consisted of routing particles from node to node in a finite difference grid. In this scheme, the direction of an individual particle is randomly selected and the prob- ability associated with the particle movement in a given direction is proportional to the fluid flux in that direction. The same method was used by Moreno et al. [1988] to investigate advective transport in a variable-aperture planar fracture.

Water Resources Research↗
USGS1999

Stochastic analysis of virus transport in aquifers

A large-scale model of virus transport in aquifers is derived using spectral perturbation analysis. The effects of spatial variability in aquifer hydraulic conductivity and virus transport (attachment, detachment, and inactivation) parameters on large-scale virus transport are evaluated. A stochastic mean model of virus transport is developed by linking a simple system of local-scale free-virus transport and attached-virus conservation equations from the current literature with a random-field representation of aquifer and virus transport properties. The resultant mean equations for free and attached viruses are found to differ considerably from the local-scale equations on which they are based and include effects such as a free-virus effective velocity that is a function of aquifer heterogeneity as well as virus transport parameters. Stochastic mean free-virus breakthrough curves are compared with local model output in order to observe the effects of spatial variability on mean one-dimensional virus transport in three-dimensionally heterogeneous porous media. Significant findings from this theoretical analysis include the following: (1) Stochastic model breakthrough occurs earlier than local model breakthrough, and this effect is most pronounced for the least conductive aquifers studied. (2) A high degree of aquifer heterogeneity can lead to virus breakthrough actually preceding that of a conservative tracer. (3) As the mean hydraulic conductivity is increased, the mean model shows less sensitivity to the variance of the natural-logarithm hydraulic conductivity and mean virus diameter. (4) Incorporation of a heterogeneous colloid filtration term results in higher predicted concentrations than a simple first-order adsorption term for a given mean attachment rate. (5) Incorporation of aquifer heterogeneity leads to a greater range of virus diameters for which significant breakthrough occurs. (6) The mean model is more sensitive to the inactivation rate of viruses associated with solid surfaces than to the inactivation rate of viruses in solution.

Water Resources Research↗
USGS2008

Comment on “Flow and tracer transport in a single fracture: A stochastic model and its relation to some field observations” by L. Moreno et al.

Moreno et al. [1988] (hereinafter referred to as MT) used a particle-tracking scheme to investigate the physics of solute movement in a variable-aperture planar fracture. The spatially heterogeneous fluid velocity was assumed to be the only mechanism of solute movement; local or pore scale dispersion and molecular diffusion were assumed to be negligible. The particle-tracking scheme used by MT consisted of routing particles from node to node in a finite difference grid. In this scheme, the direction of an individual particle is randomly selected and the probability associated with the particle movement in a given direction is proportional to the fluid flux in that direction. The same method was used by Desbarats [1990] to investigate advective transport in aquifers composed of two porous media of different hydraulic conductivities.

Water Resources Research↗
USGS2005

Virus fate and transport during artificial recharge with recycled water

A field‐scale experiment was conducted at a research site using bacterial viruses (bacteriophage) MS2 and PRD1 as surrogates for human viruses, bromide as a conservative tracer, and tertiary‐treated municipal wastewater (recycled water) to investigate the fate and transport of viruses during artificial recharge. Observed virus concentrations were fitted using a mathematical model that simulates virus transport in one‐dimensional, homogeneous, water‐saturated porous media accounting for virus sorption (or filtration), virus inactivation, and time‐dependent source concentration. The fitted time‐dependent clogging rate constants were used to estimate the collision efficiencies for bacteriophage MS2 and PRD1 during vertical fully saturated flow. Furthermore, the corresponding time‐dependent collision efficiencies for both bacteriophage asymptotically reached similar values at the various sampling locations. These results can be used to develop an optimal management scenario to maximize the amount of recycled water that can be applied to the spreading grounds while still maintaining favorable attachment conditions for virus removal.

Water Resources Research↗
USGS1977

The use of Galerkin finite-element methods to solve mass-transport equations

The partial differential equation that describes the transport and reaction of chemical solutes in porous media was solved using the Galerkin finite-element technique. These finite elements were superimposed over finite-difference cells used to solve the flow equation. Both convection and flow due to hydraulic dispersion were considered. Linear and Hermite cubic approximations (basis functions) provided satisfactory results: however, the linear functions were computationally more efficient for two-dimensional problems. Successive over relaxation (SOR) and iteration techniques using Tchebyschef polynomials were used to solve the sparce matrices generated using the linear and Hermite cubic functions, respectively. Comparisons of the finite-element methods to the finite-difference methods, and to analytical results, indicated that a high degree of accuracy may be obtained using the method outlined. The technique was applied to a field problem involving an aquifer contaminated with chloride, tritium, and strontium-90. (Woodard-USGS)

Water-Resources Investigations Report↗
USGS1983

Front-tracking model for convective transport in flowing ground water

A finite-difference numerical model that simulates the convective transport of water or tracer particles through porous media is described. It can be applied to one- or two-dimensional problems involving either steady-state or transient flow. The model tracks representative water or tracer particles, initially located along specified lines, as they move in response to the groundwater velocity filed. Aquifer properties may be both anisotropic and nonhomogeneous. Included in the report is a listing of the program along with input formats and test problem results. The front-tracking model provides a useful first approximation determining the movement of solutes in an aquifer, particularly in cases when dispersion and dilution is of minor consideration. (USGS)

Water-Resources Investigations Report↗
USGS2018

Simulation of less‐mobile porosity dynamics in contrasting sediment water interface porous media

Considering heterogeneity in porous media pore size and connectivity is essential to predicting reactive solute transport across interfaces. However, exchange with less‐mobile porosity is rarely considered in surface water/groundwater recharge studies. Previous research indicates that a combination of pore‐fluid sampling and geoelectrical measurements can be used to quantify less‐mobile porosity exchange dynamics using the time‐varying relation between fluid and bulk electrical conductivity. For this study, we use macro‐scale (10 s of cm) advection–dispersion solute transport models linked with electrical conduction in COMSOL Multiphysics to explore less‐mobile porosity dynamics in two different types of observed sediment water interface porous media. Modeled sediment textures contrast from strongly layered streambed deposits to poorly sorted lakebed sands and cobbles. During simulated ionic tracer perturbations, a lag between fluid and bulk electrical conductivity, and the resultant hysteresis, is observed for all simulations indicating differential loading of pore spaces with tracer. Less‐mobile exchange parameters are determined graphically from these tracer time series data without the need for inverse numerical model simulation. In both sediment types, effective less‐mobile porosity exchange parameters are variable in response to changes in flow direction and fluid flux. These observed flow‐dependent effects directly impact local less‐mobile residence times and associated contact time for biogeochemical reaction. The simulations indicate that for the sediment textures explored here, less‐mobile porosity exchange is dominated by variable rates of advection through the domain, rather than diffusion of solute, for typical low‐to‐moderate rate (approximately 3–40 cm/day) hyporheic fluid fluxes. Overall, our model‐based results show that less‐mobile porosity may be expected in a range of natural hyporheic sediments and that changes in flowpath orientation and magnitude will impact less‐mobile exchange parameters. These temporal dynamics can be assessed with the geoelectrical experimental tracer method applied at laboratory and field scales.

Hydrological Processes↗
USGS2021

Co-transport of biogenic nano-hydroxyapatite and Pb(II) in saturated sand columns: Controlling factors and stochastic modeling

Biogenic nano-hydroxyapatite (bio-nHAP) has recently gained great interest in many domains, especially in the remediation of heavy metal-contaminated soil, due to its high reactivity, low cost, and eco-friendly nature. The co-transport and reaction of bio-nHAP with Pb(II) in saturated porous media, however, are not well understood. This work investigated the effects of ionic strength (IS), ionic composition (IC), dissolved organic matter (DOM), and flow velocity on transport-reaction dynamics of Pb(II) and bio-nHAP by combining column breakthrough experiments and model simulations. Results showed that the mobility of Pb(II) was significantly enhanced with increasing IS/IC but less affected by flow velocity during the transport-reaction process of bio-nHAP and Pb(II) in the saturated sand column; while the transport of bio-nHAP was restricted by increasing IS/IC but facilitated by increasing velocity. IC, IS, and velocity only slightly affected the reaction kinetics between Pb(II) and bio-nHAP, likely due to the fast reaction rate between Pb(II) and bio-nHAP and precipitation of pyromorphite. The transport dynamics of bio-nHAP and Pb(II) were significantly changed by DOM, and this effect depended strongly on the type of DOM with different molecular weights. Breakthrough curves of Pb(II) and bio-nHAP exhibited apparent “anomalous”, sub-diffusive transport behaviors, which could be well quantified by a novel tempered fractional derivative bimolecular reaction equation (T-FBRE). Our findings highlighted the accurate simulation of the co-transport and reaction of bio-nHAP with Pb(II) using T-FBRE and had a great benefit for risk assessment and remediation strategy development for Pb(II) contaminated soil.

Chemosphere↗
USGS2002

Effect of growth conditions and staining procedure upon the subsurface transport and attachment behaviors of a groundwater protist

The transport and attachment behaviors of Spumella guttula (Kent), a nanoflagellate (protist) found in contaminated and uncontaminated aquifer sediments in Cape Cod, Mass., were assessed in flowthrough and static columns and in a field injection-and-recovery transport experiment involving an array of multilevel samplers. Transport of S. guttula harvested from low-nutrient (10 mg of dissolved organic carbon per liter), slightly acidic, granular (porous) growth media was compared to earlier observations involving nanoflagellates grown in a traditional high-nutrient liquid broth. In contrast to the highly retarded (retardation factor of ∼3) subsurface transport previously reported for S. guttula , the peak concentration of porous-medium-grown S. guttula traveled concomitantly with that of a conservative (bromide) tracer. About one-third of the porous-medium-grown nanoflagellates added to the aquifer were transported at least 2.8 m downgradient, compared to only ∼2% of the broth-grown nanoflagellates. Flowthrough column studies revealed that a vital (hydroethidine [HE]) staining procedure resulted in considerably less attachment (more transport) of S. guttula in aquifer sediments than did a staining-and-fixation procedure involving 4′,6′-diamidino-2-phenylindole (DAPI) and glutaraldehyde. The calculated collision efficiency (∼10 −2 for porous-medium-grown, DAPI-stained nanoflagellates) was comparable to that observed earlier for the indigenous community of unattached groundwater bacteria that serve as prey. The attachment of HE-labeled S. guttula onto aquifer sediment grains was independent of pH (over the range from pH 3 to 9) suggesting a primary attachment mechanism that may be fundamentally different from that of their prey bacteria, which exhibit sharp decreases in fractional attachment with increasing pH. The high degree of mobility of S. guttula in the aquifer sediments has important ecological implications for the protistan community within the temporally changing plume of organic contaminants in the Cape Cod aquifer.

Applied and Environmental Microbiology↗
USGS2018

The dual‐domain porosity apparatus: Characterizing dual porosity at the sediment/water interface

The characterization of pore-space connectivity in porous media at the sediment/water interface is critical to understanding contaminant transport and reactive biogeochemical processes in zones of groundwater and surface-water exchange. Previous in situ studies of dual-domain (i.e., mobile/less-mobile porosity) studies have been limited to solute tracer injections at scales of meters to 100s of meters and subsequent numerical model parameterization using fluid concentration histories. Pairing fine-scale (e.g., sub-meter) geoelectrical measurements with fluid tracer data over time alleviates dependence on flowpath-scale experiments, enabling spatially targeted characterization of shallow sediment/water interface media where biogeochemical reactivity is often high. The Dual-Domain Porosity Apparatus is a field-tested device capable of variable rate-controlled downward flow experiments. The Dual-Domain Porosity Apparatus facilitates meter-scale inference of dual-domain parameters, i.e., mobile/less-mobile exchange rate coefficient and the ratio of less mobile to mobile porosity. The Dual-Domain Porosity Apparatus experimental procedure uses water electrical conductivity as a conservative tracer of differential loading and flushing of pore spaces within the region of measurement. Variable injection rates permit the direct quantification of the flow-dependence of dual-domain parameters, which has been theorized for decades but remains challenging to assess using existing experimental methodologies.

Groundwater↗
USGS1990

Reevaluation of large-scale dispersivities for a waste chloride plume: Effects of transient flow

This paper investigates the effects of transient groundwater flow on dispersion of a waste chloride plume in the basaltic aquifer beneath the Idaho (USA) National Engineering Laboratory. In an early application of numerical modeling techniques to the two-dimensional simulation of field-scale plumes, previous investigators identified longitudinal and transverse dispersivities using an independently calibrated steady-state flow model and matching contours of observed and simulated concentrations. The unusual result of calibrated transverse dispersivity (140 m) being significantly larger than longitudinal dispersivity (90 m) has been attributed to spatial heterogeneity, the fractured nature of the aquifer, and to the use of a two-dimensional model. New calibrations of the solute-transport model are performed on point concentration observations using a flow model incorporating transient recharge conditions that cause significant fluctuations in velocity. Under transient flow conditions, lowest calibration errors are achieved with significantly larger dispersivities than previously hypothesized, and with the longitudinal component larger than the transverse component. Unfortunately, the sensitivity of the model calibration error to dispersivity is low. Incorporating transient flow in this two-dimensional porous-media model does not significantly improve our understanding of the processes controlling chloride transport at this site.

Conference Paper↗
USGS2004

Extent of reaction in open systems with multiple heterogeneous reactions

The familiar batch concept of extent of reaction is reexamined for systems of reactions occurring in open systems. Because species concentrations change as a result of transport processes as well as reactions in open systems, the extent of reaction has been less useful in practice in these applications. It is shown that by defining the extent of the equivalent batch reaction and a second contribution to the extent of reaction due to the transport processes, it is possible to treat the description of the dynamics of flow through porous media accompanied by many chemical reactions in a uniform, concise manner. This approach tends to isolate the reaction terms among themselves and away from the model partial differential equations, thereby enabling treatment of large problems involving both equilibrium and kinetically controlled reactions. Implications on the number of coupled partial differential equations necessary to be solved and on numerical algorithms for solving such problems are discussed. Examples provided illustrate the theory applied to solute transport in groundwater flow.

American Institute of Chemical Engineers↗
USGS2010

A solute flux approach to transport in heterogeneous formations: 2. Uncertainty analysis

Uncertainty in the mass flux for advection dominated solute movement in heterogeneous porous media is investigated using the Lagrangian framework developed in paper 1 by Dagan et al. (this issue). Expressions for the covariance of the mass flux and cumulative mass flux are derived as functions of the injection volume and sampling area size relative to the scale of heterogeneity. The result is illustrated for solute advection in three types of heterogeneous porous media: stratified formations, two‐ and three‐dimensional porous media; small perturbation approximation is used for the two‐ and three‐dimensional cases. Variances of the mass flux and cumulative mass flux are evaluated as functions of the injection volume (area) scale versus log‐hydraulic conductivity integral scale. The greatest decrease in coefficient of variation (CV) of the mass flux is for the source scale 1–5 times the hydraulic conductivity integral scale; further increase in the source size decreases CV comparatively less. The variance of the cumulative mass flux (or total discharge) indicates that for the source size of 20 hydraulic conductivity integral scales, the transport conditions are almost ergodic. The present results also indicate that the cumulative mass flux is a relatively robust quantity for describing field‐scale solute transport.

Water Resources Research↗
USGS2008

Electrical characterization of non‐Fickian transport in groundwater and hyporheic systems

Recent work indicates that processes controlling solute mass transfer between mobile and less mobile domains in porous media may be quantified by combining electrical geophysical methods and electrically conductive tracers. Whereas direct geochemical measurements of solute preferentially sample the mobile domain, electrical geophysical methods are sensitive to changes in bulk electrical conductivity (bulk EC) and therefore sample EC in both the mobile and immobile domains. Consequently, the conductivity difference between direct geochemical samples and remotely sensed electrical geophysical measurements may provide an indication of mass transfer rates and mobile and immobile porosities in situ. Here we present (1) an overview of a theoretical framework for determining parameters controlling mass transfer with electrical resistivity in situ; (2) a review of a case study estimating mass transfer processes in a pilot‐scale aquifer storage recovery test; and (3) an example application of this method for estimating mass transfer in watershed settings between streams and the hyporheic corridor. We demonstrate that numerical simulations of electrical resistivity studies of the stream/hyporheic boundary can help constrain volumes and rates of mobile‐immobile mass transfer. We conclude with directions for future research applying electrical geophysics to understand field‐scale transport in aquifer and fluvial systems subject to rate‐limited mass transfer.

Water Resources Research↗
USGS2003

MODFLOW-2000, the U.S. Geological Survey Modular Ground-Water Model--Documentation of the SEAWAT-2000 Version with the Variable-Density Flow Process (VDF) and the Integrated MT3DMS Transport Process (IMT)

SEAWAT-2000 is the latest release of the SEAWAT computer program for simulation of three-dimensional, variable-density, transient ground-water flow in porous media. SEAWAT-2000 was designed by combining a modified version of MODFLOW-2000 and MT3DMS into a single computer program. The code was developed using the MODFLOW-2000 concept of a process, which is defined as ?part of the code that solves a fundamental equation by a specified numerical method.? SEAWAT-2000 contains all of the processes distributed with MODFLOW-2000 and also includes the Variable-Density Flow Process (as an alternative to the constant-density Ground-Water Flow Process) and the Integrated MT3DMS Transport Process. Processes may be active or inactive, depending on simulation objectives; however, not all processes are compatible. For example, the Sensitivity and Parameter Estimation Processes are not compatible with the Variable-Density Flow and Integrated MT3DMS Transport Processes. The SEAWAT-2000 computer code was tested with the common variable-density benchmark problems and also with problems representing evaporation from a salt lake and rotation of immiscible fluids.

Open-File Report↗
USGS2012

Phast4Windows: A 3D graphical user interface for the reactive-transport simulator PHAST

Phast4Windows is a Windows® program for developing and running groundwater-flow and reactive-transport models with the PHAST simulator. This graphical user interface allows definition of grid-independent spatial distributions of model properties—the porous media properties, the initial head and chemistry conditions, boundary conditions, and locations of wells, rivers, drains, and accounting zones—and other parameters necessary for a simulation. Spatial data can be defined without reference to a grid by drawing, by point-by-point definitions, or by importing files, including ArcInfo® shape and raster files. All definitions can be inspected, edited, deleted, moved, copied, and switched from hidden to visible through the data tree of the interface. Model features are visualized in the main panel of the interface, so that it is possible to zoom, pan, and rotate features in three dimensions (3D). PHAST simulates single phase, constant density, saturated groundwater flow under confined or unconfined conditions. Reactions among multiple solutes include mineral equilibria, cation exchange, surface complexation, solid solutions, and general kinetic reactions. The interface can be used to develop and run simple or complex models, and is ideal for use in the classroom, for analysis of laboratory column experiments, and for development of field-scale simulations of geochemical processes and contaminant transport.

Groundwater↗