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 1,153 records · Page 64Linked to original sources
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↗
USGS2002

Hydrogeology and simulated effects of ground-water withdrawals from the Floridan aquifer system in Lake County and in the Ocala National Forest and vicinity, north-central Florida

The hydrogeology of Lake County and the Ocala National Forest in north-central Florida was evaluated (1995-2000), and a ground-water flow model was developed and calibrated to simulate the effects of both present day and future ground-water withdrawals in these areas and the surrounding vicinity. A predictive model simulation was performed to determine the effects of projected 2020 ground-water withdrawals on the water levels and flows in the surficial and Floridan aquifer systems. The principal water-bearing units in Lake County and the Ocala National Forest are the surficial and Floridan aquifer systems. The two aquifer systems generally are separated by the intermediate confining unit, which contains beds of lower permeability sediments that confine the water in the Florida aquifer system. The Floridan aquifer system has two major water-bearing zones (the Upper Floridan aquifer and the Lower Floridan aquifer), which generally are separated by one or two less-permeable confining units. The Floridan aquifer system is the major source of ground water in the study area. In 1998, ground-water withdrawals totaled about 115 million gallons per day in Lake County and 5.7 million gallons per day in the Ocala National Forest. Of the total ground water pumped in Lake County in 1998, nearly 50 percent was used for agricultural purposes, more than 40 percent for municipal, domestic, and recreation supplies, and less than 10 percent for commercial and industrial purposes. Fluctuations of lake stages, surficial and Floridan aquifer system water levels, and Upper Floridan aquifer springflows in the study area are highly related to cycles and distribution of rainfall. Long-term hydrographs for 9 lakes, 8 surficial aquifer system and Upper Floridan aquifer wells, and 23 Upper Floridan aquifer springs show the most significant increases in water levels and springflows following consecutive years with above-average rainfall, and significant decreases following consecutive years with below-average rainfall. Long-term (1940-2000) hydrographs of lake and ground-water levels and springflow show a slight downward trend; however, after the early 1960's, this downward trend generally is more pronounced, which corresponds with accumulating rainfall deficits and increased development. The U.S. Geological Survey three-dimensional ground-water flow model MODFLOW-2000 was used to simulate ground-water flow in the surficial and Floridan aquifer systems in Lake County, the Ocala National Forest, and adjacent areas. A steady-state calibration to average 1998 conditions was facilitated by using the inverse modeling capabilities of MODFLOW-2000. Values of hydrologic properties from the calibrated model were in reasonably close agreement with independently estimated values and results from previous modeling studies. The calibrated model generally produced simulated water levels and flows in reasonably close agreement with measured values and was used to simulate the hydrologic effects of projected 2020 conditions. Ground-water withdrawals in the model area have been projected to increase from 470 million gallons per day in 1998 to 704 million gallons per day in 2020. Significant drawdowns were simulated in Lake County from average 1998 to projected 2020 conditions: the average and maximum drawdowns, respectively, were 0.5 and 5.7 feet in the surficial aquifer system, 1.1 and 7.6 feet in the Upper Floridan aquifer, and 1.4 and 4.3 feet in the Lower Floridan aquifer. The largest drawdowns in Lake County were simulated in the southeastern corner of the County and in the vicinities of Clermont and Mount Dora. Closed-basin lakes and wetlands are more likely to be affected by future pumping in these large drawdown areas, as opposed to other areas of Lake County. However, within the Ocala National Forest, drawdowns were relatively small: the average and maximum drawdowns, respectively, were 0.1 and 1.0 feet in the surficial aquifer system, 0.2 and 0.8 feet in the Upper Floridan aquifer, and 0.3 and 0.8 feet in the Lower Floridan aquifer. Projected 2020 withdrawals from the Floridan aquifer system caused decreases from average 1998 conditions in the following simulated flows: combined rates of excess evapotranspiration and excess overland runoff (which represent evapotranspiration and overland runoff that occur in excess of their assumed minimum rates); ground-water discharge to streams, lakes, and wetlands; and springflow. The largest simulated flow decreases for first- or second-magnitude springs in Lake County were at Apopka (28 percent), Seminole (12 percent), and Bugg Springs (9 percent). The largest simulated flow decrease for first- or second-magnitude springs in the Ocala National Forest was at Juniper Springs (4 percent). Particle-tracking analyses were used to delineate areas that contribute recharge to selected springs. Based on average 1998 conditions, the contributing area for Apopka Spring covers approximately 30 square miles and has an average contributing recharge flux of 15 inches per year, and the contributing area for Alexander Springs covers approximately 76 square miles and has an average contributing recharge flux of 18 inches per year. The contributing area for Alexander Springs changed little as a result of projected 2020 conditions because relatively little pumping exists in the vicinity of the spring's contributing area. However, the size of the contributing area for Apopka Spring decreased to 26 square miles and the average contributing recharge flux decreased to 13 inches per year as a result of projected 2020 conditions.

Florida↗
USGS1998

Ground-water discharge and base-flow nitrate loads of nontidal streams, and their relation to a hydrogeomorphic classification of the Chesapeake Bay watershed, middle Atlantic Coast

Existing data on base-flow and groundwater nitrate loads were compiled and analyzed to assess the significance of groundwater discharge as a source of the nitrate load to nontidal streams of the Chesapeake Bay watershed. These estimates were then related to hydrogeomorphic settings based on lithology and physiographic province to provide insight on the areal distribution of ground-water discharge. Base-flow nitrate load accounted for 26 to about 100 percent of total-flow nitrate load, with a median value of 56 percent, and it accounted for 17 to 80 percent of total-flow total-nitrogen load, with a median value of 48 percent. Hydrograph separations were conducted on continuous streamflow records from 276 gaging stations within the watershed. The values for base flow thus calculated were considered an estimate of ground-water discharge. The ratio of base flow to total flow provided an estimate of the relative importance of ground-water discharge within a basin. Base-flow nitrate loads, total-flow nitrate loads, and total-flow total-nitrogen loads were previously computed from water-quality and discharge measurements by use of a regression model. Base-flow nitrate loads were available from 78 stations, total-flow nitrate loads were available from 86 stations, and total-flow total-nitrogen loads were available for 48 stations. The percentage of base-flow nitrate load to total-flow nitrate load could be computed for 57 stations, whereas the percentage of base-flow nitrate load to totalflow total-nitrogen load could be computed for 36 stations. These loads were divided by the basin area to obtain yields, which were used to compare the nitrate discharge from basins of different sizes. The results indicate that ground-water discharge is a significant source of water and nitrate to the total streamflow and nitrate load. Base flow accounted for 16 to 92 percent of total streamflow at the 276 sampling sites, with a median value of 54 percent. It is estimated that of the 50 billion gallons of water that reaches the Chesapeake Bay each day, nearly 27 billion gallons is base flow. Generalized lithology (siliciclastic, carbonate, crystalline, and unconsolidated) was combined with physiographic province (the Appalachian Plateau, the Valley and Ridge, the Blue Ridge, the Piedmont, including the Mesozoic Lowland section, and the Coastal Plain) to delineate 11 hydrogeomorphic regions. Areal variation of base flow and base-flow nitrate yield were assessed by means of nonparametric, one-way analysis of variance on basins grouped by the dominant hydrogeomorphic region and by correlation analysis of base flow or base-flow nitrate yield with the percentage of land area of a given hydrogeomorphic region within a basin. Base flow appeared to have a significant relation to the hydrogeomorphic regions. The highest percentages of base flow were found in areas underlain by carbonate rock, crystalline rock with relatively low relief, and unconsolidated sediments. Lower percentages were found in areas underlain by siliclastic rocks and crystalline rocks with relatively high relief. The relation between base-flow nitrate yield and hydrogeomorphic region is less clear. Although there is a relation between low nitrate yields and areas underlain by highrelief siliciclastic rocks, and a relation between high yields and carbonate rocks, much of this relation can be explained by the strong association between the hydrogeomorphic units and land use. In addition, most basins are mixtures of several hydrogeomorphic regions, so the nitrate yield from a basin depends on a large number of complex interacting factors. These unclear results indicate that the sample of available data used here may not be adequate to fully assess the relation between base-flow nitrate yield and the hydrogeomorphic setting of the basin. The results appear to show, however, that ground-water discharge is an important component of the total nontidal streamflow, and that ground-water discharge varies according to the hydrogeomorphic regions. Environmental management of the nontidal streams in the Chesapeake Bay watershed will thus have to consider the prevention of nutrient infiltration into aquifers as well as prevention of overland runoff of high-nitrogen waters.

Chesapeake Bay watershed↗
USGS2004

Geochemical characterization of ground-water flow in the Santa Fe Group aquifer system, Middle Rio Grande Basin, New Mexico

Chemical and isotopic data were obtained from ground water and surface water throughout the Middle Rio Grande Basin (MRGB), New Mexico, and supplemented with selected data from the U.S. Geological Survey (USGS) National Water Information System (NWIS) and City of Albuquerque water-quality database in an effort to refine the conceptual model of ground-water flow in the basin. The ground-water data collected as part of this study include major- and minor-element chemistry (30 elements), oxygen-18 and deuterium content of water, carbon-13 content and carbon-14 activity of dissolved inorganic carbon, sulfur-34 content of dissolved sulfate, tritium, and dissolved atmospheric gases including nitrogen, argon, helium, chlorofluorocarbons, and sulfur hexafluoride from 288 wells and springs in parts of the Santa Fe Group aquifer system. The surface-water data collected as part of this study include monthly measurements of major- and minor-element chemistry (30 elements), oxygen-18 and deuterium content of water, chlorofluorocarbons, and tritium content at 14 locations throughout the basin. Additional data include stable isotope analyses of precipitation and of ground water from City of Albuquerque production wells collected and archived from the early 1980?s, and other data on the chemical and isotopic composition of air, unsaturated zone air, plants, and carbonate minerals from throughout the basin. The data were used to identify 12 sources of water to the basin, map spatial and vertical extents of ground-water flow, map water chemistry in relation to hydrogeologic, stratigraphic, and structural properties of the basin, determine radiocarbon ages of ground water, and reconstruct paleo-environmental conditions in the basin over the past 30,000 years. The data indicate that concentrations of most elements and isotopes generally parallel the predominant north to south direction of ground-water flow. The radiocarbon ages of dissolved inorganic carbon in ground water range from modern (post-1950) to more than 30,000 years before present, and appear to be particularly well defined in the predominantly siliciclastic aquifer system. Major sources of water to the basin include (1) recharge from mountains along the north, east and southwest margins (median age 5,000-9,000 years); (2) seepage from the Rio Grande and Rio Puerco (median age 4,000-8,000 years), and from Abo and Tijeras Arroyos (median age 3,000-9,000 years); (3) inflow of saline water along the southwestern basin margin (median age 20,000 years); and (4) inflow along the northern basin margin that probably represents recharge from the Jemez Mountains during the last glacial period (median age 20,000 years). Water recharged from the Jemez Mountains during the last glacial period occurs at the water table in the central part of the basin and beneath younger recharge along the Rio Grande and the northern mountain front. In some parts of the basin, boundaries between hydrochemical zones appear to be near major faults that may affect ground-water flow. However, in other parts of the basin, such as along the east side of Albuquerque near the Sandia Fault zone, ground-water flow appears to be unaffected by major faults. Upward leakage of saline water occurs along some faults and can be a source of salinity and elevated arsenic concentrations in some ground water. A trough in the modern and predevelopment water table west of Albuquerque is centered along a zone of predominantly late Pleistocene age water through the center of the basin and is flanked and overlain along the trough boundary by water that infiltrated from the Rio Puerco on the west and the Rio Grande to the east. It is suggested that the groundwater trough is a relatively recent transient feature of the Santa Fe Group aquifer system. At Albuquerque, a distinct north-south boundary in deuterium content of ground water marks the division between recharge from the eastern mountain front and that from the Rio Grande.

New Mexico↗
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 pumping ceases and would reco ver more quickly with continued use of the infiltration basin for recharge.

Nevada↗
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.

California↗
USGS2002

Ground-water flow simulation and chemical and isotopic mixing equation analysis to determine source contributions to the Missouri River alluvial aquifer in the vicinity of the Independence, Missouri, well field

The city of Independence, Missouri, operates a well field in the Missouri River alluvial aquifer. Steady-state ground-water flow simulation, particle tracking, and the use of chemical and isotopic composition of river water, ground water, and well-field pumpage in a two-component mixing equation were used to determine the source contributions of induced inflow from the Missouri River and recharge to ground water from precipitation in well-field pumpage. Steady-state flow-budget analysis for the simulation-defined zone of contribution to the Independence well field indicates that 86.7 percent of well-field pumpage is from induced inflow from the river, and 6.7 percent is from ground-water recharge from precipitation. The 6.6 percent of flow from outside the simulation-defined zone of contribution is a measure of the uncertainty of the estimation, and occurs because model cells are too large to uniquely define the actual zone of contribution. Flow-budget calculations indicate that the largest source of water to most wells is the Missouri River. Particle-tracking techniques indicate that the Missouri River supplies 82.3 percent of the water to the Independence well field, ground-water recharge from precipitation supplies 9.7 percent, and flow from outside defined zones of contribution supplies 8.0 percent. Particle tracking was used to determine the relative amounts of source water to total well-field pumpage as a function of traveltime from the source. Well-field pumpage that traveled 1 year or less from the source was 8.8 percent, with 0.6 percent from the Missouri River, none from precipitation, and 8.2 percent between starting cells. Well-field pumpage that traveled 2 years or less from the source was 10.3 percent, with 1.8 percent from the Missouri River, 0.2 percent from precipitation, and 8.3 percent between starting cells. Well-field pumpage that traveled 5 years or less from the source was 36.5 percent, with 27.1 percent from the Missouri River, 1.1 percent from precipitation, and 8.3 percent between starting cells. Well-field pumpage that traveled 10 years or less from the source was 42.7 percent, with 32.6 percent from the Missouri River, 1.8 percent from precipitation, and 8.3 percent between starting cells. Well-field pumpage that traveled 25 years or less from the source was 71.9 percent, with 58.9 percent from the Missouri River, 4.7 percent from precipitation, and 8.3 percent between starting cells. Results of chemical (calcium, sodium, iron, and fluoride) and isotopic (oxygen and hydrogen) analyses of water samples collected from the Missouri River, selected monitoring wells around the Independence well field, and combined well-field pumpage were used in a two component mixing equation to estimate the relative amount of Missouri River water in total well-field pumpage. The relative amounts of induced inflow from the Missouri River in well-field pumpage ranged from 49 percent for sodium to 80 percent for calcium, and sensitivities ranged from 0 percent for iron to plus or minus 35 percent for naturally occurring stable isotope ( 18 O). The average of all mixing equation results indicated that 61 percent of well-field pumpage was from induced inflow from the Missouri River. All methods used in the study indicate that more than one-half of the water in well-field pumpage was inflow from the Missouri River. River inflow estimates from ground-water simulation methods are larger and error values are smaller than those using chemical and isotopic data in the mixing equation, although substantial uncertainties exist for both estimation methods. Because of the complex hydrology of the aquifer near the Independence well field, the source estimates using particle tracking probably are the most reliable of the ground-water simulation methods. Mixing equation results are less reliable than those of the ground-water simulation for this study. However, more reliable results can be obtained from the mixing equation by increasing the number of samples and collecting samples for a longer period of time, and during different flow conditions. In the absence of a calibrated ground-water flow simulation, the mixing equation can provide a reasonable estimate of the sources of water to a well field at relatively low cost, if sources of error are clearly understood.

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↗
USGS2001

Cross-hole radar scanning of two vertical, permeable, reactive-iron walls at the Massachusetts Military Reservation, Cape Cod, Massachusetts

A pilot-scale study was conducted by the U.S. Army National Guard (USANG) at the Massachusetts Military Reservation (MMR) on Cape Cod, Massachusetts, to assess the use of a hydraulic-fracturing method to create vertical, permeable walls of zero-valent iron to passively remediate ground water contaminated with chlorinated solvents. The study was conducted near the source area of the Chemical Spill-10 (CS-10) plume, a plume containing chlorinated solvents that underlies the MMR. Ground-water contamination near the source area extends from about 24 m (meters) to 35 m below land surface. The USANG designed two reactive-iron walls to be 12 m long and positioned 24 to 37 m below land surface to intersect and remediate part of the CS-10 plume.Because iron, as an electrical conductor, absorbs electromagnetic energy, the US Geological Survey used a cross-hole common-depth, radar scanning method to assess the continuity and to estimate the lateral and vertical extent of the two reactive-iron walls. The cross-hole radar surveys were conducted in boreholes on opposite sides of the iron injection zones using electric-dipole antennas with dominant center frequencies of 100 and 250 MHz. Significant decreases in the radar-pulse amplitudes observed in scans that traversed the injection zones were interpreted by comparing field data to results of two-dimensional finite-difference time-domain numerical models and laboratory-scale physical models.The numerical and physical models simulate a wall of perfectly conducting material embedded in saturated sand. Results from the numerical and physical models show that the amplitude of the radar pulse transmitted across the edge of a conductive wall is about 43 percent of the amplitude of a radar pulse transmitted across background material. The amplitude of a radar pulse transmitted through a hole in a conductive wall increases as the aperture of the hole increases. The modeling results indicate that holes with an aperture of less than 40 percent of the dominant wavelength of the radar pulse are not likely to be detected.Based on the results of the numerical and physical modeling, the decreases in radar-pulse amplitudes observed in scans traversing the injection zones are interpreted as electrically conductive zones that outline the distribution of iron. The area interpreted as iron in the northern A-wall contains two zones -- an upper zone about 10 m wide, extending from about 25 to 31 m below land surface, and a lower zone about 8 m wide, extending from 31.5 to 34.5 m below land surface. The area interpreted as iron in the southern B-wall is about 9 m wide, extending from about 27 to 34.5 m below land surface. No discrete holes were interpreted in either the A- or B-wall zones.The interpretation of the field data suggests that (1) the hydraulic-fracturing method introduced iron into the subsurface, but not in the dimensions originally proposed; (2) the iron within the treatment zones is distributed in a generally continuous manner; and (3) excluding the discontinuity in the A-wall, holes within the iron treatment zone, if any, exist at scales smaller than about 10 cm, the resolution limit of the radar antennas and acquisition geometry used for this study. The cross-hole radar method appears to have been an effective method for delineating the distribution of iron in the two walls; however, the veracity of the results cannot be ascertained without excavation or drilling into the treatment zone.

Massachusetts↗