Geology topics

Daniel J. Goode

Publications and source records attributed to Daniel J. Goode.

46 records · Page 3Linked to original sources

USGS1999

Hydrologic effects of the Pymatuning earthquake of September 25, 1998, in northwestern Pennsylvania

Within hours after the Pymatuning earthquake of September 25, 1998, in northwestern Pennsylvania, local residents reported wells becoming dry, wells beginning to flow, and the formation of new springs. About 120 household-supply wells reportedly went dry within 3 months after the earthquake. About 80 of these wells were on a ridge between Jamestown and Greenville, where water-level declines of as much as 100 feet were documented. Accompanying the decline in water levels beneath the ridge was an increase in water levels in valley wells of as much as 62 feet. One possible explanation of the observed hydrologic effects is that the earthquake increased the vertical hydraulic conductivity of shales beneath the ridge, which allowed ground water to drain from the hilltops. Computer simulations of ground-water flow beneath the ridge between Jamestown and Greenville indicate that increasing the vertical hydraulic conductivity of shale confining beds about 10 to 60 times from their pre-quake values could cause the general pattern of decreased water levels on hilltops and increased levels in valleys.

Pennsylvania↗

USGS1999

Ground-water system, estimation of aquifer hydraulic properties, and effects of pumping on ground-water flow in Triassic sedimentary rocks in and near Lansdale, Pennsylvania

Ground water in Triassic-age sedimentary fractured-rock aquifers in the area of Lansdale, Pa., is used as drinking water and for industrial supply. In 1979, ground water in the Lansdale area was found to be contaminated with trichloroethylene, tetrachloroethylene, and other man-made organic compounds, and in 1989, the area was placed on the U.S. Environmental Protection Agency's (USEPA) National Priority List as the North Penn Area 6 site. To assist the USEPA in the hydrogeological assessment of the site, the U.S. Geological Survey began a study in 1995 to describe the ground-water system and to determine the effects of changes in the well pumping patterns on the direction of ground-water flow in the Lansdale area. This determination is based on hydrologic and geophysical data collected from 1995-98 and on results of the simulation of the regional ground-water-flow system by use of a numerical model. Correlation of natural-gamma logs indicate that the sedimentary rock beds strike generally northeast and dip at angles less than 30 degrees to the northwest. The ground-water system is confined or semi-confined, even at shallow depths; depth to bedrock commonly is less than 20 feet (6 meters); and depth to water commonly is about 15 to 60 feet (5 to 18 meters) below land surface. Single-well, aquifer-interval-isolation (packer) tests indicate that vertical permeability of the sedimentary rocks is low. Multiple-well aquifer tests indicate that the system is heterogeneous and that flow appears primarily in discrete zones parallel to bedding. Preferred horizontal flow along strike was not observed in the aquifer tests for wells open to the pumped interval. Water levels in wells that are open to the pumped interval, as projected along the dipping stratigraphy, are drawn down more than water levels in wells that do not intersect the pumped interval. A regional potentiometric map based on measured water levels indicates that ground water flows from Lansdale towards discharge areas in three drainages, the Wissahickon, Towamencin, and Neshaminy Creeks. Ground-water flow was simulated for different pumping patterns representing past and current conditions. The three-dimensional numerical flow model (MODFLOW) was automatically calibrated by use of a parameter estimation program (MODFLOWP). Steady-state conditions were assumed for the calibration period of 1996. Model calibration indicates that estimated recharge is 8.2 inches (208 millimeters) and the regional anisotropy ratio for the sedimentary-rock aquifer is about 11 to 1, with permeability greatest along strike. The regional anisotropy is caused by up- and down-dip termination of high-permeability bed-oriented features, which were not explicitly simulated in the regional-scale model. The calibrated flow model was used to compare flow directions and capture zones in Lansdale for conditions corresponding to relatively high pumping rates in 1994 and to lower pumping rates in 1997. Comparison of the 1994 and 1997 simulations indicates that wells pumped at the lower 1997 rates captured less ground water from known sites of contamination than wells pumped at the 1994 rates. Ground-water flow rates away from Lansdale increased as pumpage decreased in 1997. A preliminary evaluation of the relation between ground-water chemistry and conditions favorable for the degradation of chlorinated solvents was based on measurements of dissolved-oxygen concentration and other chemical constituents in water samples from 92 wells. About 18 percent of the samples contained less than or equal to 5 milligrams per liter dissolved oxygen, a concentration that indicates reducing conditions favorable for degradation of chlorinated solvents.

Water-Resources Investigations Report↗

USGS1998

Review of Aquifer Test Results for the Lansdale Area, Montgomery County, Pennsylvania, 1980–95

Aquifer and aquifer-isolation test results in and around North Penn Area 6 Superfund site, Lansdale, Montgomery County, Pennsylvania are reviewed to provide estimated aquifer properties for use in a numerical model of ground-water flow. This review was made to support of remedial action investigations by U.S. Environmental Protection Agency (USEPA), Region III, Philadelphia. The data reviewed are from ﬁles of the U.S. Geological Survey, USEPA, and water companies, and from unpublished consultant reports prepared for USEPA and corporations in the Lansdale area. Tested wells are in fractured sedimentary rocks of the Brunswick Formation, which are Triassic-aged, dipping shales and sandstones. Review procedures include, in some cases, new analyses of drawdown during pumping and recovery by use of analytical models of ﬂow to wells. Estimated aquifer transmissivities (T) range from zero to about 1,300 m 2 /d (meters squared per day); most tests indicate T between 10 and 100 m 2 /d. Aquifer-isolation testing results indicate that most ﬂow enters wells at a few discrete zones, probably fractures or bedding-plane openings. The vertical connection between the zones in a single borehole with multiple producing zones commonly is negligible. This suggests that the formation is vertically anisotropic; the hydraulic conductivity is much larger in the horizontal direction than in the vertical direction. Some evidence of well-ﬁeld-scale horizontal anisotropy exists, with maximum transmissivity aligned with the regional northeast strike of bedding, but this evidence is weak because of the small number of observation wells, particularly wells screened in isolated depth intervals. Analysis of recovery data after constant-pumping-rate aquifer tests and of drawdown during step tests suggests that a signiﬁcant fraction, perhaps as much as 85 percent, of the drawdown in some production wells is due to well loss or skin effects in or very near the pumped well and is not caused by resistance to ﬂow in the surrounding formations.

Pennsylvania↗

USGS1998

Ground-water age and atmospheric tracers: Simulation studies and analysis of field data from the Mirror Lake site, New Hampshire

The use of environmental tracers in characterization of ground-water systems is investigated through mathematical modeling of ground-water age and atmospheric tracer transport, and by a field study at the Mirror Lake site, New Hampshire. Theory is presented for modeling ground-water age using the advective-dispersive transport equation. The transport equation includes a zero-order source of unit strength, corresponding to the rate of aging, and can accommodate matrix diffusion and other exchange processes. The effect of temperature fluctuations and layered soils on transport of atmospheric gases to the water table is investigated using a one-dimensional numerical model of chlorofluorocarbon (CFC-11) transport. The nonlinear relation between temperature and Henry's Law coefficient (reflecting air/water phase partitioning) can cause the apparent recharge temperature to be elevated above the annual mean temperature where the water table is shallow. In addition, fine-grained soils can isolate the air phase in the unsaturated zone from the atmosphere. At the USGS' Mirror Lake, New Hampshire fractured-rock research site CFC concentrations near the water table are depleted where dissolved oxygen is low. CFC-11 and CFC-113 are completely absent under anaerobic conditions, while CFC-12 is as low as one-third of modern concentrations. Anaerobic biodegradation apparently consumes CFC's near the water table at this site. One area of active degradation appears to be associated with streamflow loss to ground water. Soil gas concentrations are generally close to atmospheric levels, although some spatial correlation is observed between depleted concentrations of CFC-11 and CFC-113 in soil gas and water-table samples. Results of unsaturated-zone monitoring indicate that recharge occurs throughout the year in the watershed, even during summer evapotranspiration periods, and that seasonal temperature fluctuations occur as much as 5 meters below land surface. Application of ground-water age and CFC-11 transport models to the large-scale ground-water system at Mirror Lake illustrates the similarities between age and chemical transport. Generally, bedrock porosities required to match observed apparent ages from CFC concentrations are high relative to porosities measured on cores. Although matrix diffusion has no effect on steady-state age, it can significantly reduce CFC concentrations in fractured rock in which the effective porosity is low.

New Hampshire↗

USGS1997

Numerical simulation of ground-water flow through glacial deposits and crystalline bedrock in the Mirror Lake area, Grafton County, New Hampshire

This report documents the development of a computer model to simulate steady-state (long-term average) flow of ground water in the vicinity of Mirror Lake, which lies at the eastern end of the Hubbard Brook valley in central New Hampshire. The 10-km2 study area includes Mirror Lake, the three streams that flow into Mirror Lake, Leeman's Brook, Paradise Brook, and parts of Hubbard Brook and the Pemigewasset River. The topography of the area is characterized by steep hillsides and relatively flat valleys. Major hydrogeologic units include glacial deposits, composed of till containing pockets of sand and gravel, and fractured crystalline bedrock, composed of schist intruded by granite, pegmatite, and lamprophyre. Ground water occurs in both the glacial deposits and bedrock. Precipitation and snowmelt infiltrate to the water table on the hillsides, flow downslope through the saturated glacial deposits and fractured bedrock, and discharge to streams and to Mirror Lake. The model domain includes the glacial deposits, the uppermost 150m of bedrock, Mirror Lake, the layer of organic sediments on the lake bottom, and streams and rivers within the study area. A streamflow routing package was included in the model to simulate baseflow in streams and interaction between streams and ground water. Recharge from precipitation is assumed to be areally uniform, and riparian evapotranspiration along stream banks is assumed negligible. The spatial distribution of hydraulic conductivity is represented by dividing the model domain into several zones, each having uniform hydraulic properties. Local variations in recharge and hydraulic conductivities are ignored; therefore, the simulation results characterize the general ground-water system, not local details of ground-water movement. The model was calibrated using a nonlinear regression method to match hydraulic heads measured in piezometers and wells, and baseflow in three inlet streams to Mirror Lake. Model calibration indicates that recharge from precipitation to the water table is 26 to 28 cm/year. Hydraulic conductivities are 1.7 x 10-6 to 2.7 x 10-6 m/s for glacial deposits, about 3 x 10-7 m/s for bedrock beneath lower hillsides and valleys, and about 6x10-8 m/s for bedrock beneath upper hillsides and hilltops. Analysis of parameter uncertainty indicates that the above values are well constrained, at least within the context of regression analysis. In the regression, several attributes of the ground-water flow model are assumed perfectly known. The hydraulic conductivity for bedrock beneath upper hillsides and hilltops was determined from few data, and additional data are needed to further confirm this result. Model fit was not improved by introducing a 10-to-1 ration of horizontal-to-vertical anisotropy in the hydraulic conductivity of the glacial deposits, or by varying hydraulic conductivity with depth in the modeled part (uppermost 150m) of the bedrock. The calibrated model was used to delineate the Mirror Lake ground-water basin, defined as the volumes of subsurface through which ground water flows from the water table to Mirror Lake or its inlet streams. Results indicate that Mirror Lake and its inlet streams drain an area of ground-water recharge that is about 1.5 times the area of the surface-water basin. The ground-water basin extends far up the hillside on the northwestern part of the study area. Ground water from this area flows at depth under Norris Brook to discharge into Mirror Lake or its inlet streams. As a result, the Mirror Lake ground-water basin extends beneath the adjacent ground-water basin that drains into Norris Brook. Model simulation indicates that approximately 300,000 m3/year of precipitation recharges the Mirror Lake ground-water basin. About half the recharge enters the basin in areas where the simulated water table lies in glacial deposits; the other half enters the basin in areas where the simulated water table lies in be

New Hampshire↗

USGS1996

Direct simulation of groundwater age

A new method is proposed to simulate groundwater age directly, by use of an advection-dispersion transport equation with a distributed zero-order source of unit (1) strength, corresponding to the rate of aging. The dependent variable in the governing equation is the mean age, a mass-weighted average age. The governing equation is derived from residence-time-distribution concepts for the case of steady flow. For the more general case of transient flow, a transient governing equation for age is derived from mass-conservation principles applied to conceptual “age mass.” The age mass is the product of the water mass and its age, and age mass is assumed to be conserved during mixing. Boundary conditions include zero age mass flux across all noflow and inflow boundaries and no age mass dispersive flux across outflow boundaries. For transient-flow conditions, the initial distribution of age must be known. The solution of the governing transport equation yields the spatial distribution of the mean groundwater age and includes diffusion, dispersion, mixing, and exchange processes that typically are considered only through tracer-specific solute transport simulation. Traditional methods have relied on advective transport to predict point values of groundwater travel time and age. The proposed method retains the simplicity and tracer-independence of advection-only models, but incorporates the effects of dispersion and mixing on volume-averaged age. Example simulations of age in two idealized regional aquifer systems, one homogeneous and the other layered, demonstrate the agreement between the proposed method and traditional particle-tracking approaches and illustrate use of the proposed method to determine the effects of diffusion, dispersion, and mixing on groundwater age.

Water Resources Research↗

USGS1993

Concentration history during pumping from a leaky aquifer with stratified initial concentration

Analytical and numerical solutions are employed to examine the concentration history of a dissolved substance in water pumped from a leaky aquifer. Many aquifer systems are characterized by stratification, for example, a sandy layer overlain by a clay layer. To obtain information about separate hydrogeologic units, aquifer pumping tests are often conducted with a well penetrating only one of the layers. When the initial concentration distribution is also stratified (the concentration varies with elevation only), the concentration breakthrough in the pumped well may be interpreted to provide information on aquifer hydraulic and transport properties. To facilitate this interpretation, we present some simple analytical and numerical solutions for limiting cases and illustrate their application to a fractured bedrock/glacial drift aquifer system where the solute of interest is dissolved radon gas. In addition to qualitative information on water source, this method may yield estimates of effective porosity and saturated thickness (or fracture transport aperture) from a single-hole test. Little information about dispersivity is obtained because the measured concentration is not significantly affected by dispersion in the aquifer.

Conference Paper↗

USGS1991

Testing a method-of-characteristics model of three-dimensional solute transport in ground water

A new three-dimensional model of solute transport in groundwater that is based on a widely used two-dimensional method of characteristics model and is coupled to a modular finite-difference flow model is under development. The model's accuracy for ideal aquifers having homogeneous properties, uniform boundary conditions, and steady flow along a grid direction is demonstrated by comparison with conventional analytical solutions. The effect of spatially and temporally variable flow velocities is investigated by comparison with special analytical solutions. To test the performance of the model for typical hydrogeologic conditions, we compare results with those from other models as well as with results from the same model using smaller grid spacings and time steps. This model generally provides accurate results for realistic simulations, and is particularly efficient for advection-dominated transport.

Conference Paper↗

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↗

USGS1988

Groundwater velocity magnitude in radionuclide transport calculations

Analytical solutions have been developed for many conceptual models of solute transport in groundwater (Bear 1979). Although these models usually rely on assumptions too restrictive for accurate description of actual field situations, they are useful in understanding groundwater transport and in evaluating the relative importance of the subsurface processes affecting transport. In addition, these simple models are often used for generic and screening-type analyses of groundwater contamination problems (Kent et al. 1985). For example, the Nuclear Regulatory Commission assesses potential doses resulting from the disposal of very slightly contaminated material in the ground using analytical solutions for one- and two-dimensional groundwater transport (Codell and Schreiber 1979; Codell et al. 1982; Goode et al. 1986). This note presents a method for determining a "worst-case" groundwater velocity value for two conceptual models of decaying radionuclide transport, resulting in maximum calculated point concentration.

Journal of Hydraulic Engineering↗