Search USGSSearch

Geology topics

A.M. Shapiro

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

13 recordsLinked to original sources

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

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

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

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

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

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

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

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

Simulation of steady-state flow in three-dimensional fracture networks using the boundary-element method

An efficient method for simulating steady-state flow in three-dimensional fracture networks is formulated with the use of the boundary-element method. The host rock is considered to be impervious, and the fractures can be of any orientation and areal extent. The fractures are treated as surfaces where fluid movement is essentially two-dimensional. Fracture intersections are regarded as one-dimensional fluid conduits. Hence, the three-dimensional geometric characteristics of the fracture geometry is retained in solutions of coupled sets of one- and two-dimentional equations. Use of the boundary-element method to evaluate the fluid responses in the fractures precludes the need to internally discretize the areal extent of the fractures.

Advances in Water Resources