Geology topics

Daniel J. Goode

Publications and source records attributed to Daniel J. Goode.

At least 37 records · Page 2Linked to original sources

USGS2011

Water budgets and groundwater volumes for abandoned underground mines in the Western Middle Anthracite Coalfield, Schuylkill, Columbia, and Northumberland Counties, Pennsylvania: Preliminary estimates with identification of data needs

This report, prepared in cooperation with the Pennsylvania Department of Environmental Protection (PaDEP), the Eastern Pennsylvania Coalition for Abandoned Mine Reclamation, and the Dauphin County Conservation District, provides estimates of water budgets and groundwater volumes stored in abandoned underground mines in the Western Middle Anthracite Coalfield, which encompasses an area of 120 square miles in eastern Pennsylvania. The estimates are based on preliminary simulations using a groundwater-flow model and an associated geographic information system that integrates data on the mining features, hydrogeology, and streamflow in the study area. The Mahanoy and Shamokin Creek Basins were the focus of the study because these basins exhibit extensive hydrologic effects and water-quality degradation from the abandoned mines in their headwaters in the Western Middle Anthracite Coalfield. Proposed groundwater withdrawals from the flooded parts of the mines and stream-channel modifications in selected areas have the potential for altering the distribution of groundwater and the interaction between the groundwater and streams in the area. Preliminary three-dimensional, steady-state simulations of groundwater flow by the use of MODFLOW are presented to summarize information on the exchange of groundwater among adjacent mines and to help guide the management of ongoing data collection, reclamation activities, and water-use planning. The conceptual model includes high-permeability mine voids that are connected vertically and horizontally within multicolliery units (MCUs). MCUs were identified on the basis of mine maps, locations of mine discharges, and groundwater levels in the mines measured by PaDEP. The locations and integrity of mine barriers were determined from mine maps and groundwater levels. The permeability of intact barriers is low, reflecting the hydraulic characteristics of unmined host rock and coal. A steady-state model was calibrated to measured groundwater levels and stream base flow, the latter at many locations composed primarily of discharge from mines. Automatic parameter estimation used MODFLOW-2000 with manual adjustments to constrain parameter values to realistic ranges. The calibrated model supports the conceptual model of high-permeability MCUs separated by low-permeability barriers and streamflow losses and gains associated with mine infiltration and discharge. The simulated groundwater levels illustrate low groundwater gradients within an MCU and abrupt changes in water levels between MCUs. The preliminary model results indicate that the primary result of increased pumping from the mine would be reduced discharge from the mine to streams near the pumping wells. The intact barriers limit the spatial extent of mine dewatering. Considering the simulated groundwater levels, depth of mining, and assumed bulk porosity of 11 or 40 percent for the mined seams, the water volume in storage in the mines of the Western Middle Anthracite Coalfield was estimated to range from 60 to 220 billion gallons, respectively. Details of the groundwater-level distribution and the rates of some mine discharges are not simulated well using the preliminary model. Use of the model results should be limited to evaluation of the conceptual model and its simulation using porous-media flow methods, overall water budgets for the Western Middle Anthracite Coalfield, and approximate storage volumes. Model results should not be considered accurate for detailed simulation of flow within a single MCU or individual flooded mine. Although improvements in the model calibration were possible by introducing spatial variability in permeability parameters and adjusting barrier properties, more detailed parameterizations have increased uncertainty because of the limited data set. The preliminary identification of data needs includes continuous streamflow, mine discharge rate, and groundwater levels in the mines and adjacent areas. Data collected whe

Pennsylvania↗

USGS2010

Particle velocity interpolation in block-centered finite difference groundwater flow models

A block-centered, finite difference model of two-dimensional groundwater flow yields velocity values at the midpoints of interfaces between adjacent blocks. Method of characteristics, random walk and particle-tracking models of solute transport require velocities at arbitrary particle locations within the finite difference grid. Particle path lines and travel times are sensitive to the spatial interpolation scheme employed, particularly in heterogeneous aquifers. This paper briefly reviews linear and bilinear interpolation of velocity and introduces a new interpolation scheme. Linear interpolation of velocity is consistent with the numerical solution of the flow equation and preserves discontinuities in velocity caused by abrupt (blocky) changes in transmissivity or hydraulic conductivity. However, linear interpolation yields discontinuous and somewhat unrealistic velocities in homogeneous aquifers. Bilinear interpolation of velocity yields continuous and realistic velocities in homogeneous and smoothly heterogeneous aquifers but does not preserve discontinuities in velocity at abrupt transmissivity boundaries. The new scheme uses potentiometric head gradients and offers improved accuracy for nonuniform flow in heterogeneous aquifers with abrupt changes in transmissivity. The new scheme is equivalent to bilinear interpolation in homogeneous media and is equivalent to linear interpolation where gradients are uniform. Selecting the best interpolation scheme depends, in part, on the conceptualization of aquifer heterogeneity, that is, whether changes in transmissivity occur abruptly or smoothly.

Water Resources Research↗

USGS2010

Apparent dispersion in transient groundwater flow

This paper investigates the effects of large-scale temporal velocity fluctuations, particularly changes in the direction of flow, on solute spreading in a two-dimensional aquifer. Relations for apparent longitudinal and transverse dispersivity are developed through an analytical solution for dispersion in a fluctuating, quasi-steady uniform flow field, in which storativity is zero. For transient flow, spatial moments are evaluated from numerical solutions. Ignored or unknown transients in the direction of flow primarily act to increase the apparent transverse dispersivity because the longitudinal dispersivity is acting in a direction that is not the assumed flow direction. This increase is a function of the angle between the transient flow vector and the assumed steady state flow direction and the ratio of transverse to longitudinal dispersivity. The maximum effect on transverse dispersivity occurs if storativity is assumed to be zero, such that the flow field responds instantly to boundary condition changes.

Water Resources Research↗

USGS2010

Multiple well-shutdown tests and site-scale flow simulation in fractured rocks

A new method was developed for conducting aquifer tests in fractured-rock flow systems that have a pump-and-treat (P&T) operation for containing and removing groundwater contaminants. The method involves temporary shutdown of individual pumps in wells of the P&T system. Conducting aquifer tests in this manner has several advantages, including (1) no additional contaminated water is withdrawn, and (2) hydraulic containment of contaminants remains largely intact because pumping continues at most wells. The well-shutdown test method was applied at the former Naval Air Warfare Center (NAWC), West Trenton, New Jersey, where a P&T operation is designed to contain and remove trichloroethene and its daughter products in the dipping fractured sedimentary rocks underlying the site. The detailed site-scale subsurface geologic stratigraphy, a three-dimensional MODFLOW model, and inverse methods in UCODE_2005 were used to analyze the shutdown tests. In the model, a deterministic method was used for representing the highly heterogeneous hydraulic conductivity distribution and simulations were conducted using an equivalent porous media method. This approach was very successful for simulating the shutdown tests, contrary to a common perception that flow in fractured rocks must be simulated using a stochastic or discrete fracture representation of heterogeneity. Use of inverse methods to simultaneously calibrate the model to the multiple shutdown tests was integral to the effectiveness of the approach.

Groundwater↗

USGS2010

Simulation of Runoff and Reservoir Inflow for Use in a Flood-Analysis Model for the Delaware River, Pennsylvania, New Jersey, and New York, 2004-2006

A model was developed to simulate inflow to reservoirs and watershed runoff to streams during three high-flow events between September 2004 and June 2006 for the main-stem subbasin of the Delaware River draining to Trenton, N.J. The model software is a modified version of the U.S. Geological Survey (USGS) Precipitation-Runoff Modeling System (PRMS), a modular, physically based, distributed-parameter modeling system developed to evaluate the impacts of various combinations of precipitation, climate, and land use on surface-water runoff and general basin hydrology. The PRMS model simulates time periods associated with main-stem flooding that occurred in September 2004, April 2005, and June 2006 and uses both daily and hourly time steps. Output from the PRMS model was formatted for use as inflows to a separately documented reservoir and riverrouting model, the HEC-ResSim model, developed by the U.S. Army Corps of Engineers Hydrologic Engineering Center to evaluate flooding. The models were integrated through a graphical user interface. The study area is the 6,780 square-mile watershed of the Delaware River in the states of Pennsylvania, New Jersey, and New York that drains to Trenton, N.J. A geospatial database was created for use with a geographic information system to assist model discretization, determine land-surface characterization, and estimate model parameters. The USGS National Elevation Dataset at 100-meter resolution, a Digital Elevation Model (DEM), was used for model discretization into streams and hydrologic response units. In addition, geospatial processing was used to estimate initial model parameters from the DEM and other data layers, including land use. The model discretization represents the study area using 869 hydrologic response units and 452 stream segments. The model climate data for point stations were obtained from multiple sources. These sources included daily data for 22 National Weather Service (NWS) Cooperative Climate Station network stations, hourly data for 15 stations from the National Climatic Data Center, hourly data for 1 station from the NWS Middle Atlantic River Forecast Center records, and daily and hourly data for 7 stations operated by the New York City Department of Environmental Protection. The NWS Multisensor Precipitation Estimate data set for 2001-2007 was used for computing daily precipitation for the model and for computing hourly precipitation for storm simulation periods. Calibration of the PRMS model included regression and optimization algorithms, as well as manual adjustments of model parameters. The general goal of the calibration procedure was to minimize the difference between discharge measured at USGS streamgages and the corresponding discharge simulated by the model. Daily streamflow data from 35 USGS streamgages were used in model calibration. The streamflow data represent areas draining from 20.2 to 6,780 square miles. The PRMS model simulates reservoir inflow and watershed runoff for use as input into HECResSim for the purpose of evaluating and comparing the effects of different watershed conditions on main-stem flooding in the Delaware River watershed draining to Trenton, N.J. The PRMS model is useful as a planning tool to simulate the effects of land-use changes and different antecedent conditions on local runoff and reservoir inflow and, as input to the HEC-ResSim model, on flood flows in the main stem of the Delaware River.

Open-File Report↗

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↗

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↗

USGS2007

Contamination in fractured-rock aquifers: Research at the former Naval Air Warfare Center, West Trenton, New Jersey

The U.S. Geological Survey and cooperators are studying chlorinated solvents in a fractured sedimentary rock aquifer underlying the former Naval Air Warfare Center (NAWC), West Trenton, New Jersey. Fractured-rock aquifers are common in many parts of the United States and are highly susceptible to contamination, particularly at industrial sites. Compared to 'unconsolidated' aquifers, there can be much more uncertainty about the direction and rate of contaminant migration and about the processes and factors that control chemical and microbial transformations of contaminants. Research at the NAWC is improving understanding of the transport and fate of chlorinated solvents in fractured-rock aquifers and will compare the effectiveness of different strategies for contaminant remediation.

New Jersey↗

USGS2006

Ground-water contamination near a uranium tailings disposal site in Colorado

Contaminants from uranium tailings disposed of at an active mill in Colorado have seeped into the shallow ground water onsite. This ground water discharges into the Arkansas River Valley through a superposed stream channel cut in the resistant sandstone ridge at the edge of a synclinal basin. In the river valley, seasonal surface-water irrigation has a significant impact on hydrodynamics. Water levels in residential wells fluctuate up to 20 ft and concentrations of uranium, molybdenum, and other contaminants also vary seasonally, with highest concentrations in the Spring, prior to irrigation, and lowest concentrations in the Fall. Results of a simple transient mixing cell model support the hypothesis that lateral ground-water inflow, and not irrigation recharge, is the source of ground-water contamination.

Colorado↗

USGS2005

Characterizing a ground water basin in a New England mountain and valley terrain

A ground water basin is defined as the volume of subsurface through which ground water flows from the water table to a specified discharge location. Delineating the topographically defined surface water basin and extending it vertically downward does not always define the ground water basin. Instead, a ground water basin is more appropriately delineated by tracking ground water flowpaths with a calibrated, three‐dimensional ground water flow model. To determine hydrologic and chemical budgets of the basin, it is also necessary to quantify flow through each hydrogeologic unit in the basin. In particular, partitioning ground water flow through unconsolidated deposits versus bedrock is of significant interest to hillslope hydrologic studies. To address these issues, a model is developed and calibrated to simulate ground water flow through glacial deposits and fractured crystalline bedrock in the vicinity of Mirror Lake, New Hampshire. Tracking of ground water flowpaths suggests that Mirror Lake and its inlet streams drain a ground water recharge area that is about 1.5 times the area of the surface water basin. Calculation of the ground water budget suggests that, of the recharge that enters the Mirror Lake ground water basin, about 40% travels through the basin along flowpaths that stay exclusively in the glacial deposits, and about 60% travels along flowpaths that involve movement in bedrock.

New Hampshire↗

USGS2005

Modeling transport in transient ground-water flow: An unacknowledged approximation

During unsteady or transient ground-water flow, the fluid mass per unit volume of aquifer changes as the potentiometric head changes, and solute transport is affected by this change in fluid storage. Three widely applied numerical models of two-dimensional transport partially account for the effects of transient flow by removing terms corresponding to the fluid continuity equation from the transport equation, resulting in a simpler governing equation. However, fluid-storage terms remaining in the transport equation that change during transient flow are, in certain cases, held constant in time in these models. For the case of increasing heads, this approximation, which is unacknowledged in these models' documentation, leads to transport velocities that are too high, and increased concentration at fluid and solute sources. If heads are dropping in time, computed transport velocities are too low. Using parameters that somewhat exaggerate the effects of this approximation, an example numerical simulation indicates solute travel time error of about 14 percent but only minor errors due to incorrect dilution volume. For horizontal flow and transport models that assume fluid density is constant, the product of porosity and aquifer thickness changes in time: initial porosity times initial thickness plus the change in head times the storage coefficient. This formula reduces to the saturated thickness in unconfined aquifers if porosity is assumed to be constant and equal to specific yield. The computational cost of this more accurate representation is insignificant and is easily incorporated in numerical models of solute transport.

Ground Water↗

USGS2005

Composite recovery type curves in normalized time from Theis' exact solution

Type curves derived from Theis’ exact nonequilibrium well function solution are proposed for graphical estimation of aquifer hydraulic properties, transmissivity (T), and storage coefficient (S), from water‐level recovery data after cessation of a constant‐rate discharge test. Drawdown (on log scale) is plotted versus the ratio of time since pumping stopped to duration of pumping, a normalized time. Under Theis conditions, individual type curves depend on only the dimensionless pumping duration, which depends in turn on S and radial distance from the pumping well. Type curve matching, in contrast to the Theis procedure for pumping data, is performed by shifting only the drawdown axis; the time axis is fixed because it is a relative or normalized time. The match‐point for the drawdown axis is used to compute T, and S is determined from matching the curve shape, which depends on early dimensionless‐time data. Multiple well data can be plotted and matched simultaneously (a composite plot), with drawdown at different radial distances matching different curves. The ratio of dimensionless pumping durations for any two matched curves is equal to one over the squared ratio of radial distances. Application to two recovery datasets from the literature confirm the utility of these type curves in normalized time for composite estimation of T and S.

Groundwater↗

USGS2003

Case study for delineating a contributing area to a well in a fractured siliciclastic-bedrock aquifer near Lansdale, Pennsylvania

A supply well used by the North Penn Water Authority near Lansdale, Pa., was selected as a case study for delineating a contributing area in a fractured siliciclastic-bedrock aquifer. The study emphasized the importance of refining the understanding of factors that control ground-water movement to the well by conducting (1) geophysical logging and flow measurements, (2) ground-water level monitoring, (3) aquifer testing, and (4) geochemical sampling. This approach could be applicable for other wells in siliciclastic-bedrock terranes, especially those of Triassic age in southeastern Pennsylvania. The principal methods for refining the understanding of hydrology at supply well MG-1125 were aquifer testing, water-level measurements, and geophysical logging. Results of two constant-discharge aquifer tests helped estimate the transmissivity of water-producing units and evaluate the anisotropy caused by dipping beds. Results from slug tests provided estimates of transmissivity that were used to evaluate the results from the constant-discharge aquifer tests. Slug tests also showed the wide distribution of transmissivity, indicating that ground-water velocities must vary considerably in the well field. Water-level monitoring in observation wells allowed maps of the potentiometric surface near the well field to be drawn. The measurements also showed that the hydraulic gradient can change abruptly in response to pumping from nearby supply wells. Water levels measured at a broader regional scale in an earlier study also provided a useful view of the potentiometric surface for purposes of delineating the contributing area. Geophysical logging and measurements of flow within wells showed that about 60 percent of water from supply well MG-1125 probably is contributed from relatively shallow water-producing fractures from 60 to 125 feet below land surface, but measurable amounts of water are contributed by fractures to a depth of 311 feet below land surface. Chemical samples supported the evidence that shallow fractures probably contribute significant amounts of water to well MG-1125. The large contribution of water from shallow fractures indicates that the area providing part of the recharge to the well is not far removed from the wellhead. Preliminary delineations of the contributing area and the 100-day time-of travel area were computed from a water budget and time-of-travel equation. These delineations provided insight into the size (but not the shape) of the contributing areas. Three other approaches were used and results compared: (1) uniform-flow equation, (2) hydrogeologic mapping, and (3) numerical modeling. The uniform-flow equation predicted a contributing area that seemed unrealistic—extending far across the ground-water divide into an adjacent watershed. Hydrogeologic mapping, if used with the potentiometric surface and constrained by the water budget, produced contributing area that was similar to that from numerical modeling. Numerical modeling allowed the incorporation of anisotropy caused by dipping water-producing units, differing transmissivity values of geologic units, and ground-water withdrawals from nearby supply wells. The numerical modeling showed that groundwater withdrawals from nearby supply wells affected the contributing area to supply well MG-1125 but had less effect on the 100-day time-of-travel area.

Water-Resources Investigations Report↗

USGS2002

Interpretation of borehole geophysical logs, aquifer-isolation tests, and water quality, supply wells 1 and 2, Willow Grove Naval Air Station/Joint Reserve Base, Horsham Township, Montgomery County, Pennsylvania

Ground water pumped from supply wells 1 and 2 on the Willow Grove Naval Air Station/Joint Reserve Base (NAS/JRB) provides water for use at the base, including potable water for drinking. The supply wells have been contaminated by volatile organic compounds (VOC's), particularly trichloroethylene (TCE) and tetrachloroethylene (PCE), and the water is treated to remove the VOC's. The Willow Grove NAS/JRB and surrounding area are underlain by sedimentary rocks of the Triassic-age Stockton Formation, which form a complex, heterogeneous aquifer. The ground-water-flow system for the supply wells was characterized by use of borehole geophysical logs and heatpulse-flowmeter measurements. The heatpulse-flowmeter measurements showed upward and downward borehole flow under nonpumping conditions in both wells. The hydraulic and chemical properties of discrete water-bearing fractures in the supply wells were characterized by isolating each water-bearing fracture with straddle packers. Eight fractures in supply well 1 and five fractures in supply well 2 were selected for testing on the basis of the borehole geophysical logs and borehole television surveys. Water samples were collected from each isolated fracture and analyzed for VOC?s and inorganic constituents. Fractures at 50–59, 79–80, 196, 124–152, 182, 241, 256, and 350–354 ft btoc (feet below top of casing) were isolated in supply well 1. Specific capacities ranged from 0.26 to 5.7 (gal/min)/ft (gallons per minute per foot) of drawdown. The highest specific capacity was for the fracture isolated at 179.8–188 ft btoc. Specific capacity and depth of fracture were not related in either supply well. The highest concentrations of PCE were in water samples collected from fractures isolated at 236.8–245 and 249.8–258 ft btoc, which are hydraulically connected. The concentration of PCE generally increased with depth to a maximum of 39 mg/L (micrograms per liter) at a depth of 249.8? 258 ft btoc and then decreased to 21 mg/L at a depth of 345.3–389 ft btoc. Fractures at 68–74, 115, 162, 182, 205, and 314 ft btoc were isolated in supply well 2. Specific capacities ranged from 0.08 to less than 2.9 (gal/ min)/ft. The highest specific capacity was for the fracture isolated at 157–165.2 ft btoc. Concentrations of detected VOC's in water samples were 3.6 mg/L or less. Lithologic units penetrated by both supply wells were determined by correlating naturalgamma and single-point-resistance borehole geophysical logs. All lithologic units are not continuous water-bearing units because water-bearing fractures are not necessarily present in the same lithologic units in each well. Although the wells penetrate the same lithologic units, the lithologic location of only three water-bearing fractures are common to both wells. The same lithologic unit may have different hydraulic properties in each well. A regional ground-water divide is southeast of the supply wells. From this divide, ground water flows northwest toward Park Creek, a tributary to Little Neshaminy Creek. Potentiometric-surface maps were prepared from water levels measured in shallow and deep wells. For both depth intervals, the direction of ground-water flow is toward the northwest. For most well clusters, the vertical head gradient is downward from the shallow to the deeper part of the aquifer. Pumping of the supply wells at times can cause the vertical flow direction to reverse.

Water-Resources Investigations Report↗

USGS2000

Simulation of aquifer tests and ground-water flowpaths at the local scale in fractured shales and sandstones of the Brunswick Group and Lockatong Formation, Lansdale, Montgomery County, Pennsylvania

The U.S. Geological Survey, as part of technical assistance to the U.S. Environmental Protection Agency, has constructed and calibrated models of local-scale ground-water flow in and near Lansdale, Pa., where numerous sources of industrial contamination have been consolidated into the North Penn Area 6 Superfund Site. The local-scale models incorporate hydrogeologic structure of northwest-dipping beds with uniform hydraulic properties identified in previous studies. Computations associated with mapping the dipping-bed structure into the three-dimensional model grid are handled by a preprocessor using a programmed geographic information system (GIS). Hydraulic properties are identified by calibration of the models using measured water levels during pumping and recovery from aquifer tests at three sites. Reduced flow across low-permeability beds is explicitly simulated. The dipping high-permeability beds are extensive in the strike direction but are of limited extent in the dip direction. This model structure yields ground-water-flow patterns characteristic of anisotropic aquifers; preferred flow is in the strike direction. The transmissivities of high-permeability beds in the local-scale models range from 142 to 1,900 ft 2 /d (feet squared per day) (13 to 177 m 2 /d). The hydraulic conductivities of low-permeability parts of the aquifer range from 9.6 x 10 -4 to 0.26 ft/d (feet per day) (2.9 x 10 -4 to 0.079 m/d). Storage coefficients and specific storage are very low, indicating the confined nature of the aquifer system. The calibrated models are used to simulate contributing areas of wells under alternative, hypothetical ground-water-management practices. Predictive contributing areas indicate the general characteristics of ground-water flow towards wells in the Lansdale area. Recharge to wells in Lansdale generally comes from infiltration near the well and over an area that extends upgradient from the well. The contributing areas for two wells pumping at 10 gal/min (gallons per minute) extend about 1,500 ft (feet) upgradient from the wells. The contributing area is more complex at ground-water divides and can extend in more than one direction to capture recharge from more than 3,300 ft away, for pumping at a rate of 30 gal/min. Locally, all recharge in the area of the pumping well is not captured; recharge in the downgradient direction about 150 ft from the pumping well will flow to other discharge locations.

Pennsylvania↗

USGS2000

Hydrogeology and simulation of ground-water flow at the Gettysburg Elevator Plant Superfund Site, Adams County, Pennsylvania

Ground water in Triassic-age sedimentary fractured-rock aquifers in the area of Gettysburg, Pa., is used as drinking water and for industrial and commercial supply. In 1983, ground water at the Gettysburg Elevator Plant was found by the Pennsylvania Department of Environmental Resources to be contaminated with trichloroethene, 1,1,1-trichloroethane, and other synthetic organic compounds. As part of the U.S. Environmental Protection Agency?s Comprehensive Environmental Response, Compensation, and Liability Act, 1980 process, a Remedial Investigation was completed in July 1991, a method of site remediation was issued in the Record of Decision dated June 1992, and a Final Design Report was completed in May 1997. In cooperation with the U.S. Environmental Protection Agency in the hydrogeologic assessment of the site remediation, the U.S. Geological Survey began a study in 1997 to determine the effects of the onsite and offsite extraction wells on ground-water flow and contaminant migration from the Gettysburg Elevator Plant. This determination is based on hydrologic and geophysical data collected from 1991 to 1998 and on results of numerical model simulations of the local ground-water flow-system. The Gettysburg Elevator Site is underlain by red, green, gray, and black shales of the Heidlersburg Member of the Gettysburg Formation. Correlation of natural-gamma logs indicates the sedimentary rock strike about N. 23 degrees E. and dip about 23 degrees NW. Depth to bedrock onsite commonly is about 6 feet but offsite may be as deep as 40 feet. The ground-water system consists of two zones?a thin, shallow zone composed of soil, clay, and highly weathered bedrock and a thicker, nonweathered or fractured bedrock zone. The shallow zone overlies the bedrock zone and truncates the dipping beds parallel to land surface. Diabase dikes are barriers to ground-water flow in the bedrock zone. The ground-water system is generally confined or semi-confined, even at shallow depths. Depth to water can range from flowing at land surface to more than 71 feet below land surface. Potentiometric maps based on measured water levels at the Gettysburg Elevator Plant indicate ground water flows from west to east, towards Rock Creek. Multiple-well aquifer tests indicate the system is heterogeneous and flow is primarily in dipping beds that contain discrete secondary openings separated by less permeable beds. Water levels in wells open to the pumped bed, as projected along the dipping stratigraphy, are drawn down more than water levels in wells not open to the pumped bed. Ground-water flow was simulated for steady-state conditions prior to pumping and long-term average pumping conditions. The three-dimensional numerical flow model (MODFLOW) was calibrated by use of a parameter estimation program (MODFLOWP). Steady-state conditions were assumed for the calibration period of 1996. An effective areal recharge rate of 7 inches was used in model calibration. The calibrated flow model was used to evaluate the effectiveness of the current onsite and offsite extraction well system. The simulation results generally indicate that the extraction system effectively captures much of the ground-water recharge at the Gettysburg Elevator Plant and, hence, contaminated ground-water migrating from the site. Some of the extraction wells pump at low rates and have very small contributing areas. Results indicate some areal recharge onsite will move to offsite extraction wells.

Pennsylvania↗

USGS1999

Simulating contaminant attenuation, double-porosity exchange, and water age in aquifers using MOC3D

MOC3D is a general-purpose computer model developed by the U.S. Geological Survey (USGS) for simulation of three-dimensional solute transport in ground water (Konikow and others, 1996). The model is an update to the widely used USGS two-dimensional solute-transport model (MOC) and is implemented as an optional “package” for the ground-water flow model MODFLOW (Harbaugh and McDonald, 1996). Directly coupling the time-tested MOC transport algorithms with the widely used MODFLOW program makes MOC3D a powerful tool for simulation of solute transport in ground water in many hydrogeologic settings. The model simulates transport processes that include: Advection - Transport of dissolved solutes at the same rate as the average ground-water flow velocity. Diffusion - Spreading of solute from areas of high concentration to areas of low concentration, caused by “random” molecular motion Dispersion - Diffusion-like spreading of solute that is caused primarily by spatial variability in aquifer properties, which results in spatial variability in transport velocity. Retardation - Reduction in the apparent solute velocity, compared to the ground-water velocity, caused by linear equilibrium sorption on aquifer materials. Decay - Disappearance of solute caused by reactions such as radioactive decay or biodegradation that are proportional to concentration. Growth - Creation (or disappearance) of solute mass caused by reactions that proceed independent of the solute concentration, such as some cases of biodegradation Double-porosity exchange - rate-limited exchange of solute mass between mobile and immobile zones; for example, between fractures and the rock matrix.

Fact Sheet↗

USGS1999

Age, double porosity, and simple reaction modifications for the MOC3D ground-water transport model

This report documents modifications for the MOC3D ground-water transport model to simulate (a) ground-water age transport; (b) double-porosity exchange; and (c) simple but flexible retardation, decay, and zero-order growth reactions. These modifications are incorporated in MOC3D version 3.0. MOC3D simulates the transport of a single solute using the method-of-characteristics numerical procedure. The age of ground water, that is the time since recharge to the saturated zone, can be simulated using the transport model with an additional source term of unit strength, corresponding to the rate of aging. The output concentrations of the model are in this case the ages at all locations in the model. Double porosity generally refers to a separate immobile-water phase within the aquifer that does not contribute to ground-water flow but can affect solute transport through diffusive exchange. The solute mass exchange rate between the flowing water in the aquifer and the immobile-water phase is the product of the concentration difference between the two phases and a linear exchange coefficient. Conceptually, double porosity can approximate the effects of dead-end pores in a granular porous media, or matrix diffusion in a fractured-rock aquifer. Options are provided for decay and zero-order growth reactions within the immobile-water phase. The simple reaction terms here extend the original model, which included decay and retardation. With these extensions, (a) the retardation factor can vary spatially within each model layer, (b) the decay rate coefficient can vary spatially within each model layer and can be different for the dissolved and sorbed phases, and (c) a zero-order growth reaction is added that can vary spatially and can be different in the dissolved and sorbed phases. The decay and growth reaction terms also can change in time to account for changing geochemical conditions during transport. The report includes a description of the theoretical basis of the model, a detailed description of input requirements and output options, and the results of model testing and evaluation. The model tests illustrate use of these modifications and demonstrate that accurate solutions can be obtained for these simple cases. Two test cases have no dispersion, illustrating the suitability of this method-of-characteristics model for simulation of advection-dominated transport in ground water.

Water-Resources Investigations Report↗