Geology Reports⌕ Search

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 1,045 records · Page 58Linked to original sources
USGS2005

Use of a ground-water flow model to delineate contributing areas to the Puchack Well Field, Pennsauken township and vicinity, Camden county, New Jersey

The New Jersey Department of Environmental Protection (NJDEP) Well Head Protection Program, developed in response to the 1986 Federal Safe Drinking Water Act Amendments, requires delineation of Well Head Protection Areas (WHPA's), commonly called contributing areas, for all public and non-community water-supply wells in New Jersey. Typically, WHPA's for public community water-supply wells in New Jersey are delineated using a two-dimensional ground-water flow model incorporating the regional hydraulic gradient; however, NJDEP guidelines allow for the use of a three-dimensional flow model to delineate contributing areas to wells in complex hydrogeologic settings. The Puchack well field in Pennsauken Township, Camden County, N.J., is an area of strong hydraulic connection between the Lower aquifer of the Potomac-Raritan-Magothy aquifer system and the Delaware River. Interactions among and within the public-supply well fields in the area are complex. To delineate the contributing area to the Puchack well field, the U.S. Geological Survey, in cooperation with the NJDEP, developed an 11-layer ground-water flow model of the Potomac-Raritan-Magothy aquifer system in the Pennsauken Township area to simulate flow in the vicinity of the well field. The model incorporates the interaction between the aquifer system and the Delaware River, and includes boundary flows from an existing regional model of the Camden area. Recharge used in the model ranged from 4.5 to 14 inches per year, and horizontal hydraulic conductivity ranged from 50 to 250 feet per day. Values of vertical hydraulic conductivity ranging from 0.001 to 0.5 feet per day were assigned to zones created on the basis of variations in hydrogeologic conditions observed in geophysical logs from wells. A steady-state simulation was used to calibrate the model to synoptic water-level data collected in March 1998. Near the Puchack well field, simulated heads generally were within 1 foot of the measured heads in both the Middle and Lower aquifers. Simulated water-level differences across the confining units at most of the nested wells were within ? 0.5 feet of the differences calculated from measured water levels. The existing flow model was modified to meet NJDEP guidelines for delineating contributing areas in complex hydrogeologic settings. These modifications included rediscretizing the model grid to a finer grid and preparing the water-use data set for use in the rediscretized model. The contributing area to the Puchack well field was delineated by means of particle tracking. An uncertainty analysis was conducted in which 36 model-input parameters were both increased and decreased until the resulting change in simulated heads exceeded the model-calibration criterion of ? 5 feet at any model cell. Porosity most affected the size and shape of the contributing area. The distribution of withdrawals at the Morris/Delair well field and variations in recharge affected both the size and shape of contributing area to the Puchack well field and the source of water to the Puchack wells. The results of the uncertainty analysis were combined to determine the 'aggregate' contributing area to the Puchack well field--a composite of areas on the land surface that contributed flow to the Puchack well field in less than 12 years in any uncertainty simulation. The shape of the aggregate contributing area was most similar to that associated with a reduction in porosity, which indirectly affected the size and shape of the contributing areas by changing travel time.

Scientific Investigations Report↗
USGS1980

Hydrology and model study of the proposed Prosperity Reservoir, Center Creek Basin, southwestern Missouri

A dam and reservoir have been proposed for construction on Center Creek, Jasper County, in southwestern Missouri. Ground-water levels in the hills adjacent to the reservoir will rise when the impoundment is completed. One of the problems is that the proposed site of Prosperity Reservoir is a few miles upstream from the lead-zinc mining area known as the Oronogo-Duenweg belt. In this belt transmissivities are variable but appear to be higher than they are in the immediate area of the reservoir. Grove Creek lies down-gradient from the reservoir area and separates it from the mining belt. A model study indicates that inflow from the proposed reservoir to the water table could cause water level rises varying from about 20 feet near the reservoir to 0.5 to 1.0 foot in the southern part of Grove Creek drainage basin. These rises will cause significant changes to the natural ground-water flow system. Increased ground-water elevations in the reservoir area could result in increased ground-water gradients and discharge to Grove and Center Creeks. The increase in ground-water discharge to Grove Creek, and in turn Center Creek, will have the beneficial effect of diluting mine-water discharge from the Oronogo-Duenweg belt during periods of low flow. However, if Grove Creek does not act as an effective drain and if conduits extend beneath Grove Creek to transfer the increased water available to the Oronogo-Duenweg belt, the flow regimen could change in the mining belt west of Grove Creek increasing mine-water discharge to Center Creek downstream from the reservoir. Bedrock in the area is Mississippian limestone, the deeply solutioned formation that contained the ore deposits. The limestone in the mining district was greatly altered by solution prior to ore deposition while the limestone in the area of the reservoir was altered less. The extent of the alteration is related to the aquifer characteristics in that high and low values of transmissivity and storage coefficient correspond to greatly altered brecciated rocks in the mining district and less altered, less brecciated rocks in the reservoir area, respectively. The authors suggest that an ancestral east-flowing White River drained the area about Joplin in Late Mississippian time. This is based on the configuration of the contact between Meramecian and Osagean rocks of Mississippian age. A high topographic area existed in the region about Joplin in which the water table stood 200 feet below the land surface when sinkholes and caverns of that depth were formed. The large number of Pennsylvanian-filled sinkholes in the Joplin area and the smaller number to the east suggest a higher land surface to the west than that to the east. The distribution of paleokarst sinkholes supports the conclusion based on the configuration of the Meramecian-Osagean contact.

Missouri↗
USGS2015

Geologic map of the Vashon 7.5' quadrangle and selected areas, King County, Washington

This map is an interpretation of a 6-ft-resolution lidar-derived digital elevation model combined with geology by Derek B. Booth and Kathy Goetz Troost. Field work by Booth and Troost was located on the 1:24,000-scale topographic map of the Vashon and Des Moines 7.5' quadrangles that were published in 1997 and 1995, respectively. Much of the geology was interpreted from landforms portrayed on the topographic maps, supplemented by field exposures, where available. In 2001, the Puget Sound Lidar Consortium (see http://pugetsoundlidar.org/) obtained a lidar-derived digital elevation model (DEM) for Vashon Island and the Des Moines quadrangle. For a brief description of lidar and this data acquisition program, see Haugerud and others (2003). This new DEM has a horizontal resolution of 6 ft (1.83 m) and mean vertical accuracy of about 1 ft (about 0.3 m). The greater resolution and accuracy of the lidar DEM facilitated a much-improved interpretation of many aspects of the surficial geology, especially the distribution and relative age of landforms and the materials inferred to comprise them. Booth and Troost were joined by Tabor to interpret the new lidar DEM but have done no futher field work for this map. This map, the Vashon quadrangle and selected adjacent areas, encompasses most of Vashon Island, Maury Island, and Three Tree Point in the south-central Puget Sound. One small area in the Vashon quadrangle on the east side of Puget Sound is excluded from this map but included on the adjacent Seattle quadrangle (Booth and others, 2005). The map displays a wide variety of surficial geologic deposits, which reflect many geologic environments and processes. Multiple ice-sheet glaciations and intervening nonglacial intervals have constructed a complexly layered sequence of deposits that underlie both islands to a depth of more than 300 m below sea level. These deposits not only record glacial and nonglacial history but also control the flow and availability of ground water, determine the susceptibility of the slopes to landslides, and provide economic reserves of sand and gravel. The islands are surrounded by channels of Puget Sound, some as deep as the islands are high (>600 ft (~200 m)). The shorelines provide many kilometers of well-exposed coastal outcrops that reveal abundant lithologic and stratigraphic details not ordinarily displayed in the heavily vegetated Puget Lowland.

Washington↗
USGS1987

Effects of proposed highway embankment modifications on water-surface elevations in the lower Pearl River flood plain near Slidell, Louisiana

Major flooding in the lower Pearl River basin in recent years has caused extensive damage to homes and highways in the area. In 1980 and 1983, Interstate Highway 10 and U.S. Highway 190 were overtopped. In 1983, the Interstate Highway 10 crossing was seriously damaged by the flood. The U.S. Geological Survey, in cooperation with the Louisiana Department of Transportation and Development, Office of Highways, used a two-dimensional finite-element surface-water flow model to evaluate the effects the proposed embankment modifications at Interstate Highway 10 and U.S. Highway 90 on the water-surface elevations in the lower Pearl River flood plain near Slidell, Louisiana. The proposed modifications that were considered for the 1983 flood are: (1) Removal of all highway embankments, the natural condition, (2) extension of the West Pearl River bridge by 1,000 feet at U.S. Highway 90, (3) construction of a new 250-foot bridge opening in the U.S. Highways 190 and 90, west of the intersection of the highways. The proposed highway bridge modifications also incorporated lowering of ground-surface elevations under the new bridges to sea level. The modification that provided the largest reduction in backwater, about 35 percent, was a new bridge in Interstate Highway 10. The modification of the West Pearl River bridge at U.S. Highway 90 and replacement of the bridge in U.S. Highway 190 provide about a 25% reduction in backwater each. For the other modification conditions that required structural modifications, maximum backwater computed on the west side of the flood plain ranges from 0.0 to 0.8 foot and on the east side from 0.0 to 0.6 foot. Results show that although backwater is greater on the west side of the flood plain than on the east side, upstream of highway embankments, backwater decreases more rapidly in the upstream direction on the west side of the flood plain than on the east side. Analysis of the proposed modifications indicates that backwater would still occur on the east and west sides of the flood plain, but values would be less than those computed with highway embankments in place. (Author 's abstract)

Water-Resources Investigations Report↗
USGS1988

Review of mechanisms, methods, and theory for determining recharge to shallow aquifers in North Dakota

Effective management of ground-water resources requires knowledge of all components of the water budget for the aquifer of interest. Efforts to simulate ground-water flow prior to development and the effects of proposed pumping in several of North Dakota's shallow glacial aquifers have been hindered by the lack of reliable estimates of ground-water recharge. This study was done to (1) review the methods that have been used to measure recharge, (2) review the theory of unsaturated flow and the methods for characterizing the physical properties of unsaturated media, (3) consider the relative merits of a rigorous data-intensive approach versus an estimation approach to the study of recharge, and (4) review past and current agronomic research in North Dakota for applicability of the research and the data generated to the study of recharge. Direct, quantitative techniques for evaluating recharge are rarely applied. The theory for computing fluxes in unsaturated media is well established and numerous physics-based models that effectively implement the theory are available, but the data required for the models generally are lacking. Many parametric approaches have been developed to avoid the large data requirements of the physics-based approaches for analyzing flow in the unsaturated zone. However, the parametric approaches normally include fitting coefficients that must be calibrated for every study site, thereby detracting from the general utility of the parametric approach. The functional relation of matric potential to moisture content is required for physics-based soil-water models, whether analytic or numeric. Laboratory methods to determine these relations are tedious, costly, and may not give results representative of the soils as they occur in the field. Many models have been proposed to estimate the moisture-characteristic curve and hydraulic-conductivity function from basic soil properties, but none yield results that are universally satisfactory. In situ methods, because they require minimal disturbance of the soil profile and may be used repeatedly on the same soil mass, have become the preferred means for acquiring physical data, especially hydraulic conductivity. Hydro logic investigations, except for recent studies of hazardous-waste disposal sites, rarely have included physical characterizations of unsaturated media. Any of four phenomena could hinder attempts to simulate unsaturated flow in settings typical of North Dakota; variability of soil properties, hysteresis, frozen ground, and macropore development. The spatial and temporal variability of soil properties probably is the greatest complicating phenomenon and must be dealt with by detailed characterization of the properties. Hysteresis can detract from the accuracy of flow calculations for some soils under certain conditions but, for the present, our scant knowledge of soil physical properties is a greater hindrance to reliable soi1-water mode 1 ing than is the hysteresis phenomenon. A1 though seasona1ly frozen ground undoubtedly affects hydrologic processes in North Dakota, much more research is needed before meaningful quantitative treatment is possible. Finally, macropores can influence soil-water movement significantly, but macropore development may not be common on the intensively farmed, coarse-textured soils that typically overlie North Dakota's glacial aquifers. Lysimetry currently is the only reliable means of analyzing macropore flow. The soil-related research that has been conducted in North Dakota to date (1983) provides little of the type of information required to estimate ground-water recharge. Useful data could be developed by systematically evaluating the hydraulic characteristics of the prominent soil types overlying North Dakota's shallow glacial aquifers. These data would be required to enable use of a physics-based approach to estimating recharge. The size of the aquifer under study, its economic value, and the resources available for data collection should be considered when choosing between parametric or physics-based methods.

Water-Resources Investigations Report↗
USGS1996

Availability and quality of water from drift aquifers in Marshall, Pennington, Polk, and Red Lake counties, northwestern Minnesota

Sand and gravel aquifers present within glacial deposits are important sources of water in Marshall, Pennington, Polk, and Red Lake Counties in northwestern Minnesota. Saturated thicknesses of the unconfined aquifers range from 0 to 30 feet. Estimated horizontal hydraulic conductivities range from 2.5 to 600 feet per day. Transmissivity of the unconfined aquifers ranges from 33 to greater than 3,910 feet squared per day. Theoretical maximum well yields for 6 wells with specific-capacity data range from 12 to 123 gallons per minute. Saturated thicknesses of shallow confined aquifers (depth to top of the aquifer less than 100 feet below land surface) range from 0 to 150 feet. Thicknesses of intermediate, deep, and basal confined aquifers (depths to top of the aquifer from 100 to 199 feet, from 200 to 299 feet, and 300 feet or more below land surface, respectively) range from 0 to more than 126 feet. Transmissivity of the confined aquifers ranges from 2 to greater than 210,000 feet squared per day. Theoretical maximum well yields range from 3 to about 2,000 gallons per minute. Recharge to ground water is predominantly from precipitation that percolates downward to the saturated zone. Recharge to unconfined aquifers in the study area ranged from 4.5 to 12.0 inches per year during 1991 and 1992, based on hydrograph analysis. Model simulations done for this study indicate that recharge rates from 8 to 9 inches per year to unconfined aquifers produce the best matches between model-simulated and measured water levels in wells. Discharge from ground water occurs by seepage to streams, lakes and wetlands, ground-water evapotranspiration, and withdrawals through wells. In 1990, total ground-water withdrawals in the study area were 6.0 million gallons per day. All of the withdrawals were from drift aquifers. Numerical models of ground-water flow were constructed to represent two beach-ridge aquifer systems under steady-state conditions. Beach-ridge aquifer systems were simulated in Pennington, Polk, and Red Lake County. Simulated recharge from the infiltration of precipitation accounts for most of the sources of water to the beach-ridge aquifer systems and simulated evapotranspiration accounts for all of the discharge other than ground-water withdrawals. The numerical-model simulations indicate that upward movement of water from underlying confined aquifers to overlying unconfined aquifers is an important component of ground-water flow within the beach-ridge aquifer systems. Simulated long-term, steady-state yields from the unconfined aquifers are generally less than 50 gallons per minute, due to the generally low saturated thickness of the aquifers and the relatively low hydraulic conductivity of the aquifer material. Water from all the drift aquifers in the study area is very hard (more than 180 milligrams per liter of calcium carbonate). The predominant ions in water from the unconfined and shallow confined aquifers were generally calcium and bicarbonate. Water from the intermediate confined aquifers includes a variety of water types, including calcium bicarbonate, calcium sulfate, mixed calcium-sodium bicarbonate, and sodium chloride type waters. Waters from the deep confined aquifers are predominantly calcium bicarbonate, mixed calcium-sodium bicarbonate, and sodium chloride type waters. Mean concentrations of calcium and magnesium generally decreased with depth below land surface. Mean concentrations of sodium and sulfate generally increased with depth. Mean chloride concentrations were greatest for the shallow and deep confined aquifers and least for the unconfined and intermediate confined aquifers. The concentration and percentage (as percent of total cations) of sodium, and concentration of dissolved solids tend to increase from east to west along regional flow paths. Concentrations and percentages (as percent of total anions) of chloride tend to be greater in the western part of the study area than in the eastern part. These trends are probably due to longer residence time of the water in the flow system, and upward leakage of water from the underlying Cretaceous and Paleozoic strata. Waters from the drift aquifers underlying most of the study area generally are suitable for domestic consumption, crop irrigation, and most other uses. Water from 20 wells screened in unconfined and confined aquifers exceeded U.S. Environmental Protection Agency recommended limits for dissolved solids concentrations. Chemical analyses of waters from the unconfined and confined aquifers generally indicated a potentially low sodium hazard and a medium to high salinity hazard for irrigation. Water samples analyzed for nitrate had nitrate concentrations below the reporting limit (0.05 milligrams per liter) in 10 out of 23 wells. Two samples had nitrate concentrations greater than 10 milligrams per liter. Pesticide concentrations in water samples from 17 wells screened in unconfined and shallow confined aquifers were below or only slightly above laboratory reporting limits.

Minnesota↗
USGS2002

Simulation of ground-water flow and evaluation of water-management alternatives in the upper Charles River basin, eastern Massachusetts

Ground water is the primary source of drinking water for towns in the upper Charles River Basin, an area of 105 square miles in eastern Massachusetts that is undergoing rapid growth. The stratified-glacial aquifers in the basin are high yield, but also are thin, discontinuous, and in close hydraulic connection with streams, ponds, and wetlands. Water withdrawals averaged 10.1 million gallons per day in 1989?98 and are likely to increase in response to rapid growth. These withdrawals deplete streamflow and lower pond levels. A study was conducted to develop tools for evaluating water-management alternatives at the regional scale in the basin. Geologic and hydrologic data were compiled and collected to characterize the ground- and surface-water systems. Numerical flow modeling techniques were applied to evaluate the effects of increased withdrawals and altered recharge on ground-water levels, pond levels, and stream base flow. Simulation-optimization methods also were applied to test their efficacy for management of multiple water-supply and water-resource needs. Steady-state and transient ground-water-flow models were developed using the numerical modeling code MODFLOW-2000. The models were calibrated to 1989?98 average annual conditions of water withdrawals, water levels, and stream base flow. Model recharge rates were varied spatially, by land use, surficial geology, and septic-tank return flow. Recharge was changed during model calibration by means of parameter-estimation techniques to better match the estimated average annual base flow; area-weighted rates averaged 22.5 inches per year for the basin. Water withdrawals accounted for about 7 percent of total simulated flows through the stream-aquifer system and were about equal in magnitude to model-calculated rates of ground-water evapotranspiration from wetlands and ponds in aquifer areas. Water withdrawals as percentages of total flow varied spatially and temporally within an average year; maximum values were 12 to 13 percent of total annual flow in some subbasins and of total monthly flow throughout the basin in summer and early fall. Water-management alternatives were evaluated by simulating hypothetical scenarios of increased withdrawals and altered recharge for average 1989?98 conditions with the flow models. Increased withdrawals to maximum State-permitted levels would result in withdrawals of about 15 million gallons per day, or about 50 percent more than current withdrawals. Model-calculated effects of these increased withdrawals included reductions in stream base flow that were greatest (as a percentage of total flow) in late summer and early fall. These reductions ranged from less than 5 percent to more than 60 percent of model-calculated 1989?98 base flow along reaches of the Charles River and major tributaries during low-flow periods. Reductions in base flow generally were comparable to upstream increases in withdrawals, but were slightly less than upstream withdrawals in areas where septic-system return flow was simulated. Increased withdrawals also increased the proportion of wastewater in the Charles River downstream of treatment facilities. The wastewater component increased downstream from a treatment facility in Milford from 80 percent of September base flow under 1989?98 conditions to 90 percent of base flow, and from 18 to 27 percent of September base flow downstream of a treatment facility in Medway. In another set of hypothetical scenarios, additional recharge equal to the transfer of water out of a typical subbasin by sewers was found to increase model-calculated base flows by about 12 percent of model-calculated base flows. Addition of recharge equal to that available from artificial recharge of residential rooftop runoff had smaller effects, augmenting simulated September base flow by about 3 percent. Simulation-optimization methods were applied to an area near Populatic Pond and the confluence of the Mill and Charles Rivers in Franklin,

Massachusetts↗
USGS2000

Documentation of a computer program to simulate lake-aquifer interaction using the MODFLOW ground water flow model and the MOC3D solute-transport model

Heads and flow patterns in surficial aquifers can be strongly influenced by the presence of stationary surface-water bodies (lakes) that are in direct contact, vertically and laterally, with the aquifer. Conversely, lake stages can be significantly affected by the volume of water that seeps through the lakebed that separates the lake from the aquifer. For these reasons, a set of computer subroutines called the Lake Package (LAK3) was developed to represent lake/aquifer interaction in numerical simulations using the U.S. Geological Survey three-dimensional, finite-difference, modular ground-water flow model MODFLOW and the U.S. Geological Survey three-dimensional method-of-characteristics solute-transport model MOC3D. In the Lake Package described in this report, a lake is represented as a volume of space within the model grid which consists of inactive cells extending downward from the upper surface of the grid. Active model grid cells bordering this space, representing the adjacent aquifer, exchange water with the lake at a rate determined by the relative heads and by conductances that are based on grid cell dimensions, hydraulic conductivities of the aquifer material, and user-specified leakance distributions that represent the resistance to flow through the material of the lakebed. Parts of the lake may become "dry" as upper layers of the model are dewatered, with a concomitant reduction in lake surface area, and may subsequently rewet when aquifer heads rise. An empirical approximation has been encoded to simulate the rewetting of a lake that becomes completely dry. The variations of lake stages are determined by independent water budgets computed for each lake in the model grid. This lake budget process makes the package a simulator of the response of lake stage to hydraulic stresses applied to the aquifer. Implementation of a lake water budget requires input of parameters including those representing the rate of lake atmospheric recharge and evaporation, overland runoff, and the rate of any direct withdrawal from, or augmentation of, the lake volume. The lake/aquifer interaction may be simulated in both transient and steady-state flow conditions, and the user may specify that lake stages be computed explicitly, semi-implicitly, or fully-implicitly in transient simulations. The lakes, and all sources of water entering the lakes, may have solute concentrations associated with them for use in solute-transport simulations using MOC3D. The Stream Package of MODFLOW-2000 and MOC3D represents stream connections to lakes, either as inflows or outflows. Because lakes with irregular bathymetry can exist as separate pools of water at lower stages, that coalesce to become a single body of water at higher stages, logic was added to the Lake Package to allow the representation of this process as a user option. If this option is selected, a system of linked pools (sublakes) is identified in each time step and stages are equalized based on current relative sublake surface areas.

Water-Resources Investigations Report↗
USGS1999

Episodic sediment-discharge events in Cascade Springs, southern Black Hills, South Dakota

Cascade Springs is a group of artesian springs in the southern Black Hills, South Dakota, with collective flow of about 19.6 cubic feet per second. Beginning on February 28, 1992, a large discharge of red suspended sediment was observed from two of the six known discharge points. Similar events during 1906-07 and 1969 were documented by local residents and newspaper accounts. Mineralogic and grain-size analyses were performed to identify probable subsurface sources of the sediment. Geochemical modeling was performed to evaluate the geochemical evolution of water discharged from Cascade Springs. Interpretations of results provide a perspective on the role of artesian springs in the regional geohydrologic framework. X-ray diffraction mineralogic analyses of the clay fraction of the suspended sediment were compared to analyses of clay-fraction samples taken from nine geologic units at and stratigraphically below the spring-discharge points. Ongoing development of a subsurface breccia pipe(s) in the upper Minnelusa Formation and/or Opeche Shale was identified as a likely source of the suspended sediment; thus, exposed breccia pipes in lower Hell Canyon were examined. Upper Minnelusa Formation breccia pipes in lower Hell Canyon occur in clusters similar to the discrete discharge points of Cascade Springs. Grain-size analyses showed that breccia masses lack clay fractions and have coarser distributions than the wall rocks, which indicates that the red, fine-grained fractions have been carried out as suspended sediment. These findings support the hypothesis that many breccia pipes were formed as throats of abandoned artesian springs. Geochemical modeling was used to test whether geochemical evolution of ground water is consistent with this hypothesis. The evolution of water at Cascade Springs could not be suitably simulated using only upgradient water from the Minnelusa aquifer. A suitable model involved dissolution of anhydrite accompanied by dedolomitization in the upper Minnelusa Formation, which is caused by upward leakage of relatively fresh water from the Madison aquifer. The anhydrite dissolution and dedolomitization account for the net removal of minerals that would lead to breccia pipe formation by gravitational collapse. Breccia pipes in the lower Minnelusa Formation are uncommon; however, networks of interconnected breccia layers and breccia dikes are common. These networks, along with vertical fractures and faults, are likely pathways for transmitting upward leakage from the Madison aquifer. It is concluded that suspended sediment discharged at Cascade Springs probably results from episodic collapse brecciation that is caused by subsurface dissolution of anhydrite beds and cements of the upper Minnelusa Formation, accompanied by replacement of dolomite by calcite. It is further concluded that many breccia pipes probably are the throats of artesian springs that have been abandoned and exposed by erosion. The locations of artesian spring-discharge points probably have been shifting outwards from the center of the Black Hills uplift, essentially keeping pace with regional erosion over geologic time. Thus, artesian springflow probably is a factor in controlling water levels in the Madison and Minnelusa aquifers, with hydraulic head declining over geologic time, in response to development of new discharge points. Development of breccia pipes as throats of artesian springs would greatly enhance vertical hydraulic conductivity in the immediate vicinity of spring-discharge points. Horizontal hydraulic conductivity in the Minnelusa Formation also may be enhanced by dissolution processes related to upward leakage from the Madison aquifer. Potential processes could include dissolution resulting from leakage in the vicinity of breccia pipes that are abandoned spring throats, active spring discharge, development of subsurface breccias with no visible surface expression or spring discharge, as well as general areal leakage from the Madison aquifer into the Minnelusa Formation.

South Dakota↗
USGS2008

Recovery of Ground-Water Levels from 1988 to 2003 and Analysis of Effects of 2003 and Full-Allocation Withdrawals in Critical Area 2, Southern New Jersey

Water levels in the Potomac-Raritan-Magothy aquifer system within Water Supply Critical Area 2 in the southern New Jersey Coastal Plain have recovered as a result of reductions in ground-water withdrawals initiated in the early 1990s. The Critical Area consists of the depleted zone and the threatened margin. The Potomac-Raritan-Magothy aquifer system consists of the Upper, Middle, and Lower aquifers. Generally, ground-water withdrawals from these aquifers declined 5 to 10 Mgal/d (million gallons per day) and water levels recovered 0 to 40 ft (foot) from 1988 to 2003. In order to reevaluate water-allocation restrictions in Critical Area 2 in response to changes in the ground-water-flow system and demands for additional water supply due to increased development, the New Jersey Department of Environmental Protection (NJDEP) needs information about the effects of changes in those allocations. Therefore, the U.S. Geological Survey (USGS), in cooperation with the NJDEP, used an existing ground-water-flow model of the New Jersey Coastal Plain to evaluate the effects of withdrawal alternatives on hydraulic heads in the Potomac-Raritan-Magothy aquifer system in Critical Area 2. The U.S. Geological Survey Regional Aquifer System Analysis model was used to simulate steady-state ground-water flow. Two withdrawal conditions were tested by using the model to evaluate hydraulic heads and differences in heads in these aquifers: 2003 withdrawals and full-allocation withdrawals (17.4 Mgal/d greater than 2003 withdrawals). Model results are presented using head maps and head-difference maps that compare 2003 to full-allocation withdrawals. Mandated hydrologic conditions for Critical Area protection are that the simulated -30-ft head contour not extend beyond the boundary of the depleted zone and (or) be at least 5 mi (miles) updip from the 250-mg/L (milligram per liter) isochlor in all three aquifers. Simulation results indicate that, for 2003 withdrawals, the simulated -30-ft head contour in all three aquifers is generally within the boundary of the depleted zone, except in the Lower aquifer in northern Camden and northwestern Burlington Counties, and is generally 1 to 10 mi downdip from the 250-mg/L isochlor. (Corresponding observed data indicate that the -30-ft water-level contour extends beyond the southwest boundary of the depleted zone in the Upper and Middle aquifers, and is generally 5 to 20 mi downdip from the 250-mg/L isochlor in all three aquifers.) The area in which heads are below -30 ft ranges from 389 mi2 (square miles) in the Middle aquifer to 427 mi2 in the Lower aquifer. For full-allocation withdrawals, the simulated -30-ft head contour extends beyond the boundary of the depleted zone in all three aquifers in northern Camden and northwestern Burlington Counties and in the Upper aquifer in Gloucester and Salem Counties, and is generally 5 to 15 mi downdip from the 250-mg/L isochlor. The area in which heads are below -30 ft ranges from 616 mi2 in the Upper aquifer to 813 mi2 in the Lower aquifer. These results and observed data indicate that any increase in withdrawals from 2003 values would likely cause heads in the three aquifers to decline below the minimum values mandated by the NJDEP for the Critical Area.

Scientific Investigations Report↗
USGS1999

Evaluation of factors that influence estimated zones of transport for six municipal wells in Clark County, Washington

A ground-water flow model was used in conjunction with particle tracking to estimate zones of transport for six municipal well sites in Clark County, Washington. A zone of transport for a well is a three-dimensional volume within a ground-water system that contains all of the ground water that will discharge from that well within a specified time period. All of the zones of transport for a well compose the zone of contribution for the well. Zones of transport and contribution are important considerations in the delineation of wellhead-protection areas. Hydrogeologic factors, such as hydraulic conductivity and porosity, influence the shape and size of the zones of transport, and, therefore, uncertainty in these and other factors can lead to uncertainty in the delineation of the zones of transport. The sensitivity of the zones of transport to uncertainty in selected hydrogeologic factors was evaluated for the six wells. Estimates of the zones of transport were delineated by the U.S. Geological Survey program MODTOOLS from three-dimensional pathlines computed by the U.S. Geological Survey program MODPATH. Input to MODPATH came from steady-state simulations calculated by the U.S. Geological Survey modular three-dimensional finite-difference ground-water flow model, MODFLOW. Three-dimensional modeling is the best method for delineating zones of transport within stratigraphically complex, heterogeneous, anisotropic aquifers that have complex boundary conditions such as streams and multiple, simultaneously discharging wells. In this study, zones of transport were delineated by using simulated particle locations computed from the results of a three-dimensional steady-state regional model for 0-0.5, 0.5-1, 1-5, 5-10, 10-20, and 20-50 year travel times to the selected wells. Zones of transport for a well were delineated by tracking particles along pathlines in the reverse direction of ground-water flow. Sensitivity of the zones of transport to change in the discharge rate of the selected well, porosity, and hydraulic conductivity, as well as to the presence or absence of interfering wells, was evaluated at six well sites to evaluate the effect of uncertainties in these factors on the size and shape of zones of transport. Uncertainty in porosity contributed the most to the uncertainty in delineating the zones of transport. Uncertainty in other factors, such as well discharge rate and horizontal hydraulic conductivity, had measurable effects on the zones of transport, but errors introduced through these factors were less significant. Insight into the causes of the changes in the size and shape of the zones of transport to varying conditions was gained by evaluating the simulated water budget and ground-water levels in the vicinity of the well. Changes in the simulated water budget and ground-water levels provided information to better understand the effects of uncertainties in the data on simulation results.The results of this study suggest that ground-water velocity is the underlying control on the size of the zones of transport. The regional hydraulic gradient is the most significant factor controlling the shape and orientation of the zones of transport. Spatial variation in recharge, discharge, and hydraulic properties can also affect the shape of the zones of transport, however. Underestimation of porosity or overestimation of horizontal hydraulic conductivity leads to overestimation of ground-water velocity and overestimation of the size of zones of transport. Overestimation of porosity or underestimation of horizontal hydraulic conductivity leads to underestimation of ground-water velocity and underestimation of the size of zones of transport. Well discharge rate affects ground-water velocities near the well. Underestimation of discharge (and therefore velocities) will result in underestimation of the size of the zones of transport. The sensitivity of estimated zones of transport to uncertainty in parameters such as porosity and horizontal hydraulic conductivity is a function of the well discharge rate and the proximity of the well to boundaries, such as streams and rivers.

Washington↗
USGS2014

Hydroclimate of the Spring Mountains and Sheep Range, Clark County, Nevada

Precipitation, potential evapotranspiration, and actual evapotranspiration often are used to characterize the hydroclimate of a region. Quantification of these parameters in mountainous terrains is difficult because limited access often hampers the collection of representative ground data. To fulfill a need to characterize ecological zones in the Spring Mountains and Sheep Range of southern Nevada, spatially and temporally explicit estimates of these hydroclimatic parameters are determined from remote-sensing and model-based methodologies. Parameter-elevation Regressions on Independent Slopes Model (PRISM) precipitation estimates for this area ranges from about 100 millimeters (mm) in the low elevations of the study area (700 meters [m]) to more than 700 mm in the high elevations of the Spring Mountains (> 2,800 m). The PRISM model underestimates precipitation by 7–15 percent based on a comparison with four high‑elevation precipitation gages having more than 20 years of record. Precipitation at 3,000-m elevation is 50 percent greater in the Spring Mountains than in the Sheep Range. The lesser amount of precipitation in the Sheep Range is attributed to partial moisture depletion by the Spring Mountains of eastward-moving, cool-season (October–April) storms. Cool-season storms account for 66–76 percent of annual precipitation. Potential evapotranspiration estimates by the Basin Characterization Model range from about 700 mm in the high elevations of the Spring Mountains to 1,600 mm in the low elevations of the study area. The model realistically simulates lower potential evapotranspiration on northeast-to-northwest facing slopes compared to adjacent southeast-to-southwest facing slopes. Actual evapotranspiration, estimated using a Moderate Resolution Imaging Spectroradiometer based water-balance model, ranges from about 100 to 600 mm. The magnitude and spatial variation of simulated, actual evapotranspiration was validated by comparison to PRISM precipitation. Estimated groundwater recharge, computed as the residual of precipitation depleted by actual evapotranspiration, is within the range of previous estimates. A climatic water deficit dataset and aridity-index-based climate zones are derived from precipitation and evapotranspiration datasets. Climate zones range from arid in the lower elevations of the study area to humid in small pockets on north- to northeast-facing slopes in the high elevations of the Spring Mountains. Correlative analyses between hydroclimatic variables and mean ecosystem elevations indicate that the climatic water deficit is the best predictor of ecosystem distribution (R 2 = 0.92). Computed water balances indicate that substantially more recharge is generated in the Spring Mountains than in the Sheep Range. A geospatial database containing compiled and developed hydroclimatic data and other pertinent information accompanies this report.

Nevada↗
USGS1988

Hydrogeology and predevelopment flow in the Texas Gulf Coast aquifer systems

A multilayered ground-water flow system exists in the Coastal Plain sediments of Texas. The Tertiary and Quaternary clastic deposits have an area! extent of 128,000 square miles onshore and in the Gulf of Mexico. Two distinct aquifer systems are recognized for the sediments, which range in thickness from a few feet to more than 12,000 feet. The older system the Texas coastal uplands aquifer system consists of four aquifers and two confining units in the Wilcox and Claiborne Groups. It is bounded from below by the practically impermeable Midway confining unit or by the top of the geopressured zone. It is bounded from above by the poorly permeable Vicksburg-Jackson confining unit, which separates it from the younger coastal lowlands aquifer system. The coastal lowlands aquifer system consists of five permeable zones and two confining units that range in age from Oligocene to Holocene. The hydrogeologic units of both systems are exposed in bands that parallel the coastline. The units dip and thicken toward the Gulf. Quality of water in the aquifer systems varies greatly, with dissolved solids ranging from a few hundred to more than 200,000 milligrams per liter. A three-dimensional, variable-density digital model was developed to simulate predevelopment flow in the aquifer systems, for which steady-state conditions were assumed. Horizontal hydraulic conductivities of the aquifers and permeable zones in the calibrated model range from 15 feet per day for the middle Wilcox aquifer, to 170 feet per day for the Holocene-upper Pleistocene aquifer. Vertical hydraulic conductivities range from 1 x 10 -5 foot per day for the Vicksburg-Jackson confining unit, to 1 x 10 -2 foot per day for four of the aquifers and permeable zones. The simulated values of transmissivity and leakance are functions of the percent of sand that is present in each model grid block. There is a large range in precipitation across the study area, from about 21 inches per year in the west to about 56 inches per year in the east. Eastward from a line through Corpus Christi and San Antonio, average annual precipitation ranges from about 30 to about 56 inches. A few inches per year reaches the saturated zone in topographically high areas and is discharged in low areas as evapotranspiration, seepage, springflow, and stream base flow. A smaller amount of water flows through the aquifers and permeable zones downdip from the outcrop areas. This flow results in upward or downward leakage into adjacent hydrogeologic units, but is generally upward into overlying units. Westward from the line through Corpus Christi and San Antonio, average annual precipitation ranges from about 30 to about 21 inches. The general pattern of flow in the aquifers and permeable zones is similar to that in the east, but rates of flow are somewhat smaller. In contrast to the east, ground-water discharge in the west is generally not visible. Evapotranspiration is the main mechanism for ground-water discharge, with most ground water being discharged through evapotranspiration by phreatophytes. Simulated discharge and recharge rates in the combined outcrop areas of all units do not exceed 6 inches per year. The large rates occur in small, local topographically low and high areas. The average discharge rate simulated in the outcrops of the units is 0.45 inch per year. The recharge area is considerably smaller than the discharge area, and the average recharge rate over this smaller area is 0.74 inch per year. Total simulated recharge in the outcrop areas is 269 million cubic feet per day, which is offset by an equal amount of discharge in the outcrop areas. The smallest rates of leakage are across the Vicksburg-Jackson confining unit, with downward and upward rates of less than one million cubic feet per day. The greatest rate of leakage is 47 million cubic feet per day upward into the Holocene-upper Pleistocene permeable zone.

Texas↗
USGS2004

Hydrogeology and simulation of ground-water flow and land-surface subsidence in the northern part of the Gulf Coast aquifer system, Texas

As a part of the Texas Water Development Board Ground- Water Availability Modeling program, the U.S. Geological Survey developed and tested a numerical finite-difference (MODFLOW) model to simulate ground-water flow and land-surface subsidence in the northern part of the Gulf Coast aquifer system in Texas from predevelopment (before 1891) through 2000. The model is intended to be a tool that water-resource managers can use to address future ground-water-availability issues. From land surface downward, the Chicot aquifer, the Evangeline aquifer, the Burkeville confining unit, the Jasper aquifer, and the Catahoula confining unit are the hydrogeologic units of the Gulf Coast aquifer system. Withdrawals of large quantities of ground water have resulted in potentiometric surface (head) declines in the Chicot, Evangeline, and Jasper aquifers and land-surface subsidence (primarily in the Houston area) from depressurization and compaction of clay layers interbedded in the aquifer sediments. In a generalized conceptual model of the aquifer system, water enters the ground-waterflow system in topographically high outcrops of the hydrogeologic units in the northwestern part of the approximately 25,000-square-mile model area. Water that does not discharge to streams flows to intermediate and deep zones of the system southeastward of the outcrop areas where it is discharged by wells and by upward leakage in topographically low areas near the coast. The uppermost parts of the aquifer system, which include outcrop areas, are under water-table conditions. As depth increases in the aquifer system and as interbedded sand and clay accumulate, water-table conditions evolve into confined conditions. The model comprises four layers, one for each of the hydrogeologic units of the aquifer system except the Catahoula confining unit, the assumed no-flow base of the system. Each layer consists of 137 rows and 245 columns of uniformly spaced grid blocks, each block representing 1 square mile. Lateral no-flow boundaries were located on the basis of outcrop extent (northwestern), major streams (southwestern, northeastern), and downdip limit of freshwater (southeastern). The MODFLOW general-head boundary package was used to simulate recharge and discharge in the outcrops of the hydrogeologic units. Simulation of land-surface subsidence (actually, compaction of clays) and release of water from storage in the clays of the Chicot and Evangeline aquifers was accomplished using the Interbed-Storage Package designed for use with the MODFLOW model. The model was calibrated by trial-anderror adjustment of selected model input data in a series of transient simulations until the model output (potentiometric surfaces, land-surface subsidence, and selected water-budget components) reasonably reproduced field measured (or estimated) aquifer responses. Model calibration comprised four elements: The first was qualitative comparison of simulated and measured heads in the aquifers for 1977 and 2000; and quantitative comparison by computation and areal distribution of the root-mean-square error between simulated and measured heads. The second calibration element was comparison of simulated and measured hydrographs from wells in the aquifers in a number of counties throughout the modeled area. The third calibration element was comparison of simulated water-budget componentsprimarily recharge and dischargeto estimates of physically reasonable ranges of actual water-budget components. The fourth calibration element was comparison of simulated land-surface subsidence from predevelopment to 2000 to measured land surface subsidence from 1906 through 1995.

Texas↗
USGS1999

Hydrogeology and water quality of the upper Floridan aquifer, western Albany area, Georgia

Geologic, hydrologic, and water-quality data were collected to refine the hydrogeologic framework conceptual model of the Upper Floridan aquifer, and to qualitatively evaluate the potential of human activities to impact water quality in the Upper Floridan aquifer in the western Albany area, Georgia. Ground-water age dating was conducted by using chlorofluorocarbons (CFC) and tritium concentrations in water from the Upper Floridan aquifer to determine if recharge and possible contaminant migration to the aquifer is recent or occurred prior to the introduction of CFCs and tritium in the early 1950's into the global natural water system. Data were collected from core holes and wells installed during this study and previously existing wells in the Albany area. Hydrogeologic data collected during this study compare well to the regional hydrogeologic conceptual model developed during previous studies. However, the greater data density available from this study shows the dynamic and local variability in the hydrologic character of the Upper Floridan aquifer in more detail. The occurrence of sediment sizes from clay to gravel in the overburden, the absence of overburden because of erosion or sinkhole collapse, and large areas lacking surface drainage west of the Flint River provide potential areas for recharge and contaminant migration from the surface to the Upper Floridan aquifer throughout the study area. Ground-water ages generally range from 9 to 34 years, indicating that recharge consisting of 'modern' water (post early-1950's) is present in the aquifer. Ground-water ages and hydraulic heads in the Upper Floridan aquifer have an irregular distribution, indicating that localized areas of recharge to the aquifer are present in the study area. Generally, water in the Upper Floridan aquifer is calcium-bicarbonate rich, having low concentrations of magnesium, potassium, sodium, chloride, and sulfate. Water in the Upper Floridan aquifer is oxygenated, having dissolved-oxygen concentrations greater than 2 milligrams per liter. Nitrite-plus-nitrate as nitrogen, is present in the aquifer at concentrations ranging from less than 0.02 to 5.5 milligrams per liter. Areas of higher nitrate concentrations in the aquifer, coupled with widely distributed localized recharge to the aquifer indicates that suburban residential and agricultural land use in the western Albany area may affect water quality in the Upper Floridan aquifer. However, concentrations exceeding drinking water criteria were not detected in the study area. Generally, water in the Upper Floridan aquifer is calcium-bicarbonate rich, having low concentrations of magnesium, potassium, sodium, chloride, and sulfate. Water in the Upper Floridan aquifer is oxygenated, having dissolved-oxygen concentrations greater than 2 milligrams per liter. Nitrite-plus-nitrate as nitrogen, is present in the aquifer at concentrations ranging from less than 0.02 to 5.5 milligrams per liter. Areas of higher nitrate concentrations in the aquifer, coupled with widely distributed localized recharge to the aquifer indicates that suburban residential and agricultural land use in the western Albany area may affect water quality in the Upper Floridan aquifer. However, concentrations exceeding drinking water criteria were not detected in the study area.

Georgia↗
USGS1985

Ground-water flow and solute transport in the Equus beds area, south-central Kansas, 1940-79

Water levels have declined about 30 ft from 1940 to 1980 in part of the Equus beds aquifer in south-central Kansas where the city of Wichita operates a well field. A three-dimensional, finite-difference, groundwater flow model was developed to: (1) Reproduce hydrologic conditions in the flow system between the Equus beds aquifer and the underlying Wellington aquifer from 1940 to 1980, and (2) simulate the effect that future withdrawals could have on water supply in the Equus beds aquifer and on relationships between water levels in the Equus beds aquifer and the Wellington aquifer from 1980 to 2020. The model favorably reproduced both measured water levels and streamflow gains in the Equus beds aquifer for 1971 and 1980, and measured water levels in the Wellington aquifer. The flow model then was used to simulate the effects of five pumping alternatives based on rates of withdrawal by wells from 1971-79. For the first alternative, withdrawal rates were decreased by one-half. Projected saturated thickness in the aquifer and streamflow gain were the greatest among the five alternatives. For the second alternative, withdrawal rates for 1971-79 were continued. Compared to results from the first alternative, saturated thickness was projected to decline from 10 to 40 ft in some areas, and streamflow gains were maintained for the Arkansas and Little Arkansas Rivers. A two-dimensional, finite-difference, solute-transport model was developed to: (1) Reproduce the movement of chloride ion in part of the Equus beds aquifer, including the Wichita municipal well field, from 1940 to 1980, and (2) simulate the effect that future withdrawal rates could have on the concentration of chloride ion from 1980 to 2020. Sources of the chloride ion were oilfield brine disposed from 1932-43 that is moving toward the well field and water in the Arkansas River. Each simulation projected an increase in the concentration of chloride ion in the Wichita well field. The projections indicated that a continuous 1,000 mg/L source of chloride ion in streamflow losses from the Arkansas River had a greater effect on increasing chloride-ion concentrations in the Wichita well field than did the movement of residual oilfield brine. (Author 's abstract)

Water-Resources Investigations Report↗
USGS1984

Evaluation of the ground-water resources of parts of Lancaster and Berks Counties, Pennsylvania

Secondary openings in bedrock are the avenues for virtually all ground-water flow in a 626-sqare-mile area in Lancaster and Berks Counties, Pennsylvania. The number, size, and interconnection of secondary openings are functions of lithology, depth, and topography. Ground water actively circulates to depths of 150 to 300 feet below land surface. Total average annual ground-water recharge for the area is 388 million gallons per day, most of which discharges to streams from local, unconfined flow systems. A digital ground-water flow model was developed to simulate unconfined flow under several different recharge and withdrawal scenarios. On the basis of lithologic and hydrologic differences, the modeled area was sub-divided into 22 hydrogeologic units. A finite-difference grid with rectangular blocks, each 2,015 by 2,332 feet, was used. The model was calibrated under steady-state and transient conditions. The steady-state calibration was used to determine hydraulic conductivities and stream leakage coefficients and the transient calibration was used to determine specific yields. The 22 hydrogeologic units fall into four general lithologies: Carbonate rocks, metamorphic rocks, Paleozoic sedimentary rocks, and Triassic sedimentary rocks. Average hydraulic conductivity ranges from about 8.8 feet per day in carbonate units to about .5 feet per day in metamorphic units. The Stonehenge Formation (limestone) has the greatest average hydraulic conductivity--85.2 feet per day in carbonate units to about 0.11 feet per day in the greatest gaining-strem leakage coefficient--16.81 feet per day. Specific yield ranges from 0.06 to 0.09 in carbonate units, and is 0.02 to 0.015, and 0.012 in metamorphic, Paleozoic sedimentary, and Triassic sedimentary units, respectively. Transient simulations were made to determine the effects of four different combinations of natural and artificial stresses. Natural aquifer conditions (no ground-water withdrawals) and actual aquifer conditions (current ground-water withdrawals) were simulated for two years under normal seasonal and hypothetical drought (60-percent reduction in winter-spring recharge) conditions. In October, 6 months after the hypothetical drought, simulated declines in water-table altitude due to the drought occurred everywhere and ranged from a median of 3.6 feet in carbonate units to 8.7 feet in carbonate units. Simulated base flows for five major streams were reduced by 33 to 51 percent during the hypothetical drought. Also in October, maximum simulated declines in water-table altitude due to ground-water withdrawls ranged from 33 feet in carbonate units to 79 feet in Triassic sedimentary units. Simulated base flows for five major streams were reduced by the amount of ground water withdrawn. Finally, again in October, maximum simulated declines in water-table altitude due to the combination of hypothetical drought and ground-water withdrawls ranged from 38 feet in carbonate units to 109 feet in Triassic sedimentary units. Due to aquifer dewatering, simulated declines were as much as 24 feet greater than the sum of the separate simulated declines that were caused by hypothetical drought and ground-water withdrawals. Some of the greatest simulated declines were in well fields, operated by three municipalities that experienced water-supply problems during the 1980-81 drought.

Water-Resources Investigations Report↗
USGS2008

Simulations of ground-water flow and particle pathline analysis in the zone of contribution of a public-supply well in Modesto, eastern San Joaquin Valley, California

Shallow ground water in the eastern San Joaquin Valley is affected by high nitrate and uranium concentrations and frequent detections of pesticides and volatile organic compounds (VOC), as a result of ground-water development and intensive agricultural and urban land use. A single public-supply well was selected for intensive study to evaluate the dominant processes affecting the vulnerability of public-supply wells in the Modesto area. A network of 23 monitoring wells was installed, and water and sediment samples were collected within the approximate zone of contribution of the public-supply well, to support a detailed analysis of physical and chemical conditions and processes affecting the water chemistry in the well. A three-dimensional, steady-state local ground-water-flow and transport model was developed to evaluate the age of ground water reaching the well and to evaluate the vulnerability of the well to nonpoint source input of nitrate and uranium. Particle tracking was used to compute pathlines and advective travel times in the ground-water flow model. The simulated ages of particles reaching the public-supply well ranged from 9 to 30,000 years, with a median of 54 years. The age of the ground water contributed to the public-supply well increased with depth below the water table. Measured nitrate concentrations, derived primarily from agricultural fertilizer, were highest (17 milligrams per liter) in shallow ground water and decreased with depth to background concentrations of less than 2 milligrams per liter in the deepest wells. Because the movement of water is predominantly downward as a result of ground-water development, and because geochemical conditions are generally oxic, high nitrate concentrations in shallow ground water are expected to continue moving downward without significant attenuation. Simulated long-term nitrate concentrations indicate that concentrations have peaked and will decrease in the public-supply well during the next 100 years because of the low nitrate concentrations in recharge beneath the urban area and the increasing proportion of urban-derived ground water reaching the well. The apparent lag time between peak input concentrations and peak concentrations in the well is about 20 to 30 years. Measured uranium concentrations were also highest (45 micrograms per liter) in shallow ground water, and decreased with depth to background concentrations of about 0.5 microgram per liter. Naturally-occurring uranium adsorbed to aquifer sediments is mobilized by oxygen-rich, high-alkalinity water. Alkalinity increased in shallow ground water in response to agricultural development. As ground-water pumping increased in the 1940s and 1950s, this alkaline water moved downward through the ground-water flow system, mobilizing the uranium adsorbed to aquifer sediments. Ground water with high alkalinity and high uranium concentrations is expected to continue to move deeper in the system, resulting in increased uranium concentrations with depth in ground water. Because alkalinity (and correspondingly uranium) concentrations were high in shallow ground water beneath both the urban and the agricultural land, long-term uranium concentrations in the public-supply well are expected to increase as the proportion of uranium-affected water contributed to the well increases. Assuming that the alkalinity near the water table remains the same, the simulation of long-term alkalinity in the public-supply well indicates that uranium concentrations in the public-supply well will likely approach the maximum contaminant level; however, the time to reach this level is more than 100 years because of the significant proportion of old, unaffected water at depth that is contributed to the public-supply well.

California↗