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 577 records · Page 32Linked to original sources
USGS2003

Methodology for estimating times of remediation associated with monitored natural attenuation

Natural attenuation processes combine to disperse, immobilize, and biologically transform anthropogenic contaminants, such as petroleum hydrocarbons and chlorinated ethenes, in ground-water systems. The time required for these processes to lower contaminant concentrations to levels protective of human health and the environment, however, varies widely between different hydrologic systems, different chemical contaminants, and varying amounts of contaminants. This report outlines a method for estimating timeframes required for natural attenuation processes, such as dispersion, sorption, and biodegradation, to lower contaminant concentrations and mass to predetermined regulatory goals in groundwater systems. The time-of-remediation (TOR) problem described in this report is formulated as three interactive components: (1) estimating the length of a contaminant plume once it has achieved a steady-state configuration from a source area of constant contaminant concentration, (2) estimating the time required for a plume to shrink to a smaller, regulatoryacceptable configuration when source-area contaminant concentrations are lowered by engineered methods, and (3) estimating the time needed for nonaqueous phase liquid (NAPL) contaminants to dissolve, disperse, and biodegrade below predetermined levels in contaminant source areas. This conceptualization was used to develop Natural Attenuation Software (NAS), an interactive computer aquifers. NAS was designed as a screening tool and requires the input of detailed site information about hydrogeology, redox conditions, and the distribution of contaminants. Because NAS is based on numerous simplifications of hydrologic, microbial, and geochemical processes, the program may introduce unacceptable errors for highly heterogeneous hydrologic systems. In such cases, application of the TOR framework outlined in this report may require more detailed, site-specific digital modeling. The NAS software may be downloaded from the Web site http://www.cee.vt.edu/NAS/ Application of NAS illustrates several general characteristics shared by all TOR problems. First, the distance of stabilization of a contaminant plume is strongly dependent on the natural attenuation capacity of particular ground-water systems. The time that it takes a plume to reach a steady-state configuration, however, is independent of natural attenuation capacity. Rather, the time of stabilization is most strongly affected by the sorptive capacity of the aquifer, which is dependent on the organic matter content of the aquifer sediments, as well as the sorptive properties of individual contaminants. As a general rule, a high sorptive capacity retards a plume.s growth or shrinkage, and increases the time of stabilization. Finally, the time of NAPL dissolution depends largely on NAPL mass, composition, geometry, and hydrologic factors, such as ground-water flow rates. An example TOR analysis for petroleum hydrocarbon NAPL was performed for the Laurel Bay site in South Carolina. About 500 to 1,000 pounds of gasoline leaked into the aquifer at this site in 1991, and the NAS simulations suggested that TOR would be on the order of 10 years for soluble and poorly sorbed compounds, such as benzene and methyl tertiary-butyl ether (MTBE). Conversely, TOR would be on the order of 40 years for less soluble, more strongly sorbed compounds, such as toluene, ethylbenzene, and xylenes (TEX). These TOR estimates are roughly consistent with contaminant concentrations observed over 10 years of monitoring at this site where benzene and MTBE concentrations were observed to decrease rapidly and are approaching regulatory maximum concentration limits, whereas toluene, ethylbenzene, and xylene concentrations decreased at a slower rate and have remained relatively high. An example TOR analysis for petroleum hydrocarbon NAPL was performed for the Laurel Bay site in South Carolina.

South Carolina↗
USGS1980

Effect of pumpage on ground-water levels as modeled in Laramie County, Wyoming

Groundwater is being extensively developed for domestic, agricultural, and industrial use in a 2,320-square mile area in Laramie County, WY., bounded approximately by Horse Creek on the north, Nebraska on the east, Colorado on the south, and pre-Tertiary outcrops on the west. Currently (1977) about 47,300 acres of land are irrigated with groundwater. Groundwater levels are declining in some areas as much as 4 feet per year. The investigation was made to provide State water administrators with data on water-level changes resulting from present (1977) groundwater withdrawals and to provide a means of predicting the future effect of groundwater development. A digital model was developed of the hydrologic system in the post-Cretaceous rocks. The ability of the model to simulate the hydrologic system was determined by comparing the water-level changes measured at 37 observation wells located in areas of irrigation pumping with the water-level changes calculated by the model for 1971-77. Comparison of the measured and calculated changes showed agreement with a root-mean-square deviation of + or - 3.6 feet with 8 feet as the maximum deviation. It is concluded that the model adequately simulates present hydrologic conditions in the post-Cretaceous rocks and may be used to predict the effect of applied stress to the system. (USGS)

Open-File Report↗
USGS2006

Simulation and particle-tracking analysis of ground-water flow near the Savannah River site, Georgia and South Carolina, 2002, and for selected ground-water management scenarios, 2002 and 2020

Ground-water flow under 2002 hydrologic conditions was evaluated in an eight-county area in Georgia and South Carolina near the Savannah River Site (SRS), by updating boundary conditions and pumping rates in an existing U.S. Geological Survey (USGS) ground-water model. The original ground-water model, developed to simulate hydrologic conditions during 1987-92, used the quasi-three-dimensional approach by dividing the Floridan, Dublin, and Midville aquifer systems into seven aquifers. The hydrogeologic system was modeled using six active layers (A2-A7) that were separated by confining units with an overlying source-sink layer to simulate the unconfined Upper Three Runs aquifer (layer A1). Potentiometric- surface maps depicting September 2002 for major aquifers were used to update, evaluate, and modify boundary conditions used by the earlier ground-water flow model. The model was updated using the USGS finite-difference code MODFLOW-2000 for mean-annual conditions during 1987-92 and 2002. The specified heads in the source-sink layer A1 were lowered to reflect observed water-level declines during the 1998-2002 drought. These declines resulted in a decrease of 12.1 million gallons per day (Mgal/d) in simulated recharge or vertical inflow to the uppermost confined aquifer (Gordon, layer A2). Although ground-water pumpage in the study area has increased by 32 Mgal/d since 1995, most of this increase (17.5 Mgal/d) was from the unconfined Upper Three Runs aquifer (source-sink layer A1) with the remaining 14.5 Mgal/d assigned to the active layers within the model (A2-A7). The simulated water budget for 2002 shows a decrease from the 1987-92 model from 1,040 Mgal/d to 1,035 Mgal/d. The decreased ground-water inflows and increased ground-water withdrawal rates reduced the simulated ground-water outflow to river cells in the active layers of the model by 43 Mgal/d. The calibration statistics for all layers of the 2002 simulation resulted in a decrease in the root mean square (RMS) of the residuals from 10.6 to 8.0 feet (ft). The residuals indicate 83.3 percent of the values for the 2002 simulation met the calibration error criteria established in the original model, whereas 88.8 percent was within the specified range for the 1987-92 simulation. Simulated ground-water outflow to the Savannah River and its tributaries during water year 2002 was 560 cubic feet per second (ft3/s), or 86 percent of the observed gain in mean-annual streamflow between streamflow gaging stations at the Millhaven, Ga., and Augusta, Ga. At Upper Three Runs Creek, simulated ground-water discharge during 2002 was 110 ft3/s, or 83 percent of the observed streamflow at two streamflow gaging stations near the SRS. These results indicate that the constructed model calibrated to 1987-92 conditions and modified for 2002 dry conditions is still representative of the hydrologic system. The USGS particle-tracking code MODPATH was used to generate advective water-particle pathlines and their associated time-of-travel based on MODFLOW simulations for 1987-92, 2002, and each of four hypothetical ground-water management scenarios. The four hypothetical ground-water management scenarios represent hydrologic conditions for (1) reported pumping for 2002 and boundary conditions for an average year; (2) reported pumping for 2002 with SRS pumping discontinued and boundary conditions for an average year; (3) projected 2020 pumping and boundary conditions for an average year; and (4) projected 2020 pumping and boundary conditions for a dry year. The MODPATH code was used in forward-tracking mode to evaluate flowpaths from areas on the SRS and in backtracking mode to evaluate further areas of previously documented trans-river flow on the Georgia side of the Savannah River. Trans-river flow is a condition in which the local head gradients might allow migration of contaminants from the SRS into the underlying aquifers and beneath the Savannah River into Georgia. More...

Georgia, South Carolina↗
USGS2018

Evaluation of the Source and Transport of High Nitrate Concentrations in Ground Water, Warren Subbasin, California

Ground water historically has been the sole source of water supply for the Town of Yucca Valley in the Warren subbasin of the Morongo ground-water basin, California. An imbalance between ground-water recharge and pumpage caused ground-water levels in the subbasin to decline by as much as 300 feet from the late 1940s through 1994. In response, the local water district, Hi-Desert Water District, instituted an artificial recharge program in February 1995 using imported surface water to replenish the ground water. The artificial recharge program resulted in water-level recoveries of as much as 250 feet in the vicinity of the recharge ponds between February 1995 and December 2001; however, nitrate concentrations in some wells also increased from a background concentration of 10 milligrams per liter to more than the U.S. Environmental Protection Agency (USEPA) maximum contaminant level (MCL) of 44 milligrams per liter (10 milligrams per liter as nitrogen). The objectives of this study were to: (1) evaluate the sources of the high-nitrate concentrations that occurred after the start of the artificial-recharge program, (2) develop a ground-water flow and solute-transport model to better understand the source and transport of nitrates in the aquifer system, and (3) utilize the calibrated models to evaluate the possible effect of a proposed conjunctive-use project. These objectives were accomplished by collecting water-level and water-quality data for the subbasin and assessing changes that have occurred since artificial recharge began. Collected data were used to calibrate the ground-water flow and solute-transport models. Data collected for this study indicate that the areal extent of the water-bearing deposits is much smaller (about 5.5 square miles versus 19 square miles) than that of the subbasin. These water-bearing deposits are referred to in this report as the Warren ground-water basin. Faults separate the ground-water basin into five hydrogeologic units: the west, the midwest, the mideast, the east and the northeast hydrogeologic units. Water-quality analyses indicate that septage from septic tanks is the primary source of the high-nitrate concentrations measured in the Warren ground-water basin. Water-quality and stable-isotope data, collected after the start of the artificial recharge program, indicate that mixing occurs between imported water and native ground water, with the highest recorded nitrate concentrations in the midwest and the mideast hydrogeologic units. In general, the timing of the increase in measured nitrate concentrations in the midwest hydrogeologic unit is directly related to the distance of the monitoring well from a recharge site, indicating that the increase in nitrate concentrations is related to the artificial recharge program. Nitrate-to-chloride and nitrogen-isotope data indicate that septage is the source of the measured increase in nitrate concentrations in the midwest and the mideast hydrogeologic units. Samples from four wells in the Warren ground-water basin were analyzed for caffeine and selected human pharmaceutical products; these analyses suggest that septage is reaching the water table. There are two possible conceptual models that explain how high-nitrate septage reaches the water table: (1) the continued downward migration of septage through the unsaturated zone to the water table and (2) rising water levels, a result of the artificial recharge program, entraining septage in the unsaturated zone. The observations that nitrate concentrations increase in ground-water samples from wells soon after the start of the artificial recharge program in 1995 and that the largest increase in nitrate concentrations occur in the midwest and mideast hydrogeologic units where the largest increase in water levels occur indicate the validity of the second conceptual model (rising water levels). The potential nitrate concentration resulting from a water-level rise in the midwest and mideast hydrogeologic units was estimated using a simple mixing-cell model. The estimated value is within the range of concentrations measured in samples from wells, further indicating the validity of the second conceptual model. A ground-water flow model and a solute-transport model were developed for the Warren ground-water basin for the period 1956-2001. MODFLOW-96 was used for the ground-water flow model and MOC3D was used for the solute-transport model. The model cell size is about 500 feet by 500 feet and the models were discretized vertically into three layers. The models were calibrated using a trial-and-error approach using water-level and nitrate-concentration data collected between 1956-2001. In order to better match the measured data, low fault hydraulic characteristic values were required, thereby compartmentalizing the ground-water basin. In addition, it was necessary to parameterize the specific yield distribution for the top model layer where unconfined ground-water conditions occur into three homogeneous zones. Separate sets of specific- yield values were needed to simulate the drawdown and subsequent water-level recovery. In addition, the calibrated natural recharge was about 83 acre-feet per year. The entrainment of unsaturated-zone septage was simulated as recharge having an associated nitrate concentration. The volume of recharge was a function of the measured water-level rise between 1994-98 and the moisture content of the unsaturated zone. The nitrate concentration of the recharge water was a weighted function of the assumed nitrate concentration in the infiltrating water associated with the overlying land use. The simulated hydraulic head and nitrate concentration results were in good agreement with the measured data indicating that the mechanism for the increase in nitrate concentrations was rising water levels entraining high-nitrate septage in the unsaturated-zone. The calibrated models were used to simulate the possible effects of a planned conjunctive-use project in the western part of the ground-water basin. The simulated project included the addition of a new recharge pond and a new extraction well. In addition, recharge at two existing recharge ponds was increased and three existing production wells were pumped, treated in a nitrate-removal facility, and used for water supply. The simulated hydraulic heads increased in the west, the mideast, and parts of the east hydrogeologic units; however, the simulated hydraulic heads decreased in the midwest and northeast hydrogeologic units. The simulated nitrate concentrations increased to above the MCL of 44 milligrams per liter (10 milligrams per liter as nitrogen) in parts of the west as a result of the increase in simulated hydraulic head. The simulated nitrate concentrations decreased in part of the midwest hydrogeologic unit as a result of the artificial recharge and pumping from the nitrate-removal wells. The simulated nitrate concentrations increased to above the MCL of 44 milligrams per liter in part of the mideast and parts of the east hydrogeologic units beneath commercial land-use areas.

California↗
USGS2008

Ground-Water Storage Change and Land Subsidence in Tucson Basin and Avra Valley, Southeastern Arizona, 1998-2002

Gravity and land subsidence were measured annually at wells and benchmarks within two networks in Tucson Basin and Avra Valley from 1998 to 2002. Both networks are within the Tucson Active Management Area. Annual estimates of ground-water storage change, ground-water budgets, and land subsidence were made based on the data. Additionally, estimates of specific yield were made at wells within the monitored region. Increases in gravity and water-level rises followed above-average natural recharge during winter 1998 in Tucson Basin. Overall declining gravity and water-level trends from 1999 to 2002 in Tucson Basin reflected general declining ground-water storage conditions and redistribution of the recent recharge throughout a larger region of the aquifer. The volume of stored ground-water in the monitored portion of Tucson Basin increased 200,000 acre-feet from December 1997 to February 1999; however, thereafter an imbalance in ground-water pumpage in excess of recharge led to a net storage loss for the monitoring period by February 2002. Ground-water storage in Avra Valley increased 70,000 acre-feet during the monitoring period, largely as a result of artificial and incidental recharge in the monitored region. The water-budget for the combined monitored regions of Tucson Basin and Avra Valley was dominated by about 460,000 acre-feet of recharge during 1998 followed by an average-annual recharge rate of about 80,000 acre-feet per year from 1999 to 2002. Above-average recharge during winter 1998, followed by average-annual deficit conditions, resulted in an overall balanced water budget for the monitored period. Monitored variations in storage compared well with simulated average-annual conditions, except for above-average recharge from 1998 to 1999. The difference in observed and simulated conditions indicate that ground-water flow models can be improved by including climate-related variations in recharge rates rather than invariable rates of average-annual recharge. Observed land-subsidence during the monitoring period was less than 1 inch except in the central part of Tucson Basin where land subsidence was about 2-3 inches. Correlations of gravity-based storage and water-level change at 37 wells were variable and illustrate the complex nature of the aquifer system. Storage and water-level variations were insufficient to estimate specific yield at many wells. Correlations at several wells were poor, inverse, or resulted in unreasonably large values of specific yield. Causes of anomalously correlated gravity and water levels include significant storage change in thick unsaturated zones, especially near major ephemeral channels, and multiple aquifers that are poorly connected hydraulically. Good correlation of storage and water-level change at 10 wells that were not near major streams where significant changes in unsaturated zone storage occur resulted in an average specific-yield value of 0.27.

Scientific Investigations Report↗
USGS1999

Evaluation of surface-water/ground-water interactions in the Santa Clara River Valley, Ventura County, California

The interactions of surface water and ground water along the Santa Clara River in Ventura County, California, were evaluated by analyzing river-discharge and water-quality data and geohydrologic information collected by the U.S. Geological Survey between 1993 and 1995 for the Piru, Fillmore, and Santa Paula subbasins. Measurements of discharge and water quality were made at multiple locations along the Santa Clara River and its tributaries at eight different time periods during different releases from Lake Piru. Geologic, hydraulic, and water-quality data were collected from three new multiple-completion ground-water monitoring wells. These data, together with data collected as part of the U.S. Geological Survey Southern California Regional Aquifer-System Analysis (RASA) study, were analyzed in order to quantify rates and locations of ground-water recharge and discharge within the river, characterize the correlation of recharge and discharge rates with ground-water conditions and reservoir releases, and better characterize the three-dimensional ground-water flow system. Analysis of the data indicates that the largest amount of ground-water recharge from the river consistently occurs in the Piru subbasin. Some ground-water recharge from the river may occur in the upper part of the Fillmore subbasin. Increases in sulfate concentrations indicate that increases in flow at the lower ends of the Piru and Fillmore subbasins result from high-sulfate ground-water discharge. Increases in flow in the lower part of the Santa Paula subbasin are not accompanied by significant sulfate increases. Several sets of regressions indicate possible correlation between net flow changes in the river and depths to ground water and release rates from Lake Piru. These statistical relations may be of use for evaluating alternative Lake Piru release strategies. Data on the stable isotopes of hydrogen and oxygen from the ground-water monitoring wells that were installed as part of this investigation were used to distinguish between zones affected by recharge from the Santa Clara River and zones affected by recharge from local precipitation. Tritium data from a new multiple-completion monitoring site indicate that near the river in the upper Santa Paula subbasin, recent (post-1950) recharge water is not present at depths greater than about 350 feet below land surface. Water-level and lithologic data from the monitoring site indicate that the river and the Shallow aquifer have only limited hydraulic connection to the underlying aquifers at this location. Water-level data from the Shallow aquifer and from an in-stream drive point were used in an analytic model to estimate hydraulic properties governing stream?aquifer interactions in the upper Santa Paula subbasin. Hydraulic conductivities in all the USGS monitoring wells were estimated on the basis of slug tests.

Water-Resources Investigations Report↗
USGS2000

Effects of ground-water withdrawals on the Rock River and associated valley aquifer, eastern Rock County, Minnesota

A better understanding of the ground-water and surface-water resources of the Rock River Valley in southwestern Minnesota was needed due to concerns surrounding future reliable sources of water for public supply. The Rock River Valley aquifer consists of a surficial sand and gravel unit that underlies the entire Rock River Valley and a buried sand and gravel unit that is present only in the vicinity of the Luverne Municipal and Airport well fields. The surficial and buried units of the aquifer are separated by a clay and till layer ranging in thickness from 1 to 38 feet. The combined maximum saturated thickness of the aquifer is 52 feet, with a median of 22 feet. The thickness of the buried unit ranges from 3 to 17 feet. Recharge to the Rock River Valley aquifer occurs primarily by infiltration of precipitation to the saturated zone (areal recharge) and by induced infiltration from the Rock River due to withdrawals from supply wells near the river. Discharge from the aquifer occurs primarily as leakage to streams and ground-water evapotranspiration. The water budget for the calibrated steady-state simulation indicated that areal recharge accounts for 38 percent of the sources of water to the Rock River Valley aquifer and leakage from streams contributes 58.7 percent. The largest discharge from the aquifer is leakage to streams, (71.1 percent). The net stream-aquifer leakage is approximately 5 cubic feet per second from the aquifer to the streams. The simulated contributing areas for the wells in the three well fields generally extend to the aquifer boundaries on the west and are generally truncated at the Rock River. The simulated transient water budget for 1996 indicated that the principal sources of water to the aquifer were as follows: (1) winter, spring, and late summer stress periods— leakage from streams and water released from storage and (2) early summer and fall stress periods—areal recharge and leakage from streams. The principal discharges from the aquifer were leakage to streams for all stress periods, ground-water evapotranspiration for the early and late summer stress periods, and addition to storage for the early summer and fall stress periods. The herbicides atrazine, alachlor, metolachlor, acetachlor, and cyanazine, and metabolites of these herbicides, occurred in concentrations of 0.05 to 11.5 micrograms per liter in the Rock River at Luverne during major runoff events following application of herbicides in the spring. Atrazine and metabolites, alachlor ESA (a metabolite of alachlor), metolachlor and metabolites, metolachlor ESA and metolachlor OA, and acetochlor metabolites acetochlor ESA and acetochlor OA, were detected at concentrations of 0.05 to 2.8 micrograms per liter in municipal supply wells less than 500 feet from the river during November 1995 through August 1997. The Rock River is the major source of the herbicides and metabolites. However, concentrations of atrazine and metabolites, alachlor ESA, metolachlor ESA, and metolachlor OA in supply wells may also reflect sources of these herbicides and metabolites in the ground-water contributing areas to the supply wells. Nitrite plus nitrate nitrogen concentrations in supply wells and in the ground-water contributing area to the Luverne Municipal well field were generally less than 1.5 milligrams per liter. Nitrite plus nitrate nitrogen concentrations of 2.4 to 8.5 milligrams per liter in the Rock River in the Rock County Rural Water well field and 14 to 18 milligrams per liter in the ground-water contributing area to the Rock County Rural Water supply wells are not having a substantial affect on nitrite plus nitrate nitrogen concentrations in most supply wells. Isotopic mixing calculations indicate that proportions of river water withdrawn from supply wells less than 500 feet from the river range from 5 to 60 percent of total withdrawals. The Rock River is a gaining stream in most reaches, but is losing water to the aquifer in the vicinity of the Luverne Municipal and Rock County Rural Water well fields, located 150 to 1,500 feet from the river. Simulated streamflow losses due to ground-water withdrawals in the well fields were approximately 2.1 cubic feet per second. Because an average of about 1.5 cubic feet per second of the water pumped by Luverne is returned to the Rock River as wastewater discharge, the net steady-state simulated streamflow loss for the study area is 0.6 cubic feet per second. The streamflow losses as a result of ground-water withdrawals are insignificant in comparison to typical streamflow, and are likely to have a measurable effect on streamflow only during low-flow conditions of less than approximately 10 cubic feet per second. Model results indicate that the additional water withdrawn by wells due to anticipated increased ground-water withdrawals was derived from a decrease in net leakage of ground water from the aquifer to the streams. The simulations indicated that the increased ground-water withdrawals and normal precipitation resulted in an increase in induced infiltration from the Rock River of 0.1 cubic feet per second for the Luverne Municipal well field and 0.3 cubic feet per second for the Rock County Rural Water well field. Maximum drawdowns ranged from 0.5 to 1.4 feet near the three well fields. For drought conditions, the simulated streamflow losses constituted approximately 30 percent and nearly 65 percent of the flows in the Rock River for the Luverne Municipal and Rock County Rural Water well fields, respectively. Maximum drawdowns ranged from 3.8 to 7.0 feet near the three well fields. Transient simulations with anticipated increased ground-water withdrawals and drought conditions indicated declines in hydraulic heads ranging from 0.2 to 0.4 feet per year in the vicinity of the three well fields, except for near the Rock River.

Minnesota↗
USGS1971

Hydrologic data for Horseshoe Lake, Arkansas and vicinity

During the summer and fall, seepage and evaporation losses from Horseshoe Lake, an oxbow or an 'old river' lake adjacent to the Mississippi River, exceed inflow to the lake, and seasonal declines of 2.5-3.0 feet in the lake level are common. In exceptionally dry years, the minimum lake level has been as much as 4 feet below the normal seasonal low. These low levels severely affect the recreational uses of the lake. Seepage and evaporation rates at Horseshoe Lake were determined from hydrologic and meteorologic data. Analysis of these data indicates that the direction of seepage is out of the lake except for a period of about 2 months in the spring, when the stage of the Mississippi River is high. The lake can be maintained at a constant level by supplementing the inflow to the lake with surface or ground water. Contributions to the lake from local drainage can be increased, but this water contains undesirable amounts of pesticides, herbicides, and plant nutrients, and the flow is insufficient to eliminate seasonal declines in the lake level. Water from r the Mississippi River can be used to maintain a given lake level, but the bacteriological quality of water from the river makes this an undesirable source of supplemental water. Water from the Quaternary alluvium contains troublesome amounts of iron, but it probably is free of pesticides, herbicides, and coliform bacteria which are commonly found in surface water. An electric-analog model was used to determine the rate at which inflow to the lake must be supplemented to maintain various lake levels. During this investigation, the lake could have been maintained very near the normal spring level by supplementing the inflow at a maximum rate of 10,600 gallons per minute. The analog model was also used to determine the effects of pumping wells on seepage. With the exception of wells near the southeast end of the lake, wells located within one-half mile of the lake would obtain more than 50 percent of their yield from the lake after pumping for 90 days.

Open-File Report↗
USGS1999

Conceptual Model and Numerical Simulation of the Ground-Water-Flow System in the Unconsolidated Sediments of Thurston County, Washington

The demand for water in Thurston County has increased steadily in recent years because of a rapid growth in population. Surface-water resources in the county have been fully appropriated for many years and Thurston County now relies entirely on ground water for new supplies of water. Thurston County is underlain by up to 2,000 feet of unconsolidated glacial and non-glacial Quaternary sediments which overlie consolidated rocks of Tertiary age. Six geohydrologic units have been identified within the unconsolidated sediments. Between 1988 and 1990, median water levels rose 0.6 to 1.9 feet in all geohydrologic units except bedrock, in which they declined 1.4 feet. Greater wet-season precipitation in 1990 (43 inches) than in 1988 (26 inches) was the probable cause of the higher 1990 water levels. Ground-water flow in the unconsolidated sediments underlying Thurston County was simulated with a computerized numerical model (MODFLOW). The model was constructed to simulate 1988 ground-water conditions as steady state. Simulated inflow to the model area from precipitation and secondary recharge was 620,000 acre-feet per year (93 percent), leakage from streams and lakes was 38,000 acre-ft/yr (6 percent), and ground water entering the model along the Chehalis River valley was 5,800 acre-ft/yr (1 percent). Simulated outflow from the model was primarily leakage to streams, springs, lakes, and seepage faces (500,000 acre-ft/yr or 75 percent of the total outflow). Submarine seepage to Puget Sound was simulated to be 88,000 acre-ft/yr (13 percent). Simulated ground-water discharge along the Chehalis River valley was simulated to be 12,000 acreft/yr (2 percent). Simulated withdrawals by wells for all purposes was 62,000 acre-ft/yr (9 percent). The numerical model was used to simulate the possible effects of increasing ground-water withdrawals by 23,000 acre-ft/yr above the 1988 rate of withdrawal. The model indicated that the increased withdrawals would come from reduced discharge to springs, seepage faces, and offshore (total of 51 percent of increased pumping) and decreased flow to rivers (46 percent). About 3 percent would come from increased leakage from rivers. Water levels would decline more than 1 foot over most of the model area, more than 10 feet over some areas, and would be at a maximum of about 35 feet. Contributing areas for water discharging at McAllister and Abbott Springs and to pumping centers near Tumwater and Lacey were estimated using a particle-tracking post-processing computer code (MODPATH) and a MODFLOW model calibrated to steady-state (1988) conditions. Water discharging at McAllister and Abbot Springs was determined to come from water entering the ground-water system at the water table in an area of about 20 square miles (mi2) to the west and south of the springs. This water is estimated to come from recharge (both precipitation and secondary) and from leakage from Lake St. Clair and several other surface-water bodies. Southeast of Lacey, about 3,800 acre-ft of ground water were pumped from five municipal wells during 1988. The source of the pumped water was determined to be an area that covers about 1.1 mi2. The water was estimated to come from recharge (both precipitation and secondary) and leakage from surface-water bodies. Along the lower Deschutes River nearly 3,900 acre-ft/yr of ground water were pumped during 1988 from 15 wells for municipal and industrial use. The calculated source of this water was an area that covers about 1.3 mi2. Within the calculated contributing area the pumped ground water comes from recharge (both precipitation and secondary) and leakage from the Deschutes River and several other surface-water bodies.

Water-Resources Investigations Report↗
USGS2001

Optimization of ground-water withdrawal at the old O-Field area, Aberdeen Proving Ground, Maryland

The U.S. Army disposed of chemical agents, laboratory materials, and unexploded ordnance at the Old O-Field landfill at Aberdeen Proving Ground, Maryland, beginning prior to World War II and continuing until at least the 1950?s. Soil, ground water, surface water, and wetland sediments in the Old O-Field area were contaminated by the disposal of these materials. The site is in the Atlantic Coastal Plain, and is characterized by a complex series of Pleistocene and Holocene sediments formed in various fluvial, estuarine, and marine-marginal hydrogeologic environments. A previously constructed transient finite-difference ground-water-flow model was used to simulate ground-water flow and the effects of a pump-and-treat remediation system designed to prevent contaminated ground water from flowing into Watson Creek (a tidal estuary and a tributary to the Gunpowder River). The remediation system consists of 14 extraction wells located between the Old O-Field landfill and Watson Creek.Linear programming techniques were applied to the results of the flow-model simulations to identify optimal pumping strategies for the remediation system. The optimal management objective is to minimize total withdrawal from the water-table aquifer, while adhering to the following constraints: (1) ground-water flow from the landfill should be prevented from reaching Watson Creek, (2) no extraction pump should be operated at a rate that exceeds its capacity, and (3) no extraction pump should be operated at a rate below its minimum capacity, the minimum rate at which an Old O-Field pump can function. Water withdrawal is minimized by varying the rate and frequency of pumping at each of the 14 extraction wells over time. This minimizes the costs of both pumping and water treatment, thus providing the least-cost remediation alternative while simultaneously meeting all operating constraints.The optimal strategy identified using this objective and constraint set involved operating 13 of the 14 extraction wells at rates ranging from 0.4 to 4.9 gallons per minute.

Water-Resources Investigations Report↗
USGS1981

Simulative models for the analysis of ground-water flow in Vekol Valley, the Waterman Wash area, and the Bosque area, Maricopa and Pina Counties, Arizona

Simulative ground-water flow models for Vekol Valley, the Waterman Wash area, and the Bosque area were developed for use in evaluating alternatives for developing a ground-water supply for the Ak-Chin Indian Community. The hydraulic properties of the basin-fill deposits used in the models were estimated primarily from aquifer tests made by the U.S. Geological Survey. Annual recharge to Vekol Valley and the Waterman Wash area is negligible in comparison to the quantity of water in storage and the quantity proposed to be pumped. The models are based on a three-dimensional, block-centered, finite-difference scheme. The Vekol Valley model was calibrated for steady-state onditions, and the Waterman Wash area model was calibrated for steady-state and transient conditions. The sensitivity of calibrated heads to changes in transmissivity was also investigated. An uncalibrated storage-depletion model was developed for the Bosque area. Simulated water levels for steady-state conditions average within 5 feet of measured values for Vekoi Valley and within 6 feet for the Waterman Wash area. Simulated water levels for transient conditions in the Waterman Wash area average within 8 feet of measured values for 15 years of analysis and within 15 feet for 24 years. Water-level declines simulated by the Waterman Wash area model average within 17 feet of those measured during the 24-year period, 1951-75.

Arizona↗
USGS1998

Steady-state simulation of ground-water flow in the Rush Springs Aquifer, western Oklahoma

A simplified steady-state ground-water flow model was prepared for the Rush Springs aquifer in western Oklahoma. A 3-kilometer square grid was established over the area containing two layers with 674 active nodes simulated in the model. The steady-state model simulation used a mean recharge rate of 3.05 x 10-4 feet per day and a hydraulic conductivity range from 0.8 to 10 feet per day. The error at each node in the model is defined as the difference between the measured and simulated water levels.The arithmetic mean error for 170 of the 674 active nodes was -0.11 feet, the absolute value mean error was 7.55 feet, and the standard deviation of the error was 10.21 feet. A net simulated recharge of 231 cubic feet per second is balanced by a discharge to drains and seeps of 190.6 cubic feet per second about 82 percent of the total recharge. Discharge to the main stem of the Washita River is about 41 cubic feet per second about 18 percent of the recharge.

Water-Resources Investigations Report↗
USGS2001

Simulation of a valley-fill aquifer system to delineate flow paths, contributing areas, and traveltime to wellfields in southwestern Broome County, New York

A valley-fill aquifer system that extends along a 14-mile reach of the Susquehanna River valley in southwestern Broome County, N.Y., is a major source of water supply to local municipalities and industries, but is highly susceptible to contamination from human activities. Protection of ground-water supplies requires accurate delineation of the areas that are the sources of water pumped by wells. A previously developed two-layer steady-state ground-water flow model of the aquifer system was upgraded with an improved method of simulating stream-aquifer interactions, then recalibrated and coupled to a particle-tracking program. Three-dimensional, ground-water flow modeling coupled with particle tracking is the most reliable method of simulating groundwater flow paths in multiaquifer systems such as this; it also allows delineation of contributing areas to well.elds. A primary advantage of three-dimensional particle-tracking analysis is that it shows the complexities of the flow paths in each aquifer. Model and particle tracking analyses indicate that groundwater frequently follows convoluted three-dimensional flow paths. The contributing areas of individual supply wells in this aquifer system each has a unique flow pattern and shape. Results of the model simulation indicate that recharge from precipitation, rivers, and tributaries contribute 35 percent, 29 percent, and 25 percent, respectively to the aquifer system and that pumpage from supply wells accounts for 67 percent of the discharge from the aquifer system. Particle-tracking results indicate that the simulated contributing areas to the 24 supply wells includes most of the valley floor.

New York↗
USGS2006

Hydrogeology and simulation of ground-water flow in the Silurian-Devonian aquifer system, Johnson County, Iowa

Bedrock of Silurian and Devonian age (termed the “Silurian-Devonian aquifer system”) is the primary source of ground water for Johnson County in east-central Iowa. Population growth within municipal and suburban areas of the county has resulted in increased amounts of water withdrawn from this aquifer and water-level declines in some areas. A 3-year study of the hydrogeology of the Silurian-Devonian aquifer system in Johnson County was undertaken to provide a quantitative assessment of ground water resources and to construct a ground-water flow model that can be used by local governmental agencies as a management tool. Johnson County is underlain by unconsolidated deposits of Quaternary age and Paleozoic-age bedrock units. The bulk of the Quaternary deposits consists of weathered and unweathered glacial till; however, shallow alluvium and buried sand and gravel deposits also are present. Six bedrock hydrogeologic units are present in Johnson County (oldest to youngest): Maquoketa confining unit, Silurian aquifer, Wapsipinicon Group (aquifer and confining unit), Cedar Valley aquifer, Upper Devonian shale confining unit, and Cherokee confining unit. Although separate aquifers and confining units are described, the Silurian- and Devonian-age units are considered as a single aquifer system. The top of the Silurian-Devonian aquifer system is considered as the top of the Cedar Valley aquifer, where present, and the base of the aquifer system is considered as the top of the Maquoketa confining unit. The hydraulic properties of the rocks that comprise the Silurian-Devonian aquifer system are highly variable as a result of the variable composition of the rocks and the presence of solution features in some of the carbonate-rock units. For the combined Silurian-Devonian aquifer system, specific capacity averages 2.1 gallons per minute per foot of drawdown, transmissivity averages about 580 feet squared per day, and hydraulic conductivity averages 8.3 feet per day. Recharge to the Silurian-Devonian aquifer system in Johnson County is predominantly from infiltration of precipitation to the bedrock. Discharge from the aquifer is primarily to municipal, industrial, and private-development wells. Reliable measurements of the amount of recharge to or discharge from the ground-water system in Johnson County, however, are not available. Altitude of the 1996 potentiometric surface ranged from more than 750 feet above the North American Vertical Datum of 1988 (NAVD88) in northern Johnson County to less than 575 feet above NAVD88 in the central part of the county. A large cone of depression within the potentiometric surface is present in the central part of the county, between Coralville and Iowa City. A large limestone quarry is located near the center of this cone of depression. Ground water generally flows from the northern and western parts of Johnson County either toward the cone of depression in the center of the county or south out of the county. Ground water also flows toward the Cedar River in the northeastern part of the county. A ground-water divide in the northeastern part of the county roughly approximates the surface-water divide between the Iowa River and Cedar River drainages. A numerical ground-water-flow model of the Silurian-Devonian aquifer system in Johnson County was used to test concepts of ground-water flow, to assess the need for additional data, and to evaluate the potential effects of anticipated increased ground-water development and drought. The 1-layer model was calibrated to average 1996 ground-water conditions, which were assumed to approximate steady-state flow conditions. The model also was used to simulate steady-state conditions for 2004, steady-state conditions using anticipated pumping rates for 2025, and potential future drought conditions. The simulated potentiometric surface generally replicated the potentiometric surface for 1996 and 2004 conditions. The calculated root mean squared error values for the 1996 and 2004 simulations were 13.6 and 18.6 feet, respectively. The mean absolute differences between measured and simulated water levels for the 1996 and 2004 simulations were about 11 and 14 feet, respectively. Total model-calculated inflow to the ground-water system for the 1996 simulation was 19.6 million gallons per day (Mgal/d), and the largest model-calculated inflow component was areal recharge (15.1 Mgal/d). Total model-calculated outflow from the ground-water system was 19.7 Mgal/d, and the largest outflow component was discharge to wells (10.5 Mgal/d). Model-calculated water-budget components for the 2004 simulation were similar to the 1996 components. Potential future steady-state conditions were simulated using anticipated 2025 pumping rates. Pumpage both for existing wells and for assumed new wells, based on anticipated population growth in the northern part of the county and for the nearby municipalities, was included in the model. Simulated 2025 pumpage was about 1.5 Mgal/d greater than simulated 2004 pumpage. Simulated steady-state ground-water levels, using anticipated 2025 pumping rates, were lower than 2004 simulated levels throughout the county, and simulated water-level declines ranged from less than 1 foot near the county boundaries to about 11 feet. Potential future drought conditions were simulated by assuming that recharge to the Silurian-Devonian aquifer system is reduced by a factor of 0.75 and that water-supply pumpage is increased by a factor of 1.25 over the anticipated 2025 pumping rates. Overall, simulated water levels for future drought conditions were greater than 5 feet lower than simulated 2004 conditions and were a maximum of about 30 feet lower in the northeastern part of the county. The greatest limitation to the model is the lack of measured or estimated water-budget components for comparison to simulated water-budget components. Because the model is only calibrated to measured water levels, and not to water-budget components, the model results are nonunique. Other model limitations include the relatively coarse grid scale, lack of detailed information on pumpage from the quarry and from private developments and domestic wells, and the lack of separate water-level data for the Silurian- and Devonian-age rocks.

Iowa↗
USGS2008

Conceptual Model of Hydrologic and Thermal Conditions of the Eastbank Aquifer System near Rocky Reach Dam, Douglas County, Washington

The Lower and Combined Aquifers of the Eastbank Aquifer system, located in a river-terrace deposit along the Columbia River near Rocky Reach Dam, Washington, are primarily recharged by the Columbia River and provide water to the Eastbank Hatchery and the regional water system servicing the cities of Wenatchee, East Wenatchee, and parts of unincorporated Chelan and Douglas Counties. In 2006, mean annual pumpage from the aquifers by the hatchery and regional water system was about 43 and 16 cubic feet per second, respectively. Reportedly, temperatures of ground water pumped by the hatchery have been increasing, thereby making water potentially too warm for salmonid fish production. An evaluation of hourly ground-water and river temperatures from January 1991 through August 2007 indicates increasing interannual trends in temperatures in most of the Lower and Combined Aquifers from 1999 through 2006 that correspond to increasing trends in the annual mean and annual maximum river temperatures during the same period of 0.07 and 0.17?C per year, respectively. There were no trends in the annual minimum river temperatures from 1999 through 2006, and there were no trends in the annual minimum, mean, and maximum river temperatures from 1991 through 1998 and from 1991 through 2007. Increases in river temperatures from 1999 through 2006 are within the natural variability of the river temperatures. Most of the Lower and Combined Aquifers reached thermal equilibrium?defined by constant time lags between changes in river temperatures and subsequent changes in ground-water temperatures?during 1991?98. The only exceptions are the Combined Aquifer north of the well field of the regional water system, which had not reached thermal equilibrium by 2006, and the Lower Aquifer west of the well fields of the hatchery and the regional water system, which reached thermal equilibrium prior to 1991. Because most of the Lower and Combined Aquifers were in thermal equilibrium from 1999 through 2006 and seasonal pumpage patterns were relatively stable, reported trends of increasing temperatures of water pumped by the hatchery well field are most likely explained by increasing trends in river temperatures. Most of the water pumped by the hatchery well field recharges in an area west to southwest of the well field about 2 months prior to the time it is pumped from the aquifer. The northern extent of the hatchery well field may pump some colder water from a bedrock depression to the north and west of the well field. The conceptual model of hydrologic and thermal conditions is supported by analyses of historical water temperatures, water-level data collected on July 18, 2007, and dissolved-constituent and bacterial concentrations in samples collected on August 20?22, 2007.

Scientific Investigations Report↗
USGS2006

The Cerrillos Uplift, the La Bajada Constriction, and Hydrogeologic Framework of the Santo Domingo Basin, Rio Grande Rift, New Mexico

The geologic, geophysical, and hydrogeologic properties of the La Bajada constriction and Santo Domingo Basin, northern New Mexico, result from tectonic and volcanic processes of the late Tertiary and Quaternary Rio Grande rift. An integrated geologic and geophysical assessment in the La Bajada constriction allows development of a geologic framework that can provide input for regional ground-water flow models. These models then can provide better estimates of future water supplies in a region that largely subsists on aquifers in Rio Grande rift basins. The combination of surface geologic investigations (stratigraphic and structural studies; chapters A, B, C, and E), airborne geophysics (aeromagnetic and time-domain electromagnetic surveys; chapters D and F), ground geophysical measurements (gravity and magnetotelluric surveys; chapters D and F), and data from the few wells in the area (chapter G) provides new constraints on the hydrogeologic framework of this area. Summary results of our investigations are synthesized in chapter G. Through-going aquifers consisting of ancestral Rio Grande axial-river sand and gravel and of coarse western-piedmont gravel form the predominant ground-water pathways through the partly buried structural trough defining the La Bajada constriction between Española and Santo Domingo Basins. Thick, clay-rich Cretaceous marine shales of low hydraulic conductivity form a pervasive regional confining unit within the Cerrillos uplift on the southeast flank of the constriction. Numerous, dominantly north-northwest-striking, intrabasin faults that project part way across the La Bajada constriction create a matrix of laterally and vertically variable hydrogeologic compartments that locally partition and deflect ground-water flow parallel to faults.

New Mexico↗
USGS2009

Reconnaissance Assessment of the Potential for Roadside Dry Wells to Affect Water Quality on the Island of Hawai'i

The County of Hawai'i Department of Public Works (DPW) uses dry wells to dispose of stormwater runoff from roads. Recently, concern has been raised that water entering the dry wells may transport contaminants to groundwater and affect the quality of receiving waters. The DPW operates 2,052 dry wells. Compiling an inventory of these dry wells and sorting it on the basis of presence or absence of urbanization in the drainage area, distance between the bottom of the dry well and the water table, and proximity to receiving waters helps identify the dry wells having greatest potential to affect the quality of receiving waters so that future studies or mitigation efforts can focus on a smaller number of dry wells. The drainage areas of some DPW dry wells encompass urbanized areas, which could be a source of contaminants. Some dry wells penetrate close to or through the water table, eliminating or substantially reducing opportunities for contaminant attenuation between the ground surface and water table. Dry wells that have drainage areas that encompass urbanization, penetrate to near the water table, and are near the coast have the highest potential to affect the quality of coastal waters (this study did not consider specific sections of coastline that may be of greater concern than others). Some DPW dry wells, including a few that have drainage areas that encompass urbanization, lie within the areas contributing recharge (ACR) to drinking-water wells. Numerical groundwater modeling studies by previous investigators indicate that water infiltrating those dry wells could eventually be pumped at drinking-water wells. Dry wells that have a high potential for affecting coastal receiving waters or drinking-water wells can be the focus of studies to further understand the effect of the dry wells on the quality of receiving waters. Possible study approaches include sampling for contaminants at the dry well and receiving water, injecting and monitoring the movement of tracers, and numerical modeling. To fully assess whether dry wells actually pose a significant contamination threat to receiving waters, results from modeling or monitoring must be compared to limits for contaminant concentration at receiving waters. These limits are usually established by the agencies tasked with protecting those waters.

Scientific Investigations Report↗
USGS2003

Simulation of Ground-Water Flow in the Irwin Basin Aquifer System, Fort Irwin National Training Center, California

Ground-water pumping in the Irwin Basin at Fort Irwin National Training Center, California resulted in water-level declines of about 30 feet from 1941 to 1996. Since 1992, artificial recharge from wastewater-effluent infiltration and irrigation-return flow has stabilized water levels, but there is concern that future water demands associated with expansion of the base may cause a resumption of water-level declines. To address these concerns, a ground-water flow model of the Irwin Basin was developed to help better understand the aquifer system, assess the long-term availability and quality of ground water, and evaluate ground-water conditions owing to current pumping and to plan for future water needs at the base. Historical data show that ground-water-level declines in the Irwin Basin between 1941 and 1996, caused the formation of a pumping depression near the pumped wells, and that recharge from the wastewater-treatment facility and disposal area caused the formation of a recharge mound. There have been two periods of water-level recovery in the Irwin Basin since the development of ground water in this basin; these periods coincide with a period of decreased pumpage from the basin and a period of increased recharge of water imported from the Bicycle Basin beginning in 1967 and from the Langford Basin beginning in 1992. Since 1992, artificial recharge has exceeded pumpage in the Irwin Basin and has stabilized water-level declines. A two-layer ground-water flow model was developed to help better understand the aquifer system, assess the long-term availability and quality of ground water, and evaluate ground-water conditions owing to current pumping and to plan for future water needs at the base. Boundary conditions, hydraulic conductivity, altitude of the bottom of the layers, vertical conductance, storage coefficient, recharge, and discharge were determined using existing geohydrologic data. Rates and distribution of recharge and discharge were determined from existing data and estimated when unavailable. Results of predictive simulations indicate that in 50 years, if artificial recharge continues to exceed pumpage in Irwin Basin, water levels could rise as much as 65 feet beneath the pumping depression, and as much as 10 feet in the wastewater-treatment facility and disposal area. Particle-tracking simulations were used to determine the pathlines and the traveltimes of water high in dissolved solids into the main pumping area. The pathlines of particles from two areas with high dissolved-solids concentrations show that in 50 years water from these areas almost reaches the nearest pumped well.

Water-Resources Investigations Report↗