USGS Science⌕ Search

Geology topics

A.M. Shapiro

Publications and source records attributed to A.M. Shapiro.

At least 19 recordsLinked to original sources
USGS2008

In situ estimation of the effective chemical diffusion coefficient of a rock matrix in a fractured aquifer

An in situ method of estimating the effective diffusion coefficient for a chemical constituent that diffuses into the primary porosity of a rock is developed by abruptly changing the concentration of the dissolved constituent in a borehole in contact with the rock matrix and monitoring the time-varying concentration. The experiment was conducted in a borehole completed in mudstone on the campus of the University of the Free State in Bloemfontein, South Africa. Numerous tracer tests were conducted at this site, which left a residual concentration of sodium chloride in boreholes that diffused into the rock matrix over a period of years. Fresh water was introduced into a borehole in contact with the mudstone, and the time-varying increase of chloride was observed by monitoring the electrical conductivity (EC) at various depths in the borehole. Estimates of the effective diffusion coefficient were obtained by interpreting measurements of EC over 34 d. The effective diffusion coefficient at a depth of 36 m was approximately 7.8×10 −6 m 2 /d, but was sensitive to the assumed matrix porosity. The formation factor and mass flux for the mudstone were also estimated from the experiment.

Bloemfontein↗
USGS2008

The use of groundwater age as a calibration target

Groundwater age (or residence time), as estimated on the basis of concentrations of one or more environmental tracers, can provide a useful and independent calibration target for groundwater models. However, concentrations of environmental tracers are affected by the complexities and mixing inherent in groundwater flow through heterogeneous media, especially in the presence of pumping wells. An analysis of flow and age distribution in the Madison aquifer in South Dakota, USA, illustrates the additional benefits and difficulties of using age as a calibration target. Alternative numerical approaches to estimating travel time and age with backward particle tracking are assessed, and the resulting estimates are used to refine estimates of effective porosity and to help assess the adequacy and credibility of the flow model.

Conference Paper↗
USGS2005

Bacterial transport experiments in fractured crystalline bedrock

The efficiency of contaminant biodegradation in ground water depends, in part, on the transport properties of the degrading bacteria. Few data exist concerning the transport of bacteria in saturated bedrock, particularly at the field scale. Bacteria and microsphere tracer experiments were conducted in a fractured crystalline bedrock under forced-gradient conditions over a distance of 36 m. Bacteria isolated from the local ground water were chosen on the basis of physicochemical and physiological differences (shape, cell-wall type, motility), and were differentially stained so that their transport behavior could be compared. No two bacterial strains transported in an identical manner, and microspheres produced distinctly different breakthrough curves than bacteria. Although there was insufficient control in this field experiment to completely separate the effects of bacteria shape, reaction to Gram staining, cell size, and motility on transport efficiency, it was observed that (1) the nonmotile, mutant strain exhibited better fractional recovery than the motile parent strain; (2) Gram-negative rod-shaped bacteria exhibited higher fractional recovery relative to the Gram-positive rod-shaped strain of similar size; and (3) coccoidal (spherical-shaped) bacteria transported better than all but one strain of the rod-shaped bacteria. The field experiment must be interpreted in the context of the specific bacterial strains and ground water environment in which they were conducted, but experimental results suggest that minor differences in the physical properties of bacteria can lead to major differences in transport behavior at the field scale.

Ground Water↗
USGS2005

How good are estimates of transmissivity from slug tests in fractured rock?

Slug tests in fractured rock usually are interpreted with models that assume homogeneous formation properties, even though hydraulic properties of fractures can vary by many orders of magnitude over the length of boreholes. To investigate the impact of heterogeneity on the interpretation of slug tests in fractured rock, slug tests were conducted over large intervals of boreholes in crystalline rock in central New Hampshire, and interpreted using a homogeneous model. The results of the slug tests were then compared with estimates of transmissivity from fluid-injection tests conducted over shorter intervals in the same boreholes. The fluid-injection tests showed transmissivity to vary more than six orders of magnitude over the length of the boreholes; however, the sum of the transmissivities from the fluid-injection tests were within an order of magnitude of the transmissivity estimated from the slug tests. Although the two estimates of transmissivity were within an order of magnitude of each other, the water level responses during the slug tests did not exactly match the responses predicted by the homogeneous model. To investigate the effect of heterogeneity on water level responses during slug tests, a Laplace-transform solution was developed for slug tests conducted in boreholes containing multiple fractures with hydraulic properties that vary over the length of the borehole. A comparison of this solution with the homogeneous model shows no difference between the shape of water level responses in a homogeneous formation and a (layered) heterogeneous formation. Furthermore, the transmissivity estimated using a homogeneous model is within an order of magnitude of the prescribed transmissivity in the heterogeneous model. Thus, differences between responses predicted from a homogeneous model and measured water levels during slug tests can be attributed to phenomena such as nonradial flow in the vicinity of the borehole, and not heterogeneous hydraulic properties over the length of the borehole. The experimental results of this investigation show that even when conditions such as nonradial flow are present in the vicinity of the borehole, interpretations of slug tests using a homogeneous model provided order-of-magnitude estimates of transmissivity in the crystalline rock terrane under consideration.Slug tests in fractured rock are usually interpreted with models that assume homogeneous formation properties, although hydraulic properties of fractures can vary by many order of magnitude over the length of the boreholes. To investigate the impact of heterogeneity on the interpretation of slug tests in fractured rocks, slug tests were conducted over large intervals of boreholes in crystalline rock in central New Hampshire. Tests results were interpreted using a homogeneous model. A Laplace-transform solution was developed for slug tests conducted in boreholes containing multiple fractures with hydraulic properties. Finally, a comparison was made between this solution and the homogeneous model.

Ground Water↗
USGS2005

AIRSLUG: A fortran program for the computation of type curves to estimate transmissivity and storativity from prematurely terminated air-pressurized slug tests

The Fortran code AIRSLUG can be used to generate the type curves needed to analyze the recovery data from prematurely terminated air-pressurized slug tests. These type curves, when used with a graphical software package, enable the engineer or scientist to analyze field tests to estimate transmissivity and storativity. Prematurely terminating the slug test can significantly reduce the overall time needed to conduct the test, especially at low-permeability sites, thus saving time and money.The Fortran code AIRSLUG can be used to generate the type curves needed to analyze the recovery data from prematurely terminated air-pressurized slug tests. These type curves, when used with a graphical software package, enable the engineer or scientist to analyze field tests to estimate transmissivity and storativity. Prematurely terminating the slug test can significantly reduce the overall time needed to conduct the test, especially at low-permeability sites, thus saving time and money.

Ground Water↗
USGS2005

Assessing the vulnerability of a municipal well field to contamination in a karst aquifer

Proposed expansion of extractive lime-rock mines near the Miami-Dade County Northwest well field and Everglades wetland areas has garnered intense scrutiny by government, public, environmental stakeholders, and the media because of concern that mining will increase the risk of pathogen contamination. Rock mines are excavated to the same depth as the well field's primary producing zone. The underlying karst Biscayne aquifer is a triple-porosity system characterized by (1) a matrix of interparticle porosity and separate vug porosity; (2) touching-vug porosity that forms preferred, stratiform passageways; and, less commonly, (3) conduit porosity formed by thin solution pipes, bedding-plane vugs, and cavernous vugs. Existing ground-water flow and particle tracking models do not provide adequate information regarding the ability the aquifer to limit the advective movement of pathogens and other contaminants. Chemical transport and colloidal mobility properties have been delineated using conservative and microsphere-surrogate tracers for Cryptosporidium parvum. Forced-gradient tests were executed by introducing conservative tracers into injection wells located 100 m (328 ft) from a municipal-supply well. Apparent mean advective velocity between the wells is one to two orders of magnitude greater than previously measured. Touching-vug, stratiform flow zones are efficient pathways for tracer movement at the well field. The effective porosity for a continuum model between the point of injection and tracer recovery ranges from 2 to 4 percent and is an order of magnitude smaller than previously assumed. Existing well-field protection zones were established using porosity estimates based on specific yield. The effective, or kinematic, porosity of a Biscayne aquifer continuum model is lower than the total porosity, because high velocities occur along preferential flow paths that result in faster times of travel than can be represented with the ground-water flow equation. Tracer tests indicate that the relative ease of contaminant movement to municipal supply wells is much greater than previously considered.

Environmental & Engineering Geoscience↗
USGS2004

Effect of cell physicochemical characteristics and motility on bacterial transport in groundwater

The influence of physicochemical characteristics and motility on bacterial transport in groundwater were examined in flow-through columns. Four strains of bacteria isolated from a crystalline rock groundwater system were investigated, with carboxylate-modified and amidine-modified latex microspheres and bromide as reference tracers. The bacterial isolates included a gram-positive rod (ML1), a gram-negative motile rod (ML2), a nonmotile mutant of ML2 (ML2m), and a gram-positive coccoid (ML3). Experiments were repeated at two flow velocities, in a glass column packed with glass beads, and in another packed with iron-oxyhydroxide coated glass beads. Bacteria breakthrough curves were interpreted using a transport equation that incorporates a sorption model from microscopic observation of bacterial deposition in flow-cell experiments. The model predicts that bacterial desorption rate will decrease exponentially with the amount of time the cell is attached to the solid surface. Desorption kinetics appeared to influence transport at the lower flow rate, but were not discernable at the higher flow rate. Iron-oxyhydroxide coatings had a lower-than-expected effect on bacterial breakthrough and no effect on the microsphere recovery in the column experiments. Cell wall type and shape also had minor effects on breakthrough. Motility tended to increase the adsorption rate, and decrease the desorption rate. The transport model predicts that at field scale, desorption rate kinetics may be important to the prediction of bacteria transport rates. ?? 2003 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology↗
USGS2003

Interpreting tracer breakthrough tailing from different forced-gradient tracer experiment configurations in fractured bedrock

Conceptual and mathematical models are presented that explain tracer breakthrough tailing in the absence of significant matrix diffusion. Model predictions are compared to field results from radially convergent, weak-dipole, and push-pull tracer experiments conducted in a saturated crystalline bedrock. The models are based upon the assumption that flow is highly channelized, that the mass of tracer in a channel is proportional to the cube of the mean channel aperture, and the mean transport time in the channel is related to the square of the mean channel aperture. These models predict the consistent −2 straight line power law slope observed in breakthrough from radially convergent and weak-dipole tracer experiments and the variable straight line power law slope observed in push-pull tracer experiments with varying injection volumes. The power law breakthrough slope is predicted in the absence of matrix diffusion. A comparison of tracer experiments in which the flow field was reversed to those in which it was not indicates that the apparent dispersion in the breakthrough curve is partially reversible. We hypothesize that the observed breakthrough tailing is due to a combination of local hydrodynamic dispersion, which always increases in the direction of fluid velocity, and heterogeneous advection, which is partially reversed when the flow field is reversed. In spite of our attempt to account for heterogeneous advection using a multipath approach, a much smaller estimate of hydrodynamic dispersivity was obtained from push-pull experiments than from radially convergent or weak dipole experiments. These results suggest that although we can explain breakthrough tailing as an advective phenomenon, we cannot ignore the relationship between hydrodynamic dispersion and flow field geometry at this site. The design of the tracer experiment can severely impact the estimation of hydrodynamic dispersion and matrix diffusion in highly heterogeneous geologic media.

Water Resources Research↗
USGS2002

Crosswell seismic investigation of hydraulically conductive, fracture bedrock near Mirror Lake, New Hampshire

Near Mirror Lake, New Hampshire (USA), hydraulically conductive, fractured bedrock was investigated with the crosswell seismic method to determine whether this method could provide any information about hydraulic conductivity between wells. To this end, crosswell seismic data, acoustic logs from boreholes, image logs from boreholes, and single borehole hydraulic tests were analyzed. The analysis showed that, first, the P-wave velocities from the acoustic logs tended to be higher in schist than they were in granite. (Schist and granite were the dominant rock types). Second, the P-wave velocities from the acoustic logs tended to be low near fractures. Third, the hydraulic conductivity was always low (always less than to 10 −8 m/s) where no fractures intersected the borehole, but the hydraulic conductivity ranged from low to high (from less than to 10 −10 m/s to 10 −4 m/s) where one or more fractures intersected the borehole. Fourth, high hydraulic conductivities were slightly more frequent when the P-wave velocity was low (less than 5200 m/s) than when it was high (greater than or equal to 5200 m/s). The interpretation of this statistical relation was that the fractures tended to increase the hydraulic conductivity and to lower the P-wave velocity. This statistical relation was applied to a velocity tomogram to create a map showing the probability of high hydraulic conductivity; the map was consistent with results from independent hydraulic tests.

New Hampshire↗
USGS2000

Estimating formation properties from early-time oscillatory water levels in a pumped well

Hydrologists often attempt to estimate formation properties from aquifer tests for which only the hydraulic responses in a pumped well are available. Borehole storage, turbulent head losses, and borehole skin, however, can mask the hydraulic behavior of the formation inferred from the water level in the pumped well. Also, in highly permeable formations or in formations at significant depth below land surface, where there is a long column of water in the well casing, oscillatory water levels may arise during the onset of pumping to further mask formation responses in the pumped well. Usually borehole phenomena are confined to the early stages of pumping or recovery, and late-time hydraulic data can be used to estimate formation properties. In many instances, however, early-time hydraulic data provide valuable information about the formation, especially if there are interferences in the late-time data. A mathematical model and its Laplace transform solution that account for inertial influences and turbulent head losses during pumping is developed for the coupled response between the pumped borehole and the formation. The formation is assumed to be homogeneous, isotropic, of infinite areal extent, and uniform thickness, with leakage from an overlying aquifer, and the screened or open interval of the pumped well is assumed to fully penetrate the pumped aquifer. Other mathematical models of aquifer flow can also be coupled with the equations describing turbulent head losses and the inertial effects on the water column in the pumped well. The mathematical model developed in this paper is sufficiently general to consider both underdamped conditions for which oscillations arise, and overdamped conditions for which there are no oscillations. Through numerical inversion of the Laplace transform solution, type curves from the mathematical model are developed to estimate formation properties through comparison with the measured hydraulic response in the pumped well. The mathematical model is applied to estimate formation properties from a singlewell test conducted near Waialua, Oahu, Hawaii. At this site, both the drawdown and recovery showed oscillatory water levels in the pumped well, and a step-drawdown test showed that approximately 86% of the drawdown is attributed to turbulent head losses. Analyses at this site using late-time drawdown data were confounded by the noise present in the measured water levels due primarily to nearby irrigation wells and ocean tides. By analyzing the early-time oscillatory recovery data at the Waialua site, upper and lower bounds were placed on the transmissivity, T, storage coefficient, S, and the leakance of the confining unit, K′/B′. The upper and lower bounds on T differ by a factor of 2. Upper and lower bounds on S and K′/B′ are much larger, because drawdown stabilized relatively quickly after the onset of pumping.

Journal of Hydrology↗
USGS1999

Movement of road salt to a small New Hampshire lake

Runoff of road salt from an interstate highway in New Hampshire has led to contamination of a lake and a stream that flows into the lake, in spite of the construction of a diversion berm to divert road salt runoff out of the lake drainage basin. Chloride concentration in the stream has increased by over an order of magnitude during the 23 yr since the highway was opened, and chloride concentration in the lake has tripled. Road salt moves to the lake primarily via the contaminated stream, which provides 53% of all the chloride to the lake and only 3% of the total streamflow to the lake. The stream receives discharge of salty water froth leakage through the diversion berm. Uncontaminated ground water dilutes the stream downstream of the berm. However, reversals of gradient during summer months, likely caused by transpiration from deciduous trees, result in flow of contaminated stream water into the adjacent ground water along the lowest 40-m reach of the stream. This contaminated ground water then discharges into the lake along a 70-m-wide segment of lake shore. Road salt is pervasive in the bedrock between the highway and the lake, but was not detected at all of the wells in the glacial overburden. Of the 500 m of shoreline that could receive discharge of saly ground water directly from the highway, only a 50-m-long segment appears to be contaminated.

Water, Air, & Soil Pollution↗
USGS1998

Estimating formation properties from early-time recovery in wells subject to turbulent head losses

A mathematical model is developed to interpret the early-time recovering water level following the termination of pumping in wells subject to turbulent head losses. The model assumes that turbulent head losses dissipate immediately when pumping ends. In wells subject to both borehole storage and turbulent head losses, the early-time recovery exhibits a slope equal to 1/2 on log-log plots of the recovery versus time. This half-slope response should not be confused with the half-slope response associated with a linear flow regime during aquifer tests. The presence of a borehole skin due to formation damage or stimulation around the pumped well alters the early-time recovery in wells subject to turbulent head losses and gives the appearance of borehole storage, where the recovery exhibits a unit slope on log-log plots of recovery versus time. Type curves can be used to estimate the formation storafivity from the early-time recovery data. In wells that are suspected of having formation damage or stimulation, the type curves can be used to estimate the 'effective' radius of the pumped well, if an estimate of the formation storativity is available from observation wells or other information. Type curves for a homogeneous and isotropic dual-porosity aquifer are developed and applied to estimate formation properties and the effect of formation stimulation from a single-well test conducted in the Madison limestone near Rapid City, South Dakota.A mathematical model is developed to interpret the early-time recovering water level following the termination of pumping in wells subject to turbulent head losses. The model assumes that turbulent head losses dissipate immediately when pumping ends. In wells subject to both borehole storage and turbulent head losses, the early-time recovery exhibits a slope equal to 1/2 on log-log plots of the recovery versus time. This half-slope response should not be confused with the half-slope response associated with a linear flow regime during aquifer tests. The presence of a borehole skin due to formation damage or stimulation around the pumped well alters the early-time recovery in wells subject to turbulent head losses and gives the appearance of borehole storage, where the recovery exhibits a unit slope on log-log plots of recovery versus time. Type curves can be used to estimate the formation storativity from the early-time recovery data. In wells that are suspected of having formation damage or stimulation, the type curves can be used to estimate the `effective' radius of the pumped well, if an estimate of the formation storativity is available from observation wells or other information. Type curves for a homogeneous and isotropic dual-porosity aquifer are developed and applied to estimate formation properties and the effect of formation stimulation from a single-well test conducted in the Madison limestone near Rapid City, South Dakota.

Journal of Hydrology↗
USGS1991

An exact solution of solute transport by one-dimensional random velocity fields

The problem of one-dimensional transport of passive solute by a random steady velocity field is investigated. This problem is representative of solute movement in porous media, for example, in vertical flow through a horizontally stratified formation of variable porosity with a constant flux at the soil surface. Relating moments of particle travel time and displacement, exact expressions for the advection and dispersion coefficients in the Focker-Planck equation are compared with the perturbation results for large distances. The first- and second-order approximations for the dispersion coefficient are robust for a lognormal velocity field. The mean Lagrangian velocity is the harmonic mean of the Eulerian velocity for large distances. This is an artifact of one-dimensional flow where the continuity equation provides for a divergence free fluid flux, rather than a divergence free fluid velocity. ?? 1991 Springer-Verlag.

Stochastic Hydrology and Hydraulics↗
USGS1990

A comparison of two- and three-dimensional stochastic models of regional solute movement

Recent models of solute movement in porous media that are based on a stochastic description of the porous medium properties have been dedicated primarily to a three-dimensional interpretation of solute movement. In many practical problems, however, it is more convenient and consistent with measuring techniques to consider flow and solute transport as an areal, two-dimensional phenomenon. The physics of solute movement, however, is dependent on the three-dimensional heterogeneity in the formation. A comparison of two- and three-dimensional stochastic interpretations of solute movement in a porous medium having a statistically isotropic hydraulic conductivity field is investigated. To provide an equitable comparison between the two- and three-dimensional analyses, the stochastic properties of the transmissivity are defined in terms of the stochastic properties of the hydraulic conductivity. The variance of the transmissivity is shown to be significantly reduced in comparison to that of the hydraulic conductivity, and the transmissivity is spatially correlated over larger distances. These factors influence the two-dimensional interpretations of solute movement by underestimating the longitudinal and transverse growth of the solute plume in comparison to its description as a three-dimensional phenomenon. Although this analysis is based on small perturbation approximations and the special case of a statistically isotropic hydraulic conductivity field, it casts doubt on the use of a stochastic interpretation of the transmissivity in describing regional scale movement. However, by assuming the transmissivity to be the vertical integration of the hydraulic conductivity field at a given position, the stochastic properties of the hydraulic conductivity can be estimated from the stochastic properties of the transmissivity and applied to obtain a more accurate interpretation of solute movement. ?? 1990 Kluwer Academic Publishers.

Transport in Porous Media↗