Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Environmental Geology”

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 973 records · Page 54Linked to original sources
USGS2013

Evaluation of total phosphorus mass balance in the lower Boise River and selected tributaries, southwestern Idaho

he U.S. Geological Survey (USGS), in cooperation with Idaho Department of Environmental Quality, developed spreadsheet mass-balance models for total phosphorus using results from three synoptic sampling periods conducted in the lower Boise River watershed during August and October 2012, and March 2013. The modeling reach spanned 46.4 river miles (RM) along the Boise River from Veteran’s Memorial Parkway in Boise, Idaho (RM 50.2), to Parma, Idaho (RM 3.8). The USGS collected water-quality samples and measured streamflow at 14 main-stem Boise River sites, two Boise River north channel sites, two sites on the Snake River upstream and downstream of its confluence with the Boise River, and 17 tributary and return-flow sites. Additional samples were collected from treated effluent at six wastewater treatment plants and two fish hatcheries. The Idaho Department of Water Resources quantified diversion flows in the modeling reach. Total phosphorus mass-balance models were useful tools for evaluating sources of phosphorus in the Boise River during each sampling period. The timing of synoptic sampling allowed the USGS to evaluate phosphorus inputs to and outputs from the Boise River during irrigation season, shortly after irrigation ended, and soon before irrigation resumed. Results from the synoptic sampling periods showed important differences in surface-water and groundwater distribution and phosphorus loading. In late August 2012, substantial streamflow gains to the Boise River occurred from Middleton (RM 31.4) downstream to Parma (RM 3.8). Mass-balance model results indicated that point and nonpoint sources (including groundwater) contributed phosphorus loads to the Boise River during irrigation season. Groundwater exchange within the Boise River in October 2012 and March 2013 was not as considerable as that measured in August 2012. However, groundwater discharge to agricultural tributaries and drains during non-irrigation season was a large source of discharge and phosphorus in the lower Boise River in October 2012 and March 2013. Model results indicate that point sources represent the largest contribution of phosphorus to the Boise River year round, but that reductions in point and nonpoint source phosphorus loads may be necessary to achieve seasonal total phosphorus concentration targets at Parma (RM 3.8) from May 1 through September 30, as set by the 2004 Snake River-Hells Canyon Total Maximum Daily Load document. The mass-balance models do not account for biological or depositional instream processes, but are useful indicators of locations where appreciable phosphorus uptake or release by aquatic plants may occur.

Idaho↗
USGS2014

Spatiotemporal variations in estrogenicity, hormones, and endocrine-disrupting compounds in influents and effluents of selected wastewater-treatment plants and receiving streams in New York, 2008-09

Endocrine-disrupting compounds (EDCs) in wastewater effluents have been linked to changes in sex ratios, intersex (in males), behavioral modifications, and developmental abnormalities in aquatic organisms. Yet efforts to identify and regulate specific EDCs in complex mixtures are problematic because little is known about the estrogen activity (estrogenicity) levels of many common and emerging contaminants. The potential effects of EDCs on the water quality and health of biota in streams of the New York City water supply is especially worrisome because more than 150 wastewater-treatment plants (WWTPs) are permitted to discharge effluents into surface waters and groundwaters of watersheds that provide potable water to more than 9 million people. In 2008, the U.S. Geological Survey (USGS), the New York State Department of Environmental Conservation (NYSDEC), New York State Department of Health (NYSDOH), and New York City Department of Environmental Protection (NYCDEP) began a pilot study to increase the understanding of estrogenicity and EDCs in effluents and receiving streams mainly in southeastern New York. The primary goals of this study were to document and assess the spatial and temporal variability of estrogenicity levels; the effectiveness of various treatment-plant types to remove estrogenicity; the concentrations of hormones, EDCs, and pharmaceuticals, personal care products (PPCPs); and the relations between estrogenicity and concentrations of hormones, EDCs, and PPCPs. The levels of estrogenicity and selected hormones, non-hormone EDCs, and PPCPs were characterized in samples collected seasonally in effluents from 7 WWTPs, once or twice in effluents from 34 WWTPs, and once in influents to 6 WWTPs. Estrogenicity was quantified, as estradiol equivalents, using both the biological e-screen assay and a chemical model. Results generally show that (1) estrogenicity levels in effluents varied spatially and seasonally, (2) a wide range of known and unknown EDCs were present in both WWTP effluents and receiving streams, (3) some effluents may be important sources of estrogenicity in weakly diluted streams, (4) measured levels of biological estrogenicity were often higher than estimated levels of chemical estrogenicity, and (5) the type of treatment had a large effect on the removal efficacy, and consequently, the estrogenicity levels observed in treated effluents.

New York↗
USGS2016

Streamflow characteristics based on data through water year 2009 for selected streamflow-gaging stations in or near Montana: Chapter E in Montana StreamStats

Chapter E of this Scientific Investigations Report documents results from a study by the U.S. Geological Survey, in cooperation with the Montana Department of Environmental Quality and the Montana Department of Natural Resources and Conservation, to provide an update of statewide streamflow characteristics based on data through water year 2009 for streamflow-gaging stations in or near Montana. Streamflow characteristics are presented for 408 streamflow-gaging stations in Montana and adjacent areas having 10 or more years of record. Data include the magnitude and probability of annual low and high streamflow, the magnitude and probability of low streamflow for three seasons (March–June, July–October, and November–February), streamflow duration statistics for monthly and annual periods, and mean streamflows for monthly and annual periods. Streamflow is considered to be regulated at streamflow-gaging stations where dams or other large-scale human modifications affect 20 percent or more of the contributing drainage basin. Separate streamflow characteristics are presented for the unregulated and regulated periods of record for streamflow-gaging stations with sufficient data.

Montana↗
USGS2016

Effects of water-supply reservoirs on streamflow in Massachusetts

State and local water-resource managers need modeling tools to help them manage and protect water-supply resources for both human consumption and ecological needs. The U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, has developed a decision-support tool to estimate the effects of reservoirs on natural streamflow. The Massachusetts Reservoir Simulation Tool is a model that simulates the daily water balance of a reservoir. The reservoir simulation tool provides estimates of daily outflows from reservoirs and compares the frequency, duration, and magnitude of the volume of outflows from reservoirs with estimates of the unaltered streamflow that would occur if no dam were present. This tool will help environmental managers understand the complex interactions and tradeoffs between water withdrawals, reservoir operational practices, and reservoir outflows needed for aquatic habitats. A sensitivity analysis of the daily water balance equation was performed to identify physical and operational features of reservoirs that could have the greatest effect on reservoir outflows. For the purpose of this report, uncontrolled releases of water (spills or spillage) over the reservoir spillway were considered to be a proxy for reservoir outflows directly below the dam. The ratio of average withdrawals to the average inflows had the largest effect on spillage patterns, with the highest withdrawals leading to the lowest spillage. The size of the surface area relative to the drainage area of the reservoir also had an effect on spillage; reservoirs with large surface areas have high evaporation rates during the summer, which can contribute to frequent and long periods without spillage, even in the absence of water withdrawals. Other reservoir characteristics, such as variability of inflows, groundwater interactions, and seasonal demand patterns, had low to moderate effects on the frequency, duration, and magnitude of spillage. The reservoir simulation tool was used to simulate 35 single- and multiple-reservoir systems in Massachusetts over a 44-year period (water years 1961 to 2004) under two water-use scenarios. The no-pumping scenario assumes no water withdrawal pumping, and the pumping scenario incorporates average annual pumping rates from 2000 to 2004. By comparing the results of the two scenarios, the total streamflow alteration can be parsed into the portion of streamflow alteration caused by the presence of a reservoir and the additional streamflow alteration caused by the level of water use of the system. For each reservoir system, the following metrics were computed to characterize the frequency, duration, and magnitude of reservoir outflow volumes compared with unaltered streamflow conditions: (1) the median number of days per year in which the reservoir did not spill, (2) the median duration of the longest consecutive period of no-spill days per year, and (3) the lowest annual flow duration exceedance probability at which the outflows are significantly different from estimated unaltered streamflow at the 95-percent confidence level. Most reservoirs in the study do not spill during the summer months even under no-pumping conditions. The median number of days during which there was no spillage was less than 365 for all reservoirs in the study, indicating that, even under reported pumping conditions, the reservoirs refill to full volume and spill at least once during nondrought years, typically in the spring. Thirteen multiple-reservoir systems consisting of two or three hydrologically connected reservoirs were included in the study. Because operating rules used to manage multiple-reservoir systems are not available, these systems were simulated under two pumping scenarios, one in which water transfers between reservoirs are minimal and one in which reservoirs continually transferred water to intermediate or terminal reservoirs. These two scenarios provided upper and lower estimates of spillage under average pumping conditions from 2000 to 2004. For sites with insufficient data to simulate daily water balances, a proxy method to estimate the three spillage metrics was developed. A series of 4,000 Monte Carlo simulations of the reservoir water balance were run. In each simulation, streamflow, physical reservoir characteristics, and daily climate inputs were randomly varied. Tobit regression equations that quantify the relation between streamflow alteration and physical and operational characteristics of reservoirs were developed from the results of the Monte Carlo simulations and can be used to estimate each of the three spillage metrics using only the withdrawal ratio and the ratio of the surface area to the drainage area, which are available statewide for all reservoirs. A graphical user-interface for the Massachusetts Reservoir Simulation Tool was developed in a Microsoft Access environment. The simulation tool contains information for 70 reservoirs in Massachusetts and allows for simulation of additional scenarios than the ones considered in this report, including controlled releases, dam seepage and leakage, demand management plans, and alternative water withdrawal and transfer rules.

Massachusetts↗
USGS2016

Delineation of areas contributing groundwater to selected receiving surface water bodies for long-term average hydrologic conditions from 1968 to 1983 for Long Island, New York

To assist resource managers and planners in developing informed strategies to address nitrogen loading to coastal water bodies of Long Island, New York, the U.S. Geological Survey and the New York State Department of Environmental Conservation initiated a program to delineate a comprehensive dataset of groundwater recharge areas (or areas contributing groundwater), travel times, and outflows to streams and saline embayments on Long Island. A four-layer regional three-dimensional finite-difference groundwater-flow model of hydrologic conditions from 1968 to 1983 was used to provide delineations of 48 groundwater watersheds on Long Island. Sixteen particle starting points were evenly spaced within each of the 4,000- by 4,000-foot model cells that receive water-table recharge and tracked using forward particle-tracking analysis modeling software to outflow zones. For each particle, simulated travel times were grouped by age as follows: less than or equal to 10 years, greater than 10 years and less than or equal to 100 years, greater than 100 years and less than or equal to 1,000 years, and greater than 1,000 years; and simulated ending zones were grouped into 48 receiving water bodies, based on the New York State Department of Environmental Conservation Waterbody Inventory/Priority Waterbodies List. Areal delineation of travel time zones and groundwater contributing areas were generated and a table was prepared presenting the sum of groundwater outflow for each area.

New York↗
USGS2017

Evaluation of the streamgage network for estimating streamflow statistics at ungaged sites in Pennsylvania and the Susquehanna River Basin in Pennsylvania and New York

The current (2015) streamgage network in Pennsylvania and the Susquehanna River Basin in Pennsylvania and New York was evaluated in order to design a network that would meet the hydrologic needs of many partners and serve a variety of purposes and interests, including estimation of streamflow statistics at ungaged sites. This study was done by the U.S. Geological Survey, in cooperation with the Pennsylvania Department of Environmental Protection and the Susquehanna River Basin Commission. The study area includes the Commonwealth of Pennsylvania and the Susquehanna River Basin in Pennsylvania and New York. For this study, 229 streamgages were identified as reference streamgages that could be used to represent ungaged watersheds. Criteria for a reference streamgage are a minimum of 10 years of continuous record, minimally altered streamflow, and a drainage area less than 1,500 square miles. Some of the reference streamgages have been discontinued but provide historical hydrologic information valuable in the determination of streamflow characteristics of ungaged watersheds. Watersheds in the study area not adequately represented by a reference streamgage were identified by examining a range of basin characteristics, the extent of geographic coverage, and the strength of estimated streamflow correlations between gaged and ungaged sites. Basin characteristics were determined for the reference streamgage watersheds and the 1,662 12-digit hydrologic unit code (HUC12) subwatersheds in Pennsylvania and the Susquehanna River Basin using a geographic information system (GIS) spatial analysis and nationally available GIS datasets. Basin characteristics selected for this study include drainage area, mean basin elevation, mean basin slope, percentage of urbanized area, percentage of forested area, percentage of carbonate bedrock, mean annual precipitation, and soil thickness. A GIS spatial analysis was used to identify HUC12 subwatersheds outside the range of basin characteristics of the reference streamgages. There were 320 HUC12 subwatersheds, or 19 percent of the study area, with basin characteristics outside the range represented by the reference streamgage watersheds. A GIS spatial analysis was used to identify geographic gaps in the streamgage network. For each streamgage, a watershed area, called the gage statistical area (GSA), was delineated. The GSA shows the drainage area within a specific drainage-area ratio of the streamgage for transfer of streamflow statistics from that streamgage to ungaged sites on the valid statistical reach of the GSA for a streamgage. In Pennsylvania, a drainage-area ratio of 0.33–3 times the drainage area of the ungaged site was found to perform as well as, if not better than, more traditional ratios such as 0.5–1.5 (or 2) for transfer of selected streamflow statistics. A total of 1,102 HUC12 subwatersheds, or 66 percent of the study area, are outside the GSA for a reference streamgage. The USGS Baseline Streamflow Estimator (BaSE) program was used to determine how well HUC12 subwatersheds outside the streamgage GSAs are represented by the reference streamgage network in Pennsylvania, based on estimated streamflow correlation. The centroid of each HUC12 subwatershed was run through the BaSE program to determine the reference streamgage with the highest estimated streamflow correlation. There were 929 HUC12 subwatersheds in Pennsylvania, or 56 percent of the State, with an estimated correlation coefficient less than 0.96. The results from the basin characteristic, geographic, and streamflow correlation analyses were combined to identify 1,405 HUC12 subwatersheds in Pennsylvania and the Susquehanna River Basin in Pennsylvania and New York that lack a representative reference, based on at least one identified gap. Of the 1,405 HUC12 subwatersheds, 139 exhibited all three gaps, indicating a 8-percent gap in the reference streamgage network. Streamgages in areas with similar hydrologic characteristics and in close proximity to one another can potentially provide similar information (termed streamgages with high substitution potential). Streamgages were considered to have a high substitution potential with a nearby streamgage(s) if (1) the streamflow correlation coefficient was equal to or greater than 0.96, (2) the streamgages had 10 years of concurrent record, and (3) the streamgages are in the same watershed within the GSA of the streamgage. Seventy-four current (2015) streamgages with high substitution potential with at least one other streamgage were identified in the study area. Although these identified streamgages have a high substitution potential, they provide valuable streamflow information to a stakeholder. Selected primary uses of these streamgages were identified to determine the overall need for an individual streamgage.

New York, Pennsylvania↗
USGS2018

Rating stability, and frequency and magnitude of shifts, for streamgages in Virginia through water year 2013

The U.S. Geological Survey, in cooperation with the Virginia Department of Environmental Quality, has quantified several measures of rating stability and the frequency and magnitude of changes to ratings through time for 174 real-time continuous streamgages active in Virginia as of September 30, 2013. Generalized additive models (GAMs) were fitted through all available flow measurements for all the streamgages in Virginia’s real-time network as of September 30, 2013, with at least 20 flow measurements with positive flow values. For each measurement with a positive flow value, residuals from the GAM curve were calculated. Time series of these residuals were used to identify major changes to the control (the stream feature or features which control the relation between stage and flow); the measurements in the periods of equilibrium between major changes were assigned to rating families. Of the 127 rating families that were identified as being distinct at sites, documented explanations were found for 67 of them. The most common reasons for the control to change enough to warrant a new rating family are moving the streamgage (28 times), floods (26 times), and construction activities (13 times). Provisional flow data from any streamgage that has recently experienced a major flood, regardless of historical stability, are more uncertain than usual until post-flood evidence emerges that the rating is stable, or if the rating has changed, until it is known to be well defined. A direct comparison between provisional flow data (those data originally displayed on the web in near-real time) and flow data approved for publication following subsequent flow measurements and review could not be made because provisional flow data have not been archived. As a substitute, alternative flow (AltFlow) tables were constructed for periods with complete records of shifts and ratings. Alternate flows consist of $Q$ same , the flow value from the shifted rating table used to compute the daily flow value at the time of the most recent flow measurement that corresponds to the gage height of each day’s daily flow value, and $Q$ prev , the flow value from the shifted rating table in effect at the time of the previous flow measurement that corresponds to the gage height of each day’s daily flow value. Several metrics that summarize AltFlow tables were computed and evaluated; particular importance was given to how well the metrics agreed with the descriptive stability class developed from interviews with hydrographers. Of these stability metrics, at least four were determined to be meaningful and to represent different aspects of control stability that might be relevant to water managers: total root mean square error between log-transformed $Q$ prev and $Q$ same , percentage of days when the difference between $Q$ prev and $Q$ same is greater than (>) $Q$ same , the sum of absolute AltFlow error divided by total flow, and percentage of days with zero difference between $Q$ prev and $Q$ same . Three other meaningful metrics of control stability or provisional flow-data quality were computed: R 2 (coefficient of determination) of GAMs from the flow measurements, percentage of total estimated days, and percentage of estimated days in the winter. Correlations among metrics varied, indicating they responded to different aspects of control stability. Relations among the various stability metrics and quantitative basin and site characteristics were weak. Although quantitative relations between stability metrics and basin and site characteristics were all weak, some common patterns still emerged. Controls and ratings on large streams [>500-square mile (mi 2 ) drainage area] and at high elevations (>1,000 feet) were more likely to be stable than controls and ratings on small streams (less than (<) 100-mi 2 drainage area) and at low elevations (<500 feet). There were exceptions to both generalizations, and streamgages that were intermediate in both characteristics varied widely in stability. Typical timing of record computation changed during water years 1991–2013. From 1991 through 2001, the median number of days between the start date of the shift and the date it was created fluctuated between about 240 and about 300 days and decreased by about 4 months from 2001 to 2002. Only in 2012 and 2013 did one-half of the shifts have a delay of about 60 days between start date and final modification.

Virginia↗
USGS2018

Revised groundwater-flow model of the glacial aquifer system north of Aberdeen, South Dakota, through water year 2015

The city of Aberdeen, in northeastern South Dakota, requires an expanded and sustainable supply of water to meet current and future demands. Conceptual and numerical models of the glacial aquifer system in the area north of Aberdeen were developed by the U.S. Geological Survey in cooperation with the City of Aberdeen in 2012. The U.S. Geological Survey, in cooperation with the City of Aberdeen, completed a study to revise the original numerical groundwater-flow model using data through water year (WY) 2015 to aid the City of Aberdeen in their development of plans and strategies for a sustainable water supply and to increase understanding of the glacial aquifer system and groundwater-flow system near Aberdeen. The original model was revised to improve the fit between model-simulated values and observed (measured or estimated) data, provide greater insight into surface-water interactions, and improve the usefulness of the model for water-supply planning. The revised groundwater-flow model (hereafter referred to as the “revised model”) presented in this report supersedes the original model. The purpose of this report is to describe a revised groundwater-flow model including data collection, model calibration, and model results for the glacial aquifer system including the Elm, Middle James, and Deep James aquifers north of Aberdeen, South Dakota, using updated hydrologic data through WY 2015. The original numerical model was revised in several ways. The model was modified by adding four new layers, which included a surficial layer, two intervening confining layers, and a shale bedrock layer. The revised model provides an improved understanding of the groundwater-flow system in comparison to the original model. The principal aquifers of the model area include portions of the Elm, Middle James, and Deep James aquifers. The lithologic information used to define and describe the aquifers in the model area was unaltered; however, aquifer properties and boundary conditions were reviewed and updated using geological information reported by the South Dakota Department of Environmental and Natural Resources and information obtained from geophysical investigations for this study. The horizontal extent of the Elm, Middle James, and Deep James aquifers was unaltered from the original model. The thickness of the Deep James aquifer was modified based on interpretations from the geophysical investigations. In general, groundwater in the Elm aquifer flowed from northwest to southeast and locally towards rivers and streams. Similarly, in the Middle James and Deep James aquifers, groundwater also typically flowed southeast. The revisions made to the original model include use of the following MODFLOW stress packages: Recharge, Evapotranspiration, Time-Variant Specified Head, Wells, Drains, and Stream Flow Routing, all of which were updated from the original model except for the Stream Flow Routing Package, which replaced the River Package used in the original model. Model calibration is the process of estimating model parameters to minimize the differences, or residuals, between observed data and simulated values; therefore, Parameter ESTimation (PEST) software was used to optimize model input parameters by matching model-simulated values to observed data. Calibration parameters included horizontal hydraulic conductivity, vertical hydraulic conductivity, specific yield, specific storage, and vertical streambed conductance for stream and drain cells. Multipliers were used to calibrate the recharge and evapotranspiration stresses. Evapotranspiration extinction depth also was adjusted during model calibration. Comparisons to the original model are described to highlight the changes made in the revised model. In general, the revised model adequately simulates the natural system and compares favorably with observed hydrologic data. Simulated water levels were evaluated by comparing them to single water-level observations at selected well locations. The selected wells were the same wells used in the original model. The coefficient of determination value between simulated and observed water levels for the revised model was 0.89 and included simulated and observed values from October 1, 1974 (WY 1975), through September 30, 2015 (WY 2015). The coefficient of determination value for the original model was 0.94 and included simulated and observed values from October 1, 1974, through September 30, 2009. The difference may indicate that the original model could have been overfit to hydraulic head observations because base flow was not simulated. The additional data used in the revised model included some climatically wetter, more extreme periods, such as 2011, in which annual precipitation was 30.9 inches. Average annual precipitation for the original model timeframe, which included data from WYs 1975–2009, was 20.26 inches. Additional precipitation data for WYs 2010–15, included in the revised model timeframe, resulted in an average annual precipitation for WYs 1975–2015 in the model area of 20.6 inches. The larger variability in climate data coupled with the additional water-level data could explain the lower coefficient of determination for water levels in the revised model. The revised model was used to calculate various groundwater-budget components for steady-state and transient conditions for WYs 1975–2015. The time-variant specified-head cells in the revised model had the largest change when compared to the original steady-state model for inflows and outflows. Comparing the transient budget components between the original and the revised models indicated that inflow from recharge and time-variant specified-head cells had the greatest effect on groundwater inflows, and outflow from storage had the greatest effect on groundwater outflows. The simulated potentiometric contours from the revised model were compared with (1) the observed (interpreted) potentiometric surface (layer 2) and the hydraulic head values (layers 4 and 6) and (2) the simulated contours from the original model. The simulated hydraulic gradients and general direction of groundwater flow in the Elm aquifer in the revised model generally matched the observed potentiometric contours, the simulated potentiometric contours from the original model, and general flow directions interpreted to be perpendicular to the contours. Minor discrepancies between simulated potentiometric contours from the revised model and the observed potentiometric contours may be due to the lack of observed data in the model area. The revised model was designed to reduce the limitations of the original model. The revisions were validated by comparing the results of the original model with the revised model. A primary benefit of the revised model is the inclusion of the surficial deposits and the confining units as explicit layers in the model. The addition of the surficial layer was beneficial for three primary reasons: (1) more accurate representation of recharge from precipitation, (2) more accurate representation of groundwater evapotranspiration, and (3) more accurate representation of groundwater and surface-water interactions. The groundwater model is a numeric approximation of a complex physical hydrologic system, and the revised model data were interpolated in regions with sparse data. Additionally, model discretization included averaged and interpolated values for water use, withdrawal rates, and hydraulic conductivity. The revised model provides a useful estimate for hydraulic gradients, groundwater-flow directions, and aquifer response to groundwater withdrawals.

South Dakota↗
USGS2019

Estimation of base flow on ungaged, periodically measured streams in small watersheds in western Pennsylvania

A 2.5-year data collection program was undertaken by the U.S. Geological Survey, in cooperation with the Pennsylvania Department of Environmental Protection (PADEP), to quantify and estimate base flow in small watersheds in western Pennsylvania where only periodic streamflow measurements had been obtained. Twelve streamgages with watershed areas of less than 10 square miles were established in western Pennsylvania for this study, with most established within Greene and Washington Counties (an area where a type of underground coal mining known as longwall mining occurs). Data from five previously established streamgages with watershed areas ranging from 48.9 to 281 square miles were also used in the analyses for this study. The index-gage method was used to relate streamflow at one streamgage referred to as the “index streamgage” to streamflow at another site of interest (usually an ungaged site, but for this study another streamgage) using a regression technique. Streamflow regressions were developed for all newly established streamgages by using the Maintenance of Variance Extension, Type 1 (MOVE.1) method. Not all streamflow data from the newly established streamgages were used for MOVE.1 regression development; only data that have little to no influence from runoff were considered. Runoff-influenced streamflow for this study was defined as streamflow on a day that precipitation occurs plus streamflow on the following 2 days. One streamflow value per day selected from a specified schedule that captures numerous non-runoff periods was used to develop a MOVE.1 regression. Prediction limits were calculated from the regression to provide the upper and lower bounds for the regression-produced streamflow estimates. Using these data, base flow at a site can be estimated with the index-gage method. The log 10 -tranformed prediction interval width and other regression diagnostics were used as indicators of regression quality when comparing streamgage relations to determine the best index streamgage among the streamgages established for this study. It was determined that index streamgages within about 10 miles of the site of interest provided the best estimated base flow and could, in the future, be used by mine operators and the PADEP to quantify base flow and to evaluate the effects of mining on streamflow.

Pennsylvania↗
USGS2019

Regression models for estimating sediment and nutrient concentrations and loads at the Kankakee River, Shelby, Indiana, December 2015 through May 2018

The Kankakee River in northern Indiana flows through the area once known as the Grand Marsh. Beginning in the 1860s, anthropogenic changes to the river within Indiana resulted in downstream flooding and additional transport of sediment and nutrients. In 2015, the U.S. Geological Survey, in cooperation with the Indiana Department of Environmental Management, upgraded the gaging station Kankakee River at Shelby, Indiana, to include the collection of water-quality data. By relating continuously monitored water-quality data to discrete data collected from December 2015 through May 2018, linear regression was used to develop models for estimating concentrations of suspended sediment, total nitrogen, and total phosphorus. Developed regression models indicated a strong correlation between turbidity and specific conductance with suspended-sediment concentration (adjusted coefficient of determination equals 0.92, predicted residual error sum of squares equals 0.151), nitrate plus nitrite and specific conductance with total nitrogen (adjusted coefficient of determination equals 0.95, predicted residual error sum of squares equals 0.0248), and turbidity with total phosphorus (adjusted coefficient of determination equals 0.89, predicted residual error sum of squares equals 0.0103). Daily loads of suspended sediment, total nitrogen, and total phosphorus were computed as the product of daily mean regression model concentrations and daily mean streamflow. Rloadest models were used to compute daily loads of each constituent during gaps in regression model loads. For 2016 and 2017, the estimated annual suspended-sediment loads were 105,000 and 91,000 tons; estimated total nitrogen loads were 8,690 and 8,890 tons; and estimated total phosphorus loads were 265 and 236 tons, respectively.

Indiana↗
USGS2020

Quantification of trace element loading in the upper Tenmile Creek drainage basin near Rimini, Montana, September 2011

The principle sources of trace elements entering upper Tenmile Creek, Montana, during September 2011, four trace metals and the metalloid arsenic, were identified and quantified by combining and analyzing streamflow data determined from tracer injection with trace-element concentrations and related water-quality data determined from synoptic sampling. The study reach was along upper Tenmile Creek, beginning downstream from the city of Helena’s diversion and extending 5,020 feet downstream. Results from the 2011 study, completed by the U.S. Geological Survey in cooperation with the Montana Department of Environmental Quality, were compared to results from a similar study conducted in 1998 to assess the effectiveness of mine reclamation and remediation work to reduce trace-element loading to upper Tenmile Creek, which has been ongoing throughout the drainage basin. Main-stem concentrations of most trace elements analyzed were generally greater in 1998 than in 2011. However, the State of Montana human-health criteria for total-recoverable cadmium and arsenic were exceeded in parts of upper Tenmile Creek, and concentrations of cadmium and zinc exceeded the acute aquatic-life criteria at all main-stem sites during both studies. Total-recoverable copper concentrations observed in 2011 exceeded the chronic aquatic-life criterion upstream from the Lee Mountain adit, whereas, in 1998, all sites exceeded the acute aquatic-life criteria. Direct comparison of loads from the 1998 and 2011 tracer studies were complicated by the differences in hydrologic conditions. Streamflow in 1998 was about 10 percent of the 2011 streamflow. The Lee Mountain Mine and Susie Lode adit were identified as major contributors of trace elements to upper Tenmile Creek in both studies. However, trace-element loading from the Lee Mountain Mine area was substantially reduced between 1998 and 2011. Total-recoverable loads of all trace elements showed substantial loss in 1998 but increased in 2011 downstream from the Susie Lode adit to the end of the study reach. This reach was one of the primary sources of trace-element loading to upper Tenmile Creek in 2011. This difference indicated that the streambed may act as a sink or a source for trace elements, depending on hydrologic conditions.

Montana↗
USGS2020

Updating data inputs, assessing trends, and evaluating a method to estimate probable high groundwater levels in selected areas of Massachusetts

A method to estimate the probable high groundwater level in Massachusetts, excluding Cape Cod and the islands, was developed in 1981. The method uses a groundwater measurement from a test site, groundwater measurements from an index well, and a distribution of high groundwater levels from wells in similar geologic and topographic settings. The U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, conducted an update to the Frimpter method for estimating the probable high groundwater levels in Massachusetts. The study evaluated the potential changes to the method resulting from four decades of additional groundwater-level data and the expansion of the network of wells for monitoring groundwater levels. The differences and potential benefits of daily, as opposed to monthly, measurements in the application of the method were examined because of the increased availability of high-frequency (subdaily) groundwater-level data. The study also considered long-term trends in groundwater levels that may alter the accuracy of the method. Finally, the accuracy of the estimated high groundwater levels was evaluated, and improved implementation guidance was prepared. For this study, groundwater levels in 153 wells in Massachusetts and surrounding States with records with lengths of 16 to 78 years were analyzed. The highest recorded groundwater levels ranged from 1.2 feet (ft) above land surface (flooded conditions) to 45.8 ft below land surface, with a median of 4.6 ft below land surface. The maximum annual groundwater-level range was 1.4 to 17.9 ft, with a median of 5.5 ft. The within-month variation, maximum annual groundwater-level range, and highest recorded groundwater level were computed using daily mean groundwater-level values from 28 wells with continuous records. The use of daily data resulted in larger maximum annual groundwater-level ranges (0.02 to 2.94 ft larger, with a median of 0.58 ft larger) and shallower highest-recorded groundwater levels (0.0 to 1.60 ft shallower, with a median of 0.18 ft shallower) than computations based on monthly measurements in the same wells. Statistical tests showed moderate to strong evidence of trends in measurements of both high and low groundwater levels within most of the periods during which water levels were analyzed. High groundwater levels rose beneath the land surface at most sites during four of the six periods used for analysis (1966–2015, 1986–2015, 1991–2010, and 1981–2010). Low groundwater levels also increased at many sites during most of the periods evaluated, but this trend was less widespread than the similar trends in high groundwater levels, and the trend was to deeper low groundwater levels at more sites than the trend to deeper high groundwater levels. There was no clear trend in annual groundwater-level ranges at most sites during the six periods analyzed. In general, the Frimpter method predicted shallower (higher) high groundwater levels than were observed but correctly classified sites according to their suitabilities for unmounded septic systems. The mean error of the predictions (difference between the estimated and observed groundwater levels) ranged from −3.23 ft to −1.40 ft for various approaches to estimating the groundwater-level range and selecting an index well. The method correctly classified 83 to 86 percent of monitoring-well sites according to their suitability for an unmounded septic system for many approaches to estimating the annual groundwater-level range and selecting an index well. The approach selected for estimating the annual groundwater-level range and selecting an index well will depend upon the importance of an accurate estimate of the high groundwater level as compared to the importance of an estimated high groundwater level that is less likely to be exceeded.

Connecticut, Massachusetts, New Hampshire, Rhode I↗
USGS2021

Hydrologic and ecological investigations in the School Branch watershed, Hendricks County, Indiana—Water years 2016–2018

School Branch in Hendricks County in central Indiana, is a small stream with a variety of agricultural and suburban land uses that drains into the Eagle Creek Reservoir, a major source of drinking water for Indianapolis, Indiana. The School Branch watershed has become the focus of a collaborative partnership of Federal, State, and local agencies; a university research center; and agricultural producers to understand the effects of land use and management practices on water quality and water quantity in the watershed. The U.S. Geological Survey, in cooperation with the Indiana Department of Environmental Management, contributed to the School Branch partnership with the operation of three streamgages (03353415 School Branch at Maloney Road near Brownsburg, Indiana; 03353420 School Branch at County Road 750 North at Brownsburg, Indiana; and 03353430 School Branch at Noble Drive at Brownsburg, Indiana) and the operation of a continuous water-quality gage (also known as a supergage) at County Road 750 North that measured dissolved oxygen, pH, temperature, specific conductance, turbidity, nitrate, and orthophosphate. Additional efforts included the use of passive samplers to identify wastewater indicators; assessment of fish and macroinvertebrate communities and stream habitat to identify ecological impairment; sampling for nutrients and sediment to estimate loads; and using major ions, stable isotopes and nested groundwater monitoring wells at County Road 750 North to determine hydrologic connectivity between the groundwater and surface water. The objectives of this study were to collect surface and groundwater data to analyze the hydrology and water quality within the watershed. Total nitrogen yields were highest at the upstream site, Maloney Road, and indicated a mixture of nitrogen sources in the watershed. Differences found in total nitrogen loading patterns throughout the watershed may be linked to differences in hydrology and land-use management from site to site. The groundwater and surface water were shown to be highly connected, and except for some low-flow periods, the water was flowing from groundwater to the stream for most of the study period. Fish and macroinvertebrate communities show improvement from upstream to downstream, with increases in diversity, richness, and species sensitive to poor water quality and habitat. These increases were most likely due to improved habitat quality at the downstream station.

Indiana↗
USGS2023

Hydrogeology of sand and gravel aquifers in the Owasco Inlet watershed, Cayuga and Tompkins Counties, New York

This study is a continuation of a series of hydrogeologic appraisals that have been conducted since 1980, as part of a cooperative, long-term, detailed aquifer mapping program by the U.S. Geological Survey and the New York State Department of Environmental Conservation. These appraisals provide a foundation for wellhead protection programs, water-resource management and planning decisions, and groundwater remediation in upstate New York. The Owasco Inlet watershed drains north directly to Owasco Lake, one of New York’s Finger Lakes. The watershed is similar in form to other watersheds of the Finger Lakes region of New York State and is characterized by broad, smooth uplands punctuated by the 18-mile-long Owasco Inlet valley and the secondary valleys of Decker, Dresserville, and Hemlock Creeks. All these streams occupy “through valleys”—the valleys are continuous into adjacent watersheds such that watershed divides are high points within the valley bottoms, rather than in uplands. Most glacial deposits in the watershed were deposited during the Valley Heads Readvance and subsequent retreat (with at least one minor readvance). Estimates are that the Valley Heads Readvance likely peaked in the Cayuga basin about 17,000 (calendar) years ago based on dates from western New York and the Mohawk Valley in east-central New York. It was the last major ice advance of the Pleistocene Epoch in this region. This readvance covered the entire Owasco watershed. Valley Heads ice in the Owasco watershed was part of the Cayuga ice lobe, which flowed parallel to the Owasco Inlet valley as far south as Moravia, but which became progressively more eastward south of Locke, effectively raking across the valley from the west. The Owasco Inlet watershed drains north, like other Finger Lakes watersheds, and during deglaciation, northward meltwater drainage was largely blocked by the ice, which resulted in development of a series of regional proglacial lakes in which fine lacustrine sediments predominated. Lacustrine sediments constitute much of the Owasco Inlet valley fill, but it has been pointed out that there are no obvious outlet channels exiting the Owasco Inlet watershed, and it has been proposed that there was northward drainage of meltwater into the ice within the larger, neighboring Cayuga Trough. The lacustrine deposits in the Owasco Inlet valley are commonly underlain by thinner coarse-grained stratified deposits (subaqueous fans and eskers) of variable sorting and permeability. An exception occurs at Groton, N.Y., where retreating ice paused long enough for coarse-grained sediments to fill much of the valley. Smaller, higher elevation valleys have a wider variety of glacial valley-fill deposits, ranging from fine lacustrine sediments to sand and gravel to till. Groundwater is the sole source of water supply in the area; glacial sand and gravel aquifers are the primary water source in the Owasco Inlet and Decker Creek valleys, and fractured bedrock aquifers typically supply domestic wells in the remaining valleys and upland areas. Municipal supplies tap a variety of aquifer types in the Owasco Inlet valley and nearby uplands. The hamlet of Locke taps the extensive confined aquifer beneath fine-grained lacustrine deposits in the valley. The Village of Moravia taps a semiconfined aquifer that overlies proglacial lake deposits and is partly confined by overlying recent lake deposits and alluvium. Withdrawals from this aquifer may also induce water from the Owasco Inlet into the aquifer. The Village of Groton draws from two aquifers: (1) a thin upland unconfined sand and gravel aquifer tapped by an infiltration gallery and (2) the local unconfined or semiconfined aquifer in the Owasco Inlet valley. Another aquifer with potential for municipal supply is the unconfined aquifer in the Decker Creek valley near Wilson Corners. The confined aquifer in the lower Dresserville Creek valley may have some water-resource potential, but it is largely untested. Unconfined aquifers are the most susceptible to contamination from activities at land surface directly above the aquifer because precipitation and subsequent recharge can transport contaminants directly to the water table. Adjacent upland areas can also contribute contaminants. Confined aquifers are less susceptible to contamination from overlying land surface areas because confining units largely prevent downward movement of water. Recharge occurs elsewhere at unconfined upvalley locations and along valley walls where alluvial fans, ice-contact deposits, or stream incision into the valley wall may provide pathways for downward movement of groundwater.

New York↗
USGS2024

Water-quality monitoring strategy for Mount Hope Bay and the Taunton River Estuary, southeastern Massachusetts

The U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection (MassDEP), began a study in 2018 to develop a water-quality monitoring strategy (WQMS) for Mount Hope Bay and the Taunton River Estuary in southeastern Massachusetts. MassDEP is interested in water-quality data in Mount Hope Bay and the Taunton River Estuary to characterize current water-quality conditions, assess nutrient-related effects, and capture conditions before and after planned upgrades to wastewater treatment facilities. The WQMS provides an overview of the environmental setting of Mount Hope Bay and the Taunton River Estuary and of dissolved oxygen and nutrient-related water-quality issues, reviews historical and existing monitoring data, and provides recommendations for future monitoring designed to meet five MassDEP management objectives: (1) support MassDEP’s review of coastal and marine dissolved oxygen criteria, (2) assess conditions within MassDEP waterbody assessment units in Mount Hope Bay and the Taunton River Estuary with respect to selected criteria in the Massachusetts Surface Water Quality Standards, (3) assess conditions along the freshwater/saltwater interface (salt wedge) of the Taunton River Estuary and delineate the boundary between the freshwater and saltwater waterbody assessment units, (4) estimate data requirements needed to determine nutrient loads flowing into Mount Hope Bay and the Taunton River Estuary, and (5) evaluate data requirements needed to support hydrodynamic and water-quality models for Mount Hope Bay and the Taunton River Estuary. This WQMS may be used by MassDEP to develop a statewide approach for monitoring estuaries in Massachusetts.

Massachusetts↗
USGS2026

Assessment of natural gas pipeline construction on stream temperature and turbidity in southwestern Virginia, 2017–25

Despite the extensive natural gas pipeline network in the United States that intersects streams and other sensitive habitats, few case studies use a comparative upstream-to-downstream approach to evaluate potential short- and long-term effects of pipeline stream crossings from pre-construction through post-restoration. In 2017, the U.S. Geological Survey, in cooperation with the Virginia Department of Environmental Quality, deployed real-time continuous stream monitoring stations upstream and downstream from six proposed Mountain Valley Pipeline stream crossings in southwestern Virginia. Water temperature and turbidity data collected at the upstream and downstream sites were compared across three periods—before stream crossing construction, during stream crossing construction, and after stream crossing construction—to determine potential influences from the pipeline stream crossing. Additionally, the monitoring network was used to notify regulators of potentially anomalous conditions throughout the entire monitoring period. The results of this study indicate that pipeline stream crossing did not affect long-term or short-term upstream-to-downstream water temperature conditions or long-term upstream-to-downstream turbidity conditions in any of the six monitored streams. Some short-term anomalously elevated turbidity conditions were observed and attributable to pipeline stream crossing; however, the magnitudes and durations were not sufficient to alter the long-term turbidity regimes of the streams in which they were observed. The application of the monitoring network as a real-time alert system successfully alerted regulators to potentially anomalous conditions.

Virginia↗
USGS2026

Phosphorus concentrations and loads in two Mississippi Alluvial Plain watersheds following nutrient reduction efforts during 2010–22

The Mississippi Alluvial Plain ecoregion in Mississippi was historically part of the largest continuous wetland in North America but was clearcut and converted to row crop agriculture beginning in the late 1800s. To reduce the effects of agriculture on surface waters in the Mississippi Alluvial Plain, nutrient reduction best management practices have been implemented throughout the region, particularly in two small watersheds: Harris Bayou watershed (174.04 square kilometers) and Porter Bayou watershed (267.15 square kilometers). Since 2009, the U.S. Geological Survey, in cooperation with the Mississippi Department of Environmental Quality, has collected water-quality samples monthly and during high-flow runoff events, at three locations within each watershed to monitor and assess nutrient trends following best management practice implementation. Here, we have analyzed long-term and seasonal trends in discharge, total nitrogen (TN), and total phosphorus (TP) concentrations and loads using a combination of quantile-Kendall analysis, quantile regression, generalized additive models, and weighted regression on time, discharge, and season. In Harris Bayou watershed, discharge values at a ~0.4–0.5 daily nonexceedance probability increased by 15 percent per year at the uppermost monitoring site. Additionally, TN concentrations decreased throughout the watershed over time, particularly in the spring and winter months. Reductions in TP were less apparent and only observed at one sampling location in Harris Bayou watershed. Nutrient concentrations in Porter Bayou watershed remained consistent over time and increased at one sampling location. However, despite discrepancies in concentration trends, flow-normalized TN and TP loads consistently decreased over time in both watersheds. At the mouth of both watersheds, peak TN and TP loads decreased between 37.7 and 56.9 percent. Collectively, these results provide evidence that best management practices may help reduce the effect of agricultural practices in the region, particularly in Harris Bayou watershed. Additionally, management efforts resulted in greater reductions in TN concentrations than in TP concentrations.

Scientific Investigations Report↗
USGS2026

Variability of discharge, nutrients, and Escherichia coli from tile drains in a small agricultural stream, a synoptic study in the School Branch watershed, Hendricks County, Indiana

The U.S. Geological Survey, in cooperation with the Indiana Department of Environmental Management, sampled a 1-mile stream reach of School Branch in Hendricks County, Indiana. Sampling consisted of 24 tile drain sites and surface-water sites both upstream and downstream from the tile drains. This study intended to (1) assess the variability of discharge, nutrients, and Escherichia coli ( E. coli ) among tiles within a stream reach in variable conditions; (2) define how the tiles affect the stream hydrology and water quality; and (3) determine potential management implications. Samples were collected on three dates—round 1 on June 12, 2018, after fertilizer was applied; round 2 on September 11, 2018, in fall harvest conditions; and round 3 on April 16, 2019, before fertilizer was applied—to encompass differing precipitation events, multiple seasons, and field conditions. Sites were sampled for water-quality parameters, dissolved and total nutrients (ammonia, nitrate plus nitrite, total nitrogen, orthophosphate, and total phosphorus), E. coli , and instantaneous discharge. Instantaneous nutrient loads were calculated from tiles and at the upstream and downstream surface-water sites. Nutrient concentrations and instantaneous loads varied depending on the tile and the sampling date. Tiles contributed between 26.5 and 66.0 percent of increased discharge between the surface-water sites. Tiles contributed between 34.7 and 87.4 percent and between 19.7 and 65.9 percent of the increased total nitrogen and phosphorus loads, respectively, for the downstream site relative to the upstream site. Most of the total nitrogen and total phosphorus from the tiles was in the dissolved forms (nitrate plus nitrite and orthophosphate). Surface-water samples had a nitrate plus nitrite to total nitrogen ratio of 71 to 93 percent, with an average of 82 percent; tile samples ranged from 44 to nearly 100 percent, with an average of 91 percent dissolved nitrogen. Surface-water samples had an orthophosphate to total phosphorus ratio of 30 to 65 percent, with an average of 53 percent; tile samples ranged from 32 to nearly 100 percent, with an average of 85 percent dissolved phosphorus. E. coli was detected in all but 5.79 percent of tile and surface-water samples. The E. coli results indicated the possibility of human or animal fecal contamination by way of septic systems leaking or directly linked to the tiles draining into the study area. Consequently, five microbial source tracking samples were collected from the upstream and downstream surface-water sites between August 3, 2021, and May 18, 2022. Human, canine, ruminant, and avian sources of fecal contamination were evaluated. The microbial source tracking samples indicated that human-specific sources are likely consistently contributing to elevated fecal contamination. Differences in discharge, nutrient concentrations, nutrient loads, and E. coli concentrations show that many factors affect how tiles alter the hydrology and water quality of receiving waters. The timing and severity of precipitation events, season, and the antecedent conditions of the field or watershed before sampling can affect the results. Tiles with the largest discharge rate had the most potential effect on the water quality and quantity of the stream. Targeting high-flowing tiles may affect efforts to improve overall water quality.

Indiana↗