Geology ReportsSearch

SEARCH · Geology Reports

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 541 records · Page 30Linked to original sources

Statistical and simulation analysis of hydraulic-conductivity data for Bear Creek and Melton Valleys, Oak Ridge Reservation, Tennessee

A total of 338 single-well aquifer tests from Bear Creek and Melton Valley, Tennessee were statistically grouped to estimate hydraulic conductivities for the geologic formations in the valleys. A cross-sectional simulation model linked to a regression model was used to further refine the statistical estimates for each of the formations and to improve understanding of ground-water flow in Bear Creek Valley. Median hydraulic-conductivity values were used as initial values in the model. Model-calculated estimates of hydraulic conductivity were generally lower than the statistical estimates. Simulations indicate that (1) the Pumpkin Valley Shale controls groundwater flow between Pine Ridge and Bear Creek; (2) all the recharge on Chestnut Ridge discharges to the Maynardville Limestone; (3) the formations having smaller hydraulic gradients may have a greater tendency for flow along strike; (4) local hydraulic conditions in the Maynardville Limestone cause inaccurate model-calculated estimates of hydraulic conductivity; and (5) the conductivity of deep bedrock neither affects the results of the model nor does it add information on the flow system. Improved model performance would require: (1) more water level data for the Copper Ridge Dolomite; (2) improved estimates of hydraulic conductivity in the Copper Ridge Dolomite and Maynardville Limestone; and (3) more water level data and aquifer tests in deep bedrock.

Tennessee

Water use, ground-water recharge and availability, and quality of water in the Greenwich area, Fairfield County, Connecticut and Westchester County, New York, 2000-2002

Ground-water budgets were developed for 32 small basin-based zones in the Greenwich area of southwestern Connecticut, where crystalline-bedrock aquifers supply private wells, to determine the status of residential ground-water consumption relative to rates of ground-water recharge and discharge. Estimated residential ground-water withdrawals for small basins (averaging 1.7 square miles (mi 2 )) ranged from 0 to 0.16 million gallons per day per square mile (Mgal/d/mi 2 ). To develop these budgets, residential ground-water withdrawals were estimated using multiple-linear regression models that relate water use from public water supply to data on residential property characteristics. Average daily water use of households with public water supply ranged from 219 to 1,082 gallons per day (gal/d). A steady-state finite-difference ground-water- flow model was developed to track water budgets, and to estimate optimal values for hydraulic conductivity of the bedrock (0.05 feet per day) and recharge to the overlying till deposits (6.9 inches) using nonlinear regression. Estimated recharge rates to the small basins ranged from 3.6 to 7.5 inches per year (in/yr) and relate to the percentage of the basin underlain by coarse- grained glacial stratified deposits. Recharge was not applied to impervious areas to account for the effects of urbanization. Net residential ground-water consumption was estimated as ground-water withdrawals increased during the growing season, and ranged from 0 to 0.9 in/yr. Long-term average stream base flows simulated by the ground-water-flow model were compared to calculated values of average base flow and low flow to determine if base flow was substantially reduced in any of the basins studied. Three of the 32 basins studied had simulated base flows less than 3 in/yr, as a result of either ground-water withdrawals or reduced recharge due to urbanization. A water-availability criteria of the difference between the 30-day 2-year low flow and the recharge rate for each basin was explored as a method to rate the status of water consumption in each basin. Water consumption ranged from 0 to 14.3 percent of available water based on this criteria for the 32 basins studied. Base-flow water quality was related to the amount of urbanized area in each basin sampled. Concentrations of total nitrogen and phosphorus, chloride, indicator bacteria, and the number of pesticide detections increased with basin urbanization, which ranged from 18 to 63 percent of basin area.

Connecticut, New York

Simulation of effects of ground-water development on water-levels in glacial-drift aquifers in the Brooten-Belgrade area, west-central Minnesota

Ground-water flow in the confined- and unconfined-drift aquifers in the Brooten-Belgrade area of west-central Minnesota was simulated with a three-dimensional finite-difference ground-water-flow model. Model results indicate that about 96 percent of the total inflow to the modeled area is from precipitation. Discounting evapotranspiration, 63 percent of the total outflow is ground-water discharge to the East Branch Chippewa and North Fork Crow Rivers, and 34 percent is ground-water pumpage. The model was used to simulate the steady-state effects of below-normal precipitation (drought) and hypothetical increases in ground-water development. Model results indicate that reduced recharge and increased pumping during a hypothetical 3-year extended drought would lower regional water levels from 2 to 5 feet in each aquifer and as much as 20 feet in the lowermost aquifer zone; ground-water discharge to the East Branch Chippewa and North Fork Crow Rivers would be reduced by 38 percent. The addition of 10 to 20 hypothetical wells in confined aquifers, pumping 123 to 246 million gallons per year, would result in regional water-level declines of 0.1 to 0.5 feet. Simulated water-level declines in wells completed in the lower part of the system would be as much as 5.0 feet as a result of pumping 246 million gallons per year from 20 hypothetical wells. Water-level declines in overlying and underlying aquifers would range from 0.4 to 2.8 feet. Ground-water discharge to the East Branch Chippewa and North Fork Crow Rivers would be unaffected by the pumpage.

Minnesota

Potential for, and possible effects of, artificial recharge in Carson Valley, Douglas County, Nevada

Rapid population growth in Carson Valley, west- central Nevada, requires a dependable municipal water source. Artificial recharge of aquifers using available flow of the Carson River is one way to increase the amount of water in underground storage and maintain a dependable ground-water supply. Ground water can be artificially recharged by routing excess surface water or, after proper treatment, routing wastewater to infiltration basins or injection wells. Withdrawal wells would remove stored water when needed. As a first step, maps showing areas in Carson Valley with high, low, moderate and unknown potential for artificial recharge were developed on the basis of the distribution of geologic units, depth to water, specific yield, infiltration rate, and location of natural recharge and discharge. For recharge by means of infiltration, areas totaling 5,700 acres have high potential, 23,900 acres have moderate potential, and 6,200 acres have low potential. For recharge through injection, areas totaling 7,800 acres have high potential and 43,500 acres have moderate potential; 23,000 acres have unknown potential because data are lacking on subsurface conditions. A ground-water-flow model was used to assess the possible results of artificial recharge. Simulations with no accompanying ground-water withdrawal show that, when recharge by injection is simulated near the valley floor, heads in the semiconfined aquifer increase over much of the valley, floor; only about 20 percent of the recharged water is stored in the aquifer after 5 years and as much as 80 percent is lost to streamflow and evapotranspiration. When recharge is simulated on the eastern side of the valley, 80 percent of the recharged water remains in storage after 5 years. When recharge is simulated near the valley floor, more water is lost to discharge than when recharge is on the eastern side of the valley. When recharge is applied for long periods without accompanying withdrawal, recharged water moves downgradient to discharge areas. The recharge water that discharges to the surface-water system could in turn replenish base flow of the Carson River and benefit downstream users.

Water-Resources Investigations Report

SUTRA: A model for 2D or 3D saturated-unsaturated, variable-density ground-water flow with solute or energy transport

SUTRA (Saturated-Unsaturated Transport) is a computer program that simulates fluid movement and the transport of either energy or dissolved substances in a subsurface environment. This upgraded version of SUTRA adds the capability for three-dimensional simulation to the former code (Voss, 1984), which allowed only two-dimensional simulation. The code employs a two- or three-dimensional finite-element and finite-difference method to approximate the governing equations that describe the two interdependent processes that are simulated: 1) fluid density-dependent saturated or unsaturated ground-water flow; and 2) either (a) transport of a solute in the ground water, in which the solute may be subject to: equilibrium adsorption on the porous matrix, and both first-order and zero-order production or decay; or (b) transport of thermal energy in the ground water and solid matrix of the aquifer. SUTRA may also be used to simulate simpler subsets of the above processes. A flow-direction-dependent dispersion process for anisotropic media is also provided by the code and is introduced in this report. As the primary calculated result, SUTRA provides fluid pressures and either solute concentrations or temperatures, as they vary with time, everywhere in the simulated subsurface system. SUTRA flow simulation may be employed for two-dimensional (2D) areal, cross sectional and three-dimensional (3D) modeling of saturated ground-water flow systems, and for cross sectional and 3D modeling of unsaturated zone flow. Solute-transport simulation using SUTRA may be employed to model natural or man-induced chemical-species transport including processes of solute sorption, production, and decay. For example, it may be applied to analyze ground-water contaminant transport problems and aquifer restoration designs. In addition, solute-transport simulation with SUTRA may be used for modeling of variable-density leachate movement, and for cross sectional modeling of saltwater intrusion in aquifers at near-well or regional scales, with either dispersed or relatively sharp transition zones between freshwater and saltwater. SUTRA energy-transport simulation may be employed to model thermal regimes in aquifers, subsurface heat conduction, aquifer thermal-energy storage systems, geothermal reservoirs, thermal pollution of aquifers, and natural hydrogeologic convection systems. Mesh construction, which is quite flexible for arbitrary geometries, employs quadrilateral finite elements in 2D Cartesian or radial-cylindrical coordinate systems, and hexahedral finite elements in 3D systems. 3D meshes are currently restricted to be logically rectangular; in other words, they are similar to deformable finite-difference-style grids. Permeabilities may be anisotropic and may vary in both direction and magnitude throughout the system, as may most other aquifer and fluid properties. Boundary conditions, sources and sinks may be time dependent. A number of input data checks are made to verify the input data set. An option is available for storing intermediate results and restarting a simulation at the intermediate time. Output options include fluid velocities, fluid mass and solute mass or energy budgets, and time-varying observations at points in the system. Both the mathematical basis for SUTRA and the program structure are highly general, and are modularized to allow for straightforward addition of new methods or processes to the simulation. The FORTRAN-90 coding stresses clarity and modularity rather than efficiency, providing easy access for later modifications.

Water-Resources Investigations Report

Evaluation of saltwater intrusion and travel time in the Atlantic City 800-foot sand, Cape May County, New Jersey, 1992, by use of a coupled-model approach and flow-path analysis

Regional and subregional ground-water-flow models were coupled, and the output was analyzed by a particle-tracking method. The results were then used to assess the effects of ground-water withdrawals on the flow of saltwater in the Atlantic City 800-foot sand in Cape May County, New Jersey, and to estimate the travel time from areas in which the chloride concentration of the ground water exceeds 250 milligrams per liter to the county's nearest public supply wells. First, a quasi-three-dimensional finite-difference computer model of freshwater and saltwater flow that simulated regional ground-water flow through the unconsolidated materials underlying the New Jersey Coastal Plain was used to estimate flow at the boundaries of the subregional study area. The results of the regional simulation were used as input to a second quasi-three-dimensional finite-difference model that was used to simulate flow in the subregion, the Atlantic City 800-foot sand in Cape May County. The results of the simulation of flow in the subregion were analyzed by a semianalytical particle-tracking method to estimate ground-water flow paths and travel time of ground water from areas in which chloride concentrations exceed 250 milligrams per liter to public supply wells located at Stone Harbor, New Jersey. Ground-water withdrawals from the Atlantic City 800-foot sand were assumed to be equal to those reported for 1991. Results of the analysis indicate that the time required for saltwater to reach the public supply wells is on the order of hundreds of years. These results, however, are based on the assumptions that the aquifer is homogeneous. The presence of zones of high permeability in the aquifer could reduce the predicted travel times of the saltwater from its present location to the supply wells. Travel times also could be reduced if ground-water withdrawals increase.

New Jersey

Hydrogeologic setting, ground-water flow, and ground-water quality at the Lake Wheeler Road research station, 2001-03 : North Carolina Piedmont and Mountains Resource Evaluation Program

Results of a 2-year field study of the regolith-fractured bedrock ground-water system at the Lake Wheeler Road research station in Wake County, North Carolina, indicate both disconnection and interaction among components of the ground-water system. The three components of the ground-water system include (1) shallow, porous regolith; (2) a transition zone, including partially weathered rock, having both secondary (fractures) and primary porosity; and (3) deeper, fractured bedrock that has little, if any, primary porosity and is dominated by secondary fractures. The research station includes 15 wells (including a well transect from topographic high to low settings) completed in the three major components of the ground-water-flow system and a surface-water gaging station on an unnamed tributary. The Lake Wheeler Road research station is considered representative of a felsic gneiss hydrogeologic unit having steeply dipping foliation and a relatively thick overlying regolith. Bedrock foliation generally strikes N. 10? E. to N. 30? E. and N. 20? W. to N. 40? W. to a depth of about 400 feet and dips between 70? and 80? SE. and NE., respectively. From 400 to 600 feet, the foliation generally strikes N. 70? E. to N. 80? E., dipping 70? to 80? SE. Depth to bedrock locally ranges from about 67 to 77 feet below land surface. Fractures in the bedrock generally occur in two primary sets: low dip angle, stress relief fractures that cross cut foliation, and steeply dipping fractures parallel to foliation. Findings of this study generally support the conceptual models of ground-water flow from high to low topographic settings developed for the Piedmont and Blue Ridge Provinces in previous investigations, but are considered a refinement of the generalized conceptual model based on a detailed local-scale investigation. Ground water flows toward a surface-water boundary, and hydraulic gradients generally are downward in recharge areas and upward in discharge areas; however, local variations in vertical gradients are apparent. Water-quality sampling and monitoring efforts were conducted to characterize the interaction of components of the ground-water system. Elevated nitrate concentrations as high as 22 milligrams per liter were detected in shallow ground water from the regolith at the study site. These elevated nitrate concentrations likely are related to land use, which includes agricultural practices that involve animal feeding operations and crop fertilization. Continuous ground-water-quality data indicate seasonal fluctuations in field water-quality properties, differences with respect to depth, and fluctuations during recharge events. Water-quality properties recorded in the regolith well following rainfall indicate the upwelling of deeper ground water in the discharge area, likely from ground water in the transition-zone fractures. Additionally, interaction with a surface-water boundary appears likely in the ground-water discharge area, as water levels in all three ground-water zones, including the deep bedrock, mimic the surface-water rise during rainfall.

North Carolina

Variability of differences between two approaches for determining ground-water discharge and pumpage, including effects of time trends, Lower Arkansas River Basin, southeastern Colorado, 1998-2002

In the mid-1990s, the Colorado Division of Water Resources (CDWR) adopted rules governing measurement of tributary ground-water pumpage for the Arkansas River Basin. The rules allowed ground-water pumpage to be determined using one of two approaches?power conversion coefficient (PCC) or totalizing flowmeters (TFM). In addition, the rules allowed a PCC to be applied to the electrical power usage up to 4 years in the future to estimate ground-water pumpage. As a result of concerns about potential errors in applying the PCC approach forward in time, a study was done by the U.S. Geological Survey, in cooperation with CDWR and Colorado Water Conservation Board, to evaluate the variability in differences in pumpage between the two approaches, including the effects of time trends. This report compared measured ground-water pumpage using TFMs to computed ground-water pumpage using PCCs by developing statistical models of relations between explanatory variables, such as site, time, and pumping water level, and dependent variables, which are based on discharge, PCC, and pumpage. When differences in pumpage (diffP) were computed using PCC measurements and power consumption for the same year (1998-2002), the median diffP, depending on the year, ranged from +0.1 to -2.9 percent; the median diffP for the entire period was -1.5 percent. However, when diffP was computed using PCC measurements applied to the next year's power consumption, the median diffP was -0.3 percent; and when PCC measurements were applied 2, 3, or 4 years into the future, median diffPs were +1.8 percent for a 2-year forward lag and +5.3 percent for a 4-year forward lag, indicating that pumpage computed with the PCC approach, as generally applied under the ground-water pumpage measurement rules by CDWR, tended to overestimate pumpage as compared to pumpage using TFMs when PCC measurement was applied to future years of measured power consumption. Analyses were done to better understand the causes of the time trend; an estimate of the overall trend with time (uncorrected for pumping water-level changes) yielded a trend of about 2.2 percent per lag year for diffP. A separate analysis that incorporated a surface-water diversion term in the statistical model rendered the time-trend term insignificant, indicating that the time trend in the models served as a surrogate for other variables, some of which reflect underlying hydrologic conditions. A more precise explanation of the potential causes of the time trend was not obtained with the available data. However, the model results with the surface-water diversion term indicate that much of the trend of 2.2 percent per lag year in diffP resulted from applying a PCC to estimate pumpage under hydrologic conditions different from those under which the PCC was measured. Although there is no evidence to conclude that the upward time trend determined in the data for this 5-year period would hold in the future, historical static ground-water levels in the study area generally have exhibited small variations over multidecadal time scales. Therefore, the approximately 2 percent per lag year trend determined in these data is expected to be a reasonable guideline for estimating potential errors in the PCC approach resulting from temporally varying hydrologic conditions between time of PCC measurement and pumpage estimation. Comparisons also were made between total, or aggregated, pumpage for a network of wells as computed by the PCC approach and the TFM approach. For 100 wells and a lag of 4 years between PCC measurement and pumpage estimation, there was a 95-percent probability that the difference between total network pumpage measured by the PCC approach and that measured using a TFM would be between 5.2 and 14.4 percent. These estimates were based on a bias of 2.2 percent per lag year estimated for the period 1998-2002 during which hydrologic conditions were known to have changed. Using the same assumptions, the estimated d

Scientific Investigations Report

Ground-water flow model of the Boone formation at the Tar Creek superfund site, Oklahoma and Kansas

Extensive mining activities conducted at the Tar Creek Superfund site, one of the largest Superfund sites in the United States, pose substantial health and safety risks. Mining activities removed a total of about 6,000,000 tons of lead and zinc by 1949. To evaluate the effect of this mining on the ground-water flow, a MODFLOW 2000 digital model has been developed to simulate ground-water flow in the carbonate formations of Mississippian age underlying the Tar Creek Superfund site. The model consists of three layers of variable thickness and a grid of 580 rows by 680 columns of cells 164 feet (50 meters) on a side. Model flux boundary conditions are specified for rivers and general head boundaries along the northern boundary of the Boone Formation. Selected cells in layer 1 are simulated as drain cells. Model calibration has been performed to minimize the difference between simulated and observed water levels in the Boone Formation. Hydraulic conductivity values specified during calibration range from 1.3 to 35 feet per day for the Boone Formation with the larger values occurring along the axis of the Miami Syncline where horizontal anisotropy is specified as 10 to 1. Hydraulic conductivity associated with the mine void is set at 50,000 feet per day and a specific yield of 1.0 is specified to represent that the mine void is filled completely with water. Residuals (the difference between measured and simulated ground-water altitudes) has a root-mean-squared value of 8.53 feet and an absolute mean value of 7.29 feet for 17 observed values of water levels in the Boone Formation. The utility of the model for simulating and evaluating the possible consequences of remediation activities has been demonstrated. The model was used to simulate the emplacement of chat (mine waste consisting of fines and fragments of chert) back into the mine. Scenarios using 1,800,000 and 6,500,000 tons of chat were run. Hydraulic conductivity was reduced from 50,000 feet per day to 35 feet per day in the model cells corresponding to chat emplacement locations. A comparison of the simulated baseline conditions and conditions after simulated chat emplacement revealed little change in water levels, drainage and stream flux, and ground-water flow velocity. Using the calibrated flow model, particle tracks were simulated using MODPATH to evaluate the simultaneous movement of particles with water in the vicinity of four potential sites at which various volumes of chat might be emplaced in the underground mine workings as part of potential remediation efforts at the site. Particle tracks were generated to follow the rate and direction of water movement for a simulated period of 100 years. In general, chat emplacement had minimal effect on the direction and rate of movement when compared to baseline (current) flow conditions. Water-level differences between baseline and chat-emplacement scenarios showed declines as much as 2 to 3 feet in areas immediately downgradient from the chat emplacement cells and little or no head change upgradient. Chat emplacements had minimal effect on changes in surfacewater flux with the largest simulated difference in one cell between baseline and chat emplacement scenarios being about 3.5 gallons per minute.

Scientific Investigations Report

Predictive analyses of ground-water discharges in the Willow Creek Watershed, northeast Nebraska

Ground-water discharge to Willow Creek, which drains a 204 mi 2 watershed in northeast Nebraska was predicted for six combinations of conditions of climate and ground-water development. A digital model of the ground water/surface-water system was stressed with recharge and withdrawal functions determined from a linear reservoir model of the soil zone. The geohydrologic system includes a regional aquifer which containssand and gravel ranging in age from Pliocene through Pleistocene. The regional aquifer is unconfined in the western part of the watershed and confined in the eastern part. The confining layer, or blue clay, in the east is of Pleistocene age and comprises principally saturated eolian silts with very fine sand interbeds overlying a basal lacustrine clay. Where unconfined conditions exist in the regional aquifer, perennial flow of Willow Creek is sustained principally by ground-water discharge. Where confined conditions exist, the low hydraulic conductivity of the blue clay effectively isolates the regional aquifer from Willow Creek. The stream-aquifer models were tested by comparing simulated against measured water levels and streamflow during 1975 and 1976. Agreement between simulated and observed values was obtained by modifying initial estimates of aquifer hydraulic conductivity and specific storage. Future ground-water discharge to Willow Creek (base flow) was simulated by super-imposing different patterns of ground-water withdrawals upon variations from average recharge for a monthly climatic sequence identical with the period 1931 to 1954. Frequency analyses of the simulated monthly base flows showed that they would be less than 12 feet 3 per second at least 50 percent of the time under all but completely undeveloped conditions. Simulated water-level declines in the regional aquifer at the end of the pumping season after a drought period similar to that of the 1930's were less than 5 feet where unconfined conditions exist and were more than 80 feet where confined conditions exist.

Nebraska

Effects of sanitary sewers on ground-water levels and streams in Nassau and Suffolk Counties, New York; part 2: Development and application of southwest Suffolk County model

By 1990, sanitary sewers in Nassau County Sewage Disposal Districts 2 and 3 and Suffolk County Southwest Sewer District will discharge to the ocean 140 cu ft of water/sec that would otherwise be returned to the groundwater system through septic tanks and similar systems. To evaluate the effects of this loss on groundwater levels and streamflow, the U.S. Geological Survey developed a groundwater flow model that couples a fine-scale subregional model to a regional model of larger scale. The regional model generates flux boundary conditions for the subregional model, and the subregional model provides detail in the area of concern. Results indicate that the water table will decline by as much as 8 ft along the Suffolk-Nassau county line, with effects decreasing eastward. Base flow is predicted to decrease by as much as 73% in a stream along the county line, but this effect will decrease to zero just east of the sewered area. This report is one of a series describing the predicted hydrologic effects of sewers in southern Nassau and southwest Suffolk Counties. (USGS)

New York

Geohydrology, geochemistry, and ground-water simulation-optimization of the Central and West Coast Basins, Los Angeles County, California

Historical ground-water development of the Central and West Coast Basins in Los Angeles County, California through the first half of the 20th century caused large water-level declines and induced seawater intrusion. Because of this, the basins were adjudicated and numerous ground-water management activities were implemented, including increased water spreading, construction of injection barriers, increased delivery of imported water, and increased use of reclaimed water. In order to improve the scientific basis for these water management activities, an extensive data collection program was undertaken, geohydrological and geochemical analyses were conducted, and ground-water flow simulation and optimization models were developed. In this project, extensive hydraulic, geologic, and chemical data were collected from new multiple-well monitoring sites. On the basis of these data and data compiled and collected from existing wells, the regional geohydrologic framework was characterized. For the purposes of modeling, the three-dimensional aquifer system was divided into four aquifer systems?the Recent, Lakewood, Upper San Pedro, and Lower San Pedro aquifer systems. Most pumpage in the two basins is from the Upper San Pedro aquifer system. Assessment of the three-dimensional geochemical data provides insight into the sources of recharge and the movement and age of ground water in the study area. Major-ion data indicate the chemical character of water containing less than 500 mg/L dissolved solids generally grades from calcium-bicarbonate/sulfate to sodium bicarbonate. Sodium-chloride water, high in dissolved solids, is present in wells near the coast. Stable isotopes of oxygen and hydrogen provide information on sources of recharge to the basin, including imported water and water originating in the San Fernando Valley, San Gabriel Valley, and the coastal plain and surrounding hills. Tritium and carbon-14 data provide information on relative ground-water ages. Water with abundant tritium (greater than 8 tritium units) is found in and downgradient from the Montebello Forebay and near the seawater barrier projects, indicating recent recharge. Water with less than measurable tritium is present in, and downgradient from, the Los Angeles Forebay and in most wells in the West Coast Basin. Water from several deep wells was analyzed for carbon-14. Uncorrected estimates of age for these samples range from 600 to more than 20,000 years before present. Chemical and isotopic data are combined to evaluate changes in chemical character along flow paths emanating from the Montebello and Los Angeles Forebays. A four-layer ground-water flow model was developed to simulate steady-state ground-water conditions representative of those in 1971 and transient conditions for the period 1971?2000. Model results indicate increases in ground-water storage in all parts of the study area over the simulated thirty-year period. The model was used to develop a three-dimensional ground-water budget and to assess impacts of two alternative future (2001?25) ground-water development scenarios?one that assumes continued pumping at average current rates and a second that assumes increasing pumping from most wells in the Central Basin. The model simulates stable or slightly increasing water levels for the first scenario and declining water levels (25 to 50 ft in the Central Basin) in the second scenario. Model sensitivity to parameter values and to the assumed Orange County boundary condition was evaluated. Particle tracking was applied to simulate advective transport of water from the spreading ponds, the coastline, and the seawater injection barriers. Particle tracking results indicate that most flow within the Upper San Pedro aquifer system occurs within about 20 percent of the total aquifer system thickness and that virtually all water injected into the seawater barrier projects has flowed inland. The simulation model was linked with optimizatio

California

Relation of nitrate concentrations in ground water to well depth, well use, and land use in Franklin Township, Gloucester County, New Jersey, 1970-85

A water-quality data base was developed to permit the investigation of the relation of concentrations of nitrate (as nitrogen) in ground water to well depth, well use, and land use (agricultural, residential, urban nonresidential, and undeveloped) in Franklin Township. Nitrate concentrations in water from 868 wells tended to decrease with depth. A rank-order regression model of nitrate concen- trations and land-use percentages was fitted to data from 98 shallow domestic wells. The model, which explains about 25 percent of the variance in the data, indicated that nitrate concentration increased with the percentage of developed land in a well's buffer zone. Further stratification of the data based on well use (commercial, domestic, or agricultural/irrigation) indicated that elevated nitrate concentrations were more common in water from agricultural/irrigation wells than in water from domestic or commercial wells. Concentrations of nitrate were indicative of human activities in water from about one-third of the wells sampled but exceeded the U.S. Environmental Protection Agency's maximum contaminant level of 10 milligrams per liter in water from only 1 percent of the wells. A sampling strategy in which water from wells of different depths located within areas in each of the four land-use categories is sampled yearly and analyzed for nitrate and other constituents would facilitate determination of the effects of human activities on ground-water quality.

New Jersey

WTAQ: A Computer Program for Calculating Drawdowns and Estimating Hydraulic Properties for Confined and Water-Table Aquifers

The computer program WTAQ calculates hydraulic-head drawdowns in a confined or water-table aquifer that result from pumping at a well of finite or infinitesimal diameter. The program is based on an analytical model of axial-symmetric ground-water flow in a homogeneous and anisotropic aquifer. The program allows for well-bore storage and well-bore skin at the pumped well and for delayed drawdown response at an observation well; by including these factors, it is possible to accurately evaluate the specific storage of a water-table aquifer from early-time drawdown data in observation wells and piezometers. For water-table aquifers, the program allows for either delayed or instantaneous drainage from the unsaturated zone. WTAQ calculates dimensionless or dimensional theoretical drawdowns that can be used with measured drawdowns at observation points to estimate the hydraulic properties of confined and water-table aquifers. Three sample problems illustrate use of WTAQ for estimating horizontal and vertical hydraulic conductivity, specific storage, and specific yield of a water-table aquifer by type-curve methods and by an automatic parameter-estimation method.

Water-Resources Investigations Report

A finite-element model for simulating hydraulic interchange of surface and ground water

A model was developed to be useful for predicting changes in streamflow as a result of groundwater pumping. The stream aquifer model is especially useful for simulating streams that flow intermittently owing to leakage to the aquifer or diversion for irrigation or streams that become perched owing to declining hydraulic head in the aquifer. The model couples the equation of two-dimensional groundwater flow with the kinematic equations of one-dimensional open-channel flow. Darcy 's law for vertical flow through a semipermeable streambed is used to couple the groundwater flow and streamflow equations. The equations of flow are approximated numerically by the finite-element method. A listing of the Fortran program that solves the equations of flow , and a description of data-input formats are given in the report. The program can simulate a variety of hydrologic characteristics including perched streams, streamflow diversions , springs, recharge from irrigated acreage, and evapotranspiration from the water table and phreatophytes. Time-dependent boundary conditions can be simulated. The program can be modified easily to simulate unconfined aquifers and aquifers with variable directions of anisotropy. (USGS)

Water-Resources Investigations Report

Simulation analysis of the unconfined aquifer, Raft River geothermal area, Idaho-Utah

This study covers about 1,000 mi 2 (2,600 km 2 ) of the southern Raft River drainage basin in south-central Idaho and northwest Utah. The main area of interest, approximately 200 mi2 (520 km 2 ) of semiarid agricultural and rangeland in the southern Raft River Valley that includes the known Geothermal Resource Area near Bridge, Idaho, was modelled numerically to evaluate the hydrodynamics of the unconfined aquifer. Computed and estimated transmissivity values range from 1,200 feet squared per day (110 meters squared per day) to 73,500 feet squared per day (6,830 meters squared per day). Water budgets, including ground-water recharge and discharge for approximate equilibrium conditions, have been computed by several previous investigators; their estimates of available ground-water recharge range from about 46,000 acre-feet per year (57 cubic hectometers per year) to 100,000 acre-feet per year (123 cubic hectometers per year). Simulation modeling of equilibrium conditions represented by 1952 water levels suggests: (1) recharge to the water-table aquifer is about 63,000 acre-feet per year (77 cubic hectometers per year); (2) a significant volume of ground water is discharged through evapotranspiration by phreatophytes growing on the valley bottomlands; (3) the major source of recharge may be from upward leakage of water from a deeper, confined reservoir; and (4) the aquifer transmissivity probably does not exceed about 12,000 feet squared per day (3,100 meters squared per day). Additional analysis carried out by simulating transient conditions from 1952 to 1965 strongly suggests that aquifer transmissivity does not exceed about 7,700 feet squared per day (700 meters squared per day). The model was calibrated using slightly modified published pumpage data; it satisfactorily reproduced the historic water-level decline over the period 1952-65.

Idaho, Utah

A model for flow through a glacial outwash aquifer in southeast Franklin County, Ohio

A glacial aquifer of 70 square miles in the Scioto River valley southeast of Columbus, Ohio, was modeled as a potentially major source of water. The model was constructed from available hydrologic data: Records of precipitation, well hydrographs, well logs, two ground-water level surveys, and analyses of six aquifer tests. Utilizing this array of data, water levels determined from a series of steady-state simulations of different hydraulic conductivity distributions were calibrated against measured (December 1977) ground-water levels. The simulations that provided the best matches used two hydraulic conductivity distributions: One was an areally varying hydraulic conductivity distribution; the other an areally uniform hydraulic conductivity (40 feet per day) distribution. After these more probable hydraulic conductivity distributions were found, they were utilized in steady state maximal pumping simulations. The maximal well-field yield for these simulations was 20.5 million gallons per day for the areally varying hydraulic conductivity distribution, and 11.3 million gallons per day for the areally uniform hydraulic conductivity. Sensitivity of well yield to changes in well position and streambed leakance changes was investigated also.

Ohio

Geohydrology of Storage Unit III and a combined flow model of the Santa Barbara and foothill ground-water basins, Santa Barbara County, California

The city of Santa Barbara pumps most of its ground water from the Santa Barbara and Foothill ground-water basins. The Santa Barbara basin is subdivided into two storage units: Storage Unit I and Storage Unit III. The Foothill basin and Storage Unit I of the Santa Barbara basin have been studied extensively and ground-water flow models have been developed for them. In this report, the geohydrology of the Santa Barbara ground- water basin is described with a special emphasis on Storage Unit III in the southwestern part of the basin. The purposes of this study were to summarize and evaluate the geohydrology of Storage Unit III and to develop an areawide model of the Santa Barbara and Foothill basins that includes the previously unmodeled Storage Unit III. Storage Unit III is in the southwestern part of the city of Santa Barbara. It is approximately 3.5 miles long and varies in width from about 2,000 feet in the southeast to 4,000 feet in the north-west. Storage Unit III is composed of the Santa Barbara Formation and overlying alluvium. The Santa Barbara Formation (the principal aquifer) consists of Pleistocene and Pliocene(?) unconsolidated marine sand, silt, and clay, and it has a maximum saturated thickness of about 160 feet. The alluvium that overlies the Santa Barbara Formation has a maximum saturated thickness of about 140 feet. The storage unit is bounded areally by faults and low-permeability deposits and is underlain by rocks of Tertiary age. The main sources of recharge to Storage Unit III are seepage from Arroyo Burro and infiltration of precipitation. Most of the recharge occurs in the northwest part of the storage unit, and ground water flows toward the southeast along the unit's long axis. Lesser amounts of recharge may occur as subsurface flow from the Hope Ranch subbasin and as upwelling from the underlying Tertiary rocks. Discharge from Storage Unit III occurs as pumpage, flow to underground drains, underflow through alluvium in the vicinity of Arroyo Burro across the Lavigia Fault, evapotranspiration, and underflow to the Pacific Ocean. The faults that bound Storage Unit III generally are considered to be effective barriers to the flow of ground water. Interbasin ground-water flow occurs where deposits of younger alluvium along stream channels cross faults. Ground-water quality in Storage Unit III deposits varies with location and depth. Upward leakage of poor-quality water from the underlying Tertiary rocks occurs in the storage unit, and such leakage can be influenced by poor well construction or by heavy localized pumping. The highest dissolved-solids concentration (4,710 milligrams per liter) in ground water resulting from this upward leakage is found in the coastal part of the storage unit. The ground-water system was modeled as two horizontal layers. In the Foothill basin and Storage Unit I the layers are separated by a confining bed. The upper layer represents the upper producing zone and the shallow zone near the coast. The lower layer represents the lower producing zone. In general, the faults in the study area were assumed to be no-flow boundaries, except for the offshore fault that forms the southeast boundary; the southeast boundary was simulated as a general-head boundary. The Storage Unit III model was combined with the preexisting Storage Unit I and Foothill basin models, using horizontal flow barriers, to form an areawide model. The areawide model was calibrated by simulating steady-state predevelopment conditions and transient conditions for 1978-92. The nonpumping steady- state simulation was used to verify that the calibrated model yielded physically reasonable results for predevelopment conditions. The calibrated areawide model calculates water levels in Storage Unit III that are within 10 feet of measured water levels at all sites of comparison. In addition, the model adequately simulates water levels in the Storage Unit I and Foothill basin areas. A total of 33,430 acre-feet of water was pumped from the study area during the simulation period. Model results indicate that 2,833 acre-feet came from storage and 5, 332 acre-feet crossed the general-head boundary from the ocean, thus indicating that seawater intrusion could occur. A sensitivity analysis indicates that, in general, the model is most sensitive to changes in transmissivity and total recharge.

California