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 793 records · Page 44Linked to original sources

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↗

Simulation of ground-water/surface-water flow in the Santa Clara-Calleguas ground-water basin, Ventura County, California

Ground water is the main source of water in the Santa Clara-Calleguas ground-water basin that covers about 310 square miles in Ventura County, California. A steady increase in the demand for surface- and ground-water resources since the late 1800s has resulted in streamflow depletion and ground-water overdraft. This steady increase in water use has resulted in seawater intrusion, inter-aquifer flow, land subsidence, and ground-water contamination. The Santa Clara-Calleguas Basin consists of multiple aquifers that are grouped into upper- and lower-aquifer systems. The upper-aquifer system includes the Shallow, Oxnard, and Mugu aquifers. The lower-aquifer system includes the upper and lower Hueneme, Fox Canyon, and Grimes Canyon aquifers. The layered aquifer systems are each bounded below by regional unconformities that are overlain by extensive basal coarse-grained layers that are the major pathways for ground-water production from wells and related seawater intrusion. The aquifer systems are bounded below and along mountain fronts by consolidated bedrock that forms a relatively impermeable boundary to ground-water flow. Numerous faults act as additional exterior and interior boundaries to ground-water flow. The aquifer systems extend offshore where they crop out along the edge of the submarine shelf and within the coastal submarine canyons. Submarine canyons have dissected these regional aquifers, providing a hydraulic connection to the ocean through the submarine outcrops of the aquifer systems. Coastal landward flow (seawater intrusion) occurs within both the upper- and lower-aquifer systems. A numerical ground-water flow model of the Santa Clara-Calleguas Basin was developed by the U.S. Geological Survey to better define the geohydrologic framework of the regional ground-water flow system and to help analyze the major problems affecting water-resources management of a typical coastal aquifer system. Construction of the Santa Clara-Calleguas Basin model required the compilation of geographic, geologic, and hydrologic data and estimation of hydraulic properties and flows. The model was calibrated to historical surface-water and ground-water flow for the period 1891-1993. Sources of water to the regional ground-water flow system are natural and artificial recharge, coastal landward flow from the ocean (seawater intrusion), storage in the coarse-grained beds, and water from compaction of fine-grained beds (aquitards). Inflows used in the regional flow model simulation include streamflows routed through the major rivers and tributaries; infiltration of mountain-front runoff and infiltration of precipitation on bedrock outcrops and on valley floors; and artificial ground-water recharge of diverted streamflow, irrigation return flow, and treated sewage effluent. Most natural recharge occurs through infiltration (losses) of streamflow within the major rivers and tributaries and the numerous arroyos that drain the mountain fronts of the basin. Total simulated natural recharge was about 114,100 acre-feet per year (acre-ft/yr) for 1984-93: 27,800 acre-ft/yr of mountain-front and bedrock recharge, 24,100 acre-ft/yr of valley-floor recharge, and 62,200 acre-ft/yr of net streamflow recharge. Artificial recharge (spreading of diverted streamflow, irrigation return, and sewage effluent) is a major source of ground-water replenishment. During the 1984-93 simulation period, the average rate of artificial recharge at the spreading grounds was about 54,400 acre-ft/yr, 13 percent less than the simulated natural recharge rate for streamflow infiltration within the major rivers and tributaries. Estimated recharge from infiltration of irrigation return flow on the valley floors averaged about 51,000 acre-ft/yr, and treated sewage effluent averaged about 9,000 acre-ft/yr. Artificial recharge as streamflow diversion to the spreading grounds has occurred since 1929, and treated-sewage effluent has been discharged to stream channels since 1930. Under

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↗

Evaluation of possible alternatives to lower the high water table of St. Charles Mesa, Pueblo County, Colorado

St. Charles Mesa, an upland terrace southeast of Pueblo, Colorado, has become increasingly urbanized as cultivated fields have been subdivided and converted to residential areas. In some areas, the water table in the terrace alluvial aquifer underlying St. Charles Mesa is very shallow. Bessemer Ditch, which delivers irrigation water to farms on the mesa and other areas of the lower Arkansas River Valley, traverses St. Charles Mesa along its southern side and is the principal source of recharge to the terrace alluvial aquifer. The ground-water flow system was assumed to be in a state of dynamic equilibrium (steady-state condition) for this study. A steady-state ground-water flow model of the terrace alluvial aquifer was constructed and calibrated. The model was run in transient state to evaluate possible alternatives of lowering the water table. The possible alternatives evaluated were (1) reducing areal recharge by reducing recharge to irrigated areas by 25 percent, (2) lining Bessemer Ditch from (a) Aspen Street to 21st Lane; (b) Aspen Street to 23rd Lane; (c) Aspen Street to 25th Lane; and (d) Aspen Street to Nicholson Road, (3) installing two drains at a depth of 10 feet below land surface upgradient from the high water table areas, and (4) installing 22 dewatering wells within the high water table areas, each pumping at 80 gallons per minute. All alternatives evaluated were at least partly effective in lowering the water table. As the simulated extent of Bessemer Ditch lining was increased, the extent and magnitude of simulated water-table declines also increased. The maximum simulated declines in the water table were 3 feet when simulated areal recharge to irrigated areas was reduced by 25 percent, 29 feet when lining of Bessemer Ditch was simulated from Aspen Street to Nicholson Road, 6.8 feet when two drains were simulated at 10-foot depth, and 14.4 feet when 22 dewatering wells, each pumping at 80 gallons per minute, were simulated. Lining Bessemer Ditch from Aspen Street to 25th Lane and from Aspen Street to Nicholson Road both resulted in water-table declines of at least 5 feet throughout most of the area. Except for reducing recharge to irrigated areas and installation of the two drains, all the alternatives evaluated probably would lower the water table enough to diminish the ground-water supply available for at least some existing wells.

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↗

Directions and rates of ground-water movement in the vicinity of Kesterson Reservoir, San Joaquin Valley, California

A three-dimensional groundwater flow model was used to simulate groundwater flow for a 124 sq mi area in the vicinity of Kesterson Reservoir in the San Joaquin Valley, California. Available data were used to calculate a probable range of groundwater flow rates, but calibration and sensitivity analysis were not done for this model. Flow directions, as inferred from measured groundwater levels and simulated hydraulic heads from all model simulations, indicate that regional groundwater flow is from the south to the north. Kesterson Reservoir acts as a recharge mound superimposed on the regional-flow system. Groundwater moves in the horizontal and vertical direction away from Kesterson Reservoir. Mud and Salt Sloughs act as groundwater discharge areas. Simulated groundwater flow from Kesterson Reservoir did not flow beyond these sloughs. Groundwater from west of Mud Slough seems to flow west toward Los Banos Creek and east toward Mud Slough. Groundwater that travels toward Salt Slough from Kesterson Reservoir probably is lost by evapotranspiration near the surface before reaching Salt Slough. Groundwater between Salt Slough and the San Joaquin River seems to flow north and toward Salt Slough and the San Joaquin River. The canals and duck ponds generally act as sources of groundwater recharge. A method was developed for determining flow directions and distance traveled in three dimensions for discrete time increments using simulated groundwater fluxes. Simulated average horizontal pore velocities away from Kesterson range less than 0.01 to 140 ft/year. The simulated average vertical pore velocities range from 0.01 to 14.7 ft/year. (Author 's abstract)

Water-Resources Investigations Report↗

Ground-water and surface-water flow and estimated water budget for Lake Seminole, southwestern Georgia and northwestern Florida

Lake Seminole is a 37,600-acre impoundment formed at the confluence of the Flint and Chattahoochee Rivers along the Georgia?Florida State line. Outflow from Lake Seminole through Jim Woodruff Lock and Dam provides headwater to the Apalachicola River, which is a major supply of freshwater, nutrients, and detritus to ecosystems downstream. These rivers,together with their tributaries, are hydraulically connected to karst limestone units that constitute most of the Upper Floridan aquifer and to a chemically weathered residuum of undifferentiated overburden. The ground-water flow system near Lake Seminole consists of the Upper Floridan aquifer and undifferentiated overburden. The aquifer is confined below by low-permeability sediments of the Lisbon Formation and, generally, is semiconfined above by undifferentiated overburden. Ground-water flow within the Upper Floridan aquifer is unconfined or semiconfined and discharges at discrete points by springflow or diffuse leakage into streams and other surface-water bodies. The high degree of connectivity between the Upper Floridan aquifer and surface-water bodies is limited to the upper Eocene Ocala Limestone and younger units that are in contact with streams in the Lake Seminole area. The impoundment of Lake Seminole inundated natural stream channels and other low-lying areas near streams and raised the water-level altitude of the Upper Floridan aquifer near the lake to nearly that of the lake, about 77 feet. Surface-water inflow from the Chattahoochee and Flint Rivers and Spring Creek and outflow to the Apalachicola River through Jim Woodruff Lock and Dam dominate the water budget for Lake Seminole. About 81 percent of the total water-budget inflow consists of surface water; about 18 percent is ground water, and the remaining 1 percent is lake precipitation. Similarly, lake outflow consists of about 89 percent surface water, as flow to the Apalachicola River through Jim Woodruff Lock and Dam, about 4 percent ground water, and about 2 percent lake evaporation. Measurement error and uncertainty in flux calculations cause a flow imbalance of about 4 percent between inflow and outflow water-budget components. Most of this error can be attributed to errors in estimating ground-water discharge from the lake, which was calculated using a ground-water model calibrated to October 1986 conditions for the entire Apalachicola?Chattahoochee?Flint River Basin and not just the area around Lake Seminole. Evaporation rates were determined using the preferred, but mathematically complex, energy budget and five empirical equations: Priestley-Taylor, Penman, DeBruin-Keijman, Papadakis, and the Priestley-Taylor used by the Georgia Automated Environmental Monitoring Network. Empirical equations require a significant amount of data but are relatively easy to calculate and compare well to long-term average annual (April 2000?March 2001) pan evaporation, which is 65 inches. Calculated annual lake evaporation, for the study period, using the energy-budget method was 67.2 inches, which overestimated long-term average annual pan evaporation by 2.2 inches. The empirical equations did not compare well with the energy-budget method during the 18-month study period, with average differences in computed evaporation using each equation ranging from 8 to 26 percent. The empirical equations also compared poorly with long-term average annual pan evaporation, with average differences in evaporation ranging from 3 to 23 percent. Energy budget and long-term average annual pan evaporation estimates did compare well, with only a 3-percent difference between estimates. Monthly evaporation estimates using all methods ranged from 0.7 to 9.5 inches and were lowest during December 2000 and highest during May 2000. Although the energy budget is generally the preferred method, the dominance of surface water in the Lake Seminole water budget makes the method inaccurate and difficult to use, because surface water makes up m

Florida, Georgia↗

Data-base system for northern Midwest regional aquifer-system analysis

The U.S. Geological Survey is conducting a study of the Cambrian and Ordovician aquifer system of the northern Midwest as part of a national series of Regional Aquifer-Systems Analysis (RASA). An integral part of this study will be a simulation of the ground-water flow regime using the Geological Survey's three-dimensional finite-difference model. The first step in the modeling effort is the design and development of a systematic set of processes to facilitate the collection, evaluation, manipulation, and use of large quantities of information. A computerized data-base system to accomplish these goals has been completed for the northern Midwest RASA. The input to the data-base system consists of either point values or contour maps of the hydrogeologic data required by the model. Digital samples of contoured surfaces are obtained by machine digitization. Uniformly spaced model-node values are then computed from the discrete data by two-dimensional interpolation. The interpolator uses the regional (long wavelength) characteristics of the data to compute a surface with minimum total curvature. Local (short wavelength) data are merged with the regional surface to produce the model-node values at the desired spacing. Nonuniformly spaced nodal values may be obtained by fitting two-dimensional polynomials (bicubic splines) to surfaces formed by uniformly spaced points and solving for the surface values at the node locations. The uniformly or nonuniformly spaced node values constitute the output of the data-base system and in turn are the input arrays of data values for the ground-water flow model. Management of the data base is facilitated by forming a data file for each model input parameter and for other hydrogeologic data used in the study. The data files are subdivided into elements consisting of individual model-layer arrays. In this form the data base can be readily accessed and edited. In addition to the interpolation programs, software components of the system include programs: (1) to transform geographic coordinates on a Lambert conformal conic projection to cartesian coordinates and vice versa, (2) to machine contour gridded data, (3) to expedite manipulation and editing of data, and (4) to perform various interarray computations. The data-base system has the following attributes: (1) manual handling of data is minimized and machine handling of data is maximized, (2) given unequally spaced point data over the extent of a study area, model input arrays for any reasonable uniform node spacing can be rapidly computed, (3) accuracy of computed node values are generally compatible with accuracy and spatial distribution of point data, (4) a relatively large class of nonuniformly spaced node configurations can be computed, (5) data within data-base files can be easily accessed and readily edited, and (6) the occurrence of data processing errors at various stages of data-base generation is monitored by machine contouring the computed grids. Functioning of the data-base system is illustrated by the sequence of steps required to produce a model-input array of transmissivity, given raw geologic and hydraulic conductivity data.

Water-Resources Investigations Report↗

Analyses of Water-Level Differentials and Variations in Recharge between the Surficial and Upper Floridan Aquifers in East-Central and Northeast Florida

Continuous (daily) water-level data collected at 29 monitoring-well cluster sites were analyzed to document variations in recharge between the surficial (SAS) and Floridan (FAS) aquifer systems in east-central and northeast Florida. According to Darcy's law, changes in the water-level differentials (differentials) between these systems are proportional to changes in the vertical flux of water between them. Variations in FAS recharge rates are of interest to water-resource managers because changes in these rates affect sensitive water resources subject to minimum flow and water-level restrictions, such as the amount of water discharged from springs and changes in lake and wetland water levels. Mean daily differentials between 2000-2004 ranged from less than 1 foot at a site in east-central Florida to more than 114 feet at a site in northeast Florida. Sites with greater mean differentials exhibited lower percentage-based ranges in fluctuations than did sites with lower mean differentials. When averaged for all sites, differentials (and thus Upper Floridan aquifer (UFA) recharge rates) decreased by about 18 percent per site between 2000-2004. This pattern can be associated with reductions in ground-water withdrawals from the UFA that occurred after 2000 as the peninsula emerged from a 3-year drought. Monthly differentials exhibited a well-defined seasonal pattern in which UFA recharge rates were greatest during the dry spring months (8 percent above the 5-year daily mean in May) and least during the wetter summer/early fall months (4 percent below the 5-year daily mean in October). In contrast, differentials exceeded the 5-year daily mean in all but 2 months of 2000, indicative of relatively high ground-water withdrawals throughout the year. On average, the UFA received about 6 percent more recharge at the project sites in 2000 than between 2000-2004. No statistically significant correlations were detected between monthly differentials and precipitation at 27 of the 29 sites between 2000-2004. For longer periods of record, double-mass plots of differentials and precipitation indicate the UFA recharge rate increased by about 34 percent at a site in west Orange County between the periods of 1974-1983 and 1983-2004. Given the absence of a trend in rainfall, the increase can likely be attributed to ground-water development. At a site in south Lake County, double-mass plots indicate that dredging of the Palatlakaha River and other nearby drainage improvements may have reduced recharge rates to the UFA by about 30 percent from the period between 1960-1965 to 1965-1970. Water-level differentials were positively correlated with land-surface altitude. The correlation was particularly strong for the 11 sites located in physiographically-defined ridge areas (coefficient of determination (R2) = 0.89). Weaker yet statistically significant negative correlations were detected between differentials and the model-calibrated leakance and thickness of the intermediate confining unit (ICU). Recharge to the UFA decreased by about 14 percent at the Charlotte Street monitoring-well site in Seminole County between 2000-2004. The decrease can be attributed to a reduction in nearby pumpage, from 57 to 49 million gallons per day over the 5-year period, with a subsequent recovery in UFA water levels that exceeded those in the SAS. Differentials at Charlotte were influenced by system memory of both precipitation and pumpage. While not statistically correlated with monthly precipitation, monthly differentials were well correlated with the 9-month moving average of precipitation. Similarly, differentials were best correlated with the 2-month moving average of pumpage. The polynomial function that quantifies the correlation between differentials and the 2-month moving average of pumpage indicates that, in terms of UFA recharge rates, the system was closer to a steady-state condition in 2000 when pumpage rates were high, than from 2001-2004 when p

Scientific Investigations Report↗

Digital-simulation and projection of head changes in the Potomac-Raritan-Magothy aquifer system, coastal plain, New Jersey

The Potomac-Raritan-Magothy aquifer system of Cretaceous age, which is the principal source of water to the major population and industrial centers in the Coastal Plain of New Jersey, has undergone continuous and widespread reduction in head. The reduced head, already below sea level throughout most of the aquifer system, in conjunction with encroachment of salty water toward centers of pumping, threatens the continued use of the aquifer as a source of freshwater in the area. A single layer, two-dimensional finite difference digital model was used to simulate the response of the aquifer system to pumping stresses during the 18-year period, 1956-73. Model simulations were based on close agreement between (1) observed and calculated heads and head trends for 10 observation wells during the period, 1956-73; and (2) the computed nonpumping steady-state potentiometric surface and a potentiometric surface based on early water-level observations (1900±) solutions. In addition, the hydrologic budget estimated by the model appears to be reasonable for the transient (1956-73) and the steady-state (1900±) solutions. The model was used to compute projected potentiometric heads and trends to the year 2000. Three sets of conditions, all using the 1973 distribution of pumping centers, were simulated. The conditions are: (1) no increase in ground-water extractions; (2) continued growth in ground-water extractions at the rate of 1.7 and 3 percent annually; and (3) continued growth in ground-water extractions at the rate of 3 percent annually, in conjunction with the activitation of a freshwater head barrier in the fresh-salty transition zone. Under the first set of conditions, further head reduction would cease over very large regions within two years. Under the second set of conditions involving a 3 percent growth rate similar to that experienced during the simulation period, the broad cone of depression already encompassing most of the New Jersey Coastal Plain would broaden and deepen. Heads would range from 60 to 160 feet below National Geodetic Vertical Datum of 1929. The reduction of head after 1973 would approach 90 feet in some areas. The resultant steeper hydraullic gradients would accelerate the rate of movement of salty ground water toward the pumping centers. A freshwater head barrier could be established in the transition zone to prevent migration of salty ground water across a 35-mile stretch in Gloucester, Camden, and Burlington Counties. A line of injection wells would be required, with total rates of injection to the head barrier ranging from about 56 cubic feet per second in 1984 to about 95 cubic feet per second in 2000. Barrier recharge rates would be equivalent to about 20 percent of the ground water pumped in the fine-grid area in any particular year for a 3 percent increase in extractions.

New Jersey↗

Ground-water recharge in the arid and semiarid southwestern United States

Ground-water recharge in the arid and semiarid southwestern United States results from the complex interplay of climate, geology, and vegetation across widely ranging spatial and temporal scales. Present-day recharge tends to be narrowly focused in time and space. Widespread water-table declines accompanied agricultural development during the twentieth century, demonstrating that sustainable ground-water supplies are not guaranteed when part of the extracted resource represents paleorecharge. Climatic controls on ground-water recharge range from seasonal cycles of summer monsoonal and winter frontal storms to multimillennial cycles of glacial and interglacial periods. Precipitation patterns reflect global-scale interactions among the oceans, atmosphere, and continents. Large-scale climatic influences associated with El Niño and Pacific Decadal Oscillations strongly, but irregularly, control weather in the study area, so that year-to-year variations in precipitation and ground-water recharge are large and difficult to predict. Proxy data indicate geologically recent periods of naturally occurring multidecadal droughts unlike any in the modern instrumental record. Any anthropogenically induced climate change will likely reduce ground-water recharge through diminished snowpack at higher elevations. Future changes in El Niño and monsoonal patterns, both crucial to precipitation in the study area, are highly uncertain in current models. Current land-use modifications influence ground-water recharge through vegetation, irrigation, and impermeable area. High mountain ranges bounding the study area—the San Bernadino Mountains and Sierra Nevada to the west, and the Wasatch and southern Colorado Rocky Mountains to the east—provide external geologic controls on ground-water recharge. Internal geologic controls stem from tectonic processes that led to numerous, variably connected alluvial-filled basins, exposure of extensive Paleozoic aquifers in mountainous recharge areas, and distinct modes of recharge in the Colorado Plateau and Basin and Range subregions. The chapters in this professional paper present (first) an overview of climatic and hydrogeologic framework (chapter A), followed by a regional analysis of ground-water recharge across the entire study area (chapter B). These are followed by an overview of site-specific case studies representing different subareas of the geographically diverse arid and semiarid southwestern United States (chapter C); the case studies themselves follow in chapters D–K. The regional analysis includes detailed hydrologic modeling within the framework of a high-resolution geographic-information system (GIS). Results from the regional analysis are used to explore both the distribution of ground-water recharge for mean climatic conditions as well as the influence of two climatic patterns—the El Niño-Southern Oscillation and Pacific Decadal Oscillation—that impart a high degree of variability to the hydrologic cycle. Individual case studies employ a variety of geophysical and geochemical techniques to investigate recharge processes and relate the processes to local geologic and climatic conditions. All of the case studies made use of naturally occurring tracers to quantify recharge. Thermal and geophysical techniques that were developed in the course of the studies are presented in appendices. The quantification of ground-water recharge in arid settings is inherently difficult due to the generally low amount of recharge, its spatially and temporally spotty nature, and the absence of techniques for directly measuring fluxes entering the saturated zone from the unsaturated zone. Deep water tables in arid alluvial basins correspond to thick unsaturated zones that produce up to millennial time lags between changes in hydrologic conditions at the land surface and subsequent changes in recharge to underlying ground water. Recent advances in physical, chemical, isotopic, and modeling techniques have fostered new types of recharge assessments. Chemical and isotopic techniques include an increasing variety of environmental tracers that are useful and robust. Physically based techniques include the use of heat as a tracer and computationally intensive geophysical imaging tools for characterizing hydrologic conditions in the unsaturated zone. Modeling-based techniques include spatially distributed water-budget computations using high-resolution remotely sensed and ground-based geographic data. Application of these techniques to arid and semiarid settings in the southwestern United States reveals distinct patterns of recharge corresponding to geologic setting, climatic and vegetative history, and land use. Analysis of recharge patterns shows that large expanses of alluvial basin floors are drying out under current climatic conditions, with little to no recharge to underlying ground water. Ground-water recharge occurs mainly beneath upland catchments in which thin soils overlie permeable bedrock, ephemeral channels in which flow may average only several hours per year, and active agricultural areas. The chapters in this professional paper represent a coordinated attempt to develop a better understanding of one of the Nation's most critical yet difficult-to-quantify renewable resources.

Professional Paper↗

Hydrology, nutrient concentrations, and nutrient yields in nearshore areas of four lakes in northern Wisconsin, 1999-2001

The effects of shoreline development on water quality and nutrient yields in nearshore areas of four lakes in northern Wisconsin were investigated from October 1999 through September 2001. The study measured surface runoff and ground-water flows from paired developed (sites containing lawn, rooftops, sidewalks, and driveways) and undeveloped (mature and immature woods) catchments adjacent to four lakes in northern Wisconsin. Water samples from surface runoff and ground water were collected and analyzed for nutrients. Coupled with water volumes, loads and subsequent yields of selected constituents were computed for developed and undeveloped catchments. The median runoff from lawn surfaces ranged from 0.0019 to 0.059 inch over the catchment area. Median surface runoff estimates from the wooded catchments were an order of magnitude less than those from the lawn catchments. The increased water volumes from the lawn catchments resulted in greater nutrient loads and subsequent annual nutrient yields from the developed sites. Soil temperature and soil moisture were measured at two sites with mixed lawn and wooded areas. At both of these sites, the area covered with a lawn commonly was warmer than the wooded area. No consistent differences in soil moisture were found. A ground-water model was constructed to simulate the local flow systems at two of the paired catchments. Model simulations showed that much of the ground water delivered to the lake originated from distant areas that did not contribute runoff directly to the lake. Surface runoff and ground-water nutrient concentrations from the lawn and wooded catchments did not have apparent patterns. Some of the median concentrations from lawns were significantly different (at the 0.05 significance level) from those at wooded catchments. Water wells and piezometers were sampled for chemical analyses three times during the study period. Variability in the shallow ground-water chemistry over time in the lawn samples was larger than samples from the wooded areas and upgradient wells. Median nutrient yields in surface runoff from lawns always were greater than those from the wooded catchments. Runoff volumes were the most important factor in determining whether lawns or wooded catchments contribute more nutrients to the lake. The ground-water system had appreciable nutrient concentrations, and are likely an important pathway for nutrient transport to the lake. The nitrate plus nitrite nitrogen and total phosphorus yields to the ground-water system from a lawn catchment were approximately 3 to 4 times greater than those from the wooded catchment. There was no difference in the yields of dissolved inorganic phosphorus to the ground-water system from the lawn and wooded catchments. Study results demonstrate that choosing the appropriate landscape position for locating lawns in sloped areas (specifically, slopes that do not terminate at the lake or areas with intervening flat or buffer zones between lawn and lake) can help reduce the adverse effect of lawns on the shallow ground water and, ultimately, the lake. Additional information would be needed to extrapolate these results to a large drainage area of a lake.

Wisconsin↗