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.

SourceResearch connectionFrom yearThrough year

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

At least 973 records · Page 54Linked to original sources

USGS1997

Hydrogeology and water quality of a surficial aquifer underlying an urban area, Manchester, Connecticut

The quality of water along flowpaths in a surficial aquifer system in Manchester, Connecticut, was studied during 1993-95 as part of the National Water Quality Assessment program. The flowpath study examined the relations among hydrogeology, land-use patterns, and the presence of contaminants in a surficial aquifer in an urban area, and evaluated ground water as a source of contamination to surface water. A two-dimensional, finite-difference groundwater- flow model was used to estimate travel distance, which ranged from about 50 to 11,000 feet, from the source areas to the sampled observation wells. Land use, land cover, and population density were determined in the source areas delineated by the ground-water-flow simulation. Source areas to the wells contained either high- or medium-density residential areas, and population density ranged from 629 to 8,895 people per square mile. Concentrations of selected inorganic constituents, including sodium, chloride, and nitrite plus nitrate nitrogen, were higher in the flowpath study wells than in wells in undeveloped areas with similar aquifer materials. One or more of 9 volatile organic compounds were detected at 12 of 14 wells. The three most commonly detected volatile organic compounds were chloroform, methyl-tert-butyl ether, and trichloroethene. Trichloroethene was detected at concentrations greater than the maximum contaminant level for drinking water (5 micrograms per liter) in samples from one well. Four pesticides, including dichloro diphenyl dichloroethylene, dieldrin, dichloroprop, and simazine were detected at low concentrations. Concentrations of sodium and chloride were higher in samples collected from wells screened in the top of the saturated zone than in samples collected from deeper zones. Volatile organic compounds and elevated concentrations of nitrite plus nitrate as nitrogen were detected at depths of as much as 60 feet below the water table, indicating that the effects of human activities on the ground-water quality extends to the bottom of the surficial aquifer. The age of ground water, as determined by tritium and 3helium concentrations, was 0.9 to 22.6 years. pH, alkalinity, and calcium were higher and concentrations of dissolved oxygen were lower in ground-water samples with ages of 10 years or more than in samples younger than 10 years. In addition, concentrations of sodium, chloride, and nitrite plus nitrate nitrogen were low in ground-water samples with ages of 10 years or more, indicating that concentrations of these compounds may be increasing with time or that the recharge areas to these wells may have had less intensive urban land use. Methyl-tert-butyl ether was detected only in wells with ground water ages of less than 11 years, which is consistent with the date of introduction of this compound as a gasoline additive in Connecticut. Analysis of additional samples collected for analysis of stable nitrogen isotopes indicated that the most likely source of elevated concentrations of nitrate nitrogen was lawn and garden fertilizers, but other sources, including wastewater effluents, soil organic nitrogen, and atmospheric deposition, may contribute to the total. Population density was positively correlated (at the 97 percent confidence level) to concentrations of nitrite plus nitrate as nitrogen. Water quality in the Hockanum River aquifer has been degraded by human activities, and, after discharge to surface water, affects the water quality in the Hockanum River. On an annual basis, ground-water discharge from the study area to the river (as measured at a downstream continuous-record gaging station) contributes about 5 percent of the annual load of nitrite plus nitrate nitrogen, but, during low flow, contributes 11 percent of the nitrite plus nitrate nitrogen, 32 percent of the calcium, and 16 percent of the chloride to the river.

Connecticut↗

USGS1998

Analytical methods, numerical modeling, and monitoring strategies for evaluating the effects of ground-water withdrawals on unconfined aquifers in the New Jersey Coastal Plain

Analytical and numerical solutions of ground-water withdrawals in the unconfined part of the Kirkwood-Cohansey aquifer system of the Coastal Plain of New Jersey were evaluated for their usefulness in predicting the area of influence of a pumped well and in determining hydraulic characteristics of an aquifer. Additionally, simulations of ground-water withdrawal using a finite-difference model provided information on the ways in which prudent well-location strategies can disperse the local effects of withdrawal over a larger part of an aquifer system. The design of a monitoring network that is sensitive to the ground-water hydraulics of streams and wetlands of the Coastal Plain of New Jersey also was considered for its utility in providing hydrologic data necessary to establish the baseline hydrologic conditions near wetlands and streams and in signaling when ground-water levels are being adversely affected by withdrawals elsewhere in the system. The application of methods based on the Theis analytical solution to ground-water flow in unconfined aquifers can lead to erroneous estimates of the size of the area of influence generated by ground-water withdrawals. Analysis oftime-drawdown data from an unconfined aquifer system are best evaluated by means of the Neuman solution, which accounts for the effects of gravity drainage; however, the pumped well must be far enough from streams so that ground water is not drawn from nearby streams. Time-drawdown data from a test well in Winslow County, N.J., were analyzed by means of the Neuman solution. Results indicate that the aquifer has a relatively high vertical to horizontal anisotropy of 1:198, and a specific yield of 0.028, an indication that the area of influence of a pumped well at the test site would be relatively large. Results from a finite-difference ground-water-flow model of the northeastern part of the Mullica River Basin near Chesilhurst, N.J., show that the area influenced by a long-term withdrawal is best estimated from a steady state ground-water-flow analysis that includes the effects of average areal recharge. Withdrawal simulations indicate an order-of-magnitude difference between the size of the area of influence generated from a 3-day (72 hour) withdrawal and the size of the area produced under steady-state conditions. An aquifer characterized by a low specific yield will cause the area of influence to extend farther away from the pumped well. The contributing area of flow to the pumped well includes areas on the water table that would, under natural conditions, be incorporated into the contributing areas of flow to streams. Ground water that is drawn to a pumped well is diverted from nearby streams; the withdrawal decreases the size of the contributing areas of flow to streams by an amount equal to the contributing area of flow to the well. Withdrawals made from a well close to a stream divert ground water that would, under natural conditions, flow to the stream. The diverted ground water causes the area of influence of the well to be smaller than it would if the well were far from the stream. Water-table declines caused by withdrawals near streams are, to some degree, mitigated by ground-water diversion from streams. However, the withdrawals can significantly reduce ground-water seepage to nearby streams, especially along stream reaches and wetlands close to the well. Alternatively, these effects can be dispersed over a large part of the aquifer if wells are located on surface-water divides. Measurements of seasonal water-level fluctuations in the Mullica River Basin indicate that the greatest fluctuations in water levels are found in upland areas, where the average fluctuation is 3.4 feet. Fluctuations in hydraulic head in the wetland areas averages 1.3 feet. The bimodal average of ranges in water levels show that upland areas are more sensitive to recharge than lowland areas. The pattern of yearly mean water levels fluctuates irregularly about a long-term mean value. Abnormally low or high yearly average values that are brought on by periods of drought or excess recharge are short lived; over time, hydrologic conditions shift back to average levels under natural conditions. Wetland areas in the New Jersey Coastal Plain are characterized by ground-water seepage into wide, shallow depressions. Periods of inundation are longest in the deepest part of the depression, whereas inundation of areas near the fringes of wetlands due to ground-water seepage is only seasonal. The seepage face in the fringe areas expand and contract in response to seasonal variation in water-table elevation and in response to precipitation. Values of the aquifer storage coefficient and transmissivity can, in some cases, be determined by use of hydraulic head or streamflow recession analysis as an alternative to aquifer testing. The recession curves developed from hydro graphs of Middle Branch and McDonalds Branch in the New Jersey Coastal Plain indicate that the aquifer near McDonalds Branch has about 2.6 times the storage capacity of the aquifer adjacent to Middle Branch; this finding is consistent with the relatively small ranges of water-level changes measured in McDonalds Branch compared to those measured in Middle Branch.

New Jersey↗

USGS2008

Simulation of Ground-Water Flow and Effects of Ground-Water Irrigation on Base Flow in the Elkhorn and Loup River Basins, Nebraska

Irrigated agriculture is vital to the livelihood of communities in the Elkhorn and Loup River Basins in Nebraska, and ground water is used to irrigate most of the cropland. Concerns about the sustainability of ground-water and surface-water resources have prompted State and regional agencies to evaluate the cumulative effects of ground-water irrigation in this area. To facilitate understanding of the effects of ground-water irrigation, a numerical computer model was developed to simulate ground-water flow and assess the effects of ground-water irrigation (including ground-water withdrawals, hereinafter referred to as pumpage, and enhanced recharge) on stream base flow. The study area covers approximately 30,800 square miles, and includes the Elkhorn River Basin upstream from Norfolk, Nebraska, and the Loup River Basin upstream from Columbus, Nebraska. The water-table aquifer consists of Quaternary-age sands and gravels and Tertiary-age silts, sands, and gravels. The simulation was constructed using one layer with 2-mile by 2-mile cell size. Simulations were constructed to represent the ground-water system before 1940 and from 1940 through 2005, and to simulate hypothetical conditions from 2006 through 2045 or 2055. The first simulation represents steady-state conditions of the system before anthropogenic effects, and then simulates the effects of early surface-water development activities and recharge of water leaking from canals during 1895 to 1940. The first simulation ends at 1940 because before that time, very little pumpage for irrigation occurred, but after that time it became increasingly commonplace. The pre-1940 simulation was calibrated against measured water levels and estimated long-term base flow, and the 1940 through 2005 simulation was calibrated against measured water-level changes and estimated long-term base flow. The calibrated 1940 through 2005 simulation was used as the basis for analyzing hypothetical scenarios to evaluate the effects of ground-water irrigation on stream base flow for 1940 through 2005 and for 2006 through 2045. Simulated base flows were compared for scenarios that alternately did or did not include a representation of the effects of ground-water irrigation. The difference between simulated base flows for the two scenarios represents the predicted effects of ground-water irrigation on base flow. Comparison of base flows between simulations with ground-water irrigation and no ground-water irrigation indicated that ground-water irrigation has cumulatively reduced streamflows from 1940 through 2005 by 888,000 acre-feet in the Elkhorn River Basin and by 2,273,000 acre-feet in the Loup River Basin. Generally, predicted cumulative effects of ground-water irrigation on base flow were 5 to 10 times larger from 2006 through 2045 than from 1940 through 2005, and were 7,678,000 acre-feet for the Elkhorn River Basin and 14,784,000 acre-feet for the Loup River Basin. The calibrated simulation also was used to estimate base-flow depletion as a percentage of pumping volumes for a 50-year future time period, because base-flow depletion percentages are used to guide the placement of management boundaries in Nebraska. Mapped results of the base-flow depletion analysis conducted for most of the interior of the study area indicated that pumpage of one additional theoretical well simulated for a future 50-year period generally would result in more than 80 percent depletion when it was located close to the stream, except in areas where depletion was partly offset by reduced ground-water discharge to evapotranspiration in wetland areas. In many areas, depletion for the 50-year future period composed greater than 10 percent of the pumped water volume for theoretical wells placed less than 7 or 8 miles from the stream, though considerable variations existed because of the heterogeneity of the natural system represented in the simulation. For a few streams, predicted future simulated base flows dec

Scientific Investigations Report↗

USGS1978

The hydrothermal system of Long Valley Caldera, California

Long Valley caldera, an elliptical depression covering 450 km 2 on the eastern front of the Sierra Nevada in east-central California, contains a hot-water convection system with numerous hot springs and measured and estimated aquifer temperatures at depths of 180?C to 280?C. In this study we have synthesized the results of previous geologic, geophysical, geochemical, and hydrologic investigations of the Long Valley area to develop a generalized conceptual and mathematical model which describes the gross features of heat and fluid flow in the hydrothermal system. Cenozoic volcanism in the Long Valley region began about 3.2 m.y. (million years) ago and has continued intermittently until the present time. The major event that resulted in the formation of the Long Valley caldera took place about 0.7 m.y. ago with the eruption of 600 km 3 or more of Bishop Tuff of Pleistocene age, a rhyolitic ash flow, and subsequent collapse of the roof of the magma chamber along one or more steeply inclined ring fractures. Subsequent intracaldera volcanism and uplift of the west-central part of the caldera floor formed a subcircular resurgent dome about 10 km in diameter surrounded by a moat containing rhyolitic, rhyodacitic, and basaltic rocks ranging in age from 0.5 to 0.05 m.y. On the basis of gravity and seismic studies, we estimate an aver- age thickness of fill of 2.4 km above the precaldera granitic and metamorphic basement rocks. A continuous layer of densely welded Bishop Tuff overlies the basement rocks, with an average thickness of 1.4 km; the fill above the welded Bishop Tuff consists of intercalated volcanic flows and tuffs and fluvial and lacustrine deposits. Assuming the average grain density of the fill is between 2.45 and 2.65 g/cm 3 , we calculate the average bulk porosity of the total fill as from 0.11 to 0.21. Comparison of published values of porosity of the welded Bishop Tuff exposed southeast of the caldera with calculated values indicates average bulk porosity for the welded tuff (including fracture porosity) from 0.05 to 0.10. Because of its continuity and depth and the likelihood of significant fracture permeability in the more competent rocks such as the welded tuff, our model of the hydrothermal system assumes that the Bishop Tuff provides the principal hot-water reservoir. However, because very little direct information exists from drill holes below 300 m, this assumption must be considered tentative. Long Valley caldera is drained by the Owens River and several tributaries which flow into Lake Crowley in the southeast end of the caldera. Streamflow and springflow measurements for water years 1964-74 indicate a total inflow to Lake Crowley of about 10,900 L/s. In contrast, the total discharge of hot water from the hydrothermal reservoir is about 300 L/s. For modeling purposes, the ground-water system is considered as comprising a shallow subsystem in the fill above the densely welded Bishop Tuff containing relatively cold ground water, and a deep subsystem or hydrothermal reservoir in the welded tuff containing relatively hot ground water. Hydrologic, isotopic, and thermal data indicate that recharge to the hydrothermal reservoir occurs in the upper Owens River drainage basin along the western periphery of the caldera. Temperature profiles in a 2.11- km-deep test well drilled by private industry in the southeastern part of the caldera suggest that an additional flux of relatively cool ground water recharges the deep subsystem around the northeast rim. Flow in the shallow ground-water subsystem is neglected in the model except in recharge areas and along Hot Creek gorge, where approximately 80 percent of the hot-water discharge from the hydrothermal reservoir moves upward along faults toward springs in the gorge. Heat-flow data from the Long Valley region indicate that the resurgent dome overlies a residual magma chamber more circular in plan than the original magma chamber that supplied the Bishop Tuff

California↗

USGS2002

Ground-water flow and numerical simulation of recharge from streamflow infiltration near Pine Nut Creek, Douglas County, Nevada

Ground-water flow and recharge from infiltration near Pine Nut Creek, east of Gardnerville, Nevada, were simulated using a single-layer numerical finite-difference model as part of a study made by the U.S. Geological Survey in cooperation with the Carson Water Subconservancy District. The model was calibrated to 190 water-level measurements made in 27 wells in December 2000, and in 9 wells from August 1999 through April 2001. The purpose of this study was to estimate reasonable limits for the approximate volume of water that may be stored by recharge through infiltration basins, and the rate at which recharged water would dissipate or move towards the valley floor. Measured water levels in the study area show that infiltration from the Allerman Canal and reservoir has created a water-table mound beneath them that decreases the hydraulic gradient east of the canal and increases the gradient west of the canal. North of Pine Nut Creek, the mound causes ground water to flow toward the northern end of the reservoir. South of Pine Nut Creek, relatively high water levels probably are maintained by the mound beneath the Allerman Canal and possibly by greater rates of recharge from the southeast. Water-level declines near Pine Nut Creek from August 1999 through April 2001 probably are caused by dissipation of recharge from infiltration of Pine Nut Creek streamflow in the springs of 1998 and 1999. Using the calibrated model, a simulation of recharge through a hypothetical infiltration basin covering 12.4 acres near Pine Nut Creek applied 700 acre-feet per year of recharge over a six-month period, for a total of 3,500 acre-feet after 5 consecutive years. This recharge requires a diversion rate of about 2 cubic feet per second and an infiltration rate of 0.3 foot per day. The simulations showed that recharge of 3,500 acre-feet caused water levels near the basin to rise over 70 feet, approaching land surface, indicating 3,500 acre-feet is the maximum that may be stored in a 5-year period, given the basin location and surface area used in the simulations. Greater amounts probably could be stored if separate infiltration basins were installed at different locations along the Pine Nut Creek alluvial fan, applying the recharge over a larger area. The water-table mound resulting from recharge extended 7,000 feet north, west, and south of the infiltration basin. After recharge ceased, water levels near the center of the mound declined rapidly to within 20 feet of initial levels after 2 years, and within 10 feet of initial levels after 7 years. The recharge mound dissipates laterally across the modeled area at decreasing rates over time. A water-level rise of 1 foot moved westward towards the valley floor 660 feet from peak conditions after 1 year, and averaged 550 feet, 440 feet, and 330 feet per year for the periods 1-4, 4-7, and 7-10 years, respectively, after recharge ceased. Simulations of subsequent pumping from hypothetical wells near the infiltration basin were made by applying pumping near the basin beginning 1 year after recharge of 3,500 acre-feet ceased. Pumping was applied over a 6-month period for 4 years from one well at 400 acre-feet per year, withdrawing 1,600 acre-feet or 45 percent of that recharged, and from two wells totaling 800 acre-feet per year, withdrawing 3,200 acre-feet or 90 percent of that recharged. Pumping of 1,600 acre-feet caused water-levels near the infiltration basin to decline only slightly below initial levels. Pumping of 3,200 acre-feet caused water-levels near the infiltration basin to decline a maximum of 30 feet below initial levels, with smaller declines extending laterally in all directions for 4,000 feet from the pumping wells. Water-level declines are a result of pumping at a rate sufficient to withdraw the majority of the water recharged through the infiltration basin. Although the declines may affect water levels in nearby domestic wells, the simulations show that water levels recover quickly after

Water-Resources Investigations Report↗

USGS1998

Ground-water hydrology and simulation of ground-water flow at Operable Unit 3 and surrounding region, U.S. Naval Air Station, Jacksonville, Florida

The Naval Air Station, Jacksonville (herein referred to as the Station), occupies 3,800 acres adjacent to the St. Johns River in Duval County, Florida. Operable Unit 3 (OU3) occupies 134 acres on the eastern side of the Station and has been used for industrial and commercial purposes since World War II. Ground water contaminated by chlorinated organic compounds has been detected in the surficial aquifer at OU3. The U.S. Navy and U.S. Geological Survey (USGS) conducted a cooperative hydrologic study to evaluate the potential for ground water discharge to the neighboring St. Johns River. A ground-water flow model, previously developed for the area, was recalibrated for use in this study. At the Station, the surficial aquifer is exposed at land surface and forms the uppermost permeable unit. The aquifer ranges in thickness from 30 to 100 feet and consists of unconsolidated silty sands interbedded with local beds of clay. The low-permeability clays of the Hawthorn Group form the base of the aquifer. The USGS previously conducted a ground-water investigation at the Station that included the development and calibration of a 1-layer regional ground-water flow model. For this investigation, the regional model was recalibrated using additional data collected after the original calibration. The recalibrated model was then used to establish the boundaries for a smaller subregional model roughly centered on OU3. Within the subregional model, the surficial aquifer is composed of distinct upper and intermediate layers. The upper layer extends from land surface to a depth of approximately 15 feet below sea level; the intermediate layer extends from the upper layer down to the top of the Hawthorn Group. In the northern and central parts of OU3, the upper and intermediate layers are separated by a low-permeability clay layer. Horizontal hydraulic conductivities in the upper layer, determined from aquifer tests, range from 0.19 to 3.8 feet per day. The horizontal hydraulic conductivity in the intermediate layer, determined from one aquifer test, is 20 feet per day. An extensive stormwater drainage system is present at OU3 and the surrounding area. Some of the stormwater drains have been documented to be draining ground water from the upper layer of the surficial aquifer, whereas other drains are only suspected to be draining ground water. The subregional model contained 78 rows and 148 columns of square model cells that were 100 feet on each side. Vertically, the surficial aquifer was divided into two layers; layer 1 represented the upper layer and layer 2 represented the intermediate layer. Steady-state ground-water flow conditions were assumed. The model was calibrated to head data collected on October 29 and 30, 1996. After calibration, the model matched all 67 measured heads to within the calibration criterion of 1 foot; and 48 of 67 simulated heads (72 percent) were within 0.5 foot. Model simulated recharge rates ranged from 0.4 inch per year in areas that were largely paved to 13.0 inches per year in irrigated areas. Simulated hydraulic conductivities in the upper layer at OU3 ranged from 0.5 foot per day in the north to 1.0 foot per day in the south. Simulated vertical leakance between the upper and intermediate layers ranged from 1.0x10-6 per day in an area with low-permeability clays to 4.3x10-2 per day in an area that had been dredged. Simulated transmissivities in the intermediate layer ranged from 25 feet squared per day in an area of low-permeability channel-fill deposits to a high of 1,200 feet squared per day in areas covering most of OU3. Simulated riverbed conductances ranged from 4 to 60 feet squared per day and simulated bottom conductances of leaking stormwater drains ranged from 5 to 20 feet squared per day. The direction and velocity of ground-water flow was determined using particle-tracking techniques. Ground-water flow in the upper layer was generally eastward toward the St. Johns River. However, leaking stormwater drains locally modified the flow system to create small areas with flow that was diverted to the drains. The flow velocities in the upper layer at OU3 were slow, averaging about 2 feet per year. The slow velocities were primarily the result of the low horizontal hydraulic conductivity and, secondarily, the result of the low recharge rate. The simulated rate at which ground water leaked into the stormwater drains was low, averaging about 0.0011 cubic feet per second per 100 feet of stormwater drainage conduit. Ground-water flow in the intermediate layer moved eastward toward and discharged into the St. Johns River. Flow velocities were significantly higher in this layer than in the upper layer. The velocity was about 35 and 12 feet per year in the northern and southern parts of OU3, respectively.

Florida↗

USGS2004

Hydrogeology, water quality, and ecology of Anderton Branch near the Quail Hollow Landfill, Bedford County, Tennessee, 1995-99

The Quail Hollow Landfill, located in southeastern Bedford County on the Highland Rim overlooking the Central Basin karst region of Tennessee, is constructed on the gravelly, clay-rich residuum of the Fort Payne Formation of Mississippian age. A conceptual hydrologic model of the landfill indicated that Anderton Branch was at risk of being affected by the landfill. Ground water flowing beneath the landfill mixes with percolating rainwater that has passed through the landfill and discharges to the surface from numerous weeps, seeps, and springs present in the area. Anderton Branch, adjacent to the landfill site on the north and east, receives most of the discharge from these weeps, seeps, and springs. Anderton Branch also receives water from the Powell Branch drainage basin to the west and south because of diverted flow of ground water through Harrison Spring Cave. The U.S. Geological Survey, in cooperation with the Bedford County Solid Waste Authority, conducted a study to evaluate the effect of the Quail Hollow Landfill on ground- and surface-water quality. During storm runoff, specific conductance was elevated, and cadmium, iron, manganese, lead, and nickel concentrations in Anderton Branch frequently exceeded maximum contaminant levels for drinking water for the State of Tennessee. High chloride inputs to Anderton Branch were detected at two locations?a barnyard straddling the stream and a tributary draining a pond that receives water directly from the landfill. The chloride inputs probably contribute to chloride load levels that are three times higher for Anderton Branch than for the control stream Anthony Branch. Although toxic volatile organic compounds were detected in water from monitoring wells at the landfill, no organic contaminants were detected in domestic water wells adjacent to the landfill or in Anderton Branch. Sons Spring, a karst spring near the landfill, has been affected by the landfill as indicated by an increase in chloride concentrations from 4 milligrams per liter in 1974 to 59 milligrams per liter in 1996. Analysis of water samples from Sons Spring detected concentrations of nickel that exceeded primary drinking-water standards and Tennessee Department of Environment and Conservation fish and aquatic life chronic standards. Trichloroethene, 1,1-dichloroethene, and 1,1-dichloroethane also were detected at Sons Spring. The presence of these chlorinated solvents imply the landfill origin of the contaminants in Sons Spring. Continuous monitoring at Sons Spring indicated a pattern of decreased specific conductance and lower contaminant concentrations after a storm. Contaminant concentrations increased with specific conductance to pre-storm levels after several days. The benthic macroinvertebrate community in Anderton Branch adjacent to the landfill was not different from the communities at control sites upstream and in Anthony Branch. Sons Spring, however, has low abundance and numbers of benthic macroinvertebrate taxa. Toxicity studies using Ceriodaphnia dubia indicated no toxicity in the base flow or storm water in Anderton Branch or in a tributary draining a pond that receives water from the landfill and Sons Spring; however, water collected from Sons Spring resulted in 100 percent mortality to all organisms within 48 hours. High concentrations of nickel were detected in crayfish tissue from control sites and Anderton Branch. Analysis of sediment samples also indicates nickel concentrations are high at control sites upstream of the landfill. Increased levels of the biomarker metallothionein detected in crayfish from Anderton Branch likely are not caused by nickel or cadmium because the levels present in the tissue are not correlated with metallothionein levels. Despite the high levels of certain metals in Anderton Branch during storm flow, the lack of toxicity and the health of the benthic community imply no detectable negative effect from the landfill to the stream. Sons Spring, howe

Scientific Investigations Report↗

USGS2009

Simulations of Groundwater Flow and Particle Tracking Analysis in the Area Contributing Recharge to a Public-Supply Well near Tampa, Florida, 2002-05

Shallow ground water in the north-central Tampa Bay region, Florida, is affected by elevated nitrate concentrations, the presence of volatile organic compounds, and pesticides as a result of groundwater development and intensive urban land use. The region relies primarily on groundwater for drinking-water supplies. Sustainability of groundwater quality for public supply requires monitoring and understanding of the mechanisms controlling the vulnerability of public-supply wells to contamination. A single public-supply well was selected for intensive study based on the need to evaluate the dominant processes affecting the vulnerability of public-supply wells in the Upper Floridan aquifer in the City of Temple Terrace near Tampa, Florida, and the presence of a variety of chemical constituents in water from the well. A network of 29 monitoring wells was installed, and water and sediment samples were collected within the area contributing recharge to the selected public-supply well to support a detailed analysis of physical and chemical conditions and processes affecting the water chemistry in the well. A three-dimensional, steady-state groundwater flow model was developed to evaluate the age of groundwater reaching the well and to test hypotheses on the vulnerability of the well to nonpoint source input of nitrate. Particle tracking data were used to calculate environmental tracer concentrations of tritium and sulfur hexafluoride and to calibrate traveltimes and compute flow paths and advective travel times in the model area. The traveltime of particles reaching the selected public-supply well ranged from less than 1 day to 127.0 years, with a median of 13.1 years; nearly 45 percent of the simulated particle ages were less than about 10 years. Nitrate concentrations, derived primarily from residential/commercial fertilizer use and atmospheric deposition, were highest (2.4 and 6.11 milligrams per liter as nitrogen, median and maximum, respectively) in shallow groundwater from the surficial aquifer system and lowest (less than the detection level of 0.06 milligram per liter) in the deeper Upper Floridan aquifer. Denitrification occurred near the interface of the surficial aquifer system and the underlying intermediate confining unit, within the intermediate confining unit, and within the Upper Floridan aquifer because of reducing conditions in this part of the flow system. However, simulations indicate that the rapid movement of water from the surficial aquifer system to the selected public-supply well through karst features (sinkholes) and conduit layers that bypass the denitrifying zones (short-circuits), coupled with high pumping rates, allow nitrate to reach the selected public-supply well in concentrations that resemble those of the overlying surficial aquifer system. Water from the surficial aquifer system with elevated concentrations of nitrate and low concentrations of some volatile organic compounds and pesticides is expected to continue moving into the selected public-supply well, because calculated flux-weighted concentrations indicate the proportion of young affected water contributing to the well is likely to remain relatively stable over time. The calculated nitrate concentration in the selected public-supply well indicates a lag of 1 to 10 years between peak concentrations of nonpoint source contaminants in recharge and appearance in the well.

Scientific Investigations Report↗

USGS2007

Simulated effects of projected 2010 withdrawals on ground-water flow and water levels in the New Jersey coastal plain – A task of the New Jersey Water Supply Plan, 2006 revision

A ground-water flow model previously developed as part of a Regional Aquifer System Analysis (RASA) of the New Jersey Coastal Plain was used to simulate ground-water flow in eight major confined aquifers to help evaluate ground-water resources in support of the New Jersey Department of Environmental Protection's revision of the New Jersey State Water Supply Plan. This model was calibrated to 1998 steady-state and transient conditions. Withdrawals at wells in operation in 1998 were varied in three scenarios to evaluate their effects on flow directions, water levels, and water budgets in the confined aquifers. The scenarios used to predict changes in pumpage from 1998 to 2010 were based on (1) a continuation of 1990-99 trends in water use, (2) public-supply withdrawals estimated from county population projections, and (3) restricted withdrawals in Water-Supply Critical Areas. Total withdrawals in these three scenarios were approximately 366, 362, and 355 million gallons per day, respectively. The results of these simulations are used by New Jersey water-management officials to help address water-supply concerns for the State. In the revision of the New Jersey State Water Supply Plan, the eight major confined aquifers of the New Jersey Coastal Plain and their outcrop areas are divided into 41 hydrologic budget areas (HBAs). Simulation results were used to assess the effects of changing ground-water withdrawals on water levels and the flow budgets in each budget area. Simulation results for each scenario were compared with 1998 (baseline) simulated water levels and flow budgets. The 41 hydrologic budget areas are in areas of large ground-water withdrawals, water-level declines, and (or) saltwater-intrusion potential. Their boundaries are based on various hydrologic, geohydrologic, and withdrawal conditions, such as aquifer extent, location of the 250-milligram-per-liter isochlor, aquifer outcrop area, and ground-water divides. The budget areas include primarily the onshore, freshwater portions of the aquifers. A budget analysis was done for each of the hydrologic budget areas for each scenario. Ground-water withdrawals, leakage to streams, net leakage to overlying and underlying aquifers, lateral flow to adjacent budget areas, and the flow direction at the 250-milligram-per-liter isochlor were evaluated. Although three different methods were applied to predict future pumping rates, the simulated water levels for scenarios 1 and 2 were generally within 2 feet of each other in most areas in the confined aquifers, but differences of more than 2 feet occurred locally. Differences in values of flow-budget components between scenarios 1 and 2 as a percentage change from 1998 values were generally within 2 percent in most hydrologic budget areas, but values of some budget components in some hydrologic budget areas differed by more than 2 percent. Simulated water levels recovered as much as 4 feet more in northeastern Camden and northwestern Burlington Counties in the Lower Potomac-Raritan-Magothy aquifer, and as much as 3 feet more in the same area in the Upper and Middle Potomac-Raritan-Magothy aquifers when pumpage restrictions were imposed in Critical Area 2 (scenario 3). In the Wenonah-Mount-Laurel aquifer, water levels declined continually in Monmouth County (HBA 8) downdip from the outcrop (in Critical Area 1) from 1988 to 2010 in all three scenarios, although most of the water levels farther downdip from this area in Critical Area 1 are still recovering because of mandated reductions in pumpage in the 1990s. In the Englishtown aquifer system, water levels declined continually in small areas in HBA 13 in central Monmouth County (in Critical Area 1) and in western Monmouth County downdip from the outcrop from 1988 to 2010 in all three scenarios, although most of the water levels farther downdip from this area are still recovering because of the mandated reductions in pumpage. In the Upper Potomac-Raritan-Magothy aquif

New Jersey↗

USGS1997

Vulnerability of ground water to atrazine leaching in Kent County, Michigan

A steady-state model of pesticide leaching through the unsaturated zone was used with readily available hydrologic, lithologic, and pesticide characteristics to estimate the vulnerability of the near-surface aquifer to atrazine contamination from non-point sources in Kent County, Michigan. The modelcomputed fraction of atrazine remaining at the water table, RM , was used as the vulnerability criterion; time of travel to the water table also was computed. Model results indicate that the average fraction of atrazine remaining at the water table was 0.039 percent; the fraction ranged from 0 to 3.6 percent. Time of travel of atrazine from the soil surface to the water table averaged 17.7 years and ranged from 2.2 to 118 years. Three maps were generated to present three views of the same atrazine vulnerability characteristics using different metrics (nonlinear transformations of the computed fractions remaining). The metrics were chosen because of the highly (right) skewed distribution of computed fractions. The first metric, rm = RM λ (where λ was 0.0625), depicts a relatively uniform distribution of vulnerability across the county with localized areas of high and low vulnerability visible. The second metric, rm λ-0.5 , depicts about one-half the county at low vulnerability with discontinuous patterns of high vulnerability evident. In the third metric, rm λ-1.0 ( RM ), more than 95 percent of the county appears to have low vulnerability; small, distinct areas of high vulnerability are present. Aquifer vulnerability estimates in the RM metric were used with a steady-state, uniform atrazine application rate to compute a potential concentration of atrazine in leachate reaching the water table. The average estimated potential atrazine concentration in leachate at the water table was 0.16 μg/L (micrograms per liter) in the model area; estimated potential concentrations ranged from 0 to 26 μg/L. About 2 percent of the model area had estimated potential atrazine concentrations in leachate at the water table that exceeded the USEPA (U.S. Environmental Protection Agency) maximum contaminant level of 3 μg/L. Uncertainty analyses were used to assess effects of parameter uncertainty and spatial interpolation error on the variability of the estimated fractions of atrazine remaining at the water table. Results of Monte Carlo simulations indicate that parameter uncertainty is associated with a standard error of 0.0875 in the computed fractions (in the rm metric). Results of kriging analysis indicate that errors in spatial interpolation are associated with a standard error of 0.146 (in the rm metric). Thus, uncertainty in fractions remaining is primarily associated with spatial interpolation error, which can be reduced by increasing the density of points where the leaching model is applied. A sensitivity analysis indicated which of 13 hydrologic, lithologic, and pesticide characteristics were influential in determining fractions of atrazine remaining at the water table. Results indicate that fractions remaining are most sensitive to the unit changes in pesticide half life and in organic-carbon content in soils and unweathered rocks, and least sensitive to infiltration rates. The leaching model applied in this report provides an estimate of the vulnerability of the near-surface aquifer in Kent County to contamination by atrazine. The vulnerability estimate is related to water-quality criteria developed by the USEPA to help assess potential risks from atrazine to the near-surface aquifer. However, atrazine accounts for only 28 percent of the herbicide use in the county; additional potential for contamination exists from other pesticides and pesticide metabolites. Therefore, additional work is needed to develop a comprehensive understanding of the relative risks associated with specific pesticides. The modeling approach described in this report provides a technique for estimating relative vulnerabilities to specific pesticides and for helping to assess potential risks.

Michigan↗

USGS2002

Rainfall-runoff characteristics and effects of increased urban density on streamflow and infiltration in the eastern part of the San Jacinto River basin, Riverside County, California

To better understand the rainfall-runoff characteristics of the eastern part of the San Jacinto River Basin and to estimate the effects of increased urbanization on streamflow, channel infiltration, and land-surface infiltration, a long-term (1950?98) time series of monthly flows in and out of the channels and land surfaces were simulated using the Hydrologic Simulation Program- FORTRAN (HSPF) rainfall-runoff model. Channel and land-surface infiltration includes rainfall or runoff that infiltrates past the zone of evapotranspiration and may become ground-water recharge. The study area encompasses about 256 square miles of the San Jacinto River drainage basin in Riverside County, California. Daily streamflow (for periods with available data between 1950 and 1998), and daily rainfall and evaporation (1950?98) data; monthly reservoir storage data (1961?98); and estimated mean annual reservoir inflow data (for 1974 conditions) were used to calibrate the rainfall-runoff model. Measured and simulated mean annual streamflows for the San Jacinto River near San Jacinto streamflow-gaging station (North-South Fork subbasin) for 1950?91 and 1997?98 were 14,000 and 14,200 acre-feet, respectively, a difference of 1.4 percent. The standard error of the mean for measured and simulated annual streamflow in the North-South Fork subbasin was 3,520 and 3,160 acre-feet, respectively. Measured and simulated mean annual streamflows for the Bautista Creek streamflow-gaging station (Bautista Creek subbasin) for 1950?98 were 980 acre-feet and 991 acre-feet, respectively, a difference of 1.1 percent. The standard error of the mean for measured and simulated annual streamflow in the Bautista Creek subbasin was 299 and 217 acre-feet, respectively. Measured and simulated annual streamflows for the San Jacinto River above State Street near San Jacinto streamflow-gaging station (Poppet subbasin) for 1998 were 23,400 and 23,500 acre-feet, respectively, a difference of 0.4 percent. The simulated mean annual streamflow for the State Street gaging station at the outlet of the study basin and the simulated mean annual basin infiltration (combined infiltration from all the channels and land surfaces) were 8,720 and 41,600 acre-feet, respectively, for water years 1950-98. Simulated annual streamflow at the State Street gaging station ranged from 16.8 acre-feet in water year 1961 to 70,400 acre-feet in water year 1993, and simulated basin infiltration ranged from 2,770 acre-feet in water year 1961 to 149,000 acre-feet in water year 1983.The effects of increased urbanization on the hydrology of the study basin were evaluated by increasing the size of the effective impervious and non-effective impervious urban areas simulated in the calibrated rainfall-runoff model by 50 and 100 percent, respectively. The rainfall-runoff model simulated a long-term time series of monthly flows in and out of the channels and land surfaces using daily rainfall and potential evaporation data for water years 1950?98. Increasing the effective impervious and non-effective impervious urban areas by 100 percent resulted in a 5-percent increase in simulated mean annual streamflow at the State Street gaging station, and a 2.2-percent increase in simulated basin infiltration. Results of a frequency analysis of the simulated annual streamflow at the State Street gaging station showed that when effective impervious and non-effective impervious areas were increased 100 percent, simulated annual streamflow increased about 100 percent for low-flow conditions and was unchanged for high-flow conditions. The simulated increase in streamflow at the State Street gaging station potentially could infiltrate along the stream channel further downstream, outside of the model area.

Water-Resources Investigations Report↗

USGS2005

Compilation of geologic, hydrologic, and ground-water flow modeling information for the Spokane Valley-Rathdrum Prairie aquifer, Spokane County, Washington, and Bonner and Kootenai Counties, Idaho

The U.S. Geological Survey, in cooperation with the Idaho Department of Water Resources and Washington Department of Ecology compiled and described geologic, hydrologic, and ground-water flow modeling information about the Spokane Valley-Rathdrum Prairie (SVRP) aquifer in northern Idaho and northeastern Washington. Descriptions of the hydrogeologic framework, water-budget components, ground- and surface-water interactions, computer flow models, and further data needs are provided. The SVRP aquifer, which covers about 370 square miles including the Rathdrum Prairie, Idaho and the Spokane valley and Hillyard Trough, Washington, was designated a Sole Source Aquifer by the U.S. Environmental Protection Agency in 1978. Continued growth, water management issues, and potential effects on water availability and water quality in the aquifer and in the Spokane and Little Spokane Rivers have illustrated the need to better understand and manage the region's water resources. The SVRP aquifer is composed of sand, gravel, cobbles, and boulders primarily deposited by a series of catastrophic glacial outburst floods from ancient Glacial Lake Missoula. The material deposited in this high-energy environment is coarser-grained than is typical for most basin-fill deposits, resulting in an unusually productive aquifer with well yields as high as 40,000 gallons per minute. In most places, the aquifer is bounded laterally by bedrock composed of granite, metasedimentary rocks, or basalt. The lower boundary of the aquifer is largely unknown except along the margins or in shallower parts of the aquifer where wells have penetrated its entire thickness and reached bedrock or silt and clay deposits. Based on surface geophysics, the thickness of the aquifer is about 500 ft near the Washington-Idaho state line, but more than 600 feet within the Rathdrum Prairie and more than 700 feet in the Hillyard trough based on drilling records. Depth to water in the aquifer is greatest in the northern Rathdrum Prairie (about 500 feet) and least near the city of Spokane along the Spokane River (less than about 50 feet). Ground-water flow is south from near the southern end of Lake Pend Oreille and Hoodoo Valley, through the Rathdrum Prairie, then west toward Spokane. In Spokane, the aquifer splits and water moves north through the Hillyard Trough as well as west through the Trinity Trough. From the Trinity Trough water flows north along the western arm of the aquifer. The aquifer's discharge area is along the Little Spokane River and near Long Lake, Washington. A compilation of estimates of water-budget components, including recharge (precipitation, irrigation, canal leakage, septic tank effluent, inflow from tributary basins, and flow from the Spokane River) and discharge (withdrawals from wells, flow to the Spokane and Little Spokane Rivers, evapotranspiration, and underflow to Long Lake) illustrates that these estimated values should be compared with caution due to several variables including the area and time period of interest as well as methods employed in making the estimates. Numerous studies have documented the dynamic ground-water and surface-water interaction between the SVRP aquifer and the Spokane and Little Spokane Rivers. Gains and losses vary throughout the year, as well as the locations of gains and losses. September 2004 streamflow measurements indicated that the upper reach of the Spokane River between Post Falls and downstream at Flora Road lost 321 cubic feet per second. A gain of 736 cubic feet per second was measured between the Flora Road site and downstream at Green Street Bridge. A loss of 124 cubic feet per second was measured for the reach between the Green Street Bridge and the Spokane River at Spokane gaging station. The river gained about 87 cubic feet per second between the Spokane River at Spokane gaging station and the TJ Meenach Bridge. Overall, the Spokane River gained about 284 cubic feet per second between the Post Falls,

Scientific Investigations Report↗

USGS2009

Geologic Map of MTM 35337, 40337, and 45337 Quadrangles, Deuteronilus Mensae Region of Mars

Deuteronilus Mensae, first defined as an albedo feature at lat 35.0 deg N., long 5.0 deg E., by U.S. Geological Survey (USGS) and International Astronomical Union (IAU) nomenclature, is a gradational zone along the dichotomy boundary in the northern mid-latitudes of Mars. The boundary in this location includes the transition from the rugged cratered highlands of Arabia Terra to the northern lowland plains of Acidalia Planitia. Within Deuteronilus Mensae, polygonal mesas are prominent along with features diagnostic of Martian fretted terrain, including lobate debris aprons, lineated valley fill, and concentric crater fill. Lobate debris aprons, as well as the valley and crater fill deposits, are geomorphic indicators of ground ice, and their concentration in Deuteronilus Mensae is of great interest because of their potential association with Martian climate change. The paucity of impact craters on the surfaces of debris aprons and the presence of ice-cemented mantle material imply young (for example, Amazonian) surface ages that are consistent with recent climate change in this region of Mars. North of Deuteronilus Mensae are the northern lowlands, a potential depositional sink that may have had large standing bodies of water or an ocean in the past. The northern lowlands have elevations that are several kilometers below the ancient cratered highlands with significantly younger surface ages. The morphologic and topographic characteristics of the Deuteronilus Mensae region record a diverse geologic history, including significant modification of the ancient highland plateau and resurfacing of low-lying regions. Previous studies of this region have interpreted a complex array of geologic processes, including eolian, fluvial and glacial activity, coastal erosion, marine deposition, mass wasting, tectonic faulting, effusive volcanism, and hydrovolcanism. The origin and age of the Martian crustal dichotomy boundary are fundamental questions that remain unresolved at the present time. Several scenarios for its formation, including single and multiple large impact events, have been proposed and debated in the literature. Endogenic processes whereby crust is thinned by internal mantle convection and tectonic processes have also been proposed. Planetary accretion models and isotopic data from Martian meteorites suggest that the crust formed very early in Martian history. Using populations of quasi-circular depressions extracted from the topography of Mars, other studies suggest that the age difference between the highlands and lowlands could be ~100 m.y.. Furthermore, understanding the origin and age of the dichotomy boundary has been made more complicated due to significant erosion and deposition that have modified the boundary and its adjacent regions. The resulting diversity of terrains and features is likely a combined result of ancient and recent events. Detailed geologic analyses of dichotomy boundary zones are important for understanding the spatial and temporal variations in highland evolution. This information, and comparisons to other highland regions, can help elucidate the scale of potential environmental changes. Previous geomorphic and geologic mapping investigations of the Deuteronilus Mensae region have been completed at local to global scales. The regional geology was first mapped by Lucchitta (1978) at 1:5,000,000 scale using Mariner 9 data. This study concluded that high crater flux early in Martian history formed overlapping craters and basins that were later filled by voluminous lava flows that buried the impacted surface, creating the highlands. After this period of heavy bombardment, fluvial erosion of the highlands formed the canyons and valleys, followed by dissection that created the small mesas and buttes, and later, formation of the steep escarpment marking the present-day northern highland margin. After valley dissection, mass wasting and eolian processes caused lateral retreat of mesas and buttes

Scientific Investigations Map↗

USGS1996

Geohydrology and conceptual model of a ground-water-flow system near a Superfund site in Cheshire, Connecticut

Degradation of ground-water quality has been identified in an area of the north-central part of the town of Cheshire, Connecticut. An investigation by the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, was done during 1994-95 to characterize the unconsolidated glacial deposits and the sedimentary bedrock, integrate the local geohydrologic conditions with the regional geohydrologic system, and develop a conceptual understanding of ground-water flow in the study area. A regional ground-water-flow model developed for the region near the study area indicates that perennial streams, including Judd Brook and the Tenmile River, form hydrologic divides that separate the larger region into hydraulically independent flow systems. In the local study area, synoptic water-level measurements made in June 1995 indicate that ground water near the water table flows west and southwestward from the low hill on the eastern side of the area toward the pond and wetlands along Judd Brook. Water-level data indicate that there is good hydraulic connection between the unconsolidated materials and underlying fractured bedrock. Unconsolidated materials in the study area consist principally of glacial stratified deposits that are fine sand, silt, and clay of glaci- olacustrine origin; locally these overlie thin glacial till. The glacial sediments range in thickness from a few feet to about 25 ft in the eastern part of the study area and are as much as 100 ft thick in the western and southern part of the study area beneath the Judd Brook and Tenmile River valleys. Fluvial redbeds of the New Haven Arkose underlie the glacial deposits in the region; in the study area, the redbeds consist of (1) channel sandstone units, which are coarse sandstone to fine conglomerate, generally in 6- to 15-ft- thick sequences; and (2) overbank mudstone units, which are siltstone and silty sandstone with some fine sandstone, generally in 6- to 50-ft-thick sequences. Thin-bedded zones of siltstone that are particularly fissile are present locally within the mudstone units. Rock units strike northward and dip eastward at about 20. The eastward-dipping strata are cut by a consistent set of west to west-northwest dipping, high-angle fractures. These fractures are oriented perpendicular to bedding and are present mostly in the channel sandstone units, but locally extend into the mudstone units as well. Borehole-geophysical logging indicates that ground water flows along bedding planes in fissile zones and between fissile zones in high-angle fractures, which are perpendicular to bedding. The combined fracture types form an aquifer system in which ground water follows a stair-step flowpath, flowing horizontally through fissile zones and vertically through high-angle fractures. Heat-pulse flow meter measurements and borehole fluid-conductivity and temperature logs indicate that only a small subset of the fissile zones and some high-angle fractures are hydraulically significant. A generalized local-scale ground-water flow model based on a nonspecific, but realistic, rock and fracture geometry was developed for the study area. Simulations show that under nonpumping conditions at a hypothetical well located in the middle of the model, ground-water flow was separated into upper and lower zones in which flow paths differed but were generally from northeast to southwest. Several short-duration aquifer tests conducted in the study area indicate that there is good hydraulic connection in the fractures between the pumping well (CS-221) and two bedrock wells located approximately 100 ft to the north and south along bedding strike. During the short duration of the aquifer tests, there was no hydraulic connection in bedrock wells located to the east, perpendicular to the strike. A range of transmissivity of 27 to 46 ft2/d was calculated from the aquifer-test data for the fractured-bedrock aquifer at CS-221 and TH-2. Individual fracture zones identified by bo

Open-File Report↗

USGS2005

Simulated water sources and effects of pumping on surface and ground water, Sagamore and Monomoy flow lenses, Cape Cod, Massachusetts

The sandy sediments underlying Cape Cod, Massachusetts, compose an important aquifer that is the sole source of water for a region undergoing rapid development. Population increases and urbanization on Cape Cod lead to two primary environmental effects that relate directly to water supply: (1) adverse effects of land use on the quality of water in the aquifer and (2) increases in pumping that can adversely affect environmentally sensitive surface waters, such as ponds and streams. These considerations are particularly important on the Sagamore and Monomoy flow lenses, which underlie the largest and most populous areas on Cape Cod. Numerical models of the two flow lenses were developed to simulate ground-water-flow conditions in the aquifer and to (1) delineate areas at the water table contributing water to wells and (2) estimate the effects of pumping and natural changes in recharge on surface waters. About 350 million gallons per day (Mgal/d) of water recharges the aquifer at the water table in this area; most water (about 65 percent) discharges at the coast and most of the remaining water (about 28 percent) discharges into streams. A total of about 24.9 Mgal/d, or about 7 percent, of water in the aquifer is withdrawn for water supply; most pumped water is returned to the hydrologic system as return flow creating a state of near mass balance in the aquifer. Areas at the water table that contribute water directly to production wells total about 17 square miles; some water (about 10 percent) pumped from the wells flows through ponds prior to reaching the wells. Current (2003) steady-state pumping reduces simulated ground-water levels in some areas by more than 4 feet; projected (2020) pumping may reduce water levels by an additional 3 feet or more in these same areas. Current (2003) and future (2020) pumping reduces total streamflow by about 4 and 9 cubic feet per second (ft3/s), corresponding to about 5 percent and 9 percent, respectively, of total streamflow. Natural recharge varies with time, over both monthly and multiyear time scales. Monthly changes in recharge cause pond levels to vary between 1 and 2 feet in an average year; annual changes in recharge, which can be much larger than monthly variations, can cause pond levels to vary by more than 10 feet in some areas over a period of years. Streamflow, which also changes in response to changes in recharge, varies by a factor of two over an average year and can vary more over multiyear periods. On average, monthly pumping ranges from 15.8 Mgal/d in March to 45.3 Mgal/d in August. Pumping and the distribution of return flow can seasonally affect the hydrologic system by lowering ground-water and pond levels and by depleting streamflows, particularly in the summer months. Maximum drawdowns in March and August exceed 3 feet and 6 feet, respectively, for current (2003) pumping. Simulated drawdowns from projected (2020) pumping, relative to water levels representing 2003 pumping conditions, exceed 2 feet in March and 5 feet in August. Current (2003) and future (2020) pumping can decrease pond levels in some areas by more than 3 feet; drawdown generally is largest during the month of August of an average year. Over multiyear periods, seasonal pumping can lower pond levels in some areas by more than 4 feet; the effects of seasonal pumping are largest during periods of reduced recharge. Monthly streamflow depletion varies in individual streams but can exceed 2 ft3/s in some streams. The combined effects of seasonal pumping and drought can reduce pond levels by more than 10 feet below average levels. Water levels in Mary Dunn Pond, which is in an area of large current and projected pumping, are predicted (2020) to decline during drought conditions by about 10.6 feet: about 6.9 feet from lower recharge, about 2.3 feet from current (2003) pumping, and about 1.4 feet from additional future (2020) pumping. The results indicate that pumping generally does not cause substantial streamflow depletion and that the primary effect of pumping is on water levels in ponds. Natural changes in recharge account for most of the variation in pond levels; however, pumping can cause substantial declines in the levels of ponds near pumping wells. Also, the effects of pumping and recharge can combine to cause drawdowns of more than 10 feet in some areas.

Massachusetts↗

USGS1996

Hydraulic characteristics and nutrient transport and transformation beneath a rapid infiltration basin, Reedy Creek Improvement District, Orange County, Florida

The Reedy Creek Improvement District disposes of about 7.5 million gallons per day (1992) of reclaimed water through 85 1-acre rapid infiltration basins within a 1,000-acre area of sandy soils in Orange County, Florida. The U.S. Geological Survey conducted field experiments in 1992 at an individual basin to examine and better understand the hydraulic characteristics and nutrient transport and transformation of reclaimed water beneath a rapid infiltration basin. At the time, concentrations of total nitrogen and total phosphorus in reclaimed water were about 3 and 0.25 milligrams per liter, respectively. A two-dimensional, radial, unsaturated/saturated numerical flow model was applied to describe the flow system beneath a rapid infiltration basin under current and hypothetical basin loading scenarios and to estimate the hydraulic properties of the soil and sediment beneath a basin. The thicknesses of the unsaturated and saturated parts of the surficial aquifer system at the basin investigated were about 37 and 52 feet, respectively. The model successfully replicated the field-monitored infiltration rate (about 5.5 feet per day during the daily flooding periods of about 17 hours) and ground-water mounding response during basin operation. Horizontal and vertical hydraulic conductivity of the saturated part of the surficial aquifer system were estimated to be 150 and 45 feet per day, respectively. The field-saturated vertical hydraulic conductivity of the shallow soil, estimated to be about 5.1 feet per day, was considered to have been less than the full-saturation value because of the effects of air entrapment. Specific yield of the surficial aquifer was estimated to be 0.41. The upper 20 feet of the basin subsurface profile probably served as a system control on infiltration because of the relatively low field-saturated, vertical hydraulic conductivity of the sediments within this layer. The flow model indicates that, in the vicinity of the basin, flow in the deeper, saturated zone was relatively slow compared to the more vigorous flow in the shallow saturated zone. The large radial component of flow below the water table in the vicinity of the basin implies that reclaimed water moves preferentially in the shallow part of the saturated zone upon reaching the water table. Therefore, there may be some vertical stratification in the saturated zone, with recently infiltrated water overlying ambient water. The infiltration capacity at the basin would be unaffected by a small (less than 10 feet) increase in background water-table altitude, because the water table would remain below the system control on infiltration. However, water-table rises of 15 and 20 feet were estimated to reduce the infiltration capacity of the basin by 8 and 25 percent, respectively. Model simulations indicate that increasing ponded depth within the basin from 4 to 12 inches and from 4 to 24 inches would increase basin infiltration capacity by less than 6 and 11 percent, respectively. A loading strategy at the basin that relies on long, uninterrupted flooding was shown to offer the possibility of inducing a more anaerobic environment conducive to denitrification while maintaining reclaimed-water disposal capacity. Field measurements indicated that transient, elevated concentrations or "spikes" of nitrate (as high as 33 milligrams per liter as nitrogen) occurred at the leading edge of the infiltrating water and in the shallow saturated zone following a prolonged basin rest period. This phenomenon probably is the result of mineralization and nitrification of organic nitrogen retained within the subsurface during earlier basin loading events. The organic nitrogen was retained in the shallow soil (due to adsorption/straining) and the shallow saturated zone (due to deposition under slacking pore-water velocity). The magnitude of the nitrate spikes appears to be influenced by the scheduling of basin loading, with short flooding and resting periods being most favorable to minimization of nitrate spikes. Removal of nitrogen by denitrification from the percolating reclaimed water is minimal in the vicinity of the basin, probably because of the lack of reducing conditions and a relative paucity of organic carbon substrates. It is speculated that longer flooding periods could induce reducing conditions that would be favorable for nitrogen removal from the system, but probably would lead to more pronounced nitrate spiking. Phosphorus concentrations were decreased by about 90 percent from concentrations in reclaimed water after moving through the upper 15 feet of the soil profile. This most likely was a result of adsorption onto abundant iron and aluminum hydrous oxyhydroxide coatings on sand grains. However, some phosphorus (coarse fraction, organic) passes through the shallow soil and accumulates below the water table under slacking pore-water velocity. This phosphorus is immobilized by adsorption, fixation, or precipitation reactions during basin rest periods.

Florida↗

USGS2020

Storage capacity and sedimentation characteristics of the San Antonio Reservoir, California, 2018

The San Antonio Reservoir is a large water storage facility in Alameda County, California, and is a major component of the Hetch Hetchy Regional Water System (RWS). The RWS is a water-supply system owned and operated by the San Francisco Public Utilities Commission (SFPUC) and provides water for about 2.7 million people in the San Francisco, Santa Clara, Alameda, and San Mateo Counties. The San Antonio Reservoir is one of two RWS reservoirs in Alameda County and the third largest of the RWS reservoirs in the San Francisco Bay Area. The reservoir was formed by the James H. Turner Dam, which was completed in 1965. At the time of construction, the reservoir was estimated to have 50,500 acre-feet (acre-ft) of storage capacity. That early estimate was based on a 1963 pre-construction topographic map, which was drawn from aerial photographs. The capacity of the reservoir was later surveyed in 1994 and 2000. These two later surveys did not include the upper 18 feet (ft) of the reservoir, which represents roughly 30 percent of the overall storage volume. To determine the storage capacity and provide updated stage-capacity curves up to the spillway, the U.S. Geological Survey, in cooperation with the SFPUC, surveyed the bathymetry and shoreline of the reservoir in April 2018. The bathymetric survey was performed by making depth soundings using a boat-mounted, multibeam echosounder. At the time of the survey, the water level was between 13 and 14 ft below the spillway elevation. To measure capacity between the water line up to the spillway elevation, topography along most of the shoreline was surveyed from the boat using a terrestrial Light Detection and Ranging (LiDAR) scanner and in other areas by using ground-survey techniques. Location during bathymetric and topographic data collection was determined using a Global Navigation Satellite System-Real Time Network system. Vertical profiles of sound speed were collected periodically. The sound-speed profiles were used to spatially and temporally adjust the sound-speed calculations used to determine depth from the soundings. Approximately 125 kilometers (78 miles) of transects with a total of about 560 million depth soundings and topographic LiDAR points were collected (about 160 per square meter). In addition, approximately 500 topographic survey points were collected in shallow, wadable areas and on land near the upper reservoir area using a Global Navigation Satellite System receiver attached to a fixed length survey rod. Depth soundings, terrestrial LiDAR points, topographic survey points, and a digitized shoreline were merged and interpolated to generate a digital elevation model (DEM) of the reservoir. Gridded elevation data extracted from the DEM were then tabulated to determine total reservoir capacity and create reservoir stage-surface area and stage-storage capacity tables. Results of the reservoir capacity analysis indicated that the reservoir has 53,266 (plus or minus 140) acre-ft of storage capacity, which is an increase of 2,766 acre-ft (or 5.5 percent) greater than the original 1965 estimate; the increase is likely due to improved survey methods. Also, at the time of this 2018 survey, Intake #1 (the lowest intake) was not in operation. Intake #1 is estimated to be buried approximately 10 ft below the bed, whereas Intake #2 is about 20 ft above the bed. There are five intakes at different elevation levels; however, when consecutive lower intakes become inoperable due to sedimentation, the live storage capacity (capacity available for use) is reduced. At the time of this survey, the remaining live storage (above Intake #2) was approximately 52,363 acre-ft. The 2018 stage-capacity curve was compared to the original 1965 stage-capacity curve. Although overall, the changes indicate an increase in storage capacity, the change in volume at 372.7 ft North American Vertical Datum of 1988 (370 ft National Geodetic Vertical Datum of 1929, NGVD 29) shows a decrease of 733 acre-ft (the elevation of 370 ft NGVD 29 was used because it is the lowest elevation available for the 1965 stage-capacity curves). This finding agrees with the observed accumulation of sediment over Intake #1. That volume was converted to an annual sediment yield of 0.35 acre-ft per square mile (or 165 cubic meters per square kilometer), which is of the same order of magnitude as that found in other watersheds for the Coast Ranges in California. A decrease of 733 acre-ft between 1965 and 2018 thus represents a loss of 1.5 percent of the overall storage capacity in the reservoir. The updated stage-surface area and stage-capacity tables provided in this report and online ( https://doi.org/10.5066/P9KC9DU8 ) can be used by the SFPUC to improve reservoir operations and serve as an accurate baseline to monitor bathymetric changes in the future.

California↗

USGS2002

Ground-water, surface-water, and water-chemistry data, Black Mesa Area, northeastern Arizona — 2000–2001, and performance and sensitivity of the 1988 USGS numerical model of the N aquifer

The N aquifer is the major source of water in the 5,400-square-mile area of Black Mesa in northeastern Arizona. Availability of water is an important issue in this area because of continued industrial and municipal use, a growing population, and precipitation of about 6 to 14 inches per year. The monitoring program in Black Mesa has been operating since 1971 and is designed to determine the long-term effects of ground-water withdrawals from the N aquifer for industrial and municipal uses. The monitoring program includes measurements of (1) ground-water pumping, (2) ground-water levels, (3) spring discharge, (4) surface-water discharge, and (5) ground-water chemistry. In 2000, total ground-water withdrawals were 7,740 acre-feet, industrial use was 4,490 acre-feet, and municipal use was 3,250 acre-feet. From 1999 to 2000, total withdrawals increased by 9 percent, industrial use increased by 7 percent, and municipal use increased by 12 percent. From 1999 to 2001, water levels declined in 10 of 15 wells in the unconfined part of the aquifer, and the median change was -0.4 foot. Water levels declined in 8 of 16 wells in the confined part of the aquifer, and the median change was -0.2 foot. From the prestress period (prior to 1965) to 2001, the median water-level change for 33 wells was -17.2 feet. Median water-level changes were -1.2 feet for 15 wells in the unconfined part of the aquifer and -31.0 feet for 18 wells in the confined part. Discharges were measured once in 1999 and once in 2001 at four springs. Discharges decreased by 5 percent and 33 percent at two springs and increased by 3 percent and 81 percent at two springs. For about the past 10 years, discharges did not significantly change in Burro Spring, the unnamed spring near Dennehotso, and Moenkopi School Spring. The record of discharge from a consistent measuring point for Pasture Canyon Spring is too short for statistical analysis of trends. Continuous records of surface-water discharge have been collected from July 1976 to 2000 at Moenkopi Wash, July 1996 to 2000 at Laguna Creek, June 1993 to 2000 at Dinnebito Wash, and April 1994 to 2000 at Polacca Wash. Median flows for November, December, January, and February of each water year were used as an index of ground-water discharge to those streams. There is no significant trend in the median winter flows for Moenkopi Wash from 1977 to 2000. The records for the other three streams are too short for a statistical analysis of trends. The median winter flows for Dinnebito Wash and Polacca Wash, however, appear to have decreased during the last 6 years. There is no apparent trend in the median winter flows for Laguna Creek since 1997. In 2001, water samples were collected from 12 wells and 4 springs and analyzed for selected chemical constituents. Dissolved-solids concentrations ranged from 102 to 628 milligrams per liter. Water samples from 9 of the wells and from the 4 springs had less than 350 milligrams per liter of dissolved solids. Water-chemistry data with sufficient years of record for a statistical analysis of trends over time are available from 7 wells and 4 springs. From about the mid-1980s or early 1990s to 2001 there are no significant trends in the concentrations of dissolved solids, chloride, and sulfate in water samples from 6 of the 7 wells. The concentration of one tested constituent (dissolved solids) in samples from Rocky Ridge PM3 significantly increased from 1990 to 2001. From the late 1980s to 2001, there are no significant trends in the concentrations of dissolved solids, chloride, and sulfate in water samples from Burro Spring, the unnamed spring near Dennehotso, and Pasture Canyon Spring. From 1987 to 2001, concentrations of chloride and sulfate significantly increased in water samples from Moenkopi School Spring and concentrations of dissolved solids did not significantly change. The performance and sensitivity of the 1988 USGS numerical model of the N aquifer were analyzed. The overall performance of the model in steady-state conditions is reasonable for residuals of heads (difference between observed and simulated steady-state heads); 80 percent of the absolute values of residuals are less than 38 feet. Simulated flows are about 40 percent different than estimated flows at two of three discharge areas; however, this comparison is only a rough approximation of performance because the accuracy of the estimated steady-state flows is uncertain. The overall performance of the model for transient conditions is fair for residuals of changes in head (difference between observed and simulated changes in head from steady state to 1999); 80 percent of the absolute values of residuals are less than 31 feet. The model is biased in two areas. In the Tuba City area, simulated changes in head are more negative than observed changes in head; all six residuals are positive, and three residuals are between 75 and 155 feet. In the confined area of the aquifer, observed changes in head are more negative than simulated changes in head; 12 of the 17 residuals are negative, and 8 residuals are between -57 and -20 feet. Analysis of model sensitivity indicates that recharge, transmissivity, and storage coefficient are the most important parameters for estimating heads, changes in heads, and flows. A strong correlation between recharge and transmissivity and a lack of independent and reliable estimates of recharge, transmissivity, and discharge create a uniqueness problem in model calibration. Several models could be constructed and calibrated with different values of recharge or transmissivity and still have similar fits to the observed data. Information from recent data and studies and more advanced modeling techniques could be used to develop a more representative and less uncertain model. Future data collection and studies should focus on obtaining a better definition of recharge, discharge, transmissivity, and storage coefficient.

Arizona↗