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 811 records · Page 45Linked to original sources
USGS1993

Simulation of changes in water levels and ground-water flow in response to water-use alternatives in the Mud Lake area, eastern Snake River plain, eastern Idaho

Water users rely on surface and ground water to irrigate crops and maintain wildlife refuges in the 2,200-square-mile Mud Lake study area. Water managers need the ability to evaluate the effects of water-use changes on the future supply of surface and ground water. A five-layer, three-dimensional, finite-difference, numerical ground-water flow model, calibrated to assumed 1980 steady-state hydrologic conditions, was used to evaluate potential effects of seven water-use alternatives on ground-water levels and on losses from and gains to streams and lakes. The model was used to simulate steady-state water levels and ground-water flow for average 1980-90 hydrologic conditions and for seven water-use alternatives that represented changes from average 1980-90 conditions. Five alternatives represented reduced withdrawals from five different sets of wells, the sixth represented increased withdrawals in areas that could potentially support additional irrigation development, and the seventh represented reduced recharge in part of the study area where change from subirrigation to sprinkler irrigation is taking place. Simulated results from each alternative were compared with results for average 1980-90 conditions. Among the five water-use alternatives in which withdrawals from wells were reduced, simulated water levels were 0.1 to 40 feet higher than average 1980-90 conditions. Simulated stream and lake losses were as much as 4,700 acre-feet less and simulated gains were as much as 19,000 acre-feet greater in response to simulated water-level rises. Simulated underflow into the study area was as much as 8,200 acre-feet less and simulated underflow out of the study area was as much as 91,000 acre-feet greater. Simulated water-level declines were as great as 15 feet for the sixth alternative (increased withdrawals) and 10 feet for the seventh (reduced recharge). Simulated stream and lake losses were as much as 5,700 acre-feet greater and simulated gains were as much as 37,000 acre-feet less for stream and lake segments due to simulated water-level declines. Simulated underflow into the study area was as much as 7,200 acre-feet greater and simulated underflow out of the study area was as much as 23,000 acre-feet less.

Idaho↗
USGS1996

Geochemical and isotopic composition of ground water, with emphasis on sources of sulfate, in the upper Floridan aquifer and intermediate aquifer system in southwest Florida

In southwest Florida, sulfate concentrations in water from the Upper Floridan aquifer and overlying intermediate aquifer system are commonly above 250 milligrams per liter (the drinking water standard), particularly in coastal areas. Possible sources of sulfate include dissolution of gypsum from the deeper part of the Upper Floridan aquifer or the middle confining unit, saltwater in the aquifer, and saline waters from the middle confining unit and Lower Floridan aquifer. The sources of sulfate and geochemical processes controlling ground-water composition were evaluated for the Peace and Myakka River Basins and adjacent coastal areas of southwest Florida. Samples were collected from 63 wells and a saline spring, including wells finished at different depth intervals of the Upper Floridan aquifer and intermediate aquifer system at about 25 locations. Sampling focused along three ground-water flow paths (selected based on a predevelopment potentiometric-surface map). Ground water was analyzed for major ions, selected trace constituents, dissolved organic carbon, and stable isotopes (delta deuterium, oxygen-18, carbon-13 of inorganic carbon, and sulfur-34 of sulfate and sulfide); the ratio of strontium-87 to strontium-86 was analyzed for waters along one of the flow paths. Chemical and isotopic data indicate that dedolomitization reactions (gypsum and dolomite dissolution and calcite precipitation) control the chemical composition of water in the Upper Floridan aquifer in inland areas. This is confirmed by mass-balance modeling between wells in the shallowest interval in the aquifer along the flow paths. However, gypsum occurs deeper in the aquifer than these wells. Upwelling of sulfate-rich water that previously dissolved gypsum in deeper parts of the aquifer is a more likely source of sulfate than gypsum dissolution in shallow parts of the aquifer. This deep ground water moves to shallower zones in the aquifer discharge area. Saltwater from the Upper Floridan aquifer has not dissolved significant amounts of gypsum compared to fresher water in the aquifer. This is consistent with a shallow seawater source for the saltwater, rather than a deeper source from the underlying middle confining unit or Lower Floridan aquifer, which would have elevated sulfate concentrations. Ion exchange and dolomitization may be important reactions for saltwater in the aquifer. According to geochemical modeling, the freshwater end member for water in the saltwater mixing zone in the southwestern part of the study area is not upgradient water from the Upper Floridan aquifer that dissolved gypsum. Instead, this water appears to be isolated from the regional freshwater flow system and may be part of a more localized flow system. The chemical and isotopic composition of water in the intermediate aquifer system is controlled by differences in extent of reactions with aquifer minerals, upward leakage from the Upper Floridan aquifer, and saltwater mixing. In inland areas, water generally is characterized by relatively low sulfate concentrations (less than 250 milligrams per liter) and differences in extent of carbonate mineral dissolution. Some inland waters have elevated chloride concentrations, which may be related to evaporation prior to recharge. In coastal Sarasota County and in isolated inland areas, water from the intermediate aquifer system has high sulfate concentrations characteristic of dedolomitization waters from the Upper Floridan aquifer. The chemical and isotopic composition of these waters is controlled by upward leakage from the Upper Floridan aquifer, which naturally occurs in the discharge area but may be locally enhanced by pumping or interconnection of wells open to both aquifer systems. In western Charlotte County, the waters are dominated by sodium and chloride, and their compositions are consistent with mixing between saltwater and inland intermediate aquifer system water that has not been influenced by discharge from the

Florida↗
USGS2000

Sustainable-yield estimation for the Sparta Aquifer in Union County, Arkansas

Options for utilizing alternative sources of water to alleviate overdraft from the Sparta aquifer and ensure that the aquifer can continue to provide abundant water of excellent quality for the future are being evaluated by water managers in Union County. Sustainable yield is a critical element in identifying and designing viable water supply alternatives. With sustainable yield defined and a knowledge of total water demand in an area, any unmet demand can be calculated. The ground-water flow model of the Sparta aquifer was used to estimate sustainable yield using an iterative approach. The Sparta aquifer is a confined aquifer of regional importance that comprises a sequence of unconsolidated sand units that are contained within the Sparta Sand. Currently, the rate of withdrawal in some areas greatly exceeds the rate of recharge to the aquifer and considerable water-level declines have occurred. Ground-water flow model results indicate that the aquifer cannot continue to meet growing water-use demands indefinitely and that water levels will drop below the top of the primary producing sand unit in Union County (locally termed the El Dorado sand) by 2008 if current water-use trends continue. Declines of that magnitude will initiate dewatering of the El Dorado sand. The sustainable yield of the aquifer was calculated by targeting a specified minimum acceptable water level within Union County and varying Union County pumpage within the model to achieve the target water level. Selection of the minimum target water level for sustainable-yield estimation was an important criterion for the modeling effort. In keeping with the State Critical Ground-Water Area designation criteria and the desire of water managers in Union County to improve aquifer conditions and bring the area out of the Critical Ground-Water Area designation, the approximate altitude of the top of the Sparta Sand in central Union County was used as the minimum water level target for estimation of sustainable yield in the county. A specific category of sustainable yield? stabilization yield, reflecting the amount of water that the aquifer can provide while maintaining current water levels? also was determined and provides information for short-term management. The top of the primary producing sand unit (the El Dorado sand) was used as the minimum water-level target for estimating stabilization yield in the county because current minimum water levels in central Union County are near the top of the El Dorado sand. Model results show that withdrawals from the Sparta aquifer in Union County must be reduced to 28 percent of 1997 values to achieve sustainable yield and maintain water levels at the top of the Sparta Sand if future pumpage outside of Union County is assumed to increase at the rate observed from 1985-1997. Results of the simulation define a very large current unmet demand and represent a substantial reduction in the county?s current dependence upon the aquifer. If future pumpage outside of Union County is assumed to increase at double the rate observed from 1985-1997, withdrawals from the Sparta aquifer in Union County must be reduced to 25 percent of 1997 values to achieve sustainable yield. Withdrawals from the Sparta aquifer in Union County must be reduced to about 88 to 91 percent (depending on pumpage growth outside of the county) of 1997 values to stabilize water levels at the top of the El Dorado sand. This result shows that 1997 rate of withdrawal in the county is considerably greater than the rate needed to halt the rapid decline in water levels.

Water-Resources Investigations Report↗
USGS2004

Two-dimensional hydrodynamic simulation of surface-water flow and transport to Florida Bay through the Southern Inland and Coastal Systems (SICS)

Successful restoration of the southern Florida ecosystem requires extensive knowledge of the physical characteristics and hydrologic processes controlling water flow and transport of constituents through extremely low-gradient freshwater marshes, shallow mangrove-fringed coastal creeks and tidal embayments, and near-shore marine waters. A sound, physically based numerical model can provide simulations of the differing hydrologic conditions that might result from various ecosystem restoration scenarios. Because hydrology and ecology are closely linked in southern Florida, hydrologic model results also can be used by ecologists to evaluate the degree of ecosystem restoration that could be achieved for various hydrologic conditions. A robust proven model, SWIFT2D, (Surface-Water Integrated Flow and Transport in Two Dimensions), was modified to simulate Southern Inland and Coastal Systems (SICS) hydrodynamics and transport conditions. Modifications include improvements to evapotranspiration and rainfall calculation and to the algorithms that describe flow through coastal creeks. Techniques used in this model should be applicable to other similar low-gradient marsh settings in southern Florida and elsewhere. Numerous investigations were conducted within the SICS area of southeastern Everglades National Park and northeastern Florida Bay to provide data and parameter values for model development and testing. The U.S. Geological Survey and the National Park Service supported investigations for quantification of evapotranspiration, vegetative resistance to flow, wind-induced flow, land elevations, vegetation classifications, salinity conditions, exchange of ground and surface waters, and flow and transport in coastal creeks and embayments. The good agreement that was achieved between measured and simulated water levels, flows, and salinities through minimal adjustment of empirical coefficients indicates that hydrologic processes within the SICS area are represented properly in the SWIFT2D model, and that the spatial and temporal resolution of these processes in the model is adequate. Sensitivity analyses were conducted to determine the effect of changes in boundary conditions and parameter values on simulation results, which aided in identifying areas of greatest uncertainty in the model. The parameter having the most uncertainty (most in need of further field study) was the flow coefficient for coastal creeks. Smaller uncertainties existed for wetlands frictional resistance and wind. Evapotranspiration and boundary inflows indicated the least uncertainty as determined by varying parameters used in their formulation and definition. Model results indicated that wind was important in reversing coastal creek flows. At Trout Creek (the major tributary connecting Taylor Slough wetlands with Florida Bay), flow in the landward direction was not simulated properly unless wind forcing was included in the simulation. Simulations also provided insight into the major influence that wind has on salinity mixing along the coast, the varying distribution of wetland flows at differing water levels, and the importance of topography in controlling flows to the coast. Slight topographic variations were shown to highly influence the routing of water. A multiple regression analysis was performed to relate inflows at the northern boundary of Taylor Slough bridge to a major pump station (S-332) north of the SICS model area. This analysis allows Taylor Slough bridge boundary conditions to be defined for the model from operating scenarios at S-332, which should facilitate use of the SICS model as an operational tool.

Florida↗
USGS2004

Simulated effects of pumping and drought on ground-water levels and the freshwater-saltwater interface on the north fork of Long Island, New York

Ground water is the sole source of freshwater on the North Fork of Long Island. Future demands for the limited freshwater supply during a prolonged drought could cause drawdowns that induce saltwater intrusion and render the supply unusable. The freshwater system on the North Fork contains several localized, hydraulically isolated aquifers bounded by salty water. The need for information on the ability of these aquifers to meet future demands prompted a 4-year study to develop a ground-water flow model to simulate several proposed pumping scenarios under long-term average conditions and during a hypothetical drought, and to delineate the resulting ground-water levels and movement of the freshwater-saltwater interface. The model code selected was SHARP, a quasi-three-dimensional finite-difference method of simulating freshwater and saltwater flow simultaneously. Two sets of four proposed pumping scenarios were evaluated. The first represented average recharge from precipitation during 2006-20; the second represented the same period and conditions except for a 5-year period of drought conditions. The average-recharge simulations used the long-term (1959-99) rate of recharge; the drought simulations applied a 20-percent reduction in recharge rate and a 20-percent increase in the 1999 rate of agricultural pumpage during 2011–15. The simulated movement of the freshwater-saltwater interface in future withdrawal and recharge scenarios indicates that the interface may rise beneath pumped wells at Inlet Drive, Brecknock Hall, Main Bayview Road, Islands End, North Road, and Alvah's Lane. Either (1) movement of the interface to within 50 feet of the well screen, (2) a large percent change in the distance between the interface and the well screen, or (3) movement of the interface through a clay layer is a cause for concern. Wellfields in which saltwater intrusion does not appear to be a cause for concern were those at Ackerly Pond, Kenney's Road, Middle Road, Rocky Point Road, and hypothetical sites where future wellfields have been proposed.

Water-Resources Investigations Report↗
USGS1988

Hydrogeology and preliminary assessment of regional flow in the upper Cretaceous and adjacent aquifers in the northern Mississippi embayment

On a regional scale, the ground-water system of the northern Mississippi embayment is composed of a series of nonindurated elastic sediments that overlie a thick sequence of Paleozoic carbonates, sandstones, and shales. Precambrian crystalline rocks form both the structural and the hydrogeologic basement throughout the northern embayment. The units that comprise the hydrogeologic framework of this study are the alluvium-lower Wilcox aquifer, the Midway confining unit, the Upper Cretaceous aquifer, the Cretaceous-Paleozoic confining unit, and the Ozark-St. Francois aquifer. The Upper Cretaceous aquifer of Late Cretaceous age is the primary focus of this investigation; the study is part of the Gulf Coast Regional Aquifer-System Analysis. A ground-water flow model was developed as the main tool to refine the concepts of deep regional flow in the northern Mississippi embayment. This four layer finite-difference model enabled testing of alternative boundary concepts and provided a refined definition of the hydrologic budget of the deep aquifers. The alluvium-lower Wilcox aquifer, the Upper Cretaceous aquifer, and the Ozark-St.Francois aquifer form layers 2 through 4, respectively. Layer 1 is an inactive layer of constant heads representing shallow water levels, which are a major control on recharge to and discharge from the regional system. A matrix of leakance values simulates each confining unit, allowing vertical interchange of water between different aquifers. The model was calibrated to 1980 conditions by using the assumption that 1980 was near steady-state conditions; it was calibrated to simulate observed heads within acceptable limits. For this preliminary model, calculated heads were found to be most sensitive to pumping, and least sensitive to the leakance. By using all available water-quality and water-level data, alternative boundary conditions were tested by comparing model simulated heads to observed heads. Simulation indicated that the major discharge zone for the Upper Cretaceous aquifer occurred along a narrow area coincident with the boundary of a buried rift. The results of the early modeling effort also contribute to a better understanding of the regional hydrologic budget, indicating that upward leakage from the Ozark-St. Francois aquifer to the Upper Cretaceous aquifer is about 43 cubic feet per second, with about 30 cubic feet per second occuring west of the western margin of the embayment. Calculations suggest upward recharge of about 68 cubic feet per second occurs to the lower Wilcox-alluvium aquifer from the Upper Cretaceous aquifer. Simulation results also indicate that the Midway is an effective regional confining unit.

northern Mississippi embayment↗
USGS2004

Delineation of Areas Contributing Water to the Dry Brook Public-Supply Well, South Hadley, Massachusetts

Areas contributing water to the Dry Brook public-supply well in South Hadley, Massachusetts, were delineated with a numerical ground-water-flow model that is based on geologic and hydrologic information for the confined sand and gravel aquifer pumped by the supply well. The study area is along the Connecticut River in central Massachusetts, about 12 miles north of Springfield, Massachusetts. Geologic units in the study area consist of Mesozoic-aged sedimentary and igneous bedrock, late-Pleistocene glaciolacustrine sediments, and recent alluvial deposits of the Connecticut River flood plain. Dry Brook Hill, immediately south of the supply well, is a large subaqueous lacustrine fan and delta formed during the last glacial retreat by sediment deposition into glacial Lake Hitchcock from a meltwater tunnel that was likely near where the Connecticut River cuts through the Holyoke Range. The sediments that compose the aquifer grade from very coarse sand and gravel along the northern flank of the hill, to medium sands in the body of the hill, and to finer-grained sediments along the southern flank of the hill. The interbedded and overlapping fine-grained lacustrine sediments associated with Dry Brook Hill include varved silt and clay deposits. These fine-grained sediments form a confining bed above the coarse-grained aquifer at the supply well and partially extend under the Connecticut River adjacent to the supply well. Ground-water flow in the aquifer supplying water to Dry Brook well was simulated with the U.S. Geological Survey ground-water-flow modeling code MODFLOW. The Dry Brook aquifer model was calibrated to drawdown data collected from 8 observation wells during an aquifer test conducted by pumping the supply well for 10 days at a rate of 122.2 cubic feet per minute (ft3/min; 914 gallons per minute) and to water levels collected from observation wells across the study area. Generally, the largest hydraulic conductivity values used in the model were in the sand and gravel aquifer near the Dry Brook well, which is consistent with the geologic information. Results of aquifer-test simulation indicated that spatially variable aquifer hydraulic properties and boundary conditions affected heads and ground-water flow near the well. A comparison and analysis of water-level fluctuations in study area wells to fluctuations in the Connecticut River indicated a hydraulic connection of the aquifer with the river, which is also consistent with geologic information. Simulated ground-water levels indicated that most ground water in the study area flowed toward and discharged to the Connecticut River and the Dry Brook well. Small amounts of ground water also discharged to smaller streams (Dry Brook and Bachelor Brook) in the study area. Areas contributing water to the well were delineated with the MODPATH particle-tracking routine. Results of the contributing-area analysis indicated that the greatest sources of water to the well were recharge in the Dry Brook Hill area and infiltration of Connecticut River water in an area beyond the extent of the confining bed where the aquifer is in hydraulic connection with the river. The amount of water entering the Dry Brook well from recharge dominated at a lower pumping rate (40.0 ft3/min); about 90 percent of the pumped water originated from recharge and boundary flow, and infiltration from the Connecticut River supplied the remaining 10 percent. At a high pumping rate (122.2 ft3/min), however, about half of the water pumped from the Dry Brook well originated from recharge and boundary flow (49 percent), and half originated from infiltration of water from the Connecticut River (51 percent). Results of a sensitivity analysis of the extent of areas contributing water to the Dry Brook well when pumped at 122.2 ft3/min indicated that the size of these areas did not substantially change when aquifer properties were varied. In contrast, however, the size of these areas changed most when the recharge

Water-Resources Investigations Report↗
USGS2004

Simulated interaction between freshwater and saltwater and effects of ground-water pumping and sea-level change, lower Cape Cod aquifer system, Massachusetts

The U.S. Geological Survey, in cooperation with the National Park Service, Massachusetts Executive Office of Environmental Affairs, Cape Cod Commission, and the Towns of Eastham, Provincetown, Truro, and Wellfleet, began an investigation in 2000 to improve the understanding of the hydrogeology of the four freshwater lenses of the Lower Cape Cod aquifer system and to assess the effects of changing ground-water pumping, recharge conditions, and sea level on ground-water flow in Lower Cape Cod, Massachusetts. A numerical flow model was developed with the computer code SEAWAT to assist in the analysis of freshwater and saltwater flow. Model simulations were used to determine water budgets, flow directions, and the position and movement of the freshwater/saltwater interface. Model-calculated water budgets indicate that approximately 68 million gallons per day of freshwater recharge the Lower Cape Cod aquifer system with about 68 percent of this water moving through the aquifer and discharging directly to the coast, 31 percent flowing through the aquifer, discharging to streams, and then reaching the coast as surface-water discharge, and the remaining 1 percent discharging to public-supply wells. The distribution of streamflow varies greatly among flow lenses and streams; in addition, the subsurface geology greatly affects the position and movement of the underlying freshwater/saltwater interface. The depth to the freshwater/saltwater interface varies throughout the study area and is directly proportional to the height of the water table above sea level. Simulated increases in sea level appear to increase water levels and streamflows throughout the Lower Cape Cod aquifer system, and yet decrease the depth to the freshwater/saltwater interface. The resulting change in water levels and in the depth to the freshwater/saltwater interface from sea-level rise varies throughout the aquifer system and is controlled largely by non-tidal freshwater streams. Pumping from large-capacity municipal-supply wells increases the potential for effects on surface-water bodies, which are affected by pumping and wastewater-disposal locations and rates. Pumping wells that are upgradient of surface-water bodies potentially capture water that would otherwise discharge to these surface-water bodies, thereby reducing streamflow and pond levels. Kettle-hole ponds, such as Duck Pond in Wellfleet, that are near the top of a freshwater flow lens, appear to be more susceptible to changing pumping and recharge conditions than kettle-hole ponds closer to the coast or near discharge boundaries, such as the Herring River.

Scientific Investigations Report↗
USGS2003

Fate and Transport Modeling of Selected Chlorinated Organic Compounds at Hangar 1000, U.S. Naval Air Station, Jacksonville, Florida

The Jacksonville Naval Air Station occupies 3,800 acres adjacent to the St. Johns River in Jacksonville, Florida. Two underground storage tanks at Hangar 1000 contained solvents from the late 1960s until they were removed in 1994. Ground-water samples at one of the tank sites had levels of trichloroethene (TCE) and total dichloroethene (DCE) of 8,710 micrograms per liter (mg/L) and 4,280 mg/L, respectively. Vinyl chloride (VC) at the site is the result of the biodegradation of DCE. Ground water beneath Hangar 1000 flows toward a storm sewer. TCE and DCE plumes travel with the ground water and presumably have reached the storm sewer, which discharges to the St. Johns River. Simulation of solute transport indicates that the traveltime from the storage tank site to the storm sewer is 16, 14, and 12 years for TCE, DCE, and VC respectively. TCE has the longest traveltime because it has the highest retardation factor at 2.5, DCE takes less time with a retardation factor of 2.0, and VC has the quickest traveltime because it has the lowest retardation factor of 1.7. Based on modeling results, the release of contaminants in the aquifer occurred more than 16 years ago. Model-derived dispersivity values at Hangar 1000 were: longitudinal 1.5 feet (ft), transverse 0.27 ft, and vertical 0.27 ft. The model-derived first order decay rates for biodegradation of TCE, DCE, and VC were 0.0002 per day (d-1), 0.0002 d-1, and 0.06 d-1, respectively. These rates are equivalent to half-lives of 13.7 years for TCE and DCE and 17 days for VC. Source area reductions in contaminant concentrations of 50 and 100 percent were modeled to simulate remediation. As expected, reducing the source concentration by 50 percent resulted in eventual TCE, DCE, and VC concentrations that were half of the original concentrations. About 16 years were needed for new steady-state TCE concentrations to develop, about 14 years for DCE, and about 12 years for VC. Reducing the source area concentrations by 100 percent in the model eventually resulted in zero concentrations of TCE, DCE, and VC. The modeled period of time for the contaminants to be removed from the aquifer once the source was removed was about 17 years for TCE, 15 years for DCE, and 13 years for VC.

Water-Resources Investigations Report↗
USGS1998

Natural attenuation of chlorinated-hydrocarbon contamination at Fort Wainwright, Alaska: A hydrogeochemical and microbiological investigation workplan

Natural attenuation processes include biological degradation, by which microorganisms break down contaminants into simpler product compounds; adsorption of contaminants to soil particles, which decreases the mass of contaminants dissolved in ground water; and dispersion, which decreases dissolved contaminant concentrations through dilution. The primary objectives of this study are to (1) assess the degree to which such natural processes are attenuating chlorinated-hydrocarbon contamination in ground water, and (2) evaluate the effects of ground-water/surface-water interactions on natural-attenuation processes in the area of the former East and West Quartermasters Fueling Systems for Fort Wainwright, Alaska. The study will include investigations of the hydrologic, geochemical, and microbiological processes occurring at this site that influence the transport and fate of chlorinated hydrocarbons in ground water. To accomplish these objectives, a data-collection program has been initiated that includes measurements of water-table elevations and the stage of the Chena River; measurements of vertical temperature profiles within the subsurface; characterization of moisture distribution and movement in the unsaturated zone; collection of ground-water samples for determination of both organic and inorganic chemical constituents; and collection of ground-water samples for enumeration of microorganisms and determination of their potential to mineralize contaminants. We will use results from the data-collection program described above to refine our conceptual model of hydrology and contaminant attenuation at this site. Measurements of water-table elevations and river stage will help us to understand the magnitude and direction of ground-water flow and how changes in the stage of the Chena River affect ground-water flow. Because ambient ground water and surface water typically have different temperature characteristics, temperature monitoring will likely provide further insight into ground-water/surface-water interactions in the subsurface. Characterization of the unsaturated zone will improve our understanding of interactions among ground water, the unsaturated zone, and the atmosphere. The interactions likely of importance to this study include the migration of water, dissolved contaminants, nutrients, and gases (oxygen, carbon dioxide, and methane) between the saturated and unsaturated zones. We will use the results of ground-water chemical analyses to determine the spatial and temporal distribution of (1) chlorinated-hydrocarbon contaminants and their degradation products, (2) oxidation-reduction indicators, (3) nutrients, and (4) major ground-water ions. These water-quality data will provide insight into ground-water flow directions, interactions between ground water and surface water, attenuation of contaminant concentrations caused by dispersion, and intrinsic microbiological processes. Microbiological analyses will indicate whether microorganisms at the site are capable of degrading the contaminants of interest, and will allow us to estimate their potential to attenuate existing contamination. Physical and chemical data interpreted as part of the analysis of ground water and surface water mixing will improve our understanding of the relationship between water quality and contaminant source mixing.

Alaska↗
USGS1986

Ground-water resource assessment of the Montauk area, Long Island, New York

The water resources of the Montauk area were investigated from October 1980 through September 1983 to assess the availability of fresh groundwater. The principal aquifer, which consists of fine- to coarse-grained stratified glacial drift, is the sole source of freshwater. The freshwater/saltwater interface lies as much as 150 ft below sea level, but the till unit and the marine-clay unit limit the thickness of the zone from which freshwater may easily be withdrawn to < 100 ft in most places. Precipitation, the sole source of freshwater in the Montauk area , averages about 42 in/yr. Direct runoff is approximately 1 in/yr, and evapotranspiration about 20 in/yr; the remaining 21 inches is the net recharge to the principal aquifer. The horizontal hydraulic conductivity, as calculated from specific-capacity data, ranges from 130 to 350 ft/d. Analysis of an aquifer pumping test indicates that horizontal hydraulic conductivity of the principal aquifer is approximately 275 ft/d and the vertical hydraulic conductivity is approximately 90 ft/d; the storage coefficient ranges from 1.1 x .001 to 2.2 x .001. A numerical two-dimensional flow model that simulates the areal flow of fresh- and saltwater was developed to evaluate the effects of the present level of groundwater development and the potential effects of increased future withdrawals. Model results indicate that the principal aquifer is capable of producing several times more than the 0.6 mil gal/d now being withdrawn for public supply. Examination of well design and pumping rates, and analysis of analog-model studies of saltwater upconing to wells, indicate that the peak pumping rate of public supply wells exceeds the limit necessary to avoid contamination from saltwater. With proper well design and distribution of pumping stress, the aquifer system would be capable of supplying enough water to meet the needs of the population projected through 1995.

New York↗
USGS1996

Estimation of the recharge areas contributing water to the South Well Field, Columbus, Ohio

The city of Columbus, Ohio, operates four radial collector wells, designed to yield 42 Mgal/d (million gallons per day), in southern Franklin County, Ohio, as part of their municipal supply of water. The collector wells are adjacent to, and designed to induce infiltration from, Big Walnut Creek and Scioto River. A previously constructed, three-dimensional, steady-state and transient ground-water-flow model of this river-aquifer system was used to estimate contributing recharge areas (CRA's) and calculate particle flowpaths in southern Franklin County. The simulations were of two steady-state periods (October 1979 and March 1986) and one 5-year transient period (March 1986---June 1991). The first simulation (1979) was of conditions before construction of the collector wells. The second simulation (1986) was of conditions when the collector wells were producing 8 Mgal/d. During the 5 years covered in the transient simulation, production at the well field averaged 18.5 Mgal/d. Under the 1979 conditions, the largest ground-water contributing areas were of the quarries and Scioto River (41 and 47 percent of the study area, respectively). During 1986, when 8 Mgal/d was withdrawn, the primary contributing areas were of the quarries (40 percent), collector wells (34 percent), and rivers (8 percent). Travel times associated with simulated particles of water tracked from cells along Big Walnut Creek to their discharge points in cells along Scioto River were about 5 to 60 years in the 1979 simulation and about 7 to 41 years in the 1986 simulation. The endpoints of these particles varied as simulated pumping rates were increased to 22 Mgal/d. The 1986, 10-year CRA's of the collector wells under 8 Mgal/d-conditions totalled about 4.5 mi2. As the pumping rate was increased to 22 Mgal/d in a predictive simulation, 10-year CRA's of the collector wells increased to 6.7mi2. Because the transient simulation encompassed only 5 years, the 10-year CRA's could not be estimated from the transient simulation. However, the size of the 1- to 5-year CRA's for the transient simulation was similar to the size of the 1- to 5-year CRA's for a steady-state predictive simulation if well-field production were 16 Mgal/d. The transient simulations predicted discontinuous CRA's, especially adjacent to the rivers, due to changes in hydrologic stresses. Analyses of the steady-state and transient models showed that sizes of CRA's were most sensitive to changes in porosity, pumping rate, riverbed conductance, and horizontal hydraulic conductivity.

Ohio↗
USGS1986

Ground-water monitoring at Santa Barbara, California: Phase 3 — Development of a three-dimensional digital ground-water flow model for storage unit I of the Santa Barbara ground-water basin

Water-bearing rocks within the 7 sq mi of Storage Unit I of the Santa Barbara Groundwater Basin, consist of unconsolidated deposits that range in thickness from < 300 ft along the north perimeter of the unit to > 1,000 ft near the Pacific Ocean. The groundwater system was simulated as two horizontal layers separated by a confining bed. The model boundaries coincide with mapped faults on all sides. The faults were considered no-flow boundaries except for the offshore fault that forms the south boundary. This boundary was simulated as a general-head boundary , which allows water to move into and out of the modeled area. The model was calibrated by simulating both steady-state conditions (approximated by July 1978 and February 1983 water levels) and transient-state conditions (represented by May 1978 through December 1979 water level changes). The calibrated model was then used to simulate the period from January 1980 through December 1983 in order to verify the model. Model results generally closely matched measured data throughout Storage Unit I. During the transient and verification simulations, 9,980 acre-ft of groundwater was pumped from Storage Unit I for municipal use. Results of the model indicate that 42% (4,190 acre-ft) of the water pumped from the system was withdrawn from storage, 33% (3,290 acre-ft) was derived from changes in underflow across the offshore fault, and 25% (2,500 acre-ft) was derived from decreased groundwater discharge to drains. The model simulated that municipal pumpage induced about 1,380 acre-ft of water to move across the offshore fault toward Storage Unit I. Several model simulations were used to estimate aquifer response to different municipal pumpage patterns that could be used as management alternatives. Results of the simulations indicate that spreading municipal pumpage more evenly throughout Storage Unit I, by increasing the number of wells while reducing the pumping rate at the individual wells to maintain the same total pumpage, significantly reduces the inflow of groundwater across the offshore fault. (Author 's abstract)

California↗
USGS1996

Estimation of evapotranspiration in the Rainbow Springs and Silver Springs basins in North-Central Florida

Estimates of evapotranspiration (ET) for the Rainbow and Silver Springs ground-water basins in north-central Florida were determined using a regional water-~budget approach and compared to estimates computed using a modified Priestley-Taylor (PT) model calibrated with eddy-correlation data. Eddy-correlation measurements of latent 0~E) and sensible (H) heat flux were made monthly for a few days at a time, and the PT model was used to estimate 3,E between times of measurement during the 1994 water year. A water-budget analysis for the two-basin area indicated that over a 30-year period (196594) annual rainfall was 51.7 inches. Of the annual rainfall, ET accounted for about 37.9 inches; springflow accounted for 13.1 inches; and the remaining 0.7 inch was accounted for by stream-flow, by ground-water withdrawals from the Floridan aquifer system, and by net change in storage. For the same 30-year period, the annual estimate of ET for the Silver Springs basin was 37.6 inches and was 38.5 inches for the Rainbow Springs basin. Wet- and dry-season estimates of ET for each basin averaged between nearly 19 inches and 20 inches, indicating that like rainfall, ET rates during the 4-month wet season were about twice the ET rates during the 8-month dry season. Wet-season estimates of ET for the Rainbow Springs and Silver Springs basins decreased 2.7 inches, and 3.4 inches, respectively, over the 30-year period; whereas, dry-season estimates for the basins decreased about 0.4 inch and1.0 inch, respectively, over the 30-year period. This decrease probably is related to the general decrease in annual rainfall and reduction in net radiation over the basins during the 30-year period. ET rates computed using the modified PT model were compared to rates computed from the water budget for the 1994 water year. Annual ET, computed using the PT model, was 32.0 inches, nearly equal to the ET water-budget estimate of 31.7 inches computed for the Rainbow Springs and Silver Springs basins. Modeled ET rates for 1994 ranged from 14.4 inches per year in January to 51.6 inches per year in May. Water-budget ET rates for 1994 ranged from 12.0 inches per year in March to 61.2 inches per year in July. Potential evapotranspiration rates for 1994 averaged 46.8 inches per year and ranged from 21.6 inches per year in January to 74.4 inches per year in May. Lake evaporation rates averaged 47.1 inches per year and ranged from 18.0 inches per year in January to 72.0 inches per year in May 1994.

Water-Resources Investigations Report↗
USGS2007

Simulated Ground-Water Withdrawals by Cabot WaterWorks from the Mississippi River Valley Alluvial Aquifer, Lonoke County, Arkansas

Cabot WaterWorks, located in Lonoke County, Arkansas, plans to increase ground-water withdrawals from the Mississippi River Valley alluvial aquifer from a 2004 rate of approximately 2.24 million gallons per day to between 4.8 and 8 million gallons per day by the end of 2049. The effects of increased pumping from several wells were simulated using a digital model of ground-water flow. The proposed additional withdrawals by Cabot WaterWorks were specified in three 1-square-mile model cells with increased pumping beginning in 2007. Increased pumping was specified at various combined rates for a period of 44 years. In addition, augmented pumping from wells owned by Grand Prairie Water Users Association, located about 2 miles from the nearest Cabot WaterWorks wells, was added to the model beginning in 2007 and continuing through to the end of 2049 in 10 of the 16 scenarios analyzed. Eight of the scenarios included reductions in pumping rates in model cells corresponding to either the Grand Prairie Water Users Association wells or to wells contained within the Grand Prairie Area Demonstration Project. Drawdown at the end of 44 years of pumping at 4.8 million gallons per day from the Cabot WaterWorks wells ranged from 15 to 25 feet in the three model cells; pumping at 8 million gallons per day resulted in water-level drawdown ranging from about 15 to 40 feet. Water levels in those cells showed no indication of leveling out at the end of the simulation period, indicating non-steady-state conditions after 44 years of pumping. From one to four new dry cells occurred in each of the scenarios by the end of 2049 when compared to a baseline scenario in which pumping was maintained at 2004 rates, even in scenarios with reduced pumping in the Grand Prairie Area Demonstration Project; however, reduced pumping produced cells that were no longer dry when compared to the baseline scenario at the end of 2049. Saturated thickness at the end of 2049 in the three Cabot WaterWorks wells ranged from about 52 to 68.5 feet for pumping rates of 4.8 million gallons per day, and from about 38 to 64 feet for pumping rates of 8 million gallons per day, the latter causing water level to fall below half the aquifer thickness in the most heavily pumped of the three cells.

Scientific Investigations Report↗
USGS1993

Hydrogeology, simulated ground-water flow, and ground-water quality, Wright-Patterson Air Force Base, Ohio

Ground water is the primary source of water in the Wright-Patterson Air Force Base area. The aquifer consists of glacial sands and gravels that fill a buried bedrock-valley system. Consolidated rocks in the area consist of poorly permeable Ordovician shale of the Richmondian stage, in the upland areas, the Brassfield Limestone of Silurian age. The valleys are filled with glacial sediments of Wisconsinan age consisting of clay-rich tills and coarse-grained outwash deposits. Estimates of hydraulic conductivity of the shales based on results of displacement/recovery tests range from 0.0016 to 12 feet per day; estimates for the glacial sediments range from less than 1 foot per day to more than 1,000 feet per day. Ground water flow from the uplands towards the valleys and the major rivers in the region, the Great Miami and the Mad Rivers. Hydraulic-head data indicate that ground water flows between the bedrock and unconsolidated deposits. Data from a gain/loss study of the Mad River System and hydrographs from nearby wells reveal that the reach of the river next to Wright-Patterson Air Force Base is a ground-water discharge area. A steady-state, three-dimensional ground-water-flow model was developed to simulate ground-water flow in the region. The model contains three layers and encompasses about 100 square miles centered on Wright-Patterson Air Force Base. Ground water enters the modeled area primarily by river leakage and underflow at the model boundary. Ground water exits the modeled area primarily by flow through the valleys at the model boundaries and through production wells. A model sensitivity analysis involving systematic changes in values of hydrologic parameters in the model indicates that the model is most sensitive to decreases in riverbed conductance and vertical conductance between the upper two layers. The analysis also indicates that the contribution of water to the buried-valley aquifer from the bedrock that forms the valley walls is about 2 to 4 percent of the total ground-water flow in the study area. Ground waters in the vicinity of Wright-Patterson Air Force Base can be classified into two compositional groups on the basis of their chemical composition: calcium magnesium bicarbonate-type and sodium chloride-type waters. Calcium magnesium bicarbonate-type waters are found in the glacial deposits and the Brassfield Limestone, whereas the sodium chloride waters are exclusively associated with the shales. Equilibrium speciation calculations indicate that ground water of the glacial drift aquifer is in equilibrium with calcite, dolomite, and chalcedony, but is undersaturated with respect to gypsum and fluorite. Waters from the shales are slightly supersaturated with respect to calcite, dolomite, and siderite but are undersaturated with respect to chalcedony. Simple-mass balance calculations treating boron as a conservative species indicate that little (< 5 percent) or no recharge from the shales to the glacial drift aquifer takes place. Data on the stable isotopes of oxygen and hydrogen indicate a meteoric origin for all ground water beneath Wright-Patterson Air Force Base, but the data were inconclusive with respect to identification of distinct isotopic differences between water collected from the glacial drift and bedrock aquifers. Tritium concentrations used to distinguish waters having a pre-and post-1953 recharge component indicate that most water entered the glacial drift aquifer after 1953. This finding indicates that recharge from shallow to deep parts (greater than 150 feet) of the aquifer takes place over time intervals of a few years or decades. However, the fact that some deep parts of the glacial aquifer did not contain measurable tritium indicates that ground-water flow from recharge zones to these parts of the aquifer takes decades or longer.

Water-Resources Investigations Report↗
USGS2008

Use of Numerical Models to Simulate Transport of Sewage-Derived Nitrate in a Coastal Aquifer, Central and Western Cape Cod, Massachusetts

The unconsolidated glacial sediments underlying Cape Cod, Massachusetts compose a regional aquifer system that is used both as a source of drinking water and as a disposal site for wastewater; in addition, the discharge of clean ground water from the aquifer system is needed for the maintenance of freshwater and marine ecosystems throughout the region. Because these uses of the aquifer conflict with one another in many areas of the Cape, local and regional planners have begun to develop sustainable wastewater plans that will facilitate the disposal of wastewater while protecting water supplies and improving the health of aquatic ecosystems. To assist local and regional planners in these efforts, the U.S. Geological Survey conducted a 2-year investigation to (1) assist local and regional planners in the evaluation of potential wastewater scenarios, (2) use results and interpretation from these analyses to develop hydrologic concepts transferable throughout the region, and (3) establish and test methods that would be of use in future evaluations. Wastewater-disposal scenarios need to be evaluated in the context of the regional ground-water-flow system. For a given rate of disposal, wastewater from sites at or near a regional ground-water divide is transported in a wider arc of flow directions, flows deeper in the system, and contaminates a larger part of the aquifer than does wastewater discharged from sites farther from the divide. Also, traveltimes of wastewater from sites near a ground-water divide to receptors are longer (as much as several hundred years) than traveltimes from sites farther from the divide. Thus, wastewater disposal at or near a divide will affect a larger part of the aquifer and likely contribute wastewater to more receptors than wastewater disposal farther from a divide; however, longer traveltimes could allow for more attenuation of wastewater-derived nitrate from those sites. Ground-water-flow models and particle tracking can be used to identify advective-transport patterns downgradient from wastewater-disposal sites and estimate traveltimes; however, these tools cannot predict the distribution of mass or concentrations of wastewater constituents, such as nitrate, in the aquifer. Flow-based particle-tracking analyses can be used to estimate mass-loading rates and time-varying concentrations at wells and ecological receptors by the accounting of mass-weighted particles discharging into the receptor of interest. This method requires no additional development beyond the flow model; however, post-modeling analyses are required. In addition, the method is based on the assumption that no mass is lost during transport, an assumption that likely is not valid in many systems. Solute-transport models simulate the subsurface transport of nitrate through the aquifer and predict the distribution of the mass of a solute in the aquifer at different transport times. This method does require additional model development beyond the flow model, but can predict timevarying concentrations at receptors. Estimates of mass-loading rates require minimal post-modeling analyses. Time-varying concentrations and mass-loading rates calculated for wells in eastern Barnstable by the two methods generally were in reasonable agreement. Inherent in the flow-based particle-tracking method is the assumption that mass is conserved along a given flow line and that there is no spreading of mass in the aquifer. Although the solute-transport models also incorporate a system-wide conservation of mass, these models allow for a spreading of mass in the aquifer, and mass is not conserved along a given flow line. As a result, estimates of concentrations and mass loading rates generally were higher in particle-tracking analyses than in solute-transport simulations. Results from the two types of simulations agreed best for wells that receive large amounts of wastewater with short traveltimes (less than 10 years) because insufficient transport

Scientific Investigations Report↗
USGS2007

Hydrogeology and simulation of ground-water flow near Mount Pleasant, South Carolina: Predevelopment, 2004, and predicted scenarios for 2030

Heavy water use from the Cretaceous Middendorf aquifer in South Carolina has created a large, regional cone of depression in the potentiometric surface of the Middendorf aquifer in Charleston and Berkeley Counties, South Carolina. Water-level declines of up to 249 feet have been observed in wells over the past 125 years and are a result of ground-water use for public-water supply, irrigation, and private industry. To address the concerns of users of the Middendorf aquifer, the U.S. Geological Survey, in cooperation with Mount Pleasant Waterworks, updated an existing ground-water flow model to incorporate additional data that have been compiled since 1989. The updated ground-water flow model incorporates water-level data collected from 349 wells in 2004, baseflow data measured at 17 streams, hydraulic property data from 265 wells, and water-use data compiled for more than 2,700 wells for the period between the early 1900s to 2004. The ground-water flow system of the Coastal Plain physiographic province of South Carolina and parts of Georgia and North Carolina was simulated using the U.S. Geological Survey finite-difference code MODFLOW-2000. The model was vertically discretized into nine layers to include the five aquifers of the surficial, the combined Floridan aquifer system and Tertiary sand aquifer, Black Creek, Middendorf, and Cape Fear, separated by four intervening confining units. Specified-head boundary conditions were used at the lateral boundaries of the model and for the lower Coastal Plain part of the surficial aquifer; no-flow boundary conditions were used at the updip and downdip extent of the model layers and at the base of the Cape Fear aquifer. Ground-water conditions for predevelopment and 2004 were simulated using steady-state and transient approximations, respectively. Simulated water levels generally matched the observed conditions, plus or minus a 20-foot calibration target, with 56.4 and 64.8 percent of the simulated values approximating the measured values for predevelopment and 2004 hydrologic conditions, respectively. The root-mean-square error of the water-level residuals for the various model layers varied between 20.2 and 34.4 feet for predevelopment and 18.2 and 36.7 feet for 2004. The general goodness of fit also was apparent in the calculation of the ratio of standard deviation of residuals to range of observations for each modeled aquifer layer. The calculated ratios for the predevelopment and 2004 hydrologic conditions were less than 0.10 for all model layers except for the Cape Fear aquifer in both predevelopment and 2004 simulations. The Mount Pleasant model was most sensitive to changes in simulated specific storage of most model layers, vertical anisotropy of the confining units above and below the Middendorf aquifer, hydraulic conductivity of the confining units, and the specified-head boundary conditions for the surficial aquifer. The model also is sensitive to horizontal hydraulic conductivity of the Floridan aquifer system and Tertiary sand aquifer and the Black Creek and Middendorf aquifers. Simulated water budgets indicate that the primary sources of water to the model are recharge and the specified-head boundaries in layers 1 and 3. More than 88 percent of the water that discharges from the model discharges from layers 1-3 through specified-head boundaries and rivers. Approximately 11 percent of the water budget was discharged through wells for the 2004 budget. In 2004, 8.11 million gallons of water per day was discharged from wells in the Mount Pleasant area. Water to these wells is provided predominantly by lateral flow within the Middendorf aquifer. Additional water is provided from aquifer storage and leakage from confining units located above and below the Middendorf aquifer. Downward flow through the Middendorf confining unit is a reversal of the predevelopment flow direction. Five predictive water-management scenarios were simulated to determine the effects on the

South Carolina↗