SEARCH · Search USGS
Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Edwards Air Force Base (EAFB) in southern California (fig. 1) has relied on ground water to meet its water-supply needs. The extraction of ground water has led to two major problems that can directly affect the mission of EAFB: declining water levels (more than 120 ft since the 1920s) and land subsidence, a gradual downward movement of the land surface (more than 4 ft since the late 1920s). As water levels decline, this valuable resource becomes depleted, thus requiring mitigating measures. Land subsidence has caused cracked (fissured) runways and accelerated erosion on Rogers lakebed. In 1988, the U.S. Geological Survey (USGS), in cooperation with the U.S. Air Force, began investigations of the effects of declining water levels and land subsidence at EAFB and possible mitigation measures, such as the injection of imported surface water into the ground-water system. The cooperative investigations included data collection and analyses, numerical simulations of ground-water flow and land subsidence, and development of a preliminary simulation-optimization model. The results of these investigations indicate that the injection of imported water may help to control land subsidence; however, the potential ground-water-quality impacts are unknown.
A digital model was developed to simulate steady-state ground-water flow in a 42,155-square-mile area of coastal Georgia and adjacent parts of South Carolina and Florida. The model was developed to (1) understand and refine the conceptual model of regional ground-water flow, (2) serve as a framework for the development of digital subregional ground-water flow and solute-transport models, and (3) serve as a tool for future evaluations of hypothetical pumping scenarios used to facilitate water management in the coastal area. Single-density ground-water flow was simulated using the U.S. Geological Survey finite-difference code MODFLOW-2000 for mean-annual conditions during predevelopment (pre?1900) and the years 1980 and 2000. The model comprises seven layers: the surficial aquifer system, the Brunswick aquifer system, the Upper Floridan aquifer, the Lower Floridan aquifer, and the intervening confining units. A combination of boundary conditions was applied, including a general-head boundary condition on the top active cells of the model and a time-variable fixed-head boundary condition along part of the southern lateral boundary. Simulated heads for 1980 and 2000 conditions indicate a good match to observed values, based on a plus-or-minus 10-foot (ft) calibration target and calibration statistics. The root-mean square of residual water levels for the Upper Floridan aquifer was 13.0 ft for the 1980 calibration and 9.94 ft for the 2000 calibration. Some spatial patterns of residuals were indicated for the 1980 and 2000 simulations, and are likely a result of model-grid cell size and insufficiently detailed hydraulic-property and pumpage data in some areas. Simulated potentiometric surfaces for predevelopment, 1980, and 2000 conditions all show major flow system features that are indicated by estimated peotentiometric maps. During 1980?2000, simulated water levels at the centers of pumping at Savannah and Brunswick rose more than 20 ft and 8 ft, respectively, in response to decreased pumping. Simulated drawdown exceeded 10 ft in the Upper Floridan aquifer across much of the western half of the model area, with drawdown exceeding 20 ft along parts of the western, northern, and southern boundaries where irrigation pumping increased during this period. From predevelopment to 2000 conditions, the simulated water budget showed an increase in inflow from, and decrease in outflow to, the general-head boundaries, and a reversal from net seaward flow to net landward flow across the coastline. Simulated changes in recharge and discharge distribution from predevelopment to 2000 conditions showed an increase in extent and magnitude of net recharge cells in the northern part of the model area, and a decrease in discharge or change to recharge in cells containing major streams and beneath major pumping centers. The model is relatively sensitive to pumping and the controlling head at the fixed-head boundary and less sensitive to the distribution of aquifer properties in general. Model limitations include: (1) its spatial scale and discretization, (2) the extent to which data are available to physically define the flow system, (3) the type of boundary conditions and controlling parameters used, (4) uncertainty in the distribution of pumping, and (5) uncertainty in field-scale hydraulic properties. The model could be improved with more accurate estimates of ground-water pumpage and better characterization of recharge and discharge.
Analysis of total dissolved nitrogen concentrations from soil water samples collected within the soil zone under septic tank leach fields in Spanish Springs Valley, Nevada, shows a median concentration of approximately 44 milligrams per liter (mg/L) from more than 300 measurements taken from four septic tank systems. Using two simple mass balance calculations, the concentration of total dissolved nitrogen potentially reaching the ground-water table ranges from 25 to 29 mg/L. This indicates that approximately 29 to 32 metric tons of nitrogen enters the aquifer every year from natural recharge and from the 2,070 houses that use septic tanks in the densely populated portion of Spanish Springs Valley. Natural recharge contributes only 0.25 metric tons because the total dissolved nitrogen concentration of natural recharge was estimated to be low (0.8 mg/L). Although there are many uncertainties in this estimate, the sensitivity of these uncertainties to the calculated load is relatively small, indicating that these values likely are accurate to within an order of magnitude. The nitrogen load calculation will be used as an input function for a ground-water flow and transport model that will be used to test management options for controlling nitrogen contamination in the basin.
In the early 1990's, two numerical models of the Death Valley regional ground-water flow system were developed by the U.S. Department of Energy. In general, the two models were based on the same basic hydrogeologic data set. In 1998, the U.S. Department of Energy requested that the U.S. Geological Survey develop and maintain a ground-water flow model of the Death Valley region in support of U.S. Department of Energy programs at the Nevada Test Site. The purpose of developing this 'second-generation' regional model was to enhance the knowledge an understanding of the ground-water flow system as new information and tools are developed. The U.S. Geological Survey also was encouraged by the U.S. Department of Energy to cooperate to the fullest extent with other Federal, State, and local entities in the region to take advantage of the benefits of their knowledge and expertise. The short-term objective of the Death Valley regional ground-water flow system project was to develop a steady-state representation of the predevelopment conditions of the ground-water flow system utilizing the two geologic interpretations used to develop the previous numerical models. The long-term objective of this project was to construct and calibrate a transient model that simulates the ground-water conditions of the study area over the historical record that utilizes a newly interpreted hydrogeologic conceptual model. This report describes the result of the predevelopment steady-state model construction and calibration. The Death Valley regional ground-water flow system is situated within the southern Great Basin, a subprovince of the Basin and Range physiographic province, bounded by latitudes 35 degrees north and 38 degrees 15 minutes north and by longitudes 115 and 118 degrees west. Hydrology in the region is a result of both the arid climatic conditions and the complex geology. Ground-water flow generally can be described as dominated by interbasinal flow and may be conceptualized as having two main components: a series of relatively shallow and localized flow paths that are superimposed on deeper regional flow paths. A significant component of the regional ground-water flow is through a thick Paleozoic carbonate rock sequence. Throughout the flow system, ground water flows through zones of high transmissivity that have resulted from regional faulting and fracturing. The conceptual model of the Death Valley regional ground-water flow system used for this study is adapted from the two previous ground-water modeling studies. The three-dimensional digital hydrogeologic framework model developed for the region also contains elements of both of the hydrogeologic framework models used in the previous investigations. As dictated by project scope, very little reinterpretation and refinement were made where these two framework models disagree; therefore, limitations in the hydrogeologic representation of the flow system exist. Despite limitations, the framework model provides the best representation to date of the hydrogeologic units and structures that control regional ground-water flow and serves as an important information source used to construct and calibrate the predevelopment, steady-state flow model. In addition to the hydrogeologic framework, a complex array of mechanisms accounts for flow into, through, and out of the regional ground-water flow system. Natural discharges from the regional ground-water flow system occur by evapotranspiration, springs, and subsurface outflow. In this study, evapotranspiration rates were adapted from a related investigation that developed maps of evapotranspiration areas and computed rates from micrometeorological data collected within the local area over a multiyear period. In some cases, historical spring flow records were used to derive ground-water discharge rates for isolated regional springs. For this investigation, a process-based, numerical model was developed to estimat
Water from wells and springs accounts for more than 90 percent of the public water supply in Calhoun County, Alabama. Springs associated with the Jacksonville Thrust Fault Complex are used for public water supply for the cities of Anniston and Jacksonville. The largest ground-water supply is Coldwater Spring, the primary source of water for Anniston, Alabama. The average discharge of Coldwater Spring is about 32 million gallons per day, and the variability of discharge is about 75 percent. Water-quality samples were collected from 6 springs and 15 wells in Calhoun County from November 2001 to January 2003. The pH of the ground water typically was greater than 6.0, and specific conductance was less than 300 microsiemens per centimeter. The water chemistry was dominated by calcium, carbonate, and bicarbonate ions. The hydrogen and oxygen isotopic composition of the water samples indicates the occurrence of a low-temperature, water-rock weathering reaction known as silicate hydrolysis. The residence time of the ground water, or ground-water age, was estimated by using analysis of chlorofluorocarbon, sulfur hexafluoride, and regression modeling. Estimated ground-water ages ranged from less than 10 to approximately 40 years, with a median age of about 18 years. The Spearman rho test was used to identify statistically significant covariance among selected physical properties and constituents in the ground water. The alkalinity, specific conductance, and dissolved solids increased as age increased; these correlations reflect common changes in ground-water quality that occur with increasing residence time and support the accuracy of the age estimates. The concentration of sodium and chloride increased as age increased; the correlation of these constituents is interpreted to indicate natural sources for chloride and sodium. The concentration of silica increased as the concentration of potassium increased; this correlation, in addition to the isotopic data, is evidence that silicate hydrolysis of clay minerals occurred. The geochemical modeling program NETPATH was used to investigate possible mixing scenarios that could yield the chemical composition of water collected from springs associated with the Jacksonville Thrust Fault Complex. The results of NETPATH modeling suggest that the primary source of water in Coldwater Spring is a deep aquifer, and only small amounts of rainwater from nearby sources are discharged from the spring. Starting with Piedmont Sports Spring and moving southwest along a conceptual ground-water flow path that parallels the Jacksonville Thrust Fault Complex, NETPATH simulated the observed water quality of each spring, in succession, by mixing rainwater and water from the spring just to the northeast of the spring being modeled. The percentage of rainwater and ground water needed to simulate the quality of water flowing from the springs ranged from 1 to 25 percent rainwater and 75 to 99 percent ground water.
An integrated borehole-geophysical and hydraulic investigation was conducted at the former landfill area near the University of Connecticut in Storrs, Connecticut, where solvents and landfill leachate have contaminated a fractured-bedrock aquifer. Borehole-geophysical techniques and hydraulic methods were used to characterize the site bedrock lithology and structure, fractures, and hydraulic properties. The geophysical and hydraulic methods included conventional logs, borehole imaging, borehole radar, flowmeter under ambient- and stressed hydraulic conditions, and discrete-zone hydraulic testing, sampling, and monitoring. The conventional geophysical-logging methods included caliper, deviation, electromagnetic induction, gamma, specific conductance, and fluid temperature. The advanced methods included optical and acoustic imaging of the borehole wall, heat-pulse flowmeter, and directional radar reflection. Borehole-geophysical methods were used to further define conductive features identified with surface-geophysical methods in the first phase of the investigation. The results of the surface- and borehole-geophysical logging were evaluated in an iterative and integrated manner to develop a conceptual model of ground-water flow at the site. The rock type, foliation, and fractures at the site were characterized from high-resolution optical televiewer (OTV) images of rocks penetrated by the boreholes and were compared to drilling logs and conventional geophysical logs. The rocks are interpreted as fine- to mediumgrained quartz-feldspar-biotite-garnet gneiss and schist with local intrusions of quartz diorite and pegmatite and minor concentrations of sulfide mineralization similar to rocks described as the Bigelow Brook Formation on regional geologic maps. Layers containing high concentrations of sulfide minerals appear as high electrical conductivity zones on electromagnetic-induction and borehole-radar logs. Foliation in the rocks generally strikes to the southwest and northeast, and dips to the northwest and southeast consistent with previous investigations in this area. The orientation of foliation, however, varies locally and with depth in some of the boreholes. These results are consistent with geologic mapping that has identified small-scale folding. The orientations of the transmissive fractures identified in the six boreholes logged for this investigation are similar to the fracture orientations mapped in a previous investigation. Many of these fractures are oriented with a north-northwest strike and have a shallow dip to the west. Other transmissive fractures have a southwest strike and dip at shallow angles to the northwest, and some strike roughly east-west and dip to the north and south. Flowmeter logging was used to identify transmissive fractures and to estimate the hydraulic properties in the boreholes. Ambient down flow was measured in one borehole, and ambient up flow and down flow were measured in another borehole. The other four bedrock boreholes did not have measurable vertical flow. Under low-rate pumping conditions (0.25 to 0.5 gallons per minute), one to three inflow zones were identified in each well. Commonly, fractures that are active under ambient conditions contribute to the well under pumping conditions. The ambient conditions were incorporated into the determination of the relative proportions of transmissivity. Specific capacity and transmissivity were determined for these open-hole low-rate pumping tests. Quasi-steady-state water levels were reached in four of the boreholes, including MW201R, MW204R, MW302R, and W202-NE. When pumped at low-rate conditions for 0.5 to 4 hours, the specific capacity ranged from 0.03 to 0.18 gallons per minute per foot. The open-hole transmissivity estimates ranged from 4.9 to 30 feet squared per day (ft2/d). Open-hole transmissivity was determined for boreholes that did not reach quasi-steady-state conditions under low-rate pumping conditions. Transmissivity was estimated for MW201R, MW202R, and MW203R using non-equilibrium methods, pumping rate, and the transient drawdown data to estimate the open-hole transmissivity. Transmissivity in these boreholes ranged from 0.98 to 3.2 ft 2 /d. The transmissivity and head of individual fractures or zones of fractures were estimated from heat-pulse flowmeter data acquired under ambient and stressed conditions. In the absence of ambient flow, data from two profiles of heat-pulse flowmeter data under two different stressed conditions were used to estimate the transmissivity and head of individual fracture zones. Only two boreholes, MW302R and W202-NE, had sufficient data for these analyses. The estimated transmissivity of individual transmissive zones ranged from 1.2 to 9.2 ft 2 /d. The transmissivity values determined by this numerical simulation method were less than the open-hole estimations, which were 15 and 30 ft 2 /d. Transmissivity also was measured directly over discrete intervals of the borehole using a straddle-packer apparatus and constant-rate pumping tests. Pumping rates were less than or equal to 0.25 gallons per minute. These discretezone single-hole pumping tests were conducted over a short period of time, usually about 30 minutes to 1 hour in duration. Pumping continued until the test zone reached a steady-state water level or until it was determined that the zone could not yield water at the pumped rate. The estimated transmissivity of individual transmissive zones ranged from about 0.21 to 11 ft 2 /d. The zone at a depth of 197 feet in W202-NE was the only zone that had discrete-interval testing with a straddle packer and sufficient heat-pulse flowmeter data for modeling the flow and estimating transmissivity and head. The two methods produced similar results. The straddle-packer method estimated a transmissivity of 4.7 ft 2 /d, and the heat-pulse flowmeter modeling results estimated a transmissivity of 6.9 ft 2 /d. A comparison of the transmissivity estimates indicate estimates typically are within an order of magnitude. The heat-pulse flowmeter methods used in this investigation to determine transmissivity of the boreholes and the individual fractures measure only the upper two or three orders of magnitude of transmissivity. Hence, other fractures in these boreholes permit the movement of water; their transmissivities, however, are lower than the detection limits of the methods that were used for this investigation and very small compared to the transmissive fractures that were studied. The data collected in this investigation were used to design discrete-zone monitoring systems for four of the boreholes used for monitoring. The results of the investigation are useful for refining the conceptual site model of ground-water flow, and for providing critical information for interpreting the results of water-quality sampling.
A computer simulation of Long Island 's regional groundwater system has been used to evaluate the effects that new-installed sewers will have on ground-water levels. Results indicate maximum water-table decliners of as much as up to 18 feet in central Nassau County and about 9 feet in Suffolk County. Total stream base flows and freshwater outflow to the south shore bay system will decrease by 22%. The regional scale of the model does not permit detailed predictions for individual streams. To quantify the effects of lowered ground-water levels on individual streams, two fine-scale sub-regional models have been designed. This report, the first in a three-part series describing the simulated effects of sewers in southern Nassau and southwestern Suffolk Counties, presents the hydrogeologic setting, pertinent literature, modeling strategy, subregional model design, and the results obtained to date from the regional ground-water model. The regional model results described will be used in the later reports to generate flux boundary conditions for the subregional models.
Water supply for the Borough of Martinsburg, Pa., is from two well fields (Wineland and Hershberger) completed in carbonate-bedrock aquifers in the Morrison Cove Valley. Water supply is plentiful; however, waters with high concentrations of nitrate are a concern. This report describes the sources of water and contaminants to the supply wells. A review of previous investigations was used to establish the aquifer framework and estimate aquifer hydraulic properties. Aquifer framework and simulation of ground-water flow in a 25-square-mile area using the MODFLOW model helped to further constrain aquifer hydraulic properties and identify water-source areas in the zone of contribution of ground water to the well fields. Flow simulation identified potential contaminant-source areas. Data on contaminants and geochemical characteristics of ground water at the well fields were compared to the results of flow simulation. The Woodbury Anticline controls the aquifer framework near the well fields and four carbonate-bedrock formations contain the primary aquifers. Three carbonate-bedrock aquifers of Ordovician age overlie the Gatesburg aquifer of Cambrian age on the flanks of the anticline. Fracture, not conduit, permeability was determined to be the dominant water-bearing characteristic of the bedrock. The horizontal hydraulic conductivity of the Gatesburg aquifer is about 36 feet per day. The other carbonate aquifers (Nittany/Stonehenge, Bellefonte/Axemann, and Coburn through Loysburg aquifers) overlying and flanking the Gatesburg aquifer have horizontal hydraulic conductivities of about 1 foot per day. Regional directions of ground-water flow are toward the major streams with Clover Creek as the major discharge point for ground water in the east. Ground-water flow to the well fields is anisotropic with a 5:1 preferential horizontal direction along strike of the axial fold of the anticline. Thus, the zone of contribution of ground water to the well fields is elongate in a north-south direction along the anticline axis, with the majority of the flow to the well fields originating from the south. Human activity in the areal extent of the zone of contribution to the well fields was the source of contaminants. The areal extent of the zone of contribution included both urban areas in the Borough and a large amount of agricultural land. By relating results of flow simulation, natural geochemistry, and analyses of anthropogenic (human-made) contaminants, the source areas for water and contaminants were determined with more confidence than by using only flow simulation. Analysis of natural geochemistry identified water sources from both limestone and dolomite aquifers. Geochemistry results also indicated fractures, not conduits, were the dominant source of water from aquifers; however, quantitative source identification was not possible. Chemical ratios of chloride and bromide were useful to show that all samples of ground water had sources with chemical contributions from land surface. Nitrogen isotope ratio analysis indicated animal manure as the possible primary source of nitrate in most ground water. Some of the nitrate in ground water had chemical fertilizer as a source. At the Wineland well field, chemical fertilizer was likely the source of nitrate. The nitrate in water from the Hershberger well field was from a mixture of fertilizer and animal-manure sources. Human sewage was ruled out as a major source of nitrate in water from the municipal wells by results showing 1) wastewater compounds in sewage were rarely detected and 2) a mass-balance calculation indicating the small contribution of nitrogen that could be attributed to septic systems.
A regional ground-water model of the Floridan aquifer was constructed for an 875-square-mile portion of the rapidly developing area north of Tampa Bay. The digital model was calibrated by comparing observed (March 1974 and May 1975) and computed potentiometric heads. A good comparison was obtained by adjusting leakance and transmissivity. Differences between the computed and measured potentiometric surface were generally less than 3 feet, with a maximum error of 15 feet. The calibrated model may be used as a predictive tool. For example, the model could be used to evaluate the regional effects of increased ground-water withdrawal on the Florida aquifer. (Woodard-USGS)
On Long Island, base flow under nonurbanized conditions constitutes 90 to 95% of total stream discharge. Base-flow data from 19 continuously gaged streams are presented as monthly mean and annual mean discharge for water years 1960-75, which includes the 1962-66 drought. The data were derived by hydrograph-separation procedures that isolate mean daily base flow from mean daily discharge. A close empirical relationship between annual mean base flow and stream discharge at the 55-% duration point facilities estimation of average base flow and can be used in place of the more time-consuming hydrograph-separation technique. These data are needed in calibration of computer models that will be used to predict the effects of hydrologic stresses, such as sewering, on the Long Island ground-water system. (USGS)
Recharge and discharge are hydrological processes that cause Everglades surface water to be exchanged for subsurface water in the peat soil and the underlying sand and limestone aquifer. These interactions are thought to be important to water budgets, water quality, and ecology in the Everglades. Nonetheless, relatively few studies of surface water and ground water interactions have been conducted in the Everglades, especially in its vast interior areas. This report is a product of a cooperative investigation conducted by the USGS and the South Florida Water Management District (SFWMD) aimed at developing and testing techniques that would provide reliable estimates of recharge and discharge in interior areas of WCA-2A (Water Conservation Area 2A) and several other sites in the central Everglades. The new techniques quantified flow from surface water to the subsurface (recharge) and the opposite (discharge) using (1) Darcy-flux calculations based on measured vertical gradients in hydraulic head and hydraulic conductivity of peat; (2) modeling transport through peat and decay of the naturally occurring isotopes 224 Ra and 223 Ra (with half-lives of 4 and 11 days, respectively); and (3) modeling transport and decay of naturally occurring and "bomb-pulse" tritium (half-life of 12.4 years) in ground water. Advantages and disadvantages of each method for quantifying recharge and discharge were compared. In addition, spatial and temporal variability of recharge and discharge were evaluated and controlling factors identified. A final goal was to develop appropriately simplified (that is, time averaged) expressions of the results that will be useful in addressing a broad range of hydrological and ecological problems in the Everglades. Results were compared with existing information about water budgets from the South Florida Water Management Model (SFWMM), a principal tool used by the South Florida Water Management District to plan many of the hydrological aspects of the Everglades restoration.
Ground-water models attempt to represent an actual ground-water system with a mathematical counterpart. The conceptualization of how and where water originates in the ground-water-flow system and how and where it leaves the system is critical to the development of an accurate model. The mathematical representation of these boundaries in the model is important because many hydrologic boundary conditions can be mathematically represented in more than one way. The determination of which mathematical representation of a boundary condition is best usually is dependent upon the objectives of the study. This report focuses on the specific aspect of describing different ways to simulate, in a numerical model, the physical features that act as hydrologic boundaries in an actual ground-water system. The ramifications, benefits, and limitations of each approach are enumerated, and descriptions of the representation of boundaries in models for Long Island, New York, and the Middle Rio Grande Basin, New Mexico, illustrate the application of some of the methods.
The Sparta aquifer is a confined aquifer of great regional importance that comprises a sequence of unconsolidated sand, silt, and clay units extending across much of eastern and southeastern Arkansas and into adjoining States. Water use from the aquifer has doubled since 1975 and continues to increase, and large water-level declines are occurring in many areas of the aquifer. To focus State attention and resources on the growing problem and to provide a mechanism for locally based education and management, the Arkansas Soil and Water Conservation Commission has designated Critical Ground-Water Areas in some counties (see page 6, ?What is a Critical Ground-Water Area??). Ground-water modeling study results show that the aquifer cannot continue to meet growing water-use demands. Dewatering of the primary producing sands is predicted to occur within 10 years in some areas if current trends continue. The predicted dewatering will cause reduced yields and damage the aquifer. Modeling also shows that a concerted ground-water conservation management plan could enable sustainable use of the aquifer. Water-conservation measures and use of alternative sources that water managers in Union County (an area of high demand and growth in Arkansas' initial five-county Critical Ground-Water Area) think to be realistic options result in considerable recovery in water levels in the aquifer during a 30-year model simulation.
The Potomac-Raritan-Magothy aquifer system is one of the primary sources of potable water in the Coastal Plain of New Jersey, particularly in heavily developed areas along the Delaware River. In Pennsauken Township, Camden County, local drinking-water supplies from this aquifer system have been contaminated by hexavalent chromium at concentrations that exceed the New Jersey maximum contaminant level. In particular, ground water at the Puchack well field has been adversely affected to the point where, since 1984, water is no longer withdrawn from this well field for public supply. The area that contains the Puchack well field was added to the National Priorities List in 1998 as a Superfund site. The U.S. Geological Survey (USGS) conducted a reconnaissance study from 1996 to 1998 during which hydrogeologic and water-quality data were collected and a ground-water-flow model was developed to describe the conditions in the aquifer system in the Pennsauken Township area. The current investigation by the USGS, in cooperation with the U.S. Environmental Protection Agency (USEPA), is an extension of the previous study. Results of the current study can be applied to a Remedial Investigation and Feasibility Study conducted at the Puchack well field Superfund site. The USGS study collected additional data on the hydrogeology and water-quality in the area. These data were incorporated into a refined model of the ground-water-flow system in the Potomac-Raritan-Magothy aquifer system. A finite-difference model was developed to simulate ground-water flow and the advective transport of chromium-contaminated ground water in the aquifers of the Potomac-Raritan-Magothy aquifer system in the Pennsauken Township area. An 11-layer model was used to represent the complex hydrogeologic framework. The model was calibrated using steady-state water-level data from March 1998, April 1998, and April 2001. Water-level recovery during the shutdown of Puchack 1 during March to April 1998 was simulated to evaluate model performance in relation to changing stresses. The Delaware River contributes appreciable-flow to the ground-water system from areas where the Middle and Lower aquifers crop out beneath the river. A transient simulation of an aquifer test near the Delaware River was run to help characterize the hydraulic conductivity of the riverbed sediments represented in the model. Vertical flow across confining units between the aquifers is highly variable and is important in the movement of water and associated contaminants through the flow system. The model was imbedded within a regional model of the Potomac-Raritan-Magothy aquifer system in Camden County. In general, a simulation of baseline conditions, which can provide a representation on which simulations of various alternatives can be based for the feasibility study, incorporated average conditions from 1998 to 2000. Ground-water withdrawals within the model area during this period averaged about 14 Mgal/d. Regional ground-water flow is from recharge areas and from the Delaware River to downgradient pumped wells located just east of the model area in central Camden County. Simulation results show an important connection between the Intermediate sand and the Lower aquifer of the Potomac-Raritan-Magothy aquifer system in the vicinity of the chromium-contaminated area. The Delaware River contributes nearly 10 Mgal/d to the flow system, whereas recharge contributes about 6 Mgal/d. Ground-water withdrawals within the model area account for nearly 14 Mgal/d (mostly from the Lower aquifer of the Potomac-Raritan-Magothy aquifer system).
A series of databases and an accompanying screening model were constructed by the U.S. Geological Survey, in cooperation with the National Park Service, to better understand the regional ground-water-flow system and its relation to stream drainage in the St. Croix River Basin. The St. Croix River and its tributaries drain about 8,000 square miles in northeastern Minnesota and northwestern Wisconsin. The databases contain information for the entire St. Croix River Basin pertaining to well logs, lithology, thickness of lithologic groups, ground-water levels, streamflow, and well pumpage. Maps and generalized cross sections created from the compiled data show the lithologic groups, extending from the water table to the crystalline bedrock, through which ground water flows. These lithologic groups are: fine-grained unconsolidated deposits; coarse-grained unconsolidated deposits; sandstone bedrock; carbonate bedrock; and other bedrock lithologies including shale, siltstone, conglomerate, and igneous intrusions. The steady-state screening model treats the ground-water-flow system as a single layer with transmissivity zones that reflect the distribution of lithologic groups, and with recharge zones that correspond to general areas of high or low evapotranspiration. The model includes representation of second- and higher-order streams and municipal and other high-capacity production wells. The analytic-element model code GFLOW was used to simulate the regional ground-water flow, the water-table surface across the St. Croix River Basin, and base-flow contributions from ground water to streams. In addition, the model routes tributary base flow through the stream network to the St. Croix River. The parameter-estimation inverse model UCODE was linked to the GFLOW model to select the combination of parameter values best able to match over 5,000 water-level measurements and base-flow estimates at 22 streamflow-gaging stations. Results from the calibrated screening model show ground-water contributing areas for selected stream reaches within the basin. The delineation of these areas is useful to water-resource managers concerned with protection of fisheries and other resources. The model results also identify the areas of the basin where ground-water travel time from the water table to streams and wells is relatively short (less than 50 years). Ninety percent of the simulated ground-water pathlines require travel times between 3 and 260 years. The median pathline distance traversed and the median pathline velocity were 1.7 mi and 177 ft/y, respectively. It is important to recognize the limitations of this screening model. Heterogeneities in subsurface properties and in recharge rates are considered only at a very broad scale (miles to tens of miles). No account is taken of vertical variations in properties or pumping rates, and no provision is made to account for stacked ground-water-flow systems that have different flow patterns at different depths. Small-scale (hundreds to thousands of feet) flow systems associated with minor water bodies are neglected, and as a result, the model is not useful for simulating typical site-specific problems. Despite its limitations, the model serves as a framework for understanding the regional pattern of ground-water flow and as a starting point for a generation of more targeted and detailed ground-water models that would be needed to address emerging water-supply and water-quality concerns in the St. Croix River Basin.
Increasing urbanization in the 67-square-mile Clover Creek Basin has generated interest in the effects of land-use changes on local water quality. To investigate these effects, water-quality and streamflow data were collected from 19 surface-water sites in the basin over a 16-month period from January 1991 through April 1992. These data were used to understand the effects of surficial geology, land-use practices, and wastewater disposal practices on surface-water quality within the basin. The basin was divided into four drainage subbasins with dissimilar hydrogeologic, land-use, and water-quality characteristics. In the Upper Clover Creek subbasin, the high permeability of surficial geologic materials promotes infiltration of precipitation to ground water and thus attenuates the response of streams to rainfall. Significant interaction occurs between surface and ground water in this subbasin, and nitrate concentrations and specific conductance values, similar to those found historically in local ground water, indicate that sources such as subsurface waste-disposal systems and fertilizers are affecting surface- water quality in this area. In the Spanaway subbasin, the presence of Spanaway and Tule Lakes affects water quality, primarily because of the reduced velocity and long residence time of water in the lakes. Reduced water velocity and long residence times (1) cause settling of suspended materials, thereby reducing concentrations of suspended sediment and constituents that are bound to the sediment; (2) promote biological activity, which tends to trap nutrients in the lakes; and (3) allow dispersion to attenuate peaks in discharge and water-quality constituent concentrations. In the North Fork subbasin, the low permeability of surficial geologic materials and areas of intensive land development inhibit infiltration of precipitation and thus promote surface runoff to streams. Surface pathways provide little attenuation of storm runoff and result in rapid increases in stream discharge in response to rainfall. Substantial increases in concentrations of constituents associated with surface wash off, for example, suspended sediment, ammonia, phosphorus, and fecal coliform, also were observed in this subbasin during rainfall. In the Lower Clover Creek subbasin, which is the most downstream subbasin, stream-discharge and water-quality characteristics show the integrated effects of the entire basin. The data show that further characterization of local ground water and discharge from stormwater outfalls entering Clover Creek and its tributaries would be necessary to successfully apply a numerical water-quality model to the basin.
Contaminants introduced into the subsurface of the Nevada Test Site at Rainier Mesa and Shoshone Mountain by underground nuclear testing are of concern to the U.S. Department of Energy and regulators responsible for protecting human health and safety. Although contaminants were introduced into low-permeability rocks above the regional flow system, the potential for contaminant movement away from the underground test areas and into the accessible environment is greatest by ground-water transport. The primary hydrologic control on this transport is evaluated and examined through a series of contour maps developed to represent the water-level distribution within each of the major aquifers underlying the area. Aquifers were identified and their extents delineated by merging and analyzing multiple hydrostratigraphic framework models developed by other investigators from existing geologic information. The contoured water-level distribution in each major aquifer was developed from a detailed evaluation and assessment of available water-level measurements. Multiple spreadsheets that accompany this report provide pertinent water-level and geologic data by well or drill hole. Aquifers are mapped, presented, and discussed in general terms as being one of three aquifer types?volcanic aquifer, upper carbonate aquifer, or lower carbonate aquifer. Each of these aquifer types was subdivided and mapped as independent continuous and isolated aquifers, based on the continuity of its component rock. Ground-water flow directions, as related to the transport of test-generated contaminants, were developed from water-level contours and are presented and discussed for each of the continuous aquifers. Contoured water-level altitudes vary across the study area and range from more than 5,000 feet in the volcanic aquifer beneath a recharge area in the northern part of the study area to less than 2,450 feet in the lower carbonate aquifer in the southern part of the study area. Variations in water-level altitudes within any single continuous aquifer range from a few hundred feet in a lower carbonate aquifer to just more than 1,100 feet in a volcanic aquifer. Flow directions throughout the study area are dominantly southward with minor eastward or westward deviations. Primary exceptions are westward flow in the northern part of the volcanic aquifer and eastward flow in the eastern part of the lower carbonate aquifer. Northward flow in the upper and lower carbonate aquifers in the northern part of the study area is possible but cannot be substantiated because data are lacking. Interflow between continuous aquifers is evaluated and mapped to define major flow paths. These flow paths delineate tributary flow systems, which converge to form the regional ground-water flow system. The implications of these tributary flow paths in controlling transport away from the underground test areas at Rainier Mesa and Shoshone Mountain are discussed. The obvious data gaps contributing to uncertainties in the delineation of aquifers and development of water-level contours are identified and evaluated.