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,063 records · Page 59Linked to original sources

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Simulation of Multiscale Ground-Water Flow in Part of the Northeastern San Joaquin Valley, California

The transport and fate of agricultural chemicals in a variety of environmental settings is being evaluated as part of the U.S. Geological Survey (USGS) National Water-Quality Assessment Program. One of the locations being evaluated is a 2,700-km2 (square kilometer) regional study area in the northeastern San Joaquin Valley surrounding the city of Modesto, an area dominated by irrigated agriculture in a semi-arid climate. Ground water is a key source of water for irrigation and public supply, and exploitation of this resource has altered the natural flow system. The aquifer system is predominantly alluvial, and an unconfined to semiconfined aquifer overlies a confined aquifer in the southwestern part of the study area; these aquifers are separated by the lacustrine Corcoran Clay. A regional-scale 16-layer steady-state model of ground-water flow in the aquifer system in the regional study area was developed to provide boundary conditions for an embedded 110-layer steady-state local-scale model of part of the aquifer system overlying the Corcoran Clay along the Merced River. The purpose of the local-scale model was to develop a better understanding of the aquifer system and to provide a basis for simulation of reactive transport of agricultural chemicals. The heterogeneity of aquifer materials was explicitly incorporated into the regional and local models using information from geologic and drillers? logs of boreholes. Aquifer materials were differentiated in the regional model by the percentage of coarse-grained sediments in a cell, and in the local model by four hydrofacies (sand, silty sand, silt, and clay). The calibrated horizontal hydraulic conductivity values of the coarse-grained materials in the zone above the Corcoran Clay in the regional model and of the sand hydrofacies used in the local model were about equal (30?80 m/d [meter per day]), and the vertical hydraulic conductivity values in the same zone of the regional model (median of 0.012 m/d), which is dominated by the finer-grained materials, were about an order of magnitude less than that for the clay hydrofacies in the local model. Data used for calibrating both models included long-term hourly water-level measurements in 20 short-screened wells installed by the USGS in the Modesto and Merced River areas. Additional calibration data for the regional model included water-level measurements in 11 wells upslope and 17 wells downslope from these areas. The root mean square error was 2.3 m (meter) for all wells in the regional model and 0.8 m for only the USGS wells; the associated average errors were 0.9 m and 0.3 m, respectively. The root mean square error for the 12 USGS wells along a transect in the local model area was 0.08 m; the average error was 0.0 m. Particle tracking was used with the local model to estimate the concentration of an environmental tracer, sulfur hexafluoride, in 10 USGS transect wells near the Merced River that were sampled for this constituent. Measured and estimated concentrations in the mid-depth and deepest wells, which would be most sensitive to errors in hydraulic conductivity estimates, were consistent. The combined results of particle tracking and sulfur hexafluoride analysis suggest that most water sampled from the transect wells was recharged less that 25 years ago.

Scientific Investigations Report↗

Determination of hydrologic properties needed to calculate average linear velocity and travel time of ground water in the principal aquifer underlying the southeastern part of Salt Lake Valley, Utah

A 48-square-mile area in the southeastern part of the Salt Lake Valley, Utah, was studied to determine if generalized information obtained from geologic maps, water-level maps, and drillers' logs could be used to estimate hydraulic conduc- tivity, porosity, and slope of the potentiometric surface: the three properties needed to calculate average linear velocity of ground water. Estimated values of these properties could be used by water- management and regulatory agencies to compute values of average linear velocity, which could be further used to estimate travel time of ground water along selected flow lines, and thus to determine wellhead protection areas around public- supply wells. The methods used to estimate the three properties are based on assumptions about the drillers' descriptions, the depositional history of the sediments, and the boundary con- ditions of the hydrologic system. These assump- tions were based on geologic and hydrologic infor- mation determined from previous investigations. The reliability of the estimated values for hydro- logic properties and average linear velocity depends on the accuracy of these assumptions. Hydraulic conductivity of the principal aquifer was estimated by calculating the thickness- weighted average of values assigned to different drillers' descriptions of material penetrated during the construction of 98 wells. Using these 98 control points, the study area was divided into zones representing approximate hydraulic- conductivity values of 20, 60, 100, 140, 180, 220, and 250 feet per day. This range of values is about the same range of values used in developing a ground-water flow model of the principal aquifer in the early 1980s. Porosity of the principal aquifer was estimated by compiling the range of porosity values determined or estimated during previous investigations of basin-fill sediments, and then using five different values ranging from 15 to 35 percent to delineate zones in the study area that were assumed to be underlain by similar deposits. Delineation of the zones was based on depositional history of the area and the distri- bution of sediments shown on a surficial geologic map. Water levels in wells were measured twice in 1990: during late winter when ground-water with- drawals were the least and water levels the highest, and again in late summer, when ground- water withdrawals were the greatest and water levels the lowest. These water levels were used to construct potentiometric-contour maps and subsequently to determine the variability of the slope in the potentiometric surface in the area. Values for the three properties, derived from the described sources of information, were used to produce a map showing the general distribution of average linear velocity of ground water moving through the principal aquifer of the study area. Velocity derived ranged from 0.06 to 144 feet per day with a median of about 3 feet per day. Values were slightly faster for late summer 1990 than for late winter 1990, mainly because increased with- drawal of water during the summer created slightly steeper hydraulic-head gradients between the recharge area near the mountain front and the well fields farther to the west. The fastest average linear-velocity values were located at the mouth of Little Cottonwood Canyon and south of Dry Creek near the mountain front, where the hydraulic con- ductivity was estimated to be the largest because the drillers described the sediments to be pre- dominantly clean and coarse grained. Both of these areas also had steep slopes in the potentiometric surface. Other areas where average linear velocity was fast included small areas near pumping wells where the slope in the potentiometric surface was locally steepened. No apparent relation between average linear velocity and porosity could be seen in the mapped distributions of these two properties. Calculation of travel time along a flow line to a well in the southwestern part of the study area during the sum

Water-Resources Investigations Report↗

Ground-water flow near two radioactive-waste-disposal areas at the Western New York Nuclear Service Center, Cattaraugus County, New York: Results of flow simulation

Two adjacent burial areas were excavated in a clay-rich till at a radioactive waste disposal site near West Valley in Cattaraugus County, N.Y.: (1) which contains mainly low-level radioactive wastes generated onsite by a nuclear fuel reprocessing plant, has been in operation since 1966; and (2) which contains commercial low-level radioactive wastes, was operated during 1963-75. Groundwater below the upper 3 meters of till generally moves downward through a 20- to 30-meter thick sequence of tills underlain by lacustrine and kame-delta deposits of fine sand and silt. Groundwater in the weathered, upper 3 meters of till can move laterally for several meters before either moving downward into the kame-delta deposits or discharging to the land surface. A two-dimensional finite-element model that simulates two vertical sections was used to evaluate hydrologic factors that control groundwater flow in the till. Conditions observed during March 1983 were reproduced accurately in steady-state simulations that used four isotropic units of differing hydraulic conductivity to represent two fractured and weathered till units near land surfaces, an intermediate group of isolated till zones that contain significant amounts of fine sand and silt, and a sequence of till units at depths that have been consolidated by overburden pressure. Recharge rates used in the best-fit simulation ranged from 1.4 cm/yr along smooth, sloping or compacted surfaces to 3.8 cm/yr near swampy areas. Values of hydraulic conductivity and infiltration used in the calibrated best-fit model were nearly identical to values used in a previous model analysis of the nearby commercial-waste burial area. Results of the model simulations of a burial pit assumed to be filled with water indicate that water near the bottom of the burial pit would migrate laterally in the shallow, weathered till for 5 to 6 meters before moving downward into the unweathered till, and water near the top of the pit would move laterally less than 20 meters before moving downward into the unweathered till. These results indicate that subsurface migration of radionuclides in groundwater to points of discharge to land surface is unlikely as long as the water level does not rise into the reworked cover material.

New York↗

Arsenic in ground water in selected parts of southwestern Ohio, 2002-03

Arsenic concentrations were measured in 57 domestic wells in Preble, Miami, and Shelby Counties, in southwestern Ohio. The median arsenic concentration was 7.1 ?g/L (micrograms per liter), and the maximum was 67.6 ?g/L. Thirty-seven percent of samples had arsenic concentrations greater than the U.S. Environmental Protection Agency drinking-water standard of 10 ?g/L. Elevated arsenic concentrations (>10 ?g/L) were detected over the entire range of depths sampled (42 to 221 feet) and in each of three aquifer types, Silurian carbonate bedrock, glacial buried-valley deposits, and glacial till with interbedded sand and gravel. One factor common in all samples with elevated arsenic concentrations was that iron concentrations were greater than 1,000 ?g/L. The observed correlations of arsenic with iron and alkalinity are consistent with the hypothesis that arsenic was released from iron oxides under reducing conditions (by reductive dissolution or reductive desorption). Comparisons among the three aquifer types revealed some differences in arsenic occurrence. For buried-valley deposits, the median arsenic concentration was 4.6 ?g/L, and the maximum was 67.6 ?g/L. There was no correlation between arsenic concentrations and depth; the highest concentrations were at intermediate depths (about 100 feet). Half of the buried-valley samples were estimated to be methanic. Most of the samples with elevated arsenic concentrations also had elevated concentrations of dissolved organic carbon and ammonia. For carbonate bedrock, the median arsenic concentration was 8.0 ?g/L, and the maximum was 30.7 ?g/L. Arsenic concentrations increased with depth. Elevated arsenic concentrations were detected in iron- or sulfate-reducing samples. Arsenic was significantly correled with molybdenum, strontium, fluoride, and silica, which are components of naturally ocurring minerals. For glacial till with interbedded sand and gravel, half of the samples had elevated arsenic concentrations. The median was 11.4 ?g/L, and the maximum was 27.6 ?g/L. At shallow depths (<100 feet), this aquifer type had higher arsenic and iron concentrations than carbonate bedrock. It is not known whether these observed differences among aquifer types are related to variations in (1) arsenic content of the aquifer material, (2) organic carbon content of the aquifer material, (3) mechanisms of arsenic mobilization (or uptake), or (4) rates of arsenic mobilization (or uptake). A followup study that includes solid-phase analyses and geochemical modeling was begun in 2004 in northwestern Preble County.

Scientific Investigations Report↗

Vulnerability of recently recharged ground water in the High Plains aquifer to nitrate contamination

Nitrate concentrations greater than background levels have been detected in ground water of the High Plains aquifer. Empirically based models and corresponding maps were developed that predict the vulnerability of the aquifer to nonpoint-source nitrate contamination. The models predict the probability of detecting nitrate concentrations larger than 4 milligrams per liter in ground water of the High Plains aquifer that was recharged during the last 50 years. The models were calibrated by correlating concentrations of nitrate in ground water from 336 wells that intercept recently recharged (less than 50 years) water with anthropogenic and hydrogeologic explanatory variables. Particle-tracking simulations delineated well-contributing areas and determined well-screen depths that intercept recently recharged ground water. The models were developed using multivariate logistic regression, and a map was generated from these models using a geographic information system. Two multivariate logistic regression models of vulnerability were found to have the most statistical significance and the best model fit and predictive ability. The two models represent the Northern High Plains and the combined Central and Southern High Plains, and they indicate that ground-water vulnerability of the entire High Plains aquifer is best explained by the spatial distribution of nonirrigated and irrigated agricultural lands, organic matter of the soil, depth to the regional water table, and clay content of the unsaturated zone. Vulnerability of the Northern High Plains is greater in areas that have more nonirrigated and irrigated agricultural lands and less organic matter in the soil. The vulnerability of the Central and Southern High Plains also is greater in areas that have more nonirrigated and irrigated agricultural lands and also in areas with shallow depths to water table and less clay in the unsaturated zone. The majority (53.3 percent) of the High Plains aquifer has less than a 40-percent predicted probability of nitrate concentrations larger than 4 milligrams per liter. Approximately 21.1 percent of the High Plains aquifer has a relatively high (greater than 60 percent) predicted probability of nitrate concentrations greater than or equal to 4 milligrams per liter. Areas with relatively high predicted probability are located in the southwestern, southern, and eastern areas of the Northern High Plains, in the eastern arm of the Central High Plains, and in southern areas of the Southern High Plains. Areas of the aquifer with relatively low (less than 40 percent) predicted vulnerability to nitrate concentrations greater than or equal to 4 milligrams per liter are located in the northwestern and north-central areas of the Northern High Plains, the central and southern areas of the Central High Plains, and a band across the north-central part of the Southern High Plains. Uncertainty of these vulnerability predictions was estimated by Latin hypercube sampling to address propagation of model and data errors inherently associated with estimates of model coefficients and explanatory variables. Results of the Latin hypercube sampling simulations are presented as uncertainty maps of the lower 5th and upper 95th percentile of the output probability distribution, which represents the 90-percent prediction interval that contains the true probability of detecting nitrate greater than or equal to 4 milligrams per liter. Generally, these uncertainty maps show greater prediction uncertainty in areas with relatively higher predicted vulnerability and lower uncertainty in areas with relatively lower predicted vulnerability.

Scientific Investigations Report↗

Evaluation of geohydrologic framework, recharge estimates and ground-water flow of the Joshua Tree area, San Bernardino County, California

Ground water historically has been the sole source of water supply for the community of Joshua Tree in the Joshua Tree ground-water subbasin of the Morongo ground-water basin in the southern Mojave Desert. The Joshua Basin Water District (JBWD) supplies water to the community from the underlying Joshua Tree ground-water subbasin. The JBWD is concerned with the long-term sustainability of the underlying aquifer. To help meet future demands, the JBWD plans to construct production wells in the adjacent Copper Mountain ground-water subbasin. As growth continues in the desert, there may be a need to import water to supplement the available ground-water resources. In order to manage the ground-water resources and to identify future mitigating measures, a thorough understanding of the ground-water system is needed. The purpose of this study was threefold: (1) improve the understanding of the geohydrologic framework of the Joshua Tree and Copper Mountain ground-water subbasins, (2) determine the distribution and quantity of recharge using field and numerical techniques, and (3) develop a ground-water flow model that can be used to help manage the water resources of the region. The geohydrologic framework was refined by collecting and interpreting water-level and water-quality data, geologic and electric logs, and gravity data. The water-bearing deposits in the Joshua Tree and Copper Mountain ground-water subbasins are Quarternary alluvial deposits and Tertiary sedimentary and volcanic deposits. The Quarternary alluvial deposits were divided into two aquifers (referred to as the 'upper' and the 'middle' alluvial aquifers), which are about 600 feet (ft) thick, and the Tertiary sedimentary and volcanic deposits were assigned to a single aquifer (referred to as the 'lower' aquifer), which is as thick as 1,500 ft. The ground-water quality of the Joshua Tree and Copper Mountain ground-water subbasins was defined by collecting 53 ground-water samples from 15 wells (10 in the Joshua Tree ground-water subbasin and 5 in the Copper Mountain ground-water subbasin) between 1980 and 2002 and analyzing the samples for major ions, nutrients, and selected trace elements. Selected samples also were analyzed for oxygen-18, deuterium, tritium, and carbon-14. The water-quality data indicated that dissolved solids and nitrate concentrations were below regulatory limits for potable water; however, fluoride concentrations in the lower aquifer exceeded regulatory limits. Arsenic concentrations and chromium concentrations were generally below regulatory limits; however, arsenic concentrations measured in water from wells perforated in the lower aquifer exceeded regulatory limits. The carbon-14 activities ranged from 2 to 72 percent modern carbon and are consistent with uncorrected ground-water ages (time since recharge) of about 32,300 to 2,700 years before present. The oxygen-18 and deuterium composition of water sampled from the upper aquifer is similar to the volume-weighted composition of present-day winter precipitation indicating that winter precipitation was the predominant source of ground-water recharge. Field studies, conducted during water years 2001 through 2003 to determine the distribution and quantity of recharge, included installation of instrumented boreholes in selected washes and at a nearby control site. Core material and cuttings from the boreholes were analyzed for physical, chemical, and hydraulic properties. Instruments installed in the boreholes were monitored to measure changes in matric potential and temperature. Borehole data were supplemented with temperature data collected from access tubes installed at additional sites along study washes. Streambed hydraulic properties and the response of instruments to infiltration were measured using infiltrometers. Physical and geochemical data collected away from the stream channels show that direct infiltration of precipitation to depths below the root zone and subsequent gro

Scientific Investigations Report↗

Reconnaissance study of the hydrology of American Memorial Park, Island of Saipan, Commonwealth of the Northern Mariana Islands

American Memorial Park, a unit of the National Park Service on the Island of Saipan, includes among its features a 27-acre estuarine system that has become a rarity within the Commonwealth of the Northern Mariana Islands. The estuarine system's mosaic of marshy areas interspersed with emergent wetlands and mixed wet forests provides critical habitat for various migratory and resident waterfowl, including two Federally listed endangered species: the Marianas gallinule ( Gallinula chloropus guami ) and the nightingale reed warbler ( Acrocephalus luscinia ). With sensitivity to the park's ecologic assets and the uncertainty associated with locally rapid urbanization, a need to better understand the hydrology of American Memorial Park was recognized. To address that need, a reconnaissance study of the park was undertaken during August and September 2005. The goals of the study were (1) to describe the occurrence and salinity of surface and ground water within the park; (2) to develop a hydrologic model of the park area of the island, with emphasis on the 27-acre estuarine system; and (3) to identify additional data needed to further develop this model. With regard to surface water, three freshwater inputs to the park's natural wetland are possible: direct rainfall, seaward-flowing ground water, and overland flow. Direct rainfall, which is an important source of freshwater to the wetland, commonly exceeds evapotranspiration both seasonally and per storm. The seaward flow of ground water is likely to be a source of freshwater to the wetland because ground water generally has an upward vertical component in the nearshore environment. Overland flow upgradient of the park could potentially contribute a significant input of freshwater during periods of intense rainfall, but roads that flank the park's perimeter act as a barrier to surficial inflows. During the reconnaissance, four discrete bodies, or zones, of surface water were observed within the park's natural wetland. Conductivity within these surface-water zones typically ranged from 1,540 to 4,370 microsiemens per centimeter (µS/cm) at 25°C although values as low as 829 and as high as 8,750 µ S/cm were measured. As a result of these observations, the American Memorial Park wetland area meets the definition criteria of an estuarine system that is dominantly oligohaline. Conductivity was also measured in a constructed wetland that was built within the park to augment the storm-drainage infrastructure of the village of Garapan. Reverse-osmosis facilities, in operation at hotels adjacent to the park, have historically discharged highly saline wastewater into the storm-drainage system. This collective storm and wastewater flow is routed into the constructed wetland and from there into the ocean. The conductivity of water in the constructed wetland ranged from 45,000 to 62,500 µS/cm , exceeding nominal seawater values by as much as 25 percent, with the highest conductivities recorded near discharging storm drains. With regard to ground water, the reconnaissance included installation of a ground-water-monitoring network. Data collected from this network identified the presence of freshwater underlying the park and indicated that surface water is directly connected to ground water in the natural wetland because the water levels of both surface water and ground water directly varied with the tide. Conductivities of ground-water samples from wells in the monitoring network indicated that ground-water salinity was geographically related: conductivities were lower (801-2,490 (µS/cm) in surficially dry areas, intermediate (6,090-9,180 (µS/cm) in natural-wetland areas, and higher (18,250-27,700 (µS/cm) in areas adjacent to the constructed wetland and its associated ocean-discharge channel. Synoptic water-level surveys were made to enhance understanding of the spatial expression of the water table; they were scheduled to overlap with peak and trough tidal signals to enable limited characteri

Scientific Investigations Report↗

Finite-difference interblock transmissivity for unconfined aquifers and for aquifers having smoothly varying transmissivity

More accurate alternatives to the widely used harmonic mean interblock transmissivity are proposed for block-centered finite-difference models of ground-water flow in unconfined aquifers and in aquifers having smoothly varying transmissivity. The harmonic mean is the exact interblock transmissivity for steady-state one-dimensional flow with no recharge if the transmissivity is assumed to be spatially uniform over each finite-difference block, changing abruptly at the block interface. However, the harmonic mean may be inferior to other means if transmissivity varies in a continuous or smooth manner between nodes. Alternative interblock transmissivity functions are analytically derived for the case of steady-state one-dimensional flow with no recharge. The second author has previously derived the exact interblock transmissivity, the logarithmic mean, for one-dimensional flow when transmissivity is a linear function of distance in the direction of flow. We show that the logarithmic mean transmissivity is also exact for uniform flow parallel to the direction of changing transmissivity in a two- or three-dimensional model, regardless of grid orientation relative to the flow vector. For the case of horizontal flow in a homogeneous unconfined or water-table aquifer with a horizontal bottom and with areally distributed recharge, the exact interblock transmissivity is the unweighted arithmetic mean of transmissivity at the nodes. This mean also exhibits no grid-orientation effect for unidirectional flow in a two-dimensional model. For horizontal flow in an unconfined aquifer with no recharge where hydraulic conductivity is a linear function of distance in the direction of flow the exact interblock transmissivity is the product of the arithmetic mean saturated thickness and the logarithmic mean hydraulic conductivity. For several hypothetical two- and three-dimensional cases with smoothly varying transmissivity or hydraulic conductivity, the harmonic mean is shown to yield the least accurate solution to the flow equation of the alternatives considered. Application of the alternative interblock transmissivities to a regional aquifer system model indicates that the changes in computed heads and fluxes are typically small, relative to model calibration error. For this example, the use of alternative interblock transmissivities resulted in an increase in computational effort of less than 3 percent. Numerical algorithms to compute alternative interblock transmissivity functions in a modular three-dimensional flow model are presented and documented.

Water-Resources Investigations Report↗

Occurrence of Uranium and 222Radon in Glacial and Bedrock Aquifers in the Northern United States, 1993-2003

Water-quality data collected from 1,426 wells during 1993-2003 as part of the U.S. Geological Survey National Water-Quality Assessment (NAWQA) program were evaluated to characterize the water quality in glacial and bedrock aquifers of the northern United States. One of the goals of the NAWQA program is to synthesize data from individual studies across the United States to gain regional- and national-scale information about the behavior of contaminants. This study focused on the regional occurrence and distribution of uranium and 222radon in ground water in the glacial aquifer system of the United States as well as in the Cambrian-Ordovician and the New York and New England crystalline aquifer systems that underlie the glacial aquifer system. The occurrence of uranium and 222radon in ground water has long been a concern throughout the United States. In the glacial aquifers, as well as the Cambrian-Ordovician and the New York and New England crystalline aquifer systems of the United States, concentrations of uranium and 222radon were highly variable. High concentrations of uranium and 222radon affect ground water used for drinking water and for agriculture. A combination of information or data on (1) national-scale ground-water regions, (2) regional-scale glacial depositional models, (3) regional-scale geology, and (4) national-scale terrestrial gamma-ray emissions were used to confirm and(or) refine the regions used in the analysis of the water-chemistry data. Significant differences in the occurrence of uranium and 222radon, based primarily on geologic information were observed and used in this report. In general, uranium was highest in the Columbia Plateau glacial, West-Central glacial, and the New York and New England crystalline aquifer groups (75th percentile concentrations of 22.3, 7.7, and 2.9 micrograms per liter (ug/L), respectively). In the Columbia Plateau glacial and the West-Central glacial aquifer groups, more than 10 percent of wells sampled had concentrations of uranium that exceeded the U.S. Environmental Protection Agency (USEPA) Maximum Contaminant Level of 30 ug/L; in the New York and New England crystalline aquifer group, 4 percent exceeded 30 ug/L. Ground-water samples with high concentrations of uranium were commonly linked to geologic sources rich in uranium. In eight of nine aquifer groups defined for this study, concentrations of uranium correlated significantly with concentrations of sulfate in ground water (Spearman's rho = 0.20 to 0.56; p < 0.05). In the Columbia Plateau, glacial aquifers were derived in part from basaltic lava flows, some felsic volcanic rocks, and some paleo-lake bed materials that may be rich in uranium. In the Columbia Plateau and West-Central glacial aquifer groups, uranium correlated with total dissolved solids, bicarbonate, boron, lithium, selenium, and strontium. In the West-Central glacial aquifer group, rocks such as Cretaceous marine shales, which are abundant in uranium, probably contribute to the high concentrations in ground water; in the southern part of this group, which extends into Nebraska, the glacial or glacial-related sediment may be interbedded with uranium-rich materials that originated to the north and west and in the Rocky Mountains. In New England, crystalline bedrock that is granitic, such as two-mica granites, as well as other high-grade metamorphic rocks, has abundant uranium that is soluble in the predominantly oxic to sub-oxic geochemical conditions. This appears to contribute to high uranium concentrations in ground water. The highest 222radon concentrations were present in samples from wells completed in the New York and New England crystalline aquifer group; the median value (2,122 picocurries per liter (pCi/L)) was about 10 times the median values of all other aquifer groups. More than 25 percent of the samples from the New York and New England crystalline aquifer group wells had 222radon concentrations that exceeded the USEPA Alternative

Scientific Investigations Report↗

Ground-Water Temperature, Noble Gas, and Carbon Isotope Data from the Espanola Basin, New Mexico

Ground-water samples were collected from 56 locations throughout the Espanola Basin and analyzed for general chemistry (major ions and trace elements), carbon isotopes (delta 13C and 14C activity) in dissolved inorganic carbon, noble gases (He, Ne, Ar, Kr, Xe, and 3He/4He ratio), and tritium. Temperature profiles were measured at six locations in the southeastern part of the basin. Temperature profiles suggest that ground water generally becomes warmer with distance from the mountains and that most ground-water flow occurs at depths <250 m below ground surface. The two dominant water types in the basin are Ca/CO3+HCO3 and Na/CO3+HCO3, followed by mixed-cation/CO3+HCO3. Waters generally evolve from Ca/CO3+HCO3 to Na/CO3+HCO3 with increasing residence time through Ca-Na cation exchange with clay minerals. Basin ground water can be divided into four hydrochemical zones based on chemical and isotopic composition: West, Southeast, Northeast, and Central Deep. Hydrochemical zone boundaries are roughly correlated with contacts between geologic units or lithosome transitions within the Tesuque Formation. Geochemical mass-transfer modeling was performed using NETPATH and 14C ages were adjusted accordingly. Isotopic input parameters were varied within reasonable limits to assess uncertainty in the adjusted 14C ages. For each sample, a preferred adjusted age was selected from multiple possible adjusted ages based primarily on the fit between measured and modeled delta 13C values. The range of possible age adjustments for most samples is about 6,000 years or less, indicating that the preferred adjusted age for most samples has a total range of uncertainty of <6,000 years. Preferred adjusted ages range from 0 to 35,400 years. First-order trends in the age distribution include older ages generally occurring farther from rivers on the east side of the basin and farther from the mountains, consistent with both mountain-front recharge and recharge on the basin floor in the form of stream-loss and arroyo recharge. Ages also increase with depth in the Southeast zone, the only area where discrete-depth samples could be collected. Recharge temperatures derived from noble gas concentrations were used in conjunction with an empirically derived local relationship between recharge temperature and elevation to constrain recharge elevation and to estimate fractions of mountain-block recharge (MBR) in sampled waters of Holocene age. Noble gas recharge temperatures indicate that ground water in the Southeast zone contains a significant fraction of MBR, commonly 20-50 percent or more. The same is apparently true for the Northeast zone, though only two data points could be used to evaluate the MBR fraction in this area. Recharge temperatures indicate that the upper 30 m of the regional aquifer on the Pajarito Plateau typically contain little or no MBR. Tritium concentrations and apparent 3H/3He ages indicate that water in the mountain block is dominantly <50 years old, and water in the basin-fill is dominantly >50 years old, consistent with the 14C ages. Terrigenic He (Heterr) concentrations in ground water are high (log Delta Heterr of 2 to 5) throughout much of the basin. High Heterr concentrations are probably caused by in situ production in the Tesuque Formation from locally high concentrations of U-bearing minerals (Northeast zone only), or by upward diffusive/advective transport of crustal- and mantle-sourced He possibly enhanced by basement piercing faults, or by both. The 3He/4He ratio of Heterr (Rterr) is commonly high (Rterr/Ra of 0.3-2.0, where Ra is the 3He/4He ratio in air) suggesting that Espanola Basin ground water commonly contains mantle-sourced He. The 3He/4He ratio of Heterr is generally the highest in the western and southern parts of the basin, closest to the western border fault system and the Quaternary to Miocene volcanics of the Jemez Mountains and Cerros del Rio.

Scientific Investigations Report↗

Possible changes in ground-water flow to the Pecos River caused by Santa Rosa Lake, Guadalupe County, New Mexico

In 1980 Santa Rosa Dam began impounding water on the Pecos River about 7 miles north of Santa Rosa, New Mexico, to provide flood control, sediment control, and storage for irrigation. Santa Rosa Lake has caused changes in the groundwater flow system, which may cause changes in the streamflow of the Pecos River that cannot be detected at the present streamflow gaging stations. Data collected at these stations are used to measure the amount of water available for downstream users. A three-dimensional groundwater flow model for a 950 sq mi area between Anton Chico and Puerto de Luna was used to simulate the effects of Santa Rosa Lake on groundwater flow to a gaining reach of the Pecos River for lake levels of 4,675, 4,715, 4,725, 4,750, 4,776, and 4,797 feet above sea level and durations of impoundment of 30, 90, 182, and 365 days for all levels except 4 ,797 feet. These simulations indicated that streamflow in the Pecos River could increase by as much as 2 cu ft/sec between the dam and Puerto de Luna if the lake level were maintained at 4 ,797 feet for 90 days or 4,776 feet for 1 year. About 90% of this increased streamflow would occur < 0.5 mi downstream from the dam, some of which would be measured at the streamflow gaging station located 0.2 mile downstream from the dam. Simulations also indicated that the lake will affect groundwater flow such that inflow to the study area may be decreased by as much as 1.9 cu ft/sec. This water may leave the Pecos River drainage basin or be diverted back to the Pecos River downstream from the gaging station near Puerto de Luna. In either case, this quantity represents a net loss of water upstream from Puerto de Luna. Most simulations indicated that the decrease in groundwater flow into the study area would be of about the same quantity as the simulated increase in streamflow downstream from the dam. Therefore, the net effect of the lake on the flow of the Pecos River in the study area appears to be negligible. Model simulations indicated that effect of lake levels below 4 ,750 feet on water levels in observation wells completed in the San Andres Limestone could not be distinguished from the effects of other hydrologic stresses. (Author 's abstract)

Water-Resources Investigations Report↗