Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Scientific Monthly”

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 649 records · Page 36Linked to original sources

Coeur d'Alene Lake, Idaho: Insights gained From limnological studies of 1991-92 and 2004-06

More than 100 years of mining and processing of metal-rich ores in northern Idaho’s Coeur d’Alene River basin have resulted in widespread metal contamination of the basin’s soil, sediment, water, and biota, including Coeur d’Alene Lake. Previous studies reported that about 85 percent of the bottom of Coeur d’Alene Lake is substantially enriched in antimony, arsenic, cadmium, copper, lead, mercury, silver, and zinc. Nutrients in the lake also are a major concern because they can change the lake’s trophic status—or level of biological productivity—which could result in secondary releases of metals from contaminated lakebed sediments. This report presents insights into the limnological functioning of Coeur d’Alene Lake based on information gathered during two large-scale limnological studies conducted during calendar years 1991–92 and water years 2004–06. Both limnological studies reported that longitudinal gradients exist from north to south for decreasing water column transparency, loss of dissolved oxygen, and increasing total phosphorus concentrations. Gradients also exist for total lead, total zinc, and hypolimnetic dissolved oxygen concentrations, ranging from high concentrations in the central part of the lake to lower concentrations at the northern and southern ends of the lake. In the southern end of the lake, seasonal anoxia serves as a mechanism to release dissolved constituents such as phosphorus, nitrogen, iron, and manganese from lakebed sediments and from detrital material within the water column. Nonparametric statistical hypothesis tests at a significance level of α=0.05 were used to compare analyte concentrations among stations, between lake zones, and between study periods. The highest dissolved oxygen concentrations were measured in winter in association with minimum water temperatures, and the lowest concentrations were measured in the Coeur d’Alene Lake hypolimnion during late summer or autumn as prolonged thermal stratification restricted mixing of the oxygenated upper water column and the hypolimnion, where oxygen was consumed. Large differences in median concentrations of dissolved inorganic nitrogen were measured between the euphotic zone and hypolimnion in the deep areas of the lake. These differences in nitrogen concentrations were attributable to several limnological processes, including seasonal inflow plume routing, isolation from wind-driven circulation and associated hypolimnetic enrichment, phytoplanktonic assimilation during summer months, and benthic flux. Increased chlorophyll- a and total phosphorus concentrations were measured throughout the lake in the 2004–06 study compared with results from the 1991–92 study. No significant change in hypolimnetic dissolved inorganic nitrogen concentration throughout the lake was noted even though total nitrogen loads into the lake decreased between study periods. Total zinc and total lead decreased throughout the lake from the 1991-92 study to the 2004-06 study except in the southern part of the lake, where concentrations were typically low. Median detected nitrogen-to-phosphorus ratios decreased from the 1991–92 study to the 2004–06 study. Whereas the lake was clearly phosphorus-limited in 1991–92, in 2004–06 the lake may have been much closer to the boundary value of 7.2 that separates nitrogen from phosphorus limitation. However, due to changes in analytical reporting limits in the period between the two studies, the data are insufficiently certain to draw reliable conclusions with regard to limiting nutrients. For both studies, the trophic state of the lake was classified as oligotrophic (less productive) or mesotrophic (moderately productive), depending on the constituent used for classification. Internal circulation from wind-generated waves and changes in the lake’s thermocline are important processes for distribution of water-quality constituents throughout Coeur d’Alene Lake. Surficial distribution of trace metals throughout most of the lake, including bays, is relatively uniform. Even south of the Coeur d’Alene River mouth, lakebed sediments are contaminated with trace metals. Inflow plume routing of the two primary inflow sources, the Coeur d’Alene and St. Joe Rivers, also significantly affects the fate and transport of contaminants. Most riverine inflows appear to move through the lake as overflow during summer, interflow during spring and autumn, and underflow during winter. Benthic flux is another key transport process for contaminants in Coeur d’Alene Lake. The results of in situ benthic flux chamber experiments indicated movement of dissolved metals, nutrients, and dissolved organic carbon out of the lakebed sediments. However, the lake is an overall sink for these constituents when they are associated with particulate material.

Idaho↗

Comparison of methods for estimating ground-water recharge and base flow at a small watershed underlain by fractured bedrock in the Eastern United States

This study by the U.S. Geological Survey (USGS), in cooperation with the Agricultural Research Service (ARS), U.S. Department of Agriculture, compared multiple methods for estimating ground-water recharge and base flow (as a proxy for recharge) at sites in east-central Pennsylvania underlain by fractured bedrock and representative of a humid-continental climate. This study was one of several within the USGS Ground-Water Resources Program designed to provide an improved understanding of methods for estimating recharge in the eastern United States. Recharge was estimated on a monthly and annual basis using four methods?(1) unsaturated-zone drainage collected in gravity lysimeters, (2) daily water balance, (3) water-table fluctuations in wells, and (4) equations of Rorabaugh. Base flow was estimated by streamflow-hydrograph separation using the computer programs PART and HYSEP. Estimates of recharge and base flow were compared for an 8-year period (1994-2001) coinciding with operation of the gravity lysimeters at an experimental recharge site (Masser Recharge Site) and a longer 34-year period (1968-2001), for which climate and streamflow data were available on a 2.8-square-mile watershed (WE-38 watershed). Estimates of mean-annual recharge at the Masser Recharge Site and WE-38 watershed for 1994-2001 ranged from 9.9 to 14.0 inches (24 to 33 percent of precipitation). Recharge, in inches, from the various methods was: unsaturated-zone drainage, 12.2; daily water balance, 12.3; Rorabaugh equations with PULSE, 10.2, or RORA, 14.0; and water-table fluctuations, 9.9. Mean-annual base flow from streamflow-hydrograph separation ranged from 9.0 to 11.6 inches (21-28 percent of precipitation). Base flow, in inches, from the various methods was: PART, 10.7; HYSEP Local Minimum, 9.0; HYSEP Sliding Interval, 11.5; and HYSEP Fixed Interval, 11.6. Estimating recharge from multiple methods is useful, but the inherent differences of the methods must be considered when comparing results. For example, although unsaturated-zone drainage from the gravity lysimeters provided the most direct measure of potential recharge, it does not incorporate spatial variability that is contained in watershed-wide estimates of net recharge from the Rorabaugh equations or base flow from streamflow-hydrograph separation. This study showed that water-level fluctuations, in particular, should be used with caution to estimate recharge in low-storage fractured-rock aquifers because of the variability of water-level response among wells and sensitivity of recharge to small errors in estimating specific yield. To bracket the largest range of plausible recharge, results from this study indicate that recharge derived from RORA should be compared with base flow from the Local-Minimum version of HYSEP.

Scientific Investigations Report↗

Methods for estimating selected low-flow frequency statistics and mean annual flow for ungaged locations on streams in North Georgia

The U.S. Geological Survey, in cooperation with the Georgia Department of Natural Resources, Environmental Protection Division, developed regional regression equations for estimating selected low-flow frequency and mean annual flow statistics for ungaged streams in north Georgia that are not substantially affected by regulation, diversions, or urbanization. Selected low-flow frequency statistics and basin characteristics for 56 streamgage locations within north Georgia and 75 miles beyond the State’s borders in Alabama, Tennessee, North Carolina, and South Carolina were combined to form the final dataset used in the regional regression analysis. Because some of the streamgages in the study recorded zero flow, the final regression equations were developed using weighted left-censored regression analysis to analyze the flow data in an unbiased manner, with weights based on the number of years of record. The set of equations includes the annual minimum 1- and 7-day average streamflow with the 10-year recurrence interval (referred to as 1Q10 and 7Q10), monthly 7Q10, and mean annual flow. The final regional regression equations are functions of drainage area, mean annual precipitation, and relief ratio for the selected low-flow frequency statistics and drainage area and mean annual precipitation for mean annual flow. The average standard error of estimate was 13.7 percent for the mean annual flow regression equation and ranged from 26.1 to 91.6 percent for the selected low-flow frequency equations. The equations, which are based on data from streams with little to no flow alterations, can be used to provide estimates of the natural flows for selected ungaged stream locations in the area of Georgia north of the Fall Line. The regression equations are not to be used to estimate flows for streams that have been altered by the effects of major dams, surface-water withdrawals, groundwater withdrawals (pumping wells), diversions, or wastewater discharges. The regression equations should be used only for ungaged sites with drainage areas between 1.67 and 576 square miles, mean annual precipitation between 47.6 and 81.6 inches, and relief ratios between 0.146 and 0.607; these are the ranges of the explanatory variables used to develop the equations. An attempt was made to develop regional regression equations for the area of Georgia south of the Fall Line by using the same approach used during this study for north Georgia; however, the equations resulted with high average standard errors of estimates and poorly predicted flows below 0.5 cubic foot per second, which may be attributed to the karst topography common in that area. The final regression equations developed from this study are planned to be incorporated into the U.S. Geological Survey StreamStats program. StreamStats is a Web-based geographic information system that provides users with access to an assortment of analytical tools useful for water-resources planning and management, and for engineering design applications, such as the design of bridges. The StreamStats program provides streamflow statistics and basin characteristics for U.S. Geological Survey streamgage locations and ungaged sites of interest. StreamStats also can compute basin characteristics and provide estimates of streamflow statistics for ungaged sites when users select the location of a site along any stream in Georgia.

Georgia↗

Evaluation of potential stresses and hydrologic conditions driving water-level fluctuations in well ER-5-3-2, Frenchman Flat, southern Nevada

Well ER-5-3-2 is part of a well network designed to monitor long-term water levels and radionuclide concentrations downgradient from underground nuclear tests that occurred in Frenchman Flat, an area of the U.S. Department of Energy Nevada National Security Site in southern Nevada. Interpretation of monitoring records for well ER-5-3-2 was confounded by previously unexplained water-level fluctuations in the well hydrograph. This study integrated geologic, hydrologic, and water-chemistry data to evaluate potential stresses and hydrologic conditions that likely affected the well ER-5-3-2 hydrograph. Numerical groundwater models were applied to evaluate four model scenarios: (1) wellbore leakage without recharge, (2) wellbore leakage with recharge, (3) equilibration to vertical heterogeneities between shallow (low transmissivity) and deep (higher transmissivity) carbonate zones, and (4) equilibration to lateral heterogeneities in carbonate rocks. Meteoric recharge was not the cause of the 21-foot (ft) water-level rise in well ER-5-3-2 from 2001 to 2011 or the 4-ft decline from 2012 to 2016. Based on observed water-level fluctuations in nearby wells, the water-level rise and decline from recharge for these periods was less than 3 and 1 ft, respectively. The lateral-heterogeneity scenario is based on the assumption that the 21-ft water-level rise from 2001 to 2011 was a natural water-level reequilibration following the pumping-induced depressurization of a large volume of high transmissivity and low-storage carbonate rock that is surrounded by low transmissivity and high-storage carbonate rock. The lateral-heterogeneity scenario was discounted because simulated water levels cannot match the well ER-5-3-2 hydrograph. Underground nuclear testing and temperature effects were discounted based on hydraulic connections and water-temperature data. Wellbore-leakage scenarios are based on the assumption that the water-level rise was sustained from leakage rates required to cause a localized mounding in the carbonate system near well ER-5-3-2, where the carbonate transmissivity is 530 square feet per day. Even though simulated and measured water levels compare favorably for scenarios of wellbore leakage with and without recharge, large volumes (178–184 million gallons) of groundwater from volcanic rocks would be required to leak into the carbonate system, which is not supported by water-chemistry data. An alternative conceptualization of wellbore leakage is based on the assumption that the 21-ft water-level rise from 2001 to 2011 was sustained by the hydraulic disconnection of well ER-5-3-2 from the carbonate system. The disconnection occurred several months after a constant-rate test in well ER-5-3-2 when carbonate rocks were hydraulically disconnected from the well by either (1) the shifting of sloughed fill in the open hole or (2) the encrusting of carbonate precipitate in the well screen. The hydraulic disconnection effectively sealed the well and caused a 21-ft water-level rise from wellbore leakage during 2001–11. In this case, total wellbore leakage from 2001 to 2011 was about 50 gallons. The 4-ft water-level decline from 2012 to 2016 was conceptualized to have occurred from the slow breaking of the seal and reconnection of the well to the carbonate system. This alternative conceptualization of wellbore leakage was consistent with water-chemistry analyses because the computed wellbore leakage (50 gallons) was small relative to purged volumes (30,000–40,000 gallons) for sampling, and the water chemistry would not be expected to change. The shallow-deep carbonate scenario provided another explanation for the well ER-5-3-2 hydrograph. This scenario is based on the assumption that well-construction effects and vertical heterogeneity of the carbonate system explain the ER-5-3-2 water-level trend. Well-construction effects are attributed to a temporary clogging of the open interval below the well screen that was opened during pumping events, which affected the hydraulic connection of deep transmissive carbonate rocks to the wellbore. The 21-ft water-level rise from 2001 to 2011 was a natural equilibration to shallow, low-transmissivity carbonate rocks during a period when the lower open interval was clogged. The 4-ft decline from 2012 to 2016 represents equilibration between the shallow and deep intervals, because of a partial unclogging of the connection between the two intervals. The low water levels from 2016 to 2021 resulted from pumping for sampling and an unclogging of the open interval so that the low head in the deep carbonate dominated the water level. Despite potential well-construction effects, from either a wellbore leakage or shallow-deep carbonate scenario, samples collected from well ER-5-3-2 are representative of the carbonate system.

Nevada↗

Numerical simulation of the groundwater-flow system of the Kitsap Peninsula, west-central Washington

A groundwater-flow model was developed to improve understanding of water resources on the Kitsap Peninsula. The Kitsap Peninsula is in the Puget Sound lowland of west-central Washington, is bounded by Puget Sound on the east and by Hood Canal on the west, and covers an area of about 575 square miles. The peninsula encompasses all of Kitsap County, Mason County north of Hood Canal, and part of Pierce County west of Puget Sound. The peninsula is surrounded by saltwater, and the hydrologic setting is similar to that of an island. The study area is underlain by a thick sequence of unconsolidated glacial and interglacial deposits that overlie sedimentary and volcanic bedrock units that crop out in the central part of the study area. Twelve hydrogeologic units consisting of aquifers, confining units, and an underlying bedrock unit form the basis of the groundwater-flow model. Groundwater flow on the Kitsap Peninsula was simulated using the groundwater-flow model, MODFLOW‑NWT. The finite difference model grid comprises 536 rows, 362 columns, and 14 layers. Each model cell has a horizontal dimension of 500 by 500 feet, and the model contains a total of 1,227,772 active cells. Groundwater flow was simulated for transient conditions. Transient conditions were simulated for January 1985–December 2012 using annual stress periods for 1985–2004 and monthly stress periods for 2005–2012. During model calibration, variables were adjusted within probable ranges to minimize differences between measured and simulated groundwater levels and stream baseflows. As calibrated to transient conditions, the model has a standard deviation for heads and flows of 47.04 feet and 2.46 cubic feet per second, respectively. Simulated inflow to the model area for the 2005–2012 period from precipitation and secondary recharge was 585,323 acre-feet per year (acre-ft/yr) (93 percent of total simulated inflow ignoring changes in storage), and simulated inflow from stream and lake leakage was 43,905 acre-ft/yr (7 percent of total simulated inflow). Simulated outflow from the model primarily was through discharge to streams, lakes, springs, seeps, and Puget Sound (594,595 acre-ft/yr; 95 percent of total simulated outflow excluding changes in storage) and through withdrawals from wells (30,761 acre-ft/yr; 5 percent of total simulated outflow excluding changes in storage). Six scenarios were formulated with input from project stakeholders and were simulated using the calibrated model to provide representative examples of how the model could be used to evaluate the effects on water levels and stream baseflows of potential changes in groundwater withdrawals, in consumptive use, and in recharge. These included simulations of a steady-state system, no-pumping and return flows, 15-percent increase in current withdrawals in all wells, 80-percent decrease in outdoor water to simulate effects of conservation efforts, 15-percent decrease in recharge from precipitation to simulate a drought, and particle tracking to determine flow paths. Changes in water-level altitudes and baseflow amounts vary depending on the stress applied to the system in these various scenarios. Reducing recharge by 15 percent between 2005 and 2012 had the largest effect, with water-level altitudes declining throughout the model domain and baseflow amounts decreasing by as much as 18 percent compared to baseline conditions. Changes in pumping volumes had a smaller effect on the model. Removing all pumping and resulting return flows caused increased water-level altitudes in many areas and increased baseflow amounts of between 1 and 3 percent.

Washington↗

Water resources of Monroe County, New York, water years 2000-02: Atmospheric deposition, ground water, streamflow, trends in water quality, and chemical loads in streams

This report, the fifth in a series that presents 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 above Blossom Road); and from three sites on tributaries to the Genesee River (Oatka Creek at Garbutt, Honeoye Creek at Honeoye Falls, and Black Creek at Churchville) and from the Genesee River at Charlotte Docks. It also interprets data from a site on Northrup Creek, which provides information on nutrient loads 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. Atmospheric Deposition : Average annual precipitation for 2000–02 was 33.11 in., 0.94 in. below normal. Average annual loads of some chemical constituents in atmospheric deposition for 2000–02 differed considerably from those for the previous period of record. Loads of all nutrients except ammonia decreased by amounts ranging from 28 percent (ammonia + organic nitrogen and phosphorus) to 2 percent (nitrite + nitrate), whereas ammonia loads an increased by 8 percent. Loads of dissolved sodium and total zinc in atmospheric deposition increased by 56 percent, and 54, percent respectively, over the previous period of record. Average annual loads of other constituents showed decreases ranging from 41 percent (dissolved magnesium) to 17 percent (dissolved chloride). Loads of all nutrients deposited in the Irondequoit Creek basin from atmospheric sources during 2000–02 greatly exceeded those transported by Irondequoit Creek. The ammonia load deposited in the basin was 165 times the load transported at Blossom Road (the most downstream site); the ammonia + organic nitrogen load was 2.8 times greater, orthophosphate 9.7 times greater, total phosphorus 1.2 times greater, and the nitrite + nitrate load was 1.6 times greater. Average yields of dissolved chloride and dissolved sulfate from atmosphoric sources were much less than those transported by streamflow at Blossom Road—chloride was about 1.5 percent and sulfate about 9.1 percent of the amount transported by Irondequoit Creek. Ground water : Ground-water-levels and water quality data were collected from 9 observation wells in Ellison Park in Monroe County. All wells except Mo 2 and Mo 659 are in the flood plain of Irondequoit Creek. Water levels indicate frequent reversals in direction of lateral flow toward or away from Irondequoit Creek, and all wells except Mo2 and Mo 659 respond to water level fluctuations in the Creek. Trend tests on water levels for the period of record indicate a slight upward trend in water levels at all nine wells, two of which (Mo 3 and Mo 667) were statistically significant. Concentrations of ammonia and ammonia + organic nitrogen showed a general decrease for the current period of record. Total phosphorus concentrations showed an increase at four wells and a decrease at four wells. Water quality data showed that the highest median concentrations of nutrients continues to occur in Mo 667 and the highest median concentrations of common ions was at Mo 664. Streamflow : Statistical analysis of long-term (greater than 15 years) streamflow records for unregulated streams in Monroe County indicated that annual mean flows for water years (A water year is the 12-month period from October 1 through September 30 of the following year.) 2000–02 generally were in the normal range (75th to 25th percentile), although Allen Creek continued to show a significant downward trend in mean monthly streamflow during the 1984–2002 water years. Chemical Concentration in Streams : 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 2002. Three of the four Irondequoit Creek sites (Allen Creek, Blossom Road, and Railroad Mills) showed downward trends in ammonia (4.6 to 12.0 percent per year) and ammonia + organic nitrogen (2.8 to 5.3 percent per year). Allen Creek showed downward trends in nitrite + nitrate and total phosphorus (both 1.2 percent per year), and Irondequoit Creek above Blossom Road showed an upward trend in orthophosphate (1.8 percent per year). Three Irondequoit Creek sites showed upward trends in dissolved chloride: Railroad Mills (4.8 percent per year), Allen Creek, and Blossom Road (both 1.9 percent per year). Allen Creek showed a downward trend in sulfate of 0.98 percent per year, whereas Blossom Road showed a downward trend in suspended solids of 4.0 percent per year. Volatile suspended solids showed an upward trend of 3.2 percent per year at Allen Creek and a downward trend of 2.2 percent per year at Blossom Road. Northrup Creek in western Monroe County, showed significant downward trends in concentrations of volatile suspended solids (2.5 percent per year), total phosphorus (5.3 percent per year), and orthophosphate (9.9 percent per year). The Genesee River at Charlotte Docks showed downward trends in volatile suspended solids (2.1 percent per year) and ammonia + organic nitrogen (4.5 percent per year). Oatka Creek at Garbutt showed an upward trend of 21.4 percent per year in turbidity. Chemical Loads in Streams : Mean annual yields (pounds or tons per square mile) of many constituents at the Irondequoit Creek sites were lower than those in previous reporting periods. Suspended solids and nitrite + nitrate yields were lower at three of the sites, and yields of volatile suspended solids, ammonia, and total phosphorus were lower at two of the sites. East Branch Allen Creek showed lower yields for five of the nine constituents for 2000–02, than for previous reporting periods. The decreased yields at East Branch Allen Creek are likely due to the Jefferson Road stormflow-detention basin and the much lower than normal runoff for the 2000–02 period.

New York↗

Simulated effects of alternative withdrawal strategies on groundwater flow in the unconfined Kirkwood-Cohansey aquifer system, the Rio Grande water-bearing zone, and the Atlantic City 800-foot sand in the Great Egg Harbor and Mullica River Basins, New Jersey

Groundwater is essential for water supply and plays a critical role in maintaining the environmental health of freshwater and estuarine ecosystems in the Atlantic Coastal basins of New Jersey. The unconfined Kirkwood-Cohansey aquifer system and the confined Atlantic City 800-foot sand are major sources of groundwater in the area, and each faces different water-supply concerns. The U.S. Geological Survey (USGS), in cooperation with the New Jersey Department of Environmental Protection (NJDEP), conducted a study to simulate the effects of withdrawals in the Kirkwood-Cohansey aquifer system, the Atlantic City 800-foot sand, and the Rio Grande water-bearing zone and to evaluate potential scenarios. The study area encompasses Atlantic County and parts of Burlington, Camden, Gloucester, Ocean, Cape May, and Cumberland Counties. The major hydrogeologic units affecting water supply in the study area are the surficial Kirkwood-Cohansey aquifer system, a thick diatomaceous clay confining unit in the upper part of Kirkwood Formation; the Rio Grande water-bearing zone; and the Atlantic City 800-foot sand of the Kirkwood Formation. Hydrogeologic data from 18 aquifer tests and specific capacity data from 230 wells were analyzed to provide horizontal hydraulic conductivity of the aquifers. Groundwater withdrawals are greatest from the Kirkwood-Cohansey aquifer system, and 65 percent of the water is used for public supply. Groundwater withdrawals from the Atlantic City 800-foot sand are about half those from the Kirkwood-Cohansey aquifer system. Ninety-five percent of the withdrawals from the Atlantic City 800-foot sand is used for public supply. Data from six streamgaging stations and 51 low-flow partial record sites were used to estimate base flow in the area. Base flow ranges from 60 to 92 percent of streamflow. A groundwater flow model of the Kirkwood-Cohansey aquifer system, the Rio Grande water-bearing zone, and the Atlantic City 800-foot sand was developed and calibrated using water-level data from 148 wells and base-flow data from 22 gaging or low-flow partial record stations. The Kirkwood-Cohansey aquifer system within the Great Egg Harbor River and the Mullica River Basins was simulated on a monthly basis from 1998 through 2006. An existing regional model of the New Jersey Coastal Plain was revised to provide boundary conditions for the Great Egg Harbor and Mullica River Basin model (referred to as the Great Egg-Mullica model). In the Great Egg-Mullica model, monthly groundwater recharge rates used in the model ranged from 10-15 inches per year in 2001 to 20-25 inches per year in 2005. The mean-absolute error for 10 of the 14 long-term hydrographs used in model calibration was less than 5 ft. Groundwater flow budgets for the Great Egg-Mullica model calibration periods, May 2005 and September 2006, and for the entire model calibration period 1998 to 2006, showed that nearly 70 percent of the water entering the Atlantic City 800-foot sand came from the horizontal connection with the Kirkwood-Cohansey aquifer system in updip areas. The groundwater flow model was used to simulate scenarios under three possible conditions: average 1998 to 2006 withdrawals (Average scenario), full-allocation withdrawals (Full Allocation scenario), and projected 2050-demand withdrawals (2050 Demand scenario). Withdrawals in the Full Allocation scenario are nearly twice the withdrawals from the Average scenario, primarily because of the potential for large agricultural withdrawals if all allocations are used. Withdrawals for the 2050 Demand scenario are about 50 percent greater than those for the Average scenario, primarily due to expected increases in withdrawals for public supply. Monthly base-flow depletion criteria were determined using the Low-Flow Margin method, currently under consideration by NJDEP, to estimate available water on an annual basis at the Hydrologic Unit Code 11 (HUC11) level and to determine whether a water-supply deficit exists. Simulations of various groundwater-withdrawal scenarios were made using the calibrated model, and results were compared with baseline conditions (no withdrawals) to determine where and when base-flow deficits may be occurring and may be expected to occur in the future. Scenarios were simulated to assess base-flow depletion that could occur from different groundwater-withdrawal situations. In the Average scenario, deficits occurred in 7 of the 14 subbasins. In the Full Allocation scenario, deficits occurred in 11 of the subbasins. In the 2050 Demand scenario, deficits occurred in 9 of the 14 subbasins. The largest deficits occurred in the Absecon Creek subbasin because the base-flow depletion criteria for this subbasin is small due to the surface-water diversions that are already occurring there and because existing groundwater withdrawals in the subbasin have resulted in base-flow depletion under current (1998-2006) conditions. Three adjusted scenarios, variations of the Average, Full Allocation, and 2050 Demand scenarios, were simulated; for the adjusted scenarios, the withdrawals were modified in stages with the intent to successively eliminate or minimize the base-flow deficits. Modifications included shifting withdrawals to a deeper part of the Kirkwood-Cohansey aquifer system, implementing seasonal conjunctive use of shallow and deep aquifers, and specifying reductions in withdrawals within a HUC11 subbasin in deficit. The adjusted scenarios are intended to show the relative effectiveness of each of the three approaches in reducing the deficits. Most of the deficits under the Average, Full Allocation, and 2050 Demand scenarios were eliminated by reductions in withdrawals or allocations. Shifting withdrawals to a deeper part of the Kirkwood-Cohansey aquifer system or seasonal conjunctive use did not eliminate deficits for any subbasin. Reductions in withdrawals accounted for more than 95 percent of the total reduction of deficits in all but one subbasin.

New Jersey↗

Baseline water quality of an area undergoing shale-gas development in the Muskingum River watershed, Ohio, 2015–16

In 2015–16, the U.S. Geological Survey, in cooperation with the Muskingum Watershed Conservancy District, led a study to assess baseline (2015–16) surface-water quality in six lake drainage basins within the Muskingum River watershed that are in the early years of shale-gas development. In 2015, 9 of the 10 most active counties in Ohio for oil and gas development were wholly or partially within the Muskingum River watershed. In addition to shale gas development, the area has a history of conventional oil and gas development and coal mining. In all, 30 surface-water sites were sampled: 20 in tributaries flowing to the lakes, 4 in lakes themselves, and 6 downstream of the lakes. At each of the 30 sites, 6 samples were collected to characterize surface-water chemistry throughout a range of hydrologic conditions. The sampling generally occurred during low flows (periods of greater groundwater contribution) rather than during runoff events (periods of high stream stage). Trilinear diagrams of major ion chemistry revealed three main types of water in the study area―sulfate-dominated waters, bicarbonate-dominated waters, and waters with mixed bicarbonate and chloride anions. Most sites produced samples of bicarbonate-dominated water, and 11 sites produced samples with sulfate-type waters. Mixed bicarbonate and chloride waters were found in samples from two of the six lake drainage basins studied. The baseline (2015–16) assessment of surface-water quality in the study area indicated that few water-chemistry constituents and properties occurred at concentrations or levels that would adversely affect aquatic organisms. Chemical-specific, aquatic life use criteria were not met in only three instances: two were for total dissolved solids at sites likely impacted by coal mining in their drainage basins (hereafter referred to as “mine-impacted sites”), and one was for dissolved oxygen. Mine drainage from historical coal mining in the region likely affected the quality of about one-third of the streams sampled. To simplify interpretation of water-chemistry results, 11 sites with sulfate-type water were identified as mine-impacted sites based on water-quality criteria established by Ohio Department of Natural Resources, Division of Mineral Resources Management, and separated out for subsequent statistical analysis. Concentrations or levels of bicarbonate, boron, calcium, carbonate, total dissolved solids, fluoride, magnesium, lithium, pH, potassium, sodium, specific conductance, strontium, sulfate, and suspended sediment in water were higher (significance level of 0.05) at mine-impacted stream sites than at non-mine-impacted stream sites. An accidental release of oil- and gas-related brines could increase salinity (sodium and chloride), the concentration of total dissolved solids in shallow groundwater and streams, and specific conductance. For this study, chloride concentrations in the study area ranged from 2.12 to 76.1 milligrams per liter. Sources of chloride in water samples were evaluated using binary mixing curves and ratios of chloride to bromide. These ratios indicated that 13 samples from 3 sites in the drainage basin that contained the highest density of conventional oil and gas wells in the study, as well as 4 samples collected from other drainage basins, likely contained a component of brine. Concentrations or levels of barium, bromide, chloride, iron, lithium, manganese, and sodium were significantly higher (alpha = 0.05) in samples with a component of brine than in samples without a component of brine. Benzene, toluene, ethylbenzene and xylene (BTEX), compounds that occur naturally in crude oil, made up 24 of the 45 detections (53 percent) of volatile organic compounds in the study area. The BTEX detections were not associated with sites containing a component of brine. The only volatile organic compound detected in any of the 17 samples that contained a component of brine was acetone, detected in 3 (18 percent) of these samples and in 11 percent of samples not containing a component of brine. Considering that BTEX are gasoline hydrocarbons and that most of the detections occurred during warmer months in and around the lakes, the BTEX detections likely are associated with increases in outdoor activities such as automobile and boating traffic. Radium-226 and radium-228 were included in the list of analytes for this study because production water from shale-gas drilling can contain these naturally occurring radioactive materials. Concentrations of radium-226 exceeded background levels in only two surface-water samples. Concentrations of radium-228 exceeded background levels in one surface-water sample. A brine signature potentially indicative of oil and gas contamination was detected in samples collected at two sites that contained active or plugged waste injection wells, or both. Results from the study indicated significant differences in the median concentrations of bromide, chloride, lithium, manganese, sodium, and total dissolved nitrogen between sites with and without injection wells in their drainage areas. Median concentrations of bromide, chloride, lithium, and sodium, which are common oil- and gas-related contaminants, were higher at sites with injection wells in their drainage areas compared to sites without injection wells. Historical (1960s, 1970s, and 1980s) chloride concentrations and streamflow data at or near five of the six sampling sites downstream from each lake dam were compared to current (2015–16) values. An analysis of covariance was done to test the effects of streamflow, time (decade), and the combined effects (cross product) of streamflow and time on chloride concentrations. Those analyses indicated that streamflow was not significant in explaining the variation in chloride concentration, likely because streamflow in those locations is controlled by dam operations; therefore, association between runoff-generating events and streamflow is less direct than in unregulated streams. From the 1980s to the study period (2015–16), data for three of the five lakes indicated an increase in chloride concentrations. The comparison of historical and current (2015–16) study data from samples collected at another lake indicated that chloride concentrations increased from the 1960s to the 1970s, but concentrations in the 1970s and 2015–16 were similar even though 13 samples from this lake drainage basin were classified as having a component of brine. Median chloride concentrations for the fifth lake, however, seemed to decrease from the 1980s to 2015–16.

Ohio↗

U.S. Geological Survey response to Hurricane Maria flooding in Puerto Rico and characterization of peak streamflows observed September 20–22, 2017

Hurricane Maria struck the island of Puerto Rico on September 20, 2017, as a Category 4 storm. The hurricane traversed the island from southeast to northwest and produced recorded 48-hour rainfall totals of up to 30.01 inches. Estimates of the human death toll range from 2,975 to 4,645, possibly more. The U.S. Geological Survey (USGS) hydrologic monitoring network sustained substantial wind and flood damage during the hurricane. Eighty-five of the 300 hydrologic monitoring stations operating in Puerto Rico and the U.S. Virgin Islands prior to the passage of Hurricane Maria were destroyed or damaged. During the weeks and months after the hurricane, USGS field crews in Puerto Rico prioritized repair of the hydrologic monitoring network and collected hydrologic information to characterize the magnitude of observed peak streamflows at 20 streamgage and to develop new theoretical stage-streamflow relations for 58 streamgages where stream channels were substantially altered; the theoretical stage-streamflow relations were used to estimate Hurricane Maria peak streamflows for 39 of those sites. As part of a pilot program, USGS field crews installed continuous slope-area monitoring equipment at two remote streamgages to automate the collection of high-streamflow stage data. Hurricane Maria peak streamflows and rankings were determined for 73 USGS streamgages in Puerto Rico. New rank 1 period-of-record peak streamflows occurred at 28 sites, rank 2 period-of-record peak streamflows occurred at 17 sites, and rank 3 period-of-record peak streamflows occurred at 9 sites; period-of-record peak streamflows at the remaining 19 sites either ranked from 4th to 20th or were not ranked. Annual exceedance probabilities for 53 unregulated peak streamflows ranged from greater than 50.0 percent (recurrence interval of less than 2 years) to 0.3 percent (recurrence interval of 333 years), with the majority (28 of 53) in the range of 10.0–2.1 percent (recurrence intervals of 10–48 years). A comparison of period-of-record ranks for the largest flood events that have occurred in Puerto Rico since the 1960s indicated that Hurricane Maria produced more record peak streamflows than either Hurricane Hortense in 1996 or Hurricane Georges in 1998. Limited pre-1960s hydrologic data preclude quantitative comparison with earlier storms. As part of this study, a maximum peak-streamflow envelope curve for Puerto Rico was developed using historical peak-streamflow information available through 2017. Other post-Hurricane Maria USGS activities summarized in this report include (1) Global Navigation Satellite System surveys at all stations in the USGS hydrologic monitoring network, used to tie the network to the Puerto Rico Vertical Datum of 2002; and (2) telemetered monitoring of the Lago Guajataca Dam in northwestern Puerto Rico, which was damaged and at risk of failure from October to December 2017.

Puerto Rico↗

Aquatic biological communities and associated habitats at selected sites in the Big Wood River Watershed, south-central Idaho, 2014

Assessments of streamflow (discharge) parameters, water quality, physical habitat, and biological communities were completed between May and September 2014 as part of a monitoring program in the Big Wood River watershed of south-central Idaho. The sampling was conducted by the U.S. Geological Survey in cooperation with Blaine County, Trout Unlimited, the Nature Conservancy, and the Wood River Land Trust to help identify the status of aquatic resources at selected locations in the watershed. Information in this report provides a basis with which to evaluate and monitor the long-term health of the Big Wood River and its major tributaries. Sampling sites were co-located with existing U.S. Geological Survey streamgaging stations: three on the main stem Big Wood River and four on the North Fork Big Wood River (North Fork), Warm Springs Creek (Warm Sp), Trail Creek (Trail Ck), and East Fork Big Wood River (East Fork) tributaries. The analytical results and quality-assurance information for water quality, physical habitat, and biological community samples collected at study sites during 2 weeks in September 2014 are summarized. Water-quality data include concentrations of major nutrients, suspended sediment, dissolved oxygen, and fecal-coliform bacteria. To assess the potential effects of nutrient enrichment on algal growth, concentrations of periphyton biomass (chlorophyll- a and ash free dry weight) in riffle habitats were determined at each site. Physical habitat parameters include stream channel morphology, habitat volume, instream structure, substrate composition, and riparian vegetative cover. Biological data include taxa richness, abundance, and stream-health indicator metrics for macroinvertebrate and fish communities. Statistical summaries of the water-quality, habitat, and biological data are provided along with discussion of how these findings relate to the health of aquatic resources in the Big Wood River watershed. Seasonal discharge patterns using statistical summaries of daily discharge from selected sites are reported for water years 2012–15. Results showed that annual average daily mean discharge increased from the Big Wood River near Ketchum, ID (BW Ketchum) downstream to the Big Wood River at Hailey, ID (BW Hailey), but decreased by nearly 50 percent from BW Hailey downstream to Big Wood River at Stanton Crossing near Bellevue, ID (BW Stanton). Annual variability in daily mean discharge among main-stem sites was highest at BW Stanton, suggesting that this part of the river may be subject to some level of flow alteration. Hydrologic alterations resulting in flow reduction can contribute to higher water temperature, especially during the summer months when conditions are often most stressful to fish and other stream organisms. Daily water temperature and water temperature trends from June to September 2014 are reported for select tributary and main-stem sites on the Big Wood River and can be used to assess the potential for biological impairment based on aquatic life temperature criteria for cold-water streams. The State of Idaho maximum temperature criteria for protection of cold-water aquatic life of 22 °C was exceeded at Warm Sp and BW Stanton during summer 2014, but at none of the other main-stem or tributary sites. The 13 °C critical temperature criterion for salmonid spawning was exceeded in early July 2014 at BW Ketchum and BW Hailey near the end of the rainbow trout critical spawning and rearing period. Temperature exceedances were most frequent at BW Stanton, where exceedances for rainbow trout and brown trout occurred from May through early July 2014 during most of the critical spawning and rearing period. Water quality and habitat availability did not seem to be limiting for fish or other aquatic organisms at most sites in the Big Wood River watershed. Water quality assessments in September 2014 determined no exceedances of total maximum daily load target levels. The availability and quality of habitat was limited at BW Stanton, where shallow-water habitat conditions prevailed. Macroinvertebrate community diversity was high at all sites except for BW Stanton, where low community diversity was attributed to low species richness and high abundances of a few tolerant taxa. Presence of low species diversity and high macroinvertebrate tolerance values at BW Stanton indicates that benthic community condition and stream health were reduced at that location. Fish surveys done in September 2014 did not indicate any significant reductions in native fish communities in the Big Wood River or its tributaries. Native rainbow trout ( Oncorhynchus mykiss ) and Wood River sculpin ( Cottus leiopomus ) were the dominant fish species in the drainage and were found at all tributary and main-stem sites. Non-native brown ( Salmo trutta ) and brook trout ( Salvelinus fontinalis ) were limited to lower drainage sites on the Big Wood River (BW Hailey and BW Stanton), and occurred in relatively low numbers.

Idaho↗

Surface-water quantity and quality, aquatic biology, stream geomorphology, and groundwater-flow simulation for National Guard Training Center at Fort Indiantown Gap, Pennsylvania, 2002-05

Base-line and long-term monitoring of water resources of the National Guard Training Center at Fort Indiantown Gap in south-central Pennsylvania began in 2002. Results of continuous monitoring of streamflow and turbidity and monthly and stormflow water-quality samples from two continuous-record long-term stream sites, periodic collection of water-quality samples from five miscellaneous stream sites, and annual collection of biological data from 2002 to 2005 at 27 sites are discussed. In addition, results from a stream-geomorphic analysis and classification and a regional groundwater-flow model are included. Streamflow at the facility was above normal for the 2003 through 2005 water years and extremely high-flow events occurred in 2003 and in 2004. Water-quality samples were analyzed for nutrients, sediments, metals, major ions, pesticides, volatile and semi-volatile organic compounds, and explosives. Results indicated no exceedances for any constituent (except iron) above the primary and secondary drinking-water standards or health-advisory levels set by the U.S. Environmental Protection Agency. Iron concentrations were naturally elevated in the groundwater within the watershed because of bedrock lithology. The majority of the constituents were at or below the method detection limit. Sediment loads were dominated by precipitation due to the remnants of Hurricane Ivan in September 2004. More than 60 percent of the sediment load measured during the entire study was transported past the streamgage in just 2 days during that event. Habitat and aquatic-invertebrate data were collected in the summers of 2002-05, and fish data were collected in 2004. Although 2002 was a drought year, 2003-05 were above-normal flow years. Results indicated a wide diversity in invertebrates, good numbers of taxa (distinct organisms), and on the basis of a combination of metrics, the majority of the 27 sites indicated no or slight impairment. Fish-metric data from 25 sites indicated results similar to the invertebrate data. Stream classification based on evolution of the stream channels indicates about 94 percent of the channels were considered to be in equilibrium (type B or C channels), neither aggrading nor eroding. A regional, uncalibrated groundwater-flow model indicated the surface-water and groundwater-flow divides coincided. Because of folding of rock layers, groundwater was under confined conditions and nearly all the water leaves the facility via the streams.

Scientific Investigations Report↗

Total nitrogen and suspended-sediment loads and identification of suspended-sediment sources in the Laurel Hill Creek watershed, Somerset County, Pennsylvania, water years 2010-11

Laurel Hill Creek is a watershed of 125 square miles located mostly in Somerset County, Pennsylvania, with small areas extending into Fayette and Westmoreland Counties. The upper part of the watershed is on the Pennsylvania Department of Environmental Protection 303(d) list of impaired streams because of siltation, nutrients, and low dissolved oxygen concentrations. The objectives of this study were to (1) estimate the annual sediment load, (2) estimate the annual nitrogen load, and (3) identify the major sources of fine-grained sediment using the sediment-fingerprinting approach. This study by the U.S. Geological Survey (USGS) was done in cooperation with the Somerset County Conservation District. Discharge, suspended-sediment, and nutrient data were collected at two streamflow-gaging stations—Laurel Hill Creek near Bakersville, Pa., (station 03079600) and Laurel Hill Creek at Ursina, Pa., (station 03080000)—and one ungaged stream site, Laurel Hill Creek below Laurel Hill Creek Lake at Trent (station 03079655). Concentrations of nutrients generally were low. Concentrations of ammonia were less than 0.2 milligrams per liter (mg/L), and concentrations of phosphorus were less than 0.3 mg/L. Most concentrations of phosphorus were less than the detection limit of 0.02 mg/L. Most water samples had concentrations of nitrate plus nitrite less than 1.0 mg/L. At the Bakersville station, concentrations of total nitrogen ranged from 0.63 to 1.3 mg/L in base-flow samples and from 0.57 to 1.5 mg/L in storm composite samples. Median concentrations were 0.88 mg/L in base-flow samples and 1.2 mg/L in storm composite samples. At the Ursina station, concentrations of total nitrogen ranged from 0.25 to 0.92 mg/L in base-flow samples; the median concentration was 0.57 mg/L. The estimated total nitrogen load at the Bakersville station was 262 pounds (lb) for 11 months of the 2010 water year (November 2009 to September 2010) and 266 lb for the 2011 water year. Most of the total nitrogen loading was from stormflows. The stormflow load accounted for 76.6 percent of the total load for the 2010 water year and 80.6 percent of the total load for the 2011 water year. The estimated monthly total nitrogen loads were higher during the winter and spring (December through May) than during the summer (June through August). For the Bakersville station, the estimated suspended-sediment load (SSL) was 17,700 tons for 11 months of the 2010 water year (November 2009 to September 2010). The storm beginning January 24, 2010, provided 34.4 percent of the annual SSL, and the storm beginning March 10, 2010, provided 31.9 percent of the annual SSL. Together, these two winter storms provided 66 percent of the annual SSL for the 2010 water year. For the 2011 water year, the estimated annual SSL was 13,500 tons. For the 2011 water year, the SSLs were more evenly divided among storms than for the 2010 water year. Seven of 37 storms with the highest SSLs provided a total of 65.7 percent of the annual SSL for the 2011 water year; each storm provided from 4.6 to 12.3 percent of the annual SSL. The highest cumulative SSL for the 2010 and 2011 water years generally occurred during the late winter. Stormflows with the highest peak discharges generally carried the highest SSL. The sediment-fingerprinting approach was used to quantify sources of fine-grained suspended sediment in the watershed draining to the Laurel Hill Creek near Bakersville streamflow-gaging station. Sediment source samples were collected from five source types: 20 from cropland, 9 from pasture, 18 from forested areas, 20 from unpaved roads, and 23 from streambanks. At the Bakersville station, 10 suspended-sediment samples were collected during 6 storms for sediment-source analysis. Thirty-five tracers from elemental analysis and 4 tracers from stable isotope analysis were used to fingerprint the source of sediment for the 10 storm samples. Statistical analysis determined that cropland and pasture could not be discriminated by the set of tracers and were combined into one source group—agriculture. Stepwise discriminant function analysis determined that 11 tracers best described the 4 sources. An "unmixing" model applied to the 11 tracers showed that agricultural land (cropland and pasture) was the major source of sediment, contributing an average of 53 percent of the sediment for the 10 storm samples. Streambanks, unpaved roads, and forest contributions for the 10 storm samples averaged 30, 17, and 0 percent, respectively. Agriculture was the major contributor of sediment during the highest sampled stormflows. The highest stormflows also produced the highest total nitrogen and suspended-sediment loads.

Pennsylvania↗

Groundwater simulation and management models for the upper Klamath Basin, Oregon and California

The upper Klamath Basin encompasses about 8,000 square miles, extending from the Cascade Range east to the Basin and Range geologic province in south-central Oregon and northern California. The geography of the basin is dominated by forested volcanic uplands separated by broad interior basins. Most of the interior basins once held broad shallow lakes and extensive wetlands, but most of these areas have been drained or otherwise modified and are now cultivated. Major parts of the interior basins are managed as wildlife refuges, primarily for migratory waterfowl. The permeable volcanic bedrock of the upper Klamath Basin hosts a substantial regional groundwater system that provides much of the flow to major streams and lakes that, in turn, provide water for wildlife habitat and are the principal source of irrigation water for the basin's agricultural economy. Increased allocation of surface water for endangered species in the past decade has resulted in increased groundwater pumping and growing interest in the use of groundwater for irrigation. The potential effects of increased groundwater pumping on groundwater levels and discharge to springs and streams has caused concern among groundwater users, wildlife and Tribal interests, and State and Federal resource managers. To provide information on the potential impacts of increased groundwater development and to aid in the development of a groundwater management strategy, the U.S. Geological Survey, in collaboration with the Oregon Water Resources Department and the Bureau of Reclamation, has developed a groundwater model that can simulate the response of the hydrologic system to these new stresses. The groundwater model was developed using the U.S. Geological Survey MODFLOW finite-difference modeling code and calibrated using inverse methods to transient conditions from 1989 through 2004 with quarterly stress periods. Groundwater recharge and agricultural and municipal pumping are specified for each stress period. All major streams and most major tributaries for which a substantial part of the flow comes from groundwater discharge are included in the model. Groundwater discharge to agricultural drains, evapotranspiration from aquifers in areas of shallow groundwater, and groundwater flow to and from adjacent basins also are simulated in key areas. The model has the capability to calculate the effects of pumping and other external stresses on groundwater levels, discharge to streams, and other boundary fluxes, such as discharge to drains. Historical data indicate that the groundwater system in the upper Klamath Basin fluctuates in response to decadal climate cycles, with groundwater levels and spring flows rising and declining in response to wet and dry periods. Data also show that groundwater levels fluctuate seasonally and interannually in response to groundwater pumping. The most prominent response is to the marked increase in groundwater pumping starting in 2001. The calibrated model is able to simulate observed decadal-scale climate-driven fluctuations in the groundwater system as well as observed shorter-term pumping-related fluctuations. Example model simulations show that the timing and location of the effects of groundwater pumping vary markedly depending on the pumping location. Pumping from wells close (within a few miles) to groundwater discharge features, such as springs, drains, and certain streams, can affect those features within weeks or months of the onset of pumping, and the impacts can be essentially fully manifested in several years. Simulations indicate that seasonal variations in pumping rates are buffered by the groundwater system, and peak impacts are closer to mean annual pumping rates than to instantaneous rates. Thus, pumping effects are, to a large degree, spread out over the entire year. When pumping locations are distant (more than several miles) from discharge features, the effects take many years or decades to fully impact those features, and much of the pumped water comes from groundwater storage over a broad geographic area even after two decades. Moreover, because the effects are spread out over a broad area, the impacts to individual features are much smaller than in the case of nearby pumping. Simulations show that the discharge features most affected by pumping in the area of the Bureau of Reclamation's Klamath Irrigation Project are agricultural drains, and impacts to other surface-water features are small in comparison. A groundwater management model was developed that uses techniques of constrained optimization along with the groundwater flow model to identify the optimal strategy to meet water user needs while not violating defined constraints on impacts to groundwater levels and streamflows. The coupled groundwater simulation-optimization models were formulated to help identify strategies to meet water demand in the upper Klamath Basin. The models maximize groundwater pumping while simultaneously keeping the detrimental impacts of pumping on groundwater levels and groundwater discharge within prescribed limits. Total groundwater withdrawals were calculated under alternative constraints for drawdown, reductions in groundwater discharge to surface water, and water demand to understand the potential benefits and limitations for groundwater development in the upper Klamath Basin. The simulation-optimization model for the upper Klamath Basin provides an improved understanding of how the groundwater and surface-water system responds to sustained groundwater pumping within the Bureau of Reclamation's Klamath Project. Optimization model results demonstrate that a certain amount of supplemental groundwater pumping can occur without exceeding defined limits on drawdown and stream capture. The results of the different applications of the model demonstrate the importance of identifying constraint limits in order to better define the amount and distribution of groundwater withdrawal that is sustainable.

Oregon;California↗

Maps showing estimated sediment yield from coastal landslides and active slope distribution along the Big Sur coast, Monterey and San Luis Obispo Counties, California

The 1982-83 and 1997-98 El Ni?os brought very high precipitation to California?s central coast; this precipitation resulted in raised groundwater levels, coastal flooding, and destabilized slopes throughout the region. Large landslides in the coastal mountains of Big Sur in Monterey and San Luis Obispo Counties blocked sections of California State Route 1, closing the road for months at a time. Large landslides such as these occur frequently in the winter months along the Big Sur coast due to the steep topography and weak bedrock. A large landslide in 1983 resulted in the closure of Highway 1 for over a year to repair the road and stabilize the slope. Resulting work from the 1983 landslide cost over $7 million and generated 30 million cubic yards of debris from landslide removal and excavations to re-establish the highway along the Big Sur coast. Before establishment of the Monterey Bay National Marine Sanctuary (MBNMS) in 1992, typical road opening measures involved disposal of some landslide material and excess material generated from slope stabilization onto the seaward side of the highway. It is likely that some or most of this disposed material, either directly or indirectly through subsequent erosion, was eventually transported downslope into the ocean. In addition to the landslides that initiate above the road, natural slope failures sometimes occur on the steep slopes below the road and thus deliver material to the base of the coastal mountains where it is eroded and dispersed by waves and nearshore currents. Any coastal-slope landslide, generated through natural or anthropogenic processes, can result in sediment entering the nearshore zone. The waters offshore of the Big Sur coast are part of the MBNMS. Since it was established in 1992, landslide-disposal practices came under question for two reasons. The U.S. Code of Federal Regulations, Title 15, Section 922.132 prohibits discharging or depositing, from beyond the boundary of the Sanctuary, any material or other matter that subsequently enters the Sanctuary and injures a Sanctuary resource or quality. The landslide-disposal practices previously used had the potential to alter nearshore zone habitat by converting marine habitats from rocky substrate to soft bottom. In addition, the disposal practices had the potential to increase concentrations of suspended sediment in the nearshore zone, and thereby possibly impact coastal biological communities. On the other hand, natural mass-wasting processes, including coastal cliff erosion coastal and landslides, as well as streams deliver sediment to the coast in unknown quantities and thus provide nutrients as well as source material for beaches. Current road maintenance and repair practices along the Coast Highway 1 corridor that restrict disposal of material within the MBNMS may actually reduce sediment input relative to natural processes.

Scientific Investigations Map↗

Geophysical and hydrologic analysis of an earthen dam site in southern Westchester County, New York

Ninety percent of the drinking water for New York City passes through the Hillview Reservoir facility in the City of Yonkers, Westchester County, New York. In the past, several seeps located downslope from the reservoir have flowed out from the side of the steepest slope at the southern end of the earthen embankment. One seep that has been flowing continuously was discovered during an inspection of the embankment in 1999. Efforts were made in 2001 to locate the potential sources of the continuous flowing seep. In 2005, the U.S. Geological Survey, in cooperation with the New York City Department of Environmental Protection, began a cooperative study to investigate the relevant hydrogeologic framework to characterize the local groundwater-flow system and to determine possible sources of the seeps. The two agencies used hydrologic and surface geophysical techniques to assess the earthen embankment of the Hillview Reservoir. Between April 1, 2005 and March 1, 2008, water levels were measured manually each month at 46 wells surrounding the reservoir, and flow was measured monthly at three of the five seeps on the embankment. Water levels were measured hourly in the East Basin of the reservoir, at 24 of 46 wells, and discharge was measured hourly at two of the five seeps. Slug tests were performed at 16 wells to determine the hydraulic conductivity of the geologic material surrounding the screened zone. Estimated hydraulic conductivities for 25 wells on the southern embankment ranged from 0.0063 to 1.2 feet per day and averaged 0.17 foot per day. The two-dimensional resistivity surveys indicate a subsurface mound of electrically conductive material (low-resistivity zone) beneath the terrace area (top of dam) surrounding the reservoir with a distinct elevation increase closer to the crest. Two-dimensional shear wave velocity surveys indicate a similar structure of the high shear wave velocity materials (high-velocity zone), increasing in elevation toward the crest and decreasing toward the reservoir and toward the northern part of the study area. Water-quality samples collected from 12 wells, downtake chamber 1 of the reservoir, and two seeps detected the presence of arsenic, toluene, and two trihalomethanes. Water-quality samples collected at the two seeps detected fluoride, indicating a connection with reservoir water. Shallow wells on the southern embankment exhibited the largest seasonal water-level fluctuations ranging between 6 feet and 12 feet. The embankment is constructed from reworked low-permeability glacial deposits at the site. Water-level responses in observation wells within the embankment indicate that there is a shallow (approximately the upper 45 feet of the embankment) and a deep water-bearing unit within the embankment with a large downward vertical gradient between the shallow and deep water-bearing units. Precipitation strongly affected water levels in shallow wells, whereas the basin appears to be the main control on water levels in the deep wells. Seeps on the embankment slope appear to be caused by above-average precipitation that increases water levels in the shallow water-bearing unit, but does not easily recharge the deep water-bearing unit. Based on the data that have been analyzed, source water to the seeps appears to be primarily groundwater and, to a lesser extent, water from the East Basin of the reservoir.

New York↗

Simulation of groundwater and surface-water flow in the upper Deschutes Basin, Oregon

This report describes a hydrologic model for the upper Deschutes Basin in central Oregon developed using the U.S. Geological Survey (USGS) integrated Groundwater and Surface-Water Flow model (GSFLOW). The upper Deschutes Basin, which drains much of the eastern side of the Cascade Range in Oregon, is underlain by large areas of permeable volcanic rock. That permeability, in combination with the large annual precipitation at high elevations, results in a substantial regional aquifer system and a stream system that is heavily groundwater dominated. The upper Deschutes Basin is also an area of expanding population and increasing water demand for public supply and agriculture. Surface water was largely developed for agricultural use by the mid-20th century, and is closed to additional appropriations. Consequently, water users look to groundwater to satisfy the growing demand. The well‑documented connection between groundwater and the stream system, and the institutional and legal restrictions on streamflow depletion by wells, resulted in the Oregon Water Resources Department (OWRD) instituting a process whereby additional groundwater pumping can be permitted only if the effects to streams are mitigated, for example, by reducing permitted surface-water diversions. Implementing such a program requires understanding of the spatial and temporal distribution of effects to streams from groundwater pumping. A groundwater model developed in the early 2000s by the USGS and OWRD has been used to provide insights into the distribution of streamflow depletion by wells, but lacks spatial resolution in sensitive headwaters and spring areas. The integrated model developed for this project, based largely on the earlier model, has a much finer grid spacing allowing resolution of sensitive headwater streams and important spring areas, and simulates a more complete set of surface processes as well as runoff and groundwater flow. In addition, the integrated model includes improved representation of subsurface geology and explicitly simulates the effects of hydrologically important fault zones not included in the previous model. The upper Deschutes Basin GSFLOW model was calibrated using an iterative trial and error approach using measured water-level elevations (water levels) from 800 wells, 144 of which have time series of 10 or more measurements. Streamflow was calibrated using data from 21 gage locations. At 14 locations where measured flows are heavily influenced by reservoir operations and irrigation diversions, so called “ naturalized ” flows, with the effects of reservoirs and diversion removed, developed by the Bureau of Reclamation, were used for calibration. Surface energy and moisture processes such as solar radiation, snow accumulation and melting, and evapotranspiration were calibrated using national datasets as well as data from long-term measurement sites in the basin. The calibrated Deschutes GSFLOW model requires daily precipitation, minimum and maximum air temperature data, and monthly data describing groundwater pumping and artificial recharge from leaking irrigation canals (which are a significant source of groundwater recharge). The calibrated model simulates the geographic distribution of hydraulic head over the 5,000 ft range measured in the basin, with a median absolute residual of about 53 ft. Temporal variations in head resulting from climate cycles, pumping, and canal leakage are well simulated over the model area. Simulated daily streamflow matches gaged flows or calculated naturalized flows for streams including the Crooked and Metolius Rivers, and lower parts of the mainstem Deschutes River. Seasonal patterns of runoff are less well fit in some upper basin streams. Annual water balances of streamflow are good over most of the model domain. Model fit and overall capabilities are appropriate for the objectives of the project. The integrated model results confirm findings from other studies and models indicating that most streamflow in the upper Deschutes Basin comes directly from groundwater discharge. The integrated model provides additional insights about the components of streamflow including direct groundwater discharge to streams, interflow, groundwater discharge to the land surface (Dunnian flow), and direct runoff (Hortonian flow). The new model provides improved capability for exploring the timing and distribution of streamflow capture by wells, and the hydrologic response to changes in other external stresses such as canal operation, irrigation, and drought. Because the model uses basic meteorological data as the primary input; and simulates surface energy and moisture balances, groundwater recharge and flow, and all components of streamflow; it is well suited for exploring the hydrologic response to climate change, although no such simulations are included in this report. The model was developed as a tool for future application; however, example simulations are provided in this report. In the example simulations, the model is used to explore the influence of well location and geologic structure on stream capture by pumping wells. Wells were simulated at three locations within a 12-mi area close to known groundwater discharge areas and crossed by a regional fault zone. Simulations indicate that the magnitude and timing of stream capture from pumping is largely controlled by the geographic location of the wells, but that faults can have a large influence on the propagation of pumping stresses.

Oregon↗

Relation of hydrologic processes to groundwater and surface-water levels and flow directions in a dune-beach complex at Indiana Dunes National Lakeshore and Beverly Shores, Indiana

The potential for high groundwater levels to cause wet basements (groundwater flooding) is of concern to residents of communities in northwestern Indiana. Changes in recharge from precipitation increases during 2006-9, water-level changes from restoration of nearby wetlands in the Great Marsh in 1998-2002, and changes in recharge due to the end of groundwater withdrawals for water supply since 2005 in a community at Beverly Shores, Ind., were suspected as factors in increased groundwater levels in an unconfined surficial aquifer beneath nearby parts of a dune-beach complex. Results of this study indicate that increased recharge from precipitation and snowmelt was the principal cause of raised water levels in the dune-beach complex from 2006 to 2009. Annual precipitation totals in 2006-9 ranged from 43.88 to 55.75 inches per year (in/yr) and were substantially greater than the median 1952-2009 precipitation of 36.35 in/yr. Recharge to groundwater from precipitation in 2006-9 ranged from 13.5 to 22 in/yr; it was higher than the typical 11 in/yr because of large precipitation events and precipitation amounts received during non-growing-season months. An estimated increase in net recharge from reduced groundwater use in Beverly Shores since 2005 ranged from 1.6 in/yr in 2006 to 1.9 in/yr in 2009. Surface-water levels in the wetland were as much as about 1.1 feet higher in 2007-9 (after the 1998-2002 wetland restoration) than during seasonally wet periods in 1979-89. Similar surface-water levels and ponded water were likely during winter and spring wet periods before and after wetland restoration. High water levels similar to those in 2009 were measured elsewhere in the dune-beach complex near a natural wetland during the spring months in 1991 and 1993 after receipt of near record precipitation. Recharge from similarly high precipitation amounts in 2008-9 was also a likely cause of high groundwater levels in other parts of the dune-beach complex, such as at Beverly Shores. Perennial mounding of the water table in the surficial aquifer indicates that the recharge that created the water-table mound originates within the dune-beach complex and not through flow from the adjacent hydrologic boundaries: the restored wetland, Lake Michigan, and Derby Ditch. Infiltrating precipitation causes most seasonal and episodic rises in groundwater levels beneath the dune-beach complex. Groundwater-level fluctuations lasting days to weeks in the dune-beach complex in 2008-9 were superimposed on a seasonal high water-table altitude that began with the recharge from snowmelt and rain in February 2009 and maintained through July 2009. Increases in water-table-mound altitude under the dune-beach complex recurred in 2008-9 in response to the largest rain events of 1 inch or more and to snowmelt. Smaller, shorter-term rises in water level after individual rain events persisted over hours to less than 1 week. Groundwater-level fluctuations varied over a relatively narrow range of about 2 to 3 feet, with no net fluctuations greater than 4 feet. Groundwater levels in or near low parts of the dune-beach complex were frequently within 0 to 6 feet of the land surface and indicate the potential for groundwater flooding. Groundwater-level gradients from the water-table mound to wells next to surface-water discharges increase after rainfall and snowmelt events and recede slowly as groundwater discharges from the aquifer. Evapotranspiration is responsible for part of the general pattern of decreasing water-table altitudes observed from May to August 2009. Rapid water-level rises in the restored wetland after precipitation do not likely have an effect on groundwater flooding elsewhere in the dune-beach complex. Surface-water-level fluctuations during this study generally varied over a narrower range, approximately from 1 to 1.5 feet, as compared with groundwater fluctuations, except after a very large, 10.77-inch rainfall. Time-delayed and smaller groundwater-level

Indiana↗

Comparison of water consumption in two riparian vegetation communities along the central Platte River, Nebraska, 2008–09 and 2011

The Platte River is a vital natural resource for the people, plants, and animals of Nebraska. A recent study quantified water use by riparian woodlands along central reaches of the Platte River, Nebraska, finding that water use was mainly regulated below maximum predicted levels. A comparative study was launched through a cooperative partnership between the U.S. Geological Survey, the Central Platte Natural Resources District, the Nebraska Department of Natural Resources, and the Nebraska Environmental Trust to compare water use of a riparian woodland with that of a grazed riparian grassland along the central Platte River. This report describes the results of the 3-year study by the U.S. Geological Survey to measure the evapotranspiration (ET) rates in the two riparian vegetation communities. Evapotranspiration was measured during 2008–09 and 2011 using the eddy-covariance method at a riparian woodland near Odessa, hereinafter referred to as the “woodland site,” and a riparian grassland pasture near Elm Creek, hereinafter referred to as the “grassland site.” Overall, annual ET totals at the grassland site were 90 percent of the annual ET measured at the woodland site, with averages of 653 millimeters (mm) and 726 mm, respectively. Evapotranspiration rates were similar at the grassland site and the woodland site during the spring and fall seasons, but at the woodland site ET rates were higher than those of the grassland site during the peak-growth summer months of June through August. These seasonal differences and the slightly lower ET rates at the grassland site were likely the result of differing plant communities, disturbance effects related to grazing and flooding, and climatic differences between the sites. The annual water balance was calculated for each site and indicated that the predominant factors in the water balance at both sites were ET and precipitation. Annual precipitation for the study period ranged from near to above the normal precipitation of 640 mm. Substantial precipitation fell in May and October 2008 that caused flooding along the Platte River in May of this especially wet year. There was a deficit in precipitation compared to ET at both sites in 2009 and 2011, leading to a net groundwater use of greater than 140 mm per year at the woodland site and greater than 55 mm per year at the grassland site. This indicates that the net annual groundwater use or recharge depends predominately upon the relation between ET and precipitation in these riparian areas with shallow soil layers above the groundwater table. Prior research at the woodland site provided four additional annual water balances dating back to 2002 for comparison with the study period at the woodland site. Perhaps most striking in this comparison was the 25-percent increase in annual ET for 2008–09 and 2011 despite precipitation totals and potential ET rates that were within the range of those measured in 2002–05. As a result, the water balance indicates that groundwater was discharged 2 of the 3 years of the study. This likely was caused by higher groundwater levels and a healthier plant community in 2008–09 and 2011 relative to the drought-affected years of 2002–05. As a result of these changes, the crop coefficients developed for riparian woodlands during the prior research underestimated 2008–09 and 2011 annual ET rates by an average of 35 percent. Though new crop coefficients were developed by this study, the importance of soil-moisture stress and plant community successional dynamics need to be considered when applying these coefficients at other riparian sites or into the future. Nonetheless, their development and the data on which they are based may provide improved understanding of water consumption by riparian grasslands and riparian woodlands along the central Platte River.

Nebraska↗