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,099 records · Page 61Linked to original sources
USGS2006

Sources of water to wells in updip areas of the Wenonah-Mount Laurel aquifer, Gloucester and Camden Counties, New Jersey

Since 1996, when the New Jersey Department of Environmental Protection (NJDEP) restricted ground-water withdrawals from the Potomac-Raritan-Magothy aquifer system in the southern New Jersey Coastal Plain as a result of excessive drawdown, Coastal Plain communities have been interested in developing alternate sources of water supply for their residents. The use of ground water from areas near the updip parts of the overlying confined aquifers where withdrawals are not restricted is being considered to meet the demand for drinking water. Concerns have arisen, however, regarding the potential effects of increased withdrawals from these areas on ground-water flow to streams and wetlands as well as to the deeper, confined parts of the aquifers. Therefore, the U.S. Geological Survey, in cooperation with the NJDEP, conducted a study to investigate the sources of water to currently inactive wells in the updip part of the Wenonah-Mount Laurel aquifer in Gloucester and Camden Counties, New Jersey. Of particular interest is whether the primary source of the increased withdrawals is likely to be the aquifer outcrop or the downdip, confined part of the aquifer. The outcrop of the Wenonah-Mount Laurel aquifer covers nearly 8 mi2 (square miles), or about 46 percent of Deptford Township's 17.56-mi2 area. The Deptford Township Municipal Utilities Authority owns six currently (2005) inactive wells in the Wenonah-Mount Laurel aquifer at the southeastern boundary of Deptford Township, 1.25 mi (miles) from the outcrop. For the purposes of this study, an existing ground-water-flow model of the New Jersey Coastal Plain aquifers was used to simulate ground-water-flow conditions in Gloucester and Camden Counties in 1998. Two alternative withdrawal scenarios were superimposed on the results of the 1998 simulation. In the first (the 'full-allocation' scenario), full-allocation withdrawal rates established by the NJDEP were applied to 45 existing wells in the Deptford Township area. In the second (the 'additional-withdrawal' scenario), the full-allocation scenario was modified by adding an additional withdrawal of 1.62 million gallons per day from the six inactive Deptford Township withdrawal wells. Simulated drawdown for the full-allocation scenario is zero to near zero in Deptford Township. Changes are greatest downdip from Deptford Township, where a broad area of 5- to 10-ft (feet) drawdowns is simulated; maximum drawdown at the center of the cone of depression is 20 ft. Water levels declined as much as 10 ft around individual wells whose current withdrawals are only a small percentage of their allotted allocation. Simulated drawdown for the additional-withdrawal scenario exceeds 40 ft and is centered around the six inactive Deptford Township withdrawal wells. The area in which the simulated drawdown is 5 ft extends approximately 3.75 mi downdip from the wells and 2 mi updip, into the outcrop. Water budgets based on the simulation results for the full-allocation and additional-withdrawal scenarios were calculated and compared, with particular focus on a 75-mi2 area in and around Deptford Township that includes the outcrop of the Wenonah-Mount Laurel aquifer and part of the area downdip from the outcrop (budget zone 2). The comparison of the two water budgets for zone 2 shows that 46 percent of the withdrawals from the six inactive Deptford Township wells would result from reduced stream base flow in the outcrop of the Wenonah-Mount Laurel aquifer and 35 percent would result from increased downward flow from the overlying Vincentown aquifer. Four percent would result from increased flow from the downdip areas of the Wenonah-Mount Laurel aquifer, 5 percent would result from decreased flow to the downdip areas of the Wenonah-Mount Laurel aquifer, and 5 percent would result from decreased flow to the underlying Englishtown aquifer system. The remaining 4 percent was attributed to decreased upward flow to the overlying Vincentown aquifer.

Scientific Investigations Report↗
USGS2002

Estimates of ground-water recharge from precipitation to glacial-deposit and bedrock aquifers on Lopez, San Juan, Orcas, and Shaw islands, San Juan County, Washington

An important source of fresh water on Lopez, San Juan, Orcas, and Shaw Islands in San Juan County off the northwestern coast of Washington is glacial-deposit and bedrock aquifers. Two methods were used to estimate recharge from precipitation to the water tables on the islands. A daily near-surface water-balance method, the Deep Percolation Model (DPM), was used to simulate water budgets for the period October 1, 1996, through September 30, 1998 (water years 1997-98) for six small drainage basins—three on Lopez Island and one each on San Juan, Orcas, and Shaw Islands. The calibrated soil and subsoil parameters from the DPM for each small basin were then used in island-wide applications of the DPM where the direct runoff component (which is not available on an island-wide basis) was simulated, rather than input, and calibration was not required. A spatial distribution of annual recharge was simulated for each island, with island averages of: Lopez Island, 2.49 inches per year; San Juan Island, 1.99 inches per year; Orcas Island, 1.46 inches per year; and Shaw Island, 1.44 inches per year. A chloride mass-balance method that requires measurements of atmospheric chloride deposition, precipitation, streamflow, and chloride concentrations in ground water was used to estimate recharge to the glacial-deposit aquifers of Lopez Island. Only average recharge could be estimated using this method rather than area-specific recharge. Average recharge for Lopez Island estimated by this method was only 0.63 inch per year. The range of chloride concentrations in ground-water samples from selected wells indicates that the average recharge in areas of glacial deposits is between 0.29 and 1.95 inches per year. Recharge simulated using the DPM for two drainage basins on Lopez Island overlain by glacial deposits are 2.76 and 2.64 inches per year. Sources of chloride in ground water other than from the atmosphere would cause the recharge estimated by the chloride mass-balance method to be less than the actual recharge, therefore these estimates may represent lower limits which are, at least, consistent with the higher simulated recharge from the DPM. The average island-wide recharge is most closely related to the amount of area overlain by glacial deposits. Thus, even though Lopez Island receives the least precipitation, it has the most recharge per square mile because it proportionally has the largest area overlain by glacial deposits. Recharge simulated by the DPM for areas of shallow to outcropping bedrock generally were less than 1.5 inches per year, but recharge simulated in areas of glacial deposits ranged from less than 0.5 to 3 inches per year, with recharge as high as 9 inches per year in some small areas.

Washington↗
USGS2005

Simulation of ground-water flow in the basin-fill aquifer of the Tularosa Basin, south-central New Mexico, predevelopment through 2040

The hydrology of the basin-fill aquifer in the Tularosa Basin was evaluated through construction and calibration of steady-state and transient three-dimensional ground-water-flow simulations. Simulations were made using the U.S. Geological Survey finite-difference modular ground-water-flow computer software MODFLOW-96. The transient simulation covered 1948-2040. Both steady-state and transient simulations were calibrated by matching simulation output to available ground-water-level measurements. The root-mean-square error of the steady-state calibration in the well-calibrated area of the ground-water-flow simulation was 6.3 meters, and root-mean-square errors of individual transient-calibration points ranged from 0.8 to 17.0 meters. The areal distribution of water-level measurements used in the steady-state and transient calibrations restricts the well-calibrated area of the model to the eastern side of the Tularosa Basin. Water levels in the La Luz Creek subbasin area were underestimated by both the steady-state and transient models, suggesting that the hydrology of this area is not well represented in the model. About 143,000 cubic meters per day of recharge is estimated to enter the basin-fill aquifer from subbasins that rim the Tularosa Basin. The estimated recharge is about 4-5 percent of total precipitation in most subbasins. Approximately 88 percent of total recharge left the basin-fill aquifer as evapotranspiration under predevelopment conditions. Water levels were simulated for 1948, 1995, and 2040 under scenarios of zero and maximum return flows. Estimated return flows from municipalities were calculated on the basis of data in the Tularosa Basin Regional Water Plan for 2000-2040. Agricultural return flows were estimated primarily on the basis of ground-water-withdrawal, ground-water-depletion, surface-water-withdrawal, and surface-water-depletion data for the Tularosa Basin. The ground-water-flow simulation was sensitive to the return-flow scenario in the agricultural area near Tularosa and decreasingly sensitive to the south. Declines in simulated water levels near Tularosa between 1948 and 1995 were as large as 30 meters under the zero return-flow scenario and 15 meters under the maximum return-flow scenario. Declines in simulated water levels between 1995 and 2040 were as large as 25 meters under the zero return-flow scenario and 15 meters under the maximum return-flow scenario. Comparison of water levels measured near Tularosa in 1991 and water levels simulated under the maximum return-flow scenario for 1991 suggests that declines in simulated water levels near Tularosa may be overestimated under the zero return-flow scenario. Declines in simulated water levels near the City of Alamogordo well field between 1948 and 1995 were as large as 15 meters under the zero return-flow scenario and 10 meters under the maximum return-flow scenario. Simulated declines in water levels between 1995 and 2040 were nearly 15 meters under both return-flow scenarios assuming that all projected increases in withdrawal came from existing City of Alamogordo public-supply wells and all withdrawal from the wells came from the basin-fill aquifer. Declines in simulated water levels near the Holloman Air Force Base well fields between 1948 and 1995 and between 1995 and 2040 were less than 5 meters under both the zero and maximum return-flow scenarios. In 1995 under the zero return-flow scenario, an estimated 56,000 cubic meters of water per day was removed from aquifer storage. Of the approximately 199,000 cubic meters of water per day that left the aquifer under 1995 conditions, 40 percent left the basin-fill aquifer as ground-water withdrawal, 51 percent as evapotranspiration, 7 percent by interbasin ground-water flow into the Hueco Bolson, and 2 percent by flow into creeks and springs. Generalized directions of ground-water flow were simulated for 1948, 1995, and 2040 for much of the eastern part of the Tularosa Basin. Localized

Scientific Investigations Report↗
USGS2001

Simulation of a long-term aquifer test conducted near the Rio Grande, Albuquerque, New Mexico

A long-term aquifer test was conducted near the Rio Grande in Albuquerque during January and February 1995 using 22 wells and piezometers at nine sites, with the City of Albuquerque Griegos 1 production well as the pumped well. Griegos 1 discharge averaged about 2,330 gallons per minute for 54.4 days. A three-dimensional finite-difference ground-water-flow model was used to estimate aquifer properties in the vicinity of the Griegos well field and the amount of infiltration induced into the aquifer system from the Rio Grande and riverside drains as a result of pumping during the test. The model was initially calibrated by trial-and-error adjustments of the aquifer properties. The model was recalibrated using a nonlinear least-squares regression technique. The aquifer system in the area includes the middle Tertiary to Quaternary Santa Fe Group and post-Santa Fe Group valley- and basin-fill deposits of the Albuquerque Basin. The Rio Grande and adjacent riverside drains are in hydraulic connection with the aquifer system. The hydraulic-conductivity values of the upper part of the Santa Fe Group resulting from the model calibrated by trial and error varied by zone in the model and ranged from 12 to 33 feet per day. The hydraulic conductivity of the inner-valley alluvium was 45 feet per day. The vertical to horizontal anisotropy ratio was 1:140. Specific storage was 4 x 10-6 per foot of aquifer thickness, and specific yield was 0.15 (dimensionless). The sum of squared errors between the observed and simulated drawdowns was 130 feet squared. Not all aquifer properties could be estimated using nonlinear regression because of model insensitivity to some aquifer properties at observation locations. Hydraulic conductivity of the inner-valley alluvium, middle part of the Santa Fe Group, and riverbed and riverside-drain bed and specific yield had low sensitivity values and therefore could not be estimated. Of the properties estimated, hydraulic conductivity of the upper part of the Santa Fe Group was estimated to be 12 feet per day, the vertical to horizontal anisotropy ratio was estimated to be 1:82, and specific storage was estimated to be 1.2 x 10-6 per foot of aquifer thickness. The overall sum of squared errors between the observed and simulated drawdowns was 87 feet squared, a significant improvement over the model calibrated by trial and error. At the end of aquifer-test pumping, induced infiltration from the Rio Grande and riverside drains was simulated to be 13 percent of the total amount of water pumped. The remainder was water removed from aquifer storage. After pumping stopped, induced infiltration continued to replenish aquifer storage. Simulations estimated that 5 years after pumping began (about 4.85 years after pumping stopped), 58 to 72 percent of the total amount of water pumped was replenished by induced infiltration from the Rio Grande surface-water system.

Water-Resources Investigations Report↗
USGS2018

Hydrogeologic controls and geochemical indicators of groundwater movement in the Niles Cone and southern East Bay Plain groundwater subbasins, Alameda County, California

Beginning in the 1970s, Alameda County Water District began infiltrating imported water through ponds in repurposed gravel quarries at the Quarry Lakes Regional Park, in the Niles Cone groundwater subbasin, to recharge groundwater and to minimize intrusion of saline, San Francisco Bay water into freshwater aquifers. Hydraulic connection between distinct aquifers underlying Quarry Lakes allows water to recharge the upper aquifer system to depths of 400 feet below land surface, and the Deep aquifer to depths of more than 650 feet. Previous studies of the Niles Cone and southern East Bay Plain groundwater subbasins suggested that these two subbasins may be hydraulically connected. Characterization of storage capacities and hydraulic properties of the complex aquifers and the structural and stratigraphic controls on groundwater movement aids in optimal storage and recovery of recharged water and provides information on the ability of aquifers shared by different water management agencies to fulfill competing storage and extraction demands. The movement of recharge water through the Niles Cone groundwater subbasin from Quarry Lakes and the possible hydraulic connection between the Niles Cone and the southern East Bay Plain groundwater subbasins were investigated using interferometric synthetic aperture radar (InSAR), water-chemistry, and isotopic data, including tritium/helium-3, helium-4, and carbon-14 age-dating techniques. InSAR data collected during refilling of the Quarry Lakes recharge ponds show corresponding ground-surface displacement. Maximum uplift was about 0.8 inches, reasonable for elastic expansion of sedimentary materials experiencing an increase in hydraulic head that resulted from pond refilling. Sodium concentrations increase while calcium and magnesium concentrations in groundwater decrease along groundwater flowpaths from the Niles Cone groundwater subbasin through the Deep aquifer to the northwest toward the southern East Bay Plain groundwater subbasin. Residual effects of pre-1970s intrusion of saline water from San Francisco Bay, including high chloride concentrations in groundwater, are evident in parts of the Niles Cone subbasin. Noble gas recharge temperatures indicate two primary recharge sources (Quarry Lakes and Alameda Creek) in the Niles Cone groundwater subbasin. Although recharge at Quarry Lakes affects hydraulic heads as far as the transition zone between the Niles Cone and East Bay Plain groundwater subbasins (about 5 miles), the effect of recharged water on water quality is only apparent in wells near (less than 2 miles) recharge sources. Groundwater chemistry from upper aquifer system wells near Quarry Lakes showed an evaporated signal (less negative oxygen and hydrogen isotopic values) relative to surrounding groundwater and a tritium concentration (2 tritium units) consistent with recently recharged water from a surface-water impoundment. Uncorrected carbon-14 activities measured in water sampled from wells in the Niles Cone groundwater subbasin range from 16 to 100 percent modern carbon (pmC). The geochemical reaction modeling software NETPATH was used to interpret carbon-14 ages along a flowpath from Quarry Lakes toward the East Bay Plain groundwater subbasin. Model results indicate that changes in groundwater chemistry are controlled by cation exchange on clay minerals and weathering of primary silicate minerals. Old groundwater (lower carbon-14 activities) is characterized by high dissolved silica and pH. Interpreted carbon-14 ages ranged from 830 to more than 7,000 years before present and are less than helium-4 ages that range from 2,000 to greater than 11,000 years before present. The average horizontal groundwater velocity along the studied flowpath, as calculated using interpreted carbon-14 ages, through the Deep aquifer of the Niles Cone groundwater subbasin is between 3 and 12 feet per year. The groundwater velocity decreases near the boundary of the transition zone to the southern East Bay Plain groundwater subbasin to about 0.5 feet per year. These changes may result from water recharged from different sources converging in flowpaths north of the transition zone, or a boundary to flow between the Niles Cone and southern East Bay Plain groundwater subbasins, likely owing to changes in lithology caused by depositional patterns.

California↗
USGS2002

Hydrogeology and simulation of ground-water flow and land-surface subsidence in the Chicot and Evangeline aquifers, Houston area, Texas

In November 1997, the U.S. Geological Survey, in cooperation with the City of Houston Utilities Planning Section and the City of Houston Department of Public Works & Engineering, began an investigation of the Chicot and Evangeline aquifers in the greater Houston area in Texas to better understand the hydrology, flow, and associated land-surface subsidence. The principal part of the investigation was a numerical finite-difference model (MODFLOW) developed to simulate ground-water flow and land-surface subsidence in an 18,100-square-mile area encompassing greater Houston. The focus of the study was Harris and Galveston Counties, but other counties were included to achieve the appropriate boundary conditions. The model was vertically discretized into three 103-row by 109-column layers resulting in a total of 33,681 grid cells. Layer 1 represents the water table using a specified head, layer 2 represents the Chicot aquifer, and layer 3 represents the Evangeline aquifer. Simulations were made under transient conditions for 31 ground-water-withdrawal (stress) periods spanning 1891–1996. The years 1977 and 1996 were chosen as potentiometric-surface calibration periods for the model. Simulated and measured potentiometric surfaces of the Chicot and Evangeline aquifers for 1977 match closely. Waterlevel measurements indicate that by 1977, large ground-water withdrawals in east-central and southeastern areas of Harris County had caused the potentiometric surfaces to decline as much as 250 feet below sea level in the Chicot aquifer and as much as 350 feet below sea level in the Evangeline aquifer. Simulated and measured potentiometric surfaces of the Chicot and Evangeline aquifers for 1996 also match closely. The large potentiometric-surface decline in 1977 in the southeastern Houston area showed significant recovery by 1996. The 1996 centers of potentiometric-surface decline are located much farther northwest. Potentiometric-surface declines of more than 200 feet below sea level in the Chicot aquifer and more than 350 feet below sea level in the Evangeline aquifer were measured in observation wells and simulated in the flow model. Simulation of land-surface subsidence and water released from storage in the clay layers was accomplished using the Interbed-Storage Package of the MODFLOW model. Land-surface subsidence was calibrated by comparing simulated long-term (1891–1995) and short-term (1978–95) land-surface subsidence with published maps of land-surface subsidence for about the same period until acceptable matches were achieved. Simulated 1996 Chicot aquifer flow rates indicate that a net flow of 562.5 cubic feet per second enters the Chicot aquifer in the outcrop area, and a net flow of 459.5 cubic feet per second passes through the Chicot aquifer into the Evangeline aquifer. The remaining 103.0 cubic feet per second of flow is withdrawn as pumpage, with a shortfall of about 84.9 cubic feet per second supplied to the wells from storage in sands and clays. Water simulated from storage in clays in the Chicot aquifer is about 19 percent of the total water withdrawn from the aquifer. Simulated 1996 Evangeline aquifer flow rates indicate that a net flow of 14.8 cubic feet per second enters the Evangeline aquifer in the outcrop area, and a net flow of 459.5 cubic feet per second passes through the Chicot aquifer into the Evangeline aquifer for a total inflow of 474.3 cubic feet per second. A greater amount, 528.6 cubic feet per second, is withdrawn by wells; the shortfall of about 54.8 cubic feet per second is supplied from storage in sands and clays. Water simulated from storage in clays in the Evangeline aquifer is about 10 percent of the total water withdrawn from the aquifer.

Texas↗
USGS1998

Areas contributing ground water to the Peconic Estuary, and ground-water budgets for the north and south forks and Shelter Island, eastern Suffolk County, New York

The Peconic Estuary, at the eastern end of Long Island, has been plagued by a recurrent algal bloom, locally referred to as ?Brown Tide,? that has caused the severe decline of local marine resources. Although the factors that trigger Brown Tide blooms remain uncertain, groundwater discharge has previously been shown to affect surface-water quality in the western part of the estuary. A U.S. Geological Survey groundwater- flow model of the main body of Long Island indicates that a total of about 7.5 x 106 ft3/d (cubic feet per day) of freshwater discharges to the western part of the estuary, but the model does not include the ground-water flow systems on the North and South Forks and Shelter Island, which contribute significant amounts of freshwater to the central and eastern parts of the estuary. The need for information on freshwater discharge to the entire estuary prompted the U.S. Geological Survey to evaluate ground-water discharge from the North and South Forks and Shelter Island. Source areas that contribute ground water to the Peconic Estuary were delineated, and groundwater budgets for these areas were developed, to evaluate the distribution and magnitude of ground-water discharge to the central and eastern parts of the estuary. Contributing-area boundaries that were delineated coincide with the hydraulic boundaries of the fresh ground-water-flow systems of the North and South Forks and Shelter Island; these boundaries are of two types? external (saltwater bodies) and internal (groundwater divides). Hydrologic components that were evaluated include recharge from precipitation, public-supply withdrawal and return flow, and agricultural withdrawal. Values for each of these components were calculated or estimated for the individual freshwater flow subsystems that form each ground-water-budget area, then summed to obtain the total discharge of fresh ground water to tidewater. Ground-water discharge to the Peconic Estuary is about 3.8 x 106 ft3/d from the North Fork, 11 x 106 ft3/d from the South Fork, and 1.7 x 106 ft3/d from Shelter Island. The total contribution to the estuary from these areas is about 16 x 106 ft3/d?roughly twice the total contribution from the main body of Long Island. In contrast to the freshwater contribution from the main body of Long Island, which is concentrated near the head of the estuary, the contributions from the North and South Forks and Shelter Island are distributed along the east-west length of the estuary. Changes in water-table altitude and the resulting changes in total discharge to the Peconic Estuary were estimated from the relative changes in annual mean water level at observation wells. The 1985-95 interval included 7 years (1985-88, 1991- 92, 1995) of generally below-average water-table altitudes that presumably caused similar decreases in ground-water discharge to the estuary; intense Brown Tide blooms coincided with six of these years (1985-88, 1991, 1995), and localized blooms coincided with the remaining year (1992). Watertable altitudes in the remaining 4 years of the 1985-95 interval (1989-90, 1993-94) were nearly average or above average, and presumably produced comparably near-average or increased amounts of ground-water discharge to the estuary; none of these years saw any widespread Brown Tide blooms. Fluctuations in the amounts of ground-water discharge to the estuary appear to affect the occurrence of Brown Tide blooms, although the factors that trigger the blooms have not been determined.

Water-Resources Investigations Report↗
USGS2006

Simulation of saltwater movement in the Upper Floridan aquifer in the Savannah, Georgia-Hilton Head Island, South Carolina, area, predevelopment-2004, and projected movement for 2000 pumping conditions

A digital model was developed to simulate ground-water flow and solute transport for the Upper Floridan aquifer in the Savannah, Georgia–Hilton Head Island, South Carolina, area. The model was used to (1) simulate trends of saltwater intrusion from predevelopment to the present day (1885–2004), (2) project these trends from the present day into the future, and (3) evaluate the relative influence of different assumptions regarding initial and boundary conditions and physical properties. The model is based on a regional, single-density ground-water flow model of coastal Georgia and adjacent parts of South Carolina and Florida. Variable-density ground-water flow and solute transport were simulated using the U.S. Geological Survey finite-element, variable-density solute-transport simulator SUTRA, 1885–2004. The model comprises seven layers: the surficial aquifer system, the Brunswick aquifer system, the Upper Floridan aquifer, the Lower Floridan aquifer, and the intervening confining units. The model was calibrated to September 1998 water levels, for single-density freshwater conditions, then refined using variable density and chloride concentration to give a reasonable match to the trend in the chloride distribution in the Upper Floridan aquifer inferred from field measurements of specific conductance made during 2000, 2002, 2003, and 2004. The model was modified to simulate solute transport by allowing saltwater to enter the system through localized areas near the northern end of Hilton Head Island, at Pinckney Island, and near the Colleton River, and was calibrated to match chloride concentrations inferred from field measurements of specific conductance. This simulation is called the "Base Case." Water-level residuals ranged from –5.3 to 23.4 feet for September 1998 conditions and single-density freshwater conditions. When chloride transport was simulated, water-level residuals ranged from –12.5 to 23.3 feet. The simulated chloride distribution captures the general trends in the field data. Chloride transport is sensitive to the permeabilities assigned to the confining units in the source areas and the porosity assigned to the Upper Floridan aquifer. Results of the study indicate that if present-day (year 2000) pumping conditions are maintained, plumes of saltwater in the Upper Floridan aquifer will continue to expand and move toward Savannah and across Hilton Head Island; the rate of movement of the 250-mg/L (milligram per liter) isochlor toward Savannah is between 144 feet per year and 190 feet per year and that the 250-mg/L isochlor could reach the pumping center at Savannah in 800 years; if effective porosities are lower than those used in the model, as is likely, higher rates of solute transport would result; and plumes may have occurred along the northern shore of Hilton Head Island before substantial development began in the mid-1960s, and lesser amounts of intrusion may have already occurred prior to the onset of pumping during 1885. Model limitations include uncertainty in (1) field data, (2) the conceptual model, (3) the physical properties and representation of the hydrogeologic framework, and (4) uncertainty in the boundary and initial conditions. Results of simulations projected far into the future must be interpreted with caution because they are based on an assumed future pumping distribution and fixed boundary conditions, and because these conditions may differ substantially from those for which the model is calibrated.

South Carolina↗
USGS2007

Vertical gradients in water chemistry and age in the Northern High Plains Aquifer, Nebraska, 2003

The northern High Plains aquifer is the primary source of water used for domestic, industrial, and irrigation purposes in parts of Colorado, Kansas, Nebraska, South Dakota, and Wyoming. Despite the aquifer’s importance to the regional economy, fundamental ground-water characteristics, such as vertical gradients in water chemistry and age, remain poorly defined. As part of the U.S. Geological Survey’s National Water-Quality Assessment Program, water samples from nested, short-screen monitoring wells installed in the northern High Plains aquifer were analyzed for major ions, nutrients, trace elements, dissolved organic carbon, pesticides, stable and radioactive isotopes, dissolved gases, and other parameters to evaluate vertical gradients in water chemistry and age in the aquifer. Chemical data and tritium and radiocarbon ages show that water in the aquifer was chemically and temporally stratified in the study area, with a relatively thin zone of recently recharged water (less than 50 years) near the water table overlying a thicker zone of older water (1,800 to 15,600 radiocarbon years). In areas where irrigated agriculture was an important land use, the recently recharged ground water was characterized by elevated concentrations of major ions and nitrate and the detection of pesticide compounds. Below the zone of agricultural influence, major-ion concentrations exhibited small increases with depth and distance along flow paths because of rock/water interactions. The concentration increases were accounted for primarily by dissolved calcium, sodium, bicarbonate, sulfate, and silica. In general, the chemistry of ground water throughout the aquifer was of high quality. None of the approximately 90 chemical constituents analyzed in each sample exceeded primary drinking-water standards. Mass-balance models indicate that changes in groundwater chemistry along flow paths in the aquifer can be accounted for by small amounts of feldspar and calcite dissolution; goethite and clay-mineral precipitation; organic-carbon and pyrite oxidation; oxygen reduction and denitrification; and cation exchange. Mixing with surface water affected the chemistry of ground water in alluvial sediments of the Platte River Valley. Radiocarbon ages in the aquifer, adjusted for carbon mass transfers, ranged from 1,800 to 15,600 14C years before present. These results have important implications with respect to development of ground-water resources in the Sand Hills. Most of the water in the aquifer predates modern anthropogenic activity so excessive removal of water by pumping is not likely to be replenished by natural recharge in a meaningful timeframe. Vertical gradients in ground-water age were used to estimate long-term average recharge rates in the aquifer. In most areas, the recharge rates ranged from 0.02 to 0.05 foot per year. The recharge rate was 0.2 foot per year in one part of the aquifer characterized by large downward hydraulic gradients. Nitrite plus nitrate concentrations at the water table were 0.13 to 3.13 milligrams per liter as nitrogen, and concentrations substantially decreased with depth in the aquifer. Dissolved-gas and nitrogen-isotope data indicate that denitrification in the aquifer removed 0 to 97 percent (average = 50 percent) of the nitrate originally present in recharge. The average amount of nitrate removed by denitrification in the aquifer north of the Platte River (Sand Hills) was substantially greater than the amount removed south of the river (66 as opposed to 0 percent), and the extent of nitrate removal appears to be related to the presence of thick deposits of sediment on top of the Ogallala Group in the Sand Hills that contained electron donors, such as organic carbon and pyrite, to support denitrification. Apparent rates of dissolved-oxygen reduction and denitrification were estimated on the basis of decreases in dissolved-oxygen concentrations and increases in concentrations of excess nitrogen gas and ground-water ages along flow paths from the water table to deeper wells. Median rates of dissolved-oxygen reduction and denitrification south of the Platte River were at least 10 times smaller than the median rates north of the river in the Sand Hills. The relatively large denitrification rates in the Sand Hills indicate that the aquifer in that area may have a greater capacity to attenuate nitrate contamination than the aquifer south of the river, depending on rates of ground-water movement in the two areas. Small denitrification rates south of the river indicate that nitrate contamination in that part of the aquifer would likely persist for a longer period of time.

Nebraska↗
USGS2007

Geophysical Methods for Investigating Ground-Water Recharge

While numerical modeling has revolutionized our understanding of basin-scale hydrologic processes, such models rely almost exclusively on traditional measurements?rainfall, streamflow, and water-table elevations?for calibration and testing. Model calibration provides initial estimates of ground-water recharge. Calibrated models are important yet crude tools for addressing questions about the spatial and temporal distribution of recharge. An inverse approach to recharge estimation is taken of necessity, due to inherent difficulties in making direct measurements of flow across the water table. Difficulties arise because recharging fluxes are typically small, even in humid regions, and because the location of the water table changes with time. Deep water tables in arid and semiarid regions make recharge monitoring especially difficult. Nevertheless, recharge monitoring must advance in order to improve assessments of ground-water recharge. Improved characterization of basin-scale recharge is critical for informed water-resources management. Difficulties in directly measuring recharge have prompted many efforts to develop indirect methods. The mass-balance approach of estimating recharge as the residual of generally much larger terms has persisted despite the use of increasing complex and finely gridded large-scale hydrologic models. Geophysical data pertaining to recharge rates, timing, and patterns have the potential to substantially improve modeling efforts by providing information on boundary conditions, by constraining model inputs, by testing simplifying assumptions, and by identifying the spatial and temporal resolutions needed to predict recharge to a specified tolerance in space and in time. Moreover, under certain conditions, geophysical measurements can yield direct estimates of recharge rates or changes in water storage, largely eliminating the need for indirect measures of recharge. This appendix presents an overview of physically based, geophysical methods that are currently available or under development for recharge monitoring. The material is written primarily for hydrogeologists. Uses of geophysical methods for improving recharge monitoring are explored through brief discussions and case studies. The intent is to indicate how geophysical methods can be used effectively in studying recharge processes and quantifying recharge. As such, the material constructs a framework for matching the strengths of individual geophysical methods with the manners in which they can be applied for hydrologic analyses. The appendix is organized in three sections. First, the key hydrologic parameters necessary to determine the rate, timing, and patterns of recharge are identified. Second, the basic operating principals of the relevant geophysical methods are discussed. Methods are grouped by the physical property that they measure directly. Each measured property is related to one or more of the key hydrologic properties for recharge monitoring. Third, the emerging conceptual framework for applying geophysics to recharge monitoring is presented. Examples of the application of selected geophysical methods to recharge monitoring are presented in nine case studies. These studies illustrate hydrogeophysical applications under a wide range of conditions and measurement scales, which vary from tenths of a meter to hundreds of meters. The case studies include practice-proven as well as emerging applications of geophysical methods to recharge monitoring.

Professional Paper↗
USGS2007

Chemical Characteristics, Water Sources and Pathways, and Age Distribution of Ground Water in the Contributing Recharge Area of a Public-Supply Well near Tampa, Florida, 2002-05

In 2001, the National Water-Quality Assessment (NAWQA) Program of the U.S. Geological Survey began a series of studies on the transport of anthropogenic and natural contaminants (TANC) to public-supply wells. The main goal of the TANC program was to better understand the source, transport, and receptor factors that control contaminant movement to public-supply wells in representative aquifers of the United States. Studies were first conducted at regional scales at four of the eight TANC study areas during 2002-03 and at small (local) scales during 2003-05 in California, Nebraska, Connecticut, and Florida. In the Temple Terrace study area near Tampa, Florida, multiple chemical indicators and geochemical and ground-water flow modeling techniques were used to assess the vulnerability of a public-supply well in the karstic Upper Floridan aquifer to contamination from anthropogenic and naturally occurring contaminants. During 2003-05, water samples were collected from the public-supply well and 13 surrounding monitoring wells that all tap the Upper Floridan aquifer, and from 15 monitoring wells in the overlying surficial aquifer system and the intermediate confining unit that are located within the modeled ground-water contributing recharge area of the public-supply well. Six volatile organic compounds and four pesticides were detected in trace concentrations (well below drinking-water standards) in water from the public-supply well, which had an open interval from 36 to 53 meters below land surface. These contaminants were detected more frequently in water samples from monitoring wells in the overlying clastic surficial aquifer system than in water from monitoring wells in the Upper Floridan aquifer in the study area. Likewise, nitrate-N concentrations in the public-supply well (0.72-1.4 milligrams per liter) were more similar to median concentrations in the oxic surficial aquifer system (2.1 milligrams per liter) than to median nitrate-N concentrations in the anoxic Upper Floridan aquifer (0.06 milligram per liter) under sulfate-reducing conditions. High concentrations of radon-222 and uranium in the public-supply well compared to those in monitoring wells in the Upper Floridan aquifer appear to originate from water moving downward through sands and discontinuous clay lenses that overlie the aquifer. Water samples also were collected from three overlapping depth intervals (38-53, 43-53, and 49-53 meters below land surface) in the public-supply well. The 49- to 53-meter interval was identified as a high-flow zone during geophysical logging of the wellbore. Water samples were collected from these depth intervals at a low pumping rate by placing a low-capacity submersible pump (less than 0.02 cubic meter per minute) at the top of each interval. To represent higher pumping conditions, a large-capacity portable submersible pump (1.6 cubic meters per minute) was placed near the top of the open interval; water-chemistry samples were collected using the low-capacity submersible pump. The 49- to 53-meter depth interval had distinctly different chemistry than the other two sampled intervals. Higher concentrations of nitrate-N, atrazine, radon, trichloromethane (chloroform), and arsenic (and high arsenic (V)/arsenic (III) ratios); lower concentrations of dissolved solids, strontium, iron, manganese, and lower nitrogen and sulfur isotope ratios were found in this highly transmissive zone in the limestone than in water from the two other depth intervals. Movement of water likely occurs from the overlying sands and clays of the oxic surficial aquifer system and intermediate confining unit (that contains high radon-222 and nitrate-N concentrations) into the anoxic Upper Floridan aquifer (that contains low radon-222 and nitrate-N concentrations). Differences in arsenic concentrations in water from the various depth intervals in the public-supply well (3.2-19.0 micrograms per liter) were related to pumping conditions. The high arsenic

Scientific Investigations Report↗
USGS1991

Hydrogeology and water quality of glacial-drift aquifers in the Bemidji-Bagley area, Beltrami, Clearwater, Cass, and Hubbard Counties, Minnesota

Unconfined and the upper confined aquifers in glacial drift are the primary sources of water in a 1,600 square-mile area including parts of Beltrami, Cass, Clearwater, and Hubbard Counties, Minnesota. The unconfineddrift aquifer consists of coarse sand and gravel in the center of the study area. The total area underlain by the unconfined-drift aquifer is approximately 550 square miles. The unconfined aquifer ranges in thickness from 0 to 130 feet, and is greater than 20 feet thick over an area of 280 square miles. On the basis of scant data, the transmissivity of the unconfined aquifer ranges from less than 70 feet squared per day in the south and west to greater than 8,900 feet squared per day in an area west of Bemidji. Well yields from 10 to 300 gallons per minute are possible in some areas. The unconfined and upper confined-drift aquifers are separated by a fine-grained confining unit of till or lake deposits. The thickness of the upper confined-drift aquifer ranges from 0 to 60 feet in the Bemidji area. On the basis of specific-capacity and aquifer-thickness data, and results of model simulations, the transmissivity of the upper confined-drift aquifer ranges from less than 100 feet squared per day in the south and west parts of the aquifer to about 12,800 feet squared per day in the area around Bemidji. Well yields of 10 to 2,100 gallons per minute are possible in some areas. The direction of ground-water flow in both unconfined and upper confineddrift aquifers is toward the Mississippi and Clearwater Rivers. These rivers are the major discharge points for both aquifers. Ground-water divides, which separate the ground-water flow systems that discharge to these rivers, are approximately coincidental with surface-water divides between the rivers. Water from both aquifers generally is of the calcium bicarbonate type and is very hard, averaging 309 and 267 milligrams per liter as CaCO 3 from confined and unconfined-drift aquifers, respectively. Water from both aquifers generally is suitable for drinking, crop irrigation, and most other uses. Concentrations of ammonia, boron, chromium, iron, manganese, and phenols, however, locally exceed recommended limits for drinking water (Minnesota Pollution Control Agency, 1988). Longer residence time and leakage through glacial till is believed to cause higher concentrations of common inorganic constituents in water from confined-drift aquifers than concentrations in water from the unconfined-drift aquifer. Elevated concentrations of nutrients, chloride, and phenols in the unconfined-drift aquifer may be related to land-use practices. Statistical comparisons of common chemical constituents in water from wells completed in the unconfined-drift aquifer in several land-use areas suggest that concentrations of many constituents and physical properties are generally greater for wells in areas of commercial and residential land-use than for wells in areas of agriculture or forest land-use. These constituents include ammonia plus organic nitrogen, phosphorus, calcium, sodium, potassium, chloride, sulfate, silica, dissolved solids, and specific conductances. The mean values of ammonia nitrogen, magnesium, and fluoride are generally greater for wells in commercial land-use type areas than for wells in forested and agricultural land-use type areas. The mean concentration of nitrogen (N0 2 + NO 3 , dissolved) is generally greater for wells in residential land-use type areas than for wells in forested and agricultural land-use type areas. The Kruskil-Wallis test, a nonparametric that for 12 of the 21 constituents sampled in groups in the unconfined-drift aquifer, a of these constituents and land use was found statistical technique, indicated common in all land-use type relation between the concentration to be statistically significant.

Minnesota↗
USGS2006

Water Quality, Hydrology, and Response to Changes in Phosphorus Loading of Nagawicka Lake, a Calcareous Lake in Waukesha County, Wisconsin

Nagawicka Lake is a 986-acre, usually mesotrophic, calcareous lake in southeastern Wisconsin. Because of concern over potential water-quality degradation of the lake associated with further development in its watershed, a study was conducted by the U.S. Geological Survey from 2002 to 2006 to describe the water quality and hydrology of the lake; quantify sources of phosphorus, including those associated with urban development; and determine the effects of past and future changes in phosphorus loading on the water quality of the lake. All major water and phosphorus sources were measured directly, and minor sources were estimated to construct detailed water and phosphorus budgets for the lake. The Bark River, near-lake surface inflow, precipitation, and ground water contributed 74, 8, 12, and 6 percent of the inflow, respectively. Water leaves the lake primarily through the Bark River outlet (88 percent) or by evaporation (11 percent). The water quality of Nagawicka Lake has improved dramatically since 1980 as a result of decreasing the historical loading of phosphorus to the lake. Total input of phosphorus to the lake was about 3,000 pounds in monitoring year (MY) 2003 and 6,700 pounds in MY 2004. The largest source of phosphorus entering the lake was the Bark River, which delivered about 56 percent of the total phosphorus input, compared with about 74 percent of the total water input. The next largest contributions were from the urbanized near-lake drainage area, which disproportionately accounted for 37 percent of the total phosphorus input but only about 5 percent of the total water input. Simulations with water-quality models within the Wisconsin Lakes Modeling Suite (WiLMS) indicated the response of Nagawicka Lake to 10 phosphorus-loading scenarios. These scenarios included historical (1970s) and current (base) years (MY 2003-04) for which lake water quality and loading were known, six scenarios with percentage increases or decreases in phosphorus loading from controllable sources relative to the base years 2003-04, and two scenarios corresponding to specific management actions. Because of the lake's calcareous character, the average simulated summer concentration of total phosphorus for Nagawicka Lake was about 2 times that measured in the lake. The models likely over-predict because they do not account for coprecipitation of phosphorus and dissolved organic matter with calcite, negligible release of phosphorus from the deep sediments, and external phosphorus loading with abnormally high amounts of nonavailable phosphorus. After adjusting the simulated results for the overestimation of the models, a 50-percent reduction in phosphorus loading resulted in an average predicted phosphorus concentration of 0.008 milligrams per liter (mg/L) (a decrease of 46 percent). With a 50-percent increase in phosphorus loading, the average predicted concentration was 0.020 mg/L (an increase of 45 percent). With the changes in land use under the assumed future full development conditions, the average summer total phosphorus concentration should remain similar to that measured in MY 2003-04 (approximately 0.014 mg/L). However, if stormwater and nonpoint controls are added to achieve a 50-percent reduction in loading from the urbanized near-lake drainage area, the average summer total phosphorus concentration should decrease from the present conditions (MY 2003-04) to 0.011 mg/L. Slightly more than a 25-percent reduction in phosphorus loading from that measured in MY 2003-04 would be required for the lake to be classified as oligotrophic.

Scientific Investigations Report↗
USGS1998

Preliminary estimates of residence times and apparent ages of ground water in the Chesapeake Bay watershed, and water-quality data from a survey of springs

Knowledge of the residence times of the ground-water systems in Chesapeake Bay watershed helps resource managers anticipate potential delays between implementation of land-management practices and any improvements in river and estuary water quality. This report presents preliminary estimates of ground-water residence times and apparent ages of water in the shallow aquifers of the Chesapeake Bay watershed. A simple reservoir model, published data, and analyses of spring water were used to estimate residence times and apparent ages of ground-water discharge. Ranges of aquifer hydraulic characteristics throughout the Bay watershed were derived from published literature and were used to estimate ground-water residence times on the basis of a simple reservoir model. Simple combinations of rock type and physiographic province were used to delineate hydrogeomorphic regions (HGMR's) for the study area. The HGMR's are used to facilitate organization and display of the data and analyses. Illustrations depicting the relation of aquifer characteristics and associated residence times as a continuum for each HGMR were developed. In this way, the natural variation of aquifer characteristics can be seen graphically by use of data from selected representative studies. Water samples collected in September and November 1996, from 46 springs throughout the watershed were analyzed for chlorofluorocarbons (CFC's) to estimate the apparent age of ground water. For comparison purposes, apparent ages of water from springs were calculated assuming piston flow. Additional data are given to estimate apparent ages assuming an exponential distribution of ages in spring discharge. Additionally, results from previous studies of CFC-dating of ground water from other springs and wells in the watershed were compiled. The CFC data, and the data on major ions, nutrients, and nitrogen isotopes in the water collected from the 46 springs are included in this report. The apparent ages of water discharging from 30 of the 46 springs sampled were less than 20 years, including 5 that were "modern" (0-4 years). Four samples had apparent ages of 22 to 34 years, and two others from thermal springs were 40 years or greater. The remaining ten samples were contaminated with local sources of CFC and could not be dated. Nitrate concentrations and nitrate delta 15 nitrogen (d 15 N) values in water from many springs are similar to those in shallow ground water beneath fertilized fields, and some values are high enough to indicate a probable source from animal-waste components. The nitrogen data reported here highlight the significance of the springs sampled during this study as pathways for nutrient transport in the Chesapeake Bay watershed. Ground-water samples were collected from springs during an unusually high flow period and thus may not be representative of low base-flow conditions. Residence times estimated from plausible ranges of aquifer properties and results of previous age-dating analyses generally corroborate the apparent-age analysis made in the current study and suggests that some residence times could be much longer. The shortest residence times tend to be in the Blue Ridge and northern carbonate areas; however, the data are preliminary and not appropriate for statistical tests of significance or variance. Because the age distributions in the aquifer discharging to the springs are not known, and because the apparent ages of water from the springs are based on various combinations of CFC criteria, the apparent ages and calculated residence times are compared for illustrative purposes but are considered preliminary until further work is accomplished.

Maryland, Virginia↗
USGS2014

Assessment of the spatial extent and height of flooding in Lake Champlain during May 2011, using satellite remote sensing and ground-based information

Landsat 5 and moderate resolution imaging spectro-radiometer satellite imagery were used to map the area of inundation of Lake Champlain, which forms part of the border between New York and Vermont, during May 2011. During this month, the lake’s water levels were record high values not observed in the previous 150 years. Lake inundation area determined from the satellite imagery is correlated with lake stage measured at three U.S. Geological Survey lake level gages to provide estimates of lake area at different lake levels (stage/area rating) and also compared with the levels of the high-water marks (HWMs) located on the Vermont side of the lake. The rating developed from the imagery shows a somewhat different relation than a similar stage/area rating developed from a medium-resolution digital elevation model (DEM) of the region. According to the rating derived from the imagery, the lake surface area during the peak lake level increased by about 17 percent above the average or “normal” lake level. By using a comparable rating developed from the DEM, the increase above average is estimated to be about 12 percent. The northern part of the lake (north of Burlington) showed the largest amount of flooding. Based on intersecting the inundation maps with the medium-resolution DEM, lake levels were not uniform around the lake. This is also evident from the lake level gage measurements and HWMs. The gage data indicate differences up to 0.5 feet between the northern and southern end of the lake. Additionally, the gage data show day-to-day and intradaily variation of the same range (0.5 foot). The high-water mark observations show differences up to 2 feet around the lake, with the highest level generally along the south- and west-facing shorelines. The data suggest that during most of May 2011, water levels were slightly higher and less variable in the northern part of the lake. These phenomena may be caused by wind effects as well as proximity to major river inputs to the lake. The inundation areas generated from the imagery generally coincide with flood mapping as estimated by the Federal Emergency Management Agency (FEMA) and shown on its digital flood insurance rate maps. Where areas in the flood inundation map derived from the imagery and the FEMA estimated flooded areas differ substantially, this difference may be due to differences between the flood magnitude at the time of the image and the assumed flood condition used for the FEMA modeling and mapping, wind/storage effects not accounted for by the FEMA modeling, and the resolution of the image compared to the DEM used in the FEMA mapping.

Lake Champlain↗
USGS2001

Preliminary evaluation of the importance of existing hydraulic-head observation locations to advective-transport predictions, Death Valley regional flow system, California and Nevada

When a model is calibrated by nonlinear regression, calculated diagnostic statistics and measures of uncertainty provide a wealth of information about many aspects of the system. This report presents a method of ranking the likely importance of existing observation locations using measures of prediction uncertainty. It is suggested that continued monitoring is warranted at more important locations, and unwarranted or less warranted at less important locations. The report develops the methodology and then demonstrates it using the hydraulic-head observation locations of a three-layer model of the Death Valley regional flow system. The predictions of interest are subsurface transport from beneath Yucca Mountain and 14 Underground Test Areas. The advective component of transport is considered because it is the component most affected by the system dynamics represented by the scale model being used. The problem is addressed using the capabilities of the U.S. Geological Survey computer program MODFLOW-2000, with its ADVective-Travel Observation (ADV) Package, and an additional computer program developed for this work. The methods presented in this report are used in three ways. (1) The ratings for individual observations are obtained by manipulating the measures of prediction uncertainty, and do not involve recalibrating the model. In this analysis, observation locations are each omitted individually and the resulting increase in uncertainty in the predictions is calculated. The uncertainty is quantified as standard deviations on the simulated advective transport. The increase in uncertainty is quantified as the percent increase in the standard deviations caused by omitting the one observation location from the calculation of standard deviations. In general, observation locations associated with larger increases are rated as more important. (2) Ratings for largely geographically based groups are obtained using a straightforward extension of the method used for individual observation locations. This analysis is needed where observations are clustered to determine whether the area is important to the predictions of interest. (3) Finally, the method is used to evaluate omitting a set of 100 observation locations. The locations were selected because they had low individual ratings and were not one of the few locations at which hydraulic heads from deep in the system were measured. The major results of the three analyses, when applied to the three-layer DVRFS ground-water flow system, are described in the following paragraphs. The discussion is labeled using the numbers 1 to 3 to clearly relate it to the three ways the method is used, as listed above. (1) The individual observation location analysis indicates that three observation locations are most important. They are located in Emigrant Valley, Oasis Valley, and Beatty. Of importance is that these and other observations shown to be important by this analysis are far from the travel paths considered. This displays the importance of the regional setting within which the transport occurs, the importance of including some sites throughout the area in the monitoring network, and the importance of including sites in these areas in particular. The method considered in this report indicates that the 19 observation locations that reflect hydraulic heads deeper in the system (in model layers 1, 2, and 3) are not very important. This appears to be because the locations of these observations are in the vicinity of shallow observation locations that also generally are rated as low importance, and because the model layers are hydraulically well connected vertically. The value of deep observations to testing conceptual models, however, is stressed. As a result, the deep observations are rated higher than is consistent with the results of the analysis presented, and none of these observations are omitted in the scenario discussed under (3) below. (2) The geographic grouping of the observations found one major area of importance not identified by the individual observation analysis. Five of the 49 groups are categorized as most highly important. The most important groups were those that, when omitted, produced mean increases greater than 10 percent at any UGTA site or Yucca Mountain. Four of the five groups were dominated by one individual observation. However, one group, located in Ash Meadows, had no individual observations ranked of high importance but collectively, when omitted, increased uncertainty substantially. Other groups also located in Ash Meadows, including intermediate depth observations, consistently ranked as more important than all other groups. (3) To demonstrate the importance of omitting a set of low-rated observations, one scenario is considered in which the 100 individually lowest-rated shallow and intermediate-depth observation locations are omitted. The measure of overall prediction uncertainty increased by just 0.59 percent, indicating that the wells associated with these observations probably could prudently be measured less frequently.

California, Nevada↗
USGS1998

A simulation-optimization model for water-resources management, Santa Barbara, California

In times of drought, the local water supplies of the city of Santa Barbara, California, are insufficient to satisfy water demand. In response, the city has built a seawater desalination plant and gained access to imported water in 1997. Of primary concern to the city is delivering water from the various sources at a minimum cost while satisfying water demand and controlling seawater intrusion that might result from the overpumping of ground water. A simulation-optimization model has been developed for the optimal management of Santa Barbara?s water resources. The objective is to minimize the cost of water supply while satisfying various physical and institutional constraints such as meeting water demand, maintaining minimum hydraulic heads at selected sites, and not exceeding water-delivery or pumping capacities. The model is formulated as a linear programming problem with monthly management periods and a total planning horizon of 5 years. The decision variables are water deliveries from surface water (Gibraltar Reservoir, Cachuma Reservoir, Cachuma Reservoir cumulative annual carryover, Mission Tunnel, State Water Project, and desalinated seawater) and ground water (13 production wells). The state variables are hydraulic heads. Basic assumptions for all simulations are that (1) the cost of water varies with source but is fixed over time, and (2) only existing or planned city wells are considered; that is, the construction of new wells is not allowed. The drought of 1947?51 is Santa Barbara?s worst drought on record, and simulated surface-water supplies for this period were used as a basis for testing optimal management of current water resources under drought conditions. Assumptions that were made for this base case include a head constraint equal to sea level at the coastal nodes; Cachuma Reservoir carryover of 3,000 acre-feet per year, with a maximum carryover of 8,277 acre-feet; a maximum annual demand of 15,000 acre-feet; and average monthly capacities for the Cachuma and the Gibraltar Reservoirs. The base-case results indicate that water demands can be met, with little water required from the most expensive water source (desalinated seawater), at a total cost of $5.56 million over the 5-year planning horizon. The simulation model has drains, which operate as nonlinear functions of heads and could affect the model solutions. However, numerical tests show that the drains have little effect on the optimal solution. Sensitivity analyses on the base case yield the following results: If allowable Cachuma Reservoir carryover is decreased by about 50 percent, then costs increase by about 14 percent; if the peak demand is decreased by 7 percent, then costs will decrease by about 14 percent; if the head constraints are loosened to -30 feet, then the costs decrease by about 18 percent; if the heads are constrained such that a zero hydraulic gradient condition occurs at the ocean boundary, then the optimization problem does not have a solution; if the capacity of the desalination plant is constrained to zero acre-feet, then the cost increases by about 2 percent; and if the carryover of State Water Project water is implemented, then the cost decreases by about 0.5 percent. Four additional monthly diversion distribution scenarios for the reservoirs were tested: average monthly Cachuma Reservoir deliveries with the actual (scenario 1) and proposed (scenario 2) monthly distributions of Gibraltar Reservoir water, and variable monthly Cachuma Reservoir deliveries with the actual (scenario 3) and proposed (scenario 4) monthly distributions of Gibraltar Reservoir water. Scenario 1 resulted in a total cost of about $7.55 million, scenario 2 resulted in a total cost of about $5.07 million, and scenarios 3 and 4 resulted in a total cost of about $4.53 million. Sensitivities of the scenarios 1 and 2 to desalination-plant capacity and State Water Project water carryover were tested. The scenario 1 sensitivity analysis indicated that incorpo

Water-Resources Investigations Report↗
USGS2020

Hydrogeologic characterization of the Hualapai Plateau on the western Hualapai Indian Reservation, northwestern Arizona

This study was developed to assess if groundwater from the western Hualapai Plateau could be used to supply developments in the Grand Canyon West area of the Hualapai Indian Reservation and to collect hydrogeologic data for future use in a numerical groundwater model for the reservation. Ground-based geophysical surveys; existing well, spring, and other hydrogeologic information from previous studies; and new well and spring data collected for this study were used to provide a better understanding of the hydrogeology of the western Hualapai Plateau. Surface geophysical data provided information on the depth and geologic structure of lower Paleozoic rock units and Proterozoic crystalline and metamorphic rocks that underlie the western Hualapai Plateau. The surface geophysical data and discharge information from springs were used to select a site to drill and develop the U.S. Geological Survey Hualapai Test Well. The Hualapai Test Well was drilled to understand the geophysical properties of geologic formations at depth. These data were used to verify the results of surface geophysical data and to evaluate if sufficient water was present in the Hualapai Test Well for potential groundwater development. The Hualapai Test Well was drilled to a depth of 2,468 feet and bottomed in Proterozoic granite. Water was expected in the lower part of the Muav Limestone, but water was not observed until the Tapeats Sandstone at a depth of 2,400 feet. The Tapeats Sandstone was determined to be confined with a hydrostatic head of over 900 feet. A 48-hour pumping test was conducted to determine aquifer properties. Low specific capacity indicated that although groundwater is present in the Tapeats Sandstone, well yields are likely to be small. A water-quality sample indicated the sample had a calcium, magnesium-bicarbonate water type with a total dissolved-solids concentration of 371 milligrams per liter. Alpha radioactivity of the sample, 18.3 picocuries per liter, exceeded the U.S. Environmental Protection Agency maximum contaminant level of 15 picocuries per liter for drinking water. Concentrations of iron and manganese in the water sample also exceeded the U.S. Environmental Protection Agency secondary maximum contaminant levels for drinking water. An inventory of wells and springs provided insight into the occurrence of groundwater on the western Hualapai Plateau. Data from 56 springs on and adjacent to the western Hualapai Plateau were compiled for this study, and new data were collected at 31 springs. Discharge from springs visited for this study ranged from dry to about 345 gallons per minute. The temporal data from springs, where repeat measurements were available, indicated that spring flow is highly variable and likely related to seasonal and annual precipitation. Water levels from 36 wells on and adjacent to the western Hualapai Plateau were compiled for this study, and new water levels were collected at 5 wells. The spring and well data in conjunction with the Hualapai Test Well results indicated that on the western Hualapai Plateau, bedrock aquifers have limited discrete flow paths that make extensive groundwater development unlikely.

Arizona↗