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 739 records · Page 41Linked to original sources

Groundwater levels, geochemistry, and water budget of the Tsala Apopka Lake system, west-central Florida, 2004–12

Tsala Apopka Lake is a complex system of lakes and wetlands, with intervening uplands, located in Citrus County in west-central Florida. It is located within the 2,100 square mile watershed of the Withlacoochee River, which drains north and northwest towards the Gulf of Mexico. The lake system is managed by the Southwest Florida Water Management District as three distinct “pools,” which from upstream to downstream are referred to as the Floral City Pool, Inverness Pool, and Hernando Pool. Each pool contains a mixture of deep-water lakes that remain wet year round, ephemeral (seasonal) ponds and wetlands, and dry uplands. Many of the major deep-water lakes are interconnected by canals. Flow from the Withlacoochee River, when conditions allow, can be diverted into the lake system. Flow thorough the canals can be used to control the distribution of water between the three pools. Flow in the canals is controlled using structures, such as gates and weirs. Hydrogeologic units in the study area include a surficial aquifer consisting of Quaternary-age sediments, a discontinuous intermediate confining unit consisting of Miocene- and Pliocene-age sediments, and the underlying Upper Floridan aquifer, which consists of Eocene- and Oligocene-age carbonates. The fine-grained quartz sands that constitute the surficial aquifer are generally thin, typically less than 25 feet thick, within the vicinity of Tsala Apopka Lake. A thin, discontinuous, sandy clay layer forms the intermediate confining unit. The Upper Floridan aquifer is generally unconfined in the vicinity of Tsala Apopka Lake because the intermediate confining unit is discontinuous and breached by numerous karst features. In the study area, the Upper Floridan aquifer includes the upper Avon Park Formation and Ocala Limestone. The Ocala Limestone is the primary source of drinking water and spring flow in the area. The objectives of this study are to document the interaction of Tsala Apopka Lake, the surficial aquifer, and the Upper Floridan aquifer; and to estimate an annual water budget for each pool and for the entire lake system for 2004–12. The hydrologic interactions were evaluated using hydraulic head and geochemical data. Geochemical data, including major ion, isotope, and age-tracer data, were used to evaluate sources of water and to distinguish flow paths. Hydrologic connection of the surficial environment (lakes, ponds, wetlands, and the surficial aquifer) was quantified on the basis of a conceptualized annual water-budget model. The model included the change in surface water and groundwater storage, precipitation, evapotranspiration, surface-water inflow and outflow, and net groundwater exchange with the underlying Upper Floridan aquifer. The control volume for each pool extended to the base of the surficial aquifer and covered an area defined to exceed the maximum inundated area for each pool during 2004–12 by 0.5 foot. Net groundwater flow was computed as a lumped value and was either positive or negative, with a negative value indicating downward or lateral leakage from the control volume and a positive value indicating upward leakage to the control volume. The annual water budget for Tsala Apopka Lake was calculated using a combination of field observations and remotely sensed data for each of three pools and for the composite three pool area. A digital elevation model at a 5-foot grid spacing and bathymetric survey data were used to define the land-surface elevation and volume of each pool and to calculate the changes in inundated area with change in lake stage. Continuous lake-stage and groundwater-level data were used to define the change in storage for each pool. The rainfall data used in the water-budget calculations were based on daily radar reflectance data and measured rainfall from weather stations. Evapotranspiration was computed as a function of reference evapotranspiration, adjusted to actual evapotranspiration using a monthly land-cover coefficient (based on evapotranspiration measurements at stations located in representative landscapes). Surface-water inflows and outflows were determined using stage data collected at a series of streamgages installed primarily at the water-control structures. Discharge was measured under varying flow regimes and ratings were developed for the water-control structures. The discharge data collected during the study period were used to calibrate a surface-water flow model for 2004–12. Flows predicted by the model were used in the water-budget analysis. Net groundwater flow was determined as the residual term in the water-budget equation. The results of the water-budget analysis indicate that rainfall was the largest input of water to Tsala Apopka Lake, whereas evapotranspiration was the largest output. For the 2004–12 analysis period, surface-water inflow accounted for 11 percent of the inputs, net groundwater inflow accounted for 1 percent of inputs (annual periods with positive net groundwater flow were included as inputs, while annual periods with negative net groundwater flow were counted as outputs), and rainfall accounted for the remaining 88 percent. For the same period, the outputs consisted of 2 percent surface-water outflow, 12 percent net groundwater outflow, and 86 percent evapotranspiration. Net groundwater inflows and surface-water/groundwater storage were negligible during the water-budget period but could be important components of the budget in individual years. The net groundwater flow was negative (downward) for 8 out of the 9 years modeled (2004–12), indicating that the Tsala Apopka Lake study area was primarily a recharge area for the underlying Upper Floridan aquifer during this time period. Groundwater-level elevation in paired wells (adjacent wells completed in the surficial aquifer and Upper Floridan aquifer) typically was higher in the surficial aquifer than the Upper Floridan aquifer. However, hydraulic head data indicate that the surficial aquifer often has discharge potential to the surface-water system, especially in the low lying areas near the major lakes. Surficial-aquifer water levels were often higher than lake stages, especially during wet periods, which is likely an indication of aquifer-to-lake seepage in these areas. East of the major lakes, hydraulic head data were nearly equal in the surficial aquifer and Upper Floridan aquifer, which is an indication that the Upper Floridan aquifer is unconfined. Based on deuterium and oxygen stable isotope data collected in December 2011 and December 2012, there was no evidence of recharge to the Upper Floridan aquifer from the wetlands east of the major lakes; aquifer isotopic ratios did not indicate an enriched source, which is typical of lake and wetland sources. West of the major lakes, there was evidence of enriched isotopic ratios in water samples from the Upper Floridan aquifer. Differences in hydraulic head at paired wells in the surficial aquifer and Upper Floridan aquifer indicated that the surficial aquifer has the potential to recharge the Upper Floridan aquifer in the western part of the pools and west of the major lakes.

Florida↗

Conceptualization and simulation of the Edwards aquifer, San Antonio region, Texas

A new numerical ground-water-flow model (Edwards aquifer model) that incorporates important components of the latest information and plausible conceptualization of the Edwards aquifer was developed. The model includes both the San Antonio and Barton Springs segments of the Edwards aquifer in the San Antonio region, Texas, and was calibrated for steady-state (1939–46) and transient (1947–2000) conditions, excluding Travis County. Transient simulations were conducted using monthly recharge and pumpage (withdrawal) data. The model incorporates conduits simulated as continuously connected (other than being separated in eastern Uvalde and southwestern Medina Counties), one-cell-wide (1,320 feet) zones with very large hydraulic-conductivity values (as much as 300,000 feet per day). The locations of the conduits were based on a number of factors, including major potentiometric-surface troughs in the aquifer, the presence of sinking streams, geochemical information, and geologic structures (for example, faults and grabens). The simulated directions of flow in the Edwards aquifer model are most strongly influenced by the presence of simulated conduits and barrier faults. The simulated flow in the Edwards aquifer is influenced by the locations of the simulated conduits, which tend to facilitate flow. The simulated subregional flow directions generally are toward the nearest conduit and subsequently along the conduits from the recharge zone into the confined zone and toward the major springs. Structures simulated in the Edwards aquifer model influencing groundwater flow that tend to restrict flow are barrier faults. The influence of simulated barrier faults on flow directions is most evident in northern Medina County. A water budget is an accounting of inflow to, outflow from, and storage change in the aquifer. For the Edwards aquifer model steady-state simulation, recharge (from seepage losses from streams and infiltration of rainfall) accounts for 93.5 percent of the sources of water to the Edwards aquifer, and inflow through the northern and northwestern model boundaries contributes 6.5 percent. The largest discharges are springflow (73.7 percent) and ground-water withdrawals by wells (25.7 percent). The principal source of water to the Edwards aquifer for the Edwards aquifer model transient simulation was recharge, constituting about 60 percent of the sources of water (excluding change in storage) to the Edwards aquifer during 1956, a drought period, and about 97 percent of the sources (excluding change in storage) during 1975, a period of above-normal rainfall and recharge. The principal discharges from the Edwards aquifer for the transient simulation were springflow and withdrawals by wells. During 1956, representing drought conditions, the change in storage (net water released from storage) was much greater than recharge, accounting for 75.9 percent of the total flow compared to 14.5 percent for recharge. Conversely, during 1975, representing above-normal rainfall and recharge conditions, recharge constituted 79.9 percent of the total flow, compared to 7.1 percent for the change in storage (net water added to storage). A series of sensitivity tests was made to ascertain how the model results were affected by variations greater than and less than the calibrated values of input data. Simulated hydraulic heads in the Edwards aquifer model were most sensitive to recharge, withdrawals, hydraulic conductivity of the conduit segments, and specific yield and were comparatively insensitive to spring-orifice conductance, northern boundary inflow, and specific storage. Simulated springflow in the Edwards aquifer model was most sensitive to recharge, withdrawals, hydraulic conductivity of the conduit segments, specific yield, and increases in northern boundary inflow and was comparatively insensitive to spring-orifice conductance and specific storage.

Texas↗

Hydrologic data and groundwater-flow simulations in the Brown Ditch Watershed, Indiana Dunes National Lakeshore, near Beverly Shores and Town of Pines, Indiana

The U.S. Geological Survey (USGS) collected data and simulated groundwater flow to increase understanding of the hydrology and the effects of drainage alterations on the water table in the vicinity of Great Marsh, near Beverly Shores and Town of Pines, Indiana. Prior land-management practices have modified drainage and caused changes in the distribution of open water, streams and ditches, and groundwater abundance and flow paths. Collected hydrologic data indicate that the majority of water entering Great Marsh flows from the southern dune ridge beneath Town of Pines, Indiana. Groundwater flow is intercepted by Brown Ditch in the eastern portion of the study area and Derby Ditch in the western portion of the study area. A smaller amount of groundwater from the northern dune ridge beneath Beverly Shores also contributed water to Great Marsh. Continuous groundwater-level data collected indicate that the predominant north-south groundwater-flow gradients vary during the course of the year due to increased levels of precipitation or during periods of drainage obstructions. Continuous surface-water discharge and surface-water elevation were measured at three USGS streamgages, one each on Brown, Kintzele and Derby Ditches. The monthly mean discharge statistics indicate that during the period of record— June 2012 to September 2013—streamflow in Kintzele Ditch was lowest during July 2012 and highest during April 2013. In Derby Ditch, streamflow also was lowest during July 2012 and highest during April 2013. Periods of relatively high and low groundwater levels during August 1982, March 2013, and April 2014 were examined and simulated by using MODFLOW and companion software. Results from the simulation of conditions during March 2013 include that nearly 100 percent of all water entering the area simulating Town of Pines is from recharge. Of all the water simulated to enter the eastern and western portions of Great Marsh, nearly 20 and 18 percent, respectively, flows from Town of Pines to the western and eastern portions of Great Marsh. The dune ridges beneath Town of Pines and to a lesser extent beneath Beverly Shores are a major source of recharge to the surficial aquifer and Great Marsh. Results from the simulation of the conditions of April 2014 include that, despite increases in the amount of water entering Great Marsh due to a beaver-dam-modified hydrologic condition, there is still virtually zero simulated groundwater flow from Great Marsh to Town of Pines. The volume of water simulated to be entering the zone representing Beverly Shores decreased by 0.43 cubic foot per second from the results of the March 2013 simulation. This simulated difference in water budgets can be attributed to increased simulated recharge in Great Marsh and Town of Pines. Effects of the inclusion of the beaver dam included the increase of the simulated water table and simulated inundated area upstream of the beaver dam due to the effects of ponding surface water. Results from the simulation scenario that includes six proposed pool-riffle control structures in Brown Ditch under the hydrologic conditions of March 2013 indicate areas inundated by water are larger, including areas just to the north of the entrance of Brown Ditch into Great Marsh, and areas north of the confluence of Brown and Kintzele Ditches. Results from the scenario simulating the increase of the Lake Michigan water level to the historical high of May 31, 1998, showed inundated areas of Great Marsh south of Beverly Shores enlarged on both sides of Lakeshore County Road with the greatest enlargement simulated to be southeast of the intersection of Lakeshore County Road and Beverly Drive. For the scenario simulating the decrease of the Lake Michigan water level to the historical low of December 23, 2007, results show little change from the original March 2013 inundated area. The results of this study can be used by water-resource managers to understand how surrounding ditches affect water levels in Great Marsh and other inland wetlands and residential areas. The groundwater model developed can be applied to answer questions about how alterations to the drainage system in the area affects water levels in the public and residential areas surrounding Great Marsh. The modeling methods developed in this study provide a template for other studies of groundwater flow and groundwater/surface-water interactions within the shallow surficial aquifer in northern Indiana, and in similar hydrologic settings that include surficial sand aquifers in coastal areas.

Indiana↗

Use of set blanks in reporting pesticide results at the U.S. Geological Survey National Water Quality Laboratory, 2001-15

Executive Summary Background .—Pesticide results from the U.S. Geological Survey (USGS) National Water Quality Laboratory (NWQL) are used for water-quality assessments by many agencies and organizations. The USGS is committed to providing data of the highest possible quality to the consumers of its data. A cooperator’s inquiries about specific pesticide detections in water revealed potential laboratory contamination issues for some results. Consequently, the USGS conducted an extensive evaluation of potential low-level contamination related to processing or analysis of water-quality samples at NWQL for 21 pesticide compounds of interest to the cooperator. This is the most comprehensive study of NWQL pesticide quality-control (QC) results to date. Purpose and scope .—The purpose of this study was to document protocols used by the NWQL to censor pesticide results and to determine the effects of laboratory contamination—as determined from detections in laboratory set blanks—on pesticide detections in groundwater and surface-water samples. More than 30,000 pesticide results from 113 selected batches of samples (2 percent or less of total batches) analyzed by the NWQL during the 15 years from 2001 to 2015 were reviewed. All laboratory results from the selected batches, including results from environmental (surface water and groundwater) and QC (set-blank, blind-blank, and blind-spike) samples, were evaluated. The study includes results for 21 pesticide compounds analyzed in groundwater and surface-water samples collected across the United States. Eleven pesticide compounds were analyzed by a gas chromatography/mass spectrometry method and 10 compounds by a liquid chromatography/mass spectrometry method. Objectives and methods .—The objectives of this study were to (1) determine the characteristics of laboratory contamination over time, (2) compare distributions of pesticide results in set blanks with distributions in environmental samples, (3) evaluate the potential for false-positive and false-negative reporting of results, and (4) evaluate the effects of reevaluating historical pesticide results using 2017 compound identification protocols on detections of pesticides in groundwater and surface-water samples. The 113 instrument batches selected for this study contained detections of one or more of the 21 pesticide compounds in set blanks or were among those batches with the highest pesticide detection frequencies in set blanks. As a result, the dataset for this study was targeted toward pesticides and batches with laboratory contamination. The objectives were addressed by statistically comparing environmental and set-blank results; computing moving averages of set-blank detection frequencies to identify periods of episodic contamination; and using summary statistics, tabular summaries, and graphical approaches, such as time-series plots and cumulative distribution functions. Results .— Objective 1 : Laboratory contamination, as determined by pesticide detections in set blanks, was found in 13 percent of set-blank results from the 113 targeted batches included in this study (as compared to 6 percent of set-blank results from all 7,620 batches analyzed during the study period). It is estimated that 92 percent of the laboratory contamination during the study period was episodic, meaning that it occurred during discrete periods of time. All 21 of the targeted pesticide compounds had periods of episodic contamination, with most episodes ranging in duration from about 1 to 8 months. The remaining 8 percent of laboratory contamination was random or from a known source (deterministic). Objective 2 : For some compounds, graphs of cumulative distribution functions of the entire distributions of set-blank and environmental samples overlap, suggesting that there is no difference in the distributions of the two types of samples. However, time-series graphs show that detections in set blanks often occur at different times (sometimes separated by years) than detections in environmental samples, indicating clear differences in those distributions, and indicating the importance of evaluating the timing of detections in all sample types. For most compounds detected in set-blank and environmental samples, detection frequencies were significantly greater in set blanks than in groundwater or surface-water samples ( p <0.05). There are several explanations for this finding, including that the 113 batches of samples chosen for this study targeted batches with detections in set blanks or that detections in set-blank samples were historically determined with less stringent identification criteria than for environmental samples (groundwater and surface-water samples). Objective 3 : The false-positive and false-negative rates from blind samples submitted during the study period by the USGS Quality Systems Branch generally were less than 1 and 5 percent, respectively, for the 21 pesticides. The only compound with a false-positive rate greater than 1 percent was flumetsulam (2.6 percent), indicating that there is a higher likelihood of flumetsulam being reported as a detection when it is not present in an environmental sample compared with the reporting of other compounds. Objective 4 : Altogether, for data in targeted batches, NWQL would have reported 0.1 percent of results from groundwater samples and 1.4 percent of results from surface-water samples differently if 2017 identification protocols were applied to historical pesticide results. In most of these cases, detections observed in historical results would change to nondetections. The small percentages of changes that would occur if historical data were reevaluated indicate that historical protocols used by the NWQL to identify detections in environmental samples were robust and produced results that are predominantly consistent with current [2017] practices. Conclusions .—The NWQL produces high-quality pesticide results at environmentally relevant concentrations. NWQL identification protocols and censoring practices are largely effective at minimizing the reporting of false-positive and false-negative results. Laboratory contamination, when it occurred, tended to occur in episodes; thus, evaluating the timing and magnitude of detections in set blanks relative to detections in environmental samples was determined to be an important consideration for analysis of environmental results. Because NWQL censoring practices do not address all types and occurrences of laboratory contamination, options for additional censoring practices are provided for data users with more specific or stringent data-quality objectives. The methods used to analyze the 21 compounds for this report can similarly be applied to all 173 pesticide compounds that were analyzed by the NWQL during the same time period. This study also has helped to identify potential improvements in reporting USGS data, such as conducting more frequent review of set-blank datasets.

Scientific Investigations Report↗

Evaluation of Street Sweeping as a Stormwater-Quality-Management Tool in Three Residential Basins in Madison, Wisconsin

Recent technological improvements have increased the ability of street sweepers to remove sediment and other debris from street surfaces; the effect of these technological advancements on stormwater quality is largely unknown. The U.S. Geological Survey, in cooperation with the City of Madison and the Wisconsin Department of Natural Resources, evaluated three street-sweeper technologies from 2002 through 2006. Regenerative-air, vacuum-assist, and mechanical-broom street sweepers were operated on a frequency of once per week (high frequency) in separate residential basins in Madison, Wis., to measure each sweeper's ability to not only reduce street-dirt yield but also improve the quality of stormwater runoff. A second mechanical-broom sweeper operating on a frequency of once per month (low frequency) was also evaluated to measure reductions in street-dirt yield only. A paired-basin study design was used to compare street-dirt and stormwater-quality samples during a calibration (no sweeping) and a treatment period (weekly sweeping). The basis of this paired-basin approach is that the relation between paired street-dirt and stormwater-quality loads for the control and tests basins is constant until a major change is made at one of the basins. At that time, a new relation will develop. Changes in either street-dirt and/or stormwater quality as a result of street sweeping could then be quantified by use of statistical tests. Street-dirt samples collected weekly during the calibration period and twice per week during the treatment period, once before and once after sweeping, were dried and separated into seven particle-size fractions ranging from less than 63 micrometers to greater than 2 millimeters. Street-dirt yield evaluation was based on a computed mass per unit length of pounds per curb-mile. An analysis of covariance was used to measure the significance of the effect of street sweeping at the end of the treatment period and to quantify any reduction in street-dirt yield. Both the regenerative-air and vacuum-assist sweepers produced reductions in street-dirt yield at the 5-percent significance level. Street-dirt yield was reduced by an average of 76, 63, and 20 percent in the regenerative-air, vacuum-assist, and high-frequency broom basins, respectively. The low-frequency broom basin showed no significant reductions in street-dirt yield. Sand-size particles (greater than 63 micrometers) recorded the greatest overall reduction. Street-sweeper pickup efficiency was determined by computing the difference between weekly street-dirt yields before and after sweeping cleaning. The regenerative-air and vacuum-assist sweepers had similar pickup efficiencies of 25 and 30 percent, respectively. The mechanical broom sweeper operating at high frequency was considerably less efficient, removing an average of 5 percent of street-dirt yield. The effects of street sweeping on stormwater quality were evaluated by use of statistical tests to compare event mean concentrations and loads computed for individual storms at the control and test basins. Loads were computed by multiplying the event mean concentrations by storm-runoff volumes. Only ammonia-nitrogen for the test basin with the vacuum-assist sweeper showed significant load increases over the control basin, at the 10-percent significance level, of 63 percent. Difficulty in detecting significant changes in constituent stormwater-quality loads could be due, in part, to the large amount of variability in the data. Coefficients of variation for the majority of constituent loads were greater than 1, indicating substantial variability. The ability to detect changes in constituent stormwater-quality loads was likely hampered by an inadequate number of samples in the data set. However, sediment transport in the storm-sewer pipe, sediment washing onto the street from other source areas, winter sand application, and sampling challenges were additional sources of variability within each study ba

Scientific Investigations Report↗

Hydrologic drought of water year 2011 compared to four major drought periods of the 20th century in Oklahoma

Water year 2011 (October 1, 2010, through September 30, 2011) was a year of hydrologic drought (based on streamflow) in Oklahoma and the second-driest year to date (based on precipitation) since 1925. Drought conditions worsened substantially in the summer, with the highest monthly average temperature record for all States being broken by Oklahoma in July (89.1 degrees Fahrenheit), June being the second hottest and August being the hottest on record for those months for the State since 1895. Drought conditions continued into the fall, with all of the State continuing to be in severe to exceptional drought through the end of September. In addition to effects on streamflow and reservoirs, the 2011 drought increased damage from wildfires, led to declarations of states of emergency, water-use restrictions, and outdoor burning bans; caused at least $2 billion of losses in the agricultural sector and higher prices for food and other agricultural products; caused losses of tourism and wildlife; reduced hydropower generation; and lowered groundwater levels in State aquifers. The U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, conducted an investigation to compare the severity of the 2011 drought with four previous major hydrologic drought periods during the 20th century – water years 1929–41, 1952–56, 1961–72, and 1976–81. The period of water years 1925–2011 was selected as the period of record because few continuous record streamflow-gaging stations existed before 1925, and gaps in time existed where no streamflow-gaging stations were operated before 1925. In water year 2011, statewide annual precipitation was the 2d lowest, statewide annual streamflow was 16th lowest, and statewide annual runoff was 42d lowest of those 87 years of record. Annual area-averaged precipitation totals by the nine National Weather Service climate divisions from water year 2011 were compared to those during four previous major hydrologic drought periods to show how precipitation deficits in Oklahoma varied by region. The nine climate divisions in Oklahoma had precipitation in water year 2011 ranging from 43 to 76 percent of normal annual precipitation, with the Northeast Climate Division having the closest to normal precipitation and the Southwest Climate Division having the greatest percentage of annual deficit. Based on precipitation amounts, water year 2011 ranked as the second driest of the 1925–2011 period, being exceeded only in one year of the 1952 to 1956 drought period. Regional streamflow patterns for water year 2011 indicate that streamflow in the Arkansas-White-Red water resources region, which includes all of Oklahoma, was relatively large, being only the 26th lowest since 1930, primarily because of normal or above-normal streamflow in the northern part of the region. Twelve long-term streamflow-gaging stations with periods of record ranging from 67 to 83 years were selected to show how streamflow deficits varied by region in Oklahoma. Statewide, streamflow in water year 2011 was greater than streamflows measured in years during the drought periods of 1929–41, 1952–56, 1961–72, and 1976–81. The hydrologic drought worsened going from the northeast toward the southwest in Oklahoma, ranging from 140 percent (above normal streamflow) in the northeast, to 13 percent of normal streamflow in southwestern Oklahoma. The relatively low streamflow in 2011 resulted in 83.3 percent of the statewide conservation storage being available at the end of the water year in major reservoirs, similar to conservation storage in the preceding severe drought year of 2006. The ranking of streamflow as the 16th smallest for the 1925–2011 period, despite precipitation being ranked the 2d smallest, may have been caused, in part, by the relatively large streamflow in northeastern Oklahoma during water year 2011.

Oklahoma↗

Changes in phosphorus concentrations and loads in the Assabet River, Massachusetts, October 2008 through April 2014

Treated effluent discharged from municipal wastewater-treatment plants to the Assabet River in central Massachusetts includes phosphorus, which leads to increased growth of nuisance aquatic plants that decrease the river’s water quality and aesthetics in impounded reaches during the growing season. To improve the river’s water quality and aesthetics, the U.S. Environmental Protection Agency approved a total maximum daily load for phosphorus in 2004 that directed the wastewater-treatment plants to reduce the amount of total phosphorus discharged to the river by 2012. The permitted total phosphorus monthly average of 0.75 milligrams per liter during the aquatic plant growing season (April 1 through October 31) was reduced by the total maximum daily load to a target of 0.1 milligrams per liter by 2012, and the nongrowing-season limit was unchanged at 1.0 milligrams per liter. From October 2008 through April 2014, the U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, measured streamflow and collected weekly flow-proportional, composite samples of water from the Assabet River for analysis of concentrations of total phosphorus and orthophosphate. Streamflow and concentration data were used to estimate total phosphorus and orthophosphate loads in the river. The purpose of this monitoring effort was to evaluate phosphorus concentrations and loads in the river before, during, and after the wastewater-treatment-plant upgrades and to assess the effects of seasonal differences in permitted discharges. The locations of water-quality-monitoring stations, with respect to the Hudson and Ben Smith impoundments, enabled examination of effects of phosphorus entering and leaving the impoundments. Annual median concentrations of total phosphorus in wastewater-treatment plants were reduced by more than 80 percent with the plant upgrades. Measured instream annual median concentrations of total phosphorus in the Assabet River decreased by about 38 to 50 percent at three of the four monitoring stations following the wastewater-treatment-plant upgrades. At the station farthest upstream, the median total phosphorus concentration remained unchanged throughout the study; this may be attributed to the site location and potential resuspension of particulate organic matter during periods of increased streamflow. Annual median loads from the wastewater-treatment plants were reduced by up to 91 percent following the upgrades, instream annual median total phosphorus loads at the three downstream stations decreased by 71 to 76 percent, and instream orthophosphate loads at the three downstream stations decreased by 79 to 87 percent. Seasonal fluctuations (growing versus nongrowing) of total phosphorus and orthophosphate were observed instream before the upgrades. However, after the upgrades, fluctuations in phosphorus released from the treatment plants were slight and seasonal changes were typically not observed instream. Annual loads entering and leaving the two impoundments were inconclusive in determining whether the impoundments were sources or sinks of total phosphorus during the study. Total phosphorus loads entering the Hudson impoundment were consistently greater than those leaving; however, there was uncertainty about the loads at the monitoring station upstream from this impoundment. At the Ben Smith impoundment, total phosphorus and orthophosphate loads downstream were slightly greater than those upstream from the impoundment, but the differences may reflect additions from tributaries and overland runoff. Estimated instream total phosphorus concentrations and loads indicated that the decreases in total phosphorus in wastewater-treatment-plant discharges were accompanied by reductions measured in the Assabet River. A statistical analysis which incorporates the effect of varying flow conditions demonstrated significant reductions in total phosphorus concentrations after the wastewater-treatment-plant upgrades at three of the four instream monitoring stations. No significant change was observed at the most upstream location, the Assabet River at Port Street at Hudson, Massachusetts (station number 01096835), which may have been affected by flow-related resuspension of particulate phosphorus.

Massachusetts↗

Peak streamflow trends in Illinois and their relation to changes in climate, water years 1921–2020

This report characterizes changes in peak streamflow in Illinois and the relation of these changes to climatic variability, and provides a foundation for future studies that can address nonstationarity in peak-flow frequency analysis in Illinois. Records of annual peak and daily streamflow at streamgages and gridded monthly climatic data (observed and modeled) were examined across four trend periods (100 years, water years 1921–2020; 75 years, 1946–2020; 50 years, 1971–2020; 30 years 1991–2020) for trends, change points, and other statistical properties indicative of changing conditions. Median peak streamflows generally exhibit upward trends across the State for the 100- and 75-year trend periods and in northern and southern Illinois for the 50- and 30-year trend periods. The medians of the trend magnitudes (normalized by median peak streamflow) range from a 23-percent increase during the 30-year trend period to a 41-percent increase during the 100-year trend period. Streamgages with trends in peak streamflow often also have change points, or abrupt changes, in streamflow magnitude. More than two-thirds of streamgages at the 100- and 75-year trend periods exhibit a trend and change point in median peak streamflow in the same direction. Temporally, clusters of change points are observed in the late 1960s through early 1980s for the 100- and 75-year trend periods and around 2006 for the 50- and 30-year trend periods. Trends in the 90-percent quantile of peak streamflow, which correspond to the 10-percent exceedance probability often used for the design of drainage structures, increased about the same amount as the 50-percent quantile peak streamflows, except at the 100-year trend period, where the 50-percent quantile peak flow increased more for almost all streamgages. The frequency of high flows has also increased in Illinois, with increases in peaks-over-threshold observed across much of the State for the 100- and 75-year trend periods and in northern and southern Illinois for the 50- and 30-year trend periods. Upward trends in observed temperature and observed annual precipitation dominate in all trend periods, with clusters of likely upward trends observed in northern and southern Illinois at the 50- and 30-year trend periods. As expected in response to increasing temperature, the modeled proportion of precipitation falling as snow has largely decreased in the study basins across the State, and modeled potential evapotranspiration has increased. Upward trends in modeled annual runoff, which in this report incorporates only the effects of climatic variation, are observed in the same geographic areas where there are increases in observed annual precipitation. The widespread upward trends in the magnitude of median peak streamflows and the frequency with which high flows occur across the State at the 100- and 75-year trend periods and in northern and southern Illinois at the 50- and 30-year trend periods appear to be driven largely by increases in precipitation based on spatial patterns of these changes and statistical relations between streamflow and climate metrics. Other effects not considered in this report, like urbanization, may be important drivers for certain streamgages in the State. The prevalence of nonstationarity in peak streamflow in Illinois has important implications for peak-flow frequency analysis. Average annual precipitation and the occurrence of extreme precipitation events are expected to increase across the State. If precipitation continues to increase as expected, peak-flow frequency estimates based on older records may no longer represent the hydrologic regime of today, and methods for nonstationary peak-flow frequency analysis may be needed.

Illinois↗

Estimation of selected streamflow statistics for a network of low-flow partial-record stations in areas affected by Base Realignment and Closure (BRAC) in Maryland

The U.S. Geological Survey, in cooperation with the Maryland Department of the Environment, operated a network of 20 low-flow partial-record stations during 2008 in a region that extends from southwest of Baltimore to the northeastern corner of Maryland to obtain estimates of selected streamflow statistics at the station locations. The study area is expected to face a substantial influx of new residents and businesses as a result of military and civilian personnel transfers associated with the Federal Base Realignment and Closure Act of 2005. The estimated streamflow statistics, which include monthly 85-percent duration flows, the 10-year recurrence-interval minimum base flow, and the 7-day, 10-year low flow, are needed to provide a better understanding of the availability of water resources in the area to be affected by base-realignment activities. Streamflow measurements collected for this study at the low-flow partial-record stations and measurements collected previously for 8 of the 20 stations were related to concurrent daily flows at nearby index streamgages to estimate the streamflow statistics. Three methods were used to estimate the streamflow statistics and two methods were used to select the index streamgages. Of the three methods used to estimate the streamflow statistics, two of them--the Moments and MOVE1 methods--rely on correlating the streamflow measurements at the low-flow partial-record stations with concurrent streamflows at nearby, hydrologically similar index streamgages to determine the estimates. These methods, recommended for use by the U.S. Geological Survey, generally require about 10 streamflow measurements at the low-flow partial-record station. The third method transfers the streamflow statistics from the index streamgage to the partial-record station based on the average of the ratios of the measured streamflows at the partial-record station to the concurrent streamflows at the index streamgage. This method can be used with as few as one pair of streamflow measurements made on a single streamflow recession at the low-flow partial-record station, although additional pairs of measurements will increase the accuracy of the estimates. Errors associated with the two correlation methods generally were lower than the errors associated with the flow-ratio method, but the advantages of the flow-ratio method are that it can produce reasonably accurate estimates from streamflow measurements much faster and at lower cost than estimates obtained using the correlation methods. The two index-streamgage selection methods were (1) selection based on the highest correlation coefficient between the low-flow partial-record station and the index streamgages, and (2) selection based on Euclidean distance, where the Euclidean distance was computed as a function of geographic proximity and the basin characteristics: drainage area, percentage of forested area, percentage of impervious area, and the base-flow recession time constant, t. Method 1 generally selected index streamgages that were significantly closer to the low-flow partial-record stations than method 2. The errors associated with the estimated streamflow statistics generally were lower for method 1 than for method 2, but the differences were not statistically significant. The flow-ratio method for estimating streamflow statistics at low-flow partial-record stations was shown to be independent from the two correlation-based estimation methods. As a result, final estimates were determined for eight low-flow partial-record stations by weighting estimates from the flow-ratio method with estimates from one of the two correlation methods according to the respective variances of the estimates. Average standard errors of estimate for the final estimates ranged from 90.0 to 7.0 percent, with an average value of 26.5 percent. Average standard errors of estimate for the weighted estimates were, on average, 4.3 percent less than the best average standard errors of estima

Scientific Investigations Report↗

Anthropogenic organic compounds in source water of nine community water systems that withdraw from streams, 2002-05

Source water, herein defined as stream water collected at a water-system intake prior to water treatment, was sampled at nine community water systems, ranging in size from a system serving about 3,000 people to one that serves about 2 million people. As many as 17 source-water samples were collected at each site over about a 12-month period between 2002 and 2004 for analysis of 258 anthropogenic organic compounds. Most of these compounds are unregulated in drinking water, and the compounds analyzed include pesticides and selected pesticide degradates, gasoline hydrocarbons, personal-care and domestic-use compounds, and solvents. The laboratory analytical methods used in this study have relatively low detection levels - commonly 100 to 1,000 times lower than State and Federal standards and guidelines for protecting water quality. Detections, therefore, do not necessarily indicate a concern to human health but rather help to identify emerging issues and to track changes in occurrence and concentrations over time. About one-half (134) of the compounds were detected at least once in source-water samples. Forty-seven compounds were detected commonly (in 10 percent or more of the samples), and six compounds (chloroform, atrazine, simazine, metolachlor, deethylatrazine, and hexahydrohexamethylcyclopentabenzopyran (HHCB) were detected in more than one-half of the samples. Chloroform was the most commonly detected compound - in every sample (year round) at five sites. Findings for chloroform and the fragrances HHCB and acetyl hexamethyl tetrahydronaphthalene (AHTN) indicate an association between occurrence and the presence of large upstream wastewater discharges in the watersheds. The herbicides atrazine, simazine, and metolachlor also were among the most commonly detected compounds. Degradates of these herbicides, as well as those of a few other commonly occurring herbicides, generally were detected at concentrations similar to or greater than concentrations of the parent compound. Samples typically contained mixtures of two or more compounds. The total number of compounds and their total concentration in samples generally increased with the amount of urban and agricultural land use in a watershed. Annual mean concentrations of all compounds were less than human-health benchmarks. Single-sample concentrations of anthropogenic organic compounds in source water generally were less than 0.1 microgram per liter and less than established human-health benchmarks. Human-health benchmarks used for comparison were U.S. Environmental Protection Agency (USEPA) Maximum Contaminant Levels (MCLs) for regulated compounds and U.S. Geological Survey Health-Based Screening Levels for unregulated compounds. About one-half of all detected compounds do not have human-health benchmarks or adequate toxicity information for evaluating results in a human-health context. During a second sampling phase (2004-05), source water and finished water (water that has passed through all the treatment processes but prior to distribution) were sampled at eight of the nine community water systems. Water-treatment processes differ among the systems. Specifically, treatment at five of the systems is conventional, typically including steps of coagulation, flocculation, sedimentation, filtration, and disinfection. One water system uses slow sand filtration and disinfection, a second system uses ozone as a preliminary treatment step to conventional treatment, and a third system is a direct filtration treatment plant that uses many of the steps employed in conventional treatment. Most of these treatment steps are not designed specifically to remove the compounds monitored in this study. About two-thirds of the compounds detected commonly in source water were detected at similar frequencies in finished water. Although the water-treatment steps differ somewhat among the eight water systems, the amount of change in concentration of the compounds from source- to finish

Scientific Investigations Report↗

Budgets and chemical characterization of groundwater for the Diamond Valley flow system, central Nevada, 2011–12

The Diamond Valley flow system consists of six hydraulically connected hydrographic areas in central Nevada. The general down-gradient order of the areas are southern and northern Monitor Valleys, Antelope Valley, Kobeh Valley, Stevens Basin, and Diamond Valley. Groundwater flow in the Diamond Valley flow system terminates at a large playa in the northern part of Diamond Valley. Concerns relating to continued water-resources development of the flow system resulted in a phased hydrologic investigation that began in 2005 by the U.S. Geological Survey in cooperation with Eureka County. This report presents the culmination of the phased investigation to increase understanding of the groundwater resources of the basin-fill aquifers in the Diamond Valley flow system through evaluations of groundwater chemistry and budgets. Groundwater chemistry was characterized using major ions and stable isotopes from groundwater and precipitation samples. Groundwater budgets accounted for all inflows, outflows, and changes in storage, and were developed for pre-development (pre-1950) and recent (average annual 2011&ndash;12) conditions. Major budget components include groundwater discharge by evapotranspiration and groundwater withdrawals; groundwater recharge by precipitation, and interbasin flow; and storage change. Groundwater in the basin-fill aquifer of the Diamond Valley flow system was mostly a calcium or sodium bicarbonate water type and generally within acceptable drinking-water standards. The general water type was similar among the individual hydrographic areas. Stable isotopes of oxygen-18 and deuterium from precipitation varied seasonally, such that enrichment from evaporation was greater during warmer months than cooler months. The isotopic signature of shallow groundwater was similar to cool season precipitation, indicating recharge was relatively recent (similar to recent climatic conditions) and was derived from cool season precipitation. Site-scale groundwater evapotranspiration was estimated from eddy-covariance and micrometeorological measurements collected at four sites and ranged from 0.15 feet per year in sparse, undisturbed shrubland to 1.13 feet per year in a grassland meadow. Vegetation indices calculated from satellite imagery and field mapping were used to define three evapotranspiration units (shrubland, grassland, and playa) and to extrapolate site-scale groundwater evapotranspiration rates to basin-scale estimates. Annual pre-development groundwater evapotranspiration for individual hydrographic areas ranged from 2,900 acre-feet per year (acre-ft/yr) in northern Monitor Valley to 35,000 acre-ft/yr in Diamond Valley. Total groundwater evapotranspiration from the Diamond Valley flow system under pre-development conditions was about 70,000 acre-ft/yr. Areas of irrigated land in the Diamond Valley flow system increased from less than 5,000 acres in the early 1960s to more than 25,000 acres in 2012 and are mostly for growing alfalfa in southern Diamond Valley. Annual (2011&ndash;12) net groundwater withdrawals for irrigation, assumed to be the volume of groundwater consumed by crops and pastureland, ranged from about 420 acre-ft/yr in Antelope Valley to 67,000 acre-ft/yr in Diamond Valley. Total net groundwater withdrawals for irrigation in the Diamond Valley flow system were about 69,000 acre-ft/yr (2011&ndash;12). Groundwater recharge, the largest inflow component to the Diamond Valley flow system, was determined as the sum of groundwater evapotranspiration and net subsurface outflow (subsurface outflow minus subsurface inflow). Annual groundwater recharge estimates ranged from 200 acre-ft/yr in Stevens Basin to 35,000 acre-ft/yr in Diamond Valley. Subsurface flow between hydrographic basins was evaluated using estimated transmissivity, groundwater-flow sections derived from remotely sensed imagery, and hydraulic gradients determined from 2012 water-level data. Subsurface outflow ranged from 0 acre-ft/yr for Diamond Valley to 3,400 acre-ft/yr for northern Monitor Valley into western Kobeh Valley. Subsurface inflow ranged from 0 acre-ft/yr for southern Monitor Valley to 4,200 acre-ft/yr for Kobeh Valley from northern Monitor and Antelope Valleys. The pre-development, steady state, groundwater budget for the Diamond Valley flow system was estimated at about 70,000 acre-ft/yr of inflow and outflow. During years 2011&ndash;12, inflow components of groundwater recharge from precipitation and subsurface inflow from adjacent basins totaled 70,000 acre-ft/yr for the DVFS, whereas outflow components included 64,000 acre-ft/yr of groundwater evapotranspiration and 69,000 acre-ft/yr of net groundwater withdrawals, or net pumpage. Spring discharge in northern Diamond Valley declined about 6,000 acre-ft/yr between pre-development time and years 2011&ndash;12. Assuming net groundwater withdrawals minus spring flow decline is equivalent to the storage change, the 2011&ndash;12 summation of inflow and storage change was balanced with outflow at about 133,000 acre-ft/yr.

Nevada↗

Surface-Water Hydrology and Quality at the Pike Hill Superfund Site, Corinth, Vermont, October 2004 to December 2005

The hydrology and quality of surface water in and around the Pike Hill Brook watershed, in Corinth, Vermont, was studied from October 2004 to December 2005 by the U.S. Geological Survey in cooperation with the U.S. Environmental Protection Agency (USEPA). Pike Hill was mined intermittently for copper from 1847 to 1919 and the site is known to be contributing trace elements and acidity to Pike Hill Brook and an unnamed tributary to Cookville Brook. The site has been listed as a Superfund site since 2004. Streamflow, specific conductance, pH, and water temperature were measured continuously and monthly at three sites on Pike Hill Brook to determine the variation in these parameters over an annual cycle. Synoptic water-quality sampling was done at 10 stream sites in October 2004, April 2005, and June 2005 and at 13 stream sites in August 2005 to characterize the quality of surface water in the watershed on a seasonal and spatial basis, as well as to assess the effects of wetlands on water quality. Samples for analysis of benthic macroinvertebrate populations were collected at 11 stream sites in August 2005. Water samples were analyzed for 5 major ions and 32 trace elements. Concentrations of trace elements at sites in the Pike Hill Brook watershed exceeded USEPA National Recommended Water Quality Criteria acute and chronic toxicity standards for aluminum, iron, cadmium, copper, and zinc. Concentrations of copper exceeded the chronic criteria in an unnamed tributary to Cookville Brook in one sample. Concentrations of sulfate, calcium, aluminum, iron, cadmium, copper, and zinc decreased with distance from a site directly downstream from the mine (site 1), as a result of dilution and through sorption and precipitation of the trace elements. Maximum concentrations of aluminum, iron, cadmium, copper, and zinc were observed during spring snowmelt. Concentrations of sulfate, calcium, cadmium, copper, and zinc, and instantaneous loads of calcium and aluminum were statistically different (p<0.05) among the three continuously monitored sites (sites 1, 4, and 5). Instantaneous loads of aluminum, iron, and copper decreased by one to three orders of magnitude from site 1 to a site 1.1 mi downstream (site 4). Instantaneous loads of sulfate were similar between sites 1, 4, and at a site 3 mi downstream (site 5). Instantaneous loads of cadmium and zinc were similar between sites 1 and 4, and loads of iron and copper were similar between sites 4 and 5. Loads of chemical constituents were compared at site 1 (closest to the mine waste piles) and site 5 (near the mouth of Pike Hill Brook and below a majority of the wetlands). Annually, the loads of dissolved cadmium and zinc at site 1 were about five times greater than loads at site 5, and the load of dissolved copper at site 1 was about 17 times greater than at site 5. The ratio of loads for dissolved cadmium, copper, and zinc to total cadmium, copper, and zinc at site 1 was about 1. Samples collected in Pike Hill Brook upstream and downstream from the wetlands during low flows in August 2005 showed that oxidation of ferrous iron and precipitation of iron-hydroxides were probably not affecting trace metals in the wetlands through sorption; however, a significant portion of the iron entering the wetlands was in particulate form and may have transported sorbed copper and other trace metals. Thus, aerobic activity in the wetlands was probably not affecting metal cycling in the watershed. Concentrations and loads of sulfate may be unlikely to define unequivocally the role of the wetlands with regard to anaerobic bacterial sulfate reduction; however, bacterial sulfate removal may have affected loads of sulfate. Loads of copper increased downstream from the wetlands and may reflect the reductive dissolution of ferric hydroxide particulates in anaerobic parts of the wetlands.Concentrations of dissolved iron increased downstream from the wetlands. The most apparent effects on the macroinvertebr

Scientific Investigations Report↗

Characterization of ambient groundwater quality within a statewide, fixed-station monitoring network in Pennsylvania, 2015–19

Pennsylvania leads the Nation in the number of individuals that use groundwater for private domestic water supply; more than 3 million rural and suburban Pennsylvania residents rely on private domestic supplies for drinking water. These supplies are not regulated nor routinely monitored; thus relevant groundwater-quality information is not widely available. The U.S. Geological Survey (USGS), in cooperation with the Pennsylvania Department of Environmental Protection (PaDEP) Safe Drinking Water Bureau, established a statewide, fixed-station ambient groundwater quality network in 2015. The goals for the Pennsylvania Groundwater Monitoring Network (GWMN) include characterizing ambient groundwater quality conditions in rural areas of the State and documenting potential changes in conditions over time. Seventeen wells were selected for monitoring at 6-month intervals beginning in 2015. Since then, several wells have been added to the GWMN, bringing the total number of wells sampled in the fall of 2019 to 28. Routinely monitored constituents included physical characteristics and chemical concentrations in filtered and unfiltered samples (major and trace elements, nutrients, and organic compounds). Samples for volatile organic compounds (VOCs), radionuclides, and dissolved hydrocarbon gases were collected during the first sampling event at each well. To offer insights on the quality of groundwater used for domestic supply in Pennsylvania, summary statistics for the 221 GWMN samples collected during 2015–19 are compared to U.S. Environmental Protection Agency (EPA) drinking-water standards, which are applicable to public water supplies. Results show that samples across the GWMN generally meet drinking-water standards for inorganic and organic constituents; however, a percentage of samples had concentrations that exceeded maximum contaminant level (MCL) thresholds for nitrate (3 percent) and secondary maximum contaminant level (SMCL) thresholds for iron (32 percent), manganese (36 percent), and aluminum (5 percent). Radon-222 activities, which were sampled only during the initial visit to a well, exceeded the lower proposed drinking water standard of 300 picocuries per liter (pCi/L) in 64 percent of wells in the GWMN; additionally, 7 percent of wells exceeded the higher proposed standard of 4,000 pCi/L. There were no exceedances for VOCs, but one well had a tribromomethane detection. Three wells had detectable concentrations of methane, with one sample exceeding the Pennsylvania action level of 7 milligrams per liter (mg/L). The pH and dissolved oxygen concentrations varied widely across the GWMN and were correlated with dissolved metal concentrations and other chemical characteristics of groundwater samples. Considering all samples collected for the study, the pH ranged from 4.2 to 8.3; 42 percent of pH values were either above or below the SMCL range of 6.5–8.5. The highest pH values resulted from contamination of loose grout used in the construction of one well and decreased to levels consistent with other wells in the vicinity after repeated sampling rounds. Dissolved oxygen (DO), which ranged from 0 to 13.9 mg/L, influences the mobility and prevalence of constituents with variable oxidation state, including iron, manganese, and nitrogen species. Samples with acidic pH (less than 6.5) and (or) low DO had the highest concentrations of manganese and iron, whereas those with neutral to alkaline pH values had the highest concentrations of calcium, magnesium, sodium, and other major ions. Analysis of major ions indicates that calcium/bicarbonate water types are the most common, with a few characterized as calcium/chloride or sodium/chloride, and most others as mixed water types including calcium-magnesium/bicarbonate, sodium-magnesium/bicarbonate, and sodium/bicarbonate-chloride. Nonparametric statistical methods were used to evaluate the data for spatial and temporal trends. A principal components analysis (PCA) model developed with ranked data values for the entire network resulted in three components, (1) dissolved solids, (2) redox, and (3) sodium-chloride, which explained 74.5 percent of variance in the dataset. On the basis of individual contributions to the PCA, certain wells were identified through hierarchical cluster analysis that shared relevant water-quality characteristics. The spatial distribution of sampling locations and the temporal trends of constituent concentrations indicate that hydrogeologic setting and topographic position as defined in the PCA model are important factors affecting the spatial and temporal patterns of groundwater quality in the GWMN.

Pennsylvania↗

Analysis of the Shallow Groundwater Flow System at Fire Island National Seashore, Suffolk County, New York

Fire Island National Seashore (FIIS) occupies 42 kilometers of the barrier island for which it is named that lies off the southern shore of Suffolk County, N.Y. Freshwater in the highly permeable, sandy aquifer underlying Fire Island is bounded laterally by marine surface waters and at depth by saline groundwater. Interspersed throughout FIIS are 17 pre-existing residential communities that in summer months greatly increase in population through the arrival of summer residents and vacationers; in addition, the National Park Service (NPS) has established several facilities on the island to accommodate visitors to FIIS. The 2.2 million people estimated by the NPS to visit Fire Island annually impact groundwater quality through the release of waste-derived contaminants, such as nutrients, pathogens, and organic compounds, into the environment. Waste-contaminated groundwater can move through the aquifer and threaten the ecological health of the adjacent back-barrier estuaries to which much of the groundwater ultimately discharges. In 2004, the U.S. Geological Survey (USGS), in cooperation with the NPS, began a 3-year investigation to (1) collect groundwater levels and water-quality (nutrient) samples, (2) develop a three-dimensional model of the shallow (water-table) aquifer system and adjacent marine surface waters, and (3) calculate nitrogen loads in simulated groundwater discharges from the aquifer to back-barrier estuaries and the ocean. The hydrogeology of the shallow aquifer system was characterized from the results of exploratory drilling, geophysical surveying, water-level monitoring, and water-quality sampling. The investigation focused on four areas-the communities of Kismet and Robbins Rest, the NPS Visitor Center at Watch Hill, and the undeveloped Otis Pike Fire Island High Dune Wilderness. Thirty-five observation wells were installed within FIIS to characterize subsurface hydrogeology and establish a water-table monitoring network in the four study areas. A variable-density model of the shallow aquifer system and adjacent marine surface waters was developed to simulate groundwater flow patterns and rates. Nitrogen loads from the shallow aquifer system were calculated from representative total nitrogen (TN) concentrations and simulated groundwater discharges to back-barrier estuaries and the ocean. The model simulates groundwater directions, velocities, and discharge rates under 2005 mean annual conditions. Groundwater budgets were developed for recharge areas of similar land use that contribute freshwater to back-barrier estuaries, the ocean, and subsea-discharge zones. Total freshwater discharge from the shallow aquifer system is about 43,500 cubic meters per day (m3/d) (79.8 percent) to back-barrier estuaries and about 10,200 m3/d (18.7 percent) to the ocean; about 836 m3/d (1.5 percent) may exit the system as subsea underflow. The total contribution of fresh groundwater to shoreline discharge zones amounts to about 53,700 m3/d (98.5 percent). The median age of freshwater discharged to back-barrier estuaries and the ocean was 3.4 years, and the 95th-percentile age was 20 years. The TN concentrations and loads under 2005 mean annual conditions for areas that contribute fresh groundwater to back-barrier estuaries and the ocean were calculated for the principal land uses on Fire Island. The overall TN load from the shallow aquifer system to shoreline discharge zones is about 16,200 kilograms per year (kg/yr) (82.2 percent) to back-barrier estuaries and about 3,500 kg/yr (17.8 percent) to the ocean. The overall TN load to marine surface waters amounts to about 19,700 kg/yr-roughly 6 percent of the annual TN load from shallow groundwater entering the South Shore Estuary Reserve (SSER) from the Suffolk County mainland, which is about 345,000 kg/yr. In contrast to the TN load from shallow groundwater for the SSER watershed, which annually yields about 353 kilograms per square kilometer (kg/km2), the overall TN loa

Scientific Investigations Report↗

Environmental factors affecting mercury in Camp Far West Reservoir, California, 2001-03

This report documents water quality in Camp Far West Reservoir from October 2001 through August 2003. The reservoir, located at approximately 300 feet above sea level in the foothills of the northwestern Sierra Nevada, California, is a monomictic lake characterized by extreme drawdown in the late summer and fall. Thermal stratification in summer and fall is coupled with anoxic conditions in the hypolimnion. Water-quality sampling was done at approximately 3-month intervals on eight occasions at several stations in the reservoir, including a group of three stations along a flow path in the reservoir: an upstream station in the Bear River arm (principal tributary), a mid-reservoir station in the thalweg (prereservoir river channel), and a station in the deepest part of the reservoir, in the thalweg near Camp Far West Dam. Stations in other tributary arms of the reservoir included those in the Rock Creek arm of the reservoir, a relatively low-flow tributary, and the Dairy Farm arm, a small tributary that receives acidic, metal-rich drainage seasonally from the inactive Dairy Farm Mine, which produced copper, zinc, and gold from underground workings and a surface pit. Several water-quality constituents varied significantly by season at all sampling stations, including major cations and anions, total mercury (filtered and unfiltered samples), nitrogen (ammonia plus organic), and total phosphorus. A strong seasonal signal also was observed for the sulfurisotope composition of aqueous sulfate from filtered water. Although there were some spatial differences in water quality, the seasonal variations were more profound. Concentrations of total mercury (filtered and unfiltered water) were highest during fall and winter; these concentrations decreased at most stations during spring and summer. Anoxic conditions developed in deep parts of the reservoir during summer and fall in association with thermal stratification. The highest concentrations of methylmercury in unfiltered water were observed in samples collected during summer from deepwater stations in the anoxic hypolimnion. In the shallow (less than 14 meters depth) oxic epilimnion, concentrations of methylmercury in unfiltered water were highest during the spring and lowest during the fall. The ratio of methylmercury to total mercury (MeHg/HgT) increased systematically from winter to spring to summer, largely in response to the progressive seasonal decrease in total mercury concentrations, but also to some extent because of increases in MeHg concentrations during summer. Water-quality data for Camp Far West Reservoir are used in conjunction with data from linked studies of sediment and biota to develop and refine a conceptual model for mercury methylation and bioaccumulation in the reservoir and the lower Bear River watershed. It is hypothesized that MeHg is produced by sulfate-reducing bacteria in the anoxic parts of the water column and in shallow bed sediment. Conditions were optimal for this process during late summer and fall. Previous work has indicated that Camp Far West Reservoir is a phosphate-limited system - molar ratios of inorganic nitrogen to inorganic phosphorus in filtered water were consistently greater than 16 (the Redfield ratio), sometimes by orders of magnitude. Therefore, concentrations of orthophosphate were expectedly very low or below detection at all stations during all seasons. It is further hypothesized that iron-reducing bacteria facilitate release of phosphorus from iron-rich sediments during summer and early fall, stimulating phytoplankton growth in the fall and winter, and that the MeHg produced in the hypolimnion and metalimnion is released to the entire water column in the late fall during reservoir destratification (vertical mixing). Mercury bioaccumulation factors (BAF) were computed using data from linked studies of biota spanning a range of trophic position: zooplankton, midge larvae, mayfly nymphs, crayfish, threadfin shad, bluegill,

California↗

Peak streamflow trends in Missouri and their relation to changes in climate, water years 1921–2020

This report characterizes changes in peak streamflow in Missouri and the relation of these changes to climatic variability, and provides a foundation for future studies that can address nonstationarity in peak-streamflow frequency analysis in Missouri. Records of annual peak and daily streamflow at streamgages and gridded monthly climatic data (observed and modeled) were examined across four trend periods (100 years, water years 1921–2020; 75 years, 1946–2020; 50 years, 1971–2020; and 30 years, 1991–2020) for trends, change points (abrupt changes in the streamflow time series), and other statistical properties indicative of changing conditions. Peak streamflow magnitudes generally exhibit upward trends across the State for the 100-, 75-, and 50-year trend periods and only in southern Missouri for the 30-year trend period. The medians of the trend magnitudes (normalized by median peak streamflow) range from a 10-percent increase during the 30-year trend period to a 40-percent increase during the 100-year trend period. Changes in the 90-percent quantile of peak streamflow, which correspond to the 10-percent exceedance probability often used for the design of drainage structures, are not as substantial or widespread, showing consistent increases mainly in the southern part of the State in the 50- and 30-year trend periods. Streamgages with trends in peak streamflow often also have change points, or abrupt changes, in streamflow magnitude. Change points in peak streamflows generally follow that of the peak streamflow trends, with upward change points throughout most of the State at the 100- and 75-year trend periods and in southern Missouri at the 30-year trend period. Temporally, clusters upward of change points are observed in the 1970s through 1980s for the 100-, 75-, and 50-year trend periods and around 2006 and 2007 for the 50- and 30-year trend periods. A peaks-over-threshold analysis, which evaluates changes in the frequency of peak streamflows over a certain threshold, indicates that high flows have increased in frequency at 50 to 64 percent of streamgages in the 100- and 75-year trend periods. Most streamgages in the 50- and 30-year trend periods exhibit no change. Although the frequency of high flows has increased at some streamgages and trend periods in Missouri, these increases are not as widespread as the increases in the magnitude of peak streamflow. Upward trends in observed temperature and observed annual precipitation dominate in all trend periods, with no downward trends in precipitation and only two somewhat likely downward trends in temperature for the 100-year trend period. Increases in annual precipitation mostly are limited to southern Missouri for the 30-year trend period. The proportion of precipitation falling as snow has largely decreased in the study basins across the State, which is expected in response to increasing temperature. Upward trends in modeled annual runoff, which in this study incorporates only the effects of climatic variation, are observed in the same geographic areas where there are increases in observed annual precipitation. When peak streamflow and climatic trends are considered together, widespread upward trends in peak streamflows for the 100-, 75-, and 50-year trend periods and for the 30-year trend period mainly in southern Missouri (encompassing both trends and abrupt change) appear to be driven largely by increases in precipitation based on spatial patterns and statistical relations. The prevalence of nonstationarity in peak streamflow in Missouri has important implications for peak-flow frequency analysis. Winter and spring precipitation and the occurrence of extreme precipitation events are expected to increase across the State. If precipitation continues to increase as expected, peak-flow frequency estimates based on older records may no longer represent the hydrologic regime of today, and methods for nonstationary peak-flow frequency analysis may be needed.

Missouri↗

Estimated reductions in phosphorus loads from removal of leaf litter in the Lake Champlain drainage area, Vermont

Excess nutrient loading and other factors are driving eutrophication and other negative effects on water-quality conditions in Lake Champlain and other receiving waters in Vermont. Two common best management practices were evaluated to determine how these practices can be optimized by targeting maintenance and operation to align better with seasonally driven needs, specifically to help municipalities remove a greater proportion of seasonal leaves and organic debris, reduce nutrient loading, and achieve water-quality goals. To characterize solid materials typically removed by the municipal BMPs of catch-basin (CB) cleaning and street cleaning (SC), subsamples of CB and SC materials were collected each month from nine participating municipalities in central and northwestern Vermont between September 2017 and November 2018. Monthly and seasonal composites of CB and SC samples were created from the subsamples of available materials from all municipalities. Samples were analyzed for concentrations of total organic carbon, total Kjeldahl nitrogen, and total phosphorus (P), and separated into three particle-size fractions. Distribution of particle-size fractions was similar between CB and SC as both practices generally collect the coarser fraction of solid materials (greater than 125 micrometers in diameter). In the fall, however, the range of the coarser fraction of materials increased. This is attributed to the ability of SC to collect leaves and other light organic materials that commonly pass through a CB system designed to trap heavier materials. Total organic carbon, total Kjeldahl nitrogen, and total P concentrations were highest in the catch-basin samples in the fall of 2017, and concentrations in the SC samples were highest in the fall of 2018. The collection of fewer samples in 2017 may account for some of the variability between fall 2017 and fall 2018 results. A subset of SC samples collected from piles representing specific street-cleaning routes in September and November 2018 were also analyzed. Materials collected in November were dominated by leaves, and the concentrations of the analyzed species of carbon, nitrogen, and phosphorus in some samples were more than double those in samples collected on the same street-cleaning routes in September. The Vermont Department of Environmental Conservation and the University of Vermont developed estimates of load-reduction credits for CB and SC practices based on a policy developed by the Wisconsin Department of Natural Resources that determined the potential for credits associated with leaf-removal activities. This process also considered BMPs that were initiated during the U.S. Environmental Protection Agency’s Lake Champlain Basin Total Maximum Daily Load monitoring period (2000 to 2009) and adapted the Wisconsin Department of Natural Resources policies to apply to existing SC routes in the cooperating Vermont municipalities that possessed at least 17 percent tree cover. This exercise demonstrated that applying the Wisconsin Department of Natural Resources policy to existing street-cleaning routes possessing 17 percent or more tree cover would result in reductions in total P loads up to 65 percent of mandated target reductions, and about a 25 percent reduction on average. Continuous simulations of stormwater runoff volume, and of loads of suspended sediments and total P, also were created for Englesby Brook Basin, an urbanized basin in Burlington and South Burlington that drains to Lake Champlain. Although the basin is more developed than the average of the nine cooperating municipalities, streamflow and P loading data collected by the U.S. Geological Survey were available to evaluate model performance. Simulations based on a year of average climatic conditions projected potential small reductions in total P of 0.08 to 0.10 percent as a result of CB cleaning and SC practices. Simulated weekly SC practices, however, reduced street-solid loads by as much as 7 percent. When the proportion of total P seen in fall SC materials collected in Vermont was applied to these simulated street-solid loads, estimated reductions of total P were about 29 percent. The combination of analytical results, estimated load-reduction credits, and simulated reductions indicate that targeted increases of SC activities to reduce leaf loading in the fall have the potential to reduce loading to receiving waters and could help regulated communities meet their water-quality goals.

Vermont↗

Effects of streamflows on stream-channel morphology in the eastern Niobrara National Scenic River, Nebraska, 1988–2010

The Niobrara River is an important and valuable economic and ecological resource in northern Nebraska that supports ecotourism, recreational boating, wildlife, fisheries, agriculture, and hydroelectric power. Because of its uniquely rich resources, a 122-kilometer reach of the Niobrara River was designated as a National Scenic River in 1991, which has been jointly managed by the U.S. Fish and Wildlife Service and National Park Service. To assess how the remarkable qualities of the National Scenic River may change if consumptive uses of water are increased above current levels, the U.S. Geological Survey, in cooperation with the National Park Service, initiated an investigation of how stream-channel morphology might be affected by potential decreases in summer streamflows. The study included a 65-kilometer segment in the wide, braided eastern stretch of the Niobrara National Scenic River that provides important nesting habitat for migratory bird species of concern to the Nation. The study focused on three river segments, separated at the confluences with two tributaries, Plum Creek and Long Pine Creek. With an overall temporal scope of 1988&ndash;2010 that includes a short interval preceding and a long interval following the Niobrara National Scenic River Designation Act of 1991, the study analyzed five separate time periods: 1988&ndash;93, 1994&ndash;99, 2000&ndash;3, 2004&ndash;6, and 2007&ndash;10, each of which ended with a year in which aerial photography coverage was available. Streamflow duration was analyzed for one streamgage upstream from the study area and two streamgages on tributary streams within the study area. Summer streamflows (July, August, and September) were targeted for analysis because median flows of the Niobrara River are lowest during those 3 months. In addition, peak flows during the study period were used to estimate bankfull discharge, which is one determinant of channel dimensions. Changes in channel morphology were examined using aerial photographs from 1993, 1999, 2003, 2006, and 2010 to measure channel width, area of islands, and incipient flood-plain surfaces, and to compute the braided index. Channel metrics were computed for each photography year and summarized by river segment. Additionally, at fixed-location cross sections, photography analysis identified localized geomorphic change to infer processes. Accuracy of geomorphic feature classification was estimated and the root-mean-square difference (RMSD) between aerial photographs was calculated to determine associated errors in channel metric calculations. The horizontal accuracy of boundaries delineated in the classification was estimated as 5 meters (m) for boundaries based on 1993 aerial photography and 4 m for all other aerial photography. The RMSD between aerial photography years ranged from 3.04 m to 4.16 m. The largest measurable changes in channel metrics were measured between 1993 and 1999 and between 1999 and 2003. Between 1993 and 1999, average total channel width increased by 9 m (3 percent) and 14 m (5 percent) in segments 2 and 3, respectively; average active channel width increased by 13 m (5 percent) in segment 3 and decreased by 6 m (4 percent) in segment 1; and incipient flood-plain-surface area increased by 40, 44, and 33 percent in segments 1, 2, and 3, respectively. Changes in channel metrics between 1999 and 2003 included a decrease in average total channel width of 14 m (5 percent) in segment 2; a decrease in active channel widths of 8 m (3 percent) and 6 m (2 percent) in segments 2 and 3, respectively; and an increase of 5 m (3 percent) in segment 1. Incipient flood-plain areas decreased by 22 and 33 percent in segments 1 and 2, respectively, and increased by 42 percent in segment 3. Large changes were measured between 1993 and 1999, and between 1999 and 2003, at many of the fixed-location cross sections. Large changes (that is, greater than 25 percent) in total channel width were measured in all three segments between 1993 and 1999 and again between 1999 and 2003; large increases were dominant between 1993 and 1999 and large decreases were dominant between 1999 and 2003. Segment 1 was the most susceptible to localized changes as there was only one period (between 2003 and 2006) in which the active channel width largely changed in fewer than 10 percent of the cross sections. Changes in channel metrics generally corresponded to changes in streamflow conditions, but other than changes in incipient flood-plain area, these changes were small and were not measured in all three segments simultaneously. Increases in total channel width (except in segment 1) and incipient flood-plain area between 1993 and 1999 corresponded to increases in streamflow. Channel narrowing (except in segment 1) between 1999 and 2003 corresponded to lower summer streamflows and extended durations of very low summer streamflow. Although the pattern of low summer streamflow and extended durations of very low summer streamflow continued during the 2004&ndash;6 period and at the beginning of the 2007&ndash;10 period, no further narrowing was measured. Consistent tributary summer inflows help to explain the resistance of segments 2 and 3 to further narrowing. Because segment 1 is already much narrower than segments 2 and 3, its average current velocity is likely to be swifter and, therefore, competent to offset further effects of the processes that led to its narrowness.

Nebraska↗