USGS Science⌕ Search

SEARCH · USGS Science

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 757 records · Page 42Linked to original sources

Simulation of ground-water flow and pumpage in Kings and Queens Counties, Long Island, New York

The potential effects of using ground water as a supplemental source of supply in Kings and Queens Counties were evaluated through a 4-layer finite-difference ground-water-flow model with a uniform grid spacing of 1,333 feet. Hydraulic properties and boundary conditions of an existing regional ground-water-flow model of Long Island with a uniform grid spacing of 4,000 feet were refined for use in the finer grid model of Kings and Queens Counties. The model is calibrated to average pumping stresses that correspond to presumed steady-state conditions of 1983 and 1991. A transient-state simulation of the year-by- year transition between these two conditions was also conducted. Pumping scenarios representing public-supply withdrawals of 100, 150, and 400 million gallons per day (Mgal/d) were simulated to determine the duration of sustainable pumpage, defined as the length of time before a particular pumping rate induces landward hydraulic gradients from areas of salty ground water. The simulations indicate the following hydrologically feasible scenarios: (1) Pumpage of 100 Mgal/d could be sustained for about 10 months, followed by a 46-month period of pumping at reduced (1991) rates, to allow water levels to recover to 90 percent of 1991 levels. (2) Pumpage of 150 Mgal/d could be sustained for about 6 months, followed by a 79-month period of pumping at a reduced (1991) rate. (3) Pumpage of 400 Mgal/d could be sustained for about 3 months from an initial condition of maximum aquifer storage. Each of these scenarios could be modified by injecting surplus water from upstate reservoirs, available from January to May, into the proposed wells. Injection at half the pumpage rate during the recovery period reduces the recovery period to 14 months in scenario 1, 6 months in scenario 2, and 9 months in scenario 3.

New York↗

Geophysical framework investigations influencing ground-water resources in east-central Nevada and west-central Utah, with a section on geologic and geophysical basin by basin descriptions

A geophysical investigation was undertaken as part of an effort to characterize the geologic framework influencing ground-water resources in east-central Nevada and west-central Utah. New gravity data were combined with existing aeromagnetic, drill-hole, and geologic data to help determine basin geometry, infer structural features, estimate depth to pre-Cenozoic basement rocks, and further constrain the horizontal extents of exposed and buried plutons. In addition, a three-dimensional (3D) geologic model was constructed to help illustrate the often complex geometries of individual basins and aid in assessing the connectivity of adjacent basins. In general, the thirteen major valleys within the study area have axes oriented north-south and frequently contain one or more sub-basins. These basins are often asymmetric and typically reach depths of 2 km. Analysis of gravity data helped delineate geophysical lineaments and accommodation zones. Structural complexities may further compartmentalize ground-water flow within basins and the influence of tectonics on basin sedimentation further complicates their hydrologic properties. The horizontal extent of exposed and, in particular, buried plutons was estimated over the entire study area. The location and subsurface extents of these plutons will be very important for regional water resource assessments, as these features may act as either barriers or pathways for groundwater flow. A previously identified basement gravity low strikes NW within the study area and occurs within a highly extended terrane between the Butte and Confusion synclinoria. Evidence from geophysical, geologic, and seismic reflection data suggests relatively lower density plutonic rocks may extend to moderate crustal depths and rocks of similar composition may be the source of the entire basement gravity anomaly.

Nevada, Utah↗

Hydrogeology, water quality, and ground-water-development alternatives in the upper Wood River ground-water reservoir, Rhode Island

The 72.4-square-mile Upper Wood River study area is in the Pawcatuck River basin in southern Rhode Island. Stratified drift is the only principal geologic unit capable of producing yields greater than 0.5 Mgal/d. Transmissivity of the aquifer ranges from 7,600 to 49,200 sq ft/d. Water-table conditions prevail and the aquifer is in good hydraulic connection with perennial streams and ponds. Groundwater and surface water in the study area are generally suitable for most uses. Water is soft, slightly acidic, and contains less than 150 mg/L dissolved solids. Locally, however, groundwater has been contaminated with nitrate, chloride, and volatile organic compounds. A model of the groundwater-flow system was used to evaluate the effect of alternative schemes of groundwater development on water levels, pond levels, and streamflow. Till contacts were simulated as specified-flux boundaries, drainage divides as no-flow boundaries, and streams as leaky boundaries. The areas most favorable for development of 1 Mgal/d are along the Flat and Wood Rivers. From 50 to 65 percent of the water withdrawn from wells would be derived from induced recharge. Results of simulation of development alternatives indicate that the groundwater reservoir could sustain withdrawals of 6 to 12 Mgal/d from 11 wells under long-term average annual (1942-89) and simulated drought (1963-66) conditions without causing water-level declines of greater than 25 percent of the unstressed saturated thickness of the aquifer. Pumping 12 Mgal/d, however, would reduce flow of the Wood River at the basin outlet by an amount almost equal to the 7-day, 10-yr low flow of 20.4 cu ft/s.

Water-Resources Investigations Report↗

Development, calibration, and testing of ground-water flow models for the Mississippi River Valley alluvial aquifer in eastern Arkansas using one-square-mile cells

Significant water-level declines in the Mississippi River Valley alluvial aquifer prompted the need to better understand the flow system in the aquifer which, in turn, led to the development of digital groundwater flow models of the alluvial aquifer. Two models were developed in the eastern Arkansas study area with the Arkansas River dividing the study area and functioning as a hydrologic boundary to the models. Both models simulate groundwater flow in one layer with recharge entering the aquifer from head-dependent surface infiltration through the overlying confining unit and from seepage through river beds. Digital models were used to simulate flow in the aquifer during seven stress periods between 1918 and 1987. Pumpage used in the simulations ranged from 83,400,000 to 412,000,000 cu ft/d in the north model and from 12,800,000 to 58,500,000 cu ft/d in the south model. Three different spatial and temporal pumpage scenarios were tested to simulate pumpage stress in the models. The pumpage distribution used in the calibrated model was based on a combination of all three scenarios. Several criteria were used during model development to determine how well the model simulated conditions in the aquifer. Potentiometric maps of model-computed water levels were compared to measured data to check the computed water levels and direction of flow. Hydrographs of observation wells were compared to computed water levels at corresponding model cells to assess the temporal distribution of pumpage. A root-mean-square error analysis was performed during calibration by comparing observation-well and model-computed water levels for 1972. Sensitivity analyses were performed to determine the effects of changes in input parameters on computed heads (water levels). Both models were sensitive to changes in recharge and pumpage but the south model generally was less sensitive than the north model.

Water-Resources Investigations Report↗

Hydrogeology and simulated effects of ground-water development on an unconfined aquifer in the Closed Basin Division, San Luis Valley, Colorado

Wells completed in an unconfined aquifer in the Closed Basin Division of the San Luis Valley Project, Colorado, are expected to provide about 101,800 acre-ft of groundwater/year to the Rio Grande when this project is completed. Lowering of groundwater levels in the unconfined aquifer is expected to decrease the quantity of groundwater that is lost by evapotranspiration. The aquifer system, which consists of an unconfined aquifer that is 50 to 130 ft thick, overlies a thick, leaky confined aquifer. Groundwater moves from the edge of the valley toward a topographic low near the center of the Closed Basin Division, where it is lost by evapotranspiration. A two-dimensional groundwater flow model was used to evaluate the effects of projected withdrawal of about 141 cu ft/sec by 168 wells throughout a 20-year period. The simulated pumpage resulted in a projected drawdown greater than 0.1 ft in the water-levels of the unconfined aquifer over an area of about 370 sq mi. Maximum simulated drawdown was 25 ft. Simulations indicate that about 66 % of the water to be withdrawn from the unconfined aquifer would be derived from decreases of evapotranspiration, 26% from induced leakage from an underlying confined aquifer, and 8% from storage of the unconfined aquifer. Model simulations were based only on withdrawals from wells completed in the unconfined aquifer. Pumpage from the confined aquifer was not simulated. Upward leakage from the confined aquifer predicted by the model, results from the simulated declines of the potentiometric surface in the unconfined aquifer. (USGS)

Water-Resources Investigations Report↗

Regional analysis of ground-water recharge

A modeling analysis of runoff and ground-water recharge for the arid and semiarid southwestern United States was performed to investigate the interactions of climate and other controlling factors and to place the eight study-site investigations into a regional context. A distributed-parameter water-balance model (the Basin Characterization Model, or BCM) was used in the analysis. Data requirements of the BCM included digital representations of topography, soils, geology, and vegetation, together with monthly time-series of precipitation and air-temperature data. Time-series of potential evapotranspiration were generated by using a submodel for solar radiation, taking into account topographic shading, cloudiness, and vegetation density. Snowpack accumulation and melting were modeled using precipitation and air-temperature data. Amounts of water available for runoff and ground-water recharge were calculated on the basis of water-budget considerations by using measured- and generated-meteorologic time series together with estimates of soil-water storage and saturated hydraulic conductivity of subsoil geologic units. Calculations were made on a computational grid with a horizontal resolution of about 270 meters for the entire 1,033,840 square-kilometer study area. The modeling analysis was composed of 194 basins, including the eight basins containing ground-water recharge-site investigations. For each grid cell, the BCM computed monthly values of potential evapotranspiration, soil-water storage, in-place ground-water recharge, and runoff (potential stream flow). A fixed percentage of runoff was assumed to become recharge beneath channels operating at a finer resolution than the computational grid of the BCM. Monthly precipitation and temperature data from 1941 to 2004 were used to explore climatic variability in runoff and ground-water recharge. The selected approach provided a framework for classifying study-site basins with respect to climate and dominant recharge processes. The average climate for all 194 basins ranged from hyperarid to humid, with arid and semiarid basins predominating (fig. 6, chapter A, this volume). Four of the 194 basins had an aridity index of dry subhumid; two of the basins were humid. Of the eight recharge-study sites, six were in semiarid basins, and two were in arid basins. Average-annual potential evapotranspiration showed a regional gradient from less than 1 m/yr in the northeastern part of the study area to more than 2 m/yr in the southwestern part of the study area. Average-annual precipitation was lowest in the two arid-site basins and highest in the two study-site basins in southern Arizona. The relative amount of runoff to in-place recharge varied throughout the study area, reflecting differences primarily in soil water-holding capacity, saturated hydraulic conductivity of subsoil materials, and snowpack dynamics. Climatic forcing expressed in El Niño and Pacific Decadal Oscillation indices strongly influenced the generation of precipitation throughout the study area. Positive values of both indices correlated with the highest amounts of runoff and ground-water recharge.

Professional Paper↗

Potentiometric Surface of the Ozark Aquifer near Springfield, Missouri, 2006-07

INTRODUCTION A study of the water resources of the Springfield, Missouri, area in the 1970s determined that a cone of depression, formed by ground-water pumping, had developed in the Ozark aquifer beneath the city (Emmett and others, 1978). Continued ground-water usage in the 1970s and 1980s caused concern that ground-water resources would not be sufficient to meet the future needs of Springfield, Missouri, during periods of drought. As a result, a ground-water flow model of the Springfield area was developed by the U. S. Geological Survey (USGS) to assess the future role of ground water as a water source for the area (Imes, 1989). Results of the USGS model led to a decision by the City Utilities of Springfield to primarily rely on surface water from Stockton Lake as a source of city drinking water. Municipal and industrial ground-water usage continues in Springfield, but at lower rates than previously experienced (Jim Vandike, Missouri Department of Natural Resources, written commun., 2007). Rapid growth in the area has caused commercial, industrial, and domestic water use to increase. Population growth has been especially rapid in Nixa, Ozark, and Republic, and water use in the vicinity of these cities has grown an estimated 39 percent since 1990 (Dintelmann and others, 2006). Unlike Springfield, ground water is the primary source of water for these cities. The increased stress on the Ozark aquifer, the primary aquifer in the study area, has raised new concerns about possible further water-level declines in the areas of increased ground-water use. Although there continues to be new development in the Ozark aquifer, since 1987 no new water-supply wells that produce water from the Springfield Plateau aquifer have been allowed to be constructed in most of Greene and northern Christian counties (Jim Vandike, Missouri Department of Natural Resources, written commun., 2007). There is concern that if the potentiometric surface of the Ozark aquifer continues to decline, increased leakage of contaminants into the Ozark aquifer from the overlying Springfield Plateau aquifer could occur (Jim Vandike, Missouri Department of Natural Resources, written commun., 2007). To address this concern, the USGS, in cooperation with Greene County, Missouri, the U.S. Army Corps of Engineers, and the Missouri Department of Natural Resources, constructed a map of the potentiometric surface of the Ozark aquifer for 2006?2007. The map can be compared to previously constructed potentiometric-surface maps by Emmett and others (1978) and Imes (1989) to evaluate changes in ground-water flow directions, but the comparison is beyond the scope of this report.

Scientific Investigations Map↗

Ground-Water Availability in the Wailuku Area, Maui, Hawai'i

Most of the public water supply in Maui, Hawai'i, is from a freshwater lens in the Wailuku area of the island. Because of population growth, ground-water withdrawals from wells in this area increased from less than 10 Mgal/d during 1970 to about 23 Mgal/d during 2006. In response to increased withdrawals from the freshwater lens in the Wailuku area, water levels declined, the transition zone between freshwater and saltwater became shallower, and the chloride concentrations of water pumped from wells increased. These responses led to concern over the long-term sustainability of withdrawals from existing and proposed wells. A three-dimensional numerical ground-water flow and transport model was developed to simulate the effects of selected withdrawal and recharge scenarios on water levels, on the transition zone between freshwater and saltwater, and on surface-water/ground-water interactions. The model was constructed using time-varying recharge, withdrawals, and ocean levels. Hydraulic characteristics used to construct the model were initially based on published estimates but ultimately were varied to obtain better agreement between simulated and measured water levels and salinity profiles in the modeled area during the period 1926-2006. Scenarios included ground-water withdrawal at 2006 and 1996 rates and locations with average recharge (based on 2000-04 land use and 1926-2004 rainfall) and withdrawal at redistributed rates and locations with several different recharge scenarios. Simulation results indicate that continuing 1996 and 2006 withdrawal distributions into the future results in decreased water levels, a thinner freshwater lens, increased salinity from pumped wells, and higher salinity at several current withdrawal sites. A redistributed withdrawal condition in which ground-water withdrawal was redistributed to maximize withdrawal and minimize salinities in the withdrawn water was determined. The redistributed withdrawal simulates 27.1 Mgal/d of withdrawal from 14 wells or well fields in the Wailuku area. Simulation results from the five scenarios that include redistributed withdrawal conditions indicate the following for the Wailuku Aquifer Sector: (1) withdrawal during times of average recharge rates cause average water levels to decrease 2-3 ft and the transition zone to become more than 200 ft shallower after 150 years; (2) a 5-yr drought condition similar to the 1998-2002 drought results in additional salinity increases after 30 years (12.5 years of normal recharge after drought conditions) but only one well has salinity increases of concern; (3) additional recharge from restored streamflow significantly increases water levels, thickens the freshwater body, and decreases salinity at withdrawal sites in the Waihe'e and 'Iao Aquifer Systems; and (4) a complete removal of irrigation recharge decreases water levels and increases salinity in the central isthmus where irrigation is reduced, but recharge through restored streams still significantly increases water levels, thickens the freshwater body, and decreases salinity at withdrawal sites in the Waihe'e and 'Iao Aquifer Systems.

Scientific Investigations Report↗

Nutrient Loads and Ground-Water Residence Times in an Agricultural Basin in North-Central Connecticut

Nutrient loads from ground-water discharge were studied in Broad Brook Basin, a 15.8-square mile basin in north-central Connecticut, dominated by agricultural activity. Loads were calculated, along with the travel times of ground water from recharge to discharge areas, to estimate the time required for the effects of Best Management Practices (BMPs) to be observed. Most concentrations of nitrogen and phosphorus in Broad Brook exceeded U.S. Environmental Protection Agency Ecoregion XIV nutrient criteria for streams. During the study period (1993-2004), annual loads of nitrogen from Broad Brook Basin ranged from 117,000 to 270,000 pounds (lb), and yields were about 10 times larger than those from forested basins in Connecticut. Ground-water discharge from the aquifer to the streams (base flow) during the study period was estimated with hydrograph separation and accounted for 82 percent of the total runoff from the basin. Nitrate nitrogen in base flow averaged 71 percent of the annual load of total nitrogen discharged from the basin, indicating that the largest source of nitrogen was likely from ground-water discharge. Annual loads of total phosphorus from the basin ranged from 2,330 to 14,400 lb, and yields were about five times higher than those from forested basins in Connecticut. Dissolved phosphorus averaged about 71 percent of the total phosphorus load, and ground-water discharge accounted for only as much as 40 percent of the annual load of dissolved phosphorus; therefore, phosphorus loads are dominated by stormwater-runoff events. Ground-water samples collected from 11 wells in the basin contained elevated concentrations of nitrite plus nitrate nitrogen. Dissolved gas analyses indicated that little denitrification was occurring in the aquifer. Apparent ages of the ground-water samples ranged from greater than 2 to more than 50 years based on sulfur hexafluoride, tritium, and tritium/helium-3 analyses. A three-dimensional ground-water-flow model was used in conjunction with a particle-tracking program to determine advective travel times to streams from all subareas in the basin. The model simulations indicated that about half the discharge to Broad Brook and its tributaries was recharged more than 10 years ago, and that about 8 percent of the discharge was recharged prior to 1960. The effects of changes in nitrate nitrogen loading at the water table were evaluated by applying new loading rates from urban and agricultural land and the simulated advective ground-water travel times. Five scenarios were tested: reducing estimated nitrate nitrogen concentrations in recharging ground water under urban and agricultural land areas to concentrations in forested areas; reducing estimated nitrate nitrogen concentrations under urban and agricultural land areas to the U.S. Environmental Protection Agency recommended nutrient criteria for streams; and reducing estimated nitrate nitrogen concentrations under urban and agricultural land areas by 50 percent, 10 percent, and 5 percent. Under the first two scenarios, the base-flow load of nitrate nitrogen could be reduced by 25 percent in slightly more than 5 years, although the reduction required by these scenarios is likely unrealistic. A 25-percent reduction in base-flow load of nitrate nitrogen could be achieved in about 10 years under the third scenario (where concentrations from urban and agricultural areas were reduced by 50 percent). Under this scenario, a 46-percent reduction could be achieved in about 60 years. The scenarios indicate that in this basin, and in other similar basins in Connecticut underlain by extensive glacial stratified deposits, there can be a substantial time lag between activities at the land surface and effects on the quality of water discharged to streams from ground water. This finding may be important in the expectations for water-quality improvements from land-use changes or BMPs.

Scientific Investigations Report↗

Ground-water hydrology of the Toppenish Creek basin, Yakima Indian Reservation, Washington

A groundwater flow model was constructed for the Toppenish Creek basin aquifer system in eastern Washington. Flow was simulated in three aquifer units: (1) the confined old valley fill and shallow basalt (unit 2); (2) the underlying primary basalt (unit 3); and (3) the deep basalt (unit 4). Water levels in the overlying unconfined alluvial aquifer (unit 1) were held fixed. Calibrated transmissivities ranged from 0.01 to 0.48 foot squared per second. Calibrated storage coefficients were 0.0004 to 0.006. The confining-bed leakance ranged from 2.0 x 10 to the minus 11th power to 2.5 x 10 to the minus 10th power feet per second per foot. Under steady-state conditions (1954) annual natural recharge was about 29,000 acre-ft underflow from adjacent basins. Annual pumpage increased from less than 500 acre-ft in 1954 to an average of 19,600 acre-ft for 1971 and 1972. Pumpage caused simulated declines in unit 3 of up to 95 ft for 1955-72. Projected annual declines from 1973-77 using 1971-72 pumpage were about 0.5-1.5 ft in unit 2 and 0.2-1.5 ft in unit 3. The corresponding declines from 1978-82 were 0.2-1.2 ft and 0.2-0.8 ft, respectively. Using 1971-72 pumpage plus 12 ,400 acre-ft per year from unit 3, the calculated annual declines from 1978-82 were 1-36 ft in unit 2 and 4-20 ft in unit 3. (USGS)

Water-Resources Investigations Report↗

Nutrients in streams and rivers across the nation — 1992–2001

Nutrient compounds of nitrogen and phosphorus were investigated in streams and rivers sampled as part of the U.S. Geological Survey National Water-Quality Assessment (NAWQA) Program. Nutrient data were collected in 20 NAWQA study units during 1992-95, 16 study units during 1996-98, and 15 study units during 1999-2001. To facilitate comparisons among sampling sites with variable sampling frequency, daily loads were determined by using regression models that relate constituent transport to streamflow and time. Model results were used to compute mean annual loads, yields, and concentrations of ammonia, nitrate, total nitrogen, orthophosphate, and total phosphorus, which were compared among stream and river sampling sites. Variations in the occurrence and distribution of nutrients in streams and rivers on a broad national scale reflect differences in the sources of nutrient inputs to the upstream watersheds and in watershed characteristics that affect movement of those nutrients. Sites were classified by watershed size and by land use in the upstream watershed: agriculture, urban, and undeveloped (forest or rangeland). Selection of NAWQA urban sites was intended to avoid effects of major wastewater-treatment plants and other point sources, but in some locations this was not feasible. Nutrient concentrations and yields generally increased with anthropogenic development in the watershed. Median concentrations and yields for all constituents at sites downstream from undeveloped areas were less than at sites downstream from agricultural or urban areas. Concentrations of ammonia, orthophosphate, and total phosphorus at agricultural and urban sites were not significantly different; however, concentrations of nitrate and total nitrogen were higher at agricultural than at urban sites. Total nitrogen concentrations at agricultural sites were higher in areas of high nitrogen input or enhanced transport, such as irrigation or artificial drainage that can rapidly move water from cropland to streams (Midwest, Northern Plains, and western areas of the United States). Concentrations were lower in the Southeast, where more denitrification occurs during transport of nitrogen compounds in shallow ground water. At urban sites, high concentrations of ammonia and orthophosphate were more prevalent downstream from wastewater-treatment plants. At sites with large watersheds and high mean-annual streamflow (“large-watershed” sites), concentrations of most nutrients were significantly less than at sites downstream from agricultural or urban areas. Total nitrogen concentrations at large-watershed sites were higher in Midwest agricultural areas and lower in the Western United States, where agricultural and urban development is less extensive. Total phosphorus concentrations at large-watershed sites were higher in areas of greater potential erosion and low overall runoff such as the arid areas in the West. Although not as distinct as seasonal patterns of streamflow, geographic patterns of seasonally high and low concentrations of total nitrogen and total phosphorus were identified in the data. Seasonal patterns in concentrations of total nitrogen generally mirror seasonal patterns in streamflow in the humid Eastern United States but are inverse to seasonal patterns in streamflow in the semiarid interior West. Total phosphorus concentrations typically have the opposite regional relation with streamflow; high concentrations coincide with high streamflows in the interior West. In the NAWQA Program, sites downstream from relatively undeveloped areas were selected to provide a baseline for comparison to sites with potential effects of urban development and agriculture. Concentrations of nitrate, total nitrogen, and total phosphorus at NAWQA undeveloped sites were found to be greater than values reported by other studies for conditions of essentially no development (background conditions). Concentrations at NAWQA undeveloped sites represent conditions of relatively little development and provide insight in comparison to developed areas but should not, in general, be considered to represent background status. The U.S. Environmental Protection Agency has developed nutrient criteria to assist States in setting regional water-quality standards. Regional criteria were exceeded by total nitrogen concentrations at 72 percent of NAWQA undeveloped sites and by total phosphorus concentrations at 89 percent of these sites. Exceedances were even more extensive at sites with greater anthropogenic development upstream. The nitrogen criteria were exceeded at 96 percent of NAWQA sites classified as agricultural, urban, or mixed land use, and the phosphorus criteria were exceeded at 97 percent of these sites. Nationally, outflow loads of all nutrient constituents were strongly correlated to nonpoint-source inputs in the upstream watershed. The variation in input mass explained at least 69 percent of the variation in loads. Correlations between nonpoint-source input rates and outflow yields were not quite as good; variation in input rates explained only about 22-45 percent of the variations in nutrient yields. Estimation of nutrient outflow, on the basis of these correlations, likely could be improved if nationally consistent data were available for additional watershed characteristics.

Scientific Investigations Report↗

Summary of and factors affecting pesticide concentrations in streams and shallow wells of the lower Susquehanna River basin, Pennsylvania and Maryland, 1993-95

This report presents the detection frequency of 83 analyzed pesticides, describes the concentrations of those pesticides measured in water from streams and shallow wells, and presents conceptual models of the major factors affecting seasonal and areal patterns of pesticide concentrations in water from streams and shallow wells in the Lower Susquehanna River Basin. Seasonal and areal patterns of pesticide concentrations were observed in 577 samples and nearly 40,000 pesticide analyses collected from 155 stream sites and 169 shallow wells from 1993 to 1995. For this study, shallow wells were defined as those generally less than 200 feet deep. The most commonly detected pesticides were agricultural herbicides?atrazine, metolachlor, simazine, prometon, alachlor, and cyanazine. Atrazine and metolachlor are the two most-used agricultural pesticides in the Lower Susquehanna River Basin. Atrazine was detected in 92 percent of all the samples and in 98 percent of the stream samples. Metolachlor was detected in 83 percent of all the samples and in 95 percent of the stream samples. Nearly half of all the analyzed pesticides were not detected in any sample. Of the 45 pesticides that were detected at least once, the median concentrations of 39 of the pesticides were less than the detection limit for the individual compounds, indicating that for at least 50 percent of the samples collected, those pesticides were not detected. Only 10 (less than 0.025 percent) of the measured concentrations exceeded any established drinking-water standards; 25 concentrations exceeded 2 mg/L (micrograms per liter) and 55 concentrations exceeded 1 mg/L. None of the elevated concentrations were measured in samples collected from streams that are used for public drinking-water supplies, and 8 of the 10 were measured in storm-affected samples. The timing and rate of agricultural pesticide applications affect the seasonal and areal concentration patterns of atrazine, simazine, chlorpyrifos, and diazinon observed in water from wells and streams in the Lower Susquehanna River Basin. Average annual pesticide use for agricultural purposes and nonagricultural pesticide use indicators were used to explain seasonal and areal patterns. Elevated concentrations of some pesticides in streams during base-flow and storm-affected conditions were related to the seasonality of agricultural-use applications and local climate conditions. Agricultural-use patterns affected areal concentration patterns for the high-use pesticides, but indicators of nonagricultural use were needed to explain concentration patterns of pesticides with smaller amounts used for agricultural purposes. Bedrock type influences the movement and discharge of ground water, which in turn affects concentration patterns of pesticides. The ratio of atrazine concentrations in stream base flow to concentrations in shallow wells varied among the different general rock types found in the Lower Susquehanna River Basin. Median concentrations of atrazine in well water and stream base flow tended to be similar in individual areas underlain by carbonate bedrock, indicating the connectivity of water in streams and shallow wells in these areas. In areas underlain by noncarbonate bedrock, median concentrations of atrazine tended to be significantly higher in stream base flow than in well water. This suggests a deep ground-water system that delivers water to shallow wells and a near-surficial system that supplies base-flow water to streams. In addition to the presence or absence of carbonate bedrock, pesticide leaching potential and persistence, soil infiltration capacity, and agricultural land use affected areal patterns in detection frequency and concentration differences between samples collected from streams during base-flow conditions and shallow wells.

Water-Resources Investigations Report↗

Regional rainfall-runoff relations for simulation of streamflow for watersheds in Du Page County, Illinois

Rainfall and streamflow data collected from July 1986 through September 1993 were utilized to calibrate and verify a continuous-simulation rainfall-runoff model for three watersheds (11.8--18.0 square miles in area) in Du Page County. Classification of land cover into three categories of pervious (grassland, forest/wetland, and agricultural land) and one category of impervious subareas was sufficient to accurately simulate the rainfall-runoff relations for the three watersheds. Regional parameter sets were obtained by calibrating jointly all parameters except fraction of ground-water inflow that goes to inactive ground water (DEEPFR), interflow recession constant (IRC), and infiltration (INFILT) for runoff from all three watersheds. DEEPFR and IRC varied among the watersheds because of physical differences among the watersheds. Two values of INFILT were obtained: one representing the rainfall-runoff process on the silty and clayey soils on the uplands and lake plains that characterize Sawmill Creek, St. Joseph Creek, and eastern Du Page County; and one representing the rainfall-runoff process on the silty soils on uplands that characterize Kress Creek and parts of western Du Page County. Regional rainfall-runoff relations, defined through joint calibration of the rainfall-runoff model and verified for independent periods, presented in this report, allow estimation of runoff for watersheds in Du Page County with an error in the total water balance less than 4.0 percent; an average absolute error in the annual-flow estimates of 17.1 percent with the error rarely exceeding 25 percent for annual flows; and correlation coefficients and coefficients of model-fit efficiency for monthly flows of at least 87 and 76 percent, respectively. Close reproduction of the runoff-volume duration curves was obtained. A frequency analysis of storm-runoff volume indicates a tendency of the model to undersimulate large storms, which may result from underestimation of the amount of impervious land cover in the watershed and errors in measuring rainfall for convective storms. Overall, the results of regional calibration and verification of the rainfall-runoff model indicate the simulated rainfall-runoff relations are adequate for stormwater-management planning and design for watersheds in Du Page County.

Water-Resources Investigations Report↗

Geology, ground-water hydrology, geochemistry, and ground-water simulation of the Beaumont and Banning Storage Units, San Gorgonio Pass area, Riverside County, California

Ground water has been the only source of potable water supply for residential, industrial, and agricultural users in the Beaumont and Banning storage units of the San Gorgonio Pass area, Riverside County, California. Ground-water levels in the Beaumont area have declined as much as 100 feet between the early 1920s and early 2000s, and numerous natural springs have stopped flowing. In 1961, the San Gorgonio Pass Water Agency (SGPWA) entered into a contract with the California State Department of Water Resources to receive 17,300 acre-feet per year of water to be delivered by the California State Water Project (SWP) to supplement natural recharge. Currently (2005), a pipeline is delivering SWP water into the area, and the SGPWA is artificially recharging the ground-water system using recharge ponds located along Little San Gorgonio Creek in Cherry Valley with the SWP water. In addition to artificial recharge, SGPWA is considering the direct delivery of SWP water for the irrigation of local golf courses and for agricultural supply in lieu of ground-water pumpage. To better understand the potential hydrologic effects of different water-management alternatives on ground-water levels and movement in the Beaumont and Banning storage units, existing geohydrologic and geochemical data were compiled, new data from a basin-wide ground-water level and water-quality monitoring network were collected, monitoring wells were installed near the Little San Gorgonio Creek recharge ponds, geohydrologic and geochemical analyses were completed, and a ground-water flow simulation model was developed. The San Gorgonio Pass area was divided into several storage units on the basis of mapped or inferred faults. This study addresses primarily the Beaumont and Banning storage units. The geologic units in the study area were generalized into crystalline basement rocks and sedimentary deposits. The younger sedimentary deposits and the surficial deposits are the main water-bearing deposits in the San Gorgonio Pass area. The water-bearing deposits were divided into three aquifers: (1) the perched aquifer, (2) the upper aquifer, and (3) the lower aquifer based on lithologic and downhole geophysical logs. Natural recharge in the San Gorgonio Pass area was estimated using INFILv3, a deterministic distributed- parameter precipitation-runoff model. The INFILv3 model simulated that the potential recharge of precipitation and runoff in the Beaumont and Banning storage units was about 3,710 acre-feet per year and that the potential recharge in 28 sub-drainage basins upstream of the storage units was about 6,180 acre-feet per year. The water supply for the Beaumont and Banning storage units is supplied by pumping ground water from wells in the Canyon (Edgar and Banning Canyons), Banning Bench, Beaumont, and Banning storage units. Total annual pumpage from the Beaumont and Banning storage units ranged from about 1,630 acre-feet in 1936 to about 20,000 acre-feet in 2003. Ground-water levels declined by as much as 100 feet in the Beaumont storage unit from 1926-2003 in response to ground-water pumping of about 450,160 acre-feet during this period. Since ground-water development began in the San Gorgonio Pass area, there have been several sources of artificial recharge to the basin including return flow from applied water on crops, golf courses, and landscape; septic-tank seepage; and infiltration of storm runoff diversions and imported water into recharge ponds. Return flow from applied water and septic-tank seepage was estimated to reach a maximum of about 8,100 acre-feet per year in 2003. Owing to the great depth of water in much of study area (in excess of 150 feet), the return flow and septic-tank seepage takes years to decades to reach the water table. Stable-isotope data indicate that the source of ground-water recharge was precipitation from storms passing through the San Gorgonio Pass as opposed to runoff from the higher altitudes of the San Bernar

Scientific Investigations Report↗

Net-infiltration map of the Navajo Sandstone outcrop area in western Washington County, Utah

As populations grow in the arid southwestern United States and desert bedrock aquifers are increasingly targeted for future development, understanding and quantifying the spatial variability of net infiltration and recharge becomes critically important for inventorying groundwater resources and mapping contamination vulnerability. A Geographic Information System (GIS)-based model utilizing readily available soils, topographic, precipitation, and outcrop data has been developed for predicting net infiltration to exposed and soil-covered areas of the Navajo Sandstone outcrop of southwestern Utah. The Navajo Sandstone is an important regional bedrock aquifer. The GIS model determines the net-infiltration percentage of precipitation by using an empirical equation. This relation is derived from least squares linear regression between three surficial parameters (soil coarseness, topographic slope, and downgradient distance from outcrop) and the percentage of estimated net infiltration based on environmental tracer data from excavations and boreholes at Sand Hollow Reservoir in the southeastern part of the study area. Processed GIS raster layers are applied as parameters in the empirical equation for determining net infiltration for soil-covered areas as a percentage of precipitation. This net-infiltration percentage is multiplied by average annual Parameter-elevation Regressions on Independent Slopes Model (PRISM) precipitation data to obtain an infiltration rate for each model cell. Additionally, net infiltration on exposed outcrop areas is set to 10 percent of precipitation on the basis of borehole net-infiltration estimates. Soils and outcrop net-infiltration rates are merged to form a final map. Areas of low, medium, and high potential for ground-water recharge have been identified, and estimates of net infiltration range from 0.1 to 66 millimeters per year (mm/yr). Estimated net-infiltration rates of less than 10 mm/yr are considered low, rates of 10 to 50 mm/yr are considered medium, and rates of more than 50 mm/yr are considered high. A comparison of estimated net-infiltration rates (determined from tritium data) to predicted rates (determined from GIS methods) at 12 sites in Sand Hollow and at Anderson Junction indicates an average difference of about 50 percent. Two of the predicted values were lower, five were higher, and five were within the estimated range. While such uncertainty is relatively small compared with the three order-of-magnitude range in predicted net-infiltration rates, the net-infiltration map is best suited for evaluating relative spatial distribution rather than for precise quantification of recharge to the Navajo aquifer at specific locations. An important potential use for this map is land-use zoning for protecting high net-infiltration parts of the aquifer from potential surface contamination.

Utah↗

Analysis of the effects of proposed pumping from the principal artesian aquifer, Savannah, Georgia area

A two-dimensional finite-difference model of the principal artesian aquifer in the Savannah, Georgia, area, originally developed by Counts and Krause (1976), has been expanded and refined. The model was updated and the grid redesigned to provide more current and accurate detail for ground-water resources management alternatives. Improvements in the definition of the flow system were made possible by the acquisition of additional data in the area and by recently completed regional models that include the area. The model was initially calibrated by using the estimated predevelopment potentiometric surface of 1880. The flow system under predevelopment conditions was sluggish and only 100 cubic feet per second (65 million gallons per day) flowed through the model area. It was then tested for acceptance by using the May 1980 potentiometric surface and corresponding pumping stress of approximately 85 million gallons per day in the Savannah, Georgia-Hilton Head Island, South Carolina, area. The flow through the system under 1980 conditions was about 390 cubic feet per second (250 million gallons per day) and the vertical inflow from the overlying surficial aquifer more than doubled due to formerly rejected recharge that now flows vertically into the aquifer. Calibration was accurate + or - 10 feet. The absolute error per node was 3.4 feet. A hypothetical 25-percent increase in pumpage over the entire area was used to represent a gradual growth in commercial and municipal pumpage over the next 20 to 30 years. The increase produced a maximum decline of 30 feet below the existing water level of 135 feet below sea level at the center of the cone of depression in Savannah, and a 5-foot decline at a radius of 20 miles from the center of the cone of depression. (USGS)

Georgia↗

A precipitation-runoff model for analysis of the effects of water withdrawals on streamflow, Ipswich River basin, Massachusetts

Water withdrawals from the 155-square-mile Ipswich River Basin in northeastern Massachusetts affect aquatic habitat, water quality, and recreational use of the river. To better understand the effects of these withdrawals on streamflow, particularly low flow, the Hydrological Simulation Program-FORTRAN (HSPF) was used to develop a watershed-scale precipitation-runoff model of the Ipswich River to simulate its hydrology and complex water-use patterns.An analytical solution was used to compute time series of streamflow depletions resulting from ground-water withdrawals at wells. The flow depletions caused by pumping from the wells were summed along with any surface-water withdrawals to calculate the total withdrawal along a stream reach. The water withdrawals, records of precipitation, and streamflow records on the Ipswich River at South Middleton and at Ipswich for the period 1989?93 were used to calibrate the model. Model-fit analysis indicates that the simulated flows matched observed flows over a wide range of conditions; at a minimum, the coefficient of model-fit efficiency indicates that the model explained 79 percent of the variance in the observed daily flow.Six alternative water-withdrawal and land-use scenarios were simulated with the model. Three scenarios were examined for the 1989?93 calibration period, and three scenarios were examined for the 1961?95 period to test alternative withdrawals and land use over a wider range of climatic conditions, and to compute 1-, 7-, and 30-day low-flow frequencies using a log-Pearson Type III analysis. Flow-duration curves computed from results of the 1989?93 simulations indicate that, at the South Middleton and Ipswich gaging stations, streamflows when no water withdrawals are being made are nearly identical to streamflows when no ground-water withdrawals are made. Streamflow under no water withdrawals at both stations are about an order of magnitude larger at the 99.8 percent exceedence probability than simulations with only ground-water withdrawals. Long-term simulations indicate that the differences between streamflow with no water withdrawals and average 1989?93 water withdrawals is similar to the difference between simulations for the same water-use conditions made for the 1989?93 period at both sites. The 7-day, 10-year low-flow (7Q10, a widely used regulatory statistic) at the South Middleton station was 4.1 cubic feet per second (ft3/s) with no water withdrawals and 1991 land use, 5.8 ft3/s no withdrawals and undeveloped land, and 0.54 ft3/s with average 1989?93 water withdrawals and 1991 land use. The 7Q10 at the Ipswich station was about 8.3 ft3/s for simulations with no water withdrawals for both the 1991 land use and the undeveloped land conditions, and 2.7 ft3/s for simulations with average 1989?93 water withdrawals and 1991 land use. Simulation results indicate that surface-water withdrawals have little effect on the duration and frequency of low flows, but the cumulative ground-water withdrawals substantially decrease low flows.

Water-Resources Investigations Report↗

Simulation of the shallow hydrologic system in the vicinity of Middle Genesee Lake, Wisconsin, using analytic elements and parameter estimation

Middle Genesee Lake is a ground-water flow-through lake located in a developing area in southeastern Wisconsin. Because the lake is in good connection with the shallow ground-water system, hydrologic stresses to the shallow ground-water system could adversely affect the lake system. In order to assess the effects of potential stresses on the lake, a study was completed by the U.S. Geological Survey, in cooperation with the Middle Genesee Lake Management District. The objective of the study was to identify areas that contribute ground water to the lake and estimate the hydrologic budget of the lake and hydraulic parameters affecting ground-water flow. A two-dimensional, steady-state analytic element model of the lake and surrounding area was developed using the computer code GFLOW. A parameter estimation model, UCODE, was used to optimize the calibration to measured water levels and streamflow. The calibrated model was used to evaluate the effect of three hypothetical stress scenarios on the stage of Middle Genesee Lake; the simulations were linked to UCODE, which formally incorporated parameter uncertainty into 95-percent confidence intervals around the simulated value. The scenarios included: (1) pumping from upgradient irrigation wells, (2) pumping from Lower Genesee Lake to lower lake levels, and (3) reduction in recharge resulting from development. The results of the simulations demonstrated that lake levels could be affected by hydrologic stresses in the shallow hydrologic system, with effects ranging from a 2.7 feet decline in lake stage resulting from pumping in Lower Genesee Lake to a 0.1 feet decline in lake stage from development in part of the upgradient recharge area. The range of lake stage decline increased when parameter uncertainty was included, from a decline of 3.1 feet for pumping from Lower Genesee Lake to no reduction in lake stage for the development in the recharge area. Whereas these simulated effects are within the natural variation in lake stage, they represent a systematic reduction of ground-water flow to the lake. Therefore, these hypothetical stresses are expected to establish a new, lower, baseline lake stage over which the natural variation due to climatic effects are added and subtracted.

Wisconsin↗