Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Ground Water Modeling Investigation Report”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 901 records · Page 50Linked to original sources

Hydrologic and Water-Quality Responses in Shallow Ground Water Receiving Stormwater Runoff and Potential Transport of Contaminants to Lake Tahoe, California and Nevada, 2005-07

Clarity of Lake Tahoe, California and Nevada has been decreasing due to inflows of sediment and nutrients associated with stormwater runoff. Detention basins are considered effective best management practices for mitigation of suspended sediment and nutrients associated with runoff, but effects of infiltrated stormwater on shallow ground water are not known. This report documents 2005-07 hydrogeologic conditions in a shallow aquifer and associated interactions between a stormwater-control system with nearby Lake Tahoe. Selected chemical qualities of stormwater, bottom sediment from a stormwater detention basin, ground water, and nearshore lake and interstitial water are characterized and coupled with results of a three-dimensional, finite-difference, mathematical model to evaluate responses of ground-water flow to stormwater-runoff accumulation in the stormwater-control system. The results of the ground-water flow model indicate mean ground-water discharge of 256 acre feet per year, contributing 27 pounds of phosphorus and 765 pounds of nitrogen to Lake Tahoe within the modeled area. Only 0.24 percent of this volume and nutrient load is attributed to stormwater infiltration from the detention basin. Settling of suspended nutrients and sediment, biological assimilation of dissolved nutrients, and sorption and detention of chemicals of potential concern in bottom sediment are the primary stormwater treatments achieved by the detention basins. Mean concentrations of unfiltered nitrogen and phosphorus in inflow stormwater samples compared to outflow samples show that 55 percent of nitrogen and 47 percent of phosphorus are trapped by the detention basin. Organic carbon, cadmium, copper, lead, mercury, nickel, phosphorus, and zinc in the uppermost 0.2 foot of bottom sediment from the detention basin were all at least twice as concentrated compared to sediment collected from 1.5 feet deeper. Similarly, concentrations of 28 polycyclic aromatic hydrocarbon compounds were all less than laboratory reporting limits in the deeper sediment sample, but 15 compounds were detected in the uppermost 0.2 foot of sediment. Published concentrations determined to affect benthic aquatic life also were exceeded for copper, zinc, benz[a]anthracene, phenanthrene, and pyrene in the shallow sediment sample. Isotopic composition of water (oxygen 18/16 and hydrogen 2/1 ratios) for samples of shallow ground water, lakewater, and interstitial water from Lake Tahoe indicate the lake was well mixed with a slight ground-water signature in samples collected near the lakebed. One interstitial sample from 0.8 foot beneath the lakebed was nearly all ground water and concentrations of nitrogen and phosphorus were comparable to concentrations in shallow ground-water samples. However, ammonium represented 65 percent of filtered nitrogen in this interstitial sample, but only 10 percent of the average nitrogen in ground-water samples. Nitrate was less than reporting limits in interstitial water, compared with mean nitrate concentration of 750 micrograms per liter in ground-water samples, indicating either active dissimilative nitrate reduction to ammonium by micro-organisms or hydrolysis of organic nitrogen to ammonium with concomitant nitrate reduction. The other interstitial sample falls along a mixing line between ground water and lake water and most of the nitrogen was organic nitrogen.

Scientific Investigations Report↗

Distribution and sources of nitrate, and presence of fluoride and pesticides, in parts of the Pasco Basin, Washington, 1986-88

Ground water was sampled in a 900-square-mile agricultural area in the Pasco Basin, which includes parts of eastern Benton County and western Franklin County, Washington, to determine distributions of nitrate and fluoride. Additional data were obtained to determine if fertilizers, irrigation water, septic systems, and naturally occurring nitrate are sources of nitrate in ground water. Limited sampling also was done to determine if pesticides were present in the ground water. Nitrate concentrations in ground water ranged from less than 0.1 to 100 milligrams per liter as nitrogen, and median concentrations of nitrate nitrogen in ground water were 3.2 and 6.7 milligrams per liter for Benton and Franklin Counties, respectively. In Franklin County, where a large percentage of the land is used for irrigated agriculture, nitrate nitrogen concentrations in water from 31 percent of sampled wells were equal to or greater than the maximum contaminant level for drinking water of 10 milligrams per liter. In Benton County, nitrate concentrations in water from about 10 percent of the sampled wells exceeded the maximum contaminant level. Nitrate concentrations in ground water at some locations in Franklin County have increased by as much as two orders of magnitude since the early 1950's. Historical data generally were not available to evaluate changes of nitrate concentrations in ground water in Benton County, except for the area around the town of Finley. A comparison of data collected during this study with data collected during 1976-77 indicate that nitrate concentrations in ground water of the Finley area probably have not changed over the intervening period. Applied nitrogen fertilizers are a major source of nitrate in ground water at many locations in the study area. Surface water used for irrigation does not contain sufficient nitrate to cause elevated concentrations in ground water. Instead, canal seepage, which makes up about 50 percent of the ground-water recharge in the study area, tends to dilute the nitrate present in ground water. Septic systems in the Finley area of Benton County are a source of nitrate in ground water, but analyses of data and results of a numerical model analysis of nitrate concentrations in the unconfmed ground-water system indicate that they are not the primary source of nitrate in ground water in this area. Naturally occurring nitrate may be a source of nitrate in ground water underlying Badger Coulee in Benton County. Average masses of natural nitrate per unit volume of sediment in two boreholes in Badger Coulee were equivalent to 2,590 and 964 pounds of nitrogen, respectively, in a block of sediments 50 feet thick underlying an acre of land. At most other locations in the study area, the amount of natural nitrate in ground water is probably small compared with nitrate from anthropogenic sources. Fluoride concentrations in ground water in the study area ranged from less than 0.1 to 4.7 milligrams per liter; the median concentration was 0.5 milligram per liter. The concentration of fluoride in water from only two of 143 wells equalled or exceeded 2.0 milligrams per liter, which is the secondary maximum contaminant level for drinking water. Both are deep wells open to the Saddle Mountains Basalt in Franklin County. Large concentrations of fluoride in deep ground waters of the Pasco Basin are apparently the result of natural conditions in the deeper basalt aquifers. One or more pesticide compounds were detected in 10 of 29 ground-water samples, which were analyzed for selected chlorophenoxy acid herbicides, triazine herbicides, carbamate insecticides, organophosphorus insecticides, and a few other types of pesticides. The sampling locations did not represent a random distribution, but instead, most were wells open to unconfined, shallow ground water in irrigated areas. The pesticides found include the herbicides atrazine, dicamba, metribuzin, picloram, and 2,4,5-T. Also present were aldicarb sulfone and aldicarb sulfoxide, which are degradation products of the insecticide aldicarb. Except for metribuzin, pesticide concentrations were at or near the analytical reporting limits. In all instances, the concentrations of pesticides detected were below the health advisory levels that are issued by the U.S. Environmental Protection Agency

Washington↗

Estimates of ground-water recharge to the Yakima River Basin aquifer system, Washington, for predevelopment and current land-use and land-cover conditions

Two models were used to estimate ground-water recharge to the Yakima River Basin aquifer system, Washington for predevelopment (estimate of natural conditions) and current (a multi-year, 1995-2004, composite) land-use and land-cover conditions. The models were the Precipitation-Runoff Modeling System (PRMS) and the Deep Percolation Model (DPM) that are contained in the U.S. Geological Survey's Modular Modeling System. Daily values of recharge were estimated for water years 1950-98 using previously developed PRMS-watershed models for four mainly forested upland areas, and for water years 1950-2003 using DPM applied to 17 semiarid to arid areas in the basin. The mean annual recharge under predevelopment conditions was estimated to be about 11.9 in. or 5,450 ft 3 /s (about 3.9 million acre-ft) for the 6,207 mi 2 in the modeled area. In the modeled areas, recharge ranged from 0.08 in. (1.2 ft 3 /s) to 34 in. (2,825 ft 3 /s). About 97 percent of the recharge occurred in the 3,667 mi 2 area included in the upland-area models, but much of this quantity is not available to recharge the bedrock hydrogeologic units. Only about 1.0 in., or 187 ft 3 /s (about 0.14 million acre-ft), was estimated to occur in the 2,540 mi 2 area included in the semiarid to arid lowland modeled areas. The mean annual recharge to the aquifer system under current conditions was estimated to be about 15.6 in., or 7,149 ft 3 /s (about 5.2 million acre-ft). The increase in recharge is due to the application of irrigation water to croplands. The annual quantity of irrigation was more than five times the annual precipitation for some of the modeled areas. Mean annual actual evapotranspiration was estimated to have increased from predevelopment conditions by more than 1,700 ft 3 /s (about 1.2 million acre-ft) due to irrigation.

Washington↗

Effects of pumping municipal wells at Manhattan, Kansas, on streamflow in the Big Blue and Kansas Rivers, Northeast Kansas, 1992-94

A ground-water flow model was developed to simulate the effects of municipal well pumping on streamflow in the Big Blue and Kansas Rivers near Manhattan, Kansas, from 1992 through 1994.Model simulations of the effects of municipal well pumping on streamflow in the Big Blue and Kansas Rivers indicate that well pumping decreases streamflow. Simulations of May 1993 conditions indicate that well pumping decreased simulated streamflow in the Big Blue and Kansas Rivers by 5.28 ft3/s (cubic feet per second) for the month, of which 3.22 ft3/s were contributed from the streams (induced infiltration) and 2.06 ft3/s werecontributed from ground water that would have seeped to the streams if the wells had not been pumping (intercepted base flow). Of the total 414acre-feet pumped by municipal wells during May 1993, about 48 percent was from induced infiltration, and about 31 percent was from intercepted base flow. Simulations of October 16 through November 14, 1994, conditions indicate that well pumping decreased simulated streamflow in the Big Blue and Kansas Rivers by 6.67 ft3/s for the period, of which 6.51 ft3/s was from induced infiltration and 0.61 ft3/s was from intercepted base flow. Of the 506 acre-feet pumped by municipal wells during October 16 through November 14, 1994, about 76 percent was induced from infiltration, and about 2 percent was from intercepted base flow. Steady-state simulations of hypothetical conditions were conducted to develop relations among average and minimum ground- water altitudes in the Manhattan municipal well field and precipitation, pumping, and streamflow rates.

Water-Resources Investigations Report↗

Geophysical investigations at Hidden Dam, Raymond, California — Flow simulations

Numerical flow modeling and analysis of observation-well data at Hidden Dam are carried out to supplement recent geophysical field investigations at the site (Minsley and others, 2010). This work also is complementary to earlier seepage-related studies at Hidden Dam documented by Cedergren (1980a, b). Known seepage areas on the northwest right abutment area of the downstream side of the dam was documented by Cedergren (1980a, b). Subsequent to the 1980 seepage study, a drainage blanket with a sub-drain system was installed to mitigate downstream seepage. Flow net analysis provided by Cedergren (1980a, b) suggests that the primary seepage mechanism involves flow through the dam foundation due to normal reservoir pool elevations, which results in upflow that intersects the ground surface in several areas on the downstream side of the dam. In addition to the reservoir pool elevations and downstream surface topography, flow is also controlled by the existing foundation geology as well as the presence or absence of a horizontal drain in the downstream portion of the dam. The current modeling study is aimed at quantifying how variability in dam and foundation hydrologic properties influences seepage as a function of reservoir stage. Flow modeling is implemented using the COMSOL Multiphysics software package, which solves the partially saturated flow equations in a two-dimensional (2D) cross-section of Hidden Dam that also incorporates true downstream topography. Use of the COMSOL software package provides a more quantitative approach than the flow net analysis by Cedergren (1980a, b), and allows for rapid evaluation of the influence of various parameters such as reservoir level, dam structure and geometry, and hydrogeologic properties of the dam and foundation materials. Historical observation-well data are used to help validate the flow simulations by comparing observed and predicted water levels for a range of reservoir elevations. The flow models are guided by, and discussed in the context of, the geophysical work (Minsley and others, 2010) where appropriate.

California↗

Ground-water flow in the Gulf Coast aquifer systems, south central United States — A preliminary analysis

A major objective of the Gulf Coast Regional Aquifer-System Analysis is to use digital models of regional groundwater flow systems to develop better understanding and to improve management of the resource. Modeling is used to synthesize information about the aquifer systems and to test hypotheses about the relative importance of the components of the systems. The 290,000-sq mile study area in the Gulf of Mexico Coastal Plain includes the Mississippi embayment, Gulf Coastal Plain of Texas, and the Continental Shelf that are underlain by deposits of Tertiary and younger age, which contain fresh and saline water. A 10-layer, finite-difference, variable density model, with blocks 10 miles on a side, was used to simulate groundwater flow before development and in 1980, assuming steady- state conditions. Preliminary results indicate that the major factors controlling predevelopment regional flow are the topography, land-surface outcrop pattern, and geometry of aquifers and confining units. Geologic structure and the distribution of precipitation were less significant factors. The density of saline water in the deeper parts of the aquifer system probably has a substantial effect on regional groundwater flow that extends into the freshwater part of the system. Variable water density may be a significant driving force that transports salt great distances in many directions, including updip. The distribution and rates of regional recharge and discharge have been substantially changed by development. Groundwater pumpage in 1980 was about five times the value of predevelopment regional recharge. About 80% of the pumpage was supplied from increased regional recharge. Also resistance to vertical flow caused by many fine-grained beds within the permeable zones can be as important as resistance caused by regional confining units.

Gulf Coast aquifer systems↗

Design, revision, and application of ground-water flow models for simulation of selected water-management scenarios in the coastal area of Georgia and adjacent parts of South Carolina and Florida

Ground-water flow models of the Floridan aquifer system in the coastal area of Georgia and adjacent parts of South Carolina and Florida, were revised and updated to ensure consistency among the various models used, and to facilitate evaluation of the effects of pumping on the ground-water level near areas of saltwater contamination. The revised models, developed as part of regional and areal assessments of ground-water resources in coastal Georgia, are--the Regional Aquifer-System Analysis (RASA) model, the Glynn County area (Glynn) model, and the Savannah area (Savannah) model. Changes were made to hydraulic-property arrays of the RASA and Glynn models to ensure consistency among all of the models; results of theses changes are evidenced in revised water budgets and calibration statistics. Following revision, the three models were used to simulate 32 scenarios of hypothetical changes in pumpage that ranged from about 82 million gallons per day (Mgal/d) lower to about 438 Mgal/d higher, than the May 1985 pumping rate of 308 Mgal/d. The scenarios were developed by the Georgia Department of Natural Resources, Environmental Protection Division and the Chatham County-Savannah Metropolitan Planning Commission to evaluate water-management alternatives in coastal Georgia. Maps showing simulated ground-water-level decline and diagrams presenting changes in simulated flow rates are presented for each scenario. Scenarios were grouped on the basis of pumping location--entire 24-county area, central subarea, Glynn-Wayne-Camden County subarea, and Savannah-Hilton Head Island subarea. For those scenarios that simulated decreased pumpage, the water level at both Brunswick and Hilton Head Island rose, decreasing the hydraulic gradient and reducing the potential for saltwater contamination. Conversely, in response to scenarios of increased pumpage, the water level at both locations declined, increasing the hydraulic gradient and increasing the potential for saltwater contamination. Pumpage effects on ground-water levels and related saltwater contamination at Brunswick and Hilton Head Island generally diminish with increased distance from these areas. Additional development of the Upper Floridan aquifer may be possible in parts of the coastal area without affecting saltwater contamination at Brunswick or Hilton Head Island, due to the presence of two hydrologic boundaries--the Gulf Trough, separating the northern and central subareas; and the hypothesized Satilla Line, separating the central and southern subareas. These boundaries diminish pumpage effects across them; and may enable greater ground-water withdrawal in areas north of the Gulf Trough and south of the Satilla Line without producing appreciable drawdown at Brunswick or Hilton Head Island.

Florida, Georgia, South Carolina↗

Ground-water age and quality in the High Plains Aquifer near Seward, Nebraska, 2003-04

The U.S. Geological Survey, in cooperation with the City of Seward, Nebraska, conducted a study of ground-water age and quality to improve understanding of: (1) traveltimes from recharge areas to public-supply wells, (2) the effects of geochemical reactions in the aquifer on water quality, and (3) how water quality has changed historically in response to land-use practices. Samples were collected from four supply wells in the Seward west well field and from nine monitoring wells along two approximate ground-water flow paths leading to the well field. Concentrations of three different chlorofluorocarbons (CFC-12, CFC-11, and CFC-113), sulfur hexafluoride (SF 6 ), and ratios of tritium ( 3 H) to helium-3 ( 3 He) isotope derived from radioactive decay of 3H were used to determine the apparent recharge age of ground-water samples. Age interpretations were based primarily on 3 H/ 3 He and CFC-12 data. Estimates of apparent ground-water age from tracer data were complicated by mixing of water of different ages in 10 of the 13 ground-water samples collected. Apparent recharge dates of unmixed ground-water samples or mean recharge dates of young fractions of mixed water in samples collected from monitoring wells ranged from 1985 to 2002. For monitoring-well samples containing mixed water, the fraction of the sample composed of young water ranged from 26 to 77 percent of the sample. Apparent mean recharge dates of young fractions in samples collected from four supply wells in the Seward west well field ranged from about 1980 to 1990. Estimated fractions of the samples composed of young water ranged from 39 to 54 percent. It is implicit in the mixing calculations that the remainder of the sample that is not young water is composed of water that is more than 60 years old and contains no detectable quantities of modern atmospheric tracers. Estimated fractions of the mixed samples composed of "old" water ranged from 23 to 74 percent. Although alternative mixing models can be used to interpret the results, the mean age and mixing fractions from the primary mixing models used were fairly similar. Relations of ground-water age and nitrate concentrations to depth were not consistent across the study area. In some well nests, more young water and nitrate were present near the bottom than in the middle of the aquifer. These results probably reflect pumping from irrigation and supply wells, which are screened primarily in the lower part of the aquifer, and draw younger water downward in the aquifer. Substantial mixing probably occurs because the aquifer is relatively thin (50 feet) and has a relatively high density of wells (about five pumping wells per square mile). The most reliable estimate of horizontal traveltimes based on differences in ground-water ages between a shallow monitoring well at the upgradient end of the northwest well transect and the deep well at the downgradient end of the well transect was 9 years to travel a distance of about 2 miles. The general similarity of ages at similar depths between different well nests is consistent with the fact that horizontal flow in the aquifer is relatively rapid. Concentrations of nitrate (as nitrogen) in untreated ground-water samples from supply wells in the well field were larger than the U.S. Environmental Protection Agency Maximum Contaminant Level for drinking water of 10 mg/L (milligrams per liter), ranging from 11.3 to 13.5 mg/L. It is unlikely that nitrate concentrations in the aquifer near the Seward west well field are decreased by denitrification in the aquifer due to oxic geochemical conditions that preclude this reaction. Nitrate concentrations coupled with water recharge dates were compared to historical estimated fertilizer application in an attempt to reconstruct historical trends in ground-water nitrate concentrations and their relation to land-use practices. Nitrate concentrations in young-water fractions, after adjustment for mixing, may be decreasing over apparent recharge dates of 1980 to 2002, corresponding to a period of generally decreasing nitrogen fertilizer applications.

Nebraska↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Three-dimensional model simulation of transient ground-water flow in the Albuquerque-Belen Basin, New Mexico

A three-dimensional digital model that simulates transient flow in the alluvial aquifer system underlying the Albuquerque-Belen Basin, New Mexico, was constructed as part of a regional aquifer study of the southwest alluvial basins. The model simulates hydraulic heads and changes in hydraulic heads for 1907 to 1979. Hydraulic-conductivity values used in the accepted model range from 0.25 ft/day in part of the Santa Fe Group to 50 ft/day in the fluvial deposits in the Rio Grande flood plain. The majority of the basin-fill material of the Santa Fe group of Tertiary and Quaternary age was modeled as having a horizontal hydraulic conductivity of either 30 or 40 ft/day. The simulated specific storage of the aquifer was 0.000001/ft and the simulated specific yield was 0.10. The aquifer was simulated as being vertically anisotropic with a ratio of vertical to horizontal hydraulic conductivity of 1:500. Simulations for 1976-79 indicated that of the 100,000 acre-ft of groundwater withdrawn annually from the basin-fill deposits outside of the Rio Grande flood plain, 68% was obtained from recharge around the basin margin, depletion of streams that are tributary to the Rio Grande, and the stream-aquifer system in the Rio Grande flood plain. Depletion of aquifer storage accounted for 25% of the groundwater supply to wells outside of the flood plain, and the remaining 7% was obtained by induced groundwater inflow from the Santo Domingo Basin. The model displayed an acceptable performance throughout the period of simulation. However, by the end of the simulation period, 1979, the portrayal of the Rio Grande flood-plain system as a specified hydraulic-head boundary was having adverse effects on the simulation. (Author 's abstract)

Water-Resources Investigations Report↗

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↗