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,009 records · Page 56Linked to original sources
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↗
USGS2008

Water quality, hydrology, and simulated response to changes in phosphorus loading of Butternut Lake, Price and Ashland Counties, Wisconsin, with special emphasis on the effects of internal phosphorus loading in a polymictic lake

Butternut Lake is a 393-hectare, eutrophic to hypereutrophic lake in northcentral Wisconsin. After only minor improvements in water quality were observed following several actions taken to reduce the nutrient inputs to the lake, a detailed study was conducted from 2002 to 2007 by the U.S. Geological Survey to better understand how the lake functions. The goals of this study were to describe the water quality and hydrology of the lake, quantify external and internal sources of phosphorus, and determine the effects of past and future changes in phosphorus inputs on the water quality of the lake. Since the early 1970s, the water quality of Butternut Lake has changed little in response to nutrient reductions from the watershed. The largest changes were in near-surface total phosphorus concentrations: August concentrations decreased from about 0.09 milligrams per liter (mg/L) to about 0.05 mg/L, but average summer concentrations decreased only from about 0.055—0.060 mg/L to about 0.045 mg/L. Since the early 1970s, only small changes were observed in chlorophyll a concentrations and water clarity (Secchi depths). All major water and phosphorus sources, including the internal release of phosphorus from the sediments (internal loading), were measured directly, and minor sources were estimated to construct detailed water and phosphorus budgets for the lake during monitoring years (MY) 2003 and 2004. During these years, Butternut Creek, Spiller Creek, direct precipitation, small tributaries and near-lake drainage area, and ground water contributed about 62, 20, 8, 7, and 3 percent of the inflow, respectively. The average annual load of phosphorus to the lake was 2,540 kilograms (kg), of which 1,590 kg came from external sources (63 percent) and 945 kg came from the sediments in the lake (37 percent). Of the total external sources, Butternut Creek, Spiller Creek, small tributaries and near-lake drainage area, septic systems, precipitation, and ground water contributed about 63, 23, 9, 3, 1, and 1 percent, respectively. Because of the high internal phosphorus loading, the eutrophication models used in this study were unable to simulate the observed water-quality characteristics in the lake without incorporating this source of phosphorus. However, when internal loading of phosphorus was added to the BATHTUB model, it accurately simulated the average water-quality characteristics measured in MY 2003 and 2004. Model simulations demonstrated a relatively linear response between in-lake total phosphorus concentrations and external phosphorus loading; however, the changes in concentrations were smaller than the changes in external phosphorus loadings (about 25—40 percent of the change in phosphorus loading). Changes in chlorophyll a concentrations, the percentage of days with algal blooms, and Secchi depths were nonlinear and had a greater response to reductions in phosphorus loading than to increases in phosphorus loading. A 50-percent reduction in external phosphorus loading caused an 18-percent decrease in chlorophyll a concentrations, a 41-percent decrease in the percentage of days with algal blooms, and a 12-percent increase in Secchi depth. When the additional internal phosphorus loading was removed from model simulations, all of these constituents showed a much greater response to changes in external phosphorus loading. Because of Butternut Lake's morphometry, it is polymictic, which means it mixes frequently and does not develop stable thermal stratification throughout the summer. This characteristic makes it more vulnerable than dimictic lakes, which mix in spring and fall and develop stable thermal stratification during summer, to the high internal phosphorus loading that has resulted from historically high, nonnatural, external phosphorus loading. In polymictic lakes, the phosphorus released from the sediments is mixed into the upper part of the lake throughout summer. Once Butternut Lake became hypereutrophic (very productive), it became very difficult to alter its trophic state through reductions in external phosphorus loading because the high internal loading does not respond quickly to reductions in external nutrient loading. For Butternut Lake to become significantly less productive (change to a borderline mesotrophic/eutrophic state) a combined approach to reduce or eliminate internal phosphorus loading and reduce the external phosphorus loading by about 50 percent is needed.

Wisconsin↗
USGS1997

Geochemical analyses of ground-water ages, recharge rates, and hydraulic conductivity of the N aquifer, Black Mesa area, Arizona

The Navajo Nation and Hopi Tribe of the Black Mesa area, Arizona, depend on ground water from the N aquifer to meet most tribal and industrial needs. Increasing use of this aquifer is creating concerns about possible adverse effects of increased ground-water withdrawals on the water resources of the region. A thorough understanding of the N aquifer is necessary to assess the aquifer's response to ground-water withdrawals. This study used geochemical techniques as an independent means of improving the conceptual model of ground-water flow in the N aquifer and to estimate recharge rates and hydraulic conductivity. Ground water flows in a south-southeastward direction from the recharge area around Shonto into the confined part of the N aquifer underneath Black Mesa. Ground-water flow paths diverge in the confined part of the aquifer to the northeast and south. The N aquifer thins to extinction south of Black Mesa. This discontinuity could force ground water to diverge along paths of least resistance. Ground water discharges from the confined part of the aquifer into Laguna Creek and Moenkopi Wash and from springs southwest of Kykotsmovi and southeast of Rough Rock after a residence time of about 35,000 years or more. Recent recharge along the periphery of Black Mesa mixes with older ground water that discharges from the confined part of the aquifer and flows away from Black Mesa. Dissolved-ion concentrations, ratios of dissolved ions, dissolved-gas concentrations, tritium, carbon-13, and chlorine-36 data indicate that water in the overlying D aquifer could be leaking into the confined part of the N aquifer in the southeastern part of Black Mesa. The boundary between the leaky and nonleaky zones is defined roughly by a line from Rough Rock to Second Mesa and separates ground waters that have significantly different chemistries. The Dakota Sandstone and Entrada Formation of the D aquifer could be the sources of leakage. Adjusted radiocarbon ground-water ages and data on isotopes of oxygen and hydrogen indicate that more than 90 percent of the water in the confined part of the N aquifer is older than 10,000 years and was recharged during glacial periods. Estimates of recharge rates made on the basis of ground-water ages, aquifer thicknesses, and assumed porosities indicate that the annual average recharge rate in the northwestern part of the study area during the glacial periods was about four times the average annual rate of the past 10,000 years, and that recharge rates for the past 10,000 years are less than modern recharge rates assumed in a previous study. Estimates of horizontal hydraulic conductivity were 0.95 and 1.16 feet per day for the northeast and southwest flow paths, respectively. These values are within the range of hydraulic conductivities calculated from aquifer tests, which ranged from 0.05 to 2.1 feet per day and averaged 0.65 foot per day.

Arizona↗
USGS2005

Questa baseline and pre-mining ground-water quality investigation. 12. Geochemical and reactive-transport modeling based on tracer injection-synoptic sampling studies for the Red River, New Mexico, 2001-2002

Reactive-transport processes in the Red River, downstream from the town of Red River in north-central New Mexico, were simulated using the OTEQ reactive-transport model. The simulations were calibrated using physical and chemical data from synoptic studies conducted during low-flow conditions in August 2001 and during March/April 2002. Discharge over the 20-km reach from the town of Red River to the USGS streamflow-gaging station near the town of Questa ranged from 395 to 1,180 L/s during the 2001 tracer and from 234 to 421 L/s during the 2002 tracer. The pH of the Red River ranged from 7.4 to 8.5 during the 2001 tracer and from 7.1 to 8.7 during the 2002 tracer, and seep and tributary samples had pH values of 2.8 to 9.0 during the 2001 tracer and 3.8 to 7.2 during the 2002 tracer. Mass-loading calculations allowed identification of several specific locations where elevated concentrations of potential contaminants entered the Red River . These locations, characterized by features on the north side of the Red River that are known to be sources of low-pH water containing elevated metal and sulfate concentrations, are: the initial 2.4 km of the study reach, including Bitter Creek, the stream section from 6.2 to 7.8 km, encompassing La Bobita well and the Hansen debris fan, Sulphur Gulch, at about 10.5 km, the area near Portal Springs, from 12.2 to 12.6 km, and the largest contributors of mass loading, the 13.7 to 13.9 km stream section near Cabin Springs and the 14.7 to 17.5 km stream section from Shaft Spring to Thunder Bridge, Goathill Gulch, and Capulin Canyon. Speciation and saturation index calculations indicated that although solubility limits the concentration of aluminum above pH 5.0, at pH values above 7 and aluminum concentrations below 0.3 mg/L inorganic speciation and mineral solubility controls no longer dominate and aluminum-organic complexing may occur. The August 2001 reactive-transport simulations included dissolved iron(II) oxidation, constrained using measured concentrations of dissolved iron(II) and dissolved iron(total). Both simulations included precipitation of amorphous Al(OH) 3 and hydrous ferric oxide as Fe(OH) 3 , and sorption of copper and zinc to the precipitated hydrous ferric oxide. Simulations revealed that hydrogen, iron, aluminum, copper, and zinc were non-conservative and that mineral precipitation can account for iron and aluminum concentrations. Copper and zinc concentrations can be accounted for by simulating their sorption to hydrous ferric oxide forming in the water column of the Red River , although hydrous manganese oxides also may be important sorption substrates.

New Mexico↗
USGS1985

Ground-water flow in the Prairie du Chien-Jordan aquifer related to contamination by coal-tar derivatives, St. Louis Park, Minnesota

A three-dimensional, ground-water-flow model of the Prairie du Chien-Jordan aquifer and associated hydrogeologic units was developed to evaluate the movement of coal-tar derivatives from a coal-tar distillation and wood-preserving plant in St. Louis Park, Minnesota. A finite-difference grid was superimposed on the modeled area, which includes most of eastern Hennepin County. The individual cells are 400-foot squares in the center of the grid (St. Louis Park area); the cells increase in dimension toward the outside limits of the grid. Five geologic units are represented by four layers in the model. These units include the Jordan Sandstone, the Prairie du Chien Group (dolomite and sandy dolomite), the basal confining unit of the St. Peter Sandstone (silty and sandy shale), the St. Peter Sandstone, and glacial deposits in bedrock valleys. The model was calibrated for steady-state conditions for a period before significant ground-water development (1885-1930) and for a period of significant pumping stress (winter conditions, 1970 fs). A transient calibration was accomplished by simulation of a period during which seasonal changes in potentiometric head in the Prairie du Chien-Jordan aquifer were significant (1977-80). Sensitivity testing indicated that leakage to the upper model layer and the vertical hydraulic conductivity of the basal confining unit of the St. Peter Sandstone were the model hydrologic properties which, when changed, resulted in the greatest changes in model-calculated water levels. The calibrated model generally calculates water levels that are within 10 feet of measured values. Model simulations indicate that the potentiometric surface of the Prairie du Chien-Jordan aquifer would be raised by as much as 3 feet in the area of the plant site by water introduced into the aquifer through wells open to more than one aquifer system. The cones of impression created at these wells could have a significant impact on the transport of contaminants in the Prairie du Chien-Jordan. The presence of coal-tar derivatives in the aquifer has been difficult to explain in wells located upgradient from the plant site to the north, west, and southwest. Simulations suggest that, during periods of heavy withdrawal from certain of these wells (SLP10, SLP15, and SLP5), local hydraulic gradients may have been altered, resulting in the potential for the movement of contaminants from the area of the plant site to the wells. Cones of impression at multiaquifer wells near the plant site contributed to the alteration of local gradients. Simulation of a proposed gradient-control plan, in which lateral homogeneity and isotropy of individual hydro geologic units was assumed, indicates that the actions would be effective in limiting expansion of the contaminated volume in the Prairie du Chien-Jordan aquifer. The plan includes the control of withdrawal from five wells. The simulations also show, however, that modelcalculated potentiometric surfaces are sensitive to changes in withdrawal rates at wells not intended to be under the control of the plan. Management of discharge from these wells also will be important to overall effectiveness of the remedial-action plan.

Minnesota↗
USGS2004

Water flow and nutrient flux from five estuarine rivers along the southwest coast of the Everglades National Park, Florida, 1997-2001

Discharge and nutrient fluxes for five tidally affected streams were monitored and evaluated as a part of the U.S. Geological Survey Place-Based Studies Initiative and the U.S. Department of the Interior Critical Ecosystem Studies Initiative. Locations on Lostmans Creek, and Broad, Harney, Shark, and North Rivers were selected using the criterion that a large amount of the water that flows through Shark River Slough must pass these sites. Discharge and nutrient-concentration data collection started at the Broad, Harney, and Shark River stations in January 1997 and ended in early 2001. Discharge and nutrient-concentration data collection started at the Lostmans Creek and North River stations in April 1999 and ended in early 2001. Each station was equipped with a vertically oriented acoustic-velocity sensor, water-level pressure transducer, bottom water-temperature thermistor, and specific conductance four-electrode sensor. Data collected using a vessel-mounted acoustic discharge measurement system were used to calibrate regression models of the mean river velocities and the in-situ index velocities. Information from these stations, in conjunction with data from other ongoing studies, will help to determine environmental effects on the southwest coast estuaries as changes in water management of the Everglades National Park continue. Discharges from the Lostmans Creek, and Broad, Harney, Shark, and North River stations are influenced by semidiurnal tides, meteorological events, and surface- and ground-water inflow. Each of the five rivers is usually well mixed, having no greater than 500 microSiemens per centimeter at 25? Celsius difference in specific conductance from top to bottom during flood and ebb tides. Instantaneous flood discharges (water moving upstream) are typically of greater magnitude and shorter duration than instantaneous ebb discharges (water moving downstream). Instantaneous discharge data were filtered using a low-pass filter to remove predominant tidal frequencies, and the filtered data were used to compute daily mean and monthly mean residual discharges. Lostmans Creek, and Broad, Harney and Shark Rivers each contributed from 20 to 27 percent of the total measured discharge to the Gulf of Mexico, whereas North River contributed approximately 4 percent. The main discharge region of the Shark River Slough extends from as far north as Lostmans Creek to as far south as North River. North River discharge has similar response characteristics to the other four rivers measured, but with a lesser magnitude of discharge. Comparisons of monthly mean discharges from the Tamiami Canal flow control structures S-12-A, B, C, and D located on U.S. Highway 41 (Tamiami Trail) to the five station total monthly mean discharges indicate that the discharges from the five rivers are approximately 2 to 3 times the S-12-A, B, C, D discharges, and that the measured southwest coast discharge peaks lead the S-12-A, B, C, D discharge peaks by approximately 1 month. Residual total nitrogen and total phosphorus fluxes were estimated using linear regression models of discharge and flux. Monthly mean total nitrogen residual fluxes for the five southwest coast rivers ranged from approximately 0 to 390 short tons, whereas monthly mean total phosphorus residual fluxes ranged from approximately 0 to 6 short tons. Total nitrogen and total phosphorus residual fluxes at Lostmans Creek, and Broad, Harney, and Shark Rivers were similar in magnitude, each accounting for between 20 to 29 percent of the total measured residual flux. North River contributed between 3 to 4 percent of the total nitrogen and total phosphorus residual flux from the five rivers.

Florida↗
USGS2006

Pesticide toxicity index for freshwater aquatic organisms, 2nd edition

The U.S. Geological Survey's National Water-Quality Assessment (NAWQA) Program is designed to assess current water-quality conditions, changes in water quality over time, and the effects of natural and human factors on water quality for the Nation's streams and ground-water resources. For streams, one of the most difficult parts of the assessment is to link chemical conditions to effects on aquatic biota, particularly for pesticides, which tend to occur in streams as complex mixtures with strong seasonal patterns. A Pesticide Toxicity Index (PTI) was developed that combines pesticide exposure of aquatic biota (measured concentrations of pesticides in stream water) with acute toxicity estimates (standard endpoints from laboratory bioassays) to produce a single index value for a sample or site. The development of the PTI was limited to pesticide compounds routinely measured in NAWQA studies and to toxicity data readily available from existing databases. Qualifying toxicity data were found for one or more types of test organisms for 124 of the 185 pesticide compounds measured in NAWQA samples, but with a wide range of available bioassays per compound (1 to 232). In the databases examined, there were a total of 3,669 bioassays for the 124 compounds, including 398 48-hour EC50 values (concentration at which 50 percent of test organisms exhibit a sublethal response) for freshwater cladocerans, 699 96-hour LC50 values (concentration lethal to 50 percent of test organisms) for freshwater benthic invertebrates, and 2,572 96-hour LC50 values for freshwater fish. The PTI for a particular sample is the sum of toxicity quotients (measured concentration divided by the median toxicity concentration from bioassays) for each detected pesticide, and thus, is based on the concentration addition model of pesticide toxicity. The PTI can be calculated for specific groups of pesticides and for specific taxonomic groups. Although the PTI does not determine whether water in a sample is toxic to aquatic organisms, its values can be used to rank or compare the toxicity of samples or sites on a relative basis for use in further analysis or additional assessments. The PTI approach may be useful as a basis for comparing the potential significance of pesticides in different streams on a common basis, for evaluating relations between pesticide exposure and observed biological conditions, and for prioritizing where further studies are most needed.

Scientific Investigations Report↗
USGS1994

Ground-water quality and geochemistry in Dayton, Stagecoach, and Churchill Valleys, western Nevada

The U.S. Geological Survey investigated the quality of ground water in the Dayton, Stagecoach, and Churchill Valleys as part of the Carson River Basin National Water-Quality Assessment (NAWQA) pilot study. Four aquifer systems have been de- lineated in the study area. Principal aquifers are unconsolidated deposits at altitudes of less than 4,900 feet above sea level and more than 50 feet below land surface. Shallow aquifers are at altitudes of less than 4,900 feet and less than 50 feet below land surface. Upland aquifers are above 4,900 feet and provide recharge to the principal aquifers. Thermal aquifers, defined as those having a water temperature greater than 30 degrees Celsius, are also present. Ground water used in Dayton, Stagecoach, and Churchill Valleys is pumped from principal aquifers in unconsolidated basin-fill deposits. Ground water in these aquifers originates as precipitation in the adjacent mountains and is recharged by the Carson River and by underflow from adjacent upstream valleys. Ground-water flow is generally parallel to the direction of surface-water flow in the Carson River. Ground water is discharged by pumping, evapo- transpiration, and underflow into the Carson River. The results of geochemical modeling indicate that as ground water moves from upland aquifers in mountainous recharge areas to principal aquifers in basin-fill deposits, the following processes probably occur: (1) plagioclase feldspar, sodium chloride, gypsum (or pyrite), potassium feldspar, and biotite dissolve; (2) calcite precipitates; (3) kaolinite forms; (4) small amounts of calcium and magnesium in the water exchange for potassium on aquifer minerals; and (5) carbon dioxide is gained or lost. The geochemical models are consistent with (1) phases identified in basin- fill sediments; (2) chemical activity of major cations and silica; (3) saturation indices of calcite and amorphous silica; (4) phase relations for aluminosilicate minerals indicated by activity diagrams; and (5) results of optical, X-ray diffraction, and scanning-electron microscopy examination of mineral grains in the aquifer sediments. Sulfur-isotopic composition of ground- water samples also supports the models. In general, the quality of ground water in the study area meets Nevada State drinking-water standards and is acceptable for most uses. In addition to analysis for major ions, samples were analyzed for 22 inorganic trace elements, 3 nutrients, and 4 radionuclides. Selenium in 1 sample is the only constituent that exceeded Nevada State primary drinking-water standards. Nevada State secondary- drinking water standards were exceeded for fluoride in 1 sample, for iron in 7 samples, and for manganese in 19 samples. Minor constituent con- centrations are generally the result of local redox conditions, and are primarily from minerals in volcanic and marine metasedimentary rocks, metal- oxide coatings on mineral grains, and organic matter.

Open-File Report↗