Search USGS⌕ Search

SEARCH · Search USGS

Results for “Hydrological Data”

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 865 records · Page 48Linked to original sources

Baseline hydrologic studies in the lower Elwha River prior to dam removal

After the removal of two large, long‑standing dams on the Elwha River, Washington, the additional load of sediment and wood is expected to affect the hydrology of the lower river, its estuary, and the alluvial aquifer underlying the surrounding flood plain. To better understand the surface-water and groundwater characteristics of the river and estuary before dam removal, several hydrologic data sets were collected and analyzed. An experiment using a dye tracer characterized transient storage, and it was determined that the low‑flow channel of the lower Elwha River was relatively simple; 1–6 percent of the median travel time of dye was attributed to transient‑storage processes. Water data from monitoring wells adjacent to the main‑stem river indicated a strong hydraulic connectivity between stage in the river and groundwater levels in the flood plain. Analysis of temperature data from the monitoring wells showed that changes in the groundwater temperature responded weeks or months after water temperature changed in the river. A seepage investigation indicated that water from the river was moving into the aquifer (losing reach) between 1.7 and 2.8 kilometers from the river mouth. Surface‑water measurements and temperature and salinity data collected throughout the estuary helped to characterize the magnitude and nature of water movement in and out of the estuary. Salinity and stage sensors positioned in the estuarine network showed a strong surface‑water connection between the river and estuary waters east of the river. In contrast, there was a weaker connection between the river and estuarine water bodies west of the river.

Washington↗

Determination of flow losses in the Cape Fear River between B. Everett Jordan Lake and Lillington, North Carolina, 2008-2010

During 2008-2010, the U.S. Geological Survey conducted a hydrologic investigation in cooperation with the Triangle J Council of Governments Cape Fear River Flow Study Committee and the North Carolina Division of Water Resources to collect hydrologic data in the Cape Fear River between B. Everett Jordan Lake and Lillington in central North Carolina to help determine if suspected flow losses occur in the reach. Flow loss analyses were completed by summing the daily flow releases at Jordan Lake Dam with the daily discharges at Deep River at Moncure and Buckhorn Creek near Corinth, then subtracting these values from the daily discharges at Cape Fear River at Lillington. Examination of long-term records revealed that during 10,227 days of the 1983-2010 water years, 408 days (4.0 percent) had flow loss when conditions were relatively steady with respect to the previous day's records. The flow loss that occurred on these 40 days ranged from 0.49 to 2,150 cubic feet per second with a median flow loss of 37.2 cubic feet per second. The months with the highest number of days with flow losses were June (16. percent), September (16.9 percent), and October (19.4 percent). A series of synoptic discharge measurements made on six separate days in 2009 provided "snapshots" of overall flow conditions along the study reach. The largest water diversion is just downstream from the confluence of the Haw and Deep Rivers, and discharges substantially decrease in the main stem downstream from the intake point. Downstream from Buckhorn Dam, minimal gain or loss between the dam and Raven Rock State Park was noted. Analyses of discharge measurements and ratings for two streamgages-one at Deep River at Moncure and the other at Cape Fear River at Lillington-were completed to address the accuracy of the relation between stage and discharge at these sites. The ratings analyses did not indicate a particular time during the 1982-2011 water years in which a consistent bias occurred in the computations of discharge records that would indicate false flow losses. A total of 34 measured discharges at a streamgage on the Haw River below B. Everett Jordan Lake near Moncure were compared with the reported hourly flow releases from Jordan Lake Dam. Because 28 of 34 measurements were within plus or minus 10 percent of the hourly flow releases reported by the U.S Army Corps of Engineers, use of the current discharge computation tables for reporting Jordan Lake Dam flow releases is generally supported. A stage gage was operated on the Cape Fear River at Buckhorn Dam near Corinth to collect continuous stage-only records. Throughout the study period, flow over the dam was observed along its length, and flow loss within the study reach is not attributed to river-level fluctuations at the dam. Water-use information and (or) data were obtained for five industrial facilities, a regional power utility, two municipalities, one small hydropower facility on the Deep River, and one quarry operation also adjacent to the Deep River. The largest water users are the regional power producer, a small hydropower operation, and the two municipalities. The total water-use diversions for these facilities range from almost 25.5 to 38.5 cubic feet per second (39.5 to 59.5 million gallons per day) during the winter and summer periods, respectively. This range is equivalent to 69 to 104 percent of the 37 cubic feet per second median flow loss. The Lockville hydropower station is on the Deep River about 1 mile downstream from the streamgage near Moncure. Run-of-river operations at the facility do not appear to affect flow losses in the study reach. The largest water user in the study area is a regional power producer at a coal-fired power-generation plant located immediately adjacent to the Cape Fear River just downstream from the confluence of the Haw an Deep Rivers. Comparisons of daily water withdrawals, sup-plied by the regional power producer, and discharge records at a streamgage on the diversion canal indicated many days when consumption exceeded the producer's estimates for the cooling towers. Uncertainty surrounding reasonable estimates of consumption remained in effect at the end of the study. Data concerning evaporative losses were compiled using two approaches-an analysis of available pan-evaporation data from a National Weather Service cooperative observer station in Chapel Hill, North Carolina; and a compilation of reference open-water evaporation computed by the State Climate Office of North Carolina. The potential flow loss by evaporation from the main stem and the Deep River was estimated to be in the range of 4 to 14 cubic feet per second during May through October, equivalent to 10 to 38 percent of the 37 cubic feet per second median flow loss. Daily water-use diversions and evaporation losses were compared to flow-loss occurrences during the period April 2008 through September 2010. In comparing the surface-water, water-use, and evaporation data compiled for 2008-2010, it is evident that documented water diversions combined with flow losses by open-water evaporation can exceed the net flow gain in the study area and result in flow losses from the reach. Analysis of data from a streamgage downstream from the regional power plant on the diversion canal adjacent to the Cape Fear River provided insight into the occurrence of an apparent flow loss at the streamgage at Lillington. Assessment of the daily discharges and subsequent hydrographs for the canal streamgage indicated at least 24 instances during the study when the flows suddenly changed by magnitudes of 100 to more that 200 cubic feet per second, resulting in a noted time-lag effect on the downstream discharges at the Lillington streamgage, beginning 8 to 16 hours after the sudden flow change. A fiber-optic distributed temperature-sensing survey was conducted on the Cape Fear River at the Raven Rock State Park reach August 12-14, 2009, to determine if the presence of diabase dikes were preferentially directing groundwater discharge. No temperature anomalies of colder water were measured during the survey, which indicated that at the time of the survey that particular reach of the Cape Fear River was a "no-flow" or losing stream. An aerial thermal-infrared survey was conducted on the Haw and Cape Fear Rivers on February 27, 2010, from Jordan Lake Dam to Lillington to qualitatively delineate areas of groundwater discharge on the basis of the contrast between warm groundwater discharge and cold surface-water temperatures. Dis-charge generally was noted as diffuse seepage, but in a few cases springs were detected as inflow at a discrete point of discharge. Two reaches of the Cape Fear River (regional power plant and Bradley Road reaches) were selected for groundwater monitoring with a transect of piezometers installed within the flood plain. Groundwater-level altitudes at these reaches were analyzed for 1 water year (October 1, 2009, to September 30, 2010). Data collected as part of this study represent only a brief period of time and may not represent all conditions and all years; however, the data indicate that, during the dry summer months, the Cape Fear River within the study area is losing an undetermined quantity of water through seepage. Analyses completed during this investigation indicate a study reach with complex flow patterns affected by numerous concurrent factors resulting in flow losses. The causes of flow loss could not be solely attributed to any one factor. Among the factors considered, the occurrences of water diversions and evaporative losses were determined to be sufficient on some days (particularly during the base-flow period) to exceed the net gain in flows between the upstream and downstream ends of the study area. Losses by diversions and evaporation can exceed the median flow loss of 3 cubic feet per second, which indicates that flow loss from the study reach is real. Groundwater data collected during 2009-2010 indicate the possibility of localized flow loss during the summer, particularly in the impounded reach above Buckhorn Dam. However, no indication of unusual patterns was noted that would cause substantial flow loss by groundwater and surface-water interaction at the river bottom.

North Carolina↗

Hydrologic characteristics of low-impact stormwater control measures at two sites in northeastern Ohio, 2008-13

This report updates and examines hydrologic data gathered to characterize the performance of two stormwater-control measure (SCM) sites in the Chagrin River watershed, Ohio. At the Sterncrest Drive site, roadside bioswales and rain gardens were used to alleviate drainage problems in this residential neighborhood area. At the Washington Street site, a treatment train (including a pervious-paver system, rain garden, and bioswales) was used to reduce and delay stormwater runoff at a small business development. Selected metrics were used to demonstrate SCM system performance with regard to stormwater-management objectives at each site. Rain-garden overflow-frequency data collected at the Sterncrest Drive site during 2008–13 were used to characterize system sensitivity to rainfall characteristics. Approximately 70 percent of storms exceeding 0.75 inches during 3 hours or more resulted in overflows. Drainage-design features that may restrict flow through the system were identified. Overall, the data and local observations confirmed the continued success of the SCM at the Sterncrest Drive site in preventing roadway closure due to flooding. The additional years of data collected at the Washington Street site indicated that a previous analysis of increased runoff removal, based on only the first 2 years (2009–10) of data, provided premature conclusions. With 5 years of data (2009–13) and adjusting for changes in rainfall characteristics, it appears that the percentage of runoff removed by the system is decreasing; however, the lag time (time from onset of rainfall to runoff) has remained nearly constant. The annual mean percent removal for 2010–13 ranged from 55 to 37 percent with an overall mean of 45 percent, and this does meet the project objective of reducing runoff from the business complex. One possible explanation for the combination of increased volume of runoff and no change in the timing of runoff is the preferential flow paths developed in the SCM, increasing the capacity for internal drainage. Data indicated that the SCM system at the Washington Street site had reduced functionality over time.

Ohio↗

Hydraulic model and flood-inundation maps developed for the Pee Dee National Wildlife Refuge, North Carolina

A one-dimensional step-backwater model was developed by the U.S. Geological Survey (USGS) in cooperation with the U.S. Fish and Wildlife Service, Pee Dee National Wildlife Refuge, North Carolina, to provide a means for predicting flood-plain inundation. The model was developed for selected reaches of the Pee Dee River, Brown Creek, and Rocky River, using the U.S. Army Corps of Engineers Hydrologic Engineering Center River Analysis System (HEC-RAS) software. Multiple cross sections were defined on each modeled stream, and hydrologic data were collected between August 2011 and August 2013 at selected locations on the Pee Dee River and on its tributaries Brown Creek, Rocky River, and Thoroughfare Creek. Cross-section, stage, and flow data were used to develop the model and simulate water-surface profiles at 1.0-foot increments at the USGS streamgage Pee Dee River at Pee Dee Refuge near Ansonville, N.C. The profiles were produced for 31 selected water levels that ranged from approximately 193.0 feet to 223.0 feet in elevation at the Pee Dee River at Pee Dee Refuge streamgage. A series of digital flood-inundation maps were developed on the basis of the water-surface profiles produced by the model. The inundation maps, which can be accessed through the USGS Flood Inundation Mapping Program Web site at http://water.usgs.gov/osw/flood_inundation , depict estimates of the areal extent and depth of flooding corresponding to selected water levels at the USGS streamgage Pee Dee River at Pee Dee Refuge near Ansonville, N.C. These maps, when combined with real-time water-level information from USGS streamgages, provide managers with critical information to help plan flood-response activities and resource protection efforts.

North Carolina↗

Groundwater hydrology and estimation of horizontal groundwater flux from the Rio Grande at selected locations in Albuquerque, New Mexico, 2009–10

The Albuquerque area of New Mexico has two principal sources of water: (1) groundwater from the Santa Fe Group aquifer system, and (2) surface water from the Rio Grande. From 1960 to 2002, pumping from the Santa Fe Group aquifer system caused groundwater levels to decline more than 120 feet while water-level declines along the Rio Grande in Albuquerque were generally less than 40 feet. These differences in water-level declines in the Albuquerque area have resulted in a great deal of interest in quantifying the river-aquifer interaction associated with the Rio Grande. In 2003, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation, acting as fiscal agent for the Middle Rio Grande Endangered Species Collaborative Program, and the U.S. Army Corps of Engineers, began a study to characterize the hydrogeology of the Rio Grande inner valley alluvial aquifer in the Albuquerque area of New Mexico. The study provides hydrologic data in order to enhance the understanding of rates of water leakage from the Rio Grande to the alluvial aquifer, groundwater flow through the aquifer, and discharge of water from the aquifer to riverside drains. The study area extends about 20 miles along the Rio Grande in the Albuquerque area. Piezometers and surface-water gages were installed in paired transects at eight locations. Nested piezometers, completed at various depths in the alluvial aquifer, and surface-water gages, installed in the Rio Grande and riverside drains, were instrumented with pressure transducers. Water-level and water-temperature data were collected from 2009 to 2010. Water levels from the piezometers indicated that groundwater movement was usually away from the river towards the riverside drains. Annual mean horizontal groundwater gradients in the inner valley alluvial aquifer ranged from 0.0024 (I-25 East) to 0.0144 (Pajarito East). The median hydraulic conductivity values of the inner valley alluvial aquifer, determined from slug tests, ranged from 30 feet per day (ft/d) (Montaño) to 120 ft/d (Central) for paired transects, with a median hydraulic conductivity for all transects of 50 ft/d. Daily mean groundwater fluxes from the river through the inner valley alluvial aquifer computed using Darcy’s Law and the slug test results ranged from about 0.01 ft/d (Montaño West) to between 1.0 and 2.0 ft/d (Central East). Median annual groundwater fluxes from the river through the inner valley alluvial aquifer determined using the Suzuki-Stallman method was greatest at Alameda East (0.50 ft/d) and lowest at Alameda West (0.25 ft/d). The results from both methods agreed reasonably well. Seepage investigations conducted by measuring discharge in the east and west riverside drains provided information for computing changes in flow within the drains and for evaluating results from Darcy’s Law and Suzuki-Stallman method flux calculations. Discharge measured in the east riverside drain between the Barelas Bridge and the I-25 bridge indicated that the flow in the east riverside drain increased by an average of 56.5 cubic feet per day per linear foot (ft 3 /d/ft) of drain. Discharge measured in the west riverside drain between the Central bridge and the I-25 bridge indicated that flow increased between west drain miles 0 and 4, an average of 53.8 ft 3 /d/ft of drain, and that flow increased between west drain miles 7 and 10, an average of 44.9 ft 3 /d/ft of drain. In comparison to the seepage measurement results, the groundwater fluxes from the river through the inner valley alluvial aquifer calculated from Darcy’s Law (q slug ) and by the Suzuki-Stallman method (q heat ) would account for 20–36 percent or 53–95 percent, respectively, of the total flow in the east riverside drain and 22–31 percent or 19–26 percent, respectively, of the total flow in the west drain. These results indicate that the drains likely also receive water from outside the inner valley. The spatial variability of horizontal hydraulic gradients and groundwater fluxes can be primarily attributed to variability in the distances between the river and riverside drains throughout the study area and geologic heterogeneities in the alluvial aquifer. Temporal variability in the water levels, which control the horizontal hydraulic gradients and fluxes between the Rio Grande and the riverside drains, can be primarily attributed to seasonal fluctuations in river stage and irrigation practices.

New Mexico↗

Hydrology of the Claiborne aquifer and interconnection with the Upper Floridan aquifer in southwest Georgia

The U.S. Geological Survey conducted a study, in cooperation with the Georgia Environmental Protection Division, to define the hydrologic properties of the Claiborne aquifer and evaluate its connection with the Upper Floridan aquifer in southwest Georgia. The effort involved collecting and compiling hydrologic data from the aquifer in subarea 4 of southwestern Georgia. Data collected for this study include borehole geophysical logs in 7 wells, and two 72-hour aquifer tests to determine aquifer properties. The top of the Claiborne aquifer extends from an altitude of about 200 feet above the North American Vertical Datum of 1988 (NAVD 88) in Terrell County to 402 feet below NAVD 88 in Decatur County, Georgia. The base of the aquifer extends from an altitude of about 60 feet above NAVD 88 in eastern Sumter County to about 750 feet below NAVD 88 in Decatur County. Aquifer thickness ranges from about 70 feet in eastern Early County to 400 feet in Decatur County. The transmissivity of the Claiborne aquifer, determined from two 72-hour aquifer tests, was estimated to be 1,500 and 700 feet squared per day in Mitchell and Early Counties, respectively. The storage coefficient was estimated to be 0.0006 and 0.0004 for the same sites, respectively. Aquifer test data from Mitchell County indicate a small amount of leakage occurred during the test. Groundwater-flow models suggest that the source of the leakage was the underlying Clayton aquifer, which produced about 2.5 feet of drawdown in response to pumping in the Claiborne aquifer. The vertical hydraulic conductivity of the confining unit between the Claiborne and Clayton aquifers was simulated to be about 0.02 foot per day. Results from the 72-hour aquifer tests run for this study indicated no interconnection between the Claiborne and overlying Upper Floridan aquifers at the two test sites. Additional data are needed to monitor the effects that increased withdrawals from the Claiborne aquifer may have on future water resources.

Georgia↗

Hydrologic characterization of Bushy Park Reservoir, South Carolina, 2013–15

The Bushy Park Reservoir is a relatively shallow impoundment in a semi-tropical climate and is the principal water supply for the 400,000 people of the city of Charleston, South Carolina, and the surrounding areas including the Bushy Park Industrial Complex. Although there is an adequate supply of freshwater in the reservoir, taste-and-odor water-quality issues are a concern. The U.S. Geological Survey conducted an investigation in cooperation with the Charleston Water System to study the hydrology and hydrodynamics of the Bushy Park Reservoir to identify factors affecting water-quality conditions. Specifically, five areas for monitoring and (or) analysis were addressed: (1) hydrologic monitoring of the reservoir to establish a water budget, (2) flow monitoring in the tunnels to compute flow from Bushy Park Reservoir and at critical distribution junctions, (3) water-quality sampling, profiling, and continuous monitoring to identify the causes of taste-and-odor occurrence, (4) technical evaluation of appropriate hydrodynamic and water-quality simulation models for the reservoir, and (5) preliminary evaluation of alternative reservoir operations scenarios. This report describes the hydrodynamic and hydrologic data collected from 2013 to 2015 to support the application and calibration of a three-dimensional hydrodynamic model and the water-quality monitoring and analysis to gain insight into the principal causes of the Bushy Park Reservoir taste-and-odor episodes. The existing U.S. Geological Survey real-time network on the West Branch of the Cooper River was augmented with a tidal flow gage on Durham Canal Back River, and Foster Creek. The Charleston Water System intake structure was instrumented to collect water-level, water temperature (top and bottom probes), specific conductance (top and bottom probes), wind speed and direction, and photosynthetically active radiation data. In addition to the gages attached to fixed structures, four bottom-mounted velocity profilers were deployed at six locations over different periods. The deployment period for the velocity profiler ranged from 2 weeks to 4 months. During the investigation, tidal cycle (13-hour) streamflow measurements were made at 30-minute intervals at five locations. The Williams Station is a coal-fired powerplant that withdraws water from Bushy Park Reservoir for cooling purposes. The magnitude of the withdrawal (approximately 550 million gallons per day) is the major factor controlling the circulation in the reservoir. The net flow in Durham Canal to the reservoir is comparable to the withdrawal rates of the powerplant. When the Williams Station is not withdrawing water, the net flow in Durham Canal quickly goes to zero or reverses with a net flow away from the reservoir and to the Cooper River. Plan views of the velocity vectors for the tidal cycle streamflow measurements and rose diagram of the velocity profilers created with the Williams Station withdrawing and not withdrawing water show substantial effects of the distribution of magnitude and direction of the water velocities.

South Carolina↗

The Ozark Plateaus Regional Aquifer Study—Documentation of a groundwater-flow model constructed to assess water availability in the Ozark Plateaus

Recent short-term drought conditions have emphasized the need to better understand the delicate balance between abundance, sustainability, and scarcity of groundwater in the Ozark Plateaus aquifer system. In 2014, the U.S. Geological Survey began construction of a groundwater-flow model as a tool for the assessment of groundwater availability in the Ozark Plateaus aquifer system. The model was developed to benefit concurrent and future investigations involving groundwater-pumping scenarios, optimization, particle transport, and groundwater-monitoring network analysis. The groundwater model simulates 116 years (1900–2015) of hydrologic conditions and the response of the groundwater system to changes in stress including changes in recharge and groundwater pumping for water supply. Semiseasonal stress periods were simulated from the later part of 1991 to 2015 and represent higher demand and lower recharge in the spring and summer months and lower demand and higher recharge in the fall and winter months. Groundwater pumping increases throughout the simulation period with a maximum rate of about 600 million gallons per day (Mgal/d). The process of matching historical hydrologic data for the Ozark Plateaus aquifer system model was accomplished by a combination of manual changes to parameter values and automated calibration methods. Observation data used in the development and evaluation of the model included 19,045 hydraulic-head observations from 6,683 wells within the model area. Observation data also included stream leakage estimates summed to calculate a net gain or net loss value for approximately 81 named streams. The majority (mean of over 95 percent) of the recharge component is discharged through streams simulated in the model. The total simulated discharge to streams fluctuates seasonally between 7,500 and 17,500 Mgal/d with a mean outflow of 11,500 Mgal/d. Much of the remaining balance between modeled recharge inflows and stream outflows is made up by water moving into or out of storage in the aquifer system resulting in changes in modeled groundwater levels. The goal of the model was to develop a model capable of suitable accuracy at regional scales. The intent was not to reproduce individual local-scale details, which are typically not possible given the uniform cell size of 1 square mile. Although the model may not represent each local-scale detail, the model can be applied for a better understanding of the regional flow system and to evaluate responses to changes in climate and groundwater pumping.

Ozark Plateaus aquifer system↗

Trends in groundwater levels in and near the Rosebud Indian Reservation, South Dakota, water years 1956–2017

The U.S. Geological Survey (USGS), in cooperation with the Rosebud Sioux Tribe, completed a study to characterize water-level fluctuations in observation wells to examine driving factors that affect water levels in and near the Rosebud Indian Reservation, which comprises all of Todd County. The study investigates concerns regarding potential effects of groundwater withdrawals and climate conditions on groundwater levels within an area that includes Todd County and a surrounding area that extends 10 miles north, east, and west of the county border. Characterization of water-level fluctuations in observation wells and relative driving factors was accomplished by statistical trend analysis. Two statistical methods were used for analysis of temporal trends for climatic and hydrologic data. To determine which trend analysis to use, applicable datasets were tested for statistically significant short-term persistence (STP). In the absence of significant STP, existence of statistical trends was determined using the standard Mann-Kendall test for probability values less than or equal to 0.10 (90-percent confidence level); however, a modified Mann-Kendall test was used for datasets where statistically significant STP was detected. Trend magnitudes were computed using the Sen’s slope estimator. Monthly data from the Parameter-elevation Regressions on Independent Slopes Model (PRISM) were aggregated to obtain annual and seasonal datasets for total precipitation, minimum air temperature ( T min ), and maximum air temperature ( T max ) for the study area and a surrounding buffer area. Trend tests for total precipitation, T min , and T max were completed for annual and seasonal time series for water years 1956–2017, which is about 2 years before the earliest available water-level measurements. A 2-year offset was arbitrarily selected because scrutiny of water-level and precipitation data indicated that responses of groundwater levels for many of the observation wells lagged major changes in precipitation patterns by about 2 years. Statistically significant upward trends were detected for annual precipitation and annual T min for almost all of the study area and the surrounding buffer area. Statistically significant downward trends in T max were detected for a very small part of the study area; however, the sparse spatial coverage reduces confidence that these are true trends. Spatial distributions of statistically significant trends in seasonal climate data were generally similar to the annual trends, but with substantial differences in the spatial density of the trends. Groundwater trends for 58 observation wells were analyzed for three separate water-level parameters (minimum, median, and maximum) because wells are measured sporadically and data are biased towards more frequent measurements during periods of heaviest irrigation demand. Trends in the time series of annual precipitation (from PRISM) starting 2 years earlier than for the associated water-level trend also were analyzed for the location of each individual observation well. Sen’s slope and Mann-Kendall probability values (p-values) were computed for the three water-level parameters and for the annual precipitation time series. Graphs showing results of trend analyses for each observation well also showed changes over time in the sum of licensed groundwater withdrawals within six specified radii (0.5, 1, 2, 3, 4, and 5 miles) of each well as a qualitative indicator of proximal groundwater demand. Of all 58 observation wells considered, 28 wells had significant upward trends for at least one of the three water-level parameters, 11 wells had significant downward trends for at least one water-level parameter, and 19 wells did not have any significant trends. Significant upward trends in annual precipitation were detected for 48 of the 58 wells. Results of trend analyses likely show the effects of groundwater withdrawals on water levels in the Ogallala aquifer in areas of substantial demand. Precipitation trends are significantly upward for 43 of the 48 wells completed in the Ogallala aquifer that were analyzed. Of the 48 Ogallala aquifer wells, 24 had significant upward trends for at least one water-level parameter (17 with all 3); however, 10 wells had statistically significant downward trends for at least one water-level parameter (8 with all 3 parameters). All but one of the wells with significant downward trends are located in the south-central part of the study area where licensed irrigation withdrawals are concentrated.

South Dakota↗

Hydrologic analysis of an earthen embankment dam in southern Westchester County, New York

In 2001, the New York City Department of Environmental Protection installed 25 wells on the southern embankment of the Hillview Reservoir in Westchester County in an unsuccessful attempt to locate the source of a large seep (seep A) that began flowing continuously in 1999. In 2005, the U.S. Geological Survey began a cooperative study with the NYCDEP to characterize the hydrology of the local groundwater system and identify potential sources of seep A and other seeps on the embankment. At least two groundwater-flow zones—one shallow and the other deep—overlie the bedrock at the Hillview Reservoir in southern Westchester County, New York. Analyses of slug tests of wells drilled into the southern embankment of the reservoir were used to determine the three-dimensional distribution of hydraulic conductivity of the embankment materials. The wells with the minimum and maximum hydraulic conductivity values are in the deep saturated zone on the southern embankment, where hydraulic conductivity ranges from 0.0012 to 2 feet per day. Hydraulic conductivity ranges from 0.0026 to 1 foot per day in the shallow saturated zone and from 0.021 to 0.27 foot per day in the toe of the embankment. A hydraulic conductivity of 0.016 foot per day was determined for one toe well partially screened in the crystalline-bedrock aquifer. In 2005, the U.S. Geological Survey began a cooperative study with New York City Department of Environmental Protection to characterize the local groundwater-flow system and identify potential sources of seeps on the southern embankment of the Hillview Reservoir in southern Westchester County, New York. Long-term hydrologic data indicated that water levels trended downward in 29 of 41 sites, including the reservoir basin that was monitored during the 12-year study period; data from a National Weather Service precipitation gage at Central Park indicated annual precipitation also trended downward during the same 12-year period. Of the seven wells in which water levels trended upward during the study, two of the wells are on the west side of the southern embankment, proximal to a major water supply conduit, whereas the five remaining wells are screened in the toe. These data indicate an increasing hydrostatic pressure within the deep system and the toe of the dam, which could result in future seeps on the southern embankment near these wells. Results of 11 suspended-sediment samples collected from seeps along the southern embankment at 234.1- and 221.6-feet elevation, and another drainage outflow point between 2007 and 2015 indicate a poor correlation between suspended-sediment concentration and discharge. From the flowing seep at 234.1 feet, suspended-sediment concentrations ranged from 1 milligram per liter at a flow of 2.6 gallons per minute (that is, 1 milligram per 0.26 gallons) during March 2008 to 16 milligrams per liter at 12 gallons per minute during July 2014. At about 12 gallons per minute discharge, suspended-sediment concentration from samples collected at that seep during different sampling events, ranged from 3 to 16 milligrams per liter. From the seep at 221.6 feet elevation, the suspended sediment concentration was 2 milligrams per liter at a discharge of 3.4 gallons per minute and 2 milligrams per liter at a discharge of 1.1 gallons per minute. Only one sample was collected at the drainage outflow point, for which the suspended sediment concentration was 2 milligrams per liter at a discharge of 2.4 gallons per minute. Anomalously high-water levels were recorded in deep-system wells between June 5, 2013, and January 14, 2014. The period for the increase and the decrease back to more typical water-level elevations occurred rapidly during a 13-hour period in each instance. The sudden and rapid changes, in addition to the spatial distribution of magnitude of water-level response indicate that leaky water infrastructure was the source of recharge to the affected wells. A major water supply conduit was drained for repairs between July 7 and 10, 2010. The seeps indicated an immediate response and a substantial hydraulic connection to the water supply conduit. Approximately 10.5 hours after the water supply conduit was drained, flow from a seep on the southern embankment decreased from about 20 gallons per minute to less than 1 gallon per minute. This seep is located at about the same elevation and within the vicinity of the water supply conduit. A travel-time of about 10.5 hours from the source to the seep at 234.1 feet elevation was estimated from the dewatering timeline. During the 3-month shutdown of the water supply conduit, the previously flowing seeps remained dry until precipitation resulted in discharge of about 0.7 gallon per minute at the higher elevation seep, indicating a minor contribution from precipitation to the total seepage discharge. Discharge from the seeps resumed almost immediately coincident with the refilling of the water supply conduit, supporting the hydraulic connection observations during the drainage stage. In addition, during the refilling of the water supply conduit on September 21, 2010, a new seep (I) was observed on the southern embankment. Discharge from this new seep remained relatively constant until it became inaccessible under construction stone from subsequent embankment repairs by the New York City Department of Environmental Protection. Precipitation after the refilling stage of the shutdown seemed to have induced a rise in water levels in the toe wells and an increase in discharge from the seep at 234.1 feet elevation. The post shutdown discharge was less than 12 gallons per minute, compared to a discharge of about 20 gallons per minute before the repairs. The lower discharge rate measured during the period of historically higher discharge rates for the fall season indicates that the repair of the major water supply conduit may have contributed to a reduced discharge from the seeps. There were no definitive responses to the shutdown in any of the wells near the major water supply conduit. The more transmissive deep system of the southern embankment near the major water supply conduit and its associated infrastructure seems to be the preferential flow path for leaking infrastructure. The wells screened in this system showed a response during the deep system anomaly and have some of the highest hydraulic conductivities of the tested wells. All the seeps are in the elevation range of the deep system from approximately the crystalline bedrock surface around 200 feet elevation to the contact between the deep and shallow saturated zones of the reservoir at about 250 feet elevation.

New York↗

Hydrogeologic framework and conceptual groundwater-flow model of the panhandle and northwest parts of the High Plains (Ogallala) aquifer in Oklahoma, 1998–2022

This study was conducted by the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, to update the hydrogeologic framework and conceptual flow model for the panhandle and northwest parts of the High Plains (Ogallala) aquifer in Oklahoma, which together compose the Ogallala aquifer focus area. The study included the construction of a potentiometric surface, and available geologic and hydrologic data were used to evaluate saturated thickness of the aquifer. The water budget for the updated conceptual groundwater-​flow model was based on estimated inflows and outflows for the 1998–2022 study period. Saturated thickness of the Ogallala aquifer averaged 127 and 116 feet for the panhandle and northwest parts, respectively. Groundwater withdrawals from the Ogallala aquifer for 1998–2022 averaged 422,054 and 39,645 acre-feet per year (acre-ft/yr) for the panhandle and northwest parts, respectively. Recharge, the primary inflow, was estimated at 0.63 inch per year for the 1998–2022 study period, with the panhandle part of the Ogallala aquifer receiving 175,068 acre-ft/yr and the northwest part of the Ogallala aquifer receiving 49,376 acre-ft/yr. Additional inflows included irrigation return flows, estimated at 8,111 and 642 acre-ft/yr for the panhandle and northwest parts, respectively, of the Ogallala aquifer. Net lateral groundwater flows, considered to be aquifer outflows, were estimated to account for 31,908 acre-ft/yr for the Ogallala aquifer focus area. Streambed seepage, which was an outflow of 5,535 acre-ft/yr, was only present in the northwest part of the Ogallala aquifer. Vertical leakage and saturated-zone evapotranspiration were considered negligible outflows. These findings provide a revised conceptual groundwater-flow model water budget for the Ogallala aquifer focus area in Oklahoma.

Colorado, Kansas, Oklahoma, Texas↗

Occurrence of cyanobacteria and associated cyanotoxins in the Raritan Basin Water Supply Complex, New Jersey, August 2020 to August 2021

Harmful algal blooms, particularly cyanobacteria harmful algal blooms (cyanoHABs), have emerged as a substantial global concern because of their detrimental effects on water quality and aquatic ecosystem health. CyanoHABs can produce cyanotoxins, which pose serious health risks to humans and wildlife, such as liver failure and respiratory distress. This is particularly concerning for water bodies that serve as drinking-water sources. Recent trends indicate an increase in the frequency and intensity of cyanoHABs globally. This study focuses on the Raritan Basin Water Supply Complex in New Jersey, where extensive monitoring was conducted from August 2020 to August 2021 to assess the presence of cyanobacteria and associated cyanotoxins. The research utilized a combination of discrete water-quality sampling, continuous monitoring, and solid phase adsorption toxin tracking (SPATT) to capture the dynamics of cyanotoxin occurrence and potential transport. Findings revealed a widespread presence of cyanobacteria and potential for cyanotoxin production, although actual cyanotoxin concentrations remained below drinking water and recreational thresholds. The study, conducted by the U.S. Geological Survey (USGS) in collaboration with the New Jersey Water Supply Authority (NJWSA) and the New Jersey Department of Environmental Protection (NJDEP), highlighted the limitations of traditional sampling methods and emphasized that continuous monitoring can support better understanding of how cyanoHAB conditions change over time and in different places. Genetic testing included quantitative polymerase chain reaction (qPCR) analyses, which demonstrated higher sensitivity, or increased findings of cyanobacteria compared to microscopy, indicating the potential for use in early warning systems. This research underscores that integrating various detection methods and hydrological data can enhance understanding of cyanotoxin dynamics in river systems.

New Jersey↗

Model simulation of the Manasquan water-supply system in Monmouth County, New Jersey

Model simulation of the Manasquan Water Supply System in Monmouth County, New Jersey, was completed using historic hydrologic data to evaluate the effects of operational and withdrawal alternatives on the Manasquan reservoir and pumping system. Changes in the system operations can be simulated with the model using precipitation forecasts. The Manasquan Reservoir system model operates by using daily streamflow values, which were reconstructed from historical U.S. Geological Survey streamflow-gaging station records. The model is able to run in two modes--General Risk analysis Model (GRAM) and Position Analysis Model (POSA). The GRAM simulation procedure uses reconstructed historical streamflow records to provide probability estimates of certain events, such as reservoir storage levels declining below a specific level, when given an assumed set of operating rules and withdrawal rates. POSA can be used to forecast the likelihood of specified outcomes, such as streamflows falling below statutory passing flows, associated with a specific working plan for the water-supply system over a period of months. The user can manipulate the model and generate graphs and tables of streamflows and storage, for example. This model can be used as a management tool to facilitate the development of drought warning and drought emergency rule curves and safe yield values for the water-supply system.

Water-Resources Investigations Report↗

Ground-water levels in the Floridan-Midville aquifer in the Breezy Hill area, Aiken and Edgefield Counties, South Carolina, April 1999-November 2000

The Breezy Hill area in Aiken and Edgefield Counties of west-central South Carolina is a rapidly growing region in need of increasing amounts of ground water. From 1995 to 1998, the local water utility increased ground-water withdrawals in the Breezy Hill area from 1.4 to 2.1 million gallons per day to meet water-supply demands. As development continues, future demands for ground water will likely put stress on the surfaceand ground-water resources of the area. To address this issue, the U.S. Geological Survey, in cooperation with Aiken County, compiled and interpreted geologic and hydrologic data needed to map the ground-water system in the Breezy Hill study area. The Breezy Hill study area consists of four interfluvial areas comprising the regions between Horse and Little Horse Creeks, Little Horse and Hightower Creeks, Hightower Creek and Franklin Branch, and Franklin Branch and Mims Branch. Across the interfluvial areas, the average elevation of the water-level surface ranged from 200 to 480 feet above sea level, and the average saturated thickness of the Floridan-Midville aquifer ranged from less than 20 to 70 feet thick. A water-level contour map of the surface of the Floridan-Midville aquifer indicates that recharge to the aquifer occurs mainly within the interfluves. Recharge is derived principally from precipitation, although there is some potential for ground-water recharge from underlying crystalline rocks. Ground water discharges along the flanks of the interfluves into the bounding streams where the elevations of the ground water and streams coincide. From April 1999 to November 2000, calculated long-term normal precipitation totaled about 84.0 inches; however, actual recorded precipitation totaled 69.2 inches, representing about a 17.6 percent decrease in precipitation during this period. Published estimates of annual evapotranspiration range from 30 to 35 inches. A U.S. Geological Survey surface-water gaging station located near the center of the study area on Little Horse Creek monitors runoff from a drainage area of 26.6 square miles. Average annual flow for the station for water years 1990-2000 was 33.8 cubic feet per second. From April 1999 to November 2000, the monthly average flow was less than the average monthly flow for the longterm record, excluding December 1999 to March 2000 when no data were collected. Monthly average flow for Little Horse Creek exceeded the normal monthly flow during June and July 1999. Ground water in the Breezy Hill area is principally withdrawn from the unconfined Floridan- Midville aquifer. Ground-water withdrawals by the local water utility increased 37 percent from 1989 to 2000 (315.2 to 500 million gallons, respectively). From January 1999 to December 2000, the utility exceeded the long-term monthly average groundwater withdrawals for every month except September and December 2000. Calculated long-term monthly ground-water withdrawals by the utility for a 20-month period from April 1999 to November 2000 totaled 674 million gallons; however, actual ground-water withdrawals totaled 883 million gallons, which is 31 percent more than the long-term average ground-water withdrawals for the production wells. Published estimates of average annual ground-water recharge rates for the study area range from 13 to 15 inches per year. A base-flow recession analysis of streamflow data for Little Horse Creek provided an estimated recharge rate of 14.9 inches per year for the drainage area. Using an estimated average porosity ranging from 30 to 35 percent observed in sand-aquifer cores, the average annual recharge of 13 to 15 inches would cause a 3.6- to 4.1-foot water-level change to the saturated thickness of the aquifer, if applied instantaneously. The water-level declines observed in wells from April 1999 to November 2000 approximated an average decline of 4 feet. From November 1999 to November 2000, ground-water levels in six wells near utility pumping centers declined 2 to 5 feet. Long-term waterlevel declines of 10.27 and 11.50 feet were measured in two wells between May 1992 and April 2000, respectively.

South Carolina↗

Water quality of the Tlikakila River and five major tributaries to Lake Clark, Lake Clark National Park and Preserve, Alaska, 1999-2001

The Tlikakila River Basin, located in Lake Clark National Park and Preserve, drains an area of 622 square miles. This watershed comprises about 21 percent of the Lake Clark Basin, making it one of the major tributaries to Lake Clark. Due to a sharp decline in sockeye salmon population and the lack of hydrologic data, the Tlikakila River and five other major tributaries to Lake Clark were studied during the summer runoff months (May through September) from 1999 through 2001 as part of a cooperative study with the National Park Service. Measurements of pH and dissolved oxygen concentrations of the Tlikakila River are within acceptable limits for fish survival. Water temperatures at the measurement site reach 0 ?C during the winter and this part of the Tlikakila River may not be suitable for fish. Water temperatures are within acceptable limits for fish during the summer months. The Tlikakila River is a calcium bicarbonate type water with a low buffering capacity. Concentrations of un-ionized ammonia are well below the recommended value of 0.02 mg/L for fish propagation. Annual transport of suspended sediment by the Tlikakila River into Lake Clark ranged from 0.4 to 1.5 million tons during 1999?2001. The fine sediment from the Tlikakila River disperses through the lake over the summer, affecting light transmissivity. Most runoff from the Tlikakila River occurs from mid-to-late May through September. Average discharge for these months during 1999?2001 was 6,600 ft?/s. Total annual inflow to Lake Clark from the Tlikakila River ranged from 32 to 45 percent of the total inflow. The relatively high proportion of inflow is due to the presence of glaciers, which comprise 36 percent of the watershed. Monthly measurements of flow, field water-quality parameters, alkalinity, and suspended sediment were collected on the remaining five tributaries to Lake Clark: the Chokotonk River, Currant Creek, the Kijik River, the Tanalian River and the Chulitna River. Similar to the Tlikakila River, pH and dissolved oxygen concentrations of these rivers are within acceptable limits for fish survival and the rivers have a low buffering capacity. Small amounts of suspended sediment are transported by the Kijik and Tanalian Rivers due to lakes acting as settling basins in their watersheds. The Chulitna River also transports small amounts of suspended sediment due to its flat topography and the presence of many lakes in the basin. Some suspended sediment is transported by the Chokotonk River and Currant Creek during the runoff season due to the presence of glaciers within their basins, but not as much as the Tlikakila River. The Chulitna River provides the most discharge into Lake Clark after the Tlikakila River and has the warmest water temperature of the major tributaries to Lake Clark. Water temperatures of Currant Creek and the Chokotonk River are similar to the Tlikakila River. The Kijik River and Tanalian River have similar temperatures that may be due to the presence of lakes in their basins and are characterize by slowly declining and rising temperatures. At all sites water temperature approaches 0 ?C during winter months which may not be suitable for fish survival.

Water-Resources Investigations Report↗

Simulation of the shallow aquifer in the vicinity of Silver Lake, Washington County, Wisconsin, using analytic elements

Shallow ground-water flow in the vicinity of Silver Lake, Washington County, Wisconsin, was investigated to develop an understanding of the hydrology of the shallow aquifer, define a water balance for the lake, delineate ground-water recharge areas for the lake, and to estimate solute flux toward the lake. A single-layer, steady-state, analytic-element model was used to simulate shallow ground-water flow. Regional model parameters include a recharge rate of 4 inches per year, hydraulic conductivity of 50 feet per day and a model base of 800 feet above sea level. A model inhomogeneity was added to represent deviations from these regional values for an area roughly coincident with the Kettle Moraine Area that trends through the study area. Model calibration was accomplished by varying the regional parameter values and those of the inhomogeneity through trial-and-error to determine a best-fit match between simulated and measured values for head and streamflow targets. There was no change to the regional parameter values as a result of calibration, however, the calibrated values for the inhomogeneity are: recharge rate of 12 inches per year, hydraulic conductivity of 20 feet per day, and a model base of 900 feet. These changes represent a four- to five-fold reduction in transmissivity within the inhomogeneity as compared to the regional model. A Silver Lake water budget was defined using both published hydrologic data and simulations using the calibrated model. Model simulations show that 1.08 cubic feet per second of ground water enters Silver Lake on the upgradient (primarily western) side and 0.08 cubic feet per second recharges to ground water on the downgradient (primarily eastern) side. Net precipitation (precipitation minus evaporation) on the lake is 0.04 cubic feet per second. Collectively, these water-budget terms provide a residual value of 1.04 cubic feet per second flow to Silver Creek at the north end of Silver Lake, which is a very good match to the range of measured flow (0.7 to 5.2 cubic feet per second). Ground-water recharge areas for Silver Lake are largely on the western side of the lake. The recharge area for the northern two-thirds of Silver Lake is west toward Big Cedar Lake. Assuming a porosity of 20 percent, model results indicate that the 50-year time-of-travel for recharge to Silver Lake does not extend to Big Cedar Lake. The recharge area for the southern one-third of Silver Lake is west toward Little Cedar Lake. Model results indicate that time of travel for recharge to Silver Lake from Little Cedar Lake is about 15 to 20 years. For travel times greater than 15 or 20 years, the ground-water recharge area for Little Cedar Lake and inflow from Big Cedar Lake also should be considered recharge affecting Silver Lake. Solute flux toward Silver Lake was calculated based on simulated ground-water flux and measured concentrations in the upgradient piezometers and observation wells.

Wisconsin↗

Water resources of Monroe County, New York, water years 1997-99, with emphasis on water quality in the Irondequoit Creek basin—Atmospheric deposition, ground water, streamflow, trends in water quality, and chemical loads to Irondequoit Bay

Irondequoit Creek drains 169 square miles in the eastern part of Monroe County. Over time, nutrients transported by Irondequoit Creek to Irondequoit Bay on Lake Ontario have contributed to the eutrophication of the bay. Sewage-treatment-plant effluent, a major source of nutrients to the creek and its tributaries, was eliminated from the basin in 1979 by diversion to a regional wastewater-treatment facility, but sediment and contaminants from nonpoint sources continue to enter the creek and Irondequoit Bay. This report, the fourth in a series of reports that present interpretive analyses of the hydrologic data collected in Monroe County since 1984, interprets data from four surface-water monitoring sites in the Irondequoit Creek basin—Irondequoit Creek at Railroad Mills, East Branch Allen Creek at Pittsford, Allen Creek near Rochester, and Irondequoit Creek at Blossom Road. It also interprets data from three sites in the the Genesee River basin—Oatka Creek at Garbutt, Honeoye Creek at Honeoye Falls, and Black Creek at Churchville—as well as the Genesee River at Charlotte Pump Station, and also from a site on Northrup Creek at North Greece. The Northrup Creek site drains a 23.5-square-mile basin in western Monroe County, and provides information on surface-water quality in streams west of the Genesee River and on loads of nutrients delivered to Long Pond, a small eutrophic embayment of Lake Ontario. The report also includes water-level and water-quality data from nine observation wells in Ellison Park, and atmospheric-deposition data from a collection site at Mendon Ponds County Park. Average annual loads of some chemical constituents in atmospheric deposition for 1997–99 differed considerably from those for the long-term period 1984–96. Ammonia and potassium loads for 1997-99 were 144 and 118 percent greater, respectively, than for the previous period. Sodium and ammonia + organic nitrogen loads were 87 and 60 percent greater, respectively. Average annual loads of sulfate and orthophosphate for 1997-99 were 36 and 30 percent lower, respectively, than for the previous period. Loads of all nutrients deposited on the Irondequoit basin from atmospheric sources during 1997–99 greatly exceeded those transported by Irondequoit Creek. The ammonia load deposited on the basin was 139 times the load transported at Blossom Road (the most downstream site); the ammonia + organic nitrogen load was 6.3 times greater, orthophosphate 7.5 times greater, total phosphorus 1.3 times greater and nitrite + nitrate 1.5 times greater. Average yields of dissolved chloride and dissolved sulfate from atmospheric sources were much smaller than those transported by streamflow at Blossom Road.chloride was about 2 percent and sulfate about 8 percent of the amount transported. Trends in concentration of chemical constituents in surface water generally can be attributed to changes in land use, annual and seasonal variations in streamflow, and annual variations in the application of road salt to county highways and roads. Concentrations of several constituents in streams of the Irondequoit Creek basin showed statistically significant (α=0.05) trends from the beginning of their period of record through 1999. The constituent with the greatest number of significant trends was ammonia + organic nitrogen, with downward trends ranging from 4.1 to 5.6 percent per year at Allen Creek, Irondequoit Creek at Blossom Road, and East Branch Allen Creek. Orthophosphate showed an upward trend of 4.1 percent per year at Irondequoit Creek at Railroad Mills (the most upstream site). Dissolved chloride showed upward trends at Railroad Mills, Allen Creek, and Blossom Road. No trends in volatile suspended solids were noted at any of the four Irondequoit basin sites. Northrup Creek showed significant downward trends in concentrations of ammonia + organic nitrogen (3.3 percent per year), total phosphorus (3.4 percent per year), and orthophosphate (5.5 percent per year), and an upward trend for dissolved sulfate (1.8 percent per year). The Genesee River at Charlotte Pump Station showed downward trends of 6.1 percent per year for ammonia + organic nitrogen and 0.1 percent per year for chloride, and upward trends of 1.7 percent per year for total phosphorus and 6.6 percent per year for orthophosphate. Mean annual yields (mass per unit area) of most constituents at the Irondequoit Creek basin sites were similar to those noted for the previous report period (1994–96). East Branch Allen Creek showed lower yields of all constituents during 1997–99 than previously, even though runoff during 1997–99 was greater. These lower yields are attributed to the construction of an upstream detention basin on East Branch Allen Creek in 1995. Statistical analysis of long-term (greater than 12 years) streamflow records for unregulated streams in Monroe County indicated that annual mean flows for water years 1997–99 were in the normal range (75th to 25th percentile), although Allen Creek continues to show a significant downward trend in mean monthly streamflow during the 1984–99 water years.

New York↗

Simulation of ground-water flow and evaluation of water-management alternatives in the upper Charles River basin, eastern Massachusetts

Ground water is the primary source of drinking water for towns in the upper Charles River Basin, an area of 105 square miles in eastern Massachusetts that is undergoing rapid growth. The stratified-glacial aquifers in the basin are high yield, but also are thin, discontinuous, and in close hydraulic connection with streams, ponds, and wetlands. Water withdrawals averaged 10.1 million gallons per day in 1989?98 and are likely to increase in response to rapid growth. These withdrawals deplete streamflow and lower pond levels. A study was conducted to develop tools for evaluating water-management alternatives at the regional scale in the basin. Geologic and hydrologic data were compiled and collected to characterize the ground- and surface-water systems. Numerical flow modeling techniques were applied to evaluate the effects of increased withdrawals and altered recharge on ground-water levels, pond levels, and stream base flow. Simulation-optimization methods also were applied to test their efficacy for management of multiple water-supply and water-resource needs. Steady-state and transient ground-water-flow models were developed using the numerical modeling code MODFLOW-2000. The models were calibrated to 1989?98 average annual conditions of water withdrawals, water levels, and stream base flow. Model recharge rates were varied spatially, by land use, surficial geology, and septic-tank return flow. Recharge was changed during model calibration by means of parameter-estimation techniques to better match the estimated average annual base flow; area-weighted rates averaged 22.5 inches per year for the basin. Water withdrawals accounted for about 7 percent of total simulated flows through the stream-aquifer system and were about equal in magnitude to model-calculated rates of ground-water evapotranspiration from wetlands and ponds in aquifer areas. Water withdrawals as percentages of total flow varied spatially and temporally within an average year; maximum values were 12 to 13 percent of total annual flow in some subbasins and of total monthly flow throughout the basin in summer and early fall. Water-management alternatives were evaluated by simulating hypothetical scenarios of increased withdrawals and altered recharge for average 1989?98 conditions with the flow models. Increased withdrawals to maximum State-permitted levels would result in withdrawals of about 15 million gallons per day, or about 50 percent more than current withdrawals. Model-calculated effects of these increased withdrawals included reductions in stream base flow that were greatest (as a percentage of total flow) in late summer and early fall. These reductions ranged from less than 5 percent to more than 60 percent of model-calculated 1989?98 base flow along reaches of the Charles River and major tributaries during low-flow periods. Reductions in base flow generally were comparable to upstream increases in withdrawals, but were slightly less than upstream withdrawals in areas where septic-system return flow was simulated. Increased withdrawals also increased the proportion of wastewater in the Charles River downstream of treatment facilities. The wastewater component increased downstream from a treatment facility in Milford from 80 percent of September base flow under 1989?98 conditions to 90 percent of base flow, and from 18 to 27 percent of September base flow downstream of a treatment facility in Medway. In another set of hypothetical scenarios, additional recharge equal to the transfer of water out of a typical subbasin by sewers was found to increase model-calculated base flows by about 12 percent of model-calculated base flows. Addition of recharge equal to that available from artificial recharge of residential rooftop runoff had smaller effects, augmenting simulated September base flow by about 3 percent. Simulation-optimization methods were applied to an area near Populatic Pond and the confluence of the Mill and Charles Rivers in Franklin,

Massachusetts↗