Search USGS⌕ Search

SEARCH · Search USGS

Results for “Ground Water Modeling Investigation Report”

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.

At least 955 records · Page 53Linked to original sources
USGS2000

Investigation of the geology and hydrology of the Mogollon Highlands of central Arizona: a project of the Arizona Rural Watershed Initiative

The Mogollon Highlands of east central Arizona is a region of forested plateau and mountains, deep, sheerwalled canyons, and desert valleys. Known for its scenic beauty and characterized by a generally mild climate, the area, though still sparsely populated, attracts an increasing number of tourists and summer residents. Furthermore, the permanent population is expected to nearly double over the next 50 years. Consequently, there is increased pressure on the water resources of this area for several sometimes conflicting uses. Rational management of water resources is necessary to meet increased domestic requirements while ensuring an adequate supply of water for commercial and agricultural use, for Indian lands, and for preservation of valued environmental elements, including surface waters, riparian woodlands, forest and grassland areas, and wildlife and aquatic habitat. Such management requires an understanding of the relations among different components of the hydrologic system—recharge areas, surface flows, shallow aquifers, deep aquifers, discharge areas, and the regional ground-water flow system—and how each is affected by geology, climate, topography, and human use. The U.S. Geological Survey (USGS) is conducting an assessment of the hydrogeology of the Mogollon Highlands in cooperation with the Arizona Department of Water Resources. The study, funded through the State’s Rural Watershed Initiative program, is one of three assessments being conducted by the USGS. Assessments also are underway in the Upper-Middle Verde River watershed and on the Coconino Plateau. Each study has as its objectives: (1) the collection, compilation, and evaluation of all existing geologic, hydrologic, and related data pertaining to the study area and the creation of a data base that is readily accessible to the public and (2) developing an understanding of the hydrogeologic framework, which is the relation between geologic and hydrologic properties, that can be used for water-- resources management purposes and that will support the development of an interpretive and predictive model to estimate the effects of climate and water use on the sustainability of regional water resources. Although the three contiguous areas in north-central Arizona are being studied separately, a single data base is being constructed from which data on each area can be extracted separately.

Arizona↗
USGS2000

Computer-model analysis of ground-water flow and simulated effects of contaminant remediation at Naval Weapons Industrial Reserve Plant, Dallas, Texas

In June 1993, the Department of the Navy, Southern Division Naval Facilities Engineering Command (SOUTHDIV), began a Resource Conservation and Recovery Act (RCRA) Facility Investigation (RFI) of the Naval Weapons Industrial Reserve Plant (NWIRP) in north-central Texas. The RFI has found trichloroethene, dichloroethene, vinyl chloride, as well as chromium, lead, and other metallic residuum in the shallow alluvial aquifer underlying NWIRP. These findings and the possibility of on-site or off-site migration of contaminants prompted the need for a ground-water-flow model of the NWIRP area. The resulting U.S. Geological Survey (USGS) model: (1) defines aquifer properties, (2) computes water budgets, (3) delineates major flowpaths, and (4) simulates hydrologic effects of remediation activity. In addition to assisting with particle-tracking analyses, the calibrated model could support solute-transport modeling as well as help evaluate the effects of potential corrective action. The USGS model simulates steadystate and transient conditions of ground-water flow within a single model layer. The alluvial aquifer is within fluvial terrace deposits of Pleistocene age, which unconformably overlie the relatively impermeable Eagle Ford Shale of Late Cretaceous age. Over small distances and short periods, finer grained parts of the aquifer are separated hydraulically; however, most of the aquifer is connected circuitously through randomly distributed coarser grained sediments. The top of the underlying Eagle Ford Shale, a regional confining unit, is assumed to be the effective lower limit of ground-water circulation and chemical contamination. The calibrated steady-state model reproduces long-term average water levels within +5.1 or –3.5 feet of those observed; the standard error of the estimate is 1.07 feet with a mean residual of 0.02 foot. Hydraulic conductivity values range from 0.75 to 7.5 feet per day, and average about 4 feet per day. Specific yield values range from 0.005 to 0.15 and average about 0.08. Simulated infiltration rates range from 0 to 2.5 inches per year, depending mostly on local patterns of ground cover. Computer simulation indicates that, as of December 31, 1998, remediation systems at NWIRP were removing 7,375 cubic feet of water per day from the alluvial aquifer, with 3,050 cubic feet per day coming from aquifer storage. The resulting drawdown prevented 1,800 cubic feet per day of ground water from discharging into Cottonwood Bay, as well as inducing another 1,325 cubic feet per day into the aquifer from the bay. An additional 1,200 cubic feet of water per day (compared to pre-remediation conditions) was prevented from discharging into the west lagoon, east lagoon, Mountain Creek Lake, and Mountain Creek swale. Particle-tracking simulations, assuming an aquifer porosity of 0.15, were made to delineate flowpath patterns, or contaminant “capture zones,” resulting from 2.5- and 5-year periods of remediation activity at NWIRP. The resulting flowlines indicate three such zones, or areas from which ground water is simulated to have been removed during July 1996–December 1998, as well as extended areas from which ground water would be removed during the next 2.5 years (January 1999– June 2001). Simulation indicates that, as of December 31, 1998, the recovery trench was intercepting about 827 cubic feet per day of ground water that—without the trench—would have discharged into Cottonwood Bay. During this time, the trench is simulated to have removed about 3,221 cubic feet per day of water from the aquifer, with about 934 cubic feet per day (29 percent) coming from the south (Cottonwood Bay) side of the trench.

Texas↗
USGS2005

Questa baseline and premining ground-water quality investigation. 8. Lake-sediment geochemical record from 1960 to 2002, Eagle Rock and Fawn Lakes, Taos County, New Mexico

Geochemical studies of lake sediment from Eagle Rock Lake and upper Fawn Lake were conducted to evaluate the effect of mining at the Molycorp Questa porphyry molybdenum deposit located immediately north of the Red River. Two cores were taken, one from each lake near the outlet where the sediment was thinnest, and they were sampled at 1-cm intervals to provide geochemical data at less than 1-year resolution. Samples from the core intervals were digested and analyzed for 34 elements using ICP–AES (inductively coupled plasma–atomic emission spectrometry). The activity of 137 Cs has been used to establish the beginning of sedimentation in the two lakes. Correlation of the geochemistry of heavy-mineral suites in the cores from both Fawn and Eagle Rock Lakes has been used to develop a sedimentation model to date the intervals sampled. The core from upper Fawn Lake, located upstream of the deposit, provided an annual sedimentary record of the geochemical baseline for material being transported in the Red River, whereas the core from Eagle Rock Lake, located downstream of the deposit, provided an annual record of the effect of mining at the Questa mine on the sediment in the Red River. Abrupt changes in the concentrations of many lithophile and deposit-related metals occur in the middle of the Eagle Rock Lake core, which we correlate with the major flood-of-record recorded at the Questa gage at Eagle Rock Lake in 1979. Sediment from the Red River collected at low flow in 2002 is a poor match for the geochemical data from the sediment core in Eagle Rock Lake. The change in sediment geochemistry in Eagle Rock Lake in the post-1979 interval is dramatic and requires that a new source of sediment be identified that has substantially different geochemistry from that in the pre-1979 core interval. Loss of mill tailings from pipeline breaks are most likely responsible for some of the spikes in trace-element concentrations in the Eagle Rock Lake core. Enrichment of Al 2 O 3 , Cu, and Zn occurred as a result of chemical precipitation of these metals from ground water upstream in the Red River. Comparisons of the geochemistry of the post-1979 sediment core with both mine wastes and with premining sediment from the vicinity of the Questa mine indicate that both are possible sources for this new component of sediment. Existing data have not resolved this enigma.

New Mexico↗
USGS1993

Geology and ground-water resources in the Zebulon area, Georgia

The current (1991) surface-water source of drinking-water supply for the city of Zebulon, Pike County, Georgia, no longer provides an adequate water supply and periodically does not meet water-quality standards. The hydrogeology of crystalline rocks in the Zebulon area was evaluated to assess the potential of ground-water resources as a supplemental or alternative source of water to present surface-water supplies. As part of the ground-water resource evaluation, well location and construction data were compiled, a geologic map was constructed, and ground water was sampled and analyzed. Three mappable geologic units delineated during this study provide a basic understanding of hydrogeologic settings in the Zebulon area. Rock types include a variety of aluminosilicate schists, granitic rocks, amphibolites/honblende gneisses, and gondites. Several geologic features that may enhance ground-water availability were identified in the study area. These features include contacts between contrasting rock types, where a high degree of differential weathering has occurred, and well-developed structural features, such as foliation and jointing are present. High-yielding wells (greater than 25 gallons per minute) and low-yielding wells (less than one gallon per minute) were located in all three geologic units in a variety of topographic settings. Well yields range from less than one gallon per minute to 250 gallons per minute. The variable total depths and wide ranges of casing depths of the high-yielding wells are indicative of variations in depths to water-bearing zones and regolith thicknesses, respectively. The depth of water-bearing zones is highly variable, even on a local scale. Analyses of ground-water samples indicate that the distribution of iron concentration is as variable as well yield in the study area and does not seem to be related to a particular rock type. Iron concentrations in ground-water samples ranged from 0.02 to 5.3 milligrams per liter. Both iron concentration and well yield vary substantially over a relatively small area. Implementation and Verification of a One-Dimensional, Unsteady-Flow Model for Spring Brook near Warrenville, Illinois By Mary J. Turner, Anthony P. Pulokas, and Audrey L. Ishii Abstract A one-dimensional, unsteady-flow model, Full EQuations (FEQ) model, based on de Saint-Venant equations for dynamic flow in open channels, was calibrated and verified for a 0.75-mile reach of Spring Brook, a tributary to the West Branch Du Page River, near Warrenville in northeastern Illinois. The model was used to simulate streamflow in a small urban stream reach with two short culverts, one with overbank flow around the culvert during high flows. Streamflow data were collected on the reach during three high-flow periods. Data from one period were used to calibrate the model, and data from the other two periods were used to verify the model. Stages and discharges over the periods were simulated, and the results were compared graphically with stage and discharge data collected at 10 sites in the study reach. Errors in simulated stage and discharge were small except when debris, not represented in the model, clogged the culvert. The effects of changes in physical and computational model parameters also were studied. The model was insensit'lve to replacement of measured cross sections with interpolated cross sections, especially if the measured thalweg elevation was preserved. Variation of the roughness, slope, and length of the culvert over-bank section, as well as the chosen representative measured cross section, caused only slight changes in the simulated peak stage and discharge. Changes in the modeled culvert area caused large differences in the simulated highflows in the vicinity of the culvert, whereas simulated low flows were unaffected. At all flows, the misrepresentation of the culvert area caused the simulated water-surface elevations to deviate from the measured elevations, especially on the falling

Georgia↗
USGS2007

Estimated infiltration, percolation, and recharge rates at the Rillito Creek focused recharge investigation site, Pima County, Arizona

A large fraction of ground water stored in the alluvial aquifers in the Southwest is recharged by water that percolates through ephemeral stream-channel deposits. The amount of water currently recharging many of these aquifers is insufficient to meet current and future demands. Improving the understanding of streambed infiltration and the subsequent redistribution of water within the unsaturated zone is fundamental to quantifying and forming an accurate description of streambed recharge. In addition, improved estimates of recharge from ephemeral-stream channels will reduce uncertainties in water-budget components used in current ground-water models. This chapter presents a summary of findings related to a focused recharge investigation along Rillito Creek in Tucson, Arizona. A variety of approaches used to estimate infiltration, percolation, and recharge fluxes are presented that provide a wide range of temporal- and spatial-scale measurements of recharge beneath Rillito Creek. The approaches discussed include analyses of (1) cores and cuttings for hydraulic and textural properties, (2) environmental tracers from the water extracted from the cores and cuttings, (3) seepage measurements made during sustained streamflow, (4) heat as a tracer and numerical simulations of the movement of heat through the streambed sediments, (5) water-content variations, (6) water-level responses to streamflow in piezometers within the stream channel, and (7) gravity changes in response to recharge events. Hydraulic properties of the materials underlying Rillito Creek were used to estimate long-term potential recharge rates. Seepage measurements and analyses of temperature and water content were used to estimate infiltration rates, and environmental tracers were used to estimate percolation rates through the thick unsaturated zone. The presence or lack of tritium in the water was used to determine whether or not water in the unsaturated zone infiltrated within the past 40 years. Analysis of water-level and temporal-gravity data were used to estimate recharge volumes. Data presented in this chapter were collected from 1999 though 2002. Precipitation and streamflow during this period were less than the long-term average; however, two periods of significant streamflow resulted in recharge—one in the summer of 1999 and the other in the fall/winter of 2000. Flux estimates of infiltration and recharge vary from less than 0.1 to 1.0 cubic meter per second per kilometer of streamflow. Recharge-flux estimates are larger than infiltration estimates. Larger recharge fluxes than infiltration fluxes are explained by the scale of measurements. Methods used to estimate recharge rates incorporate the largest volumetric and temporal scales and are likely to have fluxes from other nearby sources, such as unmeasured tributaries, whereas the methods used to estimate infiltration incorporate the smallest scales, reflecting infiltration rates at individual measurement sites.

Arizona↗
USGS2006

Simulation of proposed increases in ground-water withdrawals on the Atlantic City 800-foot sand, New Jersey Coastal Plain

The confined Atlantic City 800-foot sand and the unconfined Kirkwood-Cohansey aquifer system (surficial aquifer) are major sources of water for southeastern New Jersey. Because of recent concerns about streamflow depletion resulting from ground-water withdrawals and the potential ecological effects on stream habitat in the area, the focus on future withdrawals has been shifted away from the surficial aquifer to the confined Atlantic City 800-foot sand until the effects of increased withdrawals from the surficial aquifer can be investigated. A study was conducted to evaluate the effects of seven proposed increases in ground-water withdrawals from the Atlantic City 800-foot sand and the Kirkwood-Cohansey aquifer system on the Atlantic City 800-foot sand. The proposed withdrawals are increases above the 2004 allocated rates (full allocation). The effects of full-allocation ground-water withdrawals and the cumulative effect of withdrawals for each of seven proposed increases in withdrawals were simulated using three previously published ground-water flow models: the New Jersey Coastal Plain Regional Aquifer System Analysis model, the Coastal Plain Optimization model, and a model of the Atlantic City 800-foot sand in Atlantic County, New Jersey. These models were used to simulate changes in water levels, the source supplying the increased ground-water flow, and the effects on saltwater movement towards production wells in Cape May County as a result of the proposed increased withdrawals at proposed or existing wells. The results of the simulations represent the effects of the proposed increase from full-allocation withdrawals to an additional 1,825 Mgal/yr (million gallons per year) from the Atlantic City 800-foot sand and an additional 1,045 Mgal/yr from the deep part of the Kirkwood-Cohansey aquifer system near the updip limit of the Atlantic City 800-foot sand. Most of the simulated decline in water levels in Atlantic County occurred as the result of the proposed increased withdrawals simulated for the New Jersey American Water Company wells. Simulated declines in water levels in Cape May were caused mainly by the simulated increased withdrawals for the Cape May City Desalination Plant wells. The additional water to supply the proposed increases in the scenarios was primarily horizontal flow from the unconfined updip part of the Kirkwood-Cohansey aquifer system, which accounted for 63 percent of the inflow, and flow from the overlying Kirkwood-Cohansey aquifer system into the Atlantic City 800-foot sand, which supplied 27 percent of the additional water. Because the withdrawals were made from the confined aquifer and the deeper part of the unconfined aquifer, the effect on streamflow was substantially less than would have occurred had the withdrawals been made directly from the shallower parts of the unconfined aquifer. The travel times from the 250-mg/L isochlor to production wells in Stone Harbor were longer as a result of all the additional withdrawals. For some scenarios, withdrawals in Atlantic County caused the saltwater to move slightly faster towards the production wells. These effects were offset by the increase in travel time caused by the potential increased withdrawals simulated for the Cape May City desalination wells, which either diverted water towards the desalination wells or increased the travel time towards production wells.

New Jersey↗
USGS1985

The ground-water system in the LaGrange Aquifer near LaGrange, southeastern Wyoming

Groundwater is being developed from the La Grange aquifer in southeastern Wyoming. It consists of saturated permeable alluvium that is hydraulically connected with most of the underlying White River Group. In the area of principal interest east of Horse Creek, Hawk Springs Reservoir and 14 adjacent wells used to supplement surface-water supply in the reservoir are in a natural discharge area. Upgradient of the reservoir there are 28 irrigation wells in about a 6-square-mile area. In this area, water levels declined between 3 and 12 feet from 1973 to 1978 causing concern about the effects of well pumpage on the hydrologic system. A digital model was developed and used to simulate the two-dimensional groundwater flow system in the unconfined La Grange aquifer. Transient simulations were made for 1973-78 using 12 time periods and for 1978-80 using 25 time periods. The calibrated digital model was used to simulate four 6-month pumping alternatives including three hypothetical alternatives for the area of principal interest east of Horse Creek. The reservoir altitude was held constant, approximating a reservoir volume 7,000 acre-feet. Pumping alternative 1 simulated historic conditions for 1973-78 including monthly recharge from precipitation which was included in the next three pumping alternatives. For pumping alternative 1, the calculated rate of discharge at the end of the 6-month simulation was 5.2 cubic feet per second from the aquifer to the reservoir. At the end of the 6-month simulations for pumping alternatives 2 and 3, the calculated rate of discharge to the reservoir was decreased to 0.4 cubic feet per second by pumpage from the 14 wells and to 3.8 cubic feet per second by pumpage from the 28 irrigation wells. For pumping alternative 4, pumpage from the total 42 wells resulted in a 1.0 cubic feet per second loss from the reservoir to the aquifer. (USGS)

Water-Resources Investigations Report↗
USGS2004

Simulated effects of ground-water augmentation on the hydrology of Round and Halfmoon Lakes in northwestern Hillsborough County, Florida

Pumpage from the Upper Floridan aquifer in northwest Hillsborough County near Tampa, Florida, has induced downward leakage from the overlying surficial aquifer and lowered the water table in many areas. Leakage is highest where the confining layer separating the aquifers is breached, which is common beneath many of the lakes in the study area. Leakage of water to the Upper Floridan aquifer has lowered the water level in many lakes and drained many wetlands. Ground water from the Upper Floridan aquifer has been added (augmented) to some lakes in an effort to maintain lake levels, but the resulting lake-water chemistry and lake leakage patterns are substantially different from those of natural lakes. Changes in lake-water chemistry can cause changes in lake flora, fauna, and lake sediment composition, and large volumes of lake leakage are suspected to enhance the formation of sinkholes near the shoreline of augmented lakes. The leakage rate of lake water through the surficial aquifer to the Upper Floridan aquifer was estimated in this study using ground-water-flow models developed for an augmented lake (Round Lake) and non-augmented lake (Halfmoon Lake). Flow models developed with MODFLOW were calibrated through nonlinear regression with UCODE to measured water levels and monthly net ground-water-flow rates from the lakes estimated from lake-water budgets. Monthly estimates of ground-water recharge were computed using an unsaturated flow model (LEACHM) that simulated daily changes in storage of water in the soil profile, thus estimating recharge as drainage to the water table. Aquifer properties in the Round Lake model were estimated through transient-state simulations using two sets of monthly recharge rates computed during July 1996 to February 1999, which spanned both average conditions (July 1996 through October 1997), and an El Niño event (November 1997 through September 1998) when the recharge rate doubled. Aquifer properties in the Halfmoon Lake model were estimated through steady-state simulations of average conditions in July 1996. Simulated hydrographs computed by the Round and Halfmoon Lake models closely matched measured water-level fluctuations, except during El Niño, when the Halfmoon Lake model was unable to accurately reproduce water levels. Possibly, potential recharge during El Niño was diverted through ground-water-flow outlets that were not represented in the Halfmoon Lake model, or a large part of the rainfall was diverted into runoff before it could become recharge. Solute transport simulations with MT3D indicate that leakage of lake water extended 250 to 400 feet into the surficial aquifer around Round Lake, and from 75 to 150 feet around Halfmoon Lake before flowing to the underlying Upper Floridan aquifer. These results are in agreement with concentrations of stable isotopes of oxygen-18 (d18O) and deuterium (dD) in the surficial aquifer. Schedules of monthly augmentation rates to maintain constant stages in Round and Halfmoon Lakes were computed using an equation that accounted for changes in the Upper Floridan aquifer head and the deviation from the mean recharge rate. Resulting lake stages were nearly constant during the first half of the study, but increased above target lake stages during El Niño; modifying the computation of augmentation rates to account for the higher recharge rate during El Niño resulted in lake stages that were closer to the target lake stage. Substantially more lake leakage flows to the Upper Floridan aquifer from Round Lake than from Halfmoon Lake, because the estimated vertical hydraulic conductivities of lake and confining layer sediments and breaches in the confining layer beneath Round Lake are much greater. Augmentation rates required to maintain the low guidance stages in Round Lake (53 feet) and Halfmoon Lake (42 feet) under average Upper Floridan aquifer heads are estimated as 33,850 cubic feet per day and 1,330 to 10,000 cubic feet per day, respectively. These rates equate to 26 inches per month of water applied to the entire surface of Round Lake and 0.34 to 2.5 inches per month of water applied to Halfmoon Lake.

Florida↗
USGS2017

A review of surface energy balance models for estimating actual evapotranspiration with remote sensing at high spatiotemporal resolution over large extents

Many approaches have been developed for measuring or estimating actual evapotranspiration ( ET a ), and research over many years has led to the development of remote sensing methods that are reliably reproducible and effective in estimating ET a . Several remote sensing methods can be used to estimate ET a at the high spatial resolution of agricultural fields and the large extent of river basins. More complex remote sensing methods apply an analytical approach to ET a estimation using physically based models of varied complexity that require a combination of ground-based and remote sensing data, and are grounded in the theory behind the surface energy balance model. This report, funded through cooperation with the International Joint Commission, provides an overview of selected remote sensing methods used for estimating water consumed through ET a and focuses on Mapping Evapotranspiration at High Resolution with Internalized Calibration (METRIC) and Operational Simplified Surface Energy Balance (SSEBop), two energy balance models for estimating ET a that are currently applied successfully in the United States. The METRIC model can produce maps of ET a at high spatial resolution (30 meters using Landsat data) for specific areas smaller than several hundred square kilometers in extent, an improvement in practice over methods used more generally at larger scales. Many studies validating METRIC estimates of ET a against measurements from lysimeters have shown model accuracies on daily to seasonal time scales ranging from 85 to 95 percent. The METRIC model is accurate, but the greater complexity of METRIC results in greater data requirements, and the internalized calibration of METRIC leads to greater skill required for implementation. In contrast, SSEBop is a simpler model, having reduced data requirements and greater ease of implementation without a substantial loss of accuracy in estimating ET a . The SSEBop model has been used to produce maps of ET a over very large extents (the conterminous United States) using lower spatial resolution (1 kilometer) Moderate Resolution Imaging Spectroradiometer (MODIS) data. Model accuracies ranging from 80 to 95 percent on daily to annual time scales have been shown in numerous studies that validated ET a estimates from SSEBop against eddy covariance measurements. The METRIC and SSEBop models can incorporate low and high spatial resolution data from MODIS and Landsat, but the high spatiotemporal resolution of ET a estimates using Landsat data over large extents takes immense computing power. Cloud computing is providing an opportunity for processing an increasing amount of geospatial “big data” in a decreasing period of time. For example, Google Earth Engine TM has been used to implement METRIC with automated calibration for regional-scale estimates of ET a using Landsat data. The U.S. Geological Survey also is using Google Earth Engine TM to implement SSEBop for estimating ET a in the United States at a continental scale using Landsat data.

Scientific Investigations Report↗
USGS1985

Hydrogeology, water quality, and ground-water development alternatives in the Beaver-Pasquiset ground-water reservoir, Rhode Island

In a 23 sq mi study area, the Beaver-Pasquiset groundwater reservoir within the Pawcatuck River basin in southern Rhode Island, stratified drift is the only principal geologic unit capable of producing yields > 350 gal/min. Transmissivity of the aquifer ranges from 7,200 to 24,300 sq ft/day. Water table conditions prevail in the aquifer, which is in good hydraulic connection with perennial streams and ponds. A digital model of two-dimensional groundwater flow was used to simulate the interaction between surface water and groundwater, and to evaluate the impact of alternative schemes of groundwater development on groundwater levels, pond levels, and streamflow in the Beaver-Pasquiset groundwater reservoir. Transient simulations of theoretical pumpage were made for a drought period (1963-66) and a wet period (1976-78). The areas most favorable for development of high-capacity wells (350 gal/min or more) are along the Beaver River and near Pasquiset Pond. The water is soft and generally contains < 100 mg/L dissolved solids. Locally, groundwater contains elevated concentrations of iron and manganese (7.5 and 3.7 mg/L, respectively), southeast of Pasquiset Pond, and will require treatment if used for public supply. The groundwater reservoir was simulated with a two-dimensional finite-difference model using a block-centered grid consisting of 33 rows and 75 columns. Differences between measured and simulated water table altitudes for the final steady state run for 21 selected observation wells averaged +0.07 ft. Combined pumping rates for simulation of groundwater development alternatives at eight sites ranged from 3.25 to 7.00 Mgal/d. Pumping rates for individual wells ranged from 0.25 to 1.50 Mgal/d. Transient simulations suggest that the Beaver-Pasquiset groundwater reservoir is capable of sustaining a pumping rate of 4.25 Mgal/d during years of average groundwater recharge with minimal impact on groundwater levels, pond levels, and streamflow. During extreme drought periods (1965 and 1966) it would be necessary to reduce pumpage below 3.25 Mgal/d to maintain flow in both the Beaver River and Pasquiset Brook. (Author 's abstract)

Water-Resources Investigations Report↗
USGS1988

Simulation of ground-water flow in the lower sand unit of the Potomac-Raritan-Magothy aquifer system, Philadelphia, Pennsylvania

Ground-water flow in the lower sand unit of the Potomac-Raritan-Magothy aquifer system in Philadelphia was simulated with a two-dimensional finite- difference ground-water model. The modeled 133-square-mile area also included parts of Delaware County, Pennsylvania, and Camden and Gloucester Counties, New Jersey. The lower sand unit is Cretaceous in age and consists of well- sorted coarse sand and fine gravel that grades upward into medium to fine sand containing a few thin beds of clay. The modeled aquifer consists of the lower sand unit in Philadelphia and the lowermost sand unit of the Potomac-Raritan-Magothy aquifer system in New Jersey. Throughout most of the area, the lower sand unit is overlain by a clay confining unit. Where the clay is absent, the lower sand unit is unconfined. A hydraulic conductivity of 1.6 x 10-3 foot per second and a storage coefficient of 3.0 x 10-4 was assigned to the lower sand unit based on 15 aquifer tests, and a hydraulic conductivity of 4.0 x 10-8 foot per second was assigned to the upper confining unit based on transient-flow sensitivity analysis. Water levels were not sensitive to changes in the value for specific storage of the upper confining unit, indicating that most vertical leakage occurs as steady leakage. Changes in the potentiometric surface of the lower sand unit for 1904-78 simulated. Differences between simulated and observed head generally were less than 10 feet. Simulations were made to determine the effects on hydraulic head of increases in industrial pumpage of 5 and 10 Mgal/d (million gallons per day) and of an emergency 60 Mgal/d municipal water supply in Philadelphia. A 5- and 10-Mgal/d increase in industrial pumpage would lower heads in the lower sand unit by as much as 33 and 66 feet, respectively. Pumping 60 Mgal/d for 30 days for an emergency municipal supply would lower heads in the lower sand unit by as much as 121 feet.

Water-Resources Investigations Report↗
USGS2005

Water budgets for selected watersheds in the Delaware River basin, eastern Pennsylvania and western New Jersey

This pilot study, done by the U.S. Geological Survey in cooperation with the Delaware River Basin Commission, developed annual water budgets using available data for five watersheds in the Delaware River Basin with different degrees of urbanization and different geological settings. A basin water budget and a water-use budget were developed for each watershed. The basin water budget describes inputs to the watershed (precipitation and imported water), outputs of water from the watershed (streamflow, exported water, leakage, consumed water, and evapotranspiration), and changes in ground-water and surface-water storage. The water-use budget describes water withdrawals in the watershed (ground-water and surface-water withdrawals), discharges of water in the watershed (discharge to surface water and ground water), and movement of water of water into and out of the watershed (imports, exports, and consumed water). The water-budget equations developed for this study can be applied to any watershed in the Delaware River Basin. Data used to develop the water budgets were obtained from available long-term meteorological and hydrological data-collection stations and from water-use data collected by regulatory agencies. In the Coastal Plain watersheds, net ground-water loss from unconfined to confined aquifers was determined by using ground-water-flow-model simulations. Error in the water-budget terms is caused by missing data, poor or incomplete measurements, overestimated or underestimated quantities, measurement or reporting errors, and the use of point measurements, such as precipitation and water levels, to estimate an areal quantity, particularly if the watershed is hydrologically or geologically complex or the data-collection station is outside the watershed. The complexity of the water budgets increases with increasing watershed urbanization and interbasin transfer of water. In the Wissahickon Creek watershed, for example, some ground water is discharged to streams in the watershed, some is exported as wastewater, and some is exported for public supply. In addition, ground water withdrawn outside the watershed is imported for public supply or imported as wastewater for treatment and discharge in the watershed. A GIS analysis was necessary to quantify many of the water-budget components. The 89.9-square mile East Branch Brandywine Creek watershed in Pennsylvania is a rural watershed with reservoir storage that is underlain by fractured rock. Water budgets were developed for 1977-2001. Average annual precipitation, streamflow, and evapotranspiration were 46.89, 21.58, and 25.88 inches, respectively. Some water was imported (average of 0.68 inches) into the watershed for public-water supply and as wastewater for treatment and discharge; these imports resulted in a net gain of water to the watershed. More water was discharged to East Branch Brandywine Creek than was withdrawn from it; the net discharge resulted in an increase in streamflow. Most ground water was withdrawn (average of 0.25 inches) for public-water supply. Surface water was withdrawn (average of 0.58 inches) for public-water and industrial supply. Discharge of water by sewage-treatment plants and industries (average of 1.22 inches) and regulation by Marsh Creek Reservoir caused base flow to appear an average of 7.2 percent higher than it would have been without these additional sources. On average, 67 percent of the difference was caused by sewage-treatment-plant and industrial discharges, and 33 percent was caused by regulation of the Marsh Creek Reservoir. Water imports, withdrawals, and discharges have been increasing as the watershed becomes increasingly urbanized. The 64-square mile Wissahickon Creek watershed in Pennsylvania is an urban watershed underlain by fractured rock. Water budgets were developed for 1987-98. Average annual precipitation, streamflow, and evapotranspiration were 47.23, 22.24, and 23.12 inches, respectively. The watershed is highly u

Pennsylvania↗
USGS1993

A modular finite-element model (MODFE) for areal and axisymmetric ground-water-flow problems, Part 3: Design philosophy and programming details

A MODular Finite-Element, digital-computer program (MODFE) was developed to simulate steady or unsteady-state, two-dimensional or axisymmetric ground-water-flow. The modular structure of MODFE places the computationally independent tasks that are performed routinely by digital-computer programs simulating ground-water flow into separate subroutines, which are executed from the main program by control statements. Each subroutine consists of complete sets of computations, or modules, which are identified by comment statements, and can be modified by the user without affecting unrelated computations elsewhere in the program. Simulation capabilities can be added or modified by either adding or modifying subroutines that perform specific computational tasks, and the modular-program structure allows the user to create versions of MODFE that contain only the simulation capabilities that pertain to the ground-water problem of interest. MODFE is written in a Fortran programming language that makes it virtually device independent and compatible with desk-top personal computers and large mainframes. MODFE uses computer storage and execution time efficiently by taking advantage of symmetry and sparseness within the coefficient matrices of the finite-element equations. Parts of the matrix coefficients are computed and stored as single-subscripted variables, which are assembled into a complete coefficient just prior to solution. Computer storage is reused during simulation to decrease storage requirements. Descriptions of subroutines that execute the computational steps of the modular-program structure are given in tables that cross reference the subroutines with particular versions of MODFE. Programming details of linear and nonlinear hydrologic terms are provided. Structure diagrams for the main programs show the order in which subroutines are executed for each version and illustrate some of the linear and nonlinear versions of MODFE that are possible. Computational aspects of changing stresses and boundary conditions with time and of mass-balance and error terms are given for each hydrologic feature. Program variables are listed and defined according to their occurrence in the main programs and in subroutines. Listings of the main programs and subroutines are given.

Techniques of Water-Resources Investigations↗
USGS1976

Availability of ground water near Carmel, Hamilton County, Indiana

A study of the hydraulic characteristics of the unconsolidated glacial deposits near the city of Carmel in central Indiana shows that 21.3 million gallons per day (933 litres per second) of additional water could be withdrawn from the aquifer for an indefinite period of time . This pumpage is approximately 5 million gallons per day (219 litres per second) above the projected water needs of Carmel for 1990. Saturated thickness, transmissivity , and storage coefficient of the outwash aquifer along the White River east of Carmel were determined , using available data supplemented by test drilling. The saturated thickness of the aquifer ranges from 10 to 110 feet (3 to 34 metres); transmissivity ranges from 1 , 000 feet squared per day (93 metres squared per day) to 24,000 feet squared per day (2,230 metres squared per day); and the average storage coefficient is 0 . 11. Seepage from the aquifer into the White River was estimated in November 1974, using data from u.S. Geological Survey gaging stations. Water- level information was obtained from a network of observation wells at that same time. Flow in the unconsolidated glacial deposits near the city of Carmel in central Indiana was simulated by a digital-computer model in a study of hydraulic characteristics of the deposits. The study shows that 21.3 million gallons per day (933 litres per second) of additional water could be withdrawn from the aquifer for an indefinite period of time. This pumpage is approximately 5 million gallons per day (219 1itres per second) above the projected water needs of Carmel for 1990. Saturated thickness, transmissivity, and storage coefficient of the outwash aquifer along the White River east of Carmel were determined, using available data supplemented by test drilling . The saturated thickness of the aquifer ranges f r om 10 to 110 feet 0 to 34 me tres); transmissivity ranges from 1,000 feet squared per day (93 metres squared per day) to 24 ,000 feet squared per day (2 ,230 metres squared per day); and the average storage coefficient is 0.11.

Indiana↗
USGS2004

Delineation of areas contributing recharge to selected public-supply wells in Glacial Valley-Fill and Wetland Settings, Rhode Island

Areas contributing recharge and sources of water to one proposed and seven present public-supply wells, screened in sand and gravel deposits and clustered in three study areas, were determined on the basis of calibrated, steady-state ground-water-flow models representing average hydrologic conditions. The area contributing recharge to a well is defined as the surface area where water recharges the ground water and then flows toward and discharges to the well. In Cumberland and Lincoln, public-supply well fields on opposite sides of the Blackstone River are in a narrow valley bordered by steep hillslopes. Ground-water-level and river-stage measurements indicated that river water was infiltrating the aquifer and flowing toward the wells during pumping conditions. Simulated areas contributing recharge to the Cumberland well field operating alone for both average (324 gallons per minute) and maximum (1,000 gallons per minute) pumping rates extend on both sides of the river to the lateral model boundaries, which is the contact between the valley and uplands. The area contributing recharge at the average pumping rate is about 0.05 square mile and the well field derives 72 percent of pumped water from upland runoff. At the maximum pumping rate, the area contributing recharge extends farther up and down the valley to 0.12 square mile and the primary source of water to the well field was infiltrated river water (53 percent). Upland areas draining toward the areas contributing recharge encompass 0.58 and 0.66 square mile for the average and maximum rates, respectively. By incorporating the backup Lincoln well-field withdrawals (2,083 gallons per minute) into the model, the area contributing recharge to the Cumberland well field operating at its maximum rate is reduced to 0.08 square mile; part of the simulated area which contributes recharge to the Cumberland well field when it is operating alone contributes instead to the Lincoln well field when both well fields are pumped. The Cumberland well field compensates by increasing the percentage of water it withdraws from the river by 11 percent. The upland area draining toward the Cumberland contributing area is 0.55 square mile. The area contributing recharge to the Lincoln well field is 0.08 square mile and infiltrated river water contributes 88 percent of the total water; the upland area draining toward the contributing area is 0.34 square mile. In North Smithfield, a public-supply well in a valley-fill setting is close to Trout Brook Pond, which is an extension of the Lower Slatersville Reservoir. A comparison of water levels from the pond and underlying sediments indicates that water is not infiltrated from Trout Brook Pond when the supply well is pumped at its maximum rate of 200 gallons per minute. Simulated areas contributing recharge for the maximum pumping rate and for the estimated maximum yield, 500 gallons per minute, of a proposed replacement well extend to the ground-water divides on both sides of Trout Brook Pond. For the 200 gallons-per-minute rate, the area contributing recharge is 0.23 square mile; the well derives almost all of its water from intercepted ground water that normally discharges to surface-water bodies. For the pumping rate of 500 gallons per minute, the area contributing recharge is 0.45 square mile. The increased pumping rate is balanced by additional intercepted ground water and by inducing 25 percent of the total withdrawn water from surface water. In Westerly, one public-supply well is in a watershed where the primarily hydrologic feature is a wetland. Water levels in piezometers surrounding the well site indicated a downward vertical gradient and the potential for water in the wetland to infiltrate the underlying aquifer. The simulated area contributing recharge for the average pumping rate (240 gallons per minute) and for the maximum pumping rate (700 gallons per minute) extends to the surrounding uplands (surficial materials not covered by t

Scientific Investigations Report↗
USGS1995

Physical and hydrologic properties of outcrop samples from a nonwelded to welded tuff transition, Yucca Mountain, Nevada

Quantitative material-property data are needed to describe lateral and vertical spatial variability of physical and hydrologic properties and to model ground-water flow and radionuclide transport at the potential Yucca Mountain nuclear-waste repository site in Nevada. As part of ongoing site characterization studies of Yucca Mountain directed toward this understanding of spatial variability, laboratory measurements of porosity, bull* and particle density, saturated hydraulic conductivity, and sorptivity have been obtained for a set of outcrop samples that form a systematic,two dimensional grid that covers a large exposure of the basal Tiva Canyon Tuff of the Paintbrush Group of Miocene age at Yucca Mountain. The samples form a detailed vertical grid roughly parallel to the transport direction of the parent ash flows, and they exhibit material-property varia- tions in an interval of major lithologic change overlying a potential nuclear-waste repository at Yucca Mountain. The observed changes in hydrologic properties were systematic and consistent with the changes expected for the nonwelded to welded transition at the base of a major ash-flow sequence. Porosity, saturated hydraulic conductivity, and sorptivity decreased upward from the base of the Tiva Canyon Tuff, indicating the progressive compaction of ash- rich volcanic debris and the onset of welding with increased overburden pressure from the accumulating ash-flow sheet. The rate of decrease in the values of these material properties varied with vertical position within the transition interval. In contrast, bulk-density values increased upward, a change that also is consistent with progressive compaction and the onset of welding. Particle-density values remained almost constant throughout the transition interval, probably indicating compositional (chemical) homogeneity.

Water-Resources Investigations Report↗
USGS1987

Relation between ground water and surface water in the Hillsborough River basin, west-central Florida

The relation between groundwater and surface water in the Hillsborough River basin was defined through the use of: seismic-reflection profiling along selected reaches of the Hillsborough River, and evaluation of streamflow, rainfall, groundwater levels, water quality, and geologic data. Major municipal well fields in the basin are Morris Bridge and Cypress Creek where an averages of 15.3 and 30.0 million gal/day (mgd), respectively, were pumped in 1980. Mean annual rainfall for the study area is 53.7 inches. Average rainfall for 1980, determined from eight rainfall stations, was 49.7 inches. Evapotranspiration, corrected for the 5% of the basin that is standing water, was 35.7 in/year. The principal geohydrologic units in the basin are the surficial aquifer, the intermediate aquifer and confining beds, the Upper Floridan aquifer, the middle confining unit, and the Lower Floridan aquifer. Total pumpage of groundwater in 1980 was 98.18 mgd. The surficial aquifer and the intermediate aquifer are not used for major groundwater supply in the basin. Continuous marine seismic-reflection data collected along selected reaches of the Hillsborough River were interpreted to define the riverbed profile, the thickness of surficial deposits, and the top of persistent limestone. Major areas of groundwater discharge near the Hillsborough River and its tributaries are the wetlands adjacent to the river between the Zephyrhills gaging stations and Fletcher Avenue and the wetlands adjacent to Cypress Creek. An estimated 20 mgd seeps upward from the Upper Floridan aquifer within those wetland areas. The runoff/sq mi is greater at the Zephyrhills station than at Morris Bridge. However, results of groundwater flow models and potentiometric-surface maps indicate that groundwater is flowing upward along the Hillsborough River between the Zephyrhills gage and the Morris Bridge gage. This upward leakage is lost to evapotranspiration. An aquifer test conducted in 1978 at the Morris Bridge well field was evaluated by using an anisotropic method. Analytical results matched observed water levels within 0.1 ft. Analysis of aquifer test results indicates that withdrawals of up to 28 mgd would have a negligible effect on the river stage or flow. (Author 's abstract)

Water-Resources Investigations Report↗
USGS1993

Geohydrology, water quality, and estimation of ground-water recharge in San Francisco, California, 1987-92

The city of San Francisco is considering further development of local groundwater resources as a supplemental source of water for potable or nonpotable use. By the year 2010, further water demand is projected to exceed the delivery capacity of the existing supply system, which is fed by surface-water sources; thus supplies are susceptible to drought conditions and damage to conveyance lines by earthquakes. The primary purpose of this study is to describe local geohydrology and water quality and to estimate groundwater recharge in the area of the city of San Francisco. Seven groundwater basins were identified in San Francisco on the basis of geologic and geophysical data. Basins on the east side of the city are relatively thin and contain a greater percentage of fine-grained sediments than those on the west side. The relatively small capacity of the basins and greater potential for contamination from sewer sources may limit the potential for groundwater development on the east side. Basins on the west side of the city have a relatively large capacity and low density sewer network. Water-level data indicate that the southern part of the largest basin on the west side of the city (Westside basin) probably cannot accommodate additional groundwater development without adversely affecting water levels and water quality in Lake Merced; however, the remainder of the basin, which is largely undeveloped, could be developed further. A hydrologic routing model was developed for estimating groundwater recharge throughout San Francisco. The model takes into account climatic factors, land and water use, irrigation, leakage from underground pipes, rainfall runoff, evapotranspiration, and other factors associated with an urban environment. Results indicate that area recharge rates for water years 1987-88 for the 7 groundwater basins ranged from 0.32 to 0.78 feet per year. Recharge for the Westside basin was estimated at 0.51 feet per year. Average annual groundwater recharge represents the maximum annual long-term yield of the basin. Attainable yield may be less than the volume of groundwater recharge because interception of all discharge from the basin may not be feasible without inducing seawater intrusion or causing other undesirable effects.

Water-Resources Investigations Report↗